Search for evolutionary changes in Cepheid periods using the Harvard plate collection: NSV 9159

We have obtained 530 photographic magnitude estimates for the long-period classical Cepheid NSV 9159 (P = 39^d) in the plate collections of the Harvard Observatory and the Sternberg Astronomical Institute. Together with the currently available CCD observations from the ASAS-3 catalog, our data have allowed us to construct an O-C diagram spanning a time interval of 119 years. The O-C diagram has the shape of a parabola, which has made it possible to determine for the first time the quadratic light elements and to calculate the rate of evolutionary decrease in the period, 314.4 (±7.3) s yr^?1, in agreement with the results of theoretical calculations for the second crossing of the instability strip. The available data reduced by the Eddington-Plakidis method do not reveal any noticeable random fluctuations in the period.     

Search for evolutionary changes in the periods of Cepheids: U Sgr

For the classical Cepheid U Sgr, we have constructed an O - C diagram spanning a time interval of 144 years. The O - C diagram has the shape of a parabola, which has made it possible to determine for the first time the quadratic light elements and to calculate the rate of evolutionary increase in the period, dP/dt = 0.39 (±0.10) s yr^?1, in agreement with the results of theoretical calculations for the third crossing of the instability strip. The available data reduced by the Eddington-Plakidis method reveal small random period fluctuations that do not distort the evolutionary trend in the O - C residuals.     

The Cepheid impostor HD 18391 and its anonymous parent cluster

