An update about recent developments of the PHITS code

PHITS (Particle and Heavy-Ion Transport code System) is a general-purpose three-dimensional Monte Carlo code, developed and maintained by RIST, JAEA and KEK in Japan together with Sihver et al. at Chalmers in Sweden. PHITS can deal with the transports of all varieties of hadrons and heavy ions with energies up to around 100 GeV/nucleon, and in this paper the current status of PHITS is presented. We introduce a relativistically covariant version of JQMD, called R-JQMD, that features an improved ground state initialization algorithm, and we will present the introduction of electron and photon transport in PHITS using EGS5, which have increased the energy region for the photon and energy transport from up to around 3 GeV to up to several hundred GeV depending on the atomic number of the target. We show how the accuracy in dose and fluence calculations can be improved by using tabulated cross sections. Benchmarking of shielding and irradiation effects of high energy protons in different materials relevant for shielding of accelerator facilities is also presented. In particular, we show that PHITS can be used for estimating the dose received by aircrews and personnel in space. In recent years, many countries have issued regulations or recommendations to set annual dose limitations for aircrews. Since estimation of cosmic-ray spectra in the atmosphere is an essential issue for the evaluation of aviation doses, we have calculated these spectra using PHITS. The accuracy of the atmospheric propagation simulation of cosmic-ray performed by PHITS has been well verified by experimental cosmic-ray spectra taken under various conditions. Based on a comprehensive analysis of the simulation results, an analytical model called “PARMA” has been proposed for instantaneously estimating the atmospheric cosmic-ray spectra below the altitude of 20 km. We have also performed preliminary simulations of long-term dose distribution measurements at the ISS performed with the joint ESA-FSA experiment MATROSHKA-R (MTR-R) led by the Russian Federation Institute of Biomedical Problems (IMBP) and the ESA supported experiment MATROSHKA (MTR), led by the German Aerospace Center (DLR). For the purpose of examining the applicability of PHITS to the shielding design in space, the absorbed doses in a tissue equivalent water phantom inside an imaginary space vessel has been estimated for different shielding materials of different thicknesses. The results confirm previous results which indicate that PHITS is a suitable tool when performing shielding design studies of spacecrafts.     

Physics of ion acceleration in the solar flare on 2005 September 7 determines γ-ray and neutron production

Relativistic neutrons were observed by the neutron monitors at Mt. Chacaltaya and Mexico City and by the solar neutron telescopes at Chacaltaya and Mt. Sierra Negra in association with an X17.0 flare on 2005 September 7. The neutron signal continued for more than 20 min with high statistical significance. Intense emissions of γ-rays were also registered by INTEGRAL, and during the decay phase by RHESSI. We analyzed these data using the solar-flare magnetic-loop transport and interaction model of Hua et al. [Hua, X.-M., Kozlovsky, B., Lingenfelter, R.E. et al. Angular and energy-dependent neutron emission from solar flare magnetic loops, Astrophys. J. Suppl. Ser. 140, 563–579, 2002], and found that the model could successfully fit the data with intermediate values of loop magnetic convergence and pitch-angle scattering parameters. These results indicate that solar neutrons were produced at the same time as the γ-ray line emission and that ions were continuously accelerated at the emission site.     

The mass spectrum analyzer (MSA) onboard BEPI COLOMBO MMO: Scientific objectives and prototype results

BEPI COLOMBO is a joint mission between ESA and JAXA that is scheduled for launch in 2014 and arrival at Mercury in 2020. A comprehensive set of ion sensors will be flown onboard the two probes that form BEPI COLOMBO. These ion sensors combined with electron analyzers will allow a detailed investigation of the structure and dynamics of the charged particle environment at Mercury. Among the ion sensors, the Mass Spectrum Analyzer (MSA) is the experiment dedicated to composition analysis onboard the Mercury Magnetospheric Orbiter (MMO). It consists of a top-hat for energy analysis followed by a Time-Of-Flight (TOF) section to derive the ion mass. A notable feature of MSA is that the TOF section is polarized with a linear electric field that provides an enhanced mass resolution, a capability that is of importance at Mercury since a variety of species originating from the planet surface and exosphere is expected. MSA exhibits two detection planes: (i) one with moderate mass resolution but a high count rate making MSA appropriate for plasma analysis, (ii) another with a high (above 40) mass resolution though a low count rate making it appropriate for planetology science. Taking advantage of the spacecraft rotation, MSA will provide three-dimensional distribution functions of magnetospheric ions, from energies characteristic of exospheric populations (a few eVs or a few tens of eVs) up to the plasma sheet energy range (up to ∼40 keV/q) in one spin (4 s).     

Correlation between cell death and induction of non-rejoining PCC breaks by carbon-ion beams.

We have shown a correlation between cell death and induction of non-rejoining chromatin breaks in two normal human cells and three human tumor cell lines irradiated by carbon-ion beams and X rays. Non-rejoining chromatin breaks were measured by counting the number of remaining chromatin fragments detected by the premature chromosome condensation (PCC) technique. Carbon-ion beams were accelerated by the Heavy Ion Medical Accelerator in Chiba (HIMAC). The cells were irradiated by two different mono-LET beams (LET = 13 keV/micrometer and 77 keV/micrometer ) and 200 kV X rays. The RBE values of cell death for carbon-ion beams relative to X rays were 1.1 to 1.4 for 13 keV/micrometer beams and 2.5 to 2.9 for 77 keV/micrometer beams. The induction rate of non-rejoining PCC breaks per cell per Gy was found to be highest for the 77 keV/micrometer beams for all of the cell lines.The results found in this study show that there is a good correlation between cell death and induction of non-rejoining PCC breaks for these human cell lines.     

