Galactic cosmic-ray anisotropy and its heliospheric modulation, inferred from the sidereal semidiurnal variations observed in the rigidity range 300–600 GV with multidirectional muon telescope at Sakashita underground station

The existence of sidereal semidiurnal variation of cosmic-ray intensity in a rigidity region 10²-10³ GV has been reported by many researchers, but there is no consensus of opinion on its origin. In this paper, using the observed semidiurnal variations in a rigidity range (300–600 GV) with 10 directional muon telescopes at Sakashita underground station (geog. lat. = 36°, long. = 138°E, DEPTH = 80 m.w.e.), the authors determine the magnitudes (η₁, η₂) and directions (a₁, a₂) of the first- and second-order anisotropies in the following galactic cosmic-ray intensity distribution (j)

jdp = j₀{1 + η₁P₁(cos χ₁) + η₂P₂(cos χ₂)}dp

, where P_nis the nth order spherical function and χ_n is the pitch angle of cosmic rays with respect to a_n. For the determination, the influence of cosmic-ray's heliomagnetospheric modulation, geomagnetic deflection and nuclear interaction with the terrestrial material and also of the geometric configuration of the telescopes are taken into account. Usually, the semidiurnal variation is produced by the second-order anisotropy. The present observation, however, requires also the first-order anisotropy which usually produces only the diurnal variation, but can produce also the semidiurnal variation as a result of the heliospheric modulation. The first- and second-order anisotropies are characterized with η₁) > 0 and η₂ < 0 have almost the same direction (a₁ a₂) specified by the right ascension ( 0.75 h) and declination (δ 50°S) and, therefore, they can be expressed, as a whole, by an axis-symmetric anisotropy of loss-cone type (i.e. deficit intensities in a cone). It is noteworthy that this anisotropy approximately coincides with that inferred from the air shower observation at Mt Norikura in the rigidity region 10⁴ GV.     

Kp index correlations with solar-wind parameters during the first and second stages of a recurrent geomagnetic storm

The dependence of geomagnetic activity during a recurrent magnetic storm on the solar-wind magnetic field and plasma parameters has been studied. According to variations of solar-wind magnetic field strength B, a recurrent magnetic storm is divided into two stages: the first proceeding during the peak of B, and the second proceeding after the return of B to quiet level. The K_p index vs solar-wind parameters scattering diagrams for stages I and II differ significantly. In particular, the random scattering for stage I is much larger than for stage II. It was found that for stage I the K_p index correlates with B, with the sign and value of northsouth field component B_z and with the magnitude ΔB of field fluctuations, the situation being similar to that during sporadic magnetic storms, though the scale of the event is smaller. For stage II, the K_p index does not correlate with B, but strongly correlates with ΔB and weaker—with B_z. So geomagnetic activity at stage II is supported mainly by solar-wind magnetic field fluctuations. The dependence of the K_p index on plasma parameters (concentration of protons n, bulk velocity v and temperature T) is weak for both stages.     

The effect of an electric field induced by a time-dependent ring current on the particle drift motion

The effect of an electric field induced by a rapidly decaying ring current on the motion of charged particles in the magnetosphere has been investigated using Euler potentials. For a model consisting of the earth dipole and the symmetric ring current, the electric field satisfies the condition E . B = 0.

Under this circumstance, the E × B drift of the particle can be identified as the motion of the magnetic field lines and vice versa. The time dependent electric field induced can be evaluated in a Spherical polar coordinate system by the formula

where and β are Euler potentials.

A model calculation on the particle drift velocity v_D = E × B/B² shows that the radial component of the drift velocity is in good agreement with those deduced from whistler duct studies.     

Effect of field-aligned currents in the polar cap at the recovery phase of a magnetic storm

Unusually great fluctuations in the ΔB module of the geomagnetic field have been observed in the polar cap from the satellite Cosmos-321. They are explained by small-scale two-sheet field-aligned current systems which exist during the periods when magnetic fields having a considerable northward B_z(B_z 10 nT) component are observed in interplanetary space.     

Electro-hydromagnetic waves in a fully ionized gas—I

Transport equations are used to determine coefficients which are generalizations for any frequency of electric field of the parallel, Pedersen and Hall conductivities in a fully ionized gas.

These coefficients are used in an investigation of the propagation of weak electromagnetic and hydromagnetic waves of all frequencies across a homogeneous and constant magnetic field in a rarefied fully ionized gas. For propagation perpendicular to the magnetic field it is found for all frequencies

(i)

(ii)

where V² = H²/4π and v, h are the perturbations of the velocity, magnetic field. Similar relationships are deduced for propagation at any angle to the field for frequencies greater than about 10 times the gyrofrequency of electrons.

