Liquid Water through Density-Functional Molecular Dynamics: Plane-Wave vs Atomic-Orbital Basis Sets

Wavefunction-Based Electrostatic-Embedding QM/MM Using CFOUR through MiMiC

We present an interface of the wavefunction-based quantum chemical software CFOUR to the multiscale modeling framework MiMiC. Electrostatic embedding of the quantum mechanical (QM) part is achieved by analytic evaluation of one-electron integrals in CFOUR, while the rest of the QM/molecular mechanical (MM) operations are treated according to the previous MiMiC-based QM/MM implementation. Long-range electrostatic interactions are treated by a multipole expansion of the potential from the QM electron density to reduce the computational cost without loss of accuracy. Testing on model water/water systems, we verified that the CFOUR interface to MiMiC is robust, guaranteeing fast convergence of the self-consistent field cycles and optimal conservation of the energy during the integration of the equations of motion. Finally, we verified that the CFOUR interface to MiMiC is compatible with the use of a QM/QM multiple time-step algorithm, which effectively reduces the cost of ab initio MD (AIMD) or QM/MM-MD simulations using higher level wavefunction-based approaches compared to cheaper density functional theory-based ones. The new wavefunction-based AIMD and QM/MM-MD implementations were tested and validated for a large number of wavefunction approaches, including Hartree-Fock and post-Hartree-Fock methods like Møller-Plesset, coupled-cluster, and complete active space self-consistent field.

Combined Computational and Experimental Investigation on the Nature of Hydrated Iodoplumbate Complexes: Insights into the Dual Role of Water in Perovskite Precursor Solutions

Water is generally considered an enemy of metal halide perovskites, being responsible for their rapid degradation and, consequently, undermining the long-term stability of perovskite-based solar cells. However, beneficial effects of liquid water have been surprisingly observed, and synthetic routes including water treatments have shown to improve the quality of perovskite films. This suggests that the interactions of water with perovskites and their precursors are far from being completely understood, as water appears to play a puzzling dual role in perovskite precursor solutions. In this context, studying the basic interactions between perovskite precursors in the aqueous environment can provide a deeper comprehension of this conundrum. In this context, it is fundamental to understand how water impacts the chemistry of iodoplumbate perovskite precursor species, PbI_x^2–x. Here, we investigate the chemistry of these complexes using a combined experimental and theoretical strategy to unveil their peculiar structural and optical properties and eventually to assign the species present in the solution. Our study indicates that iodide-rich iodoplumbates, which are generally key to the formation of lead halide perovskites, are not easily formed in aqueous solutions because of the competition between iodide and solvent molecules in coordinating Pb²⁺ ions, explaining the difficulty of depositing lead iodide perovskites from aqueous solutions. We postulate that the beneficial effect of water when used as an additive is then motivated by its behavior being similar to high coordinative polar aprotic solvents usually employed as additives in one-step perovskite depositions.

Quantifying the hydration structure of sodium and potassium ions: taking additional steps on Jacob's Ladder

The ability to reproduce the experimental structure of water around the sodium and potassium ions is a key test of the quality of interaction potentials due to the central importance of these ions in a wide range of important phenomena.

Aqueous solvation of the chloride ion revisited with density functional theory: impact of correlation and exchange approximations

Aqueous chloride is simulated using PBE-D3, PBE0-D3, and SCAN to investigate the impact of exchange and correlation approximations; we find the exact exchange fraction strongly impacts the energetics and polarizability of solvated chloride.

pH-Dependent Surface Chemistry from First Principles: Application to the BiVO4(010)-Water Interface

