Short-period variations in the proper motions of sunspots from SOHO (MDI) observations

Short-period (1–60 min) variations in the coordinates of the centers of gravity of isolated sunspots are analyzed. The sunspot coordinated were determined using two sets of observational data—magnetograms and intensities—obtained by SOHO (MDI) on December 6, 1998, from 01:00 to 21:57 UT with temporal resolution 60 s and spatial resolution 0.6″/pixel. A slow drift in the sunspot coordinates was removed using a low-frequency filter with a 61-min integration window. The guiding errors (RMS～0.014″) were determined by analyzing correlated motions in pairs of sunspots, and were removed from the time series before determining the sunspot proper motions. Based on the calculated power spectra for the sunspot proper motions, two period intervals containing appreciable power were identified. One coincides with the well-known 5-min acoustic solar oscillations. The concentration of power in this interval is greater for the coordinate variations derived the magnetograms than those derived from the intensities; the harmonic amplitude for some peaks reaches ～±30 km. The other spectral interval corresponds to periods exceeding 30 min. Overall, the rms short-period variations in the sunspot proper motions are 9.9±2.2 and 16.7±7.6 km (0.014″±0.003″ and 0.024″±0.010″) for the magnetogram and intensity data, respectively.     

Physical parameters and dynamical properties of the multiple system ι UMa (ADS 7114)

We analyze the physical parameters, orbital elements, and dynamic stability of the multiple system ?? UMa (HD 76644 = ADS 7114). We have used the positions from the WDS catalog and our own observations on the 6-m telescope of the Special Astrophysical Observatory of the Russian Academy of Sciences and the 1.5-m Russian-Turkish Telescope (Antalya, Turkey). We have obtained more precise orbital parameters of the subsystems, and spectral types, absolute magnitudes, and masses of the components. The primary has Sp = F0 V?CIV, M = 1.7 ± 0.1M _??, T _eff = 7260 ± 70 K, and log g = 4.30 ± 0.07. The companion in the close Aa subsystem is most likely a white dwarf with a mass of approximately 1.0 ± 0.3M _??. The spectral types and masses of the components in the BC subsystem are M3V, M4V and 0.35 ± 0.05M _??, 0.30 ± 0.05M _??, respectively. The total mass is 3.4 ± 0.4M _??. The Aa subsystem probably has an orbital period of 4470^d = 12.2^y and an eccentricity of approximately 0.6. The outer subsystem seems to have a period of approximately 2084 yrs and an eccentricity of approximately 0.9. We have carried out simulations using the stability criteria and shown that for all possible variations in the component parameters, the multiple system is unstable on a time scale of less than 10⁶ years with a probability exceeding 0.98. Possible reasons for this instability are discussed.     

Spot activity of the binary CG Cyg: Observations during maximum spottedness

We present the results of new multicolor photometry of the chromospherically active binary CG Cyg acquired in 2005–2009 (136 hours of observations). The light curves for each season reveal rotational brightness modulation due to spots that varies in amplitude and phase from season to season. We have determined the longitudes of spotted areas: for each season, they were located on the primary, close to the line joining the centers of the components. The longitude distribution of spots was analyzed for 44 years of observations of CG Cyg using all data available in the literature. The active longitudes of CG Cyg are not fixed at the quadratures, as was believed earlier: most of the time (1965–2003), the spots were concentrated at two active longitudes at the quadratures, at orbital phases 0.28±0.06 and 0.70±0.08, but, during a shorter time after 2004, they were located along the line joining the component centers, at orbital phases 0.50 ± 0.04 and 0.93 ± 0.03. We detected a switch of the active longitude by 180° during 1.5 months in 2008, accompanied by an increase in the amplitude of the rotational brightness modulation and only a slight increase in the star’s spotted area (by 5%). Our analysis of archive data reveals that switches of the active longitude during time intervals of 1–1.5 months were observed three times during the entire observational history of CG Cyg (in 1991, as well as in our observations of 2003 and 2008). All these switches were accompanied by similar phenomena: an increase in the amplitude of the rotational brightness modulation (by 0.06 ^m , 0.02 ^m , and 0.04 ^m ) and an increase in the spotted area (by 79%, 11%, and 5%). We used a zonal spot model to reconstruct the parameters of the spotted regions on CG Cyg. At all our observing epochs, the spots were located at low latitudes, in a region that was symmetric about the equator, 10° to 14° wide on either side. The spots are cooler than the surrounding photosphere by 2000 K. The spotted area varied only slightly from season to season, comprising 13%–15% of the surface area of the star, close to the historic spottedness maximum for the CG Cyg system.     

