IR photometry and dust-shell models for two carbon stars

We present JHKLM photometry of the carbon stars ST And and T Lyn acquired in 2000–2010. Along with brightness variations due to pulsations, changes on timescales of 2000–3000 days are also observed. Our combined light curves can be satisfactorily represented with light elements derived from visual observations, but the maxima are delayed relative to the calculated times. A color-index analysis demonstrates that the dust shell of ST And is fairly weak, and is manifest only episodically, while the presence of hot dust was always detected for T Lyn. These results confirm models of spherically symmetric stellar dust shells based on mean-flux data, supplemented with observations in the intermediate IR from the IRAS and AKARI satellites. The visual optical depth of the relatively cool dust shell of ST And assuming a dust temperature at the inner edge of T ₁ = 510 K is very low: τ _V = 0.047. The dust shell of T Lyn is considerably hotter (T ₁ = 940 K), with τ _V = 0.95. We estimate the mass-loss rate to be 1.8 × 10⁻⁷ M _⊙/year for ST And and 3.7 × 10⁻⁷ M _⊙/year for T Lyn.     

The influence of the decay of OB associations on the evolution of dwarf galaxies

The ejection of stars from spheroidal and disk dwarf galaxies resulting from the decay of OB associations is studied. This has substantial observational consequences for disk galaxies with escape velocities up to 20 km/s, or dynamical masses up to 10⁸ M _⊙. The ejection of stars can (i) reduce the abundances of the products of Type Ia supernovae and, to a lesser degree, Type II supernovae, in disk stars, (ii) chemically enrich the galactic halo and intergalactic medium, (iii) lead to the loss of 50% of the stellar mass in galaxies with masses ∼10⁷ M _⊙ and the loss of all stars in systems with masses ≲10⁵ M _⊙, (iv) increase the mass-to-luminosity ratio of the galaxy.     

Results of infrared photometry and a model for the dust envelope of the protoplanetary nebula candidate V1027 Cyg

Results of JHKLM photometry for the protoplanetary nebula candidate V1027 Cyg obtained in 1991–2008 are reported. In all bands, the brightness variations did not exceed 0.2 ^m . Estimated linear trends demonstrate no significant changes in the mean brightness and color indices of the object, with the possible exception of the L-M color index, which showed a small decrease. A search for possible periodicities in the brightness variations yielded the most probable period of 237^d. A model for a spherically-symmetric dust shell has been calculated based on the photometric results supplemented with data on the mid- and far-infrared fluxes. The estimated mass-loss rate of the star is 1.3 × 10⁻⁵ M _⊙/year.     

Map of the 6.7 GHz class II methanol maser emission of the protoplanetary disk G23.01-0.41

We present images of the star-forming regionG23.01–0.41 at 6.7GHz in the Class II methanol maser transition 5₁–6₀ A ⁺, produced from archival observations on the European VLBI Network. Our map of the source and its maser spots contains 24 maser components. The data for each spot—absolute coordinates, coordinates relative to the calibration feature, peak flux and flux integrated over the spot, size, position angle, velocity along the line of sight, and line full width at half-maximum—are collected in tabular form. The spatial region occupied by the maser spots is approximated by a 200×130 milliarcsec ellipse in position angle PA = −0.40°, centered on the absolute coordinates α ₀ = 18^h34^m40.282^s, δ ₀ = −09°00′38.27″ (J2000). If the source is a protoplanetary disk, then, for the distance estimate derived from trigonometric parallax, its diameter is 1800 AU, and the mass of the central protostar is 23.5M _⊙.     

Search for and study of hot circumstellar dust envelopes

