Hybrid Quantum Classical Molecular Dynamics Simulation of the Proton-Transfer Reaction of HO- with HBr in Aqueous Clusters

PyDFT-QMMM: A modular, extensible software framework for DFT-based QM/MM molecular dynamics

PyDFT-QMMM is a Python-based package for performing hybrid quantum mechanics/molecular mechanics (QM/MM) simulations at the density functional level of theory. The program is designed to treat short-range and long-range interactions through user-specified combinations of electrostatic and mechanical embedding procedures within periodic simulation domains, providing necessary interfaces to external quantum chemistry and molecular dynamics software. To enable direct embedding of long-range electrostatics in periodic systems, we have derived and implemented force terms for our previously described QM/MM/PME approach [Pederson and McDaniel, J. Chem. Phys. 156, 174105 (2022)]. Communication with external software packages Psi4 and OpenMM is facilitated through Python application programming interfaces (APIs). The core library contains basic utilities for running QM/MM molecular dynamics simulations, and plug-in entry-points are provided for users to implement custom energy/force calculation and integration routines, within an extensible architecture. The user interacts with PyDFT-QMMM primarily through its Python API, allowing for complex workflow development with Python scripting, for example, interfacing with PLUMED for free energy simulations. We provide benchmarks of forces and energy conservation for the QM/MM/PME and alternative QM/MM electrostatic embedding approaches. We further demonstrate a simple example use case for water solute in a water solvent system, for which radial distribution functions are computed from 100 ps QM/MM simulations; in this example, we highlight how the solvation structure is sensitive to different basis-set choices due to under- or over-polarization of the QM water molecule's electron density.

Chemical Reactivity and Spectroscopy Explored From QM/MM Molecular Dynamics Simulations Using the LIO Code

In this work we present the current advances in the development and the applications of LIO, a lab-made code designed for density functional theory calculations in graphical processing units (GPU), that can be coupled with different classical molecular dynamics engines. This code has been thoroughly optimized to perform efficient molecular dynamics simulations at the QM/MM DFT level, allowing for an exhaustive sampling of the configurational space. Selected examples are presented for the description of chemical reactivity in terms of free energy profiles, and also for the computation of optical properties, such as vibrational and electronic spectra in solvent and protein environments.

A quantum-mechanics molecular-mechanics scheme for extended systems

We introduce and discuss a hybrid quantum-mechanics molecular-mechanics (QM-MM) approach for Car-Parrinello DFT simulations with pseudopotentials and planewaves basis, designed for the treatment of periodic systems. In this implementation the MM atoms are considered as additional QM ions having fractional charges of either sign, which provides conceptual and computational simplicity by exploiting the machinery already existing in planewave codes to deal with electrostatics in periodic boundary conditions. With this strategy, both the QM and MM regions are contained in the same supercell, which determines the periodicity for the whole system. Thus, while this method is not meant to compete with non-periodic QM-MM schemes able to handle extremely large but finite MM regions, it is shown that for periodic systems of a few hundred atoms, our approach provides substantial savings in computational times by treating classically a fraction of the particles. The performance and accuracy of the method is assessed through the study of energetic, structural, and dynamical aspects of the water dimer and of the aqueous bulk phase. Finally, the QM-MM scheme is applied to the computation of the vibrational spectra of water layers adsorbed at the TiO2 anatase (1 0 1) solid-liquid interface. This investigation suggests that the inclusion of a second monolayer of H2O molecules is sufficient to induce on the first adsorbed layer, a vibrational dynamics similar to that taking place in the presence of an aqueous environment. The present QM-MM scheme appears as a very interesting tool to efficiently perform molecular dynamics simulations of complex condensed matter systems, from solutions to nanoconfined fluids to different kind of interfaces.

Influence of Thermal Fluctuations on Interfacial Electron Transfer in Functionalized TiO2 Semiconductors

The influence of thermal fluctuations on the dynamics of interfacial electron transfer in sensitized TiO2-anatase semiconductors is investigated by combining ab initio DFT molecular dynamics simulations and quantum dynamics propagation of transient electronic excitations. It is shown that thermal nuclear fluctuations speed up the underlying interfacial electron transfer dynamics by introducing nonadiabatic transitions between electron acceptor states, localized in the vicinity of the photoexcited adsorbate, and delocalized states extended throughout the semiconductor material, creating additional relaxation pathways for carrier diffusion. Furthermore, it is shown that room-temperature thermal fluctuations reduce the anisotropic character of charge diffusion along different directions in the anatase crystal and make similar the rates for electron injection from adsorbate states of different character. The reported results are particularly relevant to the understanding of temperature effects on surface charge separation mechanisms in molecular-based photo-optic devices.

