Propagation in 3D spiral-arm cosmic-ray source distribution models and secondary particle production using Picard

We study the impact of possible spiral-arm distributions of Galactic cosmic-ray sources on the flux of various cosmic-ray nuclei throughout our Galaxy. We investigate model cosmic-ray spectra at the nominal position of the sun and at different positions within the Galaxy. The modelling is performed using the recently introduced numerical cosmic ray propagation code Picard. Assuming non-axisymmetric cosmic-ray source distributions yields new insights on the behaviour of primary versus secondary nuclei.We find that primary cosmic rays are more strongly confined to the vicinity of the sources, while the distribution of secondary cosmic rays is much more homogeneous compared to the primaries. This leads to stronger spatial variation in secondary to primary ratios when compared to axisymmetric source distribution models. A good fit to the cosmic-ray data at Earth can be accomplished in different spiral-arm models, although leading to decisively different spatial distributions of the cosmic-ray flux. These lead to different cosmic ray anisotropies, where even reproducing the data becomes possible. Consequently, we advocate directions to seek best fit propagation parameters that take into account the higher complexity introduced by the spiral-arm structure on the cosmic-ray distribution. We specifically investigate whether the flux at Earth is representative for a large fraction of the Galaxy. The variance among possible spiral-arm models allows us to quantify the spatial variation of the cosmic-ray flux within the Galaxy in presence of non-axisymmetric source distributions.     

Simulations of Galactic Cosmic Ray impacts onto the Herschel/PACS photoconductor arrays with Geant4 code

We present results of simulations performed with the Geant4 software code of the effects of Galactic Cosmic Ray impacts on the photoconductor arrays of the PACS instrument. This instrument is part of the ESA-Herschel payload, which will be launched in 2008 and will operate at the Lagrangian L2 point of the Sun-Earth system. Both the Satellite plus the cryostat (the shield) and the detector act as source of secondary events, affecting the detector performance. Secondary event rates originated within the detector and from the shield are of comparable intensity. The impacts deposit energy on each photoconductor pixel but do not affect the behaviour of nearby pixels. These latter are hit with a probability always lower than 7%. The energy deposited produces a spike which can be hundreds times larger than the noise. We then compare our simulations with proton irradiation tests carried out for one of the detector modules and follow the detector behaviour under ‘real’ conditions.     

Ultra High Energy Cosmic Rays: Anisotropies and spectrum

The recent results of the Pierre Auger Observatory on the possible correlation of Ultra High Energy Cosmic Rays events and several nearby discrete sources could be the starting point of a new era with charged particles astronomy. In this paper we introduce a simple model to determine the effects of any local distribution of sources on the expected flux. We consider two populations of sources: faraway sources uniformly distributed and local point sources. We study the effects on the expected flux of the local distribution of sources, referring also to the set of astrophysical objects whose correlation with the Auger events is experimentally claimed.     

On the prospects of Ultra-High Energy Cosmic Rays detection by high altitude antennas

Radio emission from Ultra-High Energy Cosmic Rays (UHECR) showers detected after specular reflection off the Antarctic ice surface has been recently demonstrated by the ANITA balloon-borne experiment. An antenna observing a large area of ice or water from a mountaintop, a balloon or a satellite may be competitive with more conventional techniques. We present an estimate of the exposure of a high altitude antenna, which provides insight on the prospects of this technique for UHECR detection. We find that a satellite antenna may reach a significantly larger exposure than existing UHECR observatories, but an experimental characterization of the radio reflected signal is required to establish the potential of this approach. A balloon-borne or a mountaintop antenna are found not to be competitive under any circumstances.     

The Lateral Trigger Probability function for the Ultra-High Energy Cosmic Ray showers detected by the Pierre Auger Observatory

