The cosmic ray energy spectrum between 10 and 10 eV

The energy spectrum of cosmic rays with primary energies between 10¹⁴ eV and 10¹⁶ eV has been studied with the CASA-MIA air shower array. The measured differential energy spectrum is a power law (dj/dE ∝ E^−y) with spectral indices γ of 2.66±0.02 below approximately 10¹⁵ eV and 3.00±0.05 above. A new method is used for measuring primary energy derived from ground-based data in a compositionally insensitive way. In contrast with some previous reports, the “knee” of the energy spectrum does not appear sharp, but rather a smooth transition over energies from 10¹⁵ eV to 3.0 × 10¹⁵ eV.     

Constraining cosmological models with cluster power spectra

Using extensive N-body simulations we estimate redshift space power spectra of clusters of galaxies for different cosmological models (SCDM, TCDM, CHDM, ΛCDM, OCDM, BSI, τCDM) and compare the results with observational data for Abell–ACO clusters. Our mock samples of galaxy clusters have the same geometry and selection functions as the observational sample which contains 417 clusters of galaxies in a double cone of galactic latitude |b|>30° up to a depth of 240 h⁻¹ Mpc. The power spectrum has been estimated for wave numbers k in the range 0.03k0.2 h Mpc⁻¹. For k>k_max0.05 h Mpc⁻¹ the power spectrum of the Abell–ACO clusters has a power-law shape, P(k)∝kⁿ, with n≈−1.9, while it changes sharply to a positive slope at k<k_max. By comparison with the mock catalogues SCDM, TCDM (n=0.9), and also OCDM with Ω₀=0.35 are rejected. Better agreement with observation can be found for the ΛCDM model with Ω₀=0.35 and h=0.7 and the CHDM model with two degenerate neutrinos and Ω_HDM=0.2 as well as for a CDM model with broken scale invariance (BSI) and the τCDM model. As for the peak in the Abell–ACO cluster power spectrum, we find that it does not represent a very unusual finding within the set of mock samples extracted from our simulations.     

The Cr and Sr sources in BOREXINO as tools for neutrino magnetic moment searches

We calculate the event rates induced by a ⁵¹Cr ν_e source and by a ⁹⁰Sr---⁹⁰Y source in BOREXINO through elastic scattering on electrons, assuming a nonzero neutrino magnetic moment μ_ν. We consider a source activity of about 2 MCi and estimate the solar ν (“source-off”) background for various oscillation scenarios. It is shown that values of μ_ν as low as 0.5 × 10⁻¹⁰μ_B ( 0.2 × 10⁻¹⁰μ_B) can be proved with the ⁵¹Cr source (⁹⁰Sr source) in about 100 days of data taking.     

Small-scale auroral arc distortions

Small-scale spatially periodic distortions of auroral forms have been studied utilizing low-light level television observations made at various locations in the Northern and Southern Hemispheres. The most commonly observed features were folds and vortex-like curl formations. The curls, identified here with the Kelvin-Helmholtz instability due to fluid shear, invariably had a counterclockwise rotational shape and motion when viewed in a direction anti-parallel to the Earth's magnetic field. The typical measured wavelength (5 km) and measured growth rate (4.2 sec⁻¹) were used to evaluate the Kelvin-Helmholtz dispersion relation for the apparent shear ω_s = ∂ ν_x/ ∂y (28 sec⁻¹). The apparent horizontal velocities of both folds (0–5 km/sec) and curls (0–22 km/sec) were invariably observed to be counterclockwise with respect to the multiple arc centre when viewed antiparallel to B. Consistent agreement between rotational shape and rotational motion suggests that the apparent growth rate and the apparent horizontal velocities closely approximate the actual values. If the shear results from E×B drifts in a space charge field, the calculated value for ω_s, implies an unneutralized electron density 0–1 cm⁻³ and a ΔE across the arc element 500mV/m. The velocity measurements indicate that the ΔE values for individual elements can combine to produce transient electric fields at the edges of multiple arcs as high as 1000 mV/m.     

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.     

Nature of solar-cycle and heliomagnetic-polarity dependence of cosmic rays, inferred from their correlation with heliomagnetic spherical surface harmonics in the period 1976–1985

Correlation of cosmic-ray intensity (I) with the solar magnetic field expanded into the spherical surface harmonics, B_ns(n 9), by Hoeksema and Scherrer has been studied using the following regression equation:

