A many-year photometric and spectroscopic variability study of planetary nebulae (1968–2008). The young planetary nebula IC 4997

We present the results of our many-year (1968–2008) photoelectric photometry and absolute spectrophotometry of the young variable planetary nebula IC 4997, performed in uniform systems. Integrated line energy fluxes in the range 3727–5007 Å are tabulated, along with the integrated (star + nebula) continuum flux at 4500 Å, and time variations investigated. We study the time behavior of fluxes in the hydrogen, HeI, [OIII], [OII], [NeIII] emission lines. Variations of the line intensity ratio R = F([OIII]4363 Å)/F(Hγ) between 1938 and 2008 are presented. We estimate the nebula’s electron density N _e and electron temperature T _e for 1972–1992: the mean nebula N _e increased by the factor of five, from ～4 × 10⁵ to ～2 × 10⁶ cm^?3, while the mean T _e increased from 12 000 to 14 000 K. The color temperature of the exciting star, measured from the 4000–4900 Å continuum, increased from 37 000–40 000 K to 47 000 K during the time interval covered by our observations, as is confirmed by the growth of ionization of HeI, NeIII, and OIII (from the 4363 Å line). V-band variations of the object’s integrated (nebula + star) light with an amplitude about 0.5 ^m can essentially be fully explained as being due to changing contributions from variable fluxes in the [OIII] 5007 and 4959 Å lines. The V magnitude in 2009 (after 40 years) happened to be the same as at the start of our observations in 1968. A period of the order of 50 years can also be noted in log R(t). This may provide evidence for binarity of the central star and be related to its orbital period.     

Optical identifications and spectra of radio sources

We present classifications, optical identifications, and radio spectra for eight radio sources from three flux-density-complete samples in the following declination ranges: 4°–6° (B1950), S _3.9 > 200 mJy; 10°–12°30′ (J2000), S _4.85 > 200 mJy; 74°?75° (J2000), S _4.85 > 100 mJy. For all these samples, the right ascensions are 0^h–24^h and the Galactic latitudes, |b| > 15°. Our optical observations at 4000–7500 ° were made with the 6-m telescope of the Special Astrophysical Observatory; we also observed at 0.97–21.7 GHz with the RATAN-600 radio telescope of the Special Astrophysical Observatory. We classify four of the objects as quasars and four as galaxies. Five of the radio sources have power-law spectra at 0.97–21.7 GHz, while two objects have flat spectra. The quasar J2358+0430 virtually did not vary during 23 years.     

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.     

Structure of the radio source 3C 120 at 8.4 GHz from VLBA+ observations in 2002

Maps of the radio source 3C 120 obtained from VLBA+ observations at 8.4 GHz at five epochs in January–September 2002 are presented. The images were reconstructed using the maximum entropy method and the Pulkovo VLBImager software package for VLBI mapping. Apparent superluminal motions of the brightest jet knots have been estimated. The speeds of jet knots decreases with distance from the core, changing from (5.40±0.48)c to (2.00±0.48)c over 10 mas (where c is the speed of light) for a Hubble constant of 65 km s^?1 Mpc^?1. This can be explained by interaction of the jet with the medium through which it propagates.     

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.     

Spectroscopy of RC sources

The results of spectroscopic observations of the host galaxies of objects in the RC catalog (the “Big Trio” program) obtained using the new SCORPIO spectrograph of the Special Astrophysical Observatory are presented. The spectroscopic redshifts of the objects are compared with their photometric color redshifts, and the errors in the latter are estimated. Based on BV RI observations obtained on the 6-m telescope of the SAO, the errors for the population of powerful radio galaxies are close to those found previously for radio quiet galaxies (about 10–20%). The detection of Ly α in the B filter in RC 1626+0448 is confirmed. This object is the second spectrally studied FR II radio source from the RC catalog to have a redshift z>2.5. Star formation in its host galaxy began at a redshift z>3.3. This first use of the new SCORPIO spectrograph demonstrates its promise for studies of very distant steep-spectrum radio galaxies brighter than 23^m–24^m in V.     

Seyfert nuclei as sources of ultrahigh-energy cosmic rays