Long-term (1984–2008) JHKLM photometry for 254 objects is presented. The observations were carried out in the standard JHKLM photometric system using an original method and a modern IR photometer designed and built at the Sternberg Astronomical Institute. Our program of studies included searches for and studies of relatively hot circumstellar dust envelopes. The most important results obtained using these observations include the following. We have detected relatively hot dust envelopes in a number of objects for the first time, including the RCB star UV Cas, RX Cas, several classical symbiotic stars, etc. A model has been calculated for the dust envelope of FG Sge, which formed around the star as a result of several successive cycles of dust condensation beginning in Autumn 1992. Several dust-condensation episodes have been traced in the envelopes of symbiotic systems (CH Cyg, V1016 Cyg, HM Sge, etc.), as well as the role of the hot component in the formation of the dust envelopes. We have established from variations of the IR emission that the cool components in the symbiotic novae V1016 Cyg and HM Sge, and possibly CH Cyg, are Miras. The binarity of V1016 Cyg and HM Sge has also been firmly established. The variability of a whole series of object has been studied, including the stellar components of close binary systems and several dozen Mira and semi-regular variables. The ellipsoidality of the components in the RX Cas system (a prototype W Ser star) and the cool component in the symbiotic systems CI Cyg and BF Cyg has been firmly established. We have obtained the first IR light curve for the eclipsing system V444 Cyg (WN5+O6), and determined the wavelength dependence of the obtained parameters of the WN5 star. Analysis of the IR light curves of several novae indicate the condensation of dust envelopes in the transition periods of Cygnus 1992, Aquila 1993, and Aquila 1995. The IR light curve of R CrB has been obtained over a long period and analyzed. IR observations of the nova-like variable V4334 Sgr have been carried out over four years, over which the star passed through four stages during its motion along its post-AGB evolutionary track; the star’s bolometric flux and optical depth of its dust envelope have been estimated, and the structure and mass of the dust layer determined. We have analyzed the IR variability of the symbiotic star V407 Cyg over 14 years, and found its cool component to be a Mira with a period of 745 days. The observed pulsations and trend are associated with the luminosity and temperature variations of the Mira, as well as the optical depth of the dust envelope. The size of the dust grains and mass-loss rate of the Mira have been determined. We have obtained JHKL light curves for the Seyfert galaxy NGC 4151 over 23 years. The IR brightness of the galaxy grew from 1985 through 1996 (by ∼0.9 ^m at 1.25 μm, ∼1 ^m at 1.65 μm, ∼1.1 ^m at 2.2 μm, and ∼1.3 ^m at 3.5 μm), while the galaxy simultaneously reddened. The “cool” variable source in NGC 4151 was still in the active state in 1998, although its luminosity had decreased by approximately 15%–20%. If the “cool” component of the variable source in this galaxy is a dust envelope heated by the central “hot” source, it should be optically thin to the radiation of this source: its mean optical depth is in the range 0.05–0.15. Emission from dust particles heated to temperatures of 600–800 K was observed in the near IR at a distance of several parsecs from the nucleus during the period of activity in 1995–1998; the inferred mass of emitting dust was 5–20 M _⊙. In 1994–2003, we observed a tendency for NGC 4151 to become bluer at 1.25–1.65 μm while simultaneously reddening at 2.2–3.5 μm. Beginning in Autumn 2000, the galaxy began to emerge from a minimum, which lasted from March 2000 through April 2001 in the IR; a flare of the galactic nucleus was observed and followed in detail in the IR in this same period. We confirm the IR variability of the nucleus of the Seyfert galaxy NGC 1068, which can be located in various stages of activity. The variability of NGC 1068 is associated with a complex source. A periodic component has been detected in the J brightness variations of the oxygen Mira V2108 Oph; we have shown that this star is immersed in a fairly dense dust envelope, and have calculated amodel for this envelope. We have calculated an (axially symmetrical) dust-envelope model for the carbon semi-regular variable RWLMi with a density distribution characteristic for the “superwind” stage. This envelope model is able to reproduce the observed fluxes over the entire observable spectrum, and displays a good agreement with the observational data. The three hot supergiants V482 Cas, QZ Sge, and HD 338926 may be variable in the IR. Long-term photometry of eight planetary nebulae in the near IR (1.25–5 μm) has enabled us to firmly establish the IR brightness and color variability of these nebulae on time scales from several tens of days to six-to-eight years. We have analyzed long-term JK photometry of the X-ray binary Cyg X-1 (V1357 Cyg). Periodicity with a characteristic time scale of ∼11.5 years is visible in the JK brightness variations of Cyg X-1 in 1984–2007, possibly due to periodic variations of the temperature, radius and luminosity of the optical component of the binary with P ∼11.5 yr. Fourteen-year IR light curves of five RV Tau stars (R Sge, RV Tau, AC Her, V Vul, and R Sct) and the yellow supergiant and protoplanetary-nebula candidate V1027 Cyg have been analyzed. A spherically symmetrical dust-envelope model has been calculated for V1027 Cyg.     

