CCD photometry of RV rau stars. VII. EP lyrae

Photometric observations of the variable star EP Lyr were performed with a CCD photometer during the observing season of 2002. Analysis of these observations together with published data has confirmed the mean period of the main variability cycle P = 83.^d248 over almost 100 years. The periodicity of the variations in the main cycle is investigated on the basis of O-C diagrams. The time scale of its variations ranges from 1–2 to 8–20 thousand days.     

Total sunspot area before 1875 (reconstruction)

Empirical functions approximating the dependences of total sunspot area A on relative sunspot number W and group sunspot number GN have been found. In the function A(W), allowance for its dependence on the secular activity cycle has been made; it is shown that this allowance is not needed for the function A(GN). The yearly mean A for 1700–1874 have been reconstructed using these functions and the available W and GN time series. Having supplemented the original data with archival observations, we have been able to reconstruct the monthly mean A _W since 1821. We discuss the causes of the systematic difference between the reconstructions using W and GN.     

Comparison of the X-ray and optical variabilities in the Seyfert galaxy 3C 120

Based on our UBV RI observations and X-ray data from the RXTE satellite, we have investigated the variability of the galaxy 3C 120 over the period 1996–2008. The relative variability amplitude in the U and B bands without any subtraction of the contribution from the underlying galaxy is 23 and 22%, respectively, against 21% in the X-ray band. The autocorrelation function based on the B-band data is considerably wider than that based on the X-ray data. The structure functions on a time scale from 1 to ～100–300 days in the X-ray and optical spectral ranges have the form of a power law (SF ～ τ ^b ). However, their indices differ significantly: b = 0.42 in the X-ray band and b = 1.36 in the B band. Considering the X-ray and optical variabilities as a superposition of independent flares in a wide range of durations, we may conclude that the amplitudes of short flares in the X-ray band are higher than those in the optical one and, conversely, the relative amplitudes of long flares in the X-ray band are slightly lower than those in the optical one, i.e., short events dominate in the X-ray band. The optical flux variations in the R _c and I _c bands lag significantly behind those in the B band, by 3.9 _?0.7 ^+1.0 and 6.2 _?0.6 ^+1.1 days, respectively, if the lag is estimated from the centroid of the cross-correlation function. The X-ray variability on a time scale of about 1800 days (～5 yr) lags behind the B-band variations by 5.3 _?3.3 ^+2.7 days, but the confidence level of this estimate is only 87%. A more detailed analysis of the correlation between the X-ray and optical emissions has revealed a fairly complex picture: different degrees of correlation between the optical and X-ray fluxes are observed at different times.     

Optical variability of the active galactic nucleus 1E 0754.6+3928 and reverberation-based mass estimate for the central black hole

We analyze the spectroscopic and photometric observations of the active galactic nucleus (AGN) 1E 0754.6+3928 performed at Crimean Astrophysical Observatory in 1998–2004. Based on formal spectroscopic criteria for the optical wavelength range, we can classify this object as a narrow-line Seyfert 1 (NLS1) galaxy. Over the period of its observations, the AGN 1E 0754.6+3928 showed a very low rms flux variability amplitude in both continuum and Hβ (～3%). The Hβ time lag relative to the continuum has been found from the cross-correlation function centroid to be τ _cent = 112 _?67 ⁺²¹⁵ days. This result and the Hβ line dispersion estimate have allowed us to determine the central black hole mass by the reverberation method: 1.05 × 10⁸ M _⊙. The position of 1E 0754.6+3928 in the black hole mass-luminosity diagram agrees with the positions of other NLS1 galaxies and can be explained by an enhanced mass accretion rate in the central source.     

Photometric and spectroscopic study of the supergiant with an infrared excess V1027 Cygni