In this paper we introduce the concept of Lateral Trigger Probability (LTP) function, i.e., the probability for an Extensive Air Shower (EAS) to trigger an individual detector of a ground based array as a function of distance to the shower axis, taking into account energy, mass and direction of the primary cosmic ray. We apply this concept to the surface array of the Pierre Auger Observatory consisting of a 1.5 km spaced grid of about 1600 water Cherenkov stations. Using Monte Carlo simulations of ultra-high energy showers the LTP functions are derived for energies in the range between 10¹⁷ and 10¹⁹ eV and zenith angles up to 65°. A parametrization combining a step function with an exponential is found to reproduce them very well in the considered range of energies and zenith angles. The LTP functions can also be obtained from data using events simultaneously observed by the fluorescence and the surface detector of the Pierre Auger Observatory (hybrid events). We validate the Monte Carlo results showing how LTP functions from data are in good agreement with simulations.     

Wavelet analysis of non-Gaussian anisotropies from primordial voids in simulated maps of the Cosmic Microwave Background

Phase transitions taking place during the inflationary epoch give rise to bubbles of true vacuum embedded in the false vacuum. These bubbles can imprint a distinctive signal on the Cosmic Microwave Background (CMB). We evaluate the feasibility of detecting these signatures with wavelets in CMB maps, such as those that will be made available by the European Space Agency’s (ESA) Planck mission.     

A new parallel PM code for very large-scale cosmological simulations

