Supersonic steady accretion of a rotating cold “Iron Gas” onto a neutron star after gravitational collapse

We analytically generalize the well-known solution of steady supersonic spherically symmetric gas accretion onto a star (Bondi 1952) for an iron atmosphere with completely degenerate electrons with an arbitrary degree of relativity. This solution is used for typical physical conditions in the vicinity of protoneutron stars produced by gravitational collapse with masses M ₀=(1.4?1.8)M _⊙ and over a wide range of nonzero “iron gas” densities at infinity, ρ_∞=(10⁴?5×10⁶)g cm^?3. Under these conditions, we determine all accretion parameters, including the accretion rate, whose value is ～(10?50)M _⊙s^?1 at M ₀=1.8M _⊙ (it is a factor of 1.7 lower for M ₀=1.4M _⊙, because the accretion rate is exactly ∝M ₀ ² ). We take into account the effect of accreting-gas rotation in a quasi-one-dimensional approximation, which has generally proved to be marginal with respect to the accretion rate.     

Analysis of HST ultraviolet spectra for T Tauri stars: DR Tau

We analyze the spectra of DR Tau in the wavelength range 1200 to 3100 Å obtained with the GHRS and STIS spectrographs from the Hubble Space Telescope. The profiles for the C IV 1550 and He II 1640 emission lines and for the absorption features of some lines indicate that matter falls to the star at a velocity ～300 km s^?1. At the same time, absorption features were detected in the blue wings of the N I, Mg I, Fe II, Mg II, C II, and Si II lines, suggesting mass outflow at a velocity up to 400 km s^?1. The C II, Si II, and Al II intercombination lines exhibit symmetric profiles whose peaks have the same radial velocity as the star. This is also true for the emission features of the Fe II and H₂ lines. We believe that stellar activity is attributable to disk accretion of circumstellar matter, with matter reaching the star mainly through the disk and the boundary layer. At the time of observations, the accretion luminosity was L_ac ? 2L_⊙ at an accretion rate ?10^?7M_⊙ yr^?1. Concurrently, a small (<10%) fraction of matter falls to the star along magnetospheric magnetic field lines from a height ～R_*. Within a region of size ?3.5R_*, the disk atmosphere has a thickness ～0.1R_* and a temperature ?1.5 × 10⁴ K. We assume that disk rotation in this region significantly differs from Keplerian rotation. The molecular hydrogen lines are formed in the disk at a distance <1.4 AU from the star. Accretion is accompanied by mass outflow from the accretion-disk surface. In a region of size <10R_*, the wind gas has a temperature ～7000 K, but at the same time, almost all iron is singly ionized by H I L _α photons from inner disk regions. Where the warm-wind velocity reaches ?400 km s^?1, the gas moves at an angle of no less than 30° to the disk plane. We found no evidence of regions with a temperature above 10⁴ K in the wind and leave open the question of whether there is outflow in the H₂ line formation region. According to our estimate, the star has the following set of parameters: M_* ? 0.9M_⊙, R_* ? 1.8R_⊙, L_* ? 0.9L_⊙, and \(A_V \simeq 0\mathop .\limits^m 9\). The inclination i of the disk axis to the line of sight cannot be very small; however, i≤60°.     

The distribution of interstellar matter in IC 5146

The very young star cluster IC 5146 is studied using star counts, with a view to determining the distribution of interstellar matter in a region where star formation recently occurred. IC 5146 is embedded in a dark nebula which is very dense near its centre. The total mass of interstellar dust in the nebula is found to be about 4.5M _⊙. Comparison of radio and optical observations of the region indicates that gas and dust are not separated to any great degree by radiation from the embedded stars. A gas/dust ratio of about 150/1 by mass is found. This ratio varies with the dust grain model used.     

Kinematic control of the inertiality of the system of Tycho-2 and UCAC2 stellar proper motions

Based on the Ogorodnikov-Milne model, we analyze the proper motions of Tycho-2 and UCAC2 stars. We have established that the model component that describes the rotation of all stars under consideration around the Galactic y axis differs significantly from zero at various magnitudes. We interpret this rotation found using the most distant stars as a residual rotation of the ICRS/Tycho-2 system relative to the inertial reference frame. For the most distant (d≈900 pc) Tycho-2 and UCAC2 stars, the mean rotation around the Galactic y axis has been found to be M ₁₃ ^? =?0.37±0.04 mas yr^?1. The proper motions of UCAC2 stars with magnitudes in the range 12–15^m are shown to be distorted appreciably by the magnitude equation in μ_α cos δ, which has the strongest effect for northern-sky stars with a coefficient of ?0.60±0.05 mas yr^?1 mag^?1. We have detected no significant effect of the magnitude equation in the proper motions of UCAC2 stars brighter than ≈11^m.     

