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.     

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.     

Flux limits on ultra high energy neutrinos with AMANDA-B10

Data taken during 1997 with the AMANDA-B10 detector are searched for a diffuse flux of neutrinos of all flavors with energies above 10¹⁶ eV. At these energies the Earth is opaque to neutrinos, and thus neutrino induced events are concentrated at the horizon. The background are large muon bundles from down-going atmospheric air shower events. No excess events above the background expectation are observed and a neutrino flux following E⁻², with an equal mix of all flavors, is limited to E²Φ(10¹⁵ eV < E < 3 × 10¹⁸ eV)  0.99 × 10⁻⁶ GeV cm⁻² s⁻¹ sr⁻¹ at 90% confidence level. This is the most restrictive experimental bound placed by any neutrino detector at these energies. Bounds to specific extraterrestrial neutrino flux predictions are also presented.     

Comparison of AGASA data with CORSIKA simulation

An interpretation of Akeno giant air shower array (AGASA) data by comparing the experimental results with the simulated ones by cosmic ray simulation for KASCADE (CORSIKA) has been made. General features of the electromagnetic component and low energy muons observed by AGASA can be well reproduced by CORSIKA. The form of the lateral distribution of charged particles agrees well with the experimental one between a few hundred metres and 2000 m from the core, irrespective of the hadronic interaction model studied and the primary composition (proton or iron). It does not depend on the primary energy between 10^17.5 and 10²⁰ eV as the experiment shows. If we evaluate the particle density measured by scintillators of 5 cm thickness at 600 m from the core S₀(600), suffix 0 denotes the vertically incident shower) by taking into account the similar conditions as in the experiment, the conversion relation from S₀(600) to the primary energy is expressed as E (eV)=2.15×10¹⁷S₀(600)^1.015 within 10% uncertainty among the models and composition used, which suggests the present AGASA conversion factor is the lower limit. Although the form of the muon lateral distribution fits well to the experiment within 1000 m from the core, the absolute values change with hadronic interaction model and primary composition. The slope of the ρ_μ(600) (muon density above 1 GeV at 600 m from the core) vs. S₀(600) relation in experiment is flatter than that in simulation of any hadronic model and primary composition. As the experimental slope is constant from 10¹⁵ to 10¹⁹ eV, we need to study this relation in a wide primary energy range to infer the rate of change of chemical composition with energy.     

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.     

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.     

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.     

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.     

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.     

A prediction for three neutrino masses and mixings and similarity between quark and lepton mixings

If neutrinos have mass, we give reasons for a possible pattern of three (squaed) mass eigenvalues: m₁² (2.8−5.8) (eV)², m₂² 0.01 (eV)², m₃² (1.5−1) × 10⁻⁴ (eV)². The flavor states ν_μ and ν_e are mixtures of the eigenstates with m₂ and m₃ with a significant mixing, corresponding to an effective mixing angle of about 0.45. The ν_τ is nearly the state with m₁; the other two effective mixing angles are about an order of magnitude smaller than 0.45. There is a marked similarity to mixing in the quark sector.     

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.     

Synchrotron radiation from the Crab Nebula discriminates between models of space–time foam

It has been argued by Jacobson, Liberati and Mattingly that synchrotron radiation from the Crab Nebula imposes a stringent constraint on any modification of the dispersion relations of the electron that might be induced by quantum gravity. We supplement their analysis by deriving the spectrum of synchrotron radiation from the coupling of an electrically charged particle to an external magnetic fields in the presence of quantum-gravity effects of the general form (E/M_QG). We find that the synchrotron constraint from the Crab Nebula practically excludes 1.74 for M_QGm_P=1.2×10¹⁹ GeV. On the other hand, this analysis does not constrain any modification of the dispersion relation of the photon that might be induced by quantum gravity. We point out that such quantum-gravity effects need not obey the equivalence principle, a point exemplified by the Liouville-string D-particle model of space–time foam. This model suggests a linear modification of the dispersion relation for the photon, but not for the electron, and hence is compatible with known constraints from the Crab Nebula and elsewhere.     

Observation and Study of the CO, CO, HCO and N2H Molecular Lines Towards IRAS 02232+6138

Using the 13.7 m millimeter-wave telescope at the Qinghai Station of Purple Mountain Observatory, we have made observations of ¹³CO, C¹⁸O, HCO⁺ and N₂H⁺ molecular lines towards IRAS 02232+6138. As the excitation density of the probe molecule increases from ¹³CO to HCO⁺, the size of the cloud core associated with IRAS 02232+6138 decreases from 2.40 pc to 0.54 pc, and the virial mass of the cloud core decreases from 2.2 × 10³M to 5.1 × 10²M. A bipolar molecular outflow is found towards IRAS 02232+6138. Using the power function n(r) ∝ r⁻ to fit the spatial density structure of the cloud core, we obtain the power-law index  = 2.3 − 1.2; and we find that, as the probed density increases, the power function becomes more flat. The abundance ratio of ¹³CO to C¹⁸O is 12.4 ± 6.9, comparable with the values 11.8 ± 5.9 for dark clouds and the values 9.0–15.6 for massive cores. The abundance of N₂H⁺ molecules is 3.5 ± 2.5 × 10⁻¹⁰, consistent with the value 1.0 − 5.0 × 10⁻¹⁰ for dark cloud cores and the value 1.2 − 12.8 × 10⁻¹⁰ for massive cores. The abundance of HCO⁺ molecules is 0.9 ± 0.5 × 10⁻⁹, close to the value 1.6 − 2.4 × 10⁻⁹ for massive cores. An increase of HCO⁺ abundance in the outflow region was not found. Combining with the IRAS data, the luminosity-mass ratio of the cloud core is obtained in the range 37–163(L/M). Based on the IRAS luminosity, it is estimated that a main-sequence O7.5 star is probably embedded in the IRAS 02232+6138 cloud core.     