, where are subgroups of B_ns classified in ascending order of n, and τ_i is the time lag of I behind correlation coefficient between the observed and simulated intensities (I_obs, I_sml) in the period 1976–1985 is 0.87 and considerably better than that derived from any single index of solar activity. The lag time τ₃ is greater than others, indicating that the higher order magnetic disturbances effective to the cosmic-ray modulation have a longer lifetime in space than the lower order disturbances. The rigidity spectrum of the cosmic-ray intensity variation responsible for A_I due to the dipole moment is harder than those for others (A₂,A₃), indicating that the lowest order (i.e. largest scale) magnetic disturbances can modulate cosmic rays more effectively than the higher order disturbances. As another result of the present analysis, it has been found that the intensity depends also on the polarity of the polar magnetic field of the Sun; the residual (I_obs−I_sml) of the simulation changes its sign from positive to negative with a time lag (0–5 Carrington rotation periods) behind the directional change of the solar magnetic dipole moment from northward to southward, and has a softer rigidity spectrum than A_iS. The dependence is consistent with the result having been obtained in the previous period, 1936–1976, by one (K.N.) of the present authors. The polarity dependence can be found also in the 22-year variation of the time lags obtained every solar cycle in the period 1936–1985. The theoretical interpretation of these polarity dependences is discussed on the basis of the diffusion-convection-drift model.     

Density profiles of dark matter halos with anisotropic velocity tensors

We present density profiles, that are solutions of the spherical Jeans equation, derived under the following two assumptions: (i) the coarse grained phase-density follows a power-law of radius, ρ/σ³ ∝ r⁻, and (ii) the velocity anisotropy parameter is given by the relation β_a(r)=β₁+2β₂ (r/r_*)/[1+(r/r_*)²] where β₁, β₂ are parameters and r_* equals twice the virial radius, r_vir, of the system. These assumptions are well motivated by the results of N-body simulations. Density profiles have increasing logarithmic slopes γ, defined by γ=−d ln ρ/d ln r. The values of γ at r=10^−2.5r_vir, a distance where the systems could be resolved by large N-body simulations, lie in the range 1.0–1.6. These inner values of γ increase for increasing β₁ and for increasing concentration of the system. On the other hand, slopes at r=r_vir lie in the range 2.42–3.82. A model density profile that fits well the results at radial distances between 10⁻³r_vir and r_vir and connects kinematic and structural characteristics of spherical systems is described.     

Average high latitude magnetic field: Variation with interplanetary sector and with season—II : Comparison of disturbance levels and discussion of ionospheric currents

Average high latitude magnetic field data from northern observatories are examined for three ranges of magnetic disturbance level, Kp = 1− to 1+,2− to 3+ and ≥ 4−. Except for 0–8^h MLT, 55–78° invariant latitude, during away interplanetary magnetic field sectors, the variations between season and sector have the same characteristics at all Kp ranges. Because the amplitude of sector differences is much larger at sunlit local times than in the midnight sector, it is concluded that the current system of Svalgaard (1973) is not adequate to describe the sector variations in magnetic disturbance. Other current systems are discussed briefly. The disturbance morphology and seasonal variation at all Kp levels confirms the results of previous studies which indicate that latitudinally broad current systems, like S_q^p and nonionospheric sources are present in addition to latitudinally narrow electrojet currents. Comparison of data between Kp levels indicates that the Harang discontinuity shifts toward earlier MLT with increasing Kp level.     

A model independent approach to future solar neutrino experiments

We discuss here what model independent information about properties of neutrinos and of the sun can be obtained from future solar neutrino experiments (SNO, Super-Kamiokande). It is shown that in the general case of transitions of solar ν_e's into νμ and/or ντ the initial ⁸B neutrino flux can be measured by the observation of NC events. From the CC measurements the ν_e survival probability can be determined as a function of neutrino energy. The general case of transitions of solar ν_e's into active as well as sterile neutrinos is considered. A number of relations between measurable quantities the test of which will allow to answer the question whether there are sterile neutrinos in the solar neutrino flux on the earth are derived. Transitions of solar ν_e's into active and sterile states due to neutrino mixing and Dirac magnetic moments or into active left-handed neutrinos and active right-handed antineutrinos due to neutrino mixing and Majorana transition magnetic moments are also considered. It is shown that future solar neutrino experiments will allow to distinguish between the cases of Dirac and Majorana magnetic moments.     

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°.     

Structure in the two-dimensional distribution of Abell-ACO clusters

We study the structure in the two-dimensional distribution of Abell-ACO clusters of distance classes D = 5 and D = 6, using cluster analysis. The results show that there are real structures on scales ≥ 100h⁻¹ Mpc, thus providing additional evidence for the existence of structure on the super-large scale.     

The cosmic ray primary composition in the “knee” region through the EAS electromagnetic and muon measurements at EAS-TOP