Properties of Wolf-Rayet stars in eclipsing binary systems

The results of investigations of a number of eclipsing Wolf-Rayet binaries are presented. The ‘core’ radiuses, the ‘core’ temperatures and masses of WR stars in the eclipsing WR+OB binary systems V 444 Cyg, CX Cep, CQ Cep, and CV Ser are obtained (see Table I). The results obtained from the light curves analysis of the V 444 Cyg in the range λλ2460 Å-3.5μ give strong evidence for the Beals (1944) model of WR phenomenon. The chromospheric-coronal effects in the WN5 extended atmosphere are not observed up to a distance ofr?20R _⊙. In the Hertzsprung—Russell diagram all the WR stars lie on the left side from the main sequence between the main sequence and the sequence of uniform helium stars (see Figure 9). Their locations are close to those of the helium remnants formed as a result of mass exchange in massive close binary systems. The period variations in the systems V 444 Cyg and CQ Cep have been discovered and a reliable value of the mass loss rateM=10^?5 M _⊙yr^?1 is obtained, for the two WR stars. The results of the photometric and spectroscopic investigations of the WR stars with low mass companions (post X-ray binary stage?) are presented too (see Table II). The masses of the companions are (1–2)M _⊙, their optical luminosity is ～10³⁶, erg s^?1 which implies that these companions cannot be the normal stars. It is possible that these companions are neutron stars accreting from the stellar wind of the WR stars. Low values of the X-ray luminosities of such WR stars with low mass companions imply that the accretion of matter in such systems is distinct from the accretion process in classical X-ray binary systems. It is noted also that the parameters of low massive companions coupled with WR stars are close to those of helium stars.     

Kinematics and parameters of the gas in the vicinity of TW Hya

The following conclusions about the kinematics and parameters of the gas in the vicinity of TW Hya have been drawn from an analysis of optical and ultraviolet line profiles and intensities. The accreting matter rises in the magnetosphere to a distance z>R_* above the disk plane and falls to the star near its equator almost perpendicular to its plane. The matter outflows from a disk region with an outer radius of ≤0.5 AU. The [OI], [SII], and H2 lines originate in the disk atmosphere outside the outflow region, where the turbulent gas velocity is close to the local speed of sound. In the formation region of the forbidden lines, T?8500 K and N_e?5×10⁶ cm^?3, and the hydrogen is almost neutral: x_e<0.03. The absorption features observed in the blue wings of some of the ultraviolet lines originate in the part of the wind that moves almost perpendicular to the disk plane, i.e., in the jet of TW Hya. The V _z gas velocity component in the jet decreases with increasing distance from the jet axis from 200 to 30 km s^?1. The matter outflowing from the inner disk boundary, moves perpendicular to the disk plane in the formation region of blue absorption line components, at a distance of ～0.5 AU from the axis of symmetry of the disk. This region of the wind is collimated into the jet at a distance of <3 AU from the disk plane. The gas temperature in the formation region of absorption components is ?2×10⁴ K, and the gas density is <3×10⁶ cm^?3. This region of the jet is on the order of several AU away from the disk plane, while free recombination in the jet begins even farther from the disk. The mass-loss rate for TW Hya is \(\dot M_w < 7 \times 10^{ - 10} M_ \odot yr^{ - 1}\), which is a factor of 3lower than the mean accretion rate. The relative abundance of silicon and aluminum in the jet gas is at least an order of magnitude lower than its standard value.     

Modeling the luminosity function of galactic low-mass X-ray binaries

