A Model of the Mira–Type Star T UMi

Stellar evolution calculations were carried out from the main sequence to the final stage of the asymptotic giant branch for stars with initial masses 1 M_⊙ ≤ M_ZAMS ≤ 2 M_⊙ and metallicity Z = 0.01. Selected models of evolutionary sequences were used as initial conditions for solution of the equations of radiation hydrodynamics and time–dependent convection describing radial stellar pulsations. The study was aimed to construct the hydrodynamic models of Mira–type stars that show the secular decrease in the pulsation period Π commenced in 1970th at Π = 315 day. We show that such a condition for the period change is satisfied with evolutionary sequences 1 M_⊙ ≤ M_ZAMS ≤ 1.2 M_⊙ and the best agreement with observations is obtained for M_ZAMS = 1.2 M_⊙. The pulsation period reduction is due to both the stellar radius decrease during the thermal pulse of the helium burning shell and mode switch from the fundamental mode to the first overtone. Theoretical estimates of the fundament parameters of the star at the onset of pulsation period reduction are as follows: the mass is M = 0.93 M_⊙, the luminosity is L = 4080 L_⊙, and the radius is R = 220 R_⊙. The mode switch occurs 35 years after the onset of period reduction.     

Radial pulsations of helium stars with masses from 10 to 50<Emphasis Type="Italic">M</Emphasis><Subscript>⊙</Subscript>

We have performed hydrodynamic calculations of the radial pulsations of helium stars with masses 10M_⊙ ≤ M ≤ 50M_⊙, luminosity-to-mass ratios 5 × 10³L_⊙/M_⊙ ≤ L/M ≤ 2.5 × 10⁴L_⊙/M_⊙, and effective temperatures 2 × 10⁴ K ≤ T_eff ≤ 10⁵ K for helium and heavy-element mass fractions of Y=0.98 and Z=0.02, respectively. We show that the high-temperature boundary of the instability region for radial pulsations at L/M ? 10⁴L_⊙/M_⊙ extends to T_eff≈10⁵ K. The amplitude of the velocity variations for outer layers is several hundred km s^?1, while the brightness variations in the B band of the UBV photometric system are within the range from several hundredths to half a magnitude. At constant luminosity-to-mass ratio, the radial pulsation period is determined only by the effective temperature of the star. In the ranges of luminosity-to-mass ratios 10⁴L_⊙/M_⊙ ≤ L/M ≤ 2 × 10⁴L_⊙/M_⊙ and effective temperatures 5 × 10⁴ K ≤ T_eff ≤ 9 × 10⁴ K, the periods of the radial modes are within 6 min ?Π?10³ min.     

Pulsations of intermediate–mass stars on the asymptotic giant branch

Evolutionary tracks from the zero age main sequence to the asymptotic giant branch were computed for stars with initial masses 2 M _⊙ ≤ M _ZAMS ≤ 5 M _⊙ and metallicity Z = 0.02. Some models of evolutionary sequences were used as initial conditions for equations of radiation hydrodynamics and turbulent convection describing radial stellar pulsations. The early asymptotic giant branch stars are shown to pulsate in the fundamental mode with periods 30 day ? Π ? 400day. The rate of period change gradually increases as the star evolves but is too small to be detected (Π?/Π < 10^?5 yr^?1). Pulsation properties of thermally pulsing AGB stars are investigated on time intervals comprising 17 thermal pulses for evolutionary sequences with initial masses M _ZAMS = 2 M _⊙ and 3 M _⊙ and 6 thermal pulses for M _ZAMS = 4 M _⊙ and 5 M _⊙. Stars with initial masses M _ZAMS ≤ 3 M _⊙ pulsate either in the fundamental mode or in the first overtone, whereas more massive red giants (M _ZAMS ≥ 4 M _⊙) pulsate in the fundamental mode with periods Π ? 10³ day. Most rapid pulsation period change with rate ?0.02 yr^?1 ? Π?/Π ? ?0.01 yr^?1 occurs during decrease of the surface luminosity after the maximum of the luminosity in the helium shell source. The rate of subsequent increase of the period is Π?/Π ? 5 × 10^?3 yr^?1.     