The evolution of the cosmic ray primary composition in the energy range 10⁶–10⁷ GeV (i.e. the “knee” region) is studied by means of the e.m. and muon data of the Extensive Air Shower EAS-TOP array (Campo Imperatore, National Gran Sasso Laboratories). The measurement is performed through: (a) the correlated muon number (N_μ) and shower size (N_e) spectra, and (b) the evolution of the average muon numbers and their distributions as a function of the shower size. From analysis (a) the dominance of helium primaries at the knee, and therefore the possibility that the knee itself is due to a break in their energy spectrum (at E_k^He=(3.5±0.3)×10⁶ GeV) are deduced. Concerning analysis (b), the measurement accuracies allow the classification in terms of three mass groups: light (p,He), intermediate (CNO), and heavy (Fe). At primary energies E₀≈10⁶ GeV the results are consistent with the extrapolations of the data from direct experiments. In the knee region the obtained evolution of the energy spectra leads to: (i) an average steep spectrum of the light mass group (γ_p,He>3.1), (ii) a spectrum of the intermediate mass group harder than the one of the light component (γ_CNO2.75, possibly bending at E_k^CNO≈(6–7)×10⁶ GeV), (iii) a constant slope for the spectrum of the heavy primaries (γ_Fe2.3–2.7) consistent with the direct measurements. In the investigated energy range, the average primary mass increases from lnA=1.6–1.9 at E₀1.5×10⁶ GeV to lnA=2.8–3.1 at E₀1.5×10⁷ GeV. The result supports the standard acceleration and propagation models of galactic cosmic rays that predict rigidity dependent cut-offs for the primary spectra of the different nuclei. The uncertainties connected to the hadronic interaction model (QGSJET in CORSIKA) used for the interpretation are discussed.     

Helioseismological constraint on solar axion emission

Helioseismological sound-speed profiles severely constrain possible deviations from standard solar models, allowing us to derive new limits on anomalous solar energy losses by the Primakoff emission of axions. For an axion-photon coupling g_ay 5 × 10⁻¹⁰ GeV⁻¹, the solar model is almost indistinguishable from the standard case, while g_ay 10 × 10⁻¹⁰ GeV⁻¹ is probably excluded, corresponding to an axion luminosity of about 0.20 L. This constraint on g_ay is much weaker than the well-known globular-cluster limit, but about a factor of 3 more restrictive than previous solar limits. Our result is primarily of interest to the large number of current or proposed search experiments for solar axions because our limit defines the maximum g_ay for which it is self-consistent to use a standard solar model to calculate the axion luminosity.     

A possible gas for solar neutrino spectroscopy

We discuss the possibility of using pure CF₄ to fill a 2000 m³ Time Projection Chamber in order to detect the solar neutrinos through the elastic scattering v_ee⁻ → v_ee⁻, with the threshold of 100 keV on the kinetic energy of the recoiling electron. In a volume of 2000 m³ of CF₄ at normal pressure and room temperature, which corresponds to a mass of 7.4 ton, we expect ~ 3300 of such events per year. The detector can give the spectrum of the low energy neutrinos from the Sun and it can identify solar neutrinos of different origin: pp, ⁷Be, and, eventually, ⁸B. We find that ¹⁴C is a possible severe source of background: it is necessary to have a ratio ¹⁴C/¹²C lower than 10⁻¹⁹ in order to be able to identify the pp neutrinos.     

A derivation of the rate of spin-down of the Earth using ancient Chinese records of central eclipses

We describe a method of studying the variation in the rotation of the Earth in historical times using ancient records of central solar eclipses. From the records of central eclipses covering a period of some one thousand four hundred years (7th century B.C. to 7th century A.D.), given in the Chinese annals, we selected 88 entries and calculated the rotation variation parameter ΔT. Our result shows that the Earth is spinning down at a relative rate of − (1.88 ± 0.10) × 10^−10/yr.     

The configuration of the loop-like structure of the solar atmosphere

An integral, governing steady flows in an isolated thin magnetic flux tube in the hydrostatic plane-stratified atmosphere, has been obtained. The integral, that we named as the shape integral, is expressed as (1 − M_A²)B cos θ = const. Here M_A² is the Alfven Mach number, B is the magnetic field strength and θ is the flux tube inclination to the horizontal. The shape integral should hold for most loop models because it represents just the momentum balance laws and has no relation to any energy balance mode. Its application to the isothermal and static cases is discussed and illustrated.     

Galactic gamma halo by heavy neutrino annihilations?