The evolution of the family of binaries with a low-mass star and a compact neutron star companion (low-mass X-ray binaries (LMXBs) with neutron stars) ismodeled by the method of population synthesis. Continuous Roche-lobe filling by the optical star in LMXBs is assumed to be maintained by the removal of orbital angular momentum from the binary by a magnetic stellar wind from the optical star and the radiation of gravitational waves by the binary. The developed model of LMXB evolution has the following significant distinctions: (1) allowance for the effect of the rotational evolution of a magnetized compact remnant on themass transfer scenario in the binary, (2) amore accurate allowance for the response of the donor star to mass loss at the Roche-lobe filling stage. The results of theoretical calculations are shown to be in good agreement with the observed orbital period-X-ray luminosity diagrams for persistent Galactic LMXBs and their X-ray luminosity function. This suggests that the main elements of binary evolution, on the whole, are correctly reflected in the developed code. It is shown that most of the Galactic bulge LMXBs at luminosities L _x > 10³⁷ erg s^?1 should have a post-main-sequence Roche-lobe-filling secondary component (low-mass giants). Almost all of the models considered predict a deficit of LMXBs at X-ray luminosities near ～10^36.5 erg s^?1 due to the transition of the binary from the regime of angular momentum removal by a magnetic stellar wind to the regime of gravitational waves (analogous to the widely known period gap in cataclysmic variables, accreting white dwarfs). At low luminosities, the shape of the model luminosity function for LMXBs is affected significantly by their transient behavior-the accretion rate onto the compact companion is not always equal to the mass transfer rate due to instabilities in the accretion disk around the compact object. The best agreement with observed binaries is achieved in the models suggesting that heavy neutron stars with masses 1.4–1.9M _⊙ can be born.     

Low‐temperature primordial gas in merging halos

The thermal regime of the baryons behind shock waves arising in the process of virialization of dark matter halos is governed at certain conditions by radiation of HD lines. A small fraction of the shocked gas can cool down to the temperature of the cosmic microwave background (CMB). We estimate an upper limit for this fraction: at z = 10 it increases sharply from about q_T ∼ 10^–3 for dark halos of M = 5 × 10⁷ M_⊙ to ∼ 0.1 for halos with M = 10⁸ M_⊙. Further increase of the halo mass does not lead however to a significant growth of q_T – the asymptotic value for M ≫ 10⁸ M_⊙ is 0.3. We estimate the star formation rate associated with such shock waves, and show that they can provide a small but not negligible fraction of the star formation. We argue that extremely metal‐poor low‐mass stars in the Milky Way may have been formed from primordial gas behind such shocks. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The lithium-rotation correlation for WTTS in Taurus–Auriga

Surface lithium abundance and rotation velocity can serve as powerful and mutually complementary diagnostics of interior structure of stars. So far, the processes responsible for the lithium depletion during pre-main sequence evolution are still poorly understood. We investigate whether a correlation exists between equivalent widths of Li (EW(Li)) and rotation period (P_rot) for weak-line T Tauri stars (WTTSs). We find that rapidly rotating stars have lower EW(Li) and the fast burning of Li begins at the phase when star’s P_rot evolves towards 3 days among 0.9M_⊙ to 1.4M_⊙ WTTSs in Taurus–Auriga. Our results support the conclusion by Piau and Turch-Chiéze about a model for lithium depletion with age of the star and by Bouvier et al. in relation to rotation evolution. The turn over of the curve for the correlation between EW(Li) and P_rot is at the phase of zero-age main sequence (ZAMS). The EW(Li) decreases with decreasing P_rot before the star reaches the ZAMS, while it decreases with increasing P_rot (decreasing rotation velocity) for young low-mass main sequence stars. This result could be explained as an age effect of Li depletion and the rapid rotation does not inhibit Li destruction among low-mass PMS stars.     

Kinematic parameters of the galactic spiral pattern from data on open star clusters and OB stars

Data on the positions, radial velocities, and proper motions of open star clusters and OB stars are used to obtain the rotation curve of the Galaxy fitted by a polynomial in inverse powers of the distances from the Galactic rotation axis. We determine the locations of the corotation region and the inner and outer Lindblad resonances using a previously estimated pattern speed. Based on data for objects of the Carina-Sagittarius and Orion arms, we have determined the distortion amplitudes of the velocity field of the Galactic disk, ?_R = ?3.97±4.79 km s^?1 and f_θ=+13.27±2.57 km s^?1.     

Dust in the disk winds from young stars as a source of the circumstellar extinction