Do we have enough pieces of the jigsaw to integrate CO2 fluxes in the coastal ocean?

Annually integrated air-water CO₂ flux data in 44 coastal environments were compiled from literature. Data were gathered in 8 major ecosystems (inner estuaries, outer estuaries, whole estuarine systems, mangroves, salt marshes, coral reefs, upwelling systems, and open continental shelves), and up-scaled in the first attempt to integrate air-water CO₂ fluxes over the coastal ocean (26×10⁶ km²), taking into account its geographical and ecological diversity. Air-water CO₂ fluxes were then up-scaled in global ocean (362×10⁶ km²) using the present estimates for the coastal ocean and those from Takahashi et al. (2002) for the open ocean (336×10⁶ km²). If estuaries and salt marshes are not taken into consideration in the up-scaling, the coastal ocean behaves as a sink for atmospheric CO₂(−1.17 mol C m⁻² yr⁻¹) and the uptake of atmospheric CO₂ by the global ocean increases by 24% (−1.93 versus −1.56 Pg C yr⁻¹). The inclusion of the coastal ocean increases the estimates of CO₂ uptake by the global ocean by 57% for high latitude areas (−0.44 versus −0.28 Pg C yr⁻¹) and by 15% for temperate latitude areas (−2.36 versus −2.06 Pg C yr⁻¹) At subtropical and tropical latitudes, the contribution from the coastal ocean increases the CO₂ emission to the atmosphere from the global oceam by 13% (0.87 versus 0.77 Pg C yr⁻¹). If estuaries and salt marshes are taken into consideration in the upscaling, the coastal ocean behaves as a source for atmospheric CO₂ (0.38 mol C m⁻² yr⁻¹) and the uptake of atmospheric CO₂ from the global ocean decreases by 12% (−1.44 versus −1.56 Pg C yr⁻¹) At high and subtropical and tropical latitudes, the coastal ocean behaves as a source for atmospheric CO₂ but at temperate latitudes, it still behaves as a moderate CO₂ sink. A rigorous up-scaling of air-water CO₂ fluxes in the coastal ocean is hampered by the poorly constrained estimate of the surface area of inner estuaries. The present estimates clearly indicate the significance of this biogeochemically, highly active region of the biosphere in the global CO₂ cycle.     

Possible scenarios for the formation of Ap/Bp stars

Possible paths for the formation of Ap/Bp stars—massive main-sequence stars with strong magnetic fields—are analyzed based on modern theories for the evolution of single and binary stars. Assuming that the strong magnetic fields of these stars are the main reason for their comparatively slow axial rotation and the observed anomalies in the chemical compositions of their atmospheres, possible origins for these high magnetic fields are considered. Analysis of several possible scenarios for the formation of these stars leads to the conclusion that their surface magnetic fields are probably generated in the convective envelopes of the precursor stars. These precursors may be young, single stars with masses 1.5–3 M _⊙ that formed at the peripheries of forming star clusters and ended their accretion at the Hayashi boundary, or alternatively close binaries whose components have convective envelopes, whose merger leads to the formation of an Ap/Bp star. Arguments are presented supporting the view that the merger of close binaries is the main channel for the formation of Ap/Bp stars, and a detailed analysis of this scenario is presented. The initial major axes of the merging binary systems must be in the range 6–12 R _⊙, and the masses of their components in the range 0.7–1.5 M _⊙. When the merging components possess developed convective envelopes and fairly strong initial magnetic fields, these can generate powerful magnetic fields “inherited” by the products of the merger—Ap/Bp stars. The reason the components of the close binaries merge is a loss of angular momentum via the magnetic stellar winds of the components.     

