Numerical Integration of Regular Equations of the Planar Motion of Terrestrial Bodies

The parameters of L matrices are applied to the numerical integration of regular equations describing the motion of minor bodies in the Solar System. The problem of the optimal choice of the regularizing change of variables is formulated in the context of the numerical integration of the equations of motion using the Runge–Kutta–Fehlberg method. Arbitrary perturbations are taken into account. This problem is completely solved in the case of planar motion. The solution of the optimization problem reduces the amount of computations needed to determine the vector of perturbing accelerations. Results of numerical integrations are given.     

The effect of gradient and curvature drifts on Forbush decreases

The average profile of Forbush decreases, produced by eastern-, central- and western-region solar flares is obtained separately by superposed epoch analysis for the periods 1966–1969 (qA < 0) and 1971–1979 (qA > 0). It is observed that the recovery of an average Forbush decrease from the maximum depression level is faster for the situation qA > 0 than for the situation qA < 0. This is in accordance with expectations from the drift theory. It is also observed that the drift effect is more pronounced for western-flare Forbush decreases which, of course, have a smaller magnitude compared to eastern- and central-flare Forbush decreases.The average profiles of simple and complex type Forbush decreases are also obtained separately for three periods 1965–1979, 1971–1979, and 1981–1987. It is found that the average profiles of simple and complex type Forbush decreases observed during the period 1965–1969 and 1971– 1979 are quite in agreement with drift theory. The anomalous behavior of average Forbush-decrease profiles during the period 1981–1987, especially in simple type Forbush decreases, is also explained by a drift current sheet tilt model.     

Periodic motions around pulsating stars

The periodic motion of a test particle (dust, grain, or a larger body) around a pulsating star with a luminosity oscillation of small amplitude (featured by a small parameterB) is being studied. The perturbations of all orbital elements are determined to first order inB, by using Delaunay-type canonical variables and a method whose bases were put forth by von Zeipel. According to the value of the ratio oscillation frequency/dynamic frequency, three possible situations are pointed out: nonresonant (NR), quasi-resonant (QR), and resonant (R). The solution of motion equations shows that only in the (QR) and (R) cases there are orbital parameters (argument of periastron and mean anomaly) affected by secular perturbations. These solutions (which indicate a secularly stable motion in a first approximation) are valid over prediction times of orderB ^–1 in the (NR) case andB ^–1/2 in the (QR) and (R) cases. The theory may be applied to various astronomical situations.     

Light curves of V711 Tauri (HR 1099)

V andB observations of V711 Tauri (HR 1099) for the year 1982–1983 are presented. The light curves made with the data of this season are compared with those of the previous observations of the 1981–1982 season. The single-peaked maximum light had shifted slightly and the peak amplitude had decreased to 0.04 mag. inV. No clear sign of colour,B–V, variation in comparison with the low temperature dark-spot is detected.Paper presented at the IAU Third Asian-Pacific Regional Meeting, held in Kyoto, Japan, between 30 September–6 October, 1984.     

Time distribution of flare events on UV Cet stars

A statistical analysis of the flare data on three UV Cet stars obtained at Catania Observatory in the indicated periods-UV Cet (1968–1976), EQ Peg (1969–1975) and YZ CMi (1968–1976)-is presented. The distributions of the four variables (i) , the time elapsed from the start of a continuous observation interval to the occurrence of the first flare; (ii) , the time interval between consecutive flares; (iii)k, the number of flares occurring in time intervals of given duration; (iv) , the ratio between the number of flares observed in each continuous observing interval to that expected for a Poisson process, indicate that the flares of the three stars considered are not randomly distributed in time.     

A unified presentation of the Voigt functions

In the present paper, we have given a generalization of a unified study of the Voigt functionsK(x, y) andL(x, y) obtained by Srivastava and Miller (1987; Vol. 135, pp. 111–118) which play an important role in several diverse fields of physics-such as astrophysical spectroscopy and the theory of neutron reactions. Explicit expressions for these functions are given in terms of relatively more familiar special functions of one and two variables; indeed, each of these representations will naturally lead to various other needed properties of the Voigt functions.     

A study of non-thermal radio emission features using fine spectral BAO and high-sensitivity RATAN observations of a solar active region