Radial pulsations of helium stars with masses from 1 to 10 <Emphasis Type="Italic">M</Emphasis><Subscript>⊙</Subscript>

We present the results of our hydrodynamic calculations of radial pulsations in helium stars with masses 1 M_⊙ ≤ M ≤ 10 M_⊙, luminosity-to-mass ratios 1 × 10³L_⊙/M_⊙ ≤ L/M ≤ 2 × 10⁴L_⊙/M_⊙, and effective temperatures 2 × 10⁴ K ≤ T_eff ≤ 10⁵ K for mass fractions of helium Y=0.98 and heavy elements Z=0.02. We show that the lower boundary of the pulsation-instability region corresponds to L/M ～ 10³L_⊙/M_⊙ and that the instability region for L/M ? 5 × 10³L_⊙/M_⊙ is bounded by effective temperatures T_eff ? 3 × 10⁴ K. As the luminosity rises, the instability boundary moves into the left part of the Hertzsprung-Russell diagram and radial pulsations can arise in stars with effective temperatures T_eff ? 10⁵ K at L/M ? 7 × 10³L_⊙/M_⊙. The velocity amplitude for the outer boundary of the hydrodynamic model increases with L/M and lies within the range 200 ? ΔU ? 700 km s^?1 for the models under consideration. The periodic shock waves that accompany radial pulsations cause a significant change of the gas-density distribution in the stellar atmosphere, which is described by a dynamic scale height comparable to the stellar radius. The dynamic instability boundary that corresponds to the separation of the outer stellar atmospheric layers at a superparabolic velocity is roughly determined by a luminosity-to-mass ratio L/M ～ 3 × 10⁴L_⊙/M_⊙.     

Evolution and the period-luminosity relation for red supergiants in the Magellanic Clouds

Excitation of radial pulsations in red supergiants of Magellanic Clouds is investigated using the stellar evolution calculations and the self-consistent solution of the equations of radiation hydrodynamics and turbulent convection. The stars with initial masses 6M _⊙ ≤ M _ZAMS ≤ 28M _⊙ and the initial chemical composition X = 0.7, 0.004 ≤ Z ≤ 0.008 are shown to be unstable against fundamental mode oscillations with periods from 17 to 1200 days as they become helium burning red supergiants. The period-luminosity relation slightly depends on the mass loss rate varying with a factor of three, whereas its dependence on the metal abundance is given by δM _bol = 0.89δ log Z. In comparison with galactic red supergiants the low metal abundances in red supergiants of Magellanic Clouds are responsible for their higher effective temperatures and substantially narrower ranges of evolutionary radius change during helium burning. Therefore on the period-mass diagram the red supergiants of Magellanic Clouds are located within the strip with width of δ logM ≈ 0.09, so that the uncertainty of mass evaluation of the red supergiant with the known pulsation period is nearly 25%.     

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.     

Evolution and Period Change in RR Lyr Variables of the Globular Cluster M 3

The grid of evolutionary tracks of population II stars with initial masses 0.81 M_⊙ ≤ M_ZAMS ≤ 0.85 M_⊙ and chemical composition of the globular cluster M3 is computed. Selected models of horizontal branch stars were used as initial conditions for solution of the equations of radiation hydrodynamics and time–dependent convection describing radial stellar oscillations. The boundaries of the instability strip on the Herztsprung–Russel diagram were determined using ≈100 hydrodynamic models of RR Lyr pulsating variables. For each evolutionary track crossing the instability strip the pulsation period was determined as a function of evolutinary time. The rate of period change of most variables is shown to range within ?0.02 ≤ \(\dot{\Pi}\) ≤ 0.05 day/10⁶ yr. Theoretical estimate of the mean period change rate obtained by the population synthesis method is 〈\(\dot{\Pi}\)〉 = 6.0 × 10^?3 day/10⁶ yr and agrees well with observations of RR Lyr variables of the globular cluster M3.     

Instability of radial modes at the Wolf-Rayet evolutionary stage