An organic carbon budget for the Mississippi River turbidity plume and plume contributions to air-sea CO2 fluxes and bottom water hypoxia

We investigated seasonal variability in organic carbon (OC) budgets using a physical-biological model for the Mississippi River turbidity plume. Plume volume was calculated from mixed layer depth and area in each of four salinity subregions based on an extensive set of cruise data and satellite-derived suspended sediment distributions. These physical measurements were coupled with an existing food web model to determine seasonally dependent budgets for labile (reactive on time scales of days to weeks) OC in each salinity subregion. Autochthonous gross primary production (GPP) equaled 1.3×10¹² g C yr⁻¹ and dominated labile OC inputs (88% of the budget) because riverine OC was assumed mostly refractory (nonreactive). For perspective, riverine OC inputs amounted to 3.9×10¹² g C yr⁻¹, such that physical inputs were 3 times greater than biological inputs to the plume. Annually, microbial respiration (R) accounted for 65% of labile OC losses and net metabolism (GPP—R) for the entire plume was, autotrophic, equaling 5.1×10¹¹ g C yr⁻¹. Smaller losses of labile OC occurred via sedimentation (20%), advection (10%), and export to higher trophic levels (5%). In our present model, annual losses of labile OC are 10% higher than inputs, indicating future improvements are required. Application of our model to estimate air-sea carbon dioxide (CO₂) fluxes indicated the plume was a net sink of 2.0×10⁹ mol CO₂ yr⁻¹, of which 90% of the total drawdown was from biotic factors. In all seasons, low salinity waters were a source of CO₂ (pCO₂=560–890 μatm), and intermediate to high salinity waters were a sink of CO₂ (pCO₂=200–370 μatm). Our model was also used to calculate O₂ demand for the development, of regional hypoxia, and our spring and early summer budgets indicated that sedimentation of autochthonous OC from the immediate plume contributed 23% of the O₂ demand necessary for establishment of hypoxia in the region.     

Modeling the mass-metallicity relation for dSph dwarf galaxies of the local group

The relationship between the masses and metallicities of galaxies could be non-monotonic, due to the outflow of matter in these systems. It is shown using a simple, one-zone, chemical-dynamicalmodel that the metallicity should be a non-monotonic function of the mass for spheroidal dwarf galaxies with low masses of M ≤ 10⁸ M _⊙, and a monotonically growing function for galaxies with higher masses. This is consistent with observations.     

Thermoelasticity and high-<Emphasis Type="Italic">T</Emphasis> behaviour of anthophyllite

The thermoelastic behaviour of anthophyllite has been determined for a natural crystal with crystal-chemical formula ^ANa_0.01 ^B(Mg_1.30Mn_0.57Ca_0.09Na_0.04) ^C(Mg_4.95Fe_0.02Al_0.03) ^T(Si_8.00)O₂₂ ^W(OH)₂ using single-crystal X-ray diffraction to 973 K. The best model for fitting the thermal expansion data is that of Berman (J Petrol 29:445–522, 1988) in which the coefficient of volume thermal expansion varies linearly with T as α _V,T  = a ₁ + 2a ₂ (T − T ₀): α₂₉₈ = a ₁ = 3.40(6) × 10⁻⁵ K⁻¹, a ₂ = 5.1(1.0) × 10⁻⁹ K⁻². The corresponding axial thermal expansion coefficients for this linear model are: α _a ,₂₉₈ = 1.21(2) × 10⁻⁵ K⁻¹, a _2,a  = 5.2(4) × 10⁻⁹ K⁻²; α _b ,₂₉₈ = 9.2(1) × 10⁻⁶ K⁻¹, a _2,b  = 7(2) × 10⁻¹⁰ K⁻². α _c ,₂₉₈ = 1.26(3) × 10⁻⁵ K⁻¹, a _2,c  = 1.3(6) × 10⁻⁹ K⁻². The thermoelastic behaviour of anthophyllite differs from that of most monoclinic (C2/m) amphiboles: (a) the ε ₁ − ε ₂ plane of the unit-strain ellipsoid, which is normal to b in anthophyllite but usually at a high angle to c in monoclinic amphiboles; (b) the strain components are ε ₁ ≫ ε ₂ > ε ₃ in anthophyllite, but ε ₁ ~ ε ₂ ≫ ε ₃ in monoclinic amphiboles. The strain behaviour of anthophyllite is similar to that of synthetic C2/m ^ANa ^B(LiMg) ^CMg₅ ^TSi₈ O₂₂ ^W(OH)₂, suggesting that high contents of small cations at the B-site may be primarily responsible for the much higher thermal expansion ⊥(100). Refined values for site-scattering at M4 decrease from 31.64 epfu at 298 K to 30.81 epfu at 973 K, which couples with similar increases of those of M1 and M2 sites. These changes in site scattering are interpreted in terms of Mn ↔ Mg exchange involving M1,2 ↔ M4, which was first detected at 673 K.     

