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.     

The third-order theory of artificial satellites of the earth: A first exploration

We propose a mixed analytical and numerical method as a practical means of solving the perturbation equations based on Vinti's intermediate orbit. Hori's averaging method is used to define mean elements σ^* and short-period perturbations Δσ_s. Δσ_s accurate to the third-orbit, are found by Fourier analysis, and dσ*/dt, accurate to the fourth order, from numerical averaging. The present procedure can be used to reduce satellite laser observations with a high precision.     

Cyclotron amplification of geomagnetic micropulsations Pc1 in the magnetosphere

A numerical analysis of cyclotron instabilities is carried out by computing the dispersion relation for a three component cold plasma-beam system. Rates of growth and damping for various values of the stream density are calculated from the dispersion relation. The rates of growth and damping increase monotonically as the number density of the proton stream increases. It is found that the frequencies at the rates of maximum growth and the damping decrease slightly to lower frequencies and a sharp peak at these frequencies becomes blunt. The minimum e-folding times of an ion cyclotron wave for (a) σ_s = 10⁻⁴, σ_i = 10⁻² and (b) σ_s = 10⁻¹, σ_i = 10⁻² are about 3·84 and 0·16 sec respectively in the vicinity of the equatorial plane at 6 Re, where σ_s and σ_i are the ratios of the beam density N_s and the helium ion (H₆⁺) density N_i to the total positive ions in the plasma-beam system.     

Galactic kinematics derived from classical cepheids

On the basis of radial velocity and Hipparcos proper motion data, we have analyzed the galactic kinematics of classical Cepheids. Using the 3-D Ogorodnikov-Milne model we have determined the rotational velocity of the Galaxy to be V₀ = 240.5 ± 10.2 km/s, on assuming a glactocentric distance of the Sun of R₀ = 8.5 kpc. The results clearly indicate a contracting motion in the solar neighbourhood of (∂V_θ∂θ)/R = −2.60 ± 1.07 km s⁻¹ kpc⁻¹, along the direction of galactic rotation. Possible reason for this motion is discussed. The solar motion found here is S = 18.78 ± 0.86 km/s in the direction l = 54.4° ± 2.9° and b = +26.6° ± 2.6°.     

Electron excitation and auroral emission parameters

Auroral luminosities of the main emission lines in the aurora have been calculated for excitation by an isotopic primary electron flux with spectra of the form J(E) = AE exp (−E/E₁) + B(E₂)E exp (−E/E₁). The variation of emissions from O and N₂⁺ with height are shown, as are the variations of column integrated intensities and pertinent intensity ratios with the characteristic energy E₂, this leading to a method of estimating the electron spectrum from ground observation.     

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.     