A photometric and spectroscopic study of the B[e]/X-ray transient CI Cam

We compare U BV R photometry of CI Cam obtained using facilities at the Sternberg Astronomical Institute in 1998–2001 and medium-resolution spectroscopy obtained at the Special Astrophysical Observatory in the same period. The photometric measurements and fluxes in the Balmer and some Fe II emission lines in the post-outburst quiescent spectrum of CI Cam display cyclic variations with a period of 1100±50^d. These variations could be due to orbital motion in a widely separated binary; one component, a giant, produces a reflection effect on the surface facing the other, compact component. Our V-band photometry confirms the \(11\mathop d\limits_. 7\) period found by Miroshnichenko before the outburst but shows a lower amplitude (3%). We investigate the possibility that this photometric period is the same as the rotational period of the radio jets measured using the VLA. Based on modeling of the radio map, we estimate the inclination of the jet rotation axis to the line of sight, i=35°–40°; the angle between the rotational axis and jet ejection direction, ?=7°–10°; and the jet spatial velocity, 0.23–0.26 c. The equivalent widths and fluxes of various spectral lines show different variation amplitudes during outburst and significantly different behavior during quiescence. We have detected five patterns of behavior, which provide evidence for stratification of the gas and dust envelope surrounding the system. The brightening of the [N II] forbidden line is delayed by 50–250^d relative to the X-ray outburst maximum.     

Photometry and spectroscopy of CI camelpardalis from 1998–2005

Our long-time monitoring of the B[e] star and transient X-ray source CI Cam during quiescence following the 1998 outburst demonstrates that the complex, stratified circumstellar envelope has tended to stabilize after this structure was perturbed by the passage of a shock wave from the outburst. The star’s U BV R brightness shows slow, possibly cyclic, variations with an amplitude of about 0.2^m. We determined the spectral type of the primary, B4III-V, based on the widths of the absorption wings of high-numbered Balmer lines. A Doppler shift of 460 km/s was detected for the Hell λ4686 Å emission line. The shifts in this line yield an orbital period of 19.41 days, which is also manifested itself in the photometric data as a wave with a V amplitude of 0.034^m. The orbit is elliptical, with an eccentricity of 0.62. It is most likely that the secondary is a white dwarf surrounded by an accretion disk. The primary’s mass exceeds 12 M _⊙ . The system may be at a late stage of its evolution, after the stage of mass exchange.     

Features of the matter flows in the peculiar cataclysmic variable AE Aquarii

The structure ofmatter flows in close binary systems in which one of the components is a rapidly rotating magnetic white dwarf is studied. Themain example considered is the AEAquarii system; the period of the white dwarf’s rotation is about a factor of 1000 shorter than the orbital period, and the magnetic field on the white-dwarf surface is of order 50MG. The matter flows in this system were analyzed via numerical solution of a systemofmagnetohydrodynamical equatons. These computations show that the white dwarf’s magnetic field does not significantly influence the velocity field of the matter in its Roche lobe in the case of a laminar flow regime, so that the field does not hinder the formation of a transient disk (ring) surrounding the magnetosphere. However, the efficiency of the energy and angular-momentum exchange between the white dwarf and the surrounding matter increases considerably with the development of turbulent motions in the matter, accelerating the matter at the magnetosphere boundary and leading to a high escape rate from the system. The time scales for the system’s transition between the laminar and turbulent modes are close to those for the transition of AE Aquarii between its quiet and active phases.     

A photometric study of the eclipsing dwarf nova GY Cnc in quiescence and during an outburst

