An official website of the United States government
We present two methods that address the computational complexities arising in hydrogen transfer reactions in enzyme active sites. To address the challenge of reactive rare events, we begin with an ab initio molecular dynamics adaptation of the Caldeira–Leggett system-bath Hamiltonian and apply this approach to the study of the hydrogen transfer rate-determining step in soybean lipoxygenase-1. Through direct application of this method to compute an ensemble of classical trajectories, we discuss the critical role of isoleucine-839 in modulating the primary hydrogen transfer event in SLO-1. Notably, the formation of the hydrogen bond between isoleucine-839 and the acceptor-OH group regulates the electronegativity of the donor and acceptor groups to affect the hydrogen transfer process. Curtailing the formation of this hydrogen bond adversely affects the probability of hydrogen transfer. The second part of this paper deals with complementing the rare event sampled reaction pathways obtained from the aforementioned development through quantum nuclear wavepacket dynamics. Essentially the idea is to construct quantum nuclear dynamics on the potential surfaces obtained along the biased trajectories created as noted above. Here, while we are able to obtain critical insights on the quantum nuclear effects from wavepacket dynamics, we primarily engage in providing an improved computational approach for efficient representation of quantum dynamics data such as potential surfaces and transmission probabilities using tensor networks. We find that utilizing tensor networks yields an accurate and efficient description of time-dependent wavepackets, reduced dimensional nuclear eigenstates and associated potential energy surfaces at much reduced cost.
more »
« less
Rare Events Sampling Methods for Quantum and Classical Ab Initio Molecular Dynamics
We provide an approach to sample rare events during classical ab initio molecular dynamics and quantum wavepacket dynamics. For classical AIMD, a set of fictitious degrees of freedom are introduced that may harmonically interact with the electronic and nuclear degrees of freedom to steer the dynamics in a conservative fashion toward energetically forbidden regions. A similar approach when introduced for quantum wavepacket dynamics has the effect of biasing the trajectory of the wavepacket centroid toward the regions of the potential surface that are difficult to sample. The approach is demonstrated for a phenol-amine system, which is a prototypical problem for condensed phase-proton transfer, and for model potentials undergoing wavepacket dynamics. In all cases, the approach yields trajectories that conserve energy while sampling rare events.
more »
« less
- PAR ID:
- 10608139
- Publisher / Repository:
- American Chemical Society
- Date Published:
- Journal Name:
- The Journal of Physical Chemistry A
- Volume:
- 128
- Issue:
- 27
- ISSN:
- 1089-5639
- Page Range / eLocation ID:
- 5386 to 5397
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Canonical methods of quasiclassical dynamics make it possible to go beyond a strict background approximation for cosmological perturbations by including independent fields such as correlation degrees of freedom. New models are introduced and analyzed here for cosmological dynamics in the presence of quantum correlations between background and perturbations, as well as cross-correlations between different modes of a quantum field. Evolution equations for moments of a perturbation state reveal conditions required for inhomogeneity to build up out of an initial vacuum. A crucial role is played by quantum non-locality, formulated by canonical methods as an equivalent local theory with non-classical degrees of freedom given by moments of a quantum state.more » « less
-
Abstract The technologies used in the manipulation of light can be used to do analogue simulations of physical systems with wave-like equations of motion. This analogy is maximized by the use of all the degrees of freedom of light. The Helmholtz equation in physical optics and the Schodinger equation in quantum mechanics share the same mathematical form. We use this connection to prepare non-diffracting optical beams representing the spatial and temporal dynamics of a nonlinear physical system: the quantum pendulum. By using the propagation coordinate to represent time in the quantum problem, we are able to analogue-simulate quantum wavepacket dynamics. These manifest themselves in novel optical beams with rich three-dimensional structures, such as rotation and sloshing of the light's intensity as it propagates. Our experimental results agree very well with the predictions from quantum theory, thus demonstrating that our system can be used as a platform to simulate the quantum pendulum dynamics. This three-dimensional light-sculpting capability has the potential to impact fields such as manipulation with light and imaging.more » « less
-
Sodium hydride (NaH) in the gas phase presents a seemingly simple electronic structure making it a potentially tractable system for the detailed investigation of nonadiabatic molecular dynamics from both computational and experimental standpoints. The single vibrational degree of freedom, as well as the strong nonadiabatic coupling that arises from the excited electronic states taking on considerable ionic character, provides a realistic chemical system to test the accuracy of quasi-classical methods to model population dynamics where the results are directly comparable against quantum mechanical benchmarks. Using a simulated pump–probe type experiment, this work presents computational predictions of population transfer through the avoided crossings of NaH via symmetric quasi-classical Meyer–Miller (SQC/MM), Ehrenfest, and exact quantum dynamics on realistic, ab initio potential energy surfaces. The main driving force for population transfer arises from the ground vibrational level of the D 1 Σ + adiabatic state that is embedded in the manifold of near-dissociation C 1 Σ + vibrational states. When coupled through a sharply localized first-order derivative coupling most of the population transfers between t = 15 and t = 30 fs depending on the initially excited vibronic wavepacket. While quantum mechanical effects are expected due to the reduced mass of NaH, predictions of the population dynamics from both the SQC/MM and Ehrenfest models perform remarkably well against the quantum dynamics benchmark. Additionally, an analysis of the vibronic structure in the nonadiabatically coupled regime is presented using a variational eigensolver methodology.more » « less
-
The accurate computational study of wavepacketnuclear dynamics is considered to be a classically intractableproblem, particularly with increasing dimensions. Here, we presenttwo algorithms that, in conjunction with other methods developedby us, may result in one set of contributions for performingquantum nuclear dynamics in arbitrary dimensions. For one of thetwo algorithms discussed here, we present a direct map betweenthe Born−Oppenheimer Hamiltonian describing the nuclearwavepacket time evolution and the control parameters of a spin−lattice Hamiltonian that describes the dynamics of qubit states in anion-trap quantum computer. This map is exact for three qubits, andwhen implemented, the dynamics of the spin states emulates thoseof the nuclear wavepacket in a continuous representation. However, this map becomes approximate as the number of qubits grows.In a second algorithm, we present a general quantum circuit decomposition formalism for such problems using a method called theQuantum Shannon Decomposition. This algorithm is more robust and is exact for any number of qubits at the cost of increasedcircuit complexity. The resultant circuit is implemented on IBM’s quantum simulator (QASM) for 3−7 qubits, without using a noisemodel so as to test the intrinsic accuracy of the method. In both cases, the wavepacket dynamics is found to be in good agreementwith the classical propagation result and the corresponding vibrational frequencies obtained from the wavepacket density timeevolution are in agreement to within a few tenths of a wavenumber.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1021/acs.jpca.3c07385