Particles can be accelerated to ultrahigh energies E≈10²¹ eV in moderate Seyfert nuclei. This acceleration occurs in shock fronts in relativistic jets. The maximum energy and chemical composition of the accelerated particles depend on the magnetic field in the jet, which is not well known; fields in the range ～5–1000 G are considered in the model. The highest energies of E≈10²¹ eV are acquired by Fe nuclei when the field in the jet is B≈16 G. When B～(5–40) G, nuclei with Z<10 are accelerated to E≤10²⁰ eV, while nuclei with Z≥10 acquire energies E≥2×10²⁰ eV. Only particles with Z≥23 acquire energies E≤10²⁰ eV when B～1000 G. Protons are accelerated to E<4×10¹⁹ eV, and do not fall into the range of energies of interest for any magnetic field B. The particles lose a negligible amount of their energy in interactions with infrared photons in the accretion disk; losses in the thick gas-dust torus are also negligible if the luminosity of the galaxy is L≤10⁴⁶ erg/s and the angle between the normal to the galactic plane and the line of sight is sufficiently small, i.e., if the axial ratio of the galactic disk is comparatively high. The particles do not lose energy to curvature radiation if their deviations from the jet axis do not exceed 0.03–0.04 pc at distances from the center of R≈40–50 pc. Synchrotron losses are small, since the magnetic field frozen in the galactic wind at R≤40–50 pc is directed (as in the jet) primarily in the direction of motion. If the model considered is valid, the detected cosmic-ray protons could be either fragments of Seyfert nuclei or be accelerated in other sources. The jet magnetic fields can be estimated both from direct astronomical observations and from the energy spectrum and chemical composition of cosmic rays.     

Physical parameters of the eclipsing binary V798 Cep

The first high-accuracy CCD UBV RI(RI)_C light curves for the recently discovered eclipsing binary V798 Cep (P = 16 ^d .08, V = 11 ^m . 8) are presented; this star is included in our program of eclipsing systems with considerable eccentricities. A photometric solution for the light curves and physical characteristics of the component stars are derived. The orbital eccentricity is quite high, e = 0.437. The longitude of periastron is close to 180°, making studies of the apsidal motion difficult. V798 Cep may be a hierarchical system.     

Spectral-temporal evolution of low-frequency pulsations in the microwave radiation of solar flares

Low-frequency pulsations of 22 and 37 GHz microwave radiation detected during solar flares are analyzed. Several microwave bursts observed at the Metsähovi Radio Observatory are studied with time resolutions of 100 and 50 ms. A fast Fourier transformation with a sliding window and the Wigner-Ville method are used to obtain frequency-time diagrams for the low-frequency pulsations, which are interpreted as natural oscillations of coronal magnetic loops; the dynamical spectra of the pulsations are synthesized for the first time. Three types of low-frequency fluctuations modulating the flare microwave radiation can be distinguished in the observations. First, there are fast and slow magneto-acoustic oscillations with periods of 0.5–0.8 s and 200–280 s, respectively. The fast magneto-acoustic oscillations appear as trains of narrow-band signals with durations of 100–200 s, a positive frequency drift dν/dt=0.25 MHz/min, and frequency splitting δν=0.01–0.05 Hz. Second, there are natural oscillations of the coronal magnetic loops as equivalent electrical circuits. These oscillations have periods of 0.5–10 s and positive or negative frequency drift rates dν/dt=8×10^?3 Hz/min or dν/dt=?1.3×10^?2 Hz/min, depending on the phase of the radio outburst. Third, there are modulations of the microwave radiation by short periodic pulses with a period of 20 s. The dynamical spectra of the low-frequency pulsations supply important information about the parameters of the magnetic loops: the ratio of the loop radius to its length r/L≈0.1, the plasma parameter β≈10^?3, the ratio of the plasma densities outside and inside the loop ρ_e/ρ_i≈10^?2, and the electrical current flowing along the loop I≈10¹² A.     

Photometry and spectroscopy of the system V4641 sagittarii in quiescence

We present the results of our CCD photometric and moderate-dispersion spectroscopic observations of the binary system V4641 Sgr, which contains a black hole of mass ≈9.5M_⊙ and a normal B9III star. The photometric light curve reveals an ellipticity effect with very high amplitudes in V and R, 0.40^m and 0.37^m, and the color curve shows that the surface temperature is nonuniform. All this testifies to tidal distortion of the normal star's surface due to the massive companion and to a high inclination of the orbit to the line of sight. In June and July 2002, during quiescence, we obtained data during three flares with amplitudes up to 0.26^m. In particular, spectroscopic observations were acquired near the time of the black hole's inferior conjunction. One hour before conjunction, a depression by EW=0.5 Å was observed in the red wing of the Hα absorption line, interpreted as absorption by gas flowing in the direction from the observer toward the normal star. This flow is apparently associated with a rarefied gas disk around the black hole, and the conjunction grazes the stellar surface if the orbital inclination is close to 70.7°. The maximum velocity along a circular Keplerian orbit is 650 km/s at a distance of R=0.15–0.20a from the black hole (where a is the component separation). Thus, we find the mass of the black hole to be M _BH =7.1–9.5M_⊙, confirming the model of Orosz et al. (2001).     