Thermal expansion and dehydroxylation of phengite micas

Phengite samples (2M ₁ and 3T politypes) and a synthetic end-member muscovite specimen were studied by in situ high-temperature synchrotron radiation X-ray diffraction. The measured volume thermal expansion of 2M ₁ phengite (<α _V> ≈ 36.6 × 10⁻⁶ K⁻¹) was systematically greater than <α _V> of the 3T polytype (≈33.3 × 10⁻⁶ K⁻¹). A positive linear correlation between the average thermal expansion on (001) plane and the mean tetrahedral rotation angle at ambient condition is proposed on the ground of new measurements and literature data. Dehydroxylation processes were observed in 2M ₁, starting at 1,000 K in 3T at 800 and 945 K in synthetic muscovite. Rietveld refinements allowed a determination of structural variations upon heating of phengite samples and their dehydroxylate phases. The phengite structure expands by regularizing the tetrahedral sheet and by reducing the bond length differences between the outer and inner coordination shell of the interlayer site. The dehydroxylate phase derived from 2M ₁ is characterized by fivefold polyhedra in the low temperature form as a consequence of two OH groups reacting to form H₂O + O (residual). The dehydroxylate exhibits an increase of the cation–cation distances along the M–Or–M bonds with respect to low-temperature phengite structures. For the 3T phase, we were unable to achieve completion of dehydroxylation. The refined structural model of the dehydroxylate phase shows two hydroxyl sites, but at a short distance from one another. This result suggests that the dehydroxylation reaction did not proceed to completion. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Thermal equations of state of dioctahedral micas on the join muscovite–paragonite

 Powder diffraction measurements at simultaneous high pressure and temperature on samples of 2M1 polytype of muscovite (Ms) and paragonite (Pg) were performed at the beamline ID30 of ESRF (Grenoble), using the Paris-Edinburgh cell. The bulk moduli of Ms, calculated from the least-squares fitting of V–P data on each isotherm using a second-order Birch–Murnaghan EoS, were: 57.0(6), 55.1(7), 51.1(7) and 48.9(5) GPa on the isotherms at 298, 573, 723 and 873 K, respectively. The value of (∂K _T /∂T) was −0.0146(2) GPa K⁻¹. The thermal expansion coefficient α varied from 35.7(3) × 10⁻⁶ K⁻¹ at P ambient to 20.1(3) × 10⁻⁶ K⁻¹ at P = 4 GPa [(∂α/∂P) _T = −3.9(1) × 10⁻⁶ GPa⁻¹ K⁻¹]. The corresponding values for Pg on the isotherms at 298, 723 and 823 K were: bulk moduli 59.9(5), 55.7(6) and 53.8(7) GPa, (∂K _T /∂T) −0.0109(1) GPa K⁻¹. The thermal expansion coefficient α varied from 44.1(2) × 10⁻⁶ K⁻¹ at P ambient to 32.5(2) × 10⁻⁶ K⁻¹ at P = 4 GPa [(∂α/∂P) _T = −2.9(1) × 10⁻⁶ GPa⁻¹ K⁻¹]. Thermoelastic coefficients showed that Pg is stiffer than Ms; Ms softens more rapidly than Pg upon heating; thermal expansion is greater and its variation with pressure is smaller in Pg than in Ms. Received: 28 January 2002 / Accepted: 5 April 2002     