We present the results of our UBV and JHKLM photometry for the semiregular pulsating variable V1027 Cyg, a supergiant with an infrared excess, over the period from 1997 to 2015 (UBV) and in 2009–2015 (JHKLM). Together with the new data, we analyze the photometric observations of V1027 Cyg that we have obtained and published previously. Our search for a periodicity in the UBV brightness variations has led to several periods from P = 212^d to 320^d in different time intervals. We have found the period P = 237^d based on our infrared photometry. The variability amplitude, the lightcurve shape, and themagnitude of V1027 Cyg atmaximum light change noticeably from cycle to cycle. The deepest minimum was observed in 2011, when the amplitudes of brightness variations in the star reached the following values: ΔU = 1 _. ^m 28, ΔB = 1 _. ^m 10, ΔV = 1 _. ^m 05, ΔJ = 0 _. ^m 30, ΔH = 0 _. ^m 35, ΔK = 0 _. ^m 32, ΔL = 0 _. ^m 26, and ΔM = 0 _. ^m 10. An ambiguous correlation of the B ? V and U ? B colors with the brightness has been revealed. For example, a noticeable bluing of the star was observed during the deep 1992, 2008, and 2011 minima, while the variations with smaller amplitudes show an increase in B ? V at the photometric minima. The spectral energy distribution for V1027 Cyg from our photometry in the range 0.36 (U)–5.0 (M) μm corresponds to spectral types from G8I to K3I at different phases of the pulsation cycle. Low-resolution spectra of V1027 Cyg in the range λ4400–9200 ?A were taken during 16 nights over the period 1995–2015. At the 1995 and 2011 photometric minima the star’s spectrum exhibited molecular TiO bands whose intensity corresponded to spectral types M0–M1, while the photometric data point to a considerably earlier spectral type. We hypothesize that the TiO bands are formed in the upper layers of the extended stellar atmosphere. We have measured the equivalent widths of the strongest absorption lines, in particular, the infrared Ca II triplet in the spectrum of V1027 Cyg. The calcium triplet (Ca T) with W _λ(Ca T) = 20.3 ± 1.8 ?A as a luminosity indicator for supergiants places V1027 Cyg in the region of the brightest G–K supergiants. V1027 Cyg has been identified with the infrared source IRAS 20004+2955 and is currently believed to be a candidate for post-AGB stars. The evolutionary status of the star and its difference from other post-AGB objects are discussed.     

3D stellar reddening map from 2MASS photometry: An improved version

An improved version of the 3D stellar reddening map in a space with a radius of 1200 pc around the Sun and within 600 pc of the Galactic midplane is presented. As in the previous 2010 and 2012 versions of the map, photometry with an accuracy better than 0.05 ^m in the J and Ks bands for more than 70 million stars from the 2MASS catalogue is used in the new version. However, the data reduction technique is considerably more complicated. As before, an analysis of the distribution of stars near the main-sequence turnoff on the (J ? Ks)?Ks diagram, where they form a distribution maximum, provides a basis for the method. The shift of this maximum, i.e., the mode (J ? Ks), along (J ? Ks) and Ks, given the spatial variations of the mean dereddened color (J ? Ks)₀ of these stars, is interpreted as a growth of the reddening with increasing distance. The main distinction of the new method is that instead of the fixed mean absolute magnitude, dereddened color, distance, and reddening for each cell, the individual values of these quantities are calculated for each star by iterations when solving the system of equations relating them. This has allowed one to increase the random accuracy of the map to 0.01 ^m and its spatial resolution to 20 pc in coordinates and distance and to 1° in longitude and latitude. Comparison with other reddening estimates for the same spatial cells and Gaia DR1 TGAS stars shows that the constructed map is one of the best maps for the space under consideration. Its systematic errors have been estimated to be σ(E(J ? Ks)) = 0.025 ^m , or σ(E(B ? V)) = 0.04 ^m . The main purpose of the map is to analyze the characteristics of Galactic structures, clouds, and cloud complexes. For this purpose, the reddening map within each spatial cell has also been computed by analyzing the reddening along each line of sight.     

A model of possible variations of the galactic cosmic ray intensity over the recent billion years

Based on the analysis of published data on exposure ages of iron meteorites determined with the ⁴⁰K/K method (T _K) and ages calculated using short-lived cosmogenic radionuclides (with the half-life T _1/2 < 1 Myr) in combination with stable cosmogenic isotopes of noble gases (T_RS), the following results have been obtained. (1) The distribution of T _RS ages (106 values) has an exponential shape, similar to that for ordinary chondrites, but different from the distribution of T _K ages (80 values). The difference is most likely due to small amounts of data for meteorites with low T _K ages (less than ~200–300 Myr). The latter can be ascribed to the difficulty of measurement of small concentrations of cosmogenic potassium isotopes. This circumstance makes the selection of meteorites with ⁴⁰K/K ages nonrepresentative and casts doubt on the correctness of conclusions about the variations of the intensity of galactic cosmic rays (GCR) based on the analysis of distribution of these ages. (2) The magnitude of the known effect (systematic overestimation of T _K ages in comparison with T _RS ages) has been refined. The value k = T _K/T _RS = 1.51 ± 0.03 is acquired for the whole population of data. We have shown the inefficiency of the explanation of this effect on account of an exponential change in the GCR intensity (I _T ) with time (T) according to the relation I _T = I ₀exp(–γT) over the whole range of ages of iron meteorites. (3) In order to explain the overestimation of T _K ages in comparison with T _RS ages, a model has been proposed, according to which the GCR intensity has exponentially increased in the interval of 0–1500 Myr governed by the relation: I _T = I _{T = 1500} (1 + αexp(–βT)). For one of the variants of this model, the GCR intensity has exponentially increased by a factor of two only over the recent ~300 Myr, remaining approximately constant for the rest of the time. The data acquired with the use of this model indicate that the measured T _K ages are close to the actual time that the meteorites existed in space; the data are in agreement with the observed exponential distribution of T _RS ages.     