Tolerance of anhydrobiotic eggs of the Tardigrade Ramazzottius varieornatus to extreme environments

Tardigrades are tiny (less than 1?mm in length) invertebrate animals that have the potential to survive travel to other planets because of their tolerance to extreme environmental conditions by means of a dry ametabolic state called anhydrobiosis. While the tolerance of adult tardigrades to extreme environments has been reported, there are few reports on the tolerance of their eggs. We examined the ability of hydrated and anhydrobiotic eggs of the tardigrade Ramazzottius varieornatus to hatch after exposure to ionizing irradiation (helium ions), extremely low and high temperatures, and high vacuum. We previously reported that there was a similar pattern of tolerance against ionizing radiation between hydrated and anhydrobiotic adults. In contrast, anhydrobiotic eggs (50% lethal dose; 1690 Gy) were substantially more radioresistant than hydrated ones (50% lethal dose; 509 Gy). Anhydrobiotic eggs also have a broader temperature resistance compared with hydrated ones. Over 70% of the anhydrobiotic eggs treated at either -196°C or +50°C hatched successfully, but all the hydrated eggs failed to hatch. After exposure to high-vacuum conditions (5.3×10(-4) Pa to 6.2×10(-5) Pa), the hatchability of the anhydrobiotic eggs was comparable to that of untreated control eggs.     

Method of the HF wave absorption evaluation based on GIRO network data

High frequency (HF) communication is strongly dependent on the state of the ionosphere, which specifies the mode structure of the radio wave propagating in ionosphere. Another core factor defining the strength of the HF signal at the receiving site is the ionospheric absorption. Accurate modelling the effect of absorption is an essential part of many studies of the HF propagation in the ionosphere.This paper proposes a method for estimating the absorption. The method is based on analysis of vertical sounding ionograms. The main idea of the approach is to compare the main parameters retrieved from measured and simulated ionograms. The combination of Global Ionospheric Radio Observatory (GIRO, http://giro.uml.edu) data and ionograms modelling allows for developing the empiric absorption model available at near real-time. The ionogram simulation taking into account absorption utilizes the NIM-RT (North Ionospheric Model and Ray Tracing) software. As a result, the proposed technique provides more reliable and accurate evaluation of minimum frequency at which echoes are observed in vertical incidence ionosonde soundings. The values of these frequencies should be used in the following simulation to optimize parameters in the empirical formulae for defining absorption HF wave in ionosphere.The ultimate objective of this work is the designing the method, which allows the simulating of HF radio channel accounting for regular absorption due to UV radiation of the Sun. Eventually it could be considered as some kind of the HF propagation forecasting.     

Alpha Particle X-Ray Spectrometer (APXS) on-board Chandrayaan-2 rover

Alpha Particle X-ray Spectrometer (APXS) payload configuration for Chandrayaan-2 rover has been completed recently and fabrication of mechanical assembly, PCB layout design and fabrication are in progress. Here we present the design and performance evaluation of various subsystems developed for APXS payload. The low energy threshold of <1 keV and the energy resolution of ∼150 eV at 5.9 keV, for the Silicon Drift Detector (SDD), as measured from the developed APXS electronics is comparable to the standard spectrometers available off-the-shelf. We have also carried out experiments for measuring fluorescent X-ray spectrum from various standard samples from the USGS catalog irradiated by the laboratory X-ray source ²⁴¹Am with 1 mCi activity. It is shown that intensities of various characteristic X-ray lines are well correlated with the respective elemental concentrations.     

Is there an abnormal enhancement of atmospheric aerosol before the 2008 Wenchuan earthquake?

There are extensive reports of ionospheric disturbances before the great 2008 Wenchuan earthquake, which are possibly explained by seismogenic electric field hypotheses linked with the aerosols injected in atmosphere. This paper attempts to investigate the possible change of atmospheric aerosol optical depth (AOD) associated with this earthquake by using MODIS data from both Terra and Aqua satellites. The result shows a clear enhancement of AOD along the Longmenshan faults 7 days before the quake, which is 1 day and 4 days earlier than the reported negative and positive ionospheric disturbances, respectively, and is 1 day earlier than or quasi-synchronism with other reported atmospheric anomalies including air temperature, outgoing longwave radiation and relative humidity. Particularly, the spatial distribution of AOD enhancement is very local and it is correlated well with the active faults and surface ruptures. We suggest that this unique enhancement could be associated with the Lithosphere–Atmosphere–Ionosphere coupling process during the preparation of the Wenchuan earthquake.     

Lattice Boltzmann Flux Solver： An Efficient Approach for Numerical Simulation of Fluid Flows

A lattice Boltzmann flux solver （LBFS） is presented for simulation of fluid flows. Like the conventional computational fluid dynamics （CFD） solvers, the new solver also applies the finite volume method to discretize the governing differential equations, but the numerical flux at the cell interface is not evaluated by the smooth function approximation or Riemann solvers. Instead, it is evaluated from local solution of lattice Boltzmann equation （LBE） at cell interface. Two versions of LBFS are presented in this paper. One is to locally apply one-dimensional compressible lattice Boltzmann （LB） model along the normal direction to the cell interface for simulation of compressible inviscid flows with shock waves. The other is to locally apply multi-dimensional LB model at cell interface for simulation of incompressible viscous and inviscid flows. The present solver removes the drawbacks of conventional lattice Boltzmann method （LBM） such as limitation to uniform mesh, tie-up of mesh spacing and time interval, limitation to viscous flows. Numerical examples show that the present solver can be well applied to simulate fluid flows with non-uniform mesh and curved boundary.