The evolution of a Population-I star with an initial mass M _ZAMS = 60 M _⊙ has been calculated. At the stage when a red giant turns into an early-type helium star, the vast bulk of the stellar mass is concentrated in a compact core surrounded by an extended envelope that is unstable with respect to radial oscillations. The range of effective temperatures within which the instability arises extends to T _eff ? 10⁵ K. For the models corresponding to the Wolf-Rayet evolutionary stage (5 × 10⁴ K ≤ T _eff ≤ 1.05 × 10⁵ K), hydrodynamic calculations of self-exciting radial stellar pulsations have been performed. The pulsational instability develops in a time interval comparable to the dynamic timescale. Once the amplitude has ceased to grow, the pulsational motions are nonlinear traveling waves propagating from the core boundary to the stellar surface. The velocity amplitude of the outer layers is 500 km s^?1 < ΔU < 10³ km s^?1, depending on the effective temperature. During the evolution of a helium star, the mean ratio of the maximum expansion velocity of the outer layers to the local escape velocity decreases and lies within the range 0.25 < U _max/v _esc < 0.6 for the models considered. The nonlinearity of the stellar pulsations is responsible for the increase in the mean radius \(\bar r\) of the Lagrangian layers compared to the equilibrium radius r _eq. The effect of the increase in mean radius decreases with rising effective temperature from\(\bar r\)/r ～ 10 at T _eff = 7 × 10⁴ K to \(\bar r\)/r ≈ 2 at T _eff = 10⁵ K. The radial pulsation periods for the models considered lie within the range 0.1 day ≤ Π ≤ 1.6 day and the amplitude of the bolometric magnitude variations does not exceed 0 _. ^m 2.     

Theoretical rates of pulsation period change in the Galactic Cepheids

Theoretical estimates of the rates of radial pulsation period change in Galactic Cepheids with initial masses 5.5 M _⊙ ≤ M _ZAMS ≤ 13 M _⊙, chemical composition X = 0.7, Z = 0.02 and periods 1.5 day ≤ Π ≤ 100 day are obtained from consistent stellar evolution and nonlinear stellar pulsation computations. Pulsational instability was investigated for three crossings of the instability strip by the evolutionary track in the HR diagram. The first crossing occurs at the post-main sequence helium core gravitational contraction stage which proceeds in the Kelvin-Helmholtz timescale whereas the second and the third crossings take place at the evolutionary stage of thermonuclear core helium burning. During each crossing of the instability strip the period of radial pulsations is a quadratic function of the stellar evolution time. Theoretical rates of the pulsation period change agree with observations but the scatter of observational estimates of \(\dot \Pi\) noticeably exceeds the width of the band \(\left( {\delta \log \left| {\dot \Pi } \right| \leqslant 0.6} \right)\) confining evolutionary tracks in the period-period change rate diagram. One of the causes of the large scatter with very high values of \(\dot \Pi\) in Cepheids with increasing periods might be the stars that cross the instability strip for the first time. Their fraction ranges from 2% for M _ZAMS = 5.5 M _⊙ to 9% for M _ZAMS = 13 M _⊙ and variables α UMi and IX Cas seem to belong to such objects.     

Rotation curve and mass distribution in the Galaxy from the velocities of objects at distances up to 200 kpc

