On the influence of gradients in the angular velocity on the solar meridional motions

If fluctuations in the density are neglected, the large-scale, axisymmetric azimuthal momentum equation for the solar convection zone (SCZ) contains only the velocity correlations and where u are the turbulent convective velocities and the brackets denote a large-scale average. The angular velocity, , and meridional motions are expanded in Legendre polynomials and in these expansions only the two leading terms are retained (for example, where is the polar angle). Per hemisphere, the meridional circulation is, in consequence, the superposition of two flows, characterized by one, and two cells in latitude respectively. Two equations can be derived from the azimuthal momentum equation. The first one expresses the conservation of angular momentum and essentially determines the stream function of the one-cell flow in terms of : the convective motions feed angular momentum to the inner regions of the SCZ and in the steady state a meridional flow must be present to remove this angular momentum. The second equation contains also the integral indicative of a transport of angular momentum towards the equator.With the help of a formalism developed earlier we evaluate, for solid body rotation, the velocity correlations and for several values of an arbitrary parameter, D, left unspecified by the theory. The most striking result of these calculations is the increase of with D. Next we calculate the turbulent viscosity coefficients defined by whereC _ro ⁰ and C _o ⁰ are the velocity correlations for solid body rotation. In these calculations it was assumed that ₂ was a linear function of r. The arbitrary parameter D was chosen so that the meridional flow vanishes at the surface for the rotation laws specified below. The coefficients v _ro ⁱ and v _0o ⁱ that allow for the calculation of C _ro and C _0o for any specified rotation law (with the proviso that ₂ be linear) are the turbulent viscosity coefficients. These coefficients comply well with intuitive expectations: v _ro ¹ and –v _0o ³ are the largest in each group, and v _0o ³ is negative.The equations for the meridional flow were first solved with ₀ and ₂ two linear functions of r ( ₀ ¹ = – 2 × 10 ^–12 cm ^–1) and ( ₂ ¹ = – 6 × 10 ¹² cm ^–1). The corresponding angular velocity increases slightly inwards at the poles and decreases at the equator in broad agreement with heliosismic observations. The computed meridional motions are far too large ( 150m s^–1). Reasonable values for the meridional motions can only be obtained if _o (and in consequence ), increase sharply with depth below the surface. The calculated meridional motion at the surface consists of a weak equatorward flow for gq < 29° and of a stronger poleward flow for > 29°.In the Sun, the Taylor-Proudman balance (the Coriolis force is balanced by the pressure gradient), must be altered to include the buoyancy force. The consequences of this modification are far reaching: is not required, now, to be constant along cylinders. Instead, the latitudinal dependence of the superadiabatic gradient is determined by the rotation law. For the above rotation laws, the corresponding latitudinal variations of the convective flux are of the order of 7% in the lower SCZ.     

3He enrichments by magnetosonic waves in current-driven instabilities

A theoretical model for³He enrichments in solar energetic particles is developed. First, current-driven, electrostatic instabilities that have frequencies ( is the cyclotron frequency of³He) are investigated for a plasma consisting of H,⁴He,³He, and electrons with the density of³He much lower than those of H and⁴He. It is found that in many cases the oblique ion-acoustic waves can have positive growth rates at frequencies and, at the same time, negative growth rates at and at _H. This can occur near the marginal state of the instability. The wave damping at these frequencies is caused by the cyclotron resonances of⁴He and H. The cyclotron damping at is negligible, however, because the abundance of³He is very small. The H cyclotron waves can be unstable at for a wide region of plasma parameters; the electron-to-ion temperature ratio must beT _e /T _H 1.5. To destabilize the₄He cyclotron waves with , high⁴He density and high electron temperature are both required. Then,³He enrichments are studied on the basis of the theory of nonlinear magnetosonic waves, which can promptly accelerate ions. The current-driven electrostatic waves with can enhance fluctuation velocities of³He. Thus, in the presence of these waves, magnetosonic waves can selectively accelerate³He particles to high energies. Finally, cyclotron resonances of heavy ions with the waves or are briefly discussed.     

Simulations of cosmic-ray particle diffusion