New and existing photometry for the G0 Ia supergiant HD 18391 is analyzed in order to confirm the nature of the variability previously detected in the star, which lies off the hot edge of the Cepheid instability strip. Small‐amplitude variability at a level of δV = 0.016 ± 0.002 is indicated, with a period of P = 123^d.04 ± 0^d.06. A weaker second signal may be present at P = 177^d.84 ± 0^d.18 with δV = 0.007 ± 0.002, likely corresponding to fundamental mode pulsation if the primary signal represents overtone pulsation (123.04/177.84 = 0.69). The star, with a spectroscopic reddening of E_B–V = 1.02 ± 0.003, is associated with heavily‐reddened B‐type stars in its immediate vicinity that appear to be outlying members of an anonymous young cluster centered ∼10′ to the west and 1661 ± 73 pc distant. The cluster has nuclear and coronal radii of r_n = 3.5′ and R_c = 14′, respectively, while the parameters for HD 18391 derived from membership in the cluster with its outlying B stars are consistent with those implied by its Cepheid‐like pulsation, provided that it follows the semi‐period‐luminosity relation expected of such objects. Its inferred luminosity as a cluster member is M_V = –7.76 ± 0.10, its age (9 ± 1) × 10⁶ years, and its evolutionary mass ∼19 M_⊙. HD 18391 is not a classical Cepheid, yet it follows the Cepheid period‐luminosity relation closely, much like another Cepheid impostor, V810 Cen (© 2009 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Eccentric Eclipsing Binary V889 Aquilae

Several new times of minimum light recorded with photoelectric or CCD means have been gathered for the eccentric eclipsing binary V889 Aql (P = 11.1 days, e = 0.37). Its O–C diagram is presented, and improved elements of the apsidal motion and the light-time effect (LITE) are given. We found a long apsidal motion period of about 24 400 ± 2400 years and a period of the third body of about 52 ± 2 years.     

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.     

The vector approach to the problem of physical libration of the Moon: the linearized problem

A flexible and informative vector approach to the problem of physical libration of the rigid Moon has been developed in which three Euler differential equations are supplemented by 12 kinematic ones. A linearized system of equations can be split into an even and odd systems with respect to the reflection in the plane of the lunar equator, and rotational oscillations of the Moon are presented by superposition of librations in longitude and latitude. The former is described by three equations and consists of unrestricted oscillations with a period of T ₁ = 2.878 Julian years (amplitude of 1.855″) and forced oscillations with periods of T ₂ = 27.201 days (15.304″), one stellar year (0.008″), half a year (0.115″), and the third of a year (0.0003″) (five harmonics altogether). A zero frequency solution has also been obtained. The effect of the Sun on these oscillations is two orders of magnitude less than that of the Earth. The libration in latitude is presented by five equations and, at pertrubations from the Earth, is described by two harmonics of unrestricted oscillations (T ₅ ≈ 74.180 Julian years, T ₆ ≈ 27.347 days) and one harmonic of forced oscillations (T ₃ = 27.212 days). The motion of the true pole is presented by the same harmonics, with the maximum deviation from the Cassini pole being 45.3″. The fifth (zero) frequency yields a stationary solution with a conic precession of the rotation axis (previously unknown). The third Cassini law has been proved. The amplitudes of unrestricted oscillations have been determined from comparison with observations. For the ratio $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 $ \frac{{\sin I}} {{\sin \left( {I + i} \right)}} \approx 0.2311 , the theory gives 0.2319, which confirms the adequacy of the approach. Some statements of the previous theory are revised. Poinsot’s method is shown to be irrelevant in describing librations of the Moon. The Moon does not have free (Euler) oscillations; it has oscillations with a period of T ₅ ≈ 74.180 Julian years rather than T ≈ 148.167 Julian years.     

New CCD photometry of the eclipsing binary system V1067 Her

We present a new set of CCD photometric observations for the short period eclipsing binary 1SWASP J1743 (= V1067 Her). We have determined the available times of light minima and two new linear and quadratic ephemerides have been obtained. The photometric solutions for the system have been performed using Wilson and Devinney Code. The 3D and fill out configuration revealed that V1067 Her is an over contact W UMa binary with relatively low fill-out factor of about 16%.We investigated the period variation for the system. It showed a strong evidence of period changes by using the (O-C) residual diagram method and we have concluded long-term orbital period decrease rate dP/dt= −3.0 × 10⁻⁷ d/yr, corresponding to a time scale 8.6 × 10⁵ yr. Such period decrease in the A-type W UMa systems is usually interpreted to be due to mass transfer from the more to the less massive component.     

Measurement of the orbital period of the X-ray burster GS 1826-238 based on observations of its optical brightness variations

The variability of the optical and X-ray fluxes from the binary GS 1826-238 is investigated. An epoch-folding analysis of the optical data obtained with the RTT-150 telescope in 2003–2004 has revealed periodic brightness variations in the source with a period P _orb = 2.24940 ± 0.00015 h with a high statistical significance. When estimating the detection significance of the periodic signal, we have specially taken into account the presence of a powerful aperiodic component (“red noise”) in the source’s brightness variability. The source’s power density spectra in the frequency range ∼10⁻⁵–0.01 Hz have been obtained. We have detected a statistically significant break in the power density spectrum of GS 1826-238 at a frequency ν _br ≈ (8.48 ± 0.14) × 10⁻⁵ Hz in both optical and X-ray energy bands. We have estimated the orbital period of the binary GS 1826-238 using the correlation between the break frequency in the power density spectrum and the orbital period of binaries, P _orb ∝ 1/ν _br, found by Gilfanov and Arefiev (2005): P _orb = 3.7 ± 0.8 h and P _orb = 11.3 ± 5.9 h when using Sco X-1 and 1H 16267-273, respectively, as reference sources. It seems to us that the method for estimating the orbital periods of low-mass X-ray binaries using the correlation P _orb ∝ 1/ν _br may turn out to be very promising, especially for persistent low-luminosity X-ray binaries.     

Photometric study of two β Lyr-type binaries: DD Aqr and RR Lep

This paper presents detailed analysis of photometric observations of two eclipsing binary systems, DD Aqr and RR Lep. The V light curve of the neglected binary star DD Aqr from the All Sky Automated Survey was solved for the first time. The 1982–1987 UBV light curves of RR Lep from Vyas and Abhyankar (1989) were re-analysed. The final solutions give these two β Lyr-type binary stars as having near contact configurations in which the secondary components almost fill their Roche limiting lobes. Using O-C residuals formed by the updated minima times, orbital period changes of the systems were analysed. The O-C diagram of DD Aqr displays a cyclic variation, while that of RR Lep shows a quasi-sinusoidal variation superimposed on a downward parabolic form. The parabolic variation, which suggests a secular orbital period decrease in RR Lep, was interpreted in terms of the combined effect of mass transfer and loss. The cyclic O-C variations were interpreted in terms of the light travel time effect due to unseen components in these two systems. The absolute parameters of the components of the systems were estimated, and their present evolutionary status is also discussed.     

Light and period study of the ecipsing binary system TW Andromedae

Photoelectric light curve (LC) solutions of the close binary system TW And were obtained using the PHOEBE program (version 0.31a). Absolute parameters of the stellar components were then determined, enabling us to discuss the structure and evolutionary status of TW And. The configuration of the system based on the LCs solutions indicates that the secondary component is slightly detached from its critical Roche surface. In addition, times of minima data (“O –C curve”) were analyzed. Apart from an almost parabolic variation of the general trend of the O –C data, indicative of a secular increase in the orbital period with a rate 0.032 s yr^–1, which was attributed to a mass transfer with a rate of Δm₂ = –1.10 × 10^–10 M_⊙ yr^–1. Additionally, a sinusoidal variation with a period of 52.75 ± 1.80 yr, modulating the orbital period, was found, which we attribute to a third body orbiting the system. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The new period of the intermediate polar V709 Cas

We present the results of 10 years of photometric CCD observations of the intermediate polar V709 Cas obtained by using different instruments during 2003–2013. We detected a new variability with a period of P_new = 0.^d016449979(5) which seems to be real. The spin variability is not clearly seen in all our data, so we are unable to study any evolution of the white dwarf rotation. From the best night (in 2010) we obtained a spin period of P_spin = 311.^s8(5). We analyzed the orbital variability using (O – C) analysis. We found no variations of the orbital period on a timescale of 10 years, but the linear fit to the (O – C) diagram shows that the value of the orbital period is P_orb = 0.^d2222123(6), which is close to the earlier published values. (© 2014 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Determination of parameters of Long-Term variability of the X-ray pulsar LMC X-4

We have investigated the temporal variability of the X-ray flux measured from the high-mass X-ray binary LMCX-4 on time scales from several tens of days to tens of years, i.e., exceeding considerably the orbital period (~1.408 days). In particular, we have investigated the 30-day cycle of modulation of the X-ray emission from the source (superorbital or precessional variability) and refined the orbital period and its first derivative. We show that the precession period in the time interval 1991–2015 is near its equilibrium value P _sup = 30.370 days, while the observed historical changes in the phase of this variability can be interpreted in terms of the “red noise” model. We have obtained an analytical law from which the precession phase can be determined to within 5% in the entire time interval under consideration. Using archival data from several astrophysical observatories, we have found 43 X-ray eclipses in LMC X-4 that, together with the nine eclipses mentioned previously in the literature, have allowed the parameters of the model describing the evolution of the orbital period to be determined. As a result, the rate of change in the orbital period ? _orb/P _orb = (1.21 ± 0.07) × 10^?6 yr^?1 has been shown to be higher than has been expected previously.

     

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.     

A preliminary radial-velocity curve for the star <Emphasis Type="Italic">θ</Emphasis><Superscript>1</Superscript> Ori D

We measured the radial velocity of the star θ¹ Ori D from IUE spectra and used published observations. Based on these data, we determined the period of its radial-velocity variations, P=20.2675±0.0010 days, constructed the phase radial-velocity curve, and solved it by least squares. The spectroscopic orbital elements were found to be the following: the epoch of periastron passage E_p=JD 2430826.6±0.1, the system's center-of-mass velocity /Gg=32.4±1.0 km s^?1, K=14.3±1.5 km s^?1, Ω=3.3±0.1 rad, e=0.68±0.09, a₁ sin i = 3 × 10¹⁰ km, and f₁ = 0.0025M_⊙. Twice the period, P=40.528±0.002 days, is also consistent with the observations.     

Spectrophotometry of Saturn and its rings from 60 to 180 microns

We observed Saturn at far-infrared and submillimeter wavelengths during the Earth's March 1980 passage through the plane of Saturn's rings. Comparison with earlier spectroscopic observations by D. B. Ward [Icarus32, 437–442 (1977)], obtained at a time when the tilt angle of the rings was 21.8°, permits separation of the disk and ring contributions to the flux observed in this wavelength range. We present two main results: (1) The observed emission of the disk between 60 and 180 μm corresponds to a brightness temperature of 104 ± 2°K; (2) the brightness temperature of the rings drops approximately 20°K between 60 and 80 μm. Our data, in conjunction with the data obtained by other observers between 1 μm and 1 mm, permit us to derive an improved estimate for the total Saturnian surface brightness of (4.84 ± 0.32) × 10^?4W cm^?2 corresponding to an effective temperature of 96.1 ± 1.6°K. The ratio of radiated to incident power, P_R/P_I, is (1.46 ± 0.08)/(1 - A), where A is the Bond albedo. For A = 0.337 ± 0.029, P_R/P_I = 2.20 ± 0.15 and Saturn's intrinsic luminosity is L_S = (2.9 ± 0.5) × 10^?10L_⊙.     

Magnetic field structure and accretion regimes of the asynchronous polar by Cam

Series of photometric CCD observations of the asynchronous polar BY Cam in a low accretion state (R = 14^m–16^m) were made on the K-380 telescope at the Crimean Astrophysical Observatory (CrAO) over 100 hours in the course of 31 nights during 2004–2005. A period of P ₁ = 0.137120±0.000002 days was found for the variations in the brightness, along with less significant periods of P ₂ = 0.139759±0.000003 and P₃ = 0.138428±0.000002 days, where P₂ and P₃ are obviously the orbital and rotation periods, while the dominant period P₁ is the sideband period. A modulation in the brightness and an amplitude of 0.137 days in the oscillations at the orbital-rotational beat period (synodic cycle) of 14.568±0.003 day are found for the first time. The profile of the modulation period is four humped. This indicates that the magnetic field has a quadrupole component, which shows up well during the low brightness state. Accretion takes place simultaneously into two or three accretion zones, but at different rates. The times of the times of maxima for the main accretion zone vary with the phase of the beat period. Three types of variation of this sort are distinguished: linear, discontinuous, and chaotic, which indicate changes in the accretion regimes. At synodic phases 0.25 and 0.78 the bulk of the stream switches by 180°, and at phase 0.55, by ∼75°. At phases of 0.25–0.55 and 0.55–0.78, the O-C shift with a period of 0.1384 days, which can be explained by a retrograde shift of the main accretion zone relative to the magnetic pole and/or a change in the angle between the field lines and the surface of the white dwarf owing to the asynchronous rotation. For phases of 0.78–1.25 the motion of the accretion zone is quite chaotic. It is found that synchronization of the components occurs at a rate of less than dP_rot/P_rot∼10⁻⁹ day/day. __________ Translated from Astrofizika, Vol. 49, No. 1, pp. 121–137 (February 2006).     

Nonstationarity of the <Emphasis Type="Italic">α</Emphasis> CYG atmosphere: IV. Some peculiarities of time variation in line profiles

We use 240 CCD spectra taken in 1998–2000 with the coude echelle spectrograph of the 2-m telescope of the National Academy of Sciences of Azerbaijan to study temporal radial velocity and line profile variations of the ion, HeI, and Hβ lines in the spectrumof the α Cyg supergiant. We demonstrate that these variations are caused by pulsation-type motions in the star’s atmosphere. Ion and HeI lines oscillate in the main fundamental mode with a period of about 12.0 ± 0.5 ^d and an amplitude of 5.0 ± 0.5 km/s. These ion-line oscillations continue for about 35 days. Then the difference between the radial velocities of strong and weak ion lines results in a gradual decay of oscillations over a time interval of about 5.0 ± 1.0 ^d . Thereafter the process repeats itself. For the Hβ line we found two significant periods, two amplitudes, and three characteristic radial velocity variability behaviors for the blue and red halves of the absorption profile: with equal variability parameters (period P and amplitude A); with equal P and A, but with a phase shift between the radial velocity variations of the blue and red halves of the absorption profile; with different P and A for the two halves of the absorption profile. The star’s center of mass radial velocity as inferred from the γ-velocity is −4.5 ± 0.5 km/s. The average expansion velocity of the atmospheric layers, where the Hβ line forms, amounts to about −16.5 ± 0.5 km/s and varies temporally with an amplitude of about 3.0 km/s.     

The simulation of the magnetic field and spin period evolution of accreting neutron stars

By means of the Monte Carlo method, we simulate the evolutionary distribution of accreting neutron stars (NSs) in the magnetic field versus spin period (B‐P) diagram where the accretion induced magnetic‐field decay model is exploited. The simulated results show that by mass accretion the B‐P distribution of the accreting NS would evolve along the equilibrium period line to a region with low field and short period. The B‐P distributions of the simulated accreting NSs are consistent with those of the observed millisecond pulsars (MSPs) after accretion of ∼ 0.1–0.2 M⊙. We also test the effects of the initial magnetic field and the spin period on the evolved B‐P distribution of the accreting NSs. It is shown that the evolved distributions of the simulated samples are independent of the selection of the initial condition when the NS magnetic field decays to a value less than ∼10¹⁰ G. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Long‐term orbital period behaviors of the neglected Algol type binaries: CC Herculis and XZ Aquilae

Orbital period variations of two neglected Algol type binaries, CC Her and XZ Aql, are studied based on all available times of minima. In the case of CC Her, it is found that the O –C curve displays a tilted sinusoidal variation with an eccentricity of 0.54 ± 0.03 and a period of 52.4 ± 0.4 yr, which can be explained by the light‐time effect due to the presence of an unseen component. The course of the orbital period change in XZ Aql appears less reliable but its O –C curve can be represented by a periodic variation with a period of 36.7 ± 0.6 yr superimposed on an upward parabola. The parabolic variation indicates a secular period increase with a rate of dP /dt = 7.1 s per century. The corresponding conservative mass transfer from less massive component to the more massive one is about 3.26 × 10^–7 M_⊙ yr^–1. It is interesting to see that the O –C variation of CC Her displays no evidence (as upward parabola) on the mass transfer characteristic for Algols. The periodic change of the orbital period of XZ Aql, like CC Her, may be caused by the presence of the thirdbody. The lower limits of the masses of the hypothetical unseen components for CC Her and XZ Aql are found to be 2.69 M_⊙ and 0.47 M_⊙, respectively. The third body of CC Her should be detectable not only spectroscopically but also photoelectrically, if it exists. (© 2006 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)