Solvent Effects on Peroxynitrite Structure and Properties from QM/MM Simulations

We have investigated the structure and the vibrational spectrum of peroxynitrite anion in aqueous solution by means of combined quantum-classical (QM/MM) molecular dynamics simulations. In our QM/MM scheme, the reactant was modeled using density functional theory with a Gaussian basis set and the solvent was described using the mean-field TIP4P and the polarizable TIP4P-FQ force fields. The choice of basis sets, functionals and force field parameters has been validated by performing calculations on isolated peroxynitrite and on small peroxynitrite-water complexes. Poor values for isolated peroxynitrite structural properties and vibrational frequencies are found for most ab initio methods, particularly regarding the ON-OO(-) bond distance and the harmonic stretching frequency, probably due to the singlet-triplet instability found in the HF wave function. On the other hand, DFT methods yield results in better agreement with high level CCSD(T) ab initio calculations. We have computed the vibrational spectrum for aqueous peroxynitrite by calculating the Fourier transform of the velocity autocorrelation function extracted from the QM-MM molecular dynamics simulations. Our computational scheme, which allows for the inclusion of both anharmonicity and solvent effects, is able to clarify previous discrepancies regarding the experimental spectra assignments and to shed light on the subtle interplay between solvation and peroxynitrite structure and properties.

Solvation and Structure of LiAlH4 in Ethereal Solvents

The nature of the solute species present in ethereal solutions of LiAlH(4) is of crucial importance for understanding the mechanisms for the reduction of ketones and other functional groups by LiAlH(4). We have employed a combination of theoretical and experimental techniques to investigate the structure of LiAlH(4) in ethereal solutions. Using complexation agents, we measured the IR spectra of LiAlH(4) and AlH(4)(-) in tetrahydrofuran (THF). Hybrid quantum-classical (QM-MM) simulations have also been carried out to compute the IR spectra of associated and dissociated LiAlH(4) species and the free-energy profile for the dissociation process in solution. Our experimental estimate of the dissociation constant in THF is 0.021 +/- 0.002, while the predicted computational value corresponding to a model dimethyl ether solvent is 0.001. The free-energy profile shows only one minimum corresponding to a contact ion pair at a Li-Al separation distance of 3.0 A. These results are consistent with the fact that LiAlH(4) is essentially associated in ethereal solutions forming contact ion pairs.

Molecular aspects of halide ion hydration: the cluster approach

This review provides a historical context for our understanding of the hydration shell surrounding halide ions and illustrates how the cluster systems can be used, in combination with theory, to elucidate the behavior of water molecules in direct contact with the anion. We discuss how vibrational predissociation spectroscopy, carried out with weakly bound argon atoms, has been employed to deduce the morphology of the small water networks attached to anions in the primary steps of hydration. We emphasize the importance of charge-transfer in the binary interaction, and discuss how this process affects the structures of the larger networks. Finally, we survey how the negatively charged water clusters (H2O)n(-) are providing a molecular-level perspective on how diffuse excess electrons interact with the water networks.

Nitric oxide binding to ferric cytochrome P450: a computational study

The interaction between nitric oxide (NO) and the active site of ferric cytochrome P450 was studied by means of density functional theory (DFT), at the generalized gradient approximation level, and of the SAM1 semiempirical method. The electrostatic effects of the protein environment were included in our DFT scheme by using a hybrid quantum classical approach. The active-site model consisted of an iron(III) porphyrin, the adjacent cysteine residue, and one coordinated water molecule. For this system, spin populations and relative energies for selected spin states were computed. Interestingly, the unpaired electron density, the HOMO, and the LUMO were found to be highly localized on the iron and in an appreciable extent on the sulfur coordinated to the metal. This provides central information about the reactivity of nitric oxide with the active site. Since the substitution of a molecule of H2O by NO has been proposed as being responsible for the inhibition of the cytochrome in the presence of nitric oxide, we have analyzed the thermodynamic feasibility of the ligand exchange process. The structure of the nitrosylated active site was partially optimized using SAM1. A low-spin ground state was obtained for the nitrosyl complex, with a linear Fe-N-O angle. The trends found in Fe-N-O angles and Fe-N lengths of the higher energy spin states provided a notable insight into the electronic configuration of the complex within the framework of the Enemark and Feltham formalism. In relation to the protein environment, it was assessed that the electrostatic field has significant effects on several computed properties. However, in both vacuum and protein environments, the ligand exchange reaction turned out to be exergonic and the relative orders of spin states of the relevant species were the same.