We consider the problem of dust grain survival in the disk winds from T Tauri and Herbig Ae stars. For our analysis, we have chosen a disk wind model in which the gas component of the wind is heated through ambipolar diffusion to a temperature of ～10⁴ K. We show that the heating of dust grains through their collisions with gas atoms is inefficient compared to their heating by stellar radiation and, hence, the grains survive even in the hot wind component. As a result, the disk wind can be opaque to the ultraviolet and optical stellar radiation and is capable of absorbing an appreciable fraction of it. Calculations show that the fraction of the wind-absorbed radiation for T Tauri stars can be from 20 to 40% of the total stellar luminosity at an accretion rate ? _a = 10^?8-10^?6 M _⊙yr^?1. This means that the disk winds from T Tauri stars can play the same role as the puffed-up inner rim in current accretion disk models. In Herbig Ae stars, the inner layers of the disk wind (r ≤ 0.5 AU) are dust-free, since the dust in this region sublimates under the effect of stellar radiation. Therefore, the fraction of the radiation absorbed by the disk wind in this case is considerably smaller and can be comparable to the effect from the puffed-up inner rim only at an accretion rate of the order of or higher than 10^?6 M _⊙yr^?1. Since the disk wind is structurally inhomogeneous, its optical depth toward the observer can be variable, which should be reflected in the photometric activity of young stars. For the same reason, moving shadows from gas and dust streams with a spiral-like shape can be observed in high-angular-resolution circumstellar disk images.     

Detection of accreting circumstellar gas around weak emission-line Herbig Ae/Be stars

We present archival and recent IUE high dispersion spectra of late B stars which reveal the presence of accreting gas with velocities as high as 350 km s^–1, collisional ionization of the accreting gas to temperatures above the stellar T_eff, and column densities intermediate between those observed toward classical Herbig Ae/Be stars and the nearby proto-planetary system Pictoris. One of the stars, HD 176386, while lacking obvious optical signatures of youth, is a member of the R CrA star formation region, and with an inferred age of 2.8 Myr has not yet arrived on the zero-age main sequence (ZAMS). The other object, an isolated, field B star with pronounced IR excess due to warm, circumstellar dust, 51 Oph, exhibits only modest h emission. The combination of high velocity, accreting gas in systems with IR excesses due to circumstellar dust suggests that not only are these objects candidate proto-planetary systems, but that they may represent an extension to higher stellar masses of the weak-emission pre-main sequence (PMS) stars.Paper presented at the Conference onPlanetary Systems: Formation, Evolution, and Detection held 7–10 December, 1992 at CalTech, Pasadena, California, U.S.A.     

Non-LTE modeling of narrow emission components of He and Ca lines in optical spectra of classical T Tauri stars

Using LTE calculations of the structure of T Tauri stellar atmospheres heated by radiation from an accretion shock (Dodin and Lamzin 2012), we have calculated the spectrum of the hot spot emerging on the stellar surface by taking into account non-LTE effects for He I, He II, Ca I, and Ca II. Assuming the pre-shock gas density N ₀ and velocity V ₀ to be the same at all points of the accretion stream cross section, we have calculated the spectrum of the star+circular spot system at various N ₀, V ₀, and parameters characterizing the star and the spot. Using nine stars as an example, we show that the theoretical optical spectra reproduce well the observed veiling of photospheric absorption lines as well as the profiles and intensities of the so-called narrow components of He II and Ca I emission lines with an appropriate choice of parameters. The accreted gas density in all of the investigated stars except DK Tau has been found to be N ₀ > 10¹² cm^?3. We have managed to choose the parameters for eight stars at a calcium abundance in the accreted gas ξ _Ca equal to the solar one, but we have been able to achieve agreement between the calculations and observations for TW Hya only by assuming ξ _Ca to be approximately a factor of 3 lower than the solar one. The estimated parameters do not depend on interstellar extinction, because they have been determined from the spectra normalized to the continuum level. The calculated intensity of Ca II lines has turned out to be lower than the observed one, but this contradiction can be eliminated by assuming that, in addition to the accreted gas with a high density N ₀, a more rarefied gas also falls onto the star. The theoretical equivalent widths and relative intensities of the subordinate He I lines disagree significantly with the observations. This is apparently because non-LTE effects should be taken into account when calculating the structure of the upper layers of the hot spot, the accuracy of the cross sections for collisional processes from upper levels is insufficient, and the spot inhomogeneity should probably be taken into account.     

Kinematic control of the inertiality of ICRS catalogs

We perform a kinematic analysis of the Hipparcos and TRC proper motions of stars by using a linear Ogorodnikov-Milne model. All of the distant (r>0.2 kpc) stars of the Hipparcos catalog have been found to rotate around the Galactic y axis with an angular velocity of M ₁₃ ^? =?0.36±0.09 mas yr^?1. One of the causes of this rotation may be an uncertainty in the lunisolar precession constant adopted when constructing the ICRS. In this case, the correction to the IAU (1976) lunisolar precession constant in longitude is shown to be Δp₁=?3.26±0.10 mas yr^?1. Based on the TRC catalog, we have determined the mean Oort constants: A=14.9±1.0 and B=?10.8±0.3 km s^?1 kpc^?1. The component of the model that describes the rotation of all TRC stars around the Galactic y axis is nonzero for all magnitudes, M ₁₃ ^? =?0.86±0.11 mas yr^?1.     

