1
|
Anderson MC, Dodin A, Fay TP, Limmer DT. Coherent control from quantum commitment probabilities. J Chem Phys 2024; 161:024115. [PMID: 38995082 DOI: 10.1063/5.0213444] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2024] [Accepted: 06/24/2024] [Indexed: 07/13/2024] Open
Abstract
We introduce a general definition of a quantum committor in order to clarify reaction mechanisms and facilitate control in processes where coherent effects are important. With a quantum committor, we generalize the notion of a transition state to quantum superpositions and quantify the effect of interference on the progress of the reaction. The formalism is applicable to any linear quantum master equation supporting metastability for which absorbing boundary conditions designating the reactant and product states can be applied. We use this formalism to determine the dependence of the quantum transition state on coherences in a polaritonic system and optimize the initialization state of a conical intersection model to control reactive outcomes, achieving yields of the desired state approaching 100%. In addition to providing a practical tool, the quantum committor provides a conceptual framework for understanding reactions in cases when classical intuitions fail.
Collapse
Affiliation(s)
- Michelle C Anderson
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Amro Dodin
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Thomas P Fay
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - David T Limmer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Kavli Energy NanoSciences Institute, University of California, Berkeley, California 94720, USA
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
2
|
Tokar VI. Effective Hamiltonian approach to kinetic Ising models: Application to an infinitely long-range Husimi-Temperley model. Phys Rev E 2024; 109:044123. [PMID: 38755929 DOI: 10.1103/physreve.109.044123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/06/2023] [Accepted: 03/14/2024] [Indexed: 05/18/2024]
Abstract
The probability distribution (PD) of spin configurations in kinetic Ising models has been cast in the form of the canonical Boltzmann PD with a time-dependent effective Hamiltonian (EH). It has been argued that in systems with extensive energy EH depends linearly on the number of spins N leading to the exponential dependence of PD on the system size. In macroscopic systems the argument of the exponential function may reach values of the order of the Avogadro number which is impossible to deal with computationally, thus making unusable the linear master equation (ME) governing the PD evolution. To overcome the difficulty, it has been suggested to use instead the nonlinear ME (NLME) for the EH density per spin. It has been shown that in spatially homogeneous systems NLME contains only terms of order unity even in the thermodynamic limit. The approach has been illustrated with the kinetic Husimi-Temperley model (HTM) evolving under the Glauber dynamics. At finite N the known numerical results has been reproduced and extended to broader parameter ranges. In the thermodynamic limit an exact nonlinear partial differential equation of the Hamilton-Jacobi type for EH has been derived. It has been shown that the average magnetization in HTM evolves according to the conventional kinetic mean field equation.
Collapse
Affiliation(s)
- V I Tokar
- G. V. Kurdyumov Institute for Metal Physics of the N.A.S. of Ukraine, 36 Acad. Vernadsky Boulevard, UA-03142 Kyiv, Ukraine and Université de Strasbourg, CNRS, IPCMS, UMR 7504, F-67000 Strasbourg, France
| |
Collapse
|
3
|
Dorbath E, Gulzar A, Stock G. Log-periodic oscillations as real-time signatures of hierarchical dynamics in proteins. J Chem Phys 2024; 160:074103. [PMID: 38364004 DOI: 10.1063/5.0188220] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/20/2023] [Accepted: 01/23/2024] [Indexed: 02/18/2024] Open
Abstract
The time-dependent relaxation of a dynamical system may exhibit a power-law behavior that is superimposed by log-periodic oscillations. D. Sornette [Phys. Rep. 297, 239 (1998)] showed that this behavior can be explained by a discrete scale invariance of the system, which is associated with discrete and equidistant timescales on a logarithmic scale. Examples include such diverse fields as financial crashes, random diffusion, and quantum topological materials. Recent time-resolved experiments and molecular dynamics simulations suggest that discrete scale invariance may also apply to hierarchical dynamics in proteins, where several fast local conformational changes are a prerequisite for a slow global transition to occur. Employing entropy-based timescale analysis and Markov state modeling to a simple one-dimensional hierarchical model and biomolecular simulation data, it is found that hierarchical systems quite generally give rise to logarithmically spaced discrete timescales. By introducing a one-dimensional reaction coordinate that collectively accounts for the hierarchically coupled degrees of freedom, the free energy landscape exhibits a characteristic staircase shape with two metastable end states, which causes the log-periodic time evolution of the system. The period of the log-oscillations reflects the effective roughness of the energy landscape and can, in simple cases, be interpreted in terms of the barriers of the staircase landscape.
Collapse
Affiliation(s)
- Emanuel Dorbath
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| | - Adnan Gulzar
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| |
Collapse
|
4
|
Anderson MC, Woods EJ, Fay TP, Wales DJ, Limmer DT. On the Mechanism of Polaritonic Rate Suppression from Quantum Transition Paths. J Phys Chem Lett 2023:6888-6894. [PMID: 37494137 DOI: 10.1021/acs.jpclett.3c01188] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Polariton chemistry holds promise for facilitating mode-selective chemical reactions, but the underlying mechanism behind the rate modifications observed under strong vibrational coupling is not well understood. Using the recently developed quantum transition path theory, we have uncovered a mechanism of resonant suppression of a thermal reaction rate in a simple model polaritonic system consisting of a reactive mode in a bath confined to a lossless microcavity with a single photon mode. We observed the formation of a polariton during rate-limiting transitions on reactive pathways and identified the concomitant rate suppression as being due to hybridization between the reactive mode and the cavity mode, which inhibits bath-mediated tunneling. The transition probabilities that define the quantum master equation can be directly translated into a visualization of the corresponding polariton energy landscape. This landscape exhibits a double funnel structure with a large barrier between the initial and final states.
Collapse
Affiliation(s)
- Michelle C Anderson
- Department of Chemistry, University of California, Berkeley 94720, United States
| | - Esmae J Woods
- Department of Physics, University of Cambridge, Cambridge CB3 0HE, U.K
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, U.K
| | - Thomas P Fay
- Department of Chemistry, University of California, Berkeley 94720, United States
| | - David J Wales
- Department of Chemistry, University of Cambridge, Cambridge CB2 1EW, U.K
| | - David T Limmer
- Department of Chemistry, University of California, Berkeley 94720, United States
- Kavli Energy NanoSciences Institute, University of California, Berkeley 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley 94720, United States
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley 94720, United States
| |
Collapse
|
5
|
Nagel D, Sartore S, Stock G. Selecting Features for Markov Modeling: A Case Study on HP35. J Chem Theory Comput 2023. [PMID: 37167425 DOI: 10.1021/acs.jctc.3c00240] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/13/2023]
Abstract
Markov state models represent a popular means to interpret molecular dynamics trajectories in terms of memoryless transitions between metastable conformational states. To provide a mechanistic understanding of the considered biomolecular process, these states should reflect structurally distinct conformations and ensure a time scale separation between fast intrastate and slow interstate dynamics. Adopting the folding of villin headpiece (HP35) as a well-established model problem, here we discuss the selection of suitable input coordinates or "features", such as backbone dihedral angles and interresidue distances. We show that dihedral angles account accurately for the structure of the native energy basin of HP35, while the unfolded region of the free energy landscape and the folding process are best described by tertiary contacts of the protein. To construct a contact-based model, we consider various ways to define and select contact distances and introduce a low-pass filtering of the feature trajectory as well as a correlation-based characterization of states. Relying on input data that faithfully account for the mechanistic origin of the studied process, the states of the resulting Markov model are clearly discriminated by the features, describe consistently the hierarchical structure of the free energy landscape, and─as a consequence─correctly reproduce the slow time scales of the process.
Collapse
Affiliation(s)
- Daniel Nagel
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| | - Sofia Sartore
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| | - Gerhard Stock
- Biomolecular Dynamics, Institute of Physics, University of Freiburg, 79104 Freiburg, Germany
| |
Collapse
|
6
|
Koskin V, Kells A, Clayton J, Hartmann AK, Annibale A, Rosta E. Variational kinetic clustering of complex networks. J Chem Phys 2023; 158:104112. [PMID: 36922127 DOI: 10.1063/5.0105099] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Efficiently identifying the most important communities and key transition nodes in weighted and unweighted networks is a prevalent problem in a wide range of disciplines. Here, we focus on the optimal clustering using variational kinetic parameters, linked to Markov processes defined on the underlying networks, namely, the slowest relaxation time and the Kemeny constant. We derive novel relations in terms of mean first passage times for optimizing clustering via the Kemeny constant and show that the optimal clustering boundaries have equal round-trip times to the clusters they separate. We also propose an efficient method that first projects the network nodes onto a 1D reaction coordinate and subsequently performs a variational boundary search using a parallel tempering algorithm, where the variational kinetic parameters act as an energy function to be extremized. We find that maximization of the Kemeny constant is effective in detecting communities, while the slowest relaxation time allows for detection of transition nodes. We demonstrate the validity of our method on several test systems, including synthetic networks generated from the stochastic block model and real world networks (Santa Fe Institute collaboration network, a network of co-purchased political books, and a street network of multiple cities in Luxembourg). Our approach is compared with existing clustering algorithms based on modularity and the robust Perron cluster analysis, and the identified transition nodes are compared with different notions of node centrality.
Collapse
Affiliation(s)
- Vladimir Koskin
- Department of Chemistry, King's College London, SE1 1DB London, United Kingdom
| | - Adam Kells
- Department of Chemistry, King's College London, SE1 1DB London, United Kingdom
| | - Joe Clayton
- Department of Physics and Astronomy, University College London, WC1E 6BT London, United Kingdom
| | | | - Alessia Annibale
- Department of Mathematics, King's College London, SE11 6NJ London, United Kingdom
| | - Edina Rosta
- Department of Physics and Astronomy, University College London, WC1E 6BT London, United Kingdom
| |
Collapse
|
7
|
Anderson MC, Schile AJ, Limmer DT. Nonadiabatic transition paths from quantum jump trajectories. J Chem Phys 2022; 157:164105. [DOI: 10.1063/5.0102891] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a means of studying rare reactive pathways in open quantum systems using transition path theory and ensembles of quantum jump trajectories. This approach allows for the elucidation of reactive paths for dissipative, nonadiabatic dynamics when the system is embedded in a Markovian environment. We detail the dominant pathways and rates of thermally activated processes and the relaxation pathways and photoyields following vertical excitation in a minimal model of a conical intersection. We find that the geometry of the conical intersection affects the electronic character of the transition state as defined through a generalization of a committor function for a thermal barrier crossing event. Similarly, the geometry changes the mechanism of relaxation following a vertical excitation. Relaxation in models resulting from small diabatic coupling proceeds through pathways dominated by pure dephasing, while those with large diabatic coupling proceed through pathways limited by dissipation. The perspective introduced here for the nonadiabatic dynamics of open quantum systems generalizes classical notions of reactive paths to fundamentally quantum mechanical processes.
Collapse
Affiliation(s)
- Michelle C. Anderson
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Addison J. Schile
- Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - David T. Limmer
- Department of Chemistry, University of California, Berkeley, California 94720, USA
- Kavli Energy NanoSciences Institute, University of California, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
- Materials Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
8
|
Yang X, Lu ZY. Nanoparticle cluster formation mechanisms elucidated via Markov state modeling: Attraction range effects, aggregation pathways, and counterintuitive transition rates. J Chem Phys 2022; 156:214902. [DOI: 10.1063/5.0086110] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Nanoparticle clusters are promising candidates for developing functional materials. However, it is still a challenging task to fabricate them in a predictable and controllable way, which requires investigation of the possible mechanisms underlying cluster formation at the nanoscale. By constructing Markov state models (MSMs) at the microstate level, we find that for highly dispersed particles to form a highly aggregated cluster, there are multiple coexisting pathways, which correspond to direct aggregation, or pathways that need to pass through partially aggregated, intermediate states. Varying the range of attraction between nanoparticles is found to significantly affect pathways. As the attraction range becomes narrower, compared to direct aggregation, some pathways that need to pass through partially aggregated intermediate states become more competitive. In addition, from MSMs constructed at the macrostate level, the aggregation rate is found to be counterintuitively lower with a lower free-energy barrier, which is also discussed.
Collapse
Affiliation(s)
- Xi Yang
- Institute of Theoretical Chemistry, State Key Laboratory of Supramolecular Structure and Materials, Jilin University, Changchun 130021, China
| | - Zhong-Yuan Lu
- Institute of Theoretical Chemistry, State Key Laboratory of Supramolecular Structure and Materials, Jilin University, Changchun 130021, China
| |
Collapse
|
9
|
Pineda M, Stamatakis M. Kinetic Monte Carlo simulations for heterogeneous catalysis: Fundamentals, current status, and challenges. J Chem Phys 2022; 156:120902. [DOI: 10.1063/5.0083251] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022] Open
Abstract
Kinetic Monte Carlo (KMC) simulations in combination with first-principles (1p)-based calculations are rapidly becoming the gold-standard computational framework for bridging the gap between the wide range of length scales and time scales over which heterogeneous catalysis unfolds. 1p-KMC simulations provide accurate insights into reactions over surfaces, a vital step toward the rational design of novel catalysts. In this Perspective, we briefly outline basic principles, computational challenges, successful applications, as well as future directions and opportunities of this promising and ever more popular kinetic modeling approach.
Collapse
Affiliation(s)
- M. Pineda
- Thomas Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London WC1E 7JE, United Kingdom
| | - M. Stamatakis
- Thomas Young Centre and Department of Chemical Engineering, University College London, Roberts Building, Torrington Place, London WC1E 7JE, United Kingdom
| |
Collapse
|
10
|
Röder K, Wales DJ. The Energy Landscape Perspective: Encoding Structure and Function for Biomolecules. Front Mol Biosci 2022; 9:820792. [PMID: 35155579 PMCID: PMC8829389 DOI: 10.3389/fmolb.2022.820792] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 01/07/2022] [Indexed: 12/02/2022] Open
Abstract
The energy landscape perspective is outlined with particular reference to biomolecules that perform multiple functions. We associate these multifunctional molecules with multifunnel energy landscapes, illustrated by some selected examples, where understanding the organisation of the landscape has provided new insight into function. Conformational selection and induced fit may provide alternative routes to realisation of multifunctionality, exploiting the possibility of environmental control and distinct binding modes.
Collapse
|