We present a theoretical formulation for studying the pH-dependent interfacial coverage of semiconductor-water interfaces through ab initio electronic structure calculations, molecular dynamics simulations, and the thermodynamic integration method. This general methodology allows one to calculate the acidity of the individual adsorption sites on the surface and consequently the pH at the point of zero charge, pH_PZC, and the preferential adsorption mode of water molecules, either molecular or dissociative, at the semiconductor-water interface. The proposed method is applied to study the BiVO₄(010)-water interface, yields a pH_PZC in excellent agreement with the experimental characterization. Furthermore, from the calculated p K_a values of the individual adsorption sites, we construct an ab initio concentration diagram of all adsorbed species at the interface as a function of the pH of the aqueous solution. The diagram clearly illustrates the pH-dependent coverage of the surface and indicates that protons are found to be significantly adsorbed (∼1% of available sites) only in highly acidic conditions. The surface is found to be mostly covered by molecularly adsorbed water molecules in a wide interval of pH values ranging from 2 to 8. Hydroxyl ions are identified as the dominant adsorbed species at pH larger than 8.2.

Born-Oppenheimer Molecular Dynamics Simulations of a Bromate Ion in Water Reveal Its Dual Kosmotropic and Chaotropic Behavior

The solvation structure and dynamics of a bromate (BrO₃^-) ion in water are studied by means of Born-Oppenheimer molecular dynamics simulations at two different temperatures using the Becke-Lee-Yang-Parr functional with Grimme D3 dispersion corrections. The bromate ion possesses a pyramidal structure, and it has two types of solvation sites, namely, the bromine and oxygen atoms. We have looked at different radial and orientational distributions of water molecules around the bromate ion and also investigated their hydrogen bonding properties. The solvation structure of the bromate ion is also compared with that of the iodate (IO₃^-) ion, which is structurally rather similar to the bromate ion and was found to have some unusual solvation properties in water. It is found that the bromate ion follows a similar trend as that followed by the iodate ion as far as the solvation structure is concerned. However, the effect of the former on surrounding water is found to be much weaker than that of the latter. On the dynamical side, we have looked at diffusion, residence dynamics, and also the orientational and hydrogen bond relaxation of water molecules around the BrO₃^- ion and compared them with those of the bulk. Dynamical results are presented for both H₂O and D₂O around the BrO₃^- ion. Interpretation of the dynamical results in terms of structure-making (kosmotropic)/-breaking (chaotropic) properties of the BrO₃^- ion reveals that the bromine atom of this ion acts as a water structure breaker, whereas the three oxygens act as water structure makers. Thus, in spite of being a single ion, the bromate ion has dual characteristics and the experimentally observed kosmotropic ability of this ion is actually a trade-off between a chaotropic site (the bromine atom) and three kosmotropic sites (three oxygen atoms) that are present in the ion.

Ab initio theory and modeling of water

Water is of the utmost importance for life and technology. However, a genuinely predictive ab initio model of water has eluded scientists. We demonstrate that a fully ab initio approach, relying on the strongly constrained and appropriately normed (SCAN) density functional, provides such a description of water. SCAN accurately describes the balance among covalent bonds, hydrogen bonds, and van der Waals interactions that dictates the structure and dynamics of liquid water. Notably, SCAN captures the density difference between water and ice Ih at ambient conditions, as well as many important structural, electronic, and dynamic properties of liquid water. These successful predictions of the versatile SCAN functional open the gates to study complex processes in aqueous phase chemistry and the interactions of water with other materials in an efficient, accurate, and predictive, ab initio manner.

Ab initio molecular dynamics simulations of liquid water using high quality meta-GGA functionals

We have used ab initio molecular dynamics (AIMD) to characterize water properties using two meta-generalized gradient approximation (meta-GGA) functionals, M06-L-D3 and B97M-rV, and compared their performance against a standard GGA corrected for dispersion, revPBE-D3, at ambient conditions (298 K, and 1 g cm-3 or 1 atm). Simulations of the equilibrium density, radial distribution functions, self-diffusivity, the infrared spectrum, liquid dipole moments, and characterizations of the hydrogen bond network show that all three functionals have overcome the problem of the early AIMD simulations that erroneously found ambient water to be highly structured, but they differ substantially among themselves in agreement with experiment on this range of water properties. We show directly using water cluster data up through the pentamer that revPBE-D3 benefits from a cancellation of its intrinsic functional error by running classical trajectories, whereas the meta-GGA functionals are demonstrably more accurate and would require the simulation of nuclear quantum effects to realize better agreement with all cluster and condensed phase properties.