Bragg's Additivity Rule and Core and Bond model studied by real-time TDDFT electronic stopping simulations: the case of water vapor

Abstract

The electronic stopping power ($S_e$) of water vapor (H${}_2$O), hydrogen (H${}_2$) and oxygen (O${}_2$) gases for protons in a broad range of energies, centered in the Bragg peak, was calculated using real-time time-dependent density functional theory (rt-TDDFT) simulations with Gaussian basis sets. This was done for a kinetic energy of incident protons ($E_k$) ranging from 1.56 keV/amu to 1.6 MeV/amu. $S_e$ was calculated as the average over geometrically pre-sampled short ion trajectories. The average $S_e(E_k)$ values were found to rapidly converge with 25-30 pre-sampled, 2 nm-long ion trajectories. The rt-TDDFT $S_e(E_k)$ curves were compared to experimental and SRIM data, and used to validate the Bragg's Additivity Rule (BAR). Discrepancies were analyzed in terms of basis set effects and omitted nuclear stopping at low energies. At variance with SRIM, we found that BAR is applicable to our rt-TDDFT simulations of $2{\mathrm{H}}_2+{\mathrm{O}}_2$ $\to 2{\mathrm{H}}_2\mathrm{O}$ without scaling for $E_k>40$ keV/amu. The hydrogen and oxygen Core and Bond (CAB) contributions to electronic stopping were calculated and found to be slightly smaller than SRIM values as a result of a red-shift in our rt-TDDFT $S_e(E_k)$ curves and a re-distribution of weights due to some bond contributions being neglected in SRIM.

15 pages 5 figures, 1 table

References in corpus (4)

• Electronic stopping power in insulators from first principles
• Efficient ab initio calculation of electronic stopping in disordered systems via geometry pre-sampling: application to liquid water
• Water Radiolysis by Low-Energy Carbon Projectiles from First-Principles Molecular Dynamics
• INQ, a modern GPU-accelerated computational framework for (time-dependent) density functional theory

Cited by in corpus (1)

• Stopping cross-section for protons across different phases of water