Fragmentation cross sections and search for nuclear fragments with fractional charge in relativistic heavy ion collisions

We present measurements of fragmentation cross sections of relativistic nuclei and upper limits for the production probability of nuclear fragments with fractional charge using CR39 nuclear track detectors and an automated scanning system. The measurements of the total and partial charge changing fragmentation cross sections concern 16 GeV/nucleon oxygen ions, 14.5 GeV/nucleon silicon ions and 200 GeV/nucleon sulphur ions interacting in copper and CR39 targets. No evidence for fractionally charged fragments was found requiring a minimum track length of 7 mm in CR39 detectors placed after a 14 mm copper target. The combined upper limit for the production probability of fractionally charged fragments relative to ordinary ones is at the level of 1.2–2.3 × 10^–4 (90% C.L.). The charge resolution of the CR39 detectors for an average of 10 measurements of the same track is σ = 0.05e at Z = 6.     

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.     

Charged boson stars and vacuum instabilities

We consider charged boson stars and study their effect on the structure of the vacuum. For very compact particle like “stars”, with constituent mass m_* close to the Planck mass m_Pl, i.e. m²_* = O(m²_Pl), we argue that there is electric charge Z_c, which, primarily, is due to the formation of a pion condensate (Z_c 0.5⁻¹e, where is the fine structure constant and e is the electric charge of the positron). If the charge of the “star” is larger than Z_c we find numerical evidence for a complete screening indicating a limiting charge for a very compact object. There is also a less efficient competing charge screening mechanism due to spontaneous electron-positron pair creation in which case Z_c ⁻¹e. Astrophysical and cosmological abundances of charged compact boson stars are briefly discussed in terms of dark matter.     

Constraining the power spectrum using clusters

We analyze an extended redshift sample of Abell/ACO clusters and compare the results with those coming from numerical simulations of the cluster distribution, based on the truncated Zel'dovich approximation (TZA), for a list of eleven dark matter (DM) models. For each model we run several realizations, so that we generate a set of 48 independent mock Abell/ACO cluster samples per model, on which we estimate cosmic variance effects. Other than the standard CDM model, we consider (a) Ω₀ = 1 CDM models based on lowering the Hubble parameter and/or on tilting the primordial spectrum; (b) Ω₀ = 1 Cold + Hot DM models with 0.1 ≤Ω_ν ≤0.5; (c) low-density flat ΛCDM models with 0.3 ≤Ω₀ ≤0.5. We compare real and simulated cluster distributions by analysing correlation statistics, the probability density function, and supercluster properties from percolation analysis. We introduce a generalized definition of the spectrum shape parameter Γ in terms of σ₂₅/σ₈, where σ_ris the rms fluctuation amplitude within a sphere of radius r. As a general result, we find that the distribution of galaxy clusters provides a constraint only on the shape of the power spectrum, but not on its amplitude: a shape parameter 0.18 Γ 0.25 and an effective spectral index at the 20 h⁻¹ Mpc scale −1.1 n_eff −0.9 are required by the Abell/ACO data. In order to obtain complementary constraints on the spectrum amplitude, we consider the cluster abundance as estimated using the Press-Schechter approach, whose reliability is explicitly tested against N-body simulations. By combining results from the analysis of the distribution and the abundance of clusters we conclude that, of the cosmological models considered here, the only viable models are either Cold + Hot DM ones with 0.2 Ω_ν 0.3, better if shared between two massive ν species, and ΛCDM ones with 0.3 Ω₀0.5.     

Whistler wave propagation in magnetospheric ducts (in the equatorial region)

An analytical theory of whistler wave propagation in axially symmetric field-aligned density ducts is developed. Both enhancements and rarefactions of the density (crests and troughs) are considered. Simple equations giving the dependence of the number of modes n upon the angular frequency ω are derived. From these results it follows that in density crests n decreases when ω approaches ω_c/2 for ω < ω_c/2. The limiting frequency of the wave trapping is calculated. An analytical investigation of wave attenuation in a density crest due to wave leakage is presented. An analysis of the whistler modes in density troughs for ω < ω_c/2, ω > ω_c/2, and ω → ω_c/2 shows that the number of modes is of the same order of magnitude in these three cases.     

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.     

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.     

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.     

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.     

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.     

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.     

Whistler wave instabilities in collisionless current sheet

Instability of whistler wave in collisionless current sheet is studied with numerical solution of the general dispersion relation obtained in Ref.[4] for the physical model A. As revealed by the results, the whistler wave can be directly absorbed by collisionless current sheets. On the neutral sheet (z/d_i = 0) oblique whistler waves over a rather wide range of wave numbers can propagate, while they are basically stable. In the ionic inertial region (z/d_i < 1), the obliquely propagating whistler wave is unstable. On the edge of the ionic inertial region (z/d_i = 1), the whistler wave is still unstable, with an increase in the growth rate, and in the frequency of the unstable wave. The growth rate is larger for the whistler wave propagating towards the neutral sheet (k_zd_i < 0) than away from the neutral sheet (k_zd_i > 0).     

Full-wave VLF modes in a cylindrically symmetric enhancement of plasma density

VLF waves in the magnetosphere may be guided by ducts, i.e. enhancements of plasma density aligned with the geomagnetic field. It is required to account for the good transmission properties of these ducts as some leakage of energy is expected, with consequent attenuation. We consider a cylindrical duct, whose axis is parallel to a uniform magnetic field. The electron density of a cold plasma is a function only of radial distance r from the axis, taking constant values in inner and outer regions r < a and r > b, and varying smoothly in the duct wall a < r < b. We compute full-wave ‘trapped’ VLF modes. Results are presented for a range of values of the parameters of the model. In general we find that attenuation of leading modes is very low, agreeing with observation. Specific features are that leakage (arising from mode conversion in the duct wall) increases as the wall is made thinner, that leakage for the leading modes is less for weaker ducts, that many modes may propagate, and that, because of interference effects, some higher order modes also suffer little attenuation.     

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.     