The results of photometric observations of the dwarf nova GY Cnc in the Rc filter acquired in 2013–2015 (~3900 orbital cycles, 19 nights in total) are presented, including observations during its outburst in April 2014. The binary’s orbital elements have been refined. The orbital period has changed only insignificantly during the ~30 000P_orb since the earlier observations; no systematic O–C variations were detected, only fluctuations within 0.004^d on time scales of 1500–2000P_orb. A “combined” model is used to solve for the parameters of GY Cnc during two states of the system. The flux from the white dwarf is negligible due to the star’s small size. The temperature of the donor star, T₂ ~ 3667 K (Sp M0.2 V), varies between 3440 and 3900 K (Sp K8.8–M1.7 V). The semi-major axis of the disk is a ~ 0.22a₀, on average. In quiescence, a varies within ~40%. The disk has a considerable eccentricity (e ~ 0.2?0.3) for a < 0.2a₀. The disk shape becomes more circular (e < 0.1) with increasing a. The outburst of GY Cnc was associated with increased luminosity of the disk due to the parameter α_g (related to the viscosity of the disk material) decreasing to 0.1–0.2 and the temperature in the inner parts of the disk increasing twofold, to T_in ~ 95 000 K. These changes were apparently due to the infall of matter onto the surface of the white dwarf as the outburst developed. All parameters of the accretion disk in quiescence display considerable variations about their mean values.     

Photometric period and rotational brightness modulation of V833 Tau

Photometric data covering 1994–2009, obtained at the Crimean Astrophysical Observatory or retrieved from the ASAS and SuperWASP catalogs, were used to analyze the brightness variations of the rapidly rotating star V833 Tau, whose activity level is close to saturation. Combined with previously published results, these data represent for the first time all stages in the development of the star’s 19- year activity cycle. The photometric period and rotational-modulation amplitude for different epochs are determined, and the qualitative pattern of the spots is considered. The photometric period is close to the orbital period, but always exceeds it, indicating that surface inhomogeneities are located far from the equatorial plane. With the high spottedness of the star, reaching 28% at the cycle maximum, the rotationalmodulation amplitude is 0.05?0.1 ^m , and increases during the growth and decline phases of the 19-year cycle. The rotational modulation is due to spots with higher latitudes than in the case of the Sun, and concentrating on active longitudes.     

HD 181906: A star with two systems of active longitudes on its surface

We use continuous 156-day COROT photometric observations of the F dwarf HD 181906 to analyze temperature inhomogeneities on the stellar surface and follow their evolution. The analysis used the iPH code, which solves for the temperature inhomogeneities in a two-temperature approximation without any assumptions about the shape or number of spots on the surface. For the first time in studies of active regions and active longitudes, we find that the phases of the active longitudes on the surface of HD 181906 are concentrated close to two systems of active longitudes. In each system, the active longitudes are separated by 180°, with the shift between the systems being 100°. During the observing period, switches between the systems of active longitudes occurred quasi-periodicly on a time scale of 30–35 day about two-thirds of the time, while these switches occurred less frequently about one-third of the time. The positions of active regions switched either every 20–25 day or every 40–45 day. The periodicity of brightness and amplitude variations is of the order of 110 day. Variations in the spot coverage and changes in the active longitudes have shorter time scales (about 55–75 day). All these parameters are variable on time scales of 25–38 day. A wavelet analysis of the periodicity of the brightness variations indicates that all the above processes are quasi-periodic; activity on all time scales became less pronounced in the last third of the observing period.     

Full 3D MHD calculations of accretion flow structure in magnetic cataclysmic variables with strong,complex magnetic fields