Kinematics of B-F Stars as a Function of Their Dereddened Color from Gaia and PCRV Data

Parallaxes with an accuracy better than 10% and proper motions from the Gaia DR1 TGAS catalogue, radial velocities from the Pulkovo Compilation of Radial Velocities (PCRV), accurate Tycho-2 photometry, theoretical PARSEC, MIST, YaPSI, BaSTI isochrones, and the most accurate reddening and interstellar extinction estimates have been used to analyze the kinematics of 9543 thin-disk B-F stars as a function of their dereddened color. The stars under consideration are located on the Hertzsprung–Russell diagram relative to the isochrones with an accuracy of a few hundredths of a magnitude, i.e., at the level of uncertainty in the parallax, photometry, reddening, extinction, and the isochrones themselves. This has allowed us to choose the most plausible reddening and extinction estimates and to conclude that the reddening and extinction were significantly underestimated in some kinematic studies of other authors. Owing to the higher accuracy of TGAS parallaxes than that of Hipparcos ones, the median accuracy of the velocity components U, V, W in this study has improved to 1.7 km s^?1, although outside the range ?0.1 ^m < (B _T ? V _T )₀ < 0.5 ^m the kinematic characteristics are noticeably biased due to the incompleteness of the sample. We have confirmed the variations in the mean velocity of stars relative to the Sun and the stellar velocity dispersion as a function of their dereddened color known from the Hipparcos data. Given the age estimates for the stars under consideration from the TRILEGAL model and the Geneva–Copenhagen survey, these variations may be considered as variations as a function of the stellar age. A comparison of our results with the results of other studies of the stellar kinematics near the Sun has shown that selection and reddening underestimation explain almost completely the discrepancies between the results. The dispersions and mean velocities from the results of reliable studies fit into a ±2 km s^?1 corridor, while the ratios σ _V /σ _U and σ _W /σ _U fit into ±0.05. Based on all reliable studies in the range ?0.1 ^m < (B _T ? V _T )₀ < 0.5^m, i.e., for an age from 0.23 to 2.4 Gyr, we have found: W_⊙ = 7.15 km s^?1, \({\sigma _U} = 16.0{e^{1.29({B_T} - {V_T})o}}\), \({\sigma _V} = 10.9{e^{1.11({B_T} - {V_T})o}}\), \({\sigma _W} = 6.8{e^{1.46({B_T} - {V_T})o}}\), the stellar velocity dispersions in km s^?1 are proportional to the age in Gyr raised to the power β _U = 0.33, β _V = 0.285, and β _W = 0.37.     

Analysis of the Z distribution of young objects in the Galactic thin disk

We have obtained new estimates of the Sun’s distance from the symmetry plane Z _⊙ and the vertical disk scale height h using currently available data on stellar OB associations, Wolf–Rayet stars, HII regions, and Cepheids. Based on individual determinations, we have calculated the mean Z _⊙ = ?16 ± 2 pc. Based on the model of a self-gravitating isothermal disk for the density distribution, we have found the following vertical disk scale heights: h = 40.2 ± 2.1 pc from OB associations, h = 47.8 ± 3.9 pc from Wolf–Rayet stars, h = 48.4 ± 2.5 pc from HII regions, and h = 66.2 ± 1.6 pc from Cepheids. We have estimated the surface, Σ = 6 kpc^?2, and volume, D(Z _⊙) = 50.6 kpc^?3, densities from a sample of OB associations. We have found that there could be ～5000 OB associations in the Galaxy.     

FG Sge at the stage of dust shell ejection