IR photometry and models for the dust envelopes of two carbon stars

The results of JHKLM photometry of two carbon stars are presented: the irregular variable NQ Cas and the Mira star BD Vul. Data on the mean fluxes supplemented with mid-IR observations with the IRAS, AKARI, andWISE satellites are used to compute spherically symmetrical model dust envelopes for the stars, consisting of particles of amorphous carbon and silicon carbide. The optical depth in the visible for the comparatively cool dust envelope of BD Vul, with a dust temperature at its inner boundary T₁ = 610 K, is fairly low: τ_V = 0.13. The dust envelope of NQ Cas is appreciably hotter (T₁ = 1550 K), and has τ_V = 0.32. The estimated mass-loss rates are 1.5 × 10^?7M_⊙/yr for NQ Cas and 5.9 × 10^?7M_⊙/yr for BD Vul.     

Evolution of the velocity ellipsoids in the thin disk of the Galaxy and the radial migration of stars

Data from the revised Geneva-Copenhagen catalog are used to study the influence of radial migration of stars on the age dependences of parameters of the velocity ellipsoids for nearby stars in the thin disk of the Galaxy, assuming that the mean radii of the stellar orbits remain constant. It is demonstrated that precisely the radial migration of stars, together with the negative metallicity gradient in the thin disk, are responsible for the observed negative correlation between the metallicities and angular momenta of nearby stars, while the angular momenta of stars that were born at the same Galactocentric distances do not depend on either age or metallicity. The velocity components of the Sun relative to the Local Standard of Rest derived using data for stars born at the solar Galactocentric distance are (U _⊙, V _⊙, W _⊙) _LSR = (5.1 ± 0.4, 7.9 ± 0.5, 7.7 ± 0.2) km/s. The two coordinates of the apex of the solar motion remain equal to 〈l _⊙〉 = 70° ± 7° and 〈b _⊙〉 = 41° ± 2°, within the errors. The indices for the power-law age dependences of them ajor, middle, and minor semi-axes become 0.26±0.04, 0.32±0.03, and 0.07±0.03, respectively. As a result, with age, the velocity ellipsoid for thin-disk stars born at the solar Galactocentric distance increases only in the plane of the disk, remaining virtually constant in the perpendicular direction. Its shape remains far from equilibrium, and the direction of the major axis does not change with age: the ellipsoid vertex deviation remains constant and equal to zero within the errors (〈L〉 = 0.7° ± 0.6°, 〈B〉 = 1.9° ± 1.1°). Such a small increase in the velocity dispersion perpendicular to the Galactic plane with age can probably be explained by “heating” of the stellar system purely by spiral density waves, without a contribution from giant molecular clouds.     

Atmospheric lithium abundances of the carbon giants AQ And,HK Lyr,UX Dra,and WZ Cas

The LTE lithium abundances, logN (Li), have been determined for the atmospheres of the four carbon giants AQ And, HK Lyr, UX Dra, and WZ Cas, applying the synthetic-spectrum technique to the LiI λ670.8 nm resonance doublet. We used estimates of the effective temperatures T _eff, metallicities, elemental abundances, and carbon-isotope ratios for these stars from the literature [12C/13C]. The resulting lithium abundances depend significantly on the adopted T _eff, so that we can find each star’s range of possible logN(Li) values for a given range of effective temperatures (ΔT _eff = 200−300 K). The uncertainty in the derived logN(Li) values is 1.3 for AQ And and UX Dra, 0.5 for HK Lyr, and 0.8 for WZ Cas. Our most probable logN(Li) values (or ranges) are −1.25 for AQ And, −1.2 for HK Lyr, −1.0 ≤ logN(Li) ≤ +0.3 for UX Dra, and 4.3 ≤ logN(Li) ≤ 4.8 for WZ Cas. The derived C/O ratios for the stellar atmospheres also depend on T _eff. Possible uncertainties in other atmospheric parameters do not change logN(Li) by more than 0.2.     

