Vibrational Spectra of Highly Anharmonic Water Clusters: Molecular Dynamics and Harmonic Analysis Revisited with Constrained Nuclear-Electronic Orbital Methods

Modeling electronic absorption spectra with nuclear quantum effects in constrained nuclear-electronic orbital framework

Electronic absorption spectra serve as versatile and powerful tools in experiments. Accurate theoretical simulation of electronic absorption spectra is challenging because multiple factors such as environmental effects and nuclear quantum effects contribute to spectrum lineshapes. This work proposes a protocol to model electronic absorption spectra in the constrained nuclear-electronic orbital framework. Solvent effects, temperature effects, and particularly nuclear quantum effects can be taken into consideration in this unified framework. This protocol is applied to investigate the electronic absorption spectrum of the pyridine molecule in water. Nuclear quantum effects are found to induce a broadening and red shift of the absorption spectrum of pyridine.

Constrained Nuclear-Electronic Orbital Transition State Theory Using Energy Surfaces with Nuclear Quantum Effects

Hydrogen-atom transfer is crucial in a myriad of chemical and biological processes, yet the accurate and efficient description of hydrogen-atom transfer reactions and kinetic isotope effects remains challenging due to significant quantum effects on hydrogenic motion, especially tunneling and zero-point energy. In this paper, we combine transition state theory (TST) with the recently developed constrained nuclear-electronic orbital (CNEO) theory to propose a new transition state theory denoted CNEO-TST. We use CNEO-TST with CNEO density functional theory (CNEO-DFT) to predict reaction rate constants for two prototypical gas-phase hydrogen-atom transfer reactions and their deuterated isotopologic reactions. CNEO-TST is similar to conventional TST except that it employs constrained minimized energy surfaces to include zero-point energy and shallow tunneling effects in the effective potential. We find that the new theory predicts reaction rates quite accurately at room temperature. The effective potential surface must be generated by CNEO theory rather than by ordinary electronic structure theory, but because of the favorable computational scaling of CNEO-DFT, the cost is economical even for large systems. Our results show that dynamics calculations with this approach achieve accuracy comparable to variational TST with a semiclassical multidimensional tunneling transmission coefficient at and above room temperature. Therefore, CNEO-TST can be a useful tool for rate prediction, even for reactions involving highly quantal motion, such as many chemical and biochemical reactions involving transfers of hydrogen atoms, protons, or hydride ions.

Accuracy of Reaction Coordinate Based Rate Theories for Modelling Chemical Reactions: Insights From the Thermal Isomerization in Retinal

Modern potential energy surfaces have shifted attention to molecular simulations of chemical reactions. While various methods can estimate rate constants for conformational transitions in molecular dynamics simulations, their applicability to studying chemical reactions remains uncertain due to the high and sharp energy barriers and complex reaction coordinates involved. This study focuses on the thermal cis-trans isomerization in retinal, employing molecular simulations and comparing rate constant estimates based on one-dimensional rate theories with those based on sampling transitions and grid-based models for low-dimensional collective variable spaces. Even though each individual method to estimate the rate passes its quality tests, the rate constant estimates exhibit considerable disparities. Rate constant estimates based on one-dimensional reaction coordinates prove challenging to converge, even if the reaction coordinate is optimized. However, consistent estimates of the rate constant are achieved by sampling transitions and by multi-dimensional grid-based models.

Assessment of electron-proton correlation functionals for vibrational spectra of shared-proton systems by constrained nuclear-electronic orbital density functional theory