New photometric observations of the variable star FG Sge, a rapidly evolving planetary nebula nucleus, were performed in 2003–2008. On 230 nights, we obtained 86 UBV and 155 BV RI (or R _c , I _c ) magnitude estimates. The maximum amplitude of the V-band light variations was >8 ^m . Six deep minima and four high maxima were observed. Analysis of the light curve has shown that the pulsation period of the star remained constant since 1991 and was P = 115 days. We have studied the wavelength dependence of the extinction at various phases of the light curve. The blueing of the B-V color at deep minima is interpreted as the result of light scattering in the circumstellar dust shell of the star formed by preceding dust ejections since 1992. Our spectroscopic observations performed on nine nights in 2003–2007 with the 125-cm telescope at the Crimean Station of the Sternberg Institute have confirmed the previously detected intensity variations of the Swan bands and the sodium doublet with brightness. It is noted that the Swan bands originate in the upper atmosphere, the star’s extended envelope, while the sodium doublet originates mainly in the circumstellar shell of FG Sge. We suggest that the star is currently located in the temperature-luminosity diagram at the turning point of the horizontal track of cooling in the direction of hot stars—evolution caused by the last helium shell flash at the planetary nebula stage.     

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.     

Results of a long-term monitoring of the multiple system SZ Cam

We present the results of the reduction of our photometric and spectroscopic observations for the eclipsing binary SZ Cam performed with the telescopes at the Astronomical Observatory of the Ural Federal University and the Special Astrophysical Observatory of the Russian Academy of Sciences in 1996–2014. Based on an 11-year-long photometric monitoring of SZ Cam, we have obtained new elements of its photometric orbit and parameters of its components. We have detected low-amplitude periodic light variations in SZ Cam that are possibly related to the ellipsoidal shape of the components of the spectroscopic binary third body. Based on published data and our new spectroscopy, we have found new values for the mass ratio, q = 0.72 ± 0.01, and parameters of the radial velocity curves of the components, V ₀ = ?3.6 ± 1.7 km s^?1, K ₁ = 190.2 ± 1.9 km s^?1, and K ₂ = 263.0 ± 2.4 km s^?1. The component masses have been estimated to be M ₁ = 16.1 M _⊙ and M ₂ = 11.6 M _⊙. We have obtained new light elements and parameters of the radial velocity curves for the third body, V ₀ ^3b = 4.2 ± 0.6 km s^?1 and K ₁ ^3b = 26.6 ± 0.8 km s^?1. We have improved the period of the relative orbit of SZ Cam and the third body, P _orb = 55.6 ± 1.5 yr.     

Variability of the period of the star DU Monocerotis,an RR Lyrae variable with the Blazhko effect

In 2012–2014 we obtained 3641 CCD frames of the fields of the RR Lyrae (AB subtype, P = 0.583 days) variable DU Mon with BV I _c filters using the 76-cm telescope of the South African Astronomical Observatory (SAAO) and the 1-m telescopes of the Las Cumbres Observatory Global Telescope Network (LCOGT). Our observations confirmed the presence of the Blazhko effect that we suspected previously and allowed its period to be determined, \({P_{Bl}} = 60_ \cdot ^d52 \pm 0_ \cdot ^d03\). Using all of the available observations, we constructed an O–C diagram spanning a time interval of 86 years that revealed at least one abrupt change in the pulsation period (a decrease by 15.26 s).     

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.     

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_⊙.     

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.     

Irregular dwarf galaxy KK 230: A relatively young galaxy at the periphery of the local group?

We present some results of the photometric analysis of the stellar population of the irregular dwarf galaxy KK 230 on the basis of the archive database of the Hubble space telescope. The color index-magnitude diagram for KK 230 gets to magnitude 27 ^m in the V and I bands, and it comprises stellar populations of various ages. The age of the youngest main-sequence stars is 3.2 × 10⁷ yr. These stars are distributed along the north-south direction in the picture plane, and this fact can be linked to the observed kinematics of the neutral gas in the galaxy. Older blue and red supergiants are no less than 1.6 × 10⁸ years old, and such an age implies that the star formation was episodic over the last several hundreds of millions of years. As judged from the position of the tip of the red giant branch, the distance modulus for KK 230 is m ? M = 26.5 ^m . The corresponding distance is D = 2 Mpc. Based on the average absolute magnitude M _I,RC and color index (V ? I)_I,RC of the red clump, we conclude that the majority of KK 230 stars have an age of no more than (2–3) × 10⁹ yr, their metallicity being Z ≈ 0.0004.     

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.     

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.