Electron nuclear dynamics with plane wave basis sets: complete theory and formalism.

Electron nuclear dynamics (END) is an ab initio quantum dynamics method that adopts a time-dependent, variational, direct, and non-adiabatic approach. The simplest-level (SL) END (SLEND) version employs a classical mechanics description for nuclei and a Thouless single-determinantal wave function for electrons. A higher-level END version, END/Kohn-Sham density functional theory, improves the electron correlation description of SLEND. While both versions can simulate various types of chemical reactions, they have difficulties to simulate scattering/capture of electrons to/from the continuum due to their reliance on localized Slater-type basis functions. To properly describe those processes, we formulate END with plane waves (PWs, END/PW), basis functions able to represent both bound and unbound electrons. As extra benefits, PWs also afford fast algorithms to simulate periodic systems, parametric independence from nuclear positions and momenta, and elimination of basis set linear dependencies and orthogonalization procedures. We obtain the END/PW formalism by extending the Thouless wave function and associated electron density to periodic systems, expressing the energy terms as functionals of the latter entities, and deriving the energy gradients with respect to nuclear and electronic variables. END/ PW has a great potential to simulate electron processes in both periodic (crystal) and aperiodic (molecular) systems (the latter in a supercell approach). Following previous END studies, END/PW will be applied to electron scattering processes in proton cancer therapy reactions.

Symmetry-breaking effects on time-dependent dynamics: correct differential cross sections and other properties in H+ + C2H4 at ELab = 30 eV.

We present a computational procedure that introduces low degrees of symmetry breaking into a restricted Hartree-Fock (RHF) state in order to induce higher symmetry breaking during the state's subsequent dynamics. The symmetries herein considered are those of electronic HF states as classified by Fukutome; those symmetries affect bond dissociations and internal rotations among other phenomena. Therefore, this investigation extends a part of Fukutome's time-independent analysis of symmetry breaking to the time-dependent (dynamical) regime. The procedure is formulated in the framework of the simplest-level electron nuclear dynamics, a time-dependent, variational, on-the-fly and non-adiabatic method that employs classical dynamics for the nuclei and a Thouless single-determinantal state for the electrons. We test this procedure on the H+ + C2H4 reaction at 30 eV due to its conspicuous display of symmetry-breaking effects; this reaction is relevant in astrophysics and proton cancer therapy. Fukutome's axial spin density wave (ASDW) HF state is used to represent the symmetry-broken initial states. Through a Thouless parameter, small degrees of symmetry breaking are introduced into the initial ASDW states in a controlled manner. After starting the dynamics from those states, higher degrees of symmetry breaking emerge or not as determined by the direct-dynamics equations without external interventions. Simulations starting from symmetry-conforming states preserve symmetry features during dynamics, whereas simulations starting from symmetry-broken states display an upsurge of symmetry breaking when the reactants collide. Initial symmetry breaking increases the total integral cross sections of collision-induced fragmentations and of target-to-proton 1-electron-transfer reactions and decreases the scattering angle function and primary rainbow angle of the outgoing projectile. Remarkably, symmetry-breaking simulations reproduce the correct relative order and values of the experimental 0- and 1-electron-transfer differential cross sections, whereas symmetry-conforming simulations predict incorrect order and values. Our calculated scattering angle functions and differential cross sections also exhibit expected primary and secondary rainbow angle features that experiments fail to detect. A detailed discussion on the description of symmetry-breaking processes with the ASDW and Thouless states is included to provide a rigorous theoretical basis for this investigation.

Electron Nuclear Dynamics Simulations of Proton Cancer Therapy Reactions: Water Radiolysis and Proton- and Electron-Induced DNA Damage in Computational Prototypes.

Proton cancer therapy (PCT) utilizes high-energy proton projectiles to obliterate cancerous tumors with low damage to healthy tissues and without the side effects of X-ray therapy. The healing action of the protons results from their damage on cancerous cell DNA. Despite established clinical use, the chemical mechanisms of PCT reactions at the molecular level remain elusive. This situation prevents a rational design of PCT that can maximize its therapeutic power and minimize its side effects. The incomplete characterization of PCT reactions is partially due to the health risks associated with experimental/clinical techniques applied to human subjects. To overcome this situation, we are conducting time-dependent and non-adiabatic computer simulations of PCT reactions with the electron nuclear dynamics (END) method. Herein, we present a review of our previous and new END research on three fundamental types of PCT reactions: water radiolysis reactions, proton-induced DNA damage and electron-induced DNA damage. These studies are performed on the computational prototypes: proton + H2O clusters, proton + DNA/RNA bases and + cytosine nucleotide, and electron + cytosine nucleotide + H2O. These simulations provide chemical mechanisms and dynamical properties of the selected PCT reactions in comparison with available experimental and alternative computational results.

Exploring water radiolysis in proton cancer therapy: Time-dependent, non-adiabatic simulations of H+ + (H2O)1-6.

To elucidate microscopic details of proton cancer therapy (PCT), we apply the simplest-level electron nuclear dynamics (SLEND) method to H+ + (H2O)1-6 at ELab = 100 keV. These systems are computationally tractable prototypes to simulate water radiolysis reactions-i.e. the PCT processes that generate the DNA-damaging species against cancerous cells. To capture incipient bulk-water effects, ten (H2O)1-6 isomers are considered, ranging from quasi-planar/multiplanar (H2O)1-6 to "smallest-drop" prism and cage (H2O)6 structures. SLEND is a time-dependent, variational, non-adiabatic and direct method that adopts a nuclear classical-mechanics description and an electronic single-determinantal wavefunction in the Thouless representation. Short-time SLEND/6-31G* (n = 1-6) and /6-31G** (n = 1-5) simulations render cluster-to-projectile 1-electron-transfer (1-ET) total integral cross sections (ICSs) and 1-ET probabilities. In absolute quantitative terms, SLEND/6-31G* 1-ET ICS compares satisfactorily with alternative experimental and theoretical results only available for n = 1 and exhibits almost the same accuracy of the best alternative theoretical result. SLEND/6-31G** overestimates 1-ET ICS for n = 1, but a comparable overestimation is also observed with another theoretical method. An investigation on H+ + H indicates that electron direct ionization (DI) becomes significant with the large virtual-space quasi-continuum in large basis sets; thus, SLEND/6-31G** 1-ET ICS is overestimated by DI contributions. The solution to this problem is discussed. In relative quantitative terms, both SLEND/6-31* and /6-31G** 1-ET ICSs precisely fit into physically justified scaling formulae as a function of the cluster size; this indicates SLEND's suitability for predicting properties of water clusters with varying size. Long-time SLEND/6-31G* (n = 1-4) simulations predict the formation of the DNA-damaging radicals H, OH, O and H3O. While "smallest-drop" isomers are included, no early manifestations of bulk water PCT properties are observed and simulations with larger water clusters will be needed to capture those effects. This study is the largest SLEND investigation on water radiolysis to date.