Thermal expansion of gehlenite, Ca2Al[AlSiO7], and the related aluminates LnCaAl[Al2O7] with Ln = Tb, Sm

The thermal expansion of gehlenite, Ca₂Al[AlSiO₇], (up to T=830 K), TbCaAl[Al₂O₇] (up to T=1100 K) and SmCaAl[Al₂O₇] (up to T=1024 K) has been determined. All compounds are of the melilite structure type with space group Thermal expansion data were obtained from in situ X-ray powder diffraction experiments in-house and at HASYLAB at the Deutsches Elektronen Synchrotron (DESY) in Hamburg (Germany). The thermal expansion coefficients for gehlenite were found to be: α₁=7.2(4)×10⁻⁶×K⁻¹+3.6(7)×10⁻⁹ΔT×K⁻² and α₃=15.0(1)×10⁻⁶×K⁻¹. For TbCaAl[Al₂O₇] the respective values are: α₁=7.0(2)×10⁻⁶×K⁻¹+2.0(2)×10⁻⁹ΔT×K⁻² and α₃=8.5(2)×10⁻⁶×K⁻¹+2.0(3)×10⁻⁹ΔT×K⁻², and the thermal expansion coefficients for SmCaAl[Al₂O₇] are: α₁=6.9(2)×10⁻⁶×K⁻¹+1.7(2)×10⁻⁹ΔT×K⁻² and α₃=9.344(5)×10⁻⁶×K⁻¹. The expansion mechanisms of the three compounds are explained in terms of structural trends obtained from Rietveld refinements of the crystal structures of the compounds against the powder diffraction patterns. No structural phase transitions have been observed. While gehlenite behaves like a ‘proper’ layer structure, the aluminates show increased framework structure behavior. This is most probably explained by stronger coulombic interactions between the tetrahedral conformation and the layer-bridging cations due to the coupled substitution (Ca²⁺+Si⁴⁺)–(Ln ³⁺+Al³⁺) in the melilite-type structure. This article has been mistakenly published twice. The first and original version of it is available at .     

Rare-earth elements as source indicators of Pan-African granites from Obudu Plateau,Southeastern Nigeria

The rare-earth element （REE） concentrations of representative granite samples from the southeast of the Obudu Plateau, Nigeria, were analyzed with an attempt to determine the signatures of their source, evolutionary history and tectonic setting. Results indicated that the granites have high absolute REE concentrations （190×10^-6-1191×10^-6; av.=549×10^-6） with the chondrite-normalized REE patterns characterized by steep negative slopes and prominent to slight or no negative Eu anomalies. All the samples are also characterized by high and variable concentrations of the LREE （151×10^-6-1169×10^-6; av.= 466×10^-6）, while the HREE show low abundance （4×10^-6-107×10^-6; av.=28×10^-6）. These are consistent with the variable levels of REE fractionation, and differentiation of the granites. This is further supported by the range of REE contents, the chondrite-normalized patterns and the ratios of LaN/YbN （2.30-343.37）, CeN/YbN （5.94-716.87）, LaN/SmN （3.14-11.68） and TbN/YbN （0.58-1.65）. The general parallelism of the REE patterns, suggest that all the granites were comagmatic in origin, while the high Eu/Eu＊ ratios （0.085-2.807; av.=0.9398） indicate high fo2 at the source. Similarly, irregular variations in LaN/YbN, CeN/YbN and Eu/Eu＊ ratios and REE abundances among the samples suggest behaviors that are related to mantle and crustal sources.     

The evolutionary status of the most massive WNh stars in close binary systems. NGC 3603-A1