We have performed three-dimensional magnetohydrodynamical calculations of stream accretion in cataclysmic variable stars for which the white dwarf primary possesses a strong, complex magnetic field. These calculations were motivated by observations of polars: cataclysmic variables containing white dwarfs with magnetic fields sufficiently strong to prevent the formation of an accretion disk. In this case, an accretion stream flows from the L1 point and impacts directly onto one or more spots on the surface of the white dwarf. Observations indicate that the white dwarfs in some binaries possess complex (non-dipolar) magnetic fields. We performed simulations of ten polars, with the only variable being the azimuthal angle of the secondary with respect to the white dwarf. These calculations are also applicable to asynchronous polars, where the spin period of the white dwarf differs by a few percent from the orbital period. Our results are equivalent to calculating the structure of one asynchronous polar at ten different spin-orbit beat phases. Our models have an aligned dipolar plus quadrupolar magnetic field centered on the whitedwarf primary. We find that, with a sufficiently strong quadrupolar component, an accretion spot arises near the magnetic equator for slightly less than half our simulations, while a polar accretion zone is active for most of the remaining simulations. For two configurations, accretion at a dominant polar region and in an equatorial zone occurs simultaneously. Most polar studies assume that the magnetic field is dipolar, especially for single-pole accretors. We demonstrate that, with the orbital parameters and magnetic-field strengths typical of polars, the accretion flow patterns can vary widely in the case of a complex magnetic field. This may make it difficult formany polars to determine observationally whether the field is pure dipolar or is more complex, but there shoulid be indications for some systems. In particular, a complex magnetic field should be suspected if there is an accretion zone near the white dwarf’s equator (assumed to be in the orbital plane) or if there are two or more accretion regions that cannot be fitted by dipolar magnetic field. Magnetic-field constraints are expected to be substantially stronger for asynchronous polars, with clearer signs of complex field geometry due to changes in the accretion flow structure as a function of azimuthal angle. These indications become clearer in asynchronous polars because each azimuthal angle corresponds to a different spin-orbit beat phase.     

Giant pulses of PSR B0531+21 at 112 MHz

Three series of 111.88 MHz observations of giant pulses of PSR B0531+21 have been carried out in 2005 and 2007. The scattering of pulses observed in various series varies by a factor of 1.7: 10.6±0.5 ms in November 2005, 18±1 ms in January 2007, and 16±0.8 ms in June 2007. The cumulative probability distribution for the peak intensities of the giant pulses for each of these series shows that the distribution is stable and is a power law with a single slope (n = 2.3). This testifies to stability of the mechanism generating the giant pulses. The distribution functions for the 2005 and 2007 data can be superposed after correcting the intensities with a coefficient equal to the ratio of the effective pulse widths. Consequently, in the range of 23MHz-9GHz the energy in the pulses is conserved; i.e., the increase in the pulse intensity is proportional to the decrease in the scattering. Refractive scintillations at low frequencies in measurements with large time separation lead to variations in the number of giant pulses exceeding a given amplitude, proportional to the ratio of the mean flux densities of the pulsar in the corresponding observational series. The maximum energy of the recorded giant pulses is 2.5 × 10⁷ Jy µs. A comparison with the statistical properties of the giant pulses observed at other frequencies shows that the frequency dependence of the maximum energy of the giant pulses in the range of 23 MHz-9GHz is a power-law with index 2.2±0.2. The degree of linear polarization of the giant pulses at 112 MHz does not exceed 12%.     

IR photometry of the symbiotic star BF Cyg: Detection of the red giant’s ellipsoidal brightness variability

We present the results of our IR photometric observations of the classical symbiotic star BF Cyg acquired in 1978–2003. The variability range in the J and K bands was ～0.2^m. A periodic component in the cool star’s brightness variations is clearly visible, its period being half the orbital one and its J amplitude being ～0.15^m. This component is associated with the ellipsoidal shape of the red giant, which model calculations show fills its Roche lobe. This is required in order to reproduce ellipsoidal brightness variability with such a large amplitude: the calculated amplitude for a red giant filling 90% of its Roche lobe is half the observed value. At the same time, it was not possible to confidently chose the optimum component-mass ratio, q = M _giant /M_hot, and orbital inclination, i, from possible values in the ranges q = 2–4, i = 70°–90°. Including the contribution from the hot radiation sources (the hot component and ionized envelope), which vary with a period equal to the orbital period, has a considerable influence on the estimated parameters associated with the red giant’s ellipsoidal brightness variations, and this contribution cannot be neglected. The deviations of the observed from the calculated light curve are irregular, with the rms deviation being σ(O-C) ≈ 0.04^m.     

Orbital-period variations in the eclipsing binary Y Cam