Three three-component (bulge, disk, halo) model Galactic gravitational potentials differing by the expression for the dark matter halo are considered. The central (bulge) and disk components are described by the Miyamoto–Nagai expressions. The Allen–Santillán (I), Wilkinson–Evans (II), and Navarro–Frenk–White (III) models are used to describe the halo. A set of present-day observational data in the range of Galactocentric distances R from 0 to 200 kpc is used to refine the parameters of thesemodels. For the Allen–Santillán model, a dimensionless coefficient γ has been included as a sought-for parameter for the first time. In the traditional and modified versions, γ = 2.0 and 6.3, respectively. Both versions are considered in this paper. The model rotation curves have been fitted to the observed velocities by taking into account the constraints on the local matter density ρ _⊙ = 0.1 M _⊙ pc^?3 and the force K _z =1.1/2πG = 77 M _⊙ pc^?2 acting perpendicularly to the Galactic plane. The Galactic mass within a sphere of radius 50 kpc, M _G (R ≤ 50 kpc) ≈ (0.41 ± 0.12) × 10¹² M _⊙, is shown to satisfy all three models. The differences between the models become increasingly significant with increasing radius R. In model I, the Galactic mass within a sphere of radius 200 kpc at γ = 2.0 turns out to be greatest among the models considered, M _G (R ≤ 200 kpc) = (1.45 ±0.30)× 10¹² M _⊙, M _G (R ≤ 200 kpc) = (1.29± 0.14)× 10¹² M _⊙ at γ = 6.3, and the smallest value has been found in model II, M _G (R ≤ 200 kpc) = (0.61 ± 0.12) × 10¹² M _⊙. In our view, model III is the best one among those considered, because it ensures the smallest residual between the data and the constructed model rotation curve provided that the constraints on the local parameters hold with a high accuracy. Here, the Galactic mass is M _G (R ≤ 200 kpc) = (0.75 ± 0.19) × 10¹² M _⊙. A comparative analysis with the models by Irrgang et al. (2013), including those using the integration of orbits for the two globular clusters NGC 104 and NGC 1851 as an example, has been performed. The third model is shown to have subjected to a significant improvement.     

Eclipsing binary UU Cas: Radial-velocity curves

This paper reports the results of spectroscopic observations of UUCas obtained with the highresolution (R = 15 000) fiber-fed echelle spectrometer of the 1.2-m telescope of Kourovka Astronomical Observatory of Ural Federal University. The radial velocities of the secondary, more massive and fainter component are measured for the first time. The component mass ratio is found to be q = M ₁/M ₂ = 0.54. The component masses, M ₁ = 9.5M _⊙ and M ₂ = 17.7M _⊙, and the radius of the or bit, A = 52.7R _⊙, are computed for the published orbital inclination of i ~ 69°. Evidence is presented for a disk surrounding the more massive component and a common expanding envelope.     

The most metal-poor quadruple system of subdwarfs G89-14

The system of subdwarfs G89-14 is one of the most metal-poor multiple stars with an atmospheric metal abundance [m/H] = ?1.9. Speckle interferometry at the 6-m BTA telescope has revealed that G89-14 consists of four components. Measurements of the magnitude difference between the components and published data have allowed their masses to be estimated: M _A ≈ 0.67 M _⊙, M _B ≈ 0.24M _⊙,M _C ≈ 0.33M _⊙, andM _D ≈ 0.22M _⊙. The ratio of the orbital periods of the subsystems has been obtained, 0.52 yr: 3000 yr: 650 000 yr (1: 5769: 1 250 000), indicative of a high degree of hierarchy o fG89-14 and its internal dynamical stability. The calculated Galactic orbital elements and the low metallicity of the quadruple system suggest that it belongs to the Galactic halo.     

Metallicity of Young and Old Stars in Irregular Galaxies

Based on archived images obtained with the Hubble Space Telescope, stellar photometry for 105 irregular galaxies has been conducted. We have shown the red supergiant and giant branches in the obtained Hertzsprung-Russel diagrams. Using the TRGB method, distances to galaxies and metallicity of red giants have been determined. The color index (V ? I) of the supergiant branch at the luminosity level M_I = ?7 was chosen as the metallicity index of red supergiants. For the galaxies under study, the diagrams have been built, in which the correlation can be seen between the luminosity of galaxies (M_B) and metallicity of red giants and supergiants. The main source of variance of the results in the obtained diagrams is, in our opinion, uncertainty inmeasurements of galaxy luminosities and star-forming outburst. The relation between metallicity of young and old stars shows that main enrichment of galaxies with metals has taken place in the remote past. Deviations of some galaxies in the obtained relation can possibly be explained with the fall of the intergalactic gas on them, although, this inconsiderably affects metallicities of the stellar content.     

Secular Perion Changes and Fundamental Parameters of Long–Period Cepheids