The evolution of the components of the unique, massive, close binary system NGC 3603-A1, which consists of stars of spectral types WN6ha and WN6h, is analyzed. The component masses are estimated to be 116 and 89M _⊙, close to the highest measured stellar masses. Numerical modeling of the evolution of the components has been carried out, taking into account mass loss via the stellar winds of the two massive stars. It is shown that the maximum possible initial component masses are close to 140 and 125M _⊙. The components are currently slightly evolved main-sequence stars, with a comparative low degree of helium enrichment at their surfaces. Further evolution of the system will lead to filling of the Roche lobe of the primary and subsequent evolution in a common envelope. This may lead to the merger of the components, with the evolution of the system ending in the formation of a singlemassive black hole after the second supernova explosion. Otherwise, depending on the masses of the resulting black holes, either a binary system of two black holes or two unbound black holes may form, accompanied by gamma-ray bursts.     

Variations in the orbital periods of the Algol-type eclipsing binaries RZ Cas and Z Dra

A detailed study of variations of the orbital periods of the Algol-type eclipsing binary systems RZ Cas and Z Dra is presented. The fairly complex variations of the periods of both systems can be represented as a superposition of a secular increase of the period, slow periodic fluctuations, and quasiperiodic oscillations with a small amplitude occurring on timescales of decades. The secular increase of the period can be explained by the steady mass transfer from the less massive to the more massive component with conservation of the total angular momentum. The mass-transfer rate is 5.7 × 10^?9M_⊙/yr for RZ Cas and 3.0×10^?8M_⊙/yr for Z Dra. To explain the long-period cyclic variations of the orbital periods of RZCas and Z Dra, it must be assumed that the eclipsing binaries move in long-period orbits. RZ Cas moves with a period of 133 yr around a third body withmass M₃ > 0.55M_⊙, while Z Dra moves with a period of 60 yr around a third body with mass M₃ > 0.7M_⊙. The residual fluctuations of the periods may be due to a superposition of variations due to magnetic cycles and non-stationary ejections of matter.     

Thermal expansion of gehlenite, Ca2Al[AlSiO7], and the related aluminates LnCaAl[Al2O7] with Ln=Tb, Sm

The thermal expansion of gehlenite, Ca₂Al[AlSiO₇], (up to T=830 K), TbCaAl[Al₂O₇] (up to T=1,100 K) and SmCaAl[Al₂O₇] (up to T=1,024 K) has been determined. All compounds are of the melilite structure type with space group Thermal expansion data was obtained from in situ X-ray powder diffraction experiments in-house and at HASYLAB at the Deutsches Elektronen Synchrotron (DESY) in Hamburg (Germany). The thermal expansion coefficients for gehlenite were found to be: α₁=7.2(4)×10⁻⁶ K⁻¹+3.6(7)×10⁻⁹ΔT K⁻² and α₃=15.0(1)×10⁻⁶ K⁻¹. For TbCaAl[Al₂O₇] the respective values are: α₁=7.0(2)×10⁻⁶ K⁻¹+2.0(2)×10⁻⁹ΔT K⁻² and α₃=8.5(2)×10⁻⁶ K⁻¹+2.0(3)×10⁻⁹ΔT K⁻², and the thermal expansion coefficients for SmCaAl[Al₂O₇] are: α₁=6.9(2)× 10⁻⁶ K⁻¹+1.7(2)×10⁻⁹ΔT K⁻² and α₃=9.344(5)×10⁻⁶ K⁻¹. The expansion-mechanisms of the three compounds are explained in terms of structural trends obtained from Rietveld refinements of the crystal structures of the compounds against the powder diffraction patterns. No structural phase transitions have been observed. While gehlenite behaves like a ’proper’ layer structure, the aluminates show increased framework structure behaviour. This is most probably explained by stronger coulombic interactions between the tetrahedral conformation and the layer-bridging cations due to the coupled substitution (Ca²⁺+Si⁴⁺)-(Ln ³⁺+Al³⁺) in the melilite-type structure. Electronic Supplementary Material Supplementary material is available for this article at