Based on spectral observations of active region NOAA 8545 on 19 May 1999, we describe the processes responsible for non-thermal long-lasting radio emission and for narrow-band non-drifting bursts observed at the same time. Non-thermal long-lasting radio emission consisted of two components: short-duration (1–2 s) microbursts with fluxes about 0.001 s.f.u. and continuum emission with growing spectrum in the range of 1000–2000 MHz. Energetic electrons continuously existed in the active region for more than 2.2 hours. The nature and parameters of microbursts were discussed by Bogod, Mercier, and Yasnov (2001). Here we consider the continuum source nature. It is shown that the model, taking into account the cyclotron loss-cone instability of hot electrons and the generation of plasma waves at the upper hybrid frequency, may explain the observed continuum source parameters. For the narrow-band non-drifting bursts we consider two models: the first taking into account an excitation of weak shock waves across the magnetic field and the second with an excitation of the upper hybrid waves under the double plasma resonance. Continuum source parameters are close to the last model. Our estimations for the magnetic field strength are as follows: H=120–126 G, which is valid for the region where the electron density of background plasmas n=(1.4–1.9) ×10⁹ cm^–3; H=180–190 G for the region where n=(3.0–4.3) ×10⁹ cm^–3; H=290 G for the region where n _e=2.5×10¹⁰ cm^–3; and H=350 G for the region where n _e=3.5×10¹⁰ cm^–3. The speed of the fast electrons is about 0.10–0.14 c.     

Power spectra of the interplanetary magnetic field

Power spectra based on Pioneer 6 interplanetary magnetic field data in early 1966 exhibit a frequency dependence of f ^–2 in the range 2.8 × 10^–4 to 1.6 × 10^–2 cps for periods of both quiet and disturbed field conditions. Both the shape and power levels of these spectra are found to be due to the presence of directional discontinuities in the microstructure (< 0.01 AU) of the interplanetary magnetic field. Power spectra at lower frequencies, in the range of 2.3 × 10^–6 to 1.4 × 10^–4 cps, reflect the field macrostructure (> 0.1 AU) and exhibit a frequency dependence roughly between f ^–1 and f ^–3/2. The results are related to theories of galactic cosmic-ray modulation and are found to be consistent with recent observations of the modulation.     

Power Spectrum of Cosmic-ray Fluctuations During Consecutive Solar Minimum and Maximum Periods

Power spectral analysis of cosmic-ray intensity recorded by eight stations was carried out over a wide range of frequencies from 2.3 × 10^–8 Hz to 5.8 × 10^–6 Hz (2–500 days) during the period 1964–1995. Spectrum results of large-scale fluctuations have revealed the existence of a broad peak near 250–285 days and a narrower peak at 45–50 days during the studied epochs as a stable feature in all neutron monitors covering a wide rigidity range. The cosmic-ray power spectrum displayed significant peaks of varying amplitude with the solar rotation period (changed inversely with the particle rigidities) and its harmonics. The amplitudes of 27-day and 13.5-day fluctuations are greater during the positive-polarity epochs of the interplanetary magnetic field (qA>0) than during the qA<0 epochs. The comparison of cosmic-ray power spectra during the four successive solar activity minima have indicated that at the low-rigidity particles the spectrum differences between the qA>0 and qA<0 epochs are significantly large. Furthermore, the spectrum for even solar maximum years are higher and much harder than the odd years. There are significant differences in the individual spectra of solar maxima for different cycles.     

Methanol 72–81A+ emission and some maser series in Orion KL

Methanol 7₂–8₁ A ⁺ is mapped for the first time in Orion KL. Analysing the observed data and solving the statistical equilibrium and radiative transfer equations, it is concluded that line series ofJ ₂–(J+1)₁ A ⁺ (J=7,8,9) is in quasi-thermal emission rather than the masers in Orion KL. The maser spots of methanolJ ₂–J ₁ E (J=6,7) and 8₀–7₁ A ⁺ are distributed in the northeast part of the contour plot of 7₂–8₁ A ⁺. The physical conditions of the regions of maser seriesJ ₀–(J–1)₁ A ⁺ (J=7,8,9) are discussed. Also from the calculation results another maser seriesJ ₁–(J–1)₂ A ^– (J=10,11,12) that might coexist with maser seriesJ ₂–J ₁ E, is found. The sizes of the 2-dimension Gaussian fit plots of methanol 7₂–8₁ A ⁺ and HCOOCH₃ 10(0,10)–9(0,9)A are almost the same, and the main parts overlap each other.     

The masers of E-type methanol in orion KL and SGR B2

Using a simplified model the statistical equilibrium and radiative transfer equations of E-type-CH₃OH are solved for Orion KL and SgrB2. According to our calculation results and the observation data taken by Matsakiset al. (1980) and Morimotoet al. (1985a, b), the physical conditions of both sources are estimated. In theJ ₂-J ₁ E methanol maser region of Orion KL, the density, kinetic temperature, dust temperature, and the fractional abundance are 0.8–2×10⁶ cm^–3, 150, 30–90 K, 0.8–8×10^–6. In the 4_–1-3₀ E and 5_–1-4₀ E methanol maser region of Sgr B2 the correspondance physical conditions above are 10⁴ cm³, 45, 23 K, and 7×10^–7, respectively.     