We have developed a parallel Particle–Particle, Particle–Mesh (P³M) simulation code for the Cray T3E parallel supercomputer that is well suited to studying the time evolution of systems of particles interacting via gravity and gas forces in cosmological contexts. The parallel code is based upon the public-domain serial Adaptive P³M-SPH (http://coho.astro.uwo.ca/pub/hydra/hydra.html) code of Couchman et al. (1995)[ApJ, 452, 797]. The algorithm resolves gravitational forces into a long-range component computed by discretizing the mass distribution and solving Poisson's equation on a grid using an FFT convolution method, and a short-range component computed by direct force summation for sufficiently close particle pairs. The code consists primarily of a particle–particle computation parallelized by domain decomposition over blocks of neighbour-cells, a more regular mesh calculation distributed in planes along one dimension, and several transformations between the two distributions. The load balancing of the P³M code is static, since this greatly aids the ongoing implementation of parallel adaptive refinements of the particle and mesh systems. Great care was taken throughout to make optimal use of the available memory, so that a version of the current implementation has been used to simulate systems of up to 10⁹ particles with a 1024³ mesh for the long-range force computation. These are the largest Cosmological N-body simulations of which we are aware. We discuss these memory optimizations as well as those motivated by computational performance. Performance results are very encouraging, and, even without refinements, the code has been used effectively for simulations in which the particle distribution becomes highly clustered as well as for other non-uniform systems of astrophysical interest.     

A numerical scheme and benchmark tests for non-isothermal two-fluid ambipolar diffusion

Stars form in magnetized molecular clouds composed primarily of neutral gas with a trace amount of ions. We present a semi-implicit strategy for incorporating the equations that describe the coupled ion and neutral two-fluid equations, with a full energy equation, into the RIEMANN code that uses a TR-BDF2 algorithm to stably handle the stiffness of the source terms. We demonstrate that the numerical implementation works through the use of a suite of test problems that we catalog here. We show that reproducing the analytic dispersion analysis for the propagation of waves in a two-fluid plasma is an especially strong code test. We also present a two-fluid analogue of the Noh wall-shock problem and demonstrate the performance of the code on the Wardle instability. We also present a novel blast wave test, showing that the results reduce to the single fluid results under strong coupling, yet differing considerably when the coupling is weak. These test problems demonstrate that the numerical implementation can accurately capture the dissipation rate of waves and reproduce the structure of a C-shock.     

对脉冲星发射宇宙线能量的估计

S.Karakula等人根据Ostriker等人的脉冲星产生的宇宙线能谱依赖脉冲星年龄的理论模型,利用脉冲星发射宇宙线的最大能量:     

Bayesian foreground analysis with CMB data

The quality of CMB observations has improved dramatically in the last few years, and will continue to do so in the coming decade. Over a wide range of angular scales, the uncertainty due to instrumental noise is now small compared to the cosmic variance. One may claim with some justification that we have entered the era of precision CMB cosmology. However, some caution is still warranted: The errors due to residual foreground contamination in the CMB power spectrum and cosmological parameters remain largely unquantified, and the effect of these errors on important cosmological parameters such as the optical depth τ and spectral index n_s is not obvious. A major goal for current CMB analysis efforts must therefore be to develop methods that allow us to propagate such uncertainties from the raw data through to the final products. Here we review a recently proposed method that may be a first step towards that goal.     

Cosmic rays in the heliosphere: Observations

This contribution to the 100th commemoration of the discovery of cosmic rays (6–8 August, 2012 in Bad Saarow, Germany) is about observations of those cosmic rays that are sensitive to the structure and the dynamics of the heliosphere. This places them in the energy range of 10⁷–10¹⁰ eV. For higher energies the heliosphere becomes transparent; below this energy range the particles become strictly locked into the solar wind. Rather than give a strict chronological development, the paper is divided into distinct topics. It starts with the Pioneer/Voyager missions to the outer edges of the heliosphere, because the most recent observations indicate that a distinct boundary of the heliosphere might have been reached at the time of the meeting. Thereafter, the Ulysses mission is described as a unique one because it is still the only spacecraft that has explored the heliosphere at very high latitudes. Next, anomalous cosmic rays, discovered in 1972–1974, constitute a separate component that is ideally suited to study the acceleration and transport of energetic particles in the heliosphere. At this point the history and development of ground-based observations is discussed, with its unique contribution to supply a stable, long-term record. The last topic is about solar energetic particles with energies up to ∼10¹⁰ eV.     

PICsar: A 2.5D axisymmetric,relativistic, electromagnetic,Particle in Cell code with a radiation absorbing boundary

We present PICsar – a new Particle in Cell code geared towards efficiently simulating the magnetosphere of the aligned rotator. PICsar is a special relativistic, electromagnetic, charge conservative code that can be used to simulate arbitrary electromagnetics problems in axisymmetry. It features stretchable body-fitted coordinates that follow the surface of a sphere, simplifying the application of boundary conditions in the case of the aligned rotator; a radiation absorbing outer boundary, which allows a steady state to be set up dynamically and maintained indefinitely from transient initial conditions; and algorithms for injection of charged particles into the simulation domain. The code is parallelized using MPI and scales well to a large number of processors. We discuss the numerical methods used in PICsar and present tests of the code. In particular, we show that PICsar can accurately and efficiently simulate the magnetosphere of the aligned monopole rotator in the force-free limit. We present simulations of the aligned dipole rotator in a forthcoming paper.     

Helioseismology program for the PICARD satellite

The PICARD mission is a CNES micro‐satellite to be launched in 2009. Its goal is to better understand the Sun and the potential impact of its activity on earth climate by measuring simultaneously the solar total and spectral irradiance, diameter, shape and oscillations. We present the scientific objectives, instrumental requirements and data products of the helioseismology program of PICARD which aims to observe the low to medium l p‐mode oscillations in intensity and search for g‐mode oscillation signatures at the limb. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The Cosmic Foreground Explorer (COFE): A balloon-borne microwave polarimeter to characterize polarized foregrounds

The COsmic Foreground Explorer (COFE) is a balloon-borne microwave polarimeter designed to measure the low-frequency and low-ℓ characteristics of dominant diffuse polarized foregrounds. Short duration balloon flights from the Northern and Southern Hemispheres will allow the telescope to cover up to 80% of the sky with an expected sensitivity per pixel better than 100 μK/deg² from 10 GHz to 20 GHz. This is an important effort toward characterizing the polarized foregrounds for future CMB experiments, in particular the ones that aim to detect primordial gravity wave signatures in the CMB polarization angular power spectrum.     

X‐Ray observations of the Galactic Centre with Suzaku

We report on the diffuse X‐ray emission from the Galactic Centre (GCDX) observed with the X‐ray Imaging Spectrometer (XIS) on board the Suzaku satellite. The highly accurate energy calibration and extremely low background of the XIS provide many new facts on the GCDX. These are (1) the origin of the lines at 6.7 and 7.0 keV is collisional excitation in a hot plasma, (2) the discovery of new SNR and super‐bubble candidates, (3) most of the 6.4 keV line is X‐ray fluorescence, and (4) time variability of the 6.4 keV line is found from the Sgr B2 complex. (© 2008 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

A parallel TreeSPH code for galaxy formation

We describe a new implementation of a parallel TreeSPH code with the aim of simulating galaxy formation and evolution. The code has been parallelized using shmem , a Cray proprietary library to handle communications between the 256 processors of the Silicon Graphics T3E massively parallel supercomputer hosted by the Cineca Super-computing Center (Bologna, Italy). ¹
The code combines the smoothed particle hydrodynamics (SPH) method for solving hydrodynamical equations with the popular Barnes & Hut tree-code to perform gravity calculation with an N ×log  N scaling, and it is based on the scalar TreeSPH code developed by Carraro et al. Parallelization is achieved by distributing particles along processors according to a workload criterion.
Benchmarks, in terms of load balance and scalability, of the code are analysed and critically discussed against the adiabatic collapse of an isothermal gas sphere test using 2×10⁴ particles on 8 processors. The code results balance at more than the 95 per cent level. Increasing the number of processors, the load balance slightly worsens. The deviation from perfect scalability for increasing number of processors is almost negligible up to 32 processors. Finally, we present a simulation of the formation of an X-ray galaxy cluster in a flat cold dark matter cosmology, using 2×10⁵ particles and 32 processors, and compare our results with Evrard's P³M–SPH simulations.
Additionally we have incorporated radiative cooling, star formation, feedback from SNe of types II and Ia, stellar winds and UV flux from massive stars, and an algorithm to follow the chemical enrichment of the interstellar medium. Simulations with some of these ingredients are also presented.     

TreePM: A code for cosmological N-body simulations

We describe the TreePM method for carrying out large N-Body simulations to study formation and evolution of the large scale structure in the Universe. This method is a combination of Barnes and Hut tree code and Particle-Mesh code. It combines the automatic inclusion of periodic boundary conditions of PM simulations with the high resolution of tree codes. This is done by splitting the gravitational force into a short range and a long range component. We describe the splitting of force between these two parts. We outline the key differences between TreePM and some other N-Body methods.     

Cosmic rays and terrestrial life: A brief review

“The investigation into the possible effects of cosmic rays on living organisms will also offer great interest.” – Victor F. Hess, Nobel Lecture, December 12, 1936High-energy radiation bursts are commonplace in our Universe. From nearby solar flares to distant gamma ray bursts, a variety of physical processes accelerate charged particles to a wide range of energies, which subsequently reach the Earth. Such particles contribute to a number of physical processes occurring in the Earth system. A large fraction of the energy of charged particles gets deposited in the atmosphere, ionizing it, causing changes in its chemistry and affecting the global electric circuit. Remaining secondary particles contribute to the background dose of cosmic rays on the surface and parts of the subsurface region. Life has evolved over the past ∼3 billion years in presence of this background radiation, which itself has varied considerably during the period [1], [2], [3]. As demonstrated by the Miller–Urey experiment, lightning plays a very important role in the formation of complex organic molecules, which are the building blocks of more complex structures forming life. There is growing evidence of increase in the lightning rate with increasing flux of charged particles. Is there a connection between enhanced rate of cosmic rays and the origin of life? Cosmic ray secondaries are also known to damage DNA and cause mutations, leading to cancer and other diseases. It is now possible to compute radiation doses from secondary particles, in particular muons and neutrons. Have the variations in cosmic ray flux affected the evolution of life on earth? We describe the mechanisms of cosmic rays affecting terrestrial life and review the potential implications of the variation of high-energy astrophysical radiation on the history of life on earth.