The radio brightness distribution along the solar limb observed on the Siberian Solar Radio Telescope

The radio brightness distribution at 5.2 cm has been obtained along the entire solar limb down to a level of 0.01 of brightness temperature at the disk center T ₀ ^c during the solar minimum. The measurements are based on strip scans from the Siberian Solar Radio Telescope. Data reduction included model fitting of an actual solar scan for each position angle. The maximum limb brightness, 1.37 T ₀ ^c , in the derived distribution is at equatorial direction, where the radio radius exceeds the photosphere radius by 7%. Toward the poles, the brightness peak and radio radius gradually decrease to 1.01 T ₀ ^c and 1.005 R₀. The derived two-dimensional brightness distribution was used to calculate radio radius as a function of position angle for several wavelengths from 4 to 31.6 cm. These dependences are consistent with RATAN-600 observations at position angles 0° and 25°.     

IR photometry and models for the dust envelopes of two carbon stars

The results of JHKLM photometry of two carbon stars are presented: the irregular variable NQ Cas and the Mira star BD Vul. Data on the mean fluxes supplemented with mid-IR observations with the IRAS, AKARI, andWISE satellites are used to compute spherically symmetrical model dust envelopes for the stars, consisting of particles of amorphous carbon and silicon carbide. The optical depth in the visible for the comparatively cool dust envelope of BD Vul, with a dust temperature at its inner boundary T₁ = 610 K, is fairly low: τ_V = 0.13. The dust envelope of NQ Cas is appreciably hotter (T₁ = 1550 K), and has τ_V = 0.32. The estimated mass-loss rates are 1.5 × 10^?7M_⊙/yr for NQ Cas and 5.9 × 10^?7M_⊙/yr for BD Vul.     

On the correlation between spectra of solar microwave bursts and proton fluxes near the Earth

Studies of the extreme solar proton event of January 20, 2005 intensified the contest over of a long-standing problem: are solar cosmic rays arriving at the Earth accelerated by solar flares or by shocks preceding rapidly moving coronal mass ejections? Among the most important questions is the relationship between the energy spectra of the solar cosmic rays and the frequency spectra of flare microwave bursts. Some studies of previous solar-activity cycles have shown that such a relationship does exist, in particular, for protons with energies of tens of MeV. The present work analyzes this relation using data for 1987–2008. For flare events observed in the western half of the disk, there is a significant correlation between the index δ, which is equivalent to the power-law index of the integrated energy spectrum of 10–100 MeV protons detected near the Earth’s orbit, and radio burst parameters such as a ratio of peak fluxes S at two frequencies (for example, at 9 and 15 GHz) and a microwave peak frequency f _m . Proton fluxes with hard (flat) energy spectra (δ ≤ 1.5) correspond to hard microwave frequency spectra (S ₉/S ₁₅ ≤ 1 and f _m ≥ 15 GHz), while flares with soft radio spectra (S ₉/S ₁₅ ≥ 1.5 and f _m ≤ 5 GHz) result in proton fluxes with soft (steep) energy spectra (δ ≥ 1.5–2). It is also shown that powerful high-frequency bursts with the hardest radio spectra (f _m ≈ 30 GHz) can point at acceleration of significant proton fluxes in flares occurring in strong magnetic fields. These results argue that solar cosmic rays (or at least their initial impulses) are mainly accelerated in flares associated with impulsive and post-eruptive energy release, rather than in shocks driven by coronal mass ejections.     

Refinement of the parameters of the binary pulsar B0655+64 based on 111 MHz observations

An analysis of monitoring observations for the pulsar PSR B0655+64, which is located in a binary system, at 111 MHz during 2002–2015 are presented. The Keplerian parameters of the pulsar have been refived: the longitude of periastron ω = 276_.^°5785 ± 0_.^°0005 and the orbital semi-major axis is a_p sin i = 4.124976± 0.000003 s. The parameters of the perturbed motion have been determined: the motion of periastron ω = 0_.^°315 ± 0_.^°005/ year, and the derivative of the period of the binary system ? = (-1.66 ± 0.11) × 10^-14 s/s = (-0.524 ± 0.038) µs/year. The estimated time scale for the decay of the PSR 0655+64 system is (1.7 ± 0.1) × 10¹¹ yrs.     