Galactic kinematics from a sample of young massive stars

Based on published sources, we have created a kinematic database on 220 massive (> 10 M _⊙) young Galactic star systems located within ≤3 kpc of the Sun. Out of them, ≈100 objects are spectroscopic binary and multiple star systems whose components are massive OB stars; the remaining objects are massive Hipparcos B stars with parallax errors of no more than 10%. Based on the entire sample, we have constructed the Galactic rotation curve, determined the circular rotation velocity of the solar neighborhood around the Galactic center at R ₀ = 8kpc, V ₀ = 259±16 km s^?1, and obtained the following spiral density wave parameters: the amplitudes of the radial and azimuthal velocity perturbations f _R = ?10.8 ± 1.2 km s^?1 and f _θ = 7.9 ± 1.3 km s^?1, respectively; the pitch angle for a two-armed spiral pattern i = ?6.0° ± 0.4°, with the wavelength of the spiral density wave near the Sun being λ = 2.6 ± 0.2 kpc; and the radial phase of the Sun in χ _⊙ = ?120° ± 4°. We show that such peculiarities of the Gould Belt as the local expansion of the system, the velocity ellipsoid vertex deviation, and the significant additional rotation can be explained in terms of the density wave theory. All these effects decrease noticeably once the influence of the spiral density wave on the velocities of nearby stars has been taken into account. The influence of Gould Belt stars on the Galactic parameter estimates has also been revealed. Eliminating them from the kinematic equations has led to the following new values of the spiral density wave parameters: f _θ = 2.9 ± 2.1 km s^?1 and χ _⊙ = ?104° ± 6°.     

Relation between the X-ray and optical luminosities in binary systems with accreting nonmagnetic white dwarfs

We investigate the relation between the optical (g-band) and X-ray (0.5–10 keV) luminosities of accreting nonmagnetic white dwarfs. According to the present-day counts of the populations of star systems in our Galaxy, these systems have the highest space density among the close binary systems with white dwarfs. We show that the dependence of the optical luminosity of accreting white dwarfs on their X-ray luminosity forms a fairly narrow one-parameter curve. The typical half-width of this curve does not exceed 0.2–0.3 dex in optical and X-ray luminosities, which is essentially consistent with the amplitude of the aperiodic flux variability for these objects. At X-ray luminosities L _x ～ 10³² erg s^?1 or lower, the optical g-band luminosity of the accretion flow is shown to be related to its X-ray luminosity by a factor ～2–3. At even lower X-ray luminosities (L _x ? 10³⁰ erg s^?1), the contribution from the photosphere of the white dwarf begins to dominate in the optical spectrum of the binary system and its optical brightness does not drop below M _g ～ 13–14. Using the latter fact, we show that in current and planned X-ray sky surveys, the family of accreting nonmagnetic white dwarfs can be completely identified to the distance determined by the sensitivity of an optical sky survey in this region. For the Sloan Digital Sky Survey (SDSS) with a limiting sensitivity m _g ～ 22.5, this distance is ～400–600 pc.     

Kinematics of gas and stars in circumnuclear star-forming regions of early-type spirals

We present high-resolution (R～20,000) spectra in the blue and the far red of circumnuclear star-forming regions (CNSFRs) in three early-type spirals (NGC3351, NGC2903 and NGC3310), which have allowed the study of the kinematics of the stars and the ionized gas in these structures and, for the first time, the derivation of their dynamical masses for the first two. In some cases, these regions, about 100 to 150 pc in size, are composed of several individual star clusters with sizes between 1.5 and 4.9 pc, estimated from Hubble Space Telescope images. The stellar dispersions have been obtained from the Calcium triplet (CaT) lines at λ λ 8494, 8542, 8662 Å, while the gas velocity dispersions have been measured by means of Gaussian fits to the Hβ and [Oiii]λ 5007 Å lines in the high-dispersion spectra. Values of the stellar velocity dispersions are between 30 and 68 km?s^?1. We apply the virial theorem to estimate dynamical masses of the clusters, assuming that systems are gravitationally bounded and spherically symmetric, and using previously measured sizes. The measured values of the stellar velocity dispersions yield dynamical masses of the order of 10⁷ to 10⁸ M_⊙ for the entire CNSFRs. Stellar and gas velocity dispersions are found to differ by about 20 to 30 km?s^?1, with the Hβ emission lines being narrower than both the stellar lines and the [Oiii]λ 5007 Å lines. The douby-ionized oxygen, on the other hand, exhibits velocity dispersions comparable to those of the stars or, in some cases, even larger. We have found indications of the presence of two different kinematical components in the ionized gas of the regions. We have mapped the velocity field in the central kpc of the spiral galaxies NGC3351 and NGC2903. For the first object, the radial velocity curve shows deviations from circular motions for the ionized hydrogen consistent with its infall towards the central regions of the galaxy, at a velocity of about 25 km?s^?1. For NGC3310, we present preliminary results for the velocity dispersions for one of the two observed slit position angles, two CNSFRs and the nucleus.     

