No Arabic abstract
We implemented a website to deal with main effects on Cosmic Ray access to the Earth, i.e. the Solar Modulation and the Geomagnetic Field effect. In helmod.org the end user can easily access a web interface to results catalog of the HelMod Monte Carlo Code. This Model uses a Monte Carlo Approach to solves the Parker Transport Equation, obtaining a modulated proton flux for a period (monthly average) between January 1990 and december 2007. geomagsphere.org is instead based on GeoMag Backtracing Code, that solves the Lorentz equation with a Runge-Kutta method of 6th order, and, reversing charge sign and velocity, reconstruct particle trajectories in the Earth Magnetosphere back in time. We use last models of internal (IGRF-11) and external (Tsyganenko 1996 -T96- and 2005 -T05-) field components valid up to 2015. Particles are backtraced to the outer (magnetopause) or inner boundary to separate Primary (allowed trajectory) from Secondary (forbidden) Cosmic Rays. This code has been used both for reproducing known effects as East-West effect and rigidity cutoff calculations. In geomagsphere.org the user can choose the external field model from Tsyganenko (T96 or T05) and obtain for a fixed position and date from 1st Jan. 1968 (T96) and 1st Jan. 1995 (T05) respectively till 31$^{st}$ Dec 2012, the vertical rigidity cutoff estimation obtained with the backtracing technique with a rigidity step of 0.1 GV. For a more precise calculation (0.01 GV), requiring more CPU time, results are sent to the user by email (mail model)
Star formation in galaxies appears to be self-regulated by energetic feedback processes. Among the most promising agents of feedback are cosmic rays (CRs), the relativistic ion population of interstellar and intergalactic plasmas. In these environments, energetic CRs are virtually collisionless and interact via collective phenomena mediated by kinetic-scale plasma waves and large-scale magnetic fields. The enormous separation of kinetic and global astrophysical scales requires a hydrodynamic description. Here, we develop a new macroscopic theory for CR transport in the self-confinement picture, which includes CR diffusion and streaming. The interaction between CRs and electromagnetic fields of Alfvenic turbulence provides the main source of CR scattering, and causes CRs to stream along the magnetic field with the Alfven velocity if resonant waves are sufficiently energetic. However, numerical simulations struggle to capture this effect with current transport formalisms and adopt regularization schemes to ensure numerical stability. We extent the theory by deriving an equation for the CR momentum density along the mean magnetic field and include a transport equation for the Alfven-wave energy. We account for energy exchange of CRs and Alfven waves via the gyroresonant instability and include other wave damping mechanisms. Using numerical simulations we demonstrate that our new theory enables stable, self-regulated CR transport. The theory is coupled to magneto-hydrodynamics, conserves the total energy and momentum, and correctly recovers previous macroscopic CR transport formalisms in the steady-state flux limit. Because it is free of tunable parameters, it holds the promise to provide predictable simulations of CR feedback in galaxy formation.
We develop a model anisotropy best-fitting to the two-dimensional sky-map of multi-TeV galactic cosmic ray (GCR) intensity observed with the Tibet III air shower (AS) array. By incorporating a pair of intensity excesses in the hydrogen deflection plane (HDP) suggested by Gurnett et al., together with the uni-directional and bi-directional flows for reproducing the observed global feature, this model successfully reproduces the observed sky-map including the skewed feature of the excess intensity from the heliotail direction, whose physical origin has long remained unknown. These additional excesses are modeled by a pair of the northern and southern Gaussian distributions, each placed ~50 degree away from the heliotail direction. The amplitude of the southern excess is as large as ~0.2 %, more than twice the amplitude of the northern excess. This implies that the Tibet AS experiment discovered for the first time a clear evidence of the significant modulation of GCR intensity in the heliotail and the asymmetric heliosphere.
In dense and cold regions of the interstellar medium (ISM), molecules may be adsorbed onto dust grains to form the ice mantles. Once formed, they can be processed by ionizing radiation coming from stellar or interstellar medium leading to formation of several new molecules in the ice. Among the different kind of ionizing radiation, cosmic rays play an important role in the solid-phase chemistry because of the large amount of energy deposited in the ices. The physicochemical changes induced by the energetic processing of astrophysical ices are recorded in a intrinsic parameter of the matter called complex refractive index (CRI). In this paper, we present for the first time a catalogue containing 39 complex refractive indices (n, k) in the infrared from 2.0 - 16.6 micrometer for 13 different water-containing ices processed in laboratory by cosmic ray analogs. The calculation was done by using the NKABS (acronym of determination of N and K from ABSorbance data) code, which employs the Lambert-Beer and Kramers-Kronig equations to calculate the values of n and k. The results are also available at the website: http://www1.univap.br/gaa/nkabs-database/data.htm. As test case, a H2O:NH3:CO2:CH4 ice was employed in a radiative transfer simulation of a prototoplanetary disk to show that these data are indispensable to reproduce the spectrum of YSOs containing ices.
The cosmic rays modulation inside the heliosphere is well described by a transport equation introduced by Parker in 1965. To solve this equation several approaches were followed in the past. Recently the Monte Carlo approach becomes widely used in force of his advantages with respect to other numerical methods. In the Monte Carlo approach, the transport equation is associated to a fully equivalent set of Stochastic Differential Equations. This set is used to describe the stochastic path of a quasi-particle from a source, e.g., the interstellar medium, to a specific target, e.g., a detector at Earth. In this work, we present both the Forward-in-Time and Backward-in-Time Monte Carlo solutions. We present an implementation of both algorithms in the framework of HelMod Code showing that the difference between the two approach is below 5% that can be quoted as the systematic uncertain of the Method itself.
Local interstellar spectra (LIS) of secondary cosmic ray (CR) nuclei, lithium, beryllium, boron, and partially secondary nitrogen, are derived in the rigidity range from 10 MV to ~200 TV using the most recent experimental results combined with the state-of-the-art models for CR propagation in the Galaxy and in the heliosphere. The lithium spectrum appears somewhat flatter at high energies compared to other secondary species that may imply a primary lithium component. Two propagation packages, GALPROP and HelMod, are combined to provide a single framework that is run to reproduce direct measurements of CR species at different modulation levels, and at both polarities of the solar magnetic field. An iterative maximum-likelihood method is developed that uses GALPROP-predicted LIS as input to HelMod, which provides the modulated spectra for specific time periods of the selected experiments for the model-data comparison. The proposed LIS accommodate the low-energy interstellar spectra measured by Voyager 1, HEAO-3, and ACE/CRIS as well as the high-energy observations by PAMELA, AMS-02, and earlier experiments that are made deep in the heliosphere. The interstellar and heliospheric propagation parameters derived in this study are consistent with our earlier results for propagation of CR protons, helium, carbon, oxygen, antiprotons, and electrons.