The diffusion of charged particles in a stochastic magnetic field (strengthB) which is superimposed on a uniform magnetic fieldB ₀ k is studied. A slab model of the stochastic magnetic field is used. Many particles were released into different realizations of the magnetic field and their subsequent displacements z in the direction of the uniform magnetic field numerically computed. The particle trajectories were calculated over periods of many particle scattering times. The ensemble average was then used to find the parallel diffusion coefficient . The simulations were performed for several types of stochastic magnetic fields and for a wide range of particle gyro-radius and the parameterB/B ₀. The calculations have shown that the theory of charged particle diffusion is a good approximation even when the stochastic magnetic field is of the same strength as the uniform magnetic field.     

The spherically-symmetrical trajectories of spin particles in the Schwarzschild field

The motion of spin particles in the Schwarzschild field is examined in this paper. It is shown that Mathisson-Papapetrou equations under additional conditions , where is the Lie derivative of the Killing vector _j , permit only radial motion, motion in the equatorial plane and in the equilibrium points. All these types of motion are considered more in detail.     

Identification of Emissions in the Spectrum of Comet C/1989 VI (Skoritchenko–George) Obtained with the 6-m SAO RAS Reflector

The emission spectrum of comet Skoritchenko–George (C/1989 VI), unusual in its information content, was obtained on February 26.7 UT, 1990, with the use of a TV scanner installed on the 6-m BTA reflector of the Special Astronomical Observatory of the Russian Academy of Sciences (SAO RAS) in Nizhni Arkhyz. Detailed identification of the emission lines of this comet was made. The observed spectrum contains 311 emission lines, including those of the molecules. Among others, the lines of the negative carbon C ₂ ^- ion and the lines corresponding to the electron transition in the neutral CO molecule are discovered. The presence of a large number of lines of the neutral CO molecule (the Asundi bands and the triplet bands) in the visible region is one of the uncommon features of the emission spectrum of this comet. The triplet lines : 15–3, 13–2, 11–2, 9–1, 8–1, 7–1, 7–0, 5–0, 4–0; : 7–0, 6–0, 5–0; and a" : 11–1 (K = 3, 4); 16–4 (K= 0, 1, 2, 4); 9-0 (K= 0, 1, 2); 8–0 (K= 0) were identified for the first time. Prior to this work, the lines of CO in the visible range were observed only in the spectrum of comet C/1979 VI (Bradfield) in 1989.     

Meridional flow in thick accretion disks

We have studied the effect of the flow in the accretion disk. The specific angular momentum of the disk is assumed to be constant and the polytropic relation is used. We have solved the structure of the disk and the flow patterns of the irrotational perfect fluid.As far as the obtained results are concerned, the flow does not affect the shape of the configuration in the bulk of the disk, although the flow velocity reaches even a half of the sound velocity at the inner edge of the disk. Therefore, in order to study accretion disk models with the moderate mass accretion rate—i.e.,

[CII] Emission and Star Formation in Late-Type Galaxies