The theory is applied to discuss transmission of disturbance across the interplanetary medium, the temperature of the solar corona and the earth's outer atmosphere, the emission of non-thermal solar radio noise, cosmic radio noise and the anomalous emission of light from shock fronts.     

Surface temperature of a rotating charged body

The surface temperature of a rotating, charged body is found separately under the Kerr-Newman metric and the vector graviton metric. Particular reference is made to pulsars. It is found that, 1) under the Kerr-Newman metric, the surface temperature rises from the poles to the equator, when the radius R of the body is greater than a certain critical value, r_n. When R= r_n, the surface temperature is uniform. When R < r_n, the above gradient is reversed. For pulsars, the equatorial temperature is some 3 × 10⁴ K higher than the polar temperature. 2) Under the Vector graviton field metric, a similar temperature differential exists, but it is much smaller in size.     

Deactivation of N2A Σumolecules in the aurora

Recent rocket observations of the N₂ V-K (Vegard-Kaplan) system in the aurora have been reinterpreted using an atmospheric model based on mass spectrometer measurements in an aurora of similar intensity at the same time of year. In contrast to the original interpretation, we find that population by cascade from the C³Π_u and B³Π_g states in the A³Σ_u⁺v=0,1 levels, as calculated using recently measured electron excitation cross sections, accurately accounts for the observed relative emission rates (IV-K/12PG_0.0). In addition there is no need to change the production rate of A ³ Σ _u⁺ molecules relative to that of C³Π_uv=0 as a function of altitude in order to fit the profile of the deactivation probability to the atmospheric model. Quenching of A ³ Σ _u⁺ molecules at high altitudes is dominated by atomic oxygen. The rate constants for the v=0 and v=1 levels are 8 × 10⁻¹¹ cm³ sec⁻¹ and 1.7 × 10⁻¹⁰ cm³ sec⁻¹ respectively, as determined using the model atmosphere mentioned above. Recent observations with a helium cooled mass spectrometer suggest that conventional mass spectrometer measurements tend to underestimate the atomic oxygen relative concentration. The rate coefficients may therefore be too large by as much as a factor of 3. Below 130 Km we find that it is possible to account for the deactivation in bright auroras by invoking large nitric oxide concentrations, similar to those recently observed mass spectrometrically and using a rate constant of 8 × 10⁻¹¹ cm³ sec⁻¹ for both the v=1 levels. This rate constant is very nearly the same as that measured in the laboratory (7 × 10⁻¹¹ cm³ sec⁻¹). Molecular oxygen appears not to play a significant role in deactivating the lower A ³ Σ _u⁺ levels.     

On the role of the motion of the photospheric footpoints of the magnetic field lines

By means of a simple relation between the velocity v of the fluid particle and the velocity v_f of the photospheric footpoint of the magnetic field line v_z and B_z being respectively the components of v and the magnetic field B normal to the photospheric surface, it is shown formally that through the phtospheric surface the transport of all the quantities attributed to the magnetic field, such as the magnetic flux, the magnetic energy and the helicity, is independent of v_z, and v_f is the only kinematical quantity on which the transport depends. In addition, in the neighborhood of the neutral line the velocity v_l of the moving curve of constant B_z is found to be equal approximately to the component of v or v_f in the direction of v_l. Since v_l can be measured or extimated, so can the components of v and v_f near the neutral line.     

Low amplitude anisotropic wave trains in cosmic rays

In the present study the occurrence of an unusual class of low amplitude anisotropic wave trains in the cosmic ray neutron intensity, which is distinctly different from the average diurnal variation as well as from other recognized types of low amplitude anisotropic wave trains are noted and the directional distribution in the interplanetary space determined. The major objective of this paper is to study the first three harmonics of low amplitude anisotropic wave trains of cosmic ray intensity over the period 1981–1994 for Deep River neutron monitoring station. The significant characteristic of these events is that the low amplitude wave trains shows a maximum intensity of diurnal component in a direction earlier than 18:00 h/co-rotational direction. It is noticed that these events are not caused either by the high-speed solar wind streams or by the sources on the Sun responsible for producing these streams such as polar coronal holes. However the possibility of occurrence of these events during high-speed solar wind streams cannot be denied. The occurrence of low amplitude events is dominant for positive polarity of B_z. The disturbance storm time index i.e. D_st, remains consistently negative only for majority of the low amplitude wave train events, which is never been reported earlier. The amplitude as well as direction of first two harmonics seems to remain unaffected with the variation in the D_st and A_p-index. However, the amplitude as well as direction of third harmonic found to deviates with the increase of D_st and A_p-index. The corotating streams produce significant deviations in cosmic ray intensity as well as in solar wind speed during low amplitude anisotropic wave train events.     