The diffused gamma halo around our Galaxy recently discovered by EGRET could be produced by annihilations of heavy relic neutrinos N (of fourth generation), whose mass is within a narrow range (M_Z/2<m_N<M_Z). Neutrino annihilation in the halo may lead to either ultrarelativistic electron pairs whose Inverse Compton Scattering on infrared and optical galactic photons could be the source of observed GeV gamma rays, or prompt 100 MeV–1 GeV photons (due to neutral pion secondaries) born by reactions. The consequent gamma flux (10⁻⁷–10⁻⁶ cm⁻² s⁻¹ sr⁻¹) is well comparable to the EGRET observed one, and it is also compatible with the narrow window of neutrino mass 45 GeV <m_N<50 GeV, recently required to explain the underground DAMA signals.The presence of heavy neutrinos of fourth generation do not contribute much to solve the dark matter problem of the Universe, but may be easily detectable by outcoming LEP II data.     

The mass ratio of Nova Muscae 1991

We report intermediate resolution H spectroscopy of the black hole candidate Nova Muscae 1991 during quiescence. We classify the companion star as a K3-K4V which contributes 85±6 percent to the total flux from the binary. The photospheric absorption lines are broadened by 106±13 kms⁻¹ with respect to template field stars, leading to a system mass ratio of q =M₁/M₂ = 7.8_−2.0^+3.4. Doppler imaging of the H line shows strong emission coming from the secondary star (EW=3.1±0.6Å) which we associate with chromospheric activity. However, the hot-spot is not detected and this may indicate a lower mass transfer rate than in other X-ray transients of comparable orbital periods. The surface brightness distribution of the accretion disk in H follows a relation I∝R^−1.1, less steep than typically observed in cataclysmic variables. Updated system parameters are also presented.     

Properties of the infrared dark cloud G79.2+0.38

The MSX infrared dark cloud G79.2+0.38 has been observed over a 11′×′ region simultaneously in the J=1-0 rotational transition lines of the ¹²CO and its isotopic molecules ¹³CO and ¹⁸CO. The dense molecular cores defined by the C¹⁸O line are found to be associated with the two high-extinction patches shown in the MSX A-band image. The two dense cores have the column density N (H₂) (5 – 12) × 10²² cm⁻² and the mean number density n (3 ± 1) × 10⁴ cm⁻³. Their sizes are 1.7 and 1.2 pc in ¹³CO(1-0) line, 1.2 and 0.6 pc in C¹⁸O(1-0) line, respectively. The masses of these cloud cores are estimated to be in the range from 2 × 10² to 2 × 10³ M. The profile of radial mean density of the cloud core can be described by the exponential function ¯n(p) p^{−0.34±0.02}. Compared with the cases of typical optical dark clouds, the abundances of the CO isotopic molecules ¹³CO and C¹⁸O in this MSX infrared dark cloud appear to be depleted by a factor of 4–11, but at present there is no evidence for any obvious variation of the relative abundance ratio X_13/18 between ¹³CO and C¹⁸O with the column density.     

Deuterium enrichment in giant planets

Recently published laboratory measurements of the isotopic exchange rate constant k⁻(T) between CD₄ and H₂ are used to calculate f(z)—the isotopic enrichment factor between CH₄ and H₂—at every level in the outer atmosphere of the giant planets. The variation of f(z) with local vertical velocity, temperature and pressure has been calculated under the assumption that atmospheres are convective and uncertainties have been calculated by error propagation. Considering only the random errors—mainly the uncertainty on k⁻(T)—the f values in the observable upper atmospheres of giant planets (i.e. at z = 0, P = 1 bar) are: f(0) = 1.25 ± 0.05, 1.38 ± 0.06, 1.68 ± 0.09, and 1.61 ± 0.08 for Jupiter, Saturn, Uranus, and Neptune, respectively. Additional systematic errors due to the uncertainty in calculating the vertical velocity in the framework of the mixing length Prandtl theory lead to an overall uncertainty on f(0) of ±0.12, ±0.15, ±0.23, and ±0.21 for each planet, respectively. The D/H ratios in H₂ derived from the measured CH₃D/CH₄ ratios in the upper atmosphere of the four giant planets are then recalculated. Uranus and Neptune seem to be enriched in deuterium with respect to the protosolar nebula but depleted relative to the Standard Mean Oceanic Water on the Earth (SMOW). However calculations based on current interior models of Neptune suggest that ices which formed the core of the planet had a D/H ratio of the order of the SMOW. The deuterium abundance in proto-Uranian ices remains uncertain. The case where water is a major constituent of the fluid envelope of Neptune is discussed. It is shown that the D/H ratio of the planet would then be higher than the value measured in hydrogen. Even in this case, the D/H ratio in proto-Neptunian ices is less than the recently revised value in P/Halley and less than the value measured in water of the Semarkona meteorite. These results suggest that the ices which formed the core of Neptune did not have an interstellar origin. Similarly, the comparison of the most recent determination of the D/H ratio in the atmosphere of Titan with the value of D/H in P/Halley suggests that this atmosphere was not formed by infalling comets but more likely from grains embedded in the sub-nebula of Saturn.