Validation of Geant4 10.3 simulation of proton interaction for space radiation effects

Monte Carlo simulation of space radiation effects induced by protons is important for design of space missions. Geant4 is a well established toolkit for Monte Carlo simulation focused on high energy physics applications. In this work, a set of new validation results versus data for Geant4 electromagnetic and hadronic interaction of protons is presented and discussed. Optimal configuration of Geant4 physics for space applications is proposed.     

Proton flux and radiation dose from galactic cosmic rays in the lunar regolith and implications for organic synthesis at the poles of the Moon and Mercury

Galactic cosmic rays are a potential energy source to stimulate organic synthesis from simple ices. The recent detection of organic molecules at the polar regions of the Moon by LCROSS (Colaprete, A. et al. [2010]. Science 330, 463–468, http://dx.doi.org/10.1126/science.1186986), and possibly at the poles of Mercury (Paige, D.A. et al. [2013]. Science 339, 300–303, http://dx.doi.org/10.1126/science.1231106), introduces the question of whether the organics were delivered by impact or formed in situ. Laboratory experiments show that high energy particles can cause organic production from simple ices. We use a Monte Carlo particle scattering code (MCNPX) to model and report the flux of GCR protons at the surface of the Moon and report radiation dose rates and absorbed doses at the Moon’s surface and with depth as a result of GCR protons and secondary particles, and apply scaling factors to account for contributions to dose from heavier ions. We compare our results with dose rate measurements by the Cosmic Ray Telescope for the Effects of Radiation (CRaTER) experiment on Lunar Reconnaissance Orbiter (Schwadron, N.A. et al. [2012]. J. Geophys. Res. 117, E00H13, http://dx.doi.org/10.1029/2011JE003978) and find them in good agreement, indicating that MCNPX can be confidently applied to studies of radiation dose at and within the surface of the Moon. We use our dose rate calculations to conclude that organic synthesis is plausible well within the age of the lunar polar cold traps, and that organics detected at the poles of the Moon may have been produced in situ. Our dose rate calculations also indicate that galactic cosmic rays can induce organic synthesis within the estimated age of the dark deposits at the pole of Mercury that may contain organics.     

Dynamic Monte Carlo radiation transfer in SPH: radiation pressure force implementation

We present a new framework for radiation hydrodynamics simulations. Gas dynamics is modelled by smoothed particle hydrodynamics (SPH), whereas radiation transfer is simulated via a time-dependent Monte Carlo approach that traces photon packets. As a first step in the development of the method, in this paper we consider the momentum transfer between radiation field and gas, which is important for systems where radiation pressure is high. There is no fundamental limitation on the number of radiation sources, the geometry or the optical depth of the problems that can be studied with the method. However, as expected for any Monte Carlo transfer scheme, stochastic noise presents a serious limitation. We present a number of tests that show that the errors of the method can be estimated accurately by considering Poisson noise fluctuations in the number of photon packets that SPH particles interact with per dynamical time. It is found that, for a reasonable accuracy, the momentum carried by photon packets must be much smaller than the typical momentum of SPH particles. We discuss numerical limitations of the code, and future steps that can be taken to improve performance and applicability of the method.     

A Moon orbiting observatory for the measurement of the elemental and isotopic composition in the primary cosmic radiation

We describe a project for the measurement of elemental composition of the primary cosmic radiation to be performed by a space observatory orbiting around the Moon. The absence of atmosphere and the low intrinsic magnetic field of the Moon give access to the very low energy component of the cosmic radiation, allowing the search for rare events.The quest for antinuclei, the determination of the lunar lepton albedo and the abundance measurement of galactic radioactive clocks (Be¹⁰, C¹⁴, Al²⁶) are the major tasks of the ANTARES apparatus (ANTimatter Assessment RESearch).We report details of the instrument design, the expected performance for single detectors, their capability to accomplish the proposed measurements and the characteristics of the space mission.     

Circumstellar dust and circumstellar extinction curves

We have calculated the circumstellar extinction curves produced by dust grains which absorb and scatter the stellar radiation in the shells of pre-main-sequence stars. A Monte Carlo method was used to model the radiative transfer in non-spherical shells. The dependence on the particle size distribution and the dust shell parameters has been examined.The application of the theoretical results to explain the extinction and polarization of the Herbig Be star HD 45677 shows that the dust shell is not disk-like and that very small grains are absent in it.     

Hard X-ray emission of the Earth's atmosphere: Monte Carlo simulations

We perform Monte Carlo simulations of cosmic ray-induced hard X-ray radiation from the Earth's atmosphere. We find that the shape of the spectrum emergent from the atmosphere in the energy range 25–300 keV is mainly determined by Compton scatterings and photoabsorption, and is almost insensitive to the incident cosmic ray spectrum. We provide a fitting formula for the hard X-ray surface brightness of the atmosphere as would be measured by a satellite-borne instrument, as a function of energy, solar modulation level, geomagnetic cut-off rigidity and zenith angle. A recent measurement by the INTEGRAL observatory of the atmospheric hard X-ray flux during the occultation of the cosmic X-ray background by the Earth agrees with our prediction within 10 per cent. This suggests that Earth observations could be used for in-orbit calibration of future hard X-ray telescopes. We also demonstrate that the hard X-ray spectra generated by cosmic rays in the crusts of the Moon, Mars and Mercury should be significantly different from that emitted by the Earth's atmosphere.     

Multidimensional simulations of radiative transfer in Type Ia supernovae

A three-dimensional Monte Carlo code for modelling radiation transport in Type Ia supernovae is described. In addition to tracking Monte Carlo quanta to follow the emission, scattering and deposition of radiative energy, a scheme involving volume-based Monte Carlo estimators is used to allow properties of the emergent radiation field to be extracted for specific viewing angles in a multidimensional structure. This eliminates the need to compute spectra or light curves by angular binning of emergent quanta. The code is applied to two test problems to illustrate consequences of multidimensional structure on the modelling of light curves. First, elliptical models are used to quantify how large-scale asphericity can introduce angular dependence to light curves. Secondly, a model which incorporates complex structural inhomogeneity, as predicted by modern explosion models, is used to investigate how such structure may affect light-curve properties.     

On Choosing a Rational Flight Trajectory to the Moon

The algorithm for choosing a trajectory of spacecraft flight to the Moon is discussed. The characteristic velocity values needed for correcting the flight trajectory and a braking maneuver are estimated using the Monte Carlo method. The profile of insertion and flight to a near-circular polar orbit with an altitude of ~100 km of an artificial lunar satellite (ALS) is given. The case of two corrections applied during the flight and braking phases is considered. The flight to an ALS orbit is modeled in the geocentric geoequatorial nonrotating coordinate system with the influence of perturbations from the Earth, the Sun, and the Moon factored in. The characteristic correction costs corresponding to corrections performed at different time points are examined. Insertion phase errors, the errors of performing the needed corrections, and the errors of determining the flight trajectory parameters are taken into account.     

Detailed calculation of test-mass charging in the LISA mission

The electrostatic charging of the LISA test masses due to exposure of the spacecraft to energetic particles in the space environment has implications in the design and operation of the gravitational inertial sensors and can affect the quality of the science data. Robust predictions of charging rates and associated stochastic fluctuations are therefore required for the exposure scenarios expected throughout the mission. We report on detailed charging simulations with the Geant4 toolkit, using comprehensive geometry and physics models, for Galactic cosmic-ray protons and helium nuclei. These predict positive charging rates of 50+e/s (elementary charges per second) for solar minimum conditions, decreasing by half at solar maximum, and current fluctuations of up to 30+e/s/Hz^1/2. Charging from sporadic solar events involving energetic protons was also investigated. Using an event-size distribution model, we conclude that their impact on the LISA science data is manageable. Several physical processes hitherto unexplored as potential charging mechanisms have also been assessed. Significantly, the kinetic emission of very low-energy secondary electrons due to bombardment of the inertial sensors by primary cosmic rays and their secondaries can produce charging currents comparable with the Monte Carlo rates.     

Early scattering by Jupiter and its collision effects in the terrestrial zone

When Jupiter was on the order of three to ten Earth masses in size, there undoubtedly was a considerably larger mass of condensed matter in its zone, since Jupiter would have perturbed most of it to other parts of the solar system. Monte Carlo studies indicate a significant portion would have crossed the Earth's orbit. If the Earth and Moon had not yet fully formed, the probability of Earth-zone planetesimals being hit by this Jupiter-scattered material was high. Further Monte Carlo models of these collisions and their products indicate a significant portion of matter was heated to melting, even if less than 5% of the relative kinetic energy went into heat. The models include capture probabilities by an embryo Earth and a protolunar swarm. Because heat energy is correlated with comminution energy, and because the capture probability of the swarm is mass-dependent while the embryo's is not, the protolunar material suffered much higher heating on the average than did the proto-Earth material.     

The Application Of À Trous Wave Filtering And Monte Carlo Analysis On Secis 2001 Solar Eclipse Observations

Eight thousand images of the solar corona were captured during the June 2001 total solar eclipse. New software for the alignment of the images and an automated technique for detecting intensity oscillations using multi-scale wavelet analysis were developed. Large areas of the images covered by the Moon and the upper corona were scanned for oscillations and the statistical properties of the atmospheric effects were determined. The à Trous wavelet transform was used for noise reduction and Monte Carlo analysis as a significance test of the detections. The effectiveness of those techniques is discussed in detail.     

New Horizons: The First Reconnaissance Mission to Bodies in the Kuiper Belt

NASA has long been planning a mission of exploration to Pluto-Charon and the Kuiper Belt (e.g., Terrile et al., 1997). In 2001 NASA selected such a mission (NASA, 2001), called New Horizons, for design and development. New Horizons is now funded and planning a launch in January 2006. The mission plans to carry 8 scientific sensors and make flybys of Pluto-Charon and one or more KBOs. Statistical Monte Carlo simulations indicate that New Horizons has sufficient fuel to reach one or more KBOs with diameters exceeding 35 km. If launched as planned in 2006, the mission will use a Jovian gravity assist, arriving at Pluto-Charon in 2015 or 2016; if launched in its backup window in 2007, a Jovian gravity assist is not feasible and arrival will be later – 2019. Below we briefly summarize the New Horizons mission, concentrating on its role in Kuiper Belt exploration.     

A model of the cosmic ray induced atmospheric neutron environment

In order to optimise the design of space instruments making use of detection materials with low atomic numbers, an understanding of the atmospheric neutron environment and its dependencies on time and position is needed. To produce a simple equation based model, Monte Carlo simulations were performed to obtain the atmospheric neutron fluxes produced by charged galactic cosmic ray interactions with the atmosphere. Based on the simulation results the omnidirectional neutron environment was parametrized including dependencies on altitude, magnetic latitude and solar activity. The upward- and downward-moving component of the atmospheric neutron flux are considered separately. The energy spectra calculated using these equations were found to be in good agreement with data from a purpose built balloon-borne neutron detector, high altitude aircraft data and previously published simulation based spectra.     

Long-term evolution of Oort Cloud comets: methods and comparisons

Two long-term simulation methods for cometary orbits, a Monte Carlo method and a direct integration method, are compared with each other. The comparison is done in seven inclination and perihelion distance intervals, and shows differences in dynamical lifetime and capture probabilities for the following main reasons. We use a finite energy step approximation in the Monte Carlo method and the method considers only close approaches with the planets. The differences can be taken into account statistically and it is possible to calculate the correction factors for the capture probability and dynamical lifetime in the Monte Carlo method. Both corrections depend on the inclination and on the value of the minimum energy step. The capture probabilities of the short-period comets originating in the Oort Cloud are calculated by the corrected Monte Carlo method and compared with published results.     

The Effect of Method of Constructing the Mass Distribution on Single Stellar Populations

We use two methods of constructing the initial mass distribution, the traditional way and Monte Carlo simulation, to obtain integrated U - B, B - V, V-R and V-I colours and absorption-line indices denned by the Lick Observatory image dissector scanner (referred to as Lick/IDS), for instantaneous burst solarmetallicity single stellar populations with ages in the range 1-15 Gyr. We find that the evolutionary curves of all colours obtained by the traditional method are smoother than those by Monte Carlo simulation, that the U - B and B - V colours obtained by the two methods agree with one another, while the V - R and V - I colours by the traditional method are bluer than those by Monte Carlo simulation. A comparison of the Lick/IDS absorption-line indices shows that the variations in all the indices by the traditional method are smoother than that for the Monte Carlo simulation, and that all the indices except for TiO1 and TiO2 are consistent with those for the Monte Carlo simulation.     

Monte Carlo simulation of pair creation using petawatt lasers

Irradiating high-Z targets such as gold with ultra-intense lasers creates electron-positron pairs. In particular, the positron density in the plasma created by this procedure is higher than that obtained via other laboratory-based methods, with theoretical maximum densities exceeding 10¹⁸ cm⁻³. All of the significantly contributing processes are well-known and hence we can study this phenomenon using Monte Carlo simulation. We focus on the latter part of this procedure, the passage of high-energy electrons through the target creating pairs. In particular, we discuss the usefulness of CERN’s GEANT4 Monte Carlo code in simulating this process. Once this code is successfully calibrated, we will use it to perform parameter studies, and design future targets to optimize the positron yield.     

The effect of surface roughness on the transmission of microwave radiation through a planetary surface

To account for surface roughness, the transmission of microwave radiation through a planetary surface to an observer is treated by a Monte Carlo technique. Sizable effects are found near the limb of the planet, and they should be included in analyses of high-resolution observations and high-precision integrated disk observations.     

A new method to reconstruct the energy and determine the composition of cosmic rays from the measurement of Cherenkov light and particle densities in extensive air showers

A Monte Carlo study to reconstruct energy and mass of cosmic rays with energies above 300 TeV using ground based measurements of the electromagnetic part of showers initiated in the atmosphere is presented. The shower properties determined with two detector arrays measuring the air Cherenkov light and the particle densities as realized at the HEGRA experiment are processed to determine the energy of the primary particle without the need of any hypothesis concerning its mass. The mass of the primary particle is reconstructed coarsely from the same observables in parallel to the energy determination.     

Indirect transfer to the Earth–Moon L1 libration point

The Earth–Moon L1 libration point is proposed as a human gateway for space transportation system of the future. This paper studies indirect transfer using the perturbed stable manifold and lunar flyby to the Earth–Moon L1 libration point. Although traditional studies indicate that indirect transfer to the Earth–Moon L1 libration point does not save much fuel, this study shows that energy efficient indirect transfer using the perturbed stable manifold and lunar flyby could be constructed in an elegant way. The design process is given to construct indirect transfer to the Earth–Moon L1 libration point. Simulation results show that indirect transfer to the Earth–Moon L1 libration point saves about 420 m/s maneuver velocity compared to direct transfer, although the flight time is about 20 days longer.