We have studied the time behavior of the orbital period and the primary’s pulsation period for the eclipsing binary system Y Cam, whose secondary fills its Roche lobe and whose primary is a δ Scuti star. The times of minima available for this eclipsing binary cover 120 years. δ Scuti pulsations of the primary have been observed over the last 50 years, with the period of these pulsational brightness variations remaining virtually unchanged during the entire observed time interval. The large-amplitude cyclic variations of the orbital period of Y Cam cannot be explained solely by the presence of a third body in the system. It is possible to explain the period variations of Y Cam with magnetic oscillations or a superposition of a stationary matter flow from the lower-mass to the higher-mass component together with magnetic oscillations, similar to the case of AB Cas. A good agreement with observations is provided by a model assuming a stationary matter flow from the secondary filling its Roche lobe to the primary, at the rate of 2.85 × 10^?7 M _⊙/year, superposed with irregular period jumps that can be explained by instabilities in the matter flow. We have detected cyclic variations of the orbital period of Y Cam with an amplitude of 0.011^d, which can be understood if the binary moves in a long-period orbit (with a period of 38.6 years) around a third body with mass M ₃ s> 0.30M _⊙. These cyclic period variations of the eclipsing binary agree with the observed small period variations of the δ Scuti pulsations.     

AE Aquarii represents a new subclass of Cataclysmic Variables

We analyze properties of the unique nova-like star AE Aquarii identified with a close binary system containing a red dwarf and a very fast rotating magnetized white dwarf. It cannot be assigned to any of the three commonly adopted sub-classes of Cataclysmic Variables: Polars, Intermediate Polars, and Accreting non-magnetized White Dwarfs. Our study has shown that the white dwarf in AE Aqr is in the ejector state and its dipole magnetic moment is ???1.5 × 10³⁴ G cm³. It switched into this state due to intensive mass exchange between the system components during a previous epoch. A high rate of disk accretion onto the white dwarf surface resulted in temporary screening of its magnetic field and spin-up of the white dwarf to its present spin period. Transition of the white dwarf to the ejector state had occurred at a final stage of the spin-up epoch as its magnetic field emerged from the accreted plasma due to diffusion. In the frame of this scenario AE Aqr represents a missing link in the chain of Polars evolution and the white dwarf resembles a recycled pulsar.     

Detection of regular variations in the intensity and pulse time of arrival of the anomalous pulsar PSR B0943+10

Timing of the anomalous pulsar PSR B0943+10 during 2007–2013 was carried out on the Large Phased Array radio telescope of the Pushchino Radio Astronomy Observatory at 112 MHz. The astrometric and rotational parameters for epoch MJD=56 500 have been determined. Considerable deviations of the pulse times of arrival from the precalculated values with a characteristic period of several years due to the presence of correlated low-frequency noise in the pulsar spin phase have been detected. These deviations can be explained in a planetary model by the presence of two companions of the pulsar, whose orbital parameters have been determined. A continuous increase in the longitude of the pulse maximum within the emission window, the pulse width, and the intensity have been detected after each switch to the burst mode. Together with the changes in pulse shape, degree of linear polarization of the pulse, and drift rate of individual pulses detected earlier, this indicates that all the main parameters of the radio emission in the B mode are unstable. This distinguishes PSR B0943+10 from all other modes-witching pulsars. The origin of the observed properties of this pulsar are probably associated with the interaction of its extended magnetosphere with the surrounding medium.     

Three-Dimensional Numerical Simulation of a Flow Structure in the Asynchronous Polar CD Ind in the Approximation of an Offset Dipole Magnetic Field of a White Dwarf

The flow structure in the asynchronous polar CD Ind was investigated in the approximation of an offset dipole. Computation results made it possible to identify such system features as the drift of hot spots over the white dwarf surface, the flow structure dependence on the phase of the beat period, and the magnetic pole switching. To study the flow structure, a three-dimensional numerical MHD model based on the approximation of modified magnetic hydrodynamics was used. Numerical calculations were performed for ten phases of the beat period at a constant rate of mass transfer. In addition, for a more detailed study of the poles switching, additional calculation series were carried out in the corresponding phase ranges of the beat period. The energy release zones during the spin-orbit period shift in longitude on average by 20°, which corresponds to 0.05 of the orbital period phase. The magnetic pole switching occurs in a time that does not exceed 0.1 of the beat period. In the middle of this process, accretion occurs on both poles with the same intensity, and the flow resembles an arch. Based on the calculation results, synthetic light curves were constructed for the optical range of the spectrum.     