Mangeto-optical effects in sunspots and their vector magnetic field

The system of transfer equations of the four Stokes parameters I, Q, U, V under the action of the magneto-optical effect (i.e. the Unno-Beckers equations) are numerically solved in this paper for the magneto-sensitive lines FeI λλ 6302.499 and 5324.191 using an appropriate sunspot model. The errors in the expressions for the coefficients _r and _W in Beckers' paper [2] have been corrected for. From the results of calculations, features of the profiles of the Stokes parameters dependent on the magnetic vector have been isolated. Our computations also show that the magneto-optical effect should be taken into consideration in the measurement of the vector magnetic fields.

In the fourth section of this paper we have established a simple and convenient method for obtaining-information on the magnetic vector (including the field strength B, its inclination to the line of sight γ and its azimuth χ) from the profiles of the Stokes parameters. It consists of three steps: (1) The value of B is determined from the distance of the highest point in the V-profile from the central line. (2) γ is then found from V_max, i.e maximum value of V. (3) Lastly, the angle χ is found from Q₀, i.e. the value of Q at line centre.     

A kinematic study of the Galaxy

Using the proper motion and parallax data for 1011 O-B stars in the Hipparcos Catalogue we have derived the Oort constants, A = 17.60 ± 0.21 (km/s)/kpc, B = −14.62 ± 0.20 (km/s)/kpc, and a solar velocity V = 16.7 ± 0.10 km/s in the direction l = 45.3° ± 2.8°, b = 21.0° ± 2.3°. For a galactocentric distance of the sun of R₀ = 8.5 kpc, we then get a galactic rotational velocity of the solar neighbourhood of V_lsr = 273.9 km/s, obviously much higher than the IAU published value of 220 km/s. We have investigated the cause for this difference.     

Rotational discontinuities in an anisotropic plasma

A general theory of rotational discontinuities is developed and the changes in the components of the plasma pressure, p_| and p, and in the magnetic induction, B, are found. For a given value of λ=(p_| − p) 4πμ/B² upstream only a limited range of downstream anisotropies are possible. If λ>0.6 upstream then isotropy is not possible downstream. Some special solutions are analysed and the identification of rotational discontinuities is the solar wind is discussed.     

Changes in thermospheric molecular oxygen abundance inferred from twilight 6300 A airglow

When the local solar zenith angle, χ_L, is < 105° the 6300 A line is much stronger than expected on the basis of F region ionic recombination alone. Between 95 and 105° the additional intensity is quantitatively explained by production of O(¹D) from photolysis of O₂ in the Schumann-Runge continuum, (λλ 1300–1750 A) using current values for solar flux, atmospheric composition and quenching of O(¹D) by N₂. The Schumann-Runge (SR) component exhibits a large seasonal variation with a maximum in summer. We interpret this variation as implying a seasonal change in thermospheric O₂ abundance; the change seems largely to reflect a variation in O₂ density at the base of the diffusive regime although some contribution may come from changes in thermospheric temperature structure. Large changes in the SR component exist from day to day and with a 27 day period following a major magnetic storm. The photodissociation source becomes inadequate when x_l < 95°; at 90° more than half of the intensity comes from still another source which we identify as local photoelectron excitation of O atoms.     

A further study of the method of identifying absorption line systems in QSO spectra

The method of identifying absorption line systems in QSO spectra (Cui et al. 1983; Chen et al. 1983) is further developed here. Certain limitations of the method and their improvements are discussed. Certain other problems requiring further study are pointed out. The improved method is applied to PKS 0528-250, and gives two new absorption line systems Z_a = 0.065 and 0.0345 in addition to the four systems Z_a = 2.8110, 2.8130, 2.5275, 2.1410, consistent with the systems A₁A₂, B, C of Norton et al. (1980). However, the systems D₁, D₂, E, F and G of Chen and Norton (1984) are not recovered. The reason for this discrepancy is discussed.     

Cosmic-Ray Decreases and Interplanetary Disturbances During the Solar Magnetic Polarity Reversal in Solar Cycle 20

Using the cosmic-ray intensity data recorded with ground-based monitors at Mt. Washington and Deep River, and with cosmic-ray telescopes on Pioneer 8 and 9 spacecraft as well as the 2-hour averages of the IMF (magnitude and direction) and the solar wind bulk speed and density at 1 AU, the cosmic-ray decreases and interplanetary disturbances, that occurred during the period of solar magnetic polarity reversal in solar cycle 20, were investigated.We observed a two-step Forbush decrease on 22–23 November 1969, and a Forbush decrease on 26 November 1969, which are respectively consistent with the model of Barnden (1973), and of Parker (1963) and Barnden (1973). Only one Forbush decrease event was observed in December 1969, a period during which there was a solar magnetic polarity reversal; the Forbush decrease was attributed to a long-lived corotating high-speed solar wind stream. This is indicative that at heliolongitudes from 43° E to 70° W of S–E radial, covered by the observations, the solar magnetic polarity reversal in solar cycle 20 was not carried by, nor related to, individual transient structures, and that the reversal most probably evolved gradually.     

