Ab Initio Prediction of a Negative Barkas Coefficient for Slow Protons and Antiprotons in LiF

We report the ab initio prediction of a negative Barkas coefficient in lithium fluoride (LiF) insulator at low velocity (v<0.25  a.u., E_{kin}∼2  keV). The electronic stopping power of protons in LiF has been extensively studied both experimentally and theoretically because of a controversial threshold effect. While our time-dependent density-functional theory simulations confirm the presence of a velocity threshold below which the proton stopping power vanishes, our calculations demonstrate that the antiprotons do not experience such a threshold. The combination of those two contrasting behaviors gives rise to an unprecedented negative Barkas effect: the stopping power of antiprotons is larger than that of protons. We identify that the slow antiproton at close encounter destabilizes a p orbital of the F^{-} anion pointing toward the antiproton. This particular orbital becomes highly polarizable and hence contributes much to the stopping power.

Use of Gaussian-type functions for flux-based dose calculations in carbon ion therapy

In radiation therapy, it is very important to ensure that the radiation dose is correctly delivered to the patient. This is achieved by obtaining quantitative dose measurements for beam calibration in the treatment planning system. Dose calculations should be performed with the required accuracy to a degree of uncertainty of less than 1%. The measurement of the absorbed dose in and around body tissues irradiated with carbon ions requires careful use of materials selected from established phantom and radiation detectors. The main advantage of such materials is that when information on the energy and nature of charged particles at the desired point is incomplete or inaccurate, they can allow determination of the absorbed dose. In general, radiation interactions in a tissue representation caused by carbon ions can be characterized by calculating the linear stopping power. Carbon ions have a limited penetration depth within human tissues that depends on the energy and stopping power of these ions as they penetrate into the body. The purpose of the present study was to calculate the stopping power, range and dose to intestinal and prostate tissues of carbon ions. The stopping power values of these tissues were specified by the effective charge approach method. The 5ZaPa-NR-CV, pcemd-4 and pcSseg-4 sets of Gaussian-type functions were employed for the calculation of electronic charge density. Range calculations were made by means of the Gaussian quadrature method, making use of the continuous slowing down approximation. Flux-based dose calculations were also carried out in accordance with the Bragg-Gray theorem using the Geant4 and FLUKA simulation toolkits. The results were compared with each other and with the SRIM and CasP datasets. Then, depth-dose distributions and range values were verified by positron emission activity using the GATE toolkit. Among the different types of Gaussian functions used here, the best semi-analytical result was found for the 5ZaPa-NR-CV set. The results obtained in the present study can be used for dose verification and dose reconstruction in charged particle radiotherapy and for radiation research on the interaction of radiation with matter. The results calculated here will be useful for quantifying uncertainties associated with stopping power, range, and reconstruction of dose in charged particle therapy.

First-Principles Modeling of Electronic Stopping in Complex Matter under Ion Irradiation

Electronic stopping refers to the dynamical energy-transfer process to electrons in matter from highly energetic charged particles such as high-velocity protons. We discuss recent progress in theoretical studies of electronic stopping in condensed matter under ion irradiation, focusing on modern electronic structure theory's role in enabling the study of electronic excitation dynamics that result from the energy transfer. In the last few decades, first-principles simulation approaches based on real-time time-dependent density functional theory have greatly advanced the field. While linear response theory is widely used to study electronic stopping processes, especially for simple solids, novel first-principles dynamics approaches now allow us to study chemically complex systems and also yield detailed descriptions of electronic excitations at the molecular scale. Outstanding challenges for further advancement of electronic stopping modeling are also discussed from the viewpoint of electronic structure theory.

Temperature equilibration in a fully ionized plasma: Electron-ion mass ratio effects

Brown, Preston, and Singleton (BPS) produced an analytic calculation for energy exchange processes for a weakly to moderately coupled plasma: the electron-ion temperature equilibration rate and the charged particle stopping power. These precise calculations are accurate to leading and next-to-leading order in the plasma coupling parameter and to all orders for two-body quantum scattering within the plasma. Classical molecular dynamics can provide another approach that can be rigorously implemented. It is therefore useful to compare the predictions from these two methods, particularly since the former is theoretically based and the latter numerically. An agreement would provide both confidence in our theoretical machinery and in the reliability of the computer simulations. The comparisons can be made cleanly in the purely classical regime, thereby avoiding the arbitrariness associated with constructing effective potentials to mock up quantum effects. We present here the classical limit of the general result for the temperature equilibration rate presented in BPS. In particular, we examine the validity of the m(electron)/m(ion)-->0 limit used in BPS to obtain a very simple analytic evaluation of the long-distance collective effects in the background plasma.

Dose point kernel for boron-11 decay and the cellular S values in boron neutron capture therapy

The study of the radiobiology of boron neutron capture therapy is based on the cellular level dosimetry of boron-10's thermal neutron capture reaction 10B(n,alpha)7Li, in which one 1.47 MeV helium-4 ion and one 0.84 MeV lithium-7 ion are spawned. Because of the chemical preference of boron-10 carrier molecules, the dose is heterogeneously distributed in cells. In the present work, the (scaled) dose point kernel of boron-11 decay, called 11B-DPK, was calculated by GEANT4 Monte Carlo simulation code. The DPK curve drops suddenly at the radius of 4.26 microm, the continuous slowing down approximation (CSDA) range of a lithium-7 ion. Then, after a slight ascending, the curve decreases to near zero when the radius goes beyond 8.20 microm, which is the CSDA range of a 1.47 MeV helium-4 ion. With the DPK data, S values for nuclei and cells with the boron-10 on the cell surface are calculated for different combinations of cell and nucleus sizes. The S value for a cell radius of 10 microm and a nucleus radius of 5 microm is slightly larger than the value published by Tung et al. [Appl. Radiat. Isot. 61, 739-743 (2004)]. This result is potentially more accurate than the published value since it includes the contribution of a lithium-7 ion as well as the alpha particle.

Stopping power in insulators and metals without charge exchange

The slowing-down process of pointlike charged particles in matter has been investigated by measuring the stopping power for antiprotons in materials of qualitatively very different nature. Whereas the velocity-proportional stopping power observed for metal-like targets such as aluminum over a wide energy range of 1-50 keV is in agreement with expectations, it is surprising that the same velocity dependence is seen for a large band-gap insulator such as LiF. The validity of these observations is supported by several measurements with protons and several checks of the target properties. The observations call for both a qualitative explanation and a quantitative theoretical model.

Antiproton stopping at low energies: confirmation of velocity-proportional stopping power

The stopping power for antiprotons in various solid targets has been measured in the low-energy range of 1-100 keV. In agreement with most models, in particular free-electron gas models, the stopping power is found to be proportional to the projectile velocity below the stopping-power maximum. Although a stopping power proportional to velocity has also been observed for protons, the interpretation of such measurements is difficult due to the presence of charge exchange processes. Hence, the present measurements constitute the first unambiguous support for a velocity-proportional stopping power due to target excitations by a pointlike projectile.