A “tsunami” in the protoplanetary disk of SV Cep

The results of ~15 years of photometric observations of the UX Ori star SV Cep in the near-infrared (JHKL) are presented. They demonstrate the presence of a cyclic component with a period of ~7 years in the variations of the IR fluxes. This is clearly seen in all four IR bands, but is absent in the optical. The variation amplitude is highest in the K band: ΔK ≈ 0.68 ^m . The shape of the variations differs slightly in the transition from J to L. However, it is reproduced with good accuracy during two cycles, suggesting a periodic process is observed. If the periodic perturbations in the circumstellar disk of SV Cep are due to a companion’s orbitalmotion, the orbital semi-major axis should be ~5AU, foramass of SVCep of 2.6M _⊙. The absence of a seven-year period in the optical light curve of SV Cep means that the observed period cannot be due to variations in the circumstellar extinction. The IR brightness variations could be due to the companion’s motion along an eccentric orbit, resulting in a periodic modulation of the rate of accretion onto the star.     

Analysis of Photometric Observations of the New Cataclysmic Variable ASASSN-13cx

Photometric observations of the variable star ASASSN-13cx acquired in the course of a program of studies of cataclysmic variables and their parameters recently carried out at the Sternberg Astronomical Institute (SAI) are presented. The star was observed with the 50-cm and 60-cm telescopes of the SAI Crimean Astronomical Station and a CCD photometer (～1800 images in the V and Rc filters) during the variable’s outburst of August–September 2014 and in a period of quiescence in October–November 2016. The ASASSN-13cx system is confirmed to be a SU UMa variable. Parameters of the system are derived from eight light curves using a “composite” model that takes into account the presence of a hot spot on the lateral surface of the geometrically thick disk and of a region of enhanced energy release near the disk edge, at the base of the gas flow (the so-called “hot line”). Parameters of the system for three light curves during the outburst were obtained in the framework of a “spiral” model that additionally takes into account the presence of geometric perturbations on the accretion-disk surface. The parameters of ASASSN-13cx determined using these models provide good accuracy in reproducing the system’s light curves in both states. The basic parameters of the system have been determined for the first time: the component mass ratio q = M₁/M₂ = 7.0 ± 0.2, the orbital inclination i = 79.9°?80.1°, the distance between the components’ centers of mass a₀ = 0.821(1) R_?, and the sizes and temperatures of the stars: R₁ = 0.0124(5)a₀ = 0.0102(4) R_?, T₁ = 12 500 ± 280 K, 〈R₂〉 = 0.236(4)a₀ = 0.194(3) R_?, T₂ = 2550 ± 400 K, corresponding to M4–9V for the spectral type of the secondary. Parameters of the accretion disk have been derived for both activity states. The mass of matter in the accretion disk increased by almost a factor of two during ～400 orbital periods in quiescence.     

Unusual Photometric Activity of the Weak-line T Tauti Star V715 Per

The results of many-year photometric observations of the weak-line T Tauri star V715 Per in the optical (V RI) and infrared (JHK) photometric bands are presented. The period of the star’s lowamplitude brightness variations in the V RI bands detected earlier, 5.23d, is confirmed. This period persists during the entire 14-year time interval of the observations, but its amplitude varies with time. The same period was found in the variations of JHK infrared fluxes, with the brightness variations in the two ranges being almost synchronous. The most probable origin of these variations is periodic eclipses of the star by its own disk, perturbed by its interaction with the magnetosphere (AA Tau-type variations). In addition to the periodic brightness variations, the star also displayed sporadic Algol-like minima characteristic of UX Ori stars during several years (2005–2011). The amplitude of these variations increasedmonotonically, reaching approximately 1^m in the V band by the end of 2010, after which the deep fadings ended. During this entire time interval, the amplitude of the periodic variations was observed to decrease. After 2011, the periodic component began to become stronger again. These changes of the photometric activity of V715 Per can be explained by an increasing rate of gas accretion onto the star, decreasing the radius of warped regions of the disk and leading to accompanying Algol-like fadings. This model assumes a low inclination of the circumstellar disk in V715 Per to the line of sight.     

Photometric elements and apsidal motion of the eclipsing binary ES Lac