Adiabatic vibrational frequencies of local and nonlocal convection models of the sun

The eigen-vibrational frequencies of Xiong Da-run's nonlocal and local convection models of solar envelope are calculated and compared. The differences between the observational and theoretical vibrational frequencies are less than 1%. They can be divided into two isolated groups. For modes with l ≥ 60, all the differences between observed and theoretical eigen-vibrational frequencies are distributed in a narrow and inclined belt in the (Δv ≈ v)-diagram. This shows that the theoretical model of solar convective region can approximately reflect the intrinsic structure of the sun in the region of r = (0.70–0.95)R. The discrepancies between the theoretical and observational frequencies come from the outer layers. For modes with l < 60, the theoretical vibrational frequency is smaller than the observational one. This implies that the temperature of the upper part of the convectively unstable region is rather low. The frequency difference is more dispersed in the local convection model than in the nonlocal convection model. For the intermediate- and low-frequency ranges (v < 3000), the difference between the two models is small, while for the high-frequency range (v ≥ 3000) the frequency in the local model is higher than in the nonlocal model. This means that the temperature of the radiation region beneath the convective region is higher in the local convection model than in the nonlocal convection model. The nonlocal model is nearer to the observation than the local model.     

Separation of gamma and hadron initiated air showers with energies between 20 and 500 TeV

The discrimination between air showers initiated by γ rays and by hadrons is one of the fundamental problems in experimental cosmic-ray physics. The physics of this ‘γ/hadron separation’ is discussed in this paper. We restrict ourselves to the energy range from about 20 to 500 TeV, and take only the information contained in the lateral Čerenkov light distribution and the number of electrons at the detector level into consideration. An understanding of the differences between air showers generated by γ rays and those due to hadrons leads us to formulate suitable observables for the separation process. Angle integrating Čerenkov arrays (AICA) offer a promising new approach to ground-based γ-ray astronomy in the energy region from about 20 to 500 TeV. In order to establish this technique, an efficient suppression of the overwhelming hadronic background radiation is required. As an example for our general discussion, we present one method for γ/hadron separation in AICAs called ‘LES’. It is based on the simultaneous determination of the shower size and some characteristic parameters of the lateral distribution of the Čerenkov light. The potential inherent within this technique is demonstrated in quantitative detail for the existing ‘AIROBICC’ AICA. We also propose an objective measure of the intrinsic sensitivity of a detection scheme in ground-based γ-ray astronomy, the ‘reduced quality factor’. It is shown that AICAs may reach a sensitivity to γ-ray point sources in the high VHE range similar to that of the Čerenkov-telescope imaging technique in the low VHE region.     

XMM-Newton observations of the nearby brown dwarf LP 944-20

The nearby (d=5.0 pc) brown dwarf LP 944-20 was observed with the XMM-Newton satellite on 07 January 2001. The target was detected with the Optical Monitor (V=16.736±0.081), but it was not detected during the ≈48 ks observation with the X-ray telescopes. We determine a 3σ upper limit for the X-ray emission from this object of L_X<3.1×10²³ ergs·s⁻¹, equivalent to a luminosity ratio upper limit of log(L_X/L_bol)≤−6.28. This measurement improves by a factor of three the previous Chandra limit on the quiescent X-ray flux. This is the most sensitive limit ever obtained on the quiescent X-ray emission of a brown dwarf. Combining the XMM-Newton data with previous ROSAT and Chandra data, we derive flare duty cycles as a function of their luminosities. We find that very strong flares [Log(L_X/L_bol)>−2.5] are very rare (less than 0.7% of the time). Flares like the one detected by Chandra [Log(L_X/L_bol)=−4.1] have a duty cycle of about 6%, which is lower than the radio flare duty cycle (13%). When compared with other M dwarfs, LP 944-20 appears to be rather inactive in X-rays despite of its relative youth, fast rotation and its moderately strong activity at radio wavelengths.     

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.     

Implications of modified AGN spectra due to absorption by infrared photons

Recent observations of high energy photons from active galactic nuclei (AGN) have lead to a renewed interest in the effect of gamma-ray cascading in background radiation fields. In recent years, numerous authors have explored possible modifications to the observed high energy photon spectra from AGN due to propagation through such fields. This paper will re-examine a number of these issues. The major conclusions of this paper are: (1) The reaction γ + γ → e⁺ + e⁻ provides a highly sensitive probe of background infrared (IR) fields. However little can be inferred about the specific nature of the IR background or the effect on spectra from more distant AGN from a study of nearby sources alone. Currently, only upper limits to the background IR density can be established. (2) The contribution of secondary photons from pair-cascading off of background microwave and infrared radiation is most likely unobservable in the regime of 1 TeV unless the strength of extragalactic magnetic fields are much less than 10⁻¹³ G and the inherent source spectra continue to much higher energies. The possible contribution from this process may be conclusively ruled out through further spectral measurements of AGN in the TeV regime during high and low states. (3) There is little hope of unambiguously extracting the values of either the Hubble constant or Ω from γ-γ attenuation measurements. (4) The sensitivity of the attenuation to the density of the background IR produced prior to the epoch of the observed source suggests a future possibility to probe directly the evolution of the background IR radiation through multiple source observations.     

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