A linear correlation between the ratio of the[CII( $^{\text{2}}$ P A linear correlation between the ratio of the[CII( P → P )] line intensity to the [ CO(J:1 →0)] line emission, I /I and the equivalent width (EW) is found, over the range 2–71 ? in EW, for a sample of 21late-Type= galaxies. The latter is comprised of an optically selected sample of 12 normal Virgo Cluster spiral galaxies with [CII] detections obtained by us with ISOLWS, plus nine late-Type= galaxies with higher star formation rates (SFRs), for which [CII] data and, especially, EW data are available in the literature. As a result we infer I /I to be a reliable tracer of the current mass-normalized global SFR for non-starburst spiral galaxies. Moreover, the ratio of the [CII] line to the total far-infrared (FIR) continuum intensity, I /I , is found to decrease from ∼0.5% to ∼0.1% with decreasing SFR which we propose is due to a `[CII]-quiet' component of I from dust heated by the general interstellar radiation field (ISRF). The more `quiescent' galaxies in the sample have values of I /I different from those observed in `compact' Galactic interstellar regions. Their [CII]-emission is interpreted to be dominated by diffuse regions of the interstellar medium (ISM). For normal `star-forming' galaxies the diffuse component of the [CII] emission is estimated to account for at least 50% of the total. This revised version was published online in July 2006 with corrections to the Cover Date.     

Possible semi-annual variation of the Newtonian constant of gravitation

A possible semi-annual variation of the Newtonian constant of gravitationG is established. For the aphelion and perihelion points of the Earth's orbit we find, respectively,

Models of galactic mass distribution

Three groups of galactic mass models, each consisting of nine inhomogeneous spheroids of two kinds are described, according to three adopted values of the total density near the Sun: 0.10, 0.15 and 0.20 M pc^–3. Approximately 20% of the total mass of each model is in the halo, constructed to adequately fit recent RR Lyrae star observations. It is shown that the maxima found in the RR Lyrae star densities towards the galactic axis (Plaut, 1970) should not be interpreted as being associated with the galactic nucleus, but as the result of the greater decrease in density with increasingz over the increase in density as the galactic axis is approached. Even at the low galactic latitude of 5° (l=0°), this effect causes a 0.5 kpc correction to the distance to the galactic centre. A basic model for kpc, km s^–1, M pc^–3 is first constructed, mainly to satisfy structural conditions near the sun and in the halo. An attempt to optimize the basic model is made by scaling it so as to retain constant density and angular velocity near the sun, and to best fit kinematic data, including the recent re-examination of the 21-cm data of Simonson and Mader (1972). No unknown matter is required in the models, in accordance with the results of Weistrop (1972b), and, as pointed out earlier (Innanen, 1966b) the faintM-stars must be in a highly flattened spheroid. The optimizing indicates that an adequate fit to kinematics can be achieved for km s^–1. More detailed results are tabulated for a representative model for which . Two new galactic density functions are discussed in the Appendix.     

Gas,dust and molecules in the Galaxy

The relative abundances of cool neutral hydrogen, carbon monoxide and formaldehyde are studied using all the available observational data in the literature. The obtained mean valuesN _{H 1}/ ,N _{H 1}/N _CO,N _CO/ are approximately constant in the dark clouds of the solar neighbourhood and in the distant molecular clouds.The observed correlationsN _CO,A _v and ,A _v show that formaldehyde can also be used as an indicator of molecular hydrogen. The ratioN H₁/A _v depends on densities and decays considerably in the ranges of visual absorptions in which the molecules become detectable (A _v 2 mg); an average of /N _{H 1}10 is calculated for the dense dark clouds.Indications of systematic temperature gradiens T/A _v are found for formaldehyde and neutral hydrogen inside the dark clouds, and qualitative comparisons are made with theoretical quantum mechanics calculations.The observed carbon monoxide and formaldehyde abundances, the free electron layer in the Galaxy, the distribution of neutral hydrogen in different states are only compatible if an ionization rate of 10^–16 is accepted, provided presumably by 2 MeV protons of cosmic radiation.Three main states for neutral hydrogen and dust are identified from different kinds of observational data (21 cm line in emission, absorption in galactic radio sources and self-absorption in the hot gas background): (1) a homogeneous intercloud stratum of tenuous gas and dust with a galactic halfwidth of 350 pc and mean parametersn _H=0.2 atom cm^–3, spin temperatureT _s 10000 K andn _d 0.3 mg kpc^–1; (2) cool gas and dust concentrated in spiral features with a galactic half-width of less than 100 pc, probably forming clouds with diffuse and indefinite limits, with mean parametersn _H2 atom cm^–3,T _s <1100 K (probable average,T _s =135 K) andn _d 3 mg kpc^–1; (3) dense gas and dust clouds with a mean diameter of 7 pc and mean parametersn _H700 atom cm^–3 (90% in a molecular state),T _s 63 K andn _d 1 mg pc^–1 on which molecules as CO and H₂CO are formed.The application of the Jeans criteria for gravitational instability shows that the dense clouds are gravitationally bound while the gas in the intermediate state (2) can be protected against collapse by the total internal energy in the medium increasing due to cosmic rays and the magnetic field in the Galaxy.The observed velocity halfwidths and galacticZ-halfwidths in states (1) and (2) are compatible with a total mass density in the galactic layer of 90M pc^–2 (gas plus stars) according to the barometric equation.The relative abundancesN _{H 1}/N _CO, calculated from C¹²O and C¹³O data and comparisons with studies in the 21 cm emission line, show that the antenna temperatureT _A ⁺ in the 2.6 mm line of C¹²O is a good indicator of the cool gas densities in the Galaxy. The possible application of this for studies in galactic structure is discussed and hypothetical distributions of carbon monoxide in the zones outside the galactic planeB=0° are presented.From a synthesis based on the results obtained, a cycle is postulated for the neutral hydrogen in the Galaxy: condensation and cooling of gas molecular formation gravitational collapse and star formation gas dissipation and heating by cosmic rays and UV radiation.     

The critical and the saturation content of magnetic monopoles in rotating relativistic objects

Both the critical content _c ( N _m /N _B , whereN _m ,N _B are the total numbers of monopoles and nucleons, respectively, contained in the object), and the saturation content _s of monopoles in a rotating relativistic object are found in this paper. The results are:

Persistence of solar activity on small scales: Hurst analysis of time series coming from Hα flares

We calculated the Hurst exponent H for the daily averaged intensity Q of optical flares, an index which describes the solar activity. We found that H0.74±0.02 in the range of scales from about 20 days up to 450 days. This value is well beyond H= , expected for a stochastic Brownian process, thus indicating that the solar cycle could show persistence on small scales, in agreement with what has been found using other indices of the solar cycle.     

Effect of oblateness on the non-linear stability of L 4 in the restricted three-body problem

Non-linear stability of the libration point L ₄ of the restricted three-body problem is studied when the more massive primary is an oblate spheroid with its equatorial plane coincident with the plane of motion, Moser's conditions are utilised in this study by employing the iterative scheme of Henrard for transforming the Hamiltonian to the Birkhoff's normal form with the help of double D'Alembert's series. It is found that L ₄ is stable for all mass ratios in the range of linear stability except for the three mass ratios: $$\begin{gathered} \mu _{c1} = 0.0242{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.1790{\text{ }}...{\text{ }}A_1 , \hfill \\ \mu _{c2} = 0.0135{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.0993{\text{ }}...{\text{ }}A_1 , \hfill \\ \mu _{c3} = 0.0109{\text{ }}...{\text{ }}{}^{{\text{\_\_}}}0.0294{\text{ }}...{\text{ }}A_1 . \hfill \\ \end{gathered} $$      

Tidal effects in the earth–moon system and the earth's rotation

Analytical expressions for tidal torques induced by a tide‐arising planet which perturbs rotation of a nonrigid body are derived. Corresponding expressions both for secular and periodic perturbations of the Euler's angles are given for the case of the earth's rotation. Centennial secular rates of the nutation angle θ and of the earth's angular velocity ω, as well as the centennial logarithmic decrement ν of the Chandler wobble are evaluated:  mas, . In the Universal Time (UT) a large out‐of‐phase (sine) dissipative term with the period 18.6 years and the amplitude 2.3 ms is found. Corrections to nutation coefficients, which presumably have not been taken into account in IAU theory, are given. This revised version was published online in July 2006 with corrections to the Cover Date.     

On numerical evaluation of theH-functions of transport problems by kernel approximation for the albedo 0<ω≤1

Das Gupta represented theH-functions of transport problems for the albedo [0, 1] in the formH(z)=R(z)–S(z) (see Das Gupta, 1977) whereR(z) is a rational function ofz andS(z) is regular on [–1, 0] ^c . In this paper we have representedS(z) through a Fredholm integral equation of the second kind with a symmetric real kernelL(y, z) as . The problem is then solved as an eigenvalue problem. The kernel is converted into a degenerate kernel through finite Taylor's expansion and the integral equation forS(z) takes the form: (which is solved by the usual procedure) where _r 's are the discrete eigenvalues andF _r 's the corresponding eigenfunctions of the real symmetric kernelL(y, z).     

POWER SPECTRA OF SOLAR NETWORK AND NON-NETWORK FIELDS

We report new properties of solar magnetic fields in a quiet region as found from their magnetic power spectra. The power spectra of network and intranetwork fields (non-network fields) are separately calculated from a Big Bear magnetogram obtained with moderately high spatial resolution of 1.5 arc sec and a high sensitivity reaching 2 Mx cm^-2. The effect of seeing on the power spectrum has been corrected using Fried's (1966) Modulation Transfer Function with the seeing parameter determined in our previous analysis of the magnetogram. As a result, it is found that the two-dimensional power spectra of network and non-network fields appear in a form: ( ₁) ^-1 and ( ₁) ^-3.5. Here ₀ 0.47 Mm^-1 for network fields and ₀ 0.69 Mm^-1 for non-network fields, the latter of which corresponds to the size of mesogranulation; ₁ 3.0 Mm^-1 for both, which is about the size of a large granule. The network field spectrum below ₀ appears nearly flat, whereas that of non-network fields instead decreases towards lower wave numbers as ( ) ^1.3. The turnover behavior of magnetic field spectra around ₁ coincides with that found for the velocity power spectrum, which may justify the kinetic approach taken in previous theoretical studies of the solar magnetic power spectra.     

Earth’s precession–nutation motion: the error analysis of the theories IAU 2000 and IAU 2006 applying the VLBI data of the years 1984–2006

The long-term systematic errors of the analytical theories IAU 2000 and IAU 2006 of the Earth’s precession–nutational motion are studied making use of the VLBI data of 1984–2007. Several independent methods give indubitable evidence of the significant quadratic error in the IAU 2000 residuals of the precessional angle while the adopted value of the secular decrease /cy of the Earth’s ellipticity e (derived from Satellite Laser Ranging data) should manifest itself in the residuals of as the negative quadratic trend . The problem with the precession of the IAU 2006 theory adopted as a new international standard and based on the precession model P03 (Capitaine et al., Astron Astrophys 432:355–367, 2005) appears to be even more serious because the above mentioned quadratic term has already been incorporated into the P03 precession. Our analysis of the VLBI data demonstrates that the quadratic trend of the IAU 2006 residuals does amount to the expected value (30.0 ± 3) mas/cy². It means, first, that the theoretical precession rate of IAU 2006 should be augmented by the large secular correction and, second, that the available VLBI data have potentiality of estimating the rate . And indeed, processing these data by the numerical theory ERA of the Earth’s rotation (Krasinsky, Celest Mech Dyn Astron 96:169–217, 2006, Krasinsky and Vasilyev, Celest Mech Dyn Astron 96:219–237, 2006) yields the estimate /cy statistically in accordance with the satellite-based . On the other hand, applying IAU 2000/2006 models, the positive value /cy is found which is incompatible with the SLR estimate and, evidently, has no physical meaning. The large and steadily increasing error of the precession motion of the IAU 2006 theory makes the task of replacing IAU 2006 by a more accurate model be most pressing.     

Determination of Optical Spectra by a Modified Kramers Kronig Integral

A new modified Kramers Kronig Integral is derived and shown to produce excellent results when k data is only known over a limited range. By considering the effect of resonance features simulated using the Dirac-Delta function, the new integral is shown to be more rapidly converging than both the conventional Kramers Kronig integral and a modified (Subtractive Kramers Kronig – SKK) integral introduced by Ahrenkiel (1971). The new integral does not require extensive extrapolation of reflectance data outside the measured region in order to produce reliable results. By extending the above procedure to include n data points, it is shown that at wavelength ₀, \[ n(_0)=\sum_{i=1}^{\rm n}(-1)^{\rm n+1}\prod_{\stackrel{j=1}{j \not=i}}^{\rm n} \frac{(_j^2-_0^2)}{(_i^2- _j^2)}n(_i)+\frac{2}{\pi}P\int_{0}^{\infty}(-1)^{\rm n+1} \frac{\prod_{i=1}^{\rm n}(_i^2-_0^2)}{\prod_{i=0}^{\rm n}(^2-_i^2)} k()d \] with relative error given by, \[ R_n(_0)=\prod_{i=1}^{\rm n}\frac{_i^2- _0^2}{_^2-_i^2} . \] This n^th order expression should prove useful in establishing the internal self-consistency of data sets for which both optical coefficients have been theoretically derived.     

On the role of plasma effects in the cosmic ray propagation and isotropization in the Galaxy

Cosmic ray (c. r.) propagation in interstellar magnetic fields is often considered in the diffusion approximation, i.e. by the diffusion equation in the coordinate space. Cosmic ray momentum distribution in this case is considered isotropic when the space gradients of c.r density are absent. This approach, with the use of an unfixed effective diffusion coefficientD independent of the energyE enables one to describe all the data available However, neither the diffusion mechanism nor the limits of applicability of the diffusion approximation is clear particularly ifD is independent ofE. Furthermore, the diffusion coefficientD must be expressed through the characteristics of the interstellar medium and possibly through the flux velocity and density of c.r. etc. One of the possible approaches for the analysis of the mechanism and characteristic features of c.r. distribution and isotropization is the account taken of the plasma effects and specifically, the study of c.r. flux instability arising when c.r. are moving in the interstellar plasma. As a result of such instability c.r. may generate waves of different types (magnetohydrodynamic, high-frequency plasma and other waves). Generation of waves and scattering on them result in isotropization of cosmic rays while their propagation under certain conditions turns out similar to that under diffusion.An attempt is made here to systematically analyse the avove mentioned plasma effects and to find out to what extent they are responsible for the behaviour of c.r. in the Galaxy. It turns out that c.r. In any case this is true if this mechanism is regarded as the only c.r. isotropization mechanizm within a wide energy range from 1 to 1000 GeV. Isotropization and spatial diffusion of c.r. up toE100–1000 GeV on the waves from external sources (for example, on the waves from the supernova shells) also proved impossible if the diffusion coefficient is assumed to be independent of c.r. energy. Some new possibilities of c.r. isotropization are also considered.A List of Notations D cosmic ray (c.r.) space diffusion coefficient - degree of c.r. anyisotropy - E,E _kin total and kinetic particle energy - p,p particle momentum and its absolute value - angle between the particle momentum direction and the magnetic field direction (z-axis) - cos - v, particle velocity and its absolute value - c light velocity - f(p),f(E) momentum and energy particle distribution function - N( > E) = N( > p) = f(p) dp/(2)³ = _E f dE c.r. particle density - c.r. spectrum index,N(>E)=KE ^–+1 - n _H neutral particle density - n=n _e=n _i ion and electron density - H niagnetic field - T temperature - thermal velocities of electrons and ions - Boltzmann constant - Alfén velocity - M, m proton and electron masses - e electron charge - wave frequency - _H =eH/Mc, = _H (Mc ²/E) gyrofrequency of a plasma proton and relativistic particle - _H =eH/mc gyrofrequency of an electron - plasma frequency - v _ii,v _ei,v _en,v _in collision frequencies between ions, electrons and ions, electrons and neutrals, ions and neutrals - growth rate of wave amplitude - k,k wave vector and its absolute value - angle between the directions of the vectorsk andH - wave energy density     

A few observations of and remarks on V1500 Cygni

The development of the post-nova light curve of V1500 Cyg inUBV andHβ, for 15 nights in September and October 1975 are presented. We confirm previous reports that superimposed on the steady decline of the light curve are small amplitude cyclic variations. The times of maxima and minima are determined. These together with other published values yield the following ephemerides from JD 2 442 661 to JD 2 442 674: $$\begin{gathered} {\text{From}} 17 {\text{points:}} {\text{JD}}_{ \odot \min } = 2 442 661.4881 + 0_{^. }^{\text{d}} 140 91{\text{n}} \hfill \\ \pm 0.0027 \pm 0.000 05 \hfill \\ {\text{From}} 15 {\text{points:}} {\text{JD}}_{ \odot \max } = 2 442 661.5480 + 0_{^. }^{\text{d}} 140 89{\text{n}} \hfill \\ \pm 0.0046 \pm 0.0001 \hfill \\ \end{gathered} $$ with standard errors of the fits of ±0 _. ^d 0052 for the minima and ±0 _. ^d 0091 for the maxima. Assuming V1500 Cyg is similar to novae in M31, we foundr=750 pc and a pre-nova absolute photographic magnitude greater than 9.68.