Auroral N2(c′4Σu → aΠg emission

Inspection of recent spectra presented by Sivjee (1983) show evidence of the 0–4 and 0–5 bands of the N₂(c′₄¹Σ_u⁺ → a¹Π_g) Gaydon-Herman system. In conjunction with earlier spectra, it is now possible that this band system is a significant auroral component, with an intensity approx. 7% that of the N2 2P system. The absence in aurorae of the potentially far stronger N₂(c′₄¹Σ_u⁺ → X¹Π_g) system is discussed. It is that the O₂(A³Σ_u⁺ → X³Σ_g⁻) band system is indiscernible in Sivjee's auroral spectra, under conditio the foreground nightglow is expected to be clearly visible. On the other hand, at least one relatively strong O₂(A′³Δ_u → a¹Δ_g) band appears to be present in these spectra.     

Rocket measurements of the 6300 Å and 3914 Å dayglow features

Rocket results are presented on the OI 6300 Å line and on the N₂⁺ 3914 Å band in the dayglow. An altitude range of 78–335 km is covered. Theoretical interpretations are given, using results of simultaneous measurements of electron density and electron temperature. The apparent brightness of the 6300 Å line at the base of the emitting region is found to be 13 kR, of which 5.5 kR are ascribed to excitation through the Schumann-Runge dissociation of O₂ by the solar UV radiations, 0.55 kR to the dissociative recombination of O₂⁺ and NO⁺ ions, and 0.03 kR to the excitation of O by thermal electrons. An additional source of excitation above 280 km is suggested. The deactivation of O(¹D) by O₂(X³Σ_g⁻) is found to be appreciable below 200 km, and its rate coefficient is estimated to be 2 × 10⁻¹⁰ cm³/sec. The apparent brightness of the 3914 Å band at the base of the emitting region is found to be 6.5 kR, decreasing to 3.2 kR at 330 km. Assuming that fluorescent scattering of solar radiation is the mechanism involved the distribution of N₂⁺ ions is calculated. The rate coefficients for the loss of these ions are hence calculated.     

ETON 6: A rocket measurement of the O2 Infrared Atmospheric (0-0) band in the nightglow

Co-ordinated rocket measurements of the O₂(a¹Δ_g−X³Σ_g⁻) Infrared Atmospheric (0-0) band emission profile and the atomic oxygen densities in an undisturbed night-time atmosphere are used to investigate the processes responsible for the excitation of O₂(a¹Δ_g) in the terrestrial nightglow. It is shown that three-body recombination of atomic oxygen, and subsequent energy transfer processes, can explain only part of the observed emission profile and that at least two other sources of O₂(a¹Δ_g) emission must exist. One of these additional sources, responsible for most of the emission observed below 90km, is identified as arising from the night-time residual of the very large dayglow ¹Δ_g population. The other additional source is required to explain most of the emission observed above 95km. The processes responsible for this high altitude component cannot be identified but the vertical distribution of the required source function strongly resembles the profile of the atomic oxygen density squared and suggests that a two-body radiative recombination process may be involved. However, the measured zenith emission rates can also be explained without the high altitude source of O₂(a¹Δ_g) if optical emission at 1.27 μm was induced by the rocket as it penetrated the nightglow layer.     

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.     

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.     

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.     

Gamma ray spectrum of the crab nebula in the multi TeV region

The THEMISTOCLE array of 18 Cherenkov detectors which has a 3 TeV gamma energy threshold, has detected a signal from the Crab nebula at a 5.8 standard deviation level. Information on the energy spectrum is obtained in the range 3–15 TeV. The integrated flux can be fitted with the form, Φ (> E) = (3.7 ± 0.5) × 10^-12 (E/5)-^{−1.5 ± 0.20} cm⁻² s⁻¹ (E in TeV) compatible with the extrapolation of results at lower energies. The Crab signal is used to measure the angular resolution of the multi-telescope technique. The value obtained is 2.3 mr (0.15°) in agreement with the results of simulations, and confirms the interest of this new method for multi-TeV gamma-ray detection.     

Spectrum of the relic neutrino background from past supernovae and cosmological models

It is greatly expected that the relic neutrino background from past supernovae will be detected by Superkamiokande (SK) which is now under construction. We calculate the spectrum and the event rate at SK systematically by using the results of simulations of a supernova explosion and reasonable supernova rates. We also investigate the effect of a cosmological constant, Λ, on the spectrum, since some recent cosmological observations strongly suggest the existence of Λ. We find following results. (1) The spectrum has a peak at about 3 MeV, which is much lower than that of previous estimates (6–10 MeV). (2) The event rate at SK in the range from 10 MeV to 50 MeV, where the relic neutrinos from past supernovae are dominant, is about 25h₅₀²(R_SN/0.1 yr⁻¹)(n_Gh₅₀⁻³/0.02 Mpc⁻³) events per year, where R_SN is the supernova rate in a galaxy, n_G is the number density of galaxies, and h₅₀ = H₀/(50 km/s Mpc), where H₀ is the Hubble constant. (3) The event rate is almost insensitive to Λ. The flux increases in the low energy side (< 10 MeV) with increasing Λ, but decreases in the high energy side (> 10 MeV) in models in which the integrated number of supernovae in one galaxy is fixed.     

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.     

Observation of a thin Es-layer by the eiscat radar

High resolution E-region measurements carried out on 16 November 1983 using the EISCAT incoherent scatter radar are presented. The experiment was monostatic with a vertical radar beam, and it was based on a Barker-coded four-pulse code on one frequency channel and Barker-coded single pulses on three channels. The basic integration time was 15 s and the spatial resolution 450 m. The results reveal a short-lived but intense thin sporadic E-layer at 18:00–18:06 U.T. at an altitude of about 106 km. Both before and during the event, downward ion velocities of the order of 100 m s⁻¹ are observed above this height. A convergent null in the vertical ion speed is occasionally seen at the layer altitude. The layer occurrence is associated with auroral arcs drifting across the radar beam. Simultaneous observations of the STARE radar show an ionospheric electric field of 25–30 mV m⁻¹. The field always has a westward component, which is in accordance with the observed downward plasma flow. Most of the time when the layer is intense, the field points into the NW-sector. Theoretically, this field direction should create convergent vertical plasma motion. Therefore it is suggested that the observed E_s-layer is created by the action of the auroral electric field rather than by the wind shear mechanism.     

Analysis of the variation in inclination of Intercosmos 13 Rocket, 1975-22B

The orbit of Intercosmos 13 rocket (1975-22B) has been determined at 103 epochs between 30 April 1975 and 10 April 1980 from almost 7000 observations. One hundred and three values of inclination have been determined and corrections incoporated for the effects due to zonal harmonic, lunisolar and tesseral harmonic perturbations, precession, and solid Earth tides. The modified data have been analysed to yield values of the atmospheric rotation rate, Λ rev day⁻¹, viz. Λ = 0.94 ± 0.10 at an average height of 322 ± 6 km and Λ = 1.27 ± 0.02 at 288 km. Analysis of the inclination near 14th-order resonance has indicated lumped harmonic values 10^{9 1.0}_1.4 = − 76.13 ± 12.47, 10^{9 1,0}₁₄ = − 29.89 ± 32.64, 10^{9 −1.2}₁₄ = − 63.11 ± 15.44 10^{9 −1.2}₁₄ = − 32.52 ± 26.96, for inclination 82.952°.     

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.     

Ionospheric electric field and current distribution associated with high altitude electric field inhomogeneities

A quantitative theoretical analysis of electric field and current distributions in the ionosphere is given assuming certain time variable convection field profiles at an altitude of 1250 km. A number of idealized assumptions regarding the ionospheric characteristics are defined and discussed. A qualitative discussion of a quasi-stationary configuration with an approximately curl free electric field is also given. Geomagnetically field aligned current densities i of the order 10⁻⁵−10⁻⁴A/m² are consistent with quite reasonable assumptions about the convection field E. Oscillations in E with periods of the order of 10 sec should readily be generated when σ is large. In the quasi-stationary case there may be a mechanism that strengthens and concentrates i locally under certain conditions. It is found that a number of recent high altitude observations of convection field reversals may be consistent with large potential drops along the magnetic field lines. The solutions obtained as well as some of the basic assumptions are compared with observations.     

Comparison of high-energy galactic and atmospheric tau neutrino flux

We compare the tau neutrino flux arising from the galaxy and the earth atmosphere for 10³E/GeV10¹¹. The intrinsic and oscillated tau neutrino fluxes from both sources are calculated. The intrinsic galactic ν_τ flux (E10³ GeV) is calculated by considering the interactions of high-energy cosmic-rays with the matter present in our galaxy, whereas the oscillated galactic ν_τ flux is coming from the oscillation of the galactic ν_μ flux. For the intrinsic atmospheric ν_τ flux, we extend the validity of a previous calculation from E10⁶ GeV up to E10¹¹ GeV. The oscillated atmospheric ν_τ flux is, on the other hand, rather suppressed. We find that, for 10³E/GeV5×10⁷, the oscillated ν_τ flux along the galactic plane dominates over the maximal intrinsic atmospheric ν_τ flux, i.e., the flux along the horizontal direction. We also briefly mention the presently envisaged prospects for observing these high-energy tau neutrinos.