Fully Quantum Description of the Zundel Ion: Combining Variational Quantum Monte Carlo with Path Integral Langevin Dynamics

Identifying the acidic or basic behavior of surface water: a QM/MM-MD study

Controversies on the water surface were theoretically addressed with the help of large scale quantum mechanical molecular dynamics (QMMD) simulations on water surface model systems with and without excess hydroniums and hydroxides. It was revealed that the thermodynamic surface structures of these ions strongly depend on their location and dipole orientation. Fast hydronium diffusion by proton transfer establishes a wider kinetic depth distribution (∼6 Å) than that predicted by its thermodynamic affinity for the water surface, while slow hydroxide is shallowly trapped below the outermost molecular layer (3-4 Å). In addition, the anisotropic orientation of surface water dipole can generate a substantial magnitude of surface potential, which extends to a depth of a few molecular layers. With these distinctively different surface properties of two ions and water molecules, the seemingly contradictory observations of acidic and negatively charged water surfaces may be successfully explained. That is, the negative surface charge of neutral water mostly stems from intrinsic water properties such as water dipole orientation and electron density spillage at the surface, rather than surface OH- ions. The enhanced acidity of the water surface can be attributed in large part to the kinetic depth profile of ion density in addition to static thermodynamic origin. Furthermore, the different depth profiles of the two ions may differently affect the surface-sensitive spectroscopic observations.

Thermal dependence of the hydrated proton and optimal proton transfer in the protonated water hexamer

Water is a key ingredient for life and plays a central role as solvent in many biochemical reactions. However, the intrinsically quantum nature of the hydrogen nucleus, revealing itself in a large variety of physical manifestations, including proton transfer, gives rise to unexpected phenomena whose description is still elusive. Here we study, by a combination of state-of-the-art quantum Monte Carlo methods and path-integral molecular dynamics, the structure and hydrogen-bond dynamics of the protonated water hexamer, the fundamental unit for the hydrated proton. We report a remarkably low thermal expansion of the hydrogen bond from zero temperature up to 300 K, owing to the presence of short-Zundel configurations, characterised by proton delocalisation and favoured by the synergy of nuclear quantum effects and thermal activation. The hydrogen bond strength progressively weakens above 300 K, when localised Eigen-like configurations become relevant. Our analysis, supported by the instanton statistics of shuttling protons, reveals that the near-room-temperature range from 250 K to 300 K is optimal for proton transfer in the protonated water hexamer.

Probing anharmonic phonons by quantum correlators: A path integral approach

We devise an efficient scheme to determine vibrational properties from Path Integral Molecular Dynamics (PIMD) simulations. The method is based on zero-time Kubo-transformed correlation functions and captures the anharmonicity of the potential due to both temperature and quantum effects. Using analytical derivations and numerical calculations on toy-model potentials, we show that two different estimators built upon PIMD correlation functions fully characterize the phonon spectra and the anharmonicity strength. The first estimator is associated with the force-force quantum correlators and, in the weak anharmonic regime, yields reliable zero-point motion frequencies and thermodynamic properties of the quantum system. The second one is instead connected to displacement-displacement correlators and accurately probes the lowest-energy phonon excitations, regardless of the anharmonicity strength of the system. We also prove that the use of generalized eigenvalue equations, in place of the standard normal mode equations, leads to a significant speed-up in the PIMD phonon calculations, both in terms of faster convergence rate and smaller time step bias. Within this framework, using ab initio PIMD simulations, we compute phonon dispersions of diamond and of the high-pressure I4₁/amd phase of atomic hydrogen. We find that in the latter case, the anharmonicity is stronger than previously estimated and yields a sizeable red-shift in the vibrational spectrum of atomic hydrogen.

Toward quantum Monte Carlo forces on heavier ions: Scaling properties

Quantum Monte Carlo (QMC) forces have been studied extensively in recent decades because of their importance with spectroscopic observables and geometry optimization. Here, we benchmark the accuracy and computational cost of QMC forces. The zero-variance zero-bias (ZVZB) force estimator is used in standard variational and diffusion Monte Carlo simulations with mean-field based trial wavefunctions and atomic pseudopotentials. Statistical force uncertainties are obtained with a recently developed regression technique for heavy tailed QMC data [P. Lopez Rios and G. J. Conduit, Phys. Rev. E 99, 063312 (2019)]. By considering selected atoms and dimers with elements ranging from H to Zn (1 ≤ Z_eff ≤ 20), we assess the accuracy and the computational cost of ZVZB forces as the effective pseudopotential valence charge, Z_eff, increases. We find that the costs of QMC energies and forces approximately follow simple power laws in Z_eff. The force uncertainty grows more rapidly, leading to a best case cost scaling relationship of approximately Z_eff ^6.5(3) for diffusion Monte Carlo. We find that the accessible system size at fixed computational cost scales as Z_eff ^-2, insensitive to model assumptions or the use of the "space warp" variance-reduction technique. Our results predict the practical cost of obtaining forces for a range of materials, such as transition metal oxides where QMC forces have yet to be applied, and underscore the importance of further developing force variance-reduction techniques, particularly for atoms with high Z_eff.