Hydrodynamic computations of nonlinear Cepheid pulsation models with periods from 20 to 100 day on the evolutionary stage of core helium burning were carried out. Equations of radiation hydrodynamics and time–dependent convection were solved with initial conditions obtained from selected models of evolutionary sequences of population I stars with initial masses from 8 M_⊙ to 12.5 M_⊙. For each crossing of the instability strip the pulsation period Π and the rate of period change \(\dot \prod \) were derived as a function of evolutionary time. Comparing results of our computations with observational estimates of Π and \(\dot \prod \) we determined fundamental parameters (the age, the mass, the luminosity and the radius) of seven long–period Cepheids. Theoretical estimates of the stellar radius are shown to agree with radius measurements by the Baade–Wesselink technique within 3% for RS Pup and GY Sge whereas for SV Vul the disagreement between theory and observations does not exceed 10%.     

A model of low-mass neutron stars with a quark core

We consider an equation of state that leads to a first-order phase transition from the nucleon state to the quark state with a transition parameter λ>3/2 (λ=ρ_Q/(ρ_N+P₀/c²)) in superdense nuclear matter. Our calculations of integrated parameters for superdense stars using this equation of state show that on the stable branch of the dependence of stellar mass on central pressure dM/dP_c>0) in the range of low masses, a new local maximum with M_max=0.082_⊙ and R=1251 km appears after the formation of a toothlike kink (M=0.08M_⊙, R=205 km) attributable to quark production. For such a star, the mass and radius of the quark core are M_core=0.005M_⊙ and R_core=1.73 km, respectively. In the model under consideration, mass accretion can result in two successive transitions to a quark-core neutron star with energy release similar to a supernova explosion: initially, a low-mass star with a quark core is formed; the subsequent accretion leads to configurations with a radius of ～1000 km; and, finally, the second catastrophic restructuring gives rise to a star with a radius of ～100 km.     

Einstein–Maxwell Field Equation in Isotropic Coordinates: An Application to Modeling Superdense Star

We present a charged analogue of Pant et al. (2010, Astrophys. Space Sci., 330, 353) solution of the general relativistic field equations in isotropic coordinates by using simple form of electric intensity E that involve charge parameter K. Our solution is well behaved in all respects for all values of X lying in the range 0 <X≤ 0.11, K lying in the range 4 <K≤ 6.2 and Schwarzschild compactness parameter u lying in the range 0 <u≤ 0.247. Since our solution is well behaved for wide ranges of the parameters, we can model many different types of ultra-cold compact stars like quark stars and neutron stars. We have shown that corresponding to X = 0.077 and K = 6.13 for which u = 0.2051 and by assuming surface density ρ _b =4.6888×10¹⁴ g cm ^?3 the mass and radius are found to be 1.509M _⊙, 10.906 km respectively which match with the observed values of mass 1.51M _⊙ and radius 10.90 km of the quark star XTE J1739-217. The well behaved class of relativistic stellar models obtained in this work might have astrophysical significance in the study of more realistic internal structures of compact stars.     

Spectroscopic studies of yellow supergiants in the Cepheid instability strip

High-resolution spectra of nine yellow nonvariable supergiants (NVSs) located within the canonical Cepheid instability strip from Sandage and Tammann (1969) (α Aqr, ? Leo, μ Per, ω Gem, BD+60 2532, HD 172365, HD 187299, HD 190113, and HD 200102) were taken with the 1-m Zeiss and 6-m BTA telescopes at the Special Astrophysical Observatory of the Russian Academy of Sciences in the 1990s. These have been used to determine the atmospheric parameters, chemical composition, radial velocities, reddenings, luminosities, distances, and radii. The spectroscopic estimates of T _eff and the luminosities determined from the Hipparcos parallaxes have shown eight of the nine program NVSs on the T _eff?log(L/L _⊙) diagram to be outside the canonical Cepheid instability strip. When the edges of the Cepheid instability strip from Bono et al. (2000) are used, out of the NVSs from the list on the diagram one is within the Cepheid instability strip but closer to the red edge, two are at the red edge, three are beyond the red edge, two are at the blue edge, and one is beyond the blue edge. The evolutionary masses of the objects have been estimated. The abundances of α-elements, r- and s-process elements for all program objects have turned out to be nearly solar. The СNO, Na, Mg, and Al abundance estimates have shown that eight of the nine NVSs from the list have already passed the first dredge-up. Judging by the abundances of the key elements and its position on the T _eff?log(L/L _⊙) diagram, the lithium-rich supergiant HD 172365 is at the post-main-sequence evolutionary stage of gravitational helium core contraction and moves toward the first crossing of the Cepheid instability strip. The star ? Leo should be assigned to bright supergiants, while HD 187299 and HD 190113 may have already passed the second dredge-up and move to the asymptotic branch.     