Robust limits on Lorentz violation from gamma-ray bursts

We constrain the possibility of a non-trivial refractive index in free space corresponding to an energy-dependent velocity of light: c(E)  c₀(1 − E/M), where M is a mass scale that might represent effect of quantum-gravitational space-time foam, using the arrival times of sharp features observed in the intensities of radiation with different energies from a large sample of gamma-ray bursters (GRBs) with known redshifts. We use wavelet techniques to identify genuine features, which we confirm in simulations with artificial added noise. Using the weighted averages of the time-lags calculated using correlated features in all the GRB light curves, we find a systematic tendency for more energetic photons to arrive earlier. However, there is a very strong correlation between the parameters characterizing an intrinsic time-lag at the source and a distance-dependent propagation effect. Moreover, the significance of the earlier arrival times is less evident for a subsample of more robust spectral structures. Allowing for intrinsic stochastic time-lags in these features, we establish a statistically robust lower limit: M > 0.9 × 10¹⁶ GeV on the scale of violation of Lorentz invariance.     

Generalized limits to the number of light particle degrees of freedom from big bang nucleosynthesis

We compute the big bang nucleosynthesis limit on the number of light neutrino degrees of freedom in a model-independent likelihood analysis based on the abundances of ⁴He and ⁷Li. We use the two-dimensional likelihood functions to simultaneously constrain the baryon-to-photon ratio and the number of light neutrinos for a range of ⁴He abundances Y_p = 0.225–0.250, as well as a range in primordial ⁷Li abundances from (1.6 to 4.1) ×10⁻¹⁰. For (⁷Li/H)_p = 1.6 × 10⁻¹⁰, as can be inferred from the ⁷Li data from Population II halo stars, the upper limit to N_ν based on the current best estimate of the primordial ⁴He abundance of Y_p = 0.238 is N_ν < 4.3 and varies from N_ν < 3.3 (at 95% C.L.) when Y_p = 0.225 to N_ν < 5.3 when Y_p = 0.250. If ⁷Li is depleted in these stars the upper limit to N_ν is relaxed. Taking (⁷Li/H)_p = 4.1 × 10⁻¹⁰, the limit varies from N_ν < 3.9 when Y_p = 0.225 to N_ν 6 when Y_p = 0.250. We also consider the consequences on the upper limit to N_ν if recent observations of deuterium in high-redshift quasar absorption-line systems are confirmed.     

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.     

The HNCO ring in the Sgr B2 region

The large-scale structure associated with the 2′N HNCO peak in Sgr B2 [Minh, Y.C., Haikala, L., Hjalmarson, Å., Irvine, W.M., 1998. ApJ 498, 261 (Paper I)] has been investigated. A ring-like morphology of the HNCO emission has been found; this structure may be colliding with the Principal Cloud of Sgr B2. This “HNCO Ring” appears to be centered at (l,b) = (0.7°,−0.07°), with a radius of 5 pc and a total mass of 1.0 × 10⁵ to 1.6 × 10⁶ M. The expansion velocity of the Ring is estimated to be 30–40 km s⁻¹, which gives an expansion time scale of 1.5 × 10⁵ year. The morphology suggests that collision between the Ring and the Principal Cloud may be triggering the massive star formation in the Sgr B2 cloud sequentially, with the latest star formation taking place at the 2′N position. The chemistry related to HNCO is not certain yet, but if it forms mainly via reaction with the evaporated OCN⁻ from icy grain mantles, the observed enhancement of the HNCO abundance can be understood as resulting from shocks associated with the collision between the Principal Cloud and the expanding HNCO Ring.     

Ozone, climate and biospheric environment in the ancient oxygen-poor atmosphere

To study the climatological role of ozone in the Precambrian atmosphere and the consequences of its reduction for the ultraviolet environment of the early biosphere, a coupled one-dimensional radiative-convective and photochemical model has been developed. Oxygen levels between 10⁻⁵ and 1 time the present atmospheric level (PAL) are considered. It is shown that when the ice-albedo feedback is taken into account, relatively important temperature decreases are associated with the ozone changes linked to the progressive decrease of the oxygen level from 1 PAL to smaller values.

A similar study is performed for enhanced atmospheric CO₂ pressures (P_co2). In these conditions, the ozone column is increased at low O₂ concentrations with respect to the P_co2 = 1 PAL case. Consequently, the larger CO₂ concentration in the ancient atmosphere could have contributed to strengthen the ultraviolet screening of ozone. The surface temperature response to the ozone decrease, as well as the thermal profiles are also analyzed in these CO₂-rich models. A possible evolutionary scenario of atmospheric O₂ and CO₂ is discussed.

The consequences of these calculations for the ultraviolet environment of the primitive biosphere is discussed with a quantitative model calculating bacterial surviving rates. According to this model, the minimum ozone column being tolerable by unprotected bacteria would fall between 1 × 10¹⁸ and 4 × 10¹⁸ cm⁻², depending on the bacterial species considered and corresponding to an O₂ level somewhat lower than 10⁻² PAL. For the coccoid blue-green alga Agmenellum quadruplicatum, this minimum ozone column would be of 4.5 × 10¹⁸, a value which is only slightly less than the presently observed column in the spring time ozone hole of Antarctica.     

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.     

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