Evolution of Intergalactic Gas in the Neighborhood of Dwarf Galaxies and Its Manifestations in the HI 21 cm Line

Low-mass galaxies are known to have played the crucial role in the hydrogen reionization in the Universe. In this paper we investigate the contribution of soft x-ray radiation (E ~ 0.1–1 keV) from dwarf galaxies to hydrogen ionization during the initial reionization stages. The only possible sources of this radiation in the process of star formation in dwarf galaxies during the epochs preceding the hydrogen reionization epoch are hot intermediate-mass stars (M ~ 5–8 M_⊙) that entered the asymptotic giant branch (AGB) stage and massive x-ray binaries. We analyze the evolution of the intergalactic gas in the neighborhood of a dwarf galaxy with a total mass of 6 × 10⁸M_⊙ formed at the redshift of z ~ 15 and having constant star-formation rate of 0.01–0.1 M_⊙ yr^?1 over a starburst with a duration of up to 100 Myr. We show that the radiation from AGB stars heats intergalactic gas to above 100 K and ensures its ionization x_e ? 0.03 within about 4–10 kpc from the galaxy in the case of a star-formation rate of star formation 0.03–0.1 M_⊙ yr^?1, and that after the end of the starburst this region remains quasi-stationary over the following 200–300 Myr, i.e., until z ~ 7.5. Formation of x-ray binaries form in dwarf galaxies at z ~ 15 results in a 2–3 and 5–6 times greater size of the ionized and heated region compared to the case where ionization is produced by AGB stars exclusively, if computed with the “x-ray luminosity–star-formation rate” dependence (L_X ~ f_XSFR) factor f_X = 0.1 and f_X ~ 1, respectively. For f_X ? 0.03 the effect of x-ray binaries is smaller that that of AGB star population. Lyα emission, heating, and ionization of the intergalactic gas in the neighborhood of dwarf galaxies result in the excitation of the 21 cm HI line. We found that during the period of the starburst end at z ~11.5–12.5 the brightness temperature in the neighborhood of galaxies is 15–25 mK and the region where the brightness temperature remains close to its maximum has a size of about 12–30 kpc. Hence the epoch of the starburst end is most favorable for 21 cm HI line observations of dwarf galaxies, because at that time the size of the region of maximum brightness temperature is the greatest over the entire evolution of the dwarf galaxy. In the case of the sizes corresponding to almost 0.’1 for z ~ 12 regions with maximum emission can be detected with the Square Kilometre Array, which is currently under construction.     

Radial pulsations of red giant branch stars

We performed hydrodynamic computations of nonlinear stellar pulsations of population I stars at the evolutionary stages of the ascending red giant branch and the following luminosity drop due to the core helium flash. Red giants populating this region of the Hertzsprung–Russel diagram were found to be the fundamental mode pulsators. The pulsation period is the largest at the tip of the red giant branch and for stars with initial masses from 1.1 M _⊙ to 1.9 M _⊙ ranges from ∏ ≈ 254 day to ∏ ≈ 33 day , respectively. The rate of period change during the core helium flash is comparable with rates of secular period change in Mira type variables during the thermal pulse in the helium shell source. The period change rate is largest (∏?/∏ ≈ ?10^?2 yr^?1) in stars with initial mass M ^ZAMS = 1.1 M ^⊙ and decreases to ∏?/∏ ～ ?10^?3 yr^?1 for stars of the evolutionary sequence M ^ZAMS = 1.9 M ^⊙. Theoretical light curves of red giants pulsating with periods ∏ > 200 day show the presence of the secondary maximum similar to that observed in many Miras.