Parameters of the classical type-IIP supernova SN 1999em

Based on observations of SN 1999em, we determined the physical parameters of this supernova using hydrodynamic calculations including nonequilibrium radiative transfer. Taking the distance to SN 1999em estimated by the expanding photosphere method (EPM) to be D = 7.5 Mpc, we found the parameters of the presupernova: radius R = 450R_⊙, mass M = 15M_⊙, and explosion energy E = 7 × 10⁵⁰ erg. For the distance D = 12 Mpc determined from Cepheids, R, M, and E must be increased to the following values: R = 1000R_⊙, M = 18M_⊙, and E = 10⁵¹ erg. We show that one cannot restrict oneself to using the simple analytical formulas relating the supernova and presupernova parameters to obtain reliable parameters for type-IIP presupernovae.     

Infrared photometry for two RV Tau stars and V1027 Cyg

We present and discuss the results of our long-term JHKLM photometry for two RV Tau stars (R Sge and RV Tau) and the yellow supergiant V1027 Cyg, a candidate for protoplanetary nebulae. The amplitude of the infrared brightness variations in R Sge and RV Tau over fourteen years of observations was 0 _. ^m 9?1 ^m ; the infrared brightness variations in V1027 Cyg over eighteen years did not exceed 0 _. ^m 25. The infrared brightness and color of R Sge fluctuated about their gradually changing mean values; the infrared brightness variations agree with a period of 70.77 days. The periodic J brightness and J-H color variations in R Sge can be explained by temperature pulsations with ΔT ≤ 200 K and radial pulsations with [ΔR/R] ≤ 0.2. From 1995 to 2008, the mean J brightness of RV Tau increased, while its mean J-H color index decreased; the variations in the mean J brightness can be associated mainly with stellar temperature variations; a periodic component with P = 78.73 days is observed in the infrared brightness and color fluctuations. The variations in the mean J brightness and J-H color index of the supergiant V1027 Cyg over eighteen years of observations did not exceed a few hundredths of a magnitude; both temperature and radial pulsations may be present in the observed J brightness variations. The most probable period of the infrared brightness fluctuations in V1027 Cyg is 237 ± 2 days. The dust shell of R Sge may consist of two layers with grain temperatures of ～1000 and ～700 K; the optical depth at 1.25 µm is ～0.02 and ～0.24, respectively. The grain temperatures in the circumstellar dust shells of the supergiants RV Tau and V1027 Cyg are ～600 K (RV Tau) and ～700 K (V1027 Cyg). Their optical depths at 1.25 µm are ～0.24 (RV Tau) and ～0.008 (V1027 Cyg).     

Stars with discrepant <Emphasis Type="Italic">v</Emphasis> sin <Emphasis Type="Italic">i</Emphasis> as derived from the Ca II 3933 and Mg II 4481 Å lines. VI. HD 199892 — an SB2 spectroscopic binary

We analyzed the spectra of a well known SB1 binary HD199892 for which the projected rotational velocity v sin i, introduced in the literature, significantly differs when determined from the lines of Ca II at 3933 Å and ofMg II at 4481 Å. Contrary to the former findings, we discovered the signs of spectral lines of a companion star in the profile of Hβ as well as weak metallic lines in the high resolution high S/N spectra covering the most of the visual region. We estimated the secondary star to be a main sequence A4V star with a mass of 2.2M _⊙ and derived its radial velocity which resulted in the mass of the primary M = 4.6M _⊙. Short sections of the spectra in the Mg II 4481 Å and Ca II 3933 Å regions are analyzed as well.