Energy spectrum and time variations of hard X-rays from Cyg X-1

Experimental results on the intensity, energy spectrum and time variations in hard X-ray emission from Cyg X-1 based on a balloon observation made on 1971, April 6 from Hyderabad (India) are described. The average energy spectrum of Cyg X-1 in the 22–154 keV interval on 1971 April 6 is best represented by a power law dN/dE=(5.41±1.53)E ^{–(1.92±0.10)} photons cm^–2s^–1 keV^–1 which is in very good agreement with the spectrum of Cyg X-1 derived from an earlier observation made by us on 1969 April 16 in the 25–151 keV band and given by dN/dE=(3.54±2.44)E ^{–(1.89±0.22)} photons cm^–2s^–1 keV^–1. A thermal bremsstrahlung spectrum fails to give a good fit over the entire energy range for both the observations. Comparison with the observations of other investigators shows that almost all balloon experiments consistently give a spectrum of E ^–2, while below 20 keV the spectrum varies fromE ^–1.7 toE ^–5. There is some indication of a break in the Cyg X-1 spectrum around 20 keV. Spectral analysis of data in different time intervals for the 1971 April 6 flight demonstrates that while the source intensity varies over time scales of a few minutes, there is no appreciable variation in the spectral slope. Analysis of various hard X-ray observations for long term variations shows that over a period of about a week the intensity of Cyg X-1 varies upto a factor of four. The binary model proposed by Dolan is examined and the difficulties in explaining the observed features of Cyg X-1 by this model are pointed out.     

Luminosity function for galaxies and cosmology

Formulae for two classical cosmological tests,N(z) andN(m) are obtained with the luminosity function taken into account. From the testN(m) and Jarvis and Tyson observational data, characteristic absolute magnitude of galaxiesM ^*=–20.05 and mean extragalactic absorptionA=4×10^–4 mag Mpc^–1 have been found.     

The peaks and gaps in the redshift distributions of active galactic nuclei and quasars

The distributionsf(z) of the redshifts for active galaxies (Seyfert galaxies, radio galaxies, and quasars) have been studied. Some statistically-significant maxima and minima are observed in the distributionsf(z) for these objects. The significance of peaks and gaps increases for the brighter objects, for which the samples are more complete. The clustering of the Seyfert galaxies is significantly different from that of the nearby normal galaxies. The distributionf(z) for the radio galaxies is similar to the analogous distribution for the galaxy clusters. Three of the five peaks in the distributionf(z) for the radio quasars may be caused by the selection effects. Two peaks within the intervalsz (0.5, 0.7) and (1.0, 1.1) are probably real. The corresponding scales of the QSO clustering along the line-of-sight are about 100h ^–1 Mpc (h is the Hubble constant in the units of 75 km s^–1 Mpc^–1). The possibility of some global quasi-periodical cycles for the processes of activity is discussed. The period of a cycle for the Seyfert and radio galaxies is about 1×10⁸ years that corresponds to the distances of about 30h ^–1 Mpc between the shells.     

Seven-digit tables of Lane-Emden functions

In Table I we present seven digit numerical solutions of the Lane-Emden equation for the plane-parallel, cylindrical, and spherical case for polytropic indices ofn=–10, –5, –4, –3, –2, –1.5, –1.01, –0.9, –0.5, 0, 0.5, 1, 1.5, 2, 3, 4, 5, 6, 10, 20, ±, supplemented byn=2.5, 3.5, 4.5, and 4.99 for the spherical case.In Table II some finite boundary values of polytropic slabs, cylinders, and spheres are summarized. For polytropic spheres (N=3) we have also quoted boundary values near the minimum of the dimensionless mass -² ₁₁ occurring atn4.823 (Seidov and Kuzakhmedov, 1978).     

Low-frequency electric field fluctuations excited by ring current protons

Energetic protons haying ring type distributions are shown to generate low-frequency electrostatic waves, propagating nearly transverse to the geomagnetic field lines, in the ring current region by exciting Mode 1 arid Mode 2 nonresonant instabilities and a resonant instability. Mode 1 nonresonant instability has frequencies around ~4 Hz with transverse wavelengths of ~(8–80) km, and it is likely to occur in the region L = (7–8). Mode 2 nonresonant instability can generate frequencies ~(850–1450) Hz with transverse wavelengths ~(2–20) km. The typical frequencies and transverse wavelengths associated with the resonant instability are (950–1250) Hz and (30–65) km. Both the Mode 2 nonresonant instability and the resonant instability can occur in the ring current region with L = (4–6). The low-frequency modes driven by energetic protons could attain maximum saturation electric field amplitude varying from 0.8 mV/m to 70 mV/m. It is suggested that the turbulence produced by the low-frequency modes may cause pitch angle scattering of ring current protons in the region outside the plasmapause resulting in the ring current decay.     