We have obtained high-accuracy photoelectric measurements of ES Lac, an eclipsing binary with an elliptical orbit (B9III + B9III; P = 4.459^d, e = 0.198) in 1985–2004 at the Sternberg Astronomical Institute’s Tien Shan High-Altitude Observatory. Our detailed analysis of the 19-year uniform series of measurements has yielded the first photometric elements for this system, as well as a self-consistent set of physical and geometrical parameters for the binary. The virtually identical components (M ₁ = M ₂ = 3.0 M _⊙; R ₁ = R ₂ = 4.12 R _⊙) are appreciably separated from the main sequence, and are located on the giant branch: their age is t = (3.5 ± 0.2) × 10⁸ yrs. An analysis of our observations together with previously published times of minima has enabled a considerable refinement of the period of the apsidal motion, U = 355 ± 20 years, and a first determination of the apsidal parameter reflecting the radial density distributions for the components stars: k ₂ ^obs = 0.00213(18). This value is in a good agreement with the value expected theoretically for current evolutionary models of such stars: k ₂ ^th = 0.00257(15).     

First Data on the Age and <Emphasis Type="Italic">P</Emphasis>–<Emphasis Type="Italic">T</Emphasis> Formation Conditions of Zones of Low-Angle Dip Schistosity of the Belomorian Mobile Belt

The low-angle dip schistosity zones of the Belomorian mobile belt of northern Karelia are zones of plastic flow of thrust origin. They were formed from 1.85 to 1.90 Ga: 1879 ± 21 Ma according to ⁴⁰Ar/³⁹Ar for amphibole from amphibolites and 1857 ± 13 Ma according to the Sm–Nd isochron in amphibolites. The P–T parameters of rock metamorphism in low-angle dip schistosity zones correspond to the boundary of amphibolite and granulite facies of metamorphism: T = 640–765°C, rarely rising to 826°C; P = 8.0–11.7 kbar. The hypothesis of the two-stage Paleoproterozoic metamorphism of rocks of the Belomorian mobile belt was introduced.     

H-bonding scheme and cation partitioning in axinite: a single-crystal neutron diffraction and Mössbauer spectroscopic study

The crystal chemistry of a ferroaxinite from Colebrook Hill, Rosebery district, Tasmania, Australia, was investigated by electron microprobe analysis in wavelength-dispersive mode, inductively coupled plasma–atomic emission spectroscopy (ICP–AES), ⁵⁷Fe Mössbauer spectroscopy and single-crystal neutron diffraction at 293 K. The chemical formula obtained on the basis of the ICP–AES data is the following: \( ^{X1,X2} {\text{Ca}}_{4.03} \,^{Y} \left( {{\text{Mn}}_{0.42} {\text{Mg}}_{0.23} {\text{Fe}}^{2 + }_{1.39} } \right)_{\varSigma 2.04} \,^{Z1,Z2} \left( {{\text{Fe}}^{3 + }_{0.15} {\text{Al}}_{3.55} {\text{Ti}}_{0.12} } \right)_{\varSigma 3.82} \,^{T1,T2,T3,T4} \left( {{\text{Ti}}_{0.03} {\text{Si}}_{7.97} } \right)_{\varSigma 8} \,^{T5} {\text{B}}_{1.96} {\text{O}}_{30} \left( {\text{OH}} \right)_{2.18} \). The ⁵⁷Fe Mössbauer spectrum shows unambiguously the occurrence of Fe²⁺ and Fe³⁺ in octahedral coordination only, with Fe²⁺/Fe³⁺ = 9:1. The neutron structure refinement provides a structure model in general agreement with the previous experimental findings: the tetrahedral T1, T2, T3 and T4 sites are fully occupied by Si, whereas the T5 site is fully occupied by B, with no evidence of Si at the T5, or Al or Fe³⁺ at the T1–T5 sites. The structural and chemical data of this study suggest that the amount of B in ferroaxinite is that expected from the ideal stoichiometry: 2 a.p.f.u. (for 32 O). The atomic distribution among the X1, X2, Y, Z1 and Z2 sites obtained by neutron structure refinement is in good agreement with that based on the ICP–AES data. For the first time, an unambiguous localization of the H site is obtained, which forms a hydroxyl group with the oxygen atom at the O16 site as donor. The H-bonding scheme in axinite structure is now fully described: the O16–H distance (corrected for riding motion effect) is 0.991(1) Å and an asymmetric bifurcated bonding configuration occurs, with O5 and O13 as acceptors [i.e. with O16···O5 = 3.096(1) Å, H···O5 = 2.450(1) Å and O16–H···O5 = 123.9(1)°; O16···O13 = 2.777(1) Å, H···O13 = 1.914(1) Å and O16–H···O13 = 146.9(1)°].