A Double Wave in the Orbital-period Variations of the Eclipsing Binary RT Per

Orbital-period variations of the low-mass, semi-detached eclipsing binary RT Per are analyzed. In addition to the secular variation of the orbital period determined by the mass transfer between the components, cyclic variations are also present. Both the light-time effect and magnetic oscillations can describe the cyclic orbital-period variations of RT Per. The secular period increase can be explained by a uniform flow of matter from the lower-mass to the higher-mass component at a rate of 0.60 × 10⁻⁸M-/year, with the total angular momentum being conserved. The period variations can be represented equally accurately by either a superposition of two cyclic variations or a superposition of a secular period increase and two cyclic variations. Approximately the same parameters are derived for the lower-period (36.8 year) variations when the times of minima are fitted with a linear or a quadratic formula. For the longer-period variation, a period of 275 years and amplitude of 0.104 days are found using the linear formula, or 89 years and 0.045 days using the quadratic formula.

     

Long-term variations in the rotation of the solar corona

The characteristic time scales for variations in the differential rotation of the solar corona are determined using measurements of the intensity of the FeXIV 5303 Å coronal line made from 1939–2004. Drift waves of the variations in the rotational speed with an 11-year periodicity can be distinguished. Moving averages with time intervals from two to five years are used to identify torsional waves. In addition, longer-period variations in the rotational speed can be distinguished when longer averaging intervals are used. When the interval used for the moving average is increased to 8–12 years, a quasi-22-year rotational period appears. The low-latitude corona rotates more slowly in odd cycles than in even cycles. Increasing the duration of the averaging interval further shows that rapid rotation at low latitudes was observed in 1940–1950 and 1990–2000, while slow rotation was observed in 1960–1980, possibly suggesting the presence of a 55-year period in the rotational variations. Long-term variations are found in the rotation of polar regions. The rotational variations for high-latitude corona are in antiphase with those for the low-latitude corona. The origins of zones of anomalous coronal rotation and their dynamics in the global activity cycle are discussed.     

Photometric parameters of the dwarf nova SS Cygni in the quiescent state

The mean 1983–1996 UBV light curves of the dwarf nova SS Cyg are used to derive the binary parameters in the quiescent state. Solutions are obtained for a classical hot-spot model and a model with an energy source lying outside the accretion disk. Photometric and spectroscopic data are combined to infer the masses and radii of the binary components. The white dwarf in SS Cyg is one and a half times as massive as the red dwarf, q=M _wd /M _rd ～1.45, M _rd ～0.46M _⊙ and M _wd ～0.66M _⊙. The orbital inclination of the system is i?51°–54°. The contribution of the accretion disk to the total flux in the quiescent state is estimated to be ～47–49% and ～54% in the VU and B filters, respectively. The hot spot contributes less than ～3% to the total optical flux. In the “non-classical” hot-spot model, the disk and bulge contributions are 27 and 2–8%, respectively, depending on the orbital phase. The shape of the mean light curves of SS Cyg suggests asymmetric heating of the red-dwarf surface in the quiescent state by high-temperature radiation generated in the hot-spot region.     

PG 1316+678: A young pre-cataclysmic binary with weak reflection effects

The PG 1316+678 star is classified as a pre-cataclysmic binary, as is evidenced by its photometric and spectroscopic observations. Its orbital period is determined to be P _orb = 3.3803^d, which coincides with the photometric period. The intensities of the emission HI and HeI lines are shown to vary synchronously with the brightness of the object (Δm _V = 0.065 ^m , Δm _R = 0.08 ^m ). These variations arise as the UV radiation from the DAO white dwarf is reflected from the surface of the cold companion. The parameters of the binary are estimated and the time of its evolution after the common-envelope phase is determined to be t ≈ 240 000 years. Thus, PG 1316+678 is a young pre-cataclysmic NN Ser variable with the smallest known photometric reflection effect.