Molecular Characterization of the Surface Excess Charge Layer in Droplets

The surface excess charge layer (SECL) in droplets has often been associated with distinct chemistry. We examine the effect of the nature of ions in the composition and structure of SECL by using molecular dynamics. We find that in the presence of simple ions the thickness of SECL is invariant not only with respect to droplet size but also with respect to the nature of the ions. In the presence of simple ions, this layer has a thickness of ∼1.5-1.7 nm but in the presence of macroions it may extend to ∼2.0 nm. The proportion of ions contained in SECL depends on the nature of the ions and the droplet size. For the same droplet size, I^- and model H₃O⁺ ions show considerably higher concentration than Na⁺ and Cl^- ions. We identify the maximum ion concentration region, which, in nanodrops, may partially overlap with SECL. As the relative shape fluctuations decrease when microdrop size is approached, the overlap between SECL and maximum ion concentration region increases. We suggest the extension of the bilayer droplet structure assumed in the equilibrium partitioning model of Enke to include the maximum ion concentration region that may not coincide with SECL in nanodrops. We compute the ion concentrations in SECL, which are those that should enter the kinetic equation in the ion-evaporation mechanism, instead of the overall drop ion concentration that has been used.

Nuclear Quantum Effects Largely Influence Molecular Dissociation and Proton Transfer in Liquid Water under an Electric Field

Proton transfer in liquid water controls acid-base chemistry, crucial enzyme reactions, and the functioning of fuel cells. Externally applied static electric fields in water are capable of dissociating molecules and transferring protons across the H-bond network. However, the impact of nuclear quantum effects (NQEs) on these fundamental field-induced phenomena has not yet been reported. By comparing state-of-the-art ab initio molecular dynamics (AIMD) and path integral AIMD simulations of water under electric fields, I show that quantum delocalization of the proton lowers the molecular ionization threshold to approximately one-third. Moreover, also the water behavior as a protonic semiconductor is considerably modified by the inclusion of NQEs. In fact, when the quantum nature of the nuclei is taken into account, the proton conductivity is ∼50% larger. This work proves that NQEs sizably affect the protolysis phenomenon and proton transfer in room-temperature liquid water.

Reduced rovibrational coupling Cartesian dynamics for semiclassical calculations: Application to the spectrum of the Zundel cation

We study the vibrational spectrum of the protonated water dimer, by means of a divide-and-conquer semiclassical initial value representation of the quantum propagator, as a first step in the study of larger protonated water clusters. We use the potential energy surface from the work of Huang et al. [J. Chem. Phys. 122, 044308 (2005)]. To tackle such an anharmonic and floppy molecule, we employ fully Cartesian dynamics and carefully reduce the coupling to global rotations in the definition of normal modes. We apply the time-averaging filter and obtain clean power spectra relative to suitable reference states that highlight the spectral peaks corresponding to the fundamental excitations of the system. Our trajectory-based approach allows for the physical interpretation of the very challenging proton transfer modes. We find that it is important, for such a floppy molecule, to selectively avoid initially exciting lower energy modes, in order to obtain cleaner spectra. The estimated vibrational energies display a mean absolute error (MAE) of ∼29 cm^-1 with respect to available multiconfiguration time-dependent Hartree calculations and MAE ∼ 14 cm^-1 when compared to the optically active experimental excitations of the Ne-tagged Zundel cation. The reasonable scaling in the number of trajectories for Monte Carlo convergence is promising for applications to higher dimensional protonated cluster systems.

Sampling the thermal Wigner density via a generalized Langevin dynamics

The Wigner thermal density is a function of considerable interest in the area of approximate (linearized or semiclassical) quantum dynamics where it is employed to generate initial conditions for the propagation of appropriate sets of classical trajectories. In this paper, we propose an original approach to compute the Wigner density based on a generalized Langevin equation. The stochastic dynamics is nontrivial in that it contains a coordinate-dependent friction coefficient and a generalized force that couples momenta and coordinates. These quantities are, in general, not known analytically and have to be estimated via auxiliary calculations. The performance of the new sampling scheme is tested on standard model systems with highly nonclassical features such as relevant zero point energy effects, correlation between momenta and coordinates, and negative parts of the Wigner density. In its current brute force implementation, the algorithm, whose convergence can be systematically checked, is accurate and has only limited overhead compared to schemes with similar characteristics. We briefly discuss potential ways to further improve its numerical efficiency.

Streams, cascades, and pools: various water cluster motifs in structurally similar Ni(ii) complexes

Hydrogen bonding (HB) interactions are well known to impact the properties of water in the bulk and within hydrated materials.