skip to main content

Title: Challenges in constructing accurate methods for hydrogen transfer reactions in large biological assemblies: rare events sampling for mechanistic discovery and tensor networks for quantum nuclear effects
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 more » 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. « less
Award ID(s):
Publication Date:
Journal Name:
Faraday Discussions
Page Range or eLocation-ID:
379 to 405
Sponsoring Org:
National Science Foundation
More Like this
  1. Hydrogen tunneling plays a critical role in many biologically and chemically important processes. The nuclear–electronic orbital multistate density functional theory (NEO-MSDFT) method was developed to describe hydrogen transfer systems. In this approach, the transferring proton is treated quantum mechanically on the same level as the electrons within multicomponent DFT, and a nonorthogonal configuration interaction scheme is used to produce delocalized vibronic states from localized vibronic states. The NEO-MSDFT method has been shown to provide accurate hydrogen tunneling splittings for fixed molecular systems. Herein, the NEO-MSDFT analytical gradients for both ground and excited vibronic states are derived and implemented. The analytical gradients and semi-numerical Hessians are used to optimize and characterize equilibrium and transition state geometries and to generate minimum energy paths (MEPs), for proton transfer in the deprotonated acetylene dimer and malonaldehyde. The barriers along the resulting MEPs are lower when the transferring proton is quantized because the NEO-MSDFT method inherently includes the zero-point energy of the transferring proton. Analysis of the proton densities along the MEPs illustrates that the proton density can exhibit symmetric or asymmetric bilobal character associated with symmetric or slightly asymmetric double-well potential energy surfaces and hydrogen tunneling. Analysis of the contributions to the intrinsic reactionmore »coordinate reveals that changes in the C–O bond lengths drive proton transfer in malonaldehyde. This work provides the foundation for future reaction path studies and direct nonadiabatic dynamics simulations of a wide range of hydrogen transfer reactions.

    « less
  2. 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, predictionsmore »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.« less
  3. Ultrafast excited state processes of transition metal complexes (TMCs) are governed by complicated interplays between electronic and nuclear dynamics, which demand a detailed understanding to achieve optimal functionalities of photoactive TMC-based materials for many applications. In this work, we investigated a cyclometalated platinum( ii ) dimer known to undergo a Pt–Pt bond contraction in the metal–metal-to-ligand-charge-transfer (MMLCT) excited state using femtosecond broadband transient absorption (fs-BBTA) spectroscopy in combination with geometry optimization and normal mode calculations. Using a sub-20 fs pump and broadband probe pulses in fs-BBTA spectroscopy, we were able to correlate the coherent vibrational wavepacket (CVWP) evolution with the stimulated emission (SE) dynamics of the 1 MMLCT state. The results demonstrated that the 145 cm −1 CVWP motions with the damping times of ∼0.9 ps and ∼2 ps originate from coherent Pt–Pt stretching vibrations in the singlet and triplet MMLCT states, respectively. On the basis of excited state potential energy surface calculations in our previous work, we rationalized that the CVWP transfer from the Franck–Condon (FC) state to the 3 MMLCT state was mediated by a triplet ligand-centered ( 3 LC) intermediate state through two step intersystem crossing (ISC) on a time scale shorter than a period of themore »Pt–Pt stretching wavepacket motions. Moreover, it was found that the CVWP motion had 110 cm −1 frequency decays with the damping time of ∼0.2 ps, matching the time constant of 0.253 ps, corresponding to a redshift in the SE feature at early times. This observation indicates that the Pt–Pt bond contraction changes the stretching frequency from 110 to 145 cm −1 and stabilizes the 1 MMLCT state relative to the 3 LC state with a ∼0.2 ps time scale. Thus, the ultrafast ISC from the 1 MMLCT to the 3 LC states occurs before the Pt–Pt bond shortening. The findings herein provide insight into understanding the impact of Pt–Pt bond contraction on the ultrafast branching of the 1 MMLCT population into the direct ( 1 MMLCT → 3 MMLCT) and indirect ISC pathways ( 1 MMLCT → 3 LC → 3 MMLCT) in the Pt( ii ) dimer. These results revealed intricate excited state electronic and nuclear motions that could steer the reaction pathways with a level of detail that has not been achieved before.« less
  4. Abstract

    STIRAP (stimulated Raman adiabatic passage) is a powerful laser-based method, usually involving two photons, for efficient and selective transfer of populations between quantum states. A particularly interesting feature is the fact that the coupling between the initial and the final quantum states is via an intermediate state, even though the lifetime of the latter can be much shorter than the interaction time with the laser radiation. Nevertheless, spontaneous emission from the intermediate state is prevented by quantum interference. Maintaining the coherence between the initial and final state throughout the transfer process is crucial. STIRAP was initially developed with applications in chemical dynamics in mind. That is why the original paper of 1990 was published inThe Journal of Chemical Physics. However, from about the year 2000, the unique capabilities of STIRAP and its robustness with respect to small variations in some experimental parameters stimulated many researchers to apply the scheme to a variety of other fields of physics. The successes of these efforts are documented in this collection of articles. In Part A the experimental success of STIRAP in manipulating or controlling molecules, photons, ions or even quantum systems in a solid-state environment is documented. After a brief introduction tomore »the basic physics of STIRAP, the central role of the method in the formation of ultracold molecules is discussed, followed by a presentation of how precision experiments (measurement of the upper limit of the electric dipole moment of the electron or detecting the consequences of parity violation in chiral molecules) or chemical dynamics studies at ultralow temperatures benefit from STIRAP. Next comes the STIRAP-based control of photons in cavities followed by a group of three contributions which highlight the potential of the STIRAP concept in classical physics by presenting data on the transfer of waves (photonic, magnonic and phononic) between respective waveguides. The works on ions or ion strings discuss options for applications, e.g. in quantum information. Finally, the success of STIRAP in the controlled manipulation of quantum states in solid-state systems, which are usually hostile towards coherent processes, is presented, dealing with data storage in rare-earth ion doped crystals and in nitrogen vacancy (NV) centers or even in superconducting quantum circuits. The works on ions and those involving solid-state systems emphasize the relevance of the results for quantum information protocols. Part B deals with theoretical work, including further concepts relevant to quantum information or invoking STIRAP for the manipulation of matter waves. The subsequent articles discuss the experiments underway to demonstrate the potential of STIRAP for populating otherwise inaccessible high-lying Rydberg states of molecules, or controlling and cooling the translational motion of particles in a molecular beam or the polarization of angular-momentum states. The series of articles concludes with a more speculative application of STIRAP in nuclear physics, which, if suitable radiation fields become available, could lead to spectacular results.

    « less
  5. Photochemistry powers numerous processes from luminescence and human vision, to light harvesting. However, the elucidation of multidimensional photochemical reaction coordinates on molecular timescales remains challenging. We developed wavelength-tunable femtosecond stimulated Raman spectroscopy to simultaneously achieve pre-resonance enhancement for transient reactant and product species of the widely used photoacid pyranine undergoing excited-state proton transfer (ESPT) reaction in solution. In the low-frequency region, the 280 cm −1 ring deformation mode following 400 nm photoexcitation exhibits pronounced intensity oscillations on the sub-picosecond timescale due to anharmonic vibrational coupling to the 180 cm −1 hydrogen-bond stretching mode only in ESPT-capable solvents, indicating a primary event of functional relevance. This leads to the contact ion pair formation on the 3 ps timescale before diffusion-controlled separation. The intermolecular 180 cm −1 mode also reveals vibrational cooling time constants, ∼500 fs and 45 ps in both H 2 O and D 2 O, which differ from ESPT time constants of ∼3/8 and 90/250 ps in H 2 O/D 2 O, respectively. Spectral results using H 2 18 O further substantiate the functional role of the intermolecular 180 cm −1 mode in modulating the distance between proton donor and acceptor and forming the transient ion pair. The directmore »observation of molecular structural evolution across a wide spectral region during photochemical reactions enriches our fundamental understanding of potential energy surface and holds the key to advancing energy and biological sciences with exceptional atomic and temporal precision.« less