CCD surface photometry of an outburst of P/Schwassmann-Wachmann 1 comet

We have observed an outburst of P/Schwassmann-Wachmann 1 (SW1) on 25–26 July 1987 using a CCD imager with R-band filter. The total brightness increased from 15.5 mg (25d) to 14.4 mag (26d) during tabout 24 hrs. The southward elongated coma of 25 was detected. The radial surface brightness (B) profiles are plotted against apparent distance p from the nucleus. The logarithmic derivative k = d ln B/d ln p for the inner coma is found to have steepened from k = –1.40 (25d) to k = –1.69 (26d), whereas that for the outer coma showed no appreciable change (k = –1.19 ~ –1.22). The ellipticity of the isophotal contour of the inner coma increased about 15% fro 25d to 6d. It is concluded that the scale of this outburst was smaller than the typical ones whose magnitude change is 5–8 mag. From recent findings on the outburst natur eof SW1 including ours, a working model of the nuclues is proposed.     

Multicolor photoelectric and photographic observations of the long-period variable Y Ori

Results of observations of the long-period variable Y Ori are presented. The photoelectric observations in UBVR were made at the 60-cm telescope of the high-altitude Maydanak station of the Tashkent Astronomical Institute during the autumn of 1989. The UBVR lightcurves as well as the variations in the color indices U-B, B-V, and V-R are presented. The photographic observations were made at the 40 Schmidt telescope of the Byurakan Observatory and at the 70-cm Maksutov telescope of the Abastumani Observatory. A nebulosity was discovered around Y Ori in red light near the brightness maximum. Such a formation is observed for the first time, not only for Y Ori but also for long-period variables in general. The obtained results are discussed in this work.Translated fromAstrofizika, Vol. 38, No. 1, pp. 5–15, January–March, 1995.     

Estimation of an upper limit for the solar neutron emission during large flares

Solar neutron emission during large flares is investigated by using neutron monitor data from the mountain stations Chacaltaya (Bolivia), Mina Aguilar (Argentine), Pic-du-Midi (France) and Jungfraujoch (Switzerland). Registrations from such days on which large flares appeared around the local noon time of the monitor station are superimposed with the time of the optical flare as reference point.No positive evidence for a solar neutron emission was found with this method, However, by using an extrapolation of the neutron transport functions given by Alsmiller and Boughner a rough estimation of mean upper limits for the solar neutron flux is possible. The flux limits are compared with Lingenfelter's model calculations.From the Chacaltaya measurements it follows: N ₀2.8 × 10^{–3 N} cm^–2 s^–1 per proton flare, E > 50 MeV, if P0 = 125 MV N ₀1.4 × 10^{–2 N} cm^–2 s^–1 per proton flare, E > 50 MeV, if P ₀ = 60 MV and from Pic-du-Midi measurements: N ₀6.7 × 10^{–3 N} cm^–2 s^–1 per proton flare, E > 50 MeV, if P ₀ = 125 MV N ₀4 × 10^{–2 N} cm^–2 s^–1 per proton flare, E > 50 MeV, if P ₀ = 60 MV P ₀ = characteristic rigidity of the producing proton spectrum on the Sun.The flux limits estimated for some special proton flares are consistent with Lingenfelter's predictions for the acceleration phase but are too small for the slowing down phase. Therefore it is believed that Lingenfelter's assumption of isotropic proton emission from the flare region is not fulfilled.     

An expanding circumstellar cloud of Zeta Aurigae

Strong absorption satellite lines of CaI 6572 were found on spectrograms taken on three successive days just after the fourth contact of the 1971–72 eclipse of Zeta Aurigae. The radial velocities of the satellite lines are –88 km s^–1, –74 km s^–1, and –180 km^–1, respectively, relative to the K-type primary star (K4 Ib). These absorptions should be due to a circumstellar cloud in which the column density of neutral calcium atoms is 1×10¹⁷ cm^–2 and the turbulent velocities come to 20–50 km s^–1. It is suggested that the cloud may be formed by the rocket-effect of the Lyman quanta of the B-type component (B6 V). We estimate the density in the cloud to be 2×10¹¹ atoms cm^–3 fors=10R _K and 2×10¹⁰ atoms cm^–3 fors=10² R _K, wheres denotes the distance of the cloud from the K star andR _K the K star's radius. The mass loss rate of the K-type component is also estimated to be about 10^–7 M yr^–1, assuming that the expansion of the K star occurs isotropically.