Proton transfer plays a crucial role in various chemical and biological processes. A major theoretical challenge in simulating proton transfer arises from the quantum nature of the proton. The constrained nuclear-electronic orbital (CNEO) framework was recently developed to efficiently and accurately account for nuclear quantum effects, particularly quantum nuclear delocalization effects, in quantum chemistry calculations and molecular dynamics simulations. In this paper, we systematically investigate challenging proton transfer modes in a series of shared-proton systems using CNEO density functional theory (CNEO-DFT), focusing on evaluating existing electron-proton correlation functionals. Our results show that CNEO-DFT accurately describes proton transfer vibrational modes and significantly outperforms conventional DFT. The inclusion of the epc17-2 electron-proton correlation functional in CNEO-DFT produces similar performance to that without electron-proton correlations, while the epc17-1 functional yields less accurate results, comparable with conventional DFT. These findings hold true for both asymmetrical and symmetrical shared-proton systems. Therefore, until a more accurate electron-proton correlation functional is developed, we currently recommend performing vibrational spectrum calculations using CNEO-DFT without electron-proton correlation functionals.

Local Electronic Correlation in Multicomponent Møller-Plesset Perturbation Theory

We present in this contribution the first application of local correlation in the context of multicomponent methods. Multicomponent approaches allow for the targeted simulation of electrons together with other Fermions (most commonly protons) as quantum particles. These methods have become increasingly popular over the last years, particularly for the description of nuclear quantum effects (in strong hydrogen bonds, proton tunneling, and many more). However, most implementations are still based on canonical formulations of wave function theory, which we know for decades to be computationally inefficient for capturing dynamical correlation effects. Local correlation approaches, particularly with the use of pair natural orbitals (PNOs), enable asymptotically linear scaling of computational costs with very little impact on the overall accuracy. In this context, the efficient use of density fitting approximations in the integral calculation proves essential. We start by discussing our implementation of density-fitted NEO-MP2 and NEO-PNO-LMP2, upgrading the electronic correlation treatment up to PNO local coupled cluster level of theory. Several challenging examples are provided to benchmark the method in terms of accuracy as well as computational cost scaling. Following appropriate protocols, anharmonic corrections to localized X-H stretches can be applied routinely with little computational overhead.

Optimization of Quantum Nuclei Positions with the Adaptive Nuclear-Electronic Orbital Approach

The use of multicomponent methods has become increasingly popular over the last years. Under this framework, nuclei (commonly protons) are treated quantum mechanically on the same footing as the electronic structure problem. Under the use of atomic-centered orbitals, this can lead to some complications as the ideal location of the nuclear basis centers must be optimized. In this contribution, we propose a straightforward approach to determine the position of such centers within the self-consistent cycle of a multicomponent calculation, making use of individual proton charge centroids. We test the method on model systems including the water dimer, a protonated water tetramer, and a porphine system. Comparing to numerical gradient calculations, the adaptive nuclear-electronic orbital (NEO) procedure is able to converge the basis centers to within a few cents of an Ångström and with less than 0.1 kcal/mol differences in absolute energies. This is achieved in one single calculation and with a small added computational effort of up to 80% compared to a regular NEO- self-consistent field run. An example application for the human transketolase proton wire is also provided.

Correlations between the Structures and Spectra of Protonated Water Clusters

Badger's rule-like correlations between OH stretching frequencies and intensities and the OH bond length are used to develop a spectral mapping procedure for studies of pure and protonated water clusters. This approach utilizes the vibrationally averaged OH bond lengths, which were obtained from diffusion Monte Carlo simulations that were performed using the general potential developed by Yu and Bowman. Good agreement is achieved between the spectra obtained using this approach and previously reported spectra for H+(H2O)n clusters, with n = 3, 4, and 5, as well as their perdeuterated analogues. The analysis of the spectra obtained by this spectral mapping approach supports previous work that assigned the spectrum of H+(H2O)6 to a mixture of Eigen and Zundel-like structures. Analysis of the calculated spectra also suggests a reassignment of the frequency of one of the transitions that involves the OH stretching vibration of the OH bonds in the hydronium core in the Eigen-like structure of H+(H2O)6 from 1917 cm-1 to roughly 2100 cm-1. For D+(D2O)6, comparison of the measured spectrum to those obtained by using the spectral mapping approach suggests that the carrier of the measured spectrum is one or more of the isomers of D+(D2O)6 that contain a four-membered ring and two flanking water molecules. While there are several candidate structures, the two flanking water molecules most likely form a chain that is bound to the hydronium core.