Source author record

Daniel Muñoz-Santiburcio

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

physics.chem-ph cond-mat.mtrl-sci cond-mat.soft physics.comp-ph physics.med-ph

Catalog footprint

Research overview

3works

5topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2022arXiv

Accurate diffusion coefficients of the excess proton and hydroxide in water via extensive ab initio simulations with different schemes

[ This article may be downloaded for personal use only. Any other use requires prior permission of the author and AIP Publishing. This article appeared in The Journal of Chemical Physics, 157, 024504 (2022) and may be found at: https://doi.org/10.1063/5.0093958 ] Despite its simple molecular formula, obtaining an accurate in silico description of water is far from straightforward. Many of its very peculiar properties are quite elusive, and in particular obtaining good estimations of the diffusion coefficients of the solvated proton and hydroxide at a reasonable computational cost has been an unsolved challenge until now. Here I present extensive results of several unusually long ab initio MD simulations employing different combinations of the Born--Oppenheimer and 2 $^\mathrm{nd}$ -generation Car--Parrinello MD propagation methods with different ensembles (NVE and NVT) and thermostats, which show that these methods together with the RPBE-D3 functional provide a very accurate estimation of the diffusion coefficients of the solvated H $_3$ O $^+$ and OH $^-$ ions, together with an extremely accurate description of several properties of neutral water (such as the structure of the liquid and its diffusion and shear viscosity coefficients). In addition, I show that the estimations of $D_\mathrm{H_3O^+}$ and $D_\mathrm{OH^-}$ depend dramatically on the simulation length, being necessary to reach timescales in the order of hundreds of picoseconds to obtain reliable results.

preprint2022arXiv

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

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 $\mathrm{2H_2+O_2}$ $\mathrm{\rightarrow 2H_2O}$ 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.

preprint2020arXiv

Efficient ab initio calculation of electronic stopping in disordered systems via geometry pre-sampling: application to liquid water

Knowledge of the electronic stopping curve for swift ions, $S_e(v)$ , particularly around the Bragg peak, is important for understanding radiation damage. Experimentally, however, the determination of such feature for light ions is very challenging, especially in disordered systems such as liquid water and biological tissue. Recent developments in real-time time-dependent density functional theory (rt-TDDFT) have enabled the calculation of $S_e(v)$ along nm-sized trajectories. However, it is still a challenge to obtain a meaningful statistically averaged $S_e(v)$ that can be compared to observations. In this work, taking advantage of the correlation between the local electronic structure probed by the projectile and the distance from the projectile to the atoms in the target, we devise a trajectory pre-sampling scheme to select, geometrically, a small set of short trajectories to accelerate the convergence of the averaged $S_e(v)$ computed via rt-TDDFT. For protons in liquid water, we first calculate the reference probability distribution function (PDF) for the distance from the proton to the closest oxygen atom, $ϕ_R(r_{p{\rightarrow}O})$ , for a trajectory of a length similar to those sampled experimentally. Then, short trajectories are sequentially selected so that the accumulated PDF reproduces $ϕ_R(r_{p{\rightarrow}O})$ to increasingly high accuracy. Using these pre-sampled trajectories, we demonstrate that the averaged $S_e(v_p)$ converges in the whole velocity range with less than eight trajectories, while other averaging methods using randomly and uniformly distributed trajectories require approximately ten times the computational effort. This allows us to compare the $S_e(v_p)$ curve to experimental data, and assess widely used empirical tables based on Bragg's rule.

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.