Maximum likelihood estimation of the mean parallax and kinematic parameters of the Pleiades

A simultaneous, maximum-likelihood determination of the distance and kinematic parameters of the Pleiades is made. The results are: distance of the cluster d = 135.56 ± 0.72 pc, dispersion σ_d = 7.66 ± 0.80 pc; space velocity V = 25.94 ± 0.13 km/s, dispersion σ_v = 0.58 ± 0.09 km/s coordinates of the convergent point A = 101.95° ± 0.47°, D = −41.36° ± 0.29°.     

Catastrophic fragmentation as a stochastic process: sizes and shapes of fragments

It is rather difficult to understand theoretically and to analyse the experimental data concerning the mass and shape distributions of fragments created by catastrophic collisions. The fragmentation process is discussed as being a purely stochastical phenomenon; the size and shape distributions obtained in this way are compared with the results of laboratory experiments. The results are presented of some computer simulations of random volume fragmentation processes; they are a 3-D generalization of the numerical experiments described in Grady and Kipp (J. Appl. Phys. 58(3), 1210–1222, 1985). The features of the size distribution are discussed, comparing it with the expectations of the Mott-Linfoot and Grady-Kipp theories. In the literature the shape of fragments is defined in terms of the ratios B/A and C/A, where A, B, C are defined as the sizes of a fragment along three orthogonal axes. The definition of the shape of a fragment cannot be considered unique, since it is not obvious in which order to define the three axes when the fragments are not ellipsoidal. A few possible methods are introduced explicity, and the resulting differences are discussed. In this light, the shape results (the mean values and the distribution of the axial ratios) obtained in recent laboratory experiments are rediscussed and critically reviewed. For what concerns the stochastical modelling, the results of various simulations, corresponding to different assumptions regarding fragmentation properties are presented. It is shown that the main features of the shape distributions from laboratory experiments cannot be satisfactorily reproduced. Comparison of the results with the outcomes of the semiempirical fragmentation model by Paolicchi et al. (Icarus 121, 126–157, 1996), as well as with some results coming out from hydrodynamical simulations, shows how only a “global” and physical model, not a purely statistical one (neither global nor “local”), can afford to reproduce the observed data.     

Space–time foam and cosmic-ray interactions

It has been proposed that propagation of cosmic-rays at extreme-energy may be sensitive to Lorentz-violating metric fluctuations (“foam”). We investigate the changes in interaction thresholds for cosmic-rays and gamma-rays interacting on the CMB and IR backgrounds, for a class of stochastic models of space–time foam. The strength of the foam is characterized by the factor (E/M_P)^a, where a is a phenomenological suppression parameter. We find that there exists a critical value of a (dependent on the particular reaction: a_crit3 for cosmic-rays, 1 for gamma-rays), below which the threshold energy can only be lowered, and above which the threshold energy may be raised, but at most by a factor of two. Thus, it does not appear possible in this class of models to extend cosmic-ray spectra significantly beyond their classical absorption energies. However, the lower thresholds resulting from foam may have signatures in the cosmic-ray spectrum. In the context of this foam model, we find that cosmic-ray energies cannot exceed the fundamental Planck scale, and so set a lower bound of 10⁸ TeV for the scale of gravity. We also find that suppression of p→pπ⁰ and γ→e⁻e⁺ “decays” favors values aa_crit. Finally, we comment on the apparent non-conservation of particle energy–momentum, and speculate on its re-emergence as dark energy in the foamy vacuum.     

Modification of the mean field theory of super-dense matter

We further discussed the Modified Mean Field Theory proposed in [1,2]. We found that this model has a new solution ml = 2; adding -meson, we calculated the symmetric energy of nuclear matter a₄; with the new parameters we calculated the equations of state with and without the addition of -meson; we compared the two; we also compared with the results for ml = 4.3; we made a preliminary theoretical investigation of the form of the vector meson mean field using the two-body correlation function; lastly, using the structural equation of neutron star we calculated the maximum mass for the various equations of state. The results are:

MMFT-2 (ml = 2): M_max = 2.22M (without -meson) M_max = 2.27M (with -meson)

MMFT-1 (ml = 4.3): M_max = 1.69M (without -meson) M_max = 1.89M (with -meson)