1
|
Barbier-Chebbah A, Bénichou O, Voituriez R, Guérin T. Long-term memory induced correction to Arrhenius law. Nat Commun 2024; 15:7408. [PMID: 39198409 PMCID: PMC11358423 DOI: 10.1038/s41467-024-50938-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2024] [Accepted: 07/25/2024] [Indexed: 09/01/2024] Open
Abstract
The Kramers escape problem is a paradigmatic model for the kinetics of rare events, which are usually characterized by Arrhenius law. So far, analytical approaches have failed to capture the kinetics of rare events in the important case of non-Markovian processes with long-term memory, as occurs in the context of reactions involving proteins, long polymers, or strongly viscoelastic fluids. Here, based on a minimal model of non-Markovian Gaussian process with long-term memory, we determine quantitatively the mean FPT to a rare configuration and provide its asymptotics in the limit of a large energy barrier E. Our analysis unveils a correction to Arrhenius law, induced by long-term memory, which we determine analytically. This correction, which we show can be quantitatively significant, takes the form of a second effective energy barrierE ' < E and captures the dependence of rare event kinetics on initial conditions, which is a hallmark of long-term memory. Altogether, our results quantify the impact of long-term memory on rare event kinetics, beyond Arrhenius law.
Collapse
Affiliation(s)
- A Barbier-Chebbah
- Decision and Bayesian Computation, USR 3756 (C3BI/DBC) and Neuroscience Department CNRS UMR 3751, Institut Pasteur, Université de Paris, CNRS, 75015, Paris, France
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France
| | - O Bénichou
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France.
| | - R Voituriez
- Laboratoire de Physique Théorique de la Matière Condensée, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France
- Laboratoire Jean Perrin, CNRS/UPMC, 4 Place Jussieu, 75005, Paris, France
| | - T Guérin
- Laboratoire Ondes et Matière d'Aquitaine, CNRS/University of Bordeaux, F-33400, Talence, France
| |
Collapse
|
2
|
Farahvash A, Willard AP. A theory of phonon-induced friction on molecular adsorbates. Proc Natl Acad Sci U S A 2024; 121:e2400589121. [PMID: 39052839 PMCID: PMC11295025 DOI: 10.1073/pnas.2400589121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2024] [Accepted: 06/07/2024] [Indexed: 07/27/2024] Open
Abstract
In this manuscript, we provide a general theory for how surface phonons couple to molecular adsorbates. Our theory maps the extended dynamics of a surface's atomic vibrational motions to a generalized Langevin equation, and by doing so captures these dynamics in a single quantity: the non-Markovian friction. The different frequency components of this friction are the phonon modes of the surface slab weighted by their coupling to the adsorbate degrees of freedom. Using this formalism, we demonstrate that physisorbed species couple primarily to acoustic phonons while chemisorbed species couple to dispersionless local vibrations. We subsequently derive equations for phonon-adjusted reaction rates using transition state theory and demonstrate that these corrections improve agreement with experimental results for CO desorption rates from Pt(111).
Collapse
Affiliation(s)
- Ardavan Farahvash
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA02139
| | - Adam P. Willard
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, MA02139
| |
Collapse
|
3
|
Wiśniewski M, Spiechowicz J. Memory-induced absolute negative mobility. CHAOS (WOODBURY, N.Y.) 2024; 34:073101. [PMID: 38949530 DOI: 10.1063/5.0213706] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/12/2024] [Accepted: 06/10/2024] [Indexed: 07/02/2024]
Abstract
Non-Markovian systems form a broad area of physics that remains greatly unexplored despite years of intensive investigations. The spotlight is on memory as a source of effects that are absent in their Markovian counterparts. In this work, we dive into this problem and analyze a driven Brownian particle moving in a spatially periodic potential and exposed to correlated thermal noise. We show that the absolute negative mobility effect, in which the net movement of the particle is in the direction opposite to the average force acting on it, may be induced by the memory of the setup. To explain the origin of this phenomenon, we resort to the recently developed effective mass approach to dynamics of non-Markovian systems.
Collapse
Affiliation(s)
- M Wiśniewski
- Institute of Physics, University of Silesia, 41-500 Chorzów, Poland
| | - J Spiechowicz
- Institute of Physics, University of Silesia, 41-500 Chorzów, Poland
| |
Collapse
|
4
|
Dalton BA, Kiefer H, Netz RR. The role of memory-dependent friction and solvent viscosity in isomerization kinetics in viscogenic media. Nat Commun 2024; 15:3761. [PMID: 38704367 PMCID: PMC11069540 DOI: 10.1038/s41467-024-48016-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2023] [Accepted: 04/18/2024] [Indexed: 05/06/2024] Open
Abstract
Molecular isomerization kinetics in liquid solvent depends on a complex interplay between the solvent friction acting on the molecule, internal dissipation effects (also known as internal friction), the viscosity of the solvent, and the dihedral free energy profile. Due to the absence of accurate techniques to directly evaluate isomerization friction, it has not been possible to explore these relationships in full. By combining extensive molecular dynamics simulations with friction memory-kernel extraction techniques we consider a variety of small, isomerising molecules under a range of different viscogenic conditions and directly evaluate the viscosity dependence of the friction acting on a rotating dihedral. We reveal that the influence of different viscogenic media on isomerization kinetics can be dramatically different, even when measured at the same viscosity. This is due to the dynamic solute-solvent coupling, mediated by time-dependent friction memory kernels. We also show that deviations from the linear dependence of isomerization rates on solvent viscosity, which are often simply attributed to internal friction effects, are due to the simultaneous violation of two fundamental relationships: the Stokes-Einstein relation and the overdamped Kramers prediction for the barrier-crossing rate, both of which require explicit knowledge of friction.
Collapse
Affiliation(s)
| | - Henrik Kiefer
- Freie Universität Berlin, Fachbereich Physik, Berlin, Germany
| | - Roland R Netz
- Freie Universität Berlin, Fachbereich Physik, Berlin, Germany.
| |
Collapse
|
5
|
Fernandez L, Hess S, Klapp SHL. Nonequilibrium dynamics and entropy production of a trapped colloidal particle in a complex nonreciprocal medium. Phys Rev E 2024; 109:054129. [PMID: 38907489 DOI: 10.1103/physreve.109.054129] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2024] [Accepted: 05/02/2024] [Indexed: 06/24/2024]
Abstract
We discuss the two-dimensional motion of a Brownian particle that is confined to a harmonic trap and driven by a shear flow. The surrounding medium induces memory effects modeled by a linear, typically nonreciprocal coupling of the particle coordinates to an auxiliary (hidden) variable. The system's behavior resulting from the microscopic Langevin equations for the three variables is analyzed by means of exact moment equations derived from the Fokker-Planck representation, and numerical Brownian dynamics simulations. Increasing the shear rate beyond a critical value we observe, for suitable coupling scenarios with nonreciprocal elements, a transition from a stationary to a nonstationary state, corresponding to an escape from the trap. We analyze this behavior, analytically and numerically, in terms of the associated moments of the probability distribution, and from the perspective of nonequilibrium thermodynamics. Intriguingly, the entropy production rate remains finite when crossing the stability threshold.
Collapse
|
6
|
Tepper L, Dalton B, Netz RR. Accurate Memory Kernel Extraction from Discretized Time-Series Data. J Chem Theory Comput 2024; 20:3061-3068. [PMID: 38603471 PMCID: PMC11044577 DOI: 10.1021/acs.jctc.3c01289] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2023] [Revised: 03/25/2024] [Accepted: 04/01/2024] [Indexed: 04/13/2024]
Abstract
Memory effects emerge as a fundamental consequence of dimensionality reduction when low-dimensional observables are used to describe the dynamics of complex many-body systems. In the context of molecular dynamics (MD) data analysis, accounting for memory effects using the framework of the generalized Langevin equation (GLE) has proven efficient, accurate, and insightful, particularly when working with high-resolution time series data. However, in experimental systems, high-resolution data are often unavailable, raising questions about the impact of the data resolution on the estimated GLE parameters. This study demonstrates that direct memory extraction from time series data remains accurate when the discretization time is below the memory time. To obtain memory functions reliably, even when the discretization time exceeds the memory time, we introduce a Gaussian Process Optimization (GPO) scheme. This scheme minimizes the deviation of discretized two-point correlation functions between time series data and GLE simulations and is able to estimate accurate memory kernels as long as the discretization time stays below the longest time scale in the data, typically the barrier crossing time.
Collapse
Affiliation(s)
- Lucas Tepper
- Department of Physics, Freie
Universität Berlin, 14195 Berlin, Germany
| | - Benjamin Dalton
- Department of Physics, Freie
Universität Berlin, 14195 Berlin, Germany
| | - Roland R. Netz
- Department of Physics, Freie
Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
7
|
Zhang XZ, Shi R, Lu ZY, Qian HJ. Chemically Specific Systematic Coarse-Grained Polymer Model with Both Consistently Structural and Dynamical Properties. JACS AU 2024; 4:1018-1030. [PMID: 38559727 PMCID: PMC10976574 DOI: 10.1021/jacsau.3c00756] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 11/30/2023] [Revised: 02/23/2024] [Accepted: 02/28/2024] [Indexed: 04/04/2024]
Abstract
The coarse-grained (CG) model serves as a powerful tool for the simulation of polymer systems; its reliability depends on the accurate representation of both structural and dynamical properties. However, strong correlations between structural and dynamical properties on different scales and also a strong memory effect, enforced by chain connectivity between monomers in polymer systems, render developing a chemically specific systematic CG model a formidable task. In this study, we report a systematic CG approach that combines the iterative Boltzmann inversion (IBI) method and the generalized Langevin equation (GLE) dynamics. Structural properties are ensured by using conservative CG potentials derived from the IBI method. To retrieve the correct dynamical properties in the system, we demonstrate that using a combination of a Rouse-type delta function and a time-dependent short-time kernel in the GLE simulation is practically efficient. The former can be used to adjust the long-time diffusion dynamics, and the latter can be reconstructed from an iterative procedure according to the velocity autocorrelation function (ACF) from all-atomistic (AA) simulations. Taking the polystyrene as an example, we show that not only structural properties of radial distribution function, intramolecular bond, and angle distributions can be reproduced but also dynamical properties of mean-square displacement, velocity ACF, and force ACF resulted from our CG model have quantitative agreement with the reference AA model. In addition, reasonable agreements are observed in other collective properties between our GLE-CG model and the AA simulations as well.
Collapse
Affiliation(s)
| | | | - Zhong-Yuan Lu
- State Key Laboratory of Supramolecular
Structure and Materials, Institute of Theoretical Chemistry, College
of Chemistry, Jilin University, Changchun 130021, China
| | - Hu-Jun Qian
- State Key Laboratory of Supramolecular
Structure and Materials, Institute of Theoretical Chemistry, College
of Chemistry, Jilin University, Changchun 130021, China
| |
Collapse
|
8
|
Brünig F, Daldrop JO, Netz RR. Pair-Reaction Dynamics in Water: Competition of Memory, Potential Shape, and Inertial Effects. J Phys Chem B 2022; 126:10295-10304. [PMID: 36473702 PMCID: PMC9761671 DOI: 10.1021/acs.jpcb.2c05923] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Revised: 11/11/2022] [Indexed: 12/12/2022]
Abstract
When described by a one-dimensional reaction coordinate, pair-reaction rates in a solvent depend, in addition to the potential barrier height and the friction coefficient, on the potential shape, the effective mass, and the friction relaxation spectrum, but a rate theory that accurately accounts for all of these effects does not exist. After a review of classical reaction-rate theories, we show how to extract all parameters of the generalized Langevin equation (GLE) and, in particular, the friction memory function from molecular dynamics (MD) simulations of two prototypical pair reactions in water, the dissociation of NaCl and of two methane molecules. The memory exhibits multiple time scales and, for NaCl, pronounced oscillatory components. Simulations of the GLE by Markovian embedding techniques accurately reproduce the pair-reaction kinetics from MD simulations without any fitting parameters, which confirms the accuracy of the approximative form of the GLE and of the parameter extraction techniques. By modification of the GLE parameters, we investigate the relative importance of memory, mass, and potential shape effects. Neglect of memory slows down NaCl and methane dissociation by roughly a factor of 2; neglect of mass accelerates reactions by a similar factor, and the harmonic approximation of the potential shape gives rise to slight acceleration. This partial error cancellation explains why Kramers' theory, which neglects memory effects and treats the potential shape in harmonic approximation, describes reaction rates better than more sophisticated theories. In essence, all three effects, friction memory, inertia, and the potential shape nonharmonicity, are important to quantitatively describe pair-reaction kinetics in water.
Collapse
Affiliation(s)
- Florian
N. Brünig
- Fachbereich Physik, Freie Universität
Berlin, Arnimallee 14, 14195Berlin, Germany
| | - Jan O. Daldrop
- Fachbereich Physik, Freie Universität
Berlin, Arnimallee 14, 14195Berlin, Germany
| | - Roland R. Netz
- Fachbereich Physik, Freie Universität
Berlin, Arnimallee 14, 14195Berlin, Germany
| |
Collapse
|
9
|
Brünig FN, Hillmann P, Kim WK, Daldrop JO, Netz RR. Proton-transfer spectroscopy beyond the normal-mode scenario. J Chem Phys 2022; 157:174116. [DOI: 10.1063/5.0116686] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022] Open
Abstract
A stochastic theory is developed to predict the spectral signature of proton-transfer processes and is applied to infrared spectra computed from ab initio molecular-dynamics simulations of a single [Formula: see text] cation. By constraining the oxygen atoms to a fixed distance, this system serves as a tunable model for general proton-transfer processes with variable barrier height. Three spectral contributions at distinct frequencies are identified and analytically predicted: the quasi-harmonic motion around the most probable configuration, amenable to normal-mode analysis, the contribution due to transfer paths when the proton moves over the barrier, and a shoulder for low frequencies stemming from the stochastic transfer-waiting-time distribution; the latter two contributions are not captured by normal-mode analysis but exclusively reported on the proton-transfer kinetics. In accordance with reaction rate theory, the transfer-waiting-contribution frequency depends inversely exponentially on the barrier height, whereas the transfer-path-contribution frequency is rather insensitive to the barrier height.
Collapse
Affiliation(s)
- Florian N. Brünig
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Paul Hillmann
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Won Kyu Kim
- School of Computational Sciences, Korea Institute for Advanced Study, Seoul 02455, Republic of Korea
| | - Jan O. Daldrop
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Roland R. Netz
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| |
Collapse
|
10
|
Brünig FN, Netz RR, Kappler J. Barrier-crossing times for different non-Markovian friction in well and barrier: A numerical study. Phys Rev E 2022; 106:044133. [PMID: 36397504 DOI: 10.1103/physreve.106.044133] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2022] [Accepted: 07/21/2022] [Indexed: 06/16/2023]
Abstract
We introduce a generalized Langevin model system for different non-Markovian effects in the well and barrier regions of a potential, and use it to numerically study the barrier-crossing time. In the appropriate limits, our model interpolates between the theoretical barrier-crossing-time predictions by Grote and Hynes (GH), as well as by Pollak et al., which for a single barrier memory time can differ by several orders of magnitude. Our model furthermore allows one to test an analytic rate theory for space-inhomogeneous memory, which disagrees with our numerical results in the long well-memory regime. In this regime, we find that short barrier memory decreases the barrier-crossing time as compared to long barrier memory. This is in contrast with the short well-memory regime, where both our numerical results and the GH theory predict an acceleration of the barrier crossing time with increasing barrier memory time. Both effects, the "Markovian-barrier acceleration" and GH "non-Markovian-barrier acceleration," can be understood from a committor analysis. Our model combines finite relaxation times of orthogonal degrees of freedom with a space-inhomogeneous coupling to such degrees and represents a step towards more realistic modeling of reaction coordinates.
Collapse
Affiliation(s)
- Florian N Brünig
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Roland R Netz
- Department of Physics, Freie Universität Berlin, 14195 Berlin, Germany
| | - Julian Kappler
- Department of Applied Mathematics and Theoretical Physics, Centre for Mathematical Sciences, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom
| |
Collapse
|
11
|
Cherayil BJ. Effects of Hydrodynamic Backflow on the Transmission Coefficient of a Barrier-Crossing Brownian Particle. J Phys Chem B 2022; 126:5629-5636. [PMID: 35894587 DOI: 10.1021/acs.jpcb.2c03273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The slow power law decay of the velocity autocorrelation function of a particle moving stochastically in a condensed-phase fluid is widely attributed to the momentum that fluid molecules displaced by the particle transfer back to it during the course of its motion. The forces created by this backflow effect are known as Basset forces, and they have been found in recent analytical work and numerical simulations to be implicated in a number of interesting dynamical phenomena, including boosted particle mobility in tilted washboard potentials. Motivated by these findings, the present paper is an investigation of the role of backflow in thermally activated barrier crossing, the governing process in essentially all condensed-phase chemical reactions. More specifically, it is an exact analytical calculation, carried out within the framework of the reactive-flux formalism, of the transmission coefficient κ(t) of a Brownian particle that crosses an inverted parabola under the influence of a colored noise process originating in the Basset force and a Markovian time-local friction. The calculation establishes that κ(t) is significantly enhanced over its backflow-free limit.
Collapse
Affiliation(s)
- Binny J Cherayil
- Department of Inorganic and Physical Chemistry, Indian Institute of Science, Bangalore 560012, Karnataka, India
| |
Collapse
|
12
|
Brünig FN, Rammler M, Adams EM, Havenith M, Netz RR. Spectral signatures of excess-proton waiting and transfer-path dynamics in aqueous hydrochloric acid solutions. Nat Commun 2022; 13:4210. [PMID: 35864099 PMCID: PMC9304333 DOI: 10.1038/s41467-022-31700-x] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2021] [Accepted: 06/22/2022] [Indexed: 11/23/2022] Open
Abstract
The theoretical basis for linking spectral signatures of hydrated excess protons with microscopic proton-transfer mechanisms has so far relied on normal-mode analysis. We introduce trajectory-decomposition techniques to analyze the excess-proton dynamics in ab initio molecular-dynamics simulations of aqueous hydrochloric-acid solutions beyond the normal-mode scenario. We show that the actual proton transfer between two water molecules involves for relatively large water-water separations crossing of a free-energy barrier and thus is not a normal mode, rather it is characterized by two non-vibrational time scales: Firstly, the broadly distributed waiting time for transfer to occur with a mean value of 200-300 fs, which leads to a broad and weak shoulder in the absorption spectrum around 100 cm-1, consistent with our experimental THz spectra. Secondly, the mean duration of a transfer event of about 14 fs, which produces a rather well-defined spectral contribution around 1200 cm-1 and agrees in location and width with previous experimental mid-infrared spectra.
Collapse
Affiliation(s)
- Florian N Brünig
- Freie Universität Berlin, Department of Physics, 14195, Berlin, Germany
| | - Manuel Rammler
- Freie Universität Berlin, Department of Physics, 14195, Berlin, Germany
| | - Ellen M Adams
- Ruhr-Universität Bochum, Department of Physical Chemistry II, 44780, Bochum, Germany
| | - Martina Havenith
- Ruhr-Universität Bochum, Department of Physical Chemistry II, 44780, Bochum, Germany
| | - Roland R Netz
- Freie Universität Berlin, Department of Physics, 14195, Berlin, Germany.
| |
Collapse
|
13
|
Neumann J, Schwierz N. Artificial Intelligence Resolves Kinetic Pathways of Magnesium Binding to RNA. J Chem Theory Comput 2022; 18:1202-1212. [PMID: 35084846 PMCID: PMC8830046 DOI: 10.1021/acs.jctc.1c00752] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Abstract
Magnesium is an indispensable cofactor in countless vital processes. In order to understand its functional role, the characterization of the binding pathways to biomolecules such as RNA is crucial. Despite the importance, a molecular description is still lacking since the transition from the water-mediated outer-sphere to the direct inner-sphere coordination is on the millisecond time scale and therefore out of reach for conventional simulation techniques. To fill this gap, we use transition path sampling to resolve the binding pathways and to elucidate the role of the solvent in the binding process. The results reveal that the molecular void provoked by the leaving phosphate oxygen of the RNA is immediately filled by an entering water molecule. In addition, water molecules from the first and second hydration shell couple to the concerted exchange. To capture the intimate solute-solvent coupling, we perform a committor analysis as the basis for a machine learning algorithm that derives the optimal deep learning model from thousands of scanned architectures using hyperparameter tuning. The results reveal that the properly optimized deep network architecture recognizes the important solvent structures, extracts the relevant information, and predicts the commitment probability with high accuracy. Our results provide detailed insights into the solute-solvent coupling which is ubiquitous for kosmotropic ions and governs a large variety of biochemical reactions in aqueous solutions.
Collapse
Affiliation(s)
- Jan Neumann
- Allianz Global Investors GmbH, Bockenheimer Landstrasse 42, 60323 Frankfurt am Main, Germany
| | - Nadine Schwierz
- Department of Theoretical Biophysics, Max-Planck-Institute of Biophysics, 60438 Frankfurt am Main, Germany
| |
Collapse
|
14
|
Schleeh MM, Reiff J, García-Müller PL, Benito RM, Borondo F, Main J, Hernandez R. Mean first-passage times for solvated LiCN isomerization at intermediate to high temperatures. J Chem Phys 2022; 156:034103. [DOI: 10.1063/5.0065090] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Micha M. Schleeh
- Institut für Theoretische Physik I, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Johannes Reiff
- Institut für Theoretische Physik I, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Pablo L. García-Müller
- Departamento de Tecnología, Centro de Investigaciones Energéticas Medioambientales y Tecnológicas, Avda. Complutense 40, 28040 Madrid, Spain
| | - Rosa M. Benito
- Grupo de Sistemas Complejos, Escuela Técnica Superior de Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, 28040 Madrid, Spain
| | - Florentino Borondo
- Instituto de Ciencias Matemáticas (ICMAT), Cantoblanco, 28049 Madrid, Spain
- Departamento de Química, Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
| | - Jörg Main
- Institut für Theoretische Physik I, Universität Stuttgart, 70550 Stuttgart, Germany
| | - Rigoberto Hernandez
- Department of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, USA
- Departments of Chemical and Biomolecular Engineering and Materials Science and Engineering, Johns Hopkins University, Baltimore, Maryland 21218, USA
| |
Collapse
|
15
|
Ginot F, Caspers J, Krüger M, Bechinger C. Barrier Crossing in a Viscoelastic Bath. PHYSICAL REVIEW LETTERS 2022; 128:028001. [PMID: 35089774 DOI: 10.1103/physrevlett.128.028001] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/31/2021] [Accepted: 11/30/2021] [Indexed: 06/14/2023]
Abstract
We investigate the hopping dynamics of a colloidal particle across a potential barrier and within a viscoelastic, i.e., non-Markovian, bath and report two clearly separated timescales in the corresponding waiting time distributions. While the longer timescale exponentially depends on the barrier height, the shorter one is similar to the relaxation time of the fluid. This short timescale is a signature of the storage and release of elastic energy inside the bath that strongly increases the hopping rate. Our results are in excellent agreement with numerical simulations of a simple Maxwell model.
Collapse
Affiliation(s)
- Félix Ginot
- Fachbereich Physik, Universität Konstanz, 78457 Konstanz, Germany
| | - Juliana Caspers
- Institute for Theoretical Physics, Georg-August Universität Göttingen, 37073 Göttingen, Germany
| | - Matthias Krüger
- Institute for Theoretical Physics, Georg-August Universität Göttingen, 37073 Göttingen, Germany
| | | |
Collapse
|
16
|
Abstract
Protein-folding kinetics is often described as Markovian (i.e., memoryless) diffusion in a one-dimensional free energy landscape, governed by an instantaneous friction coefficient that is fitted to reproduce experimental or simulated folding times. For the α-helix forming polypeptide alanine9 and a specific reaction coordinate that consists of the summed native hydrogen-bond lengths, we demonstrate that the friction extracted from molecular dynamics simulations exhibits significant memory with a decay time that is in the nanosecond range and thus, of the same order as the folding and unfolding times. Our non-Markovian modeling not only reproduces the molecular dynamics simulations accurately but also demonstrates that memory friction effects lead to anomalous and drastically accelerated protein kinetics. We extract the folding free energy landscape and the time-dependent friction function, the two ingredients of the generalized Langevin equation (GLE), from explicit-water molecular dynamics (MD) simulations of the α-helix forming polypeptide alanine9 for a one-dimensional reaction coordinate based on the sum of the native H-bond distances. Folding and unfolding times from numerical integration of the GLE agree accurately with MD results, which demonstrate the robustness of our GLE-based non-Markovian model. In contrast, Markovian models do not accurately describe the peptide kinetics and in particular, cannot reproduce the folding and unfolding kinetics simultaneously, even if a spatially dependent friction profile is used. Analysis of the GLE demonstrates that memory effects in the friction significantly speed up peptide folding and unfolding kinetics, as predicted by the Grote–Hynes theory, and are the cause of anomalous diffusion in configuration space. Our methods are applicable to any reaction coordinate and in principle, also to experimental trajectories from single-molecule experiments. Our results demonstrate that a consistent description of protein-folding dynamics must account for memory friction effects.
Collapse
|
17
|
Klippenstein V, van der Vegt NFA. Cross-correlation corrected friction in (generalized) Langevin models. J Chem Phys 2021; 154:191102. [PMID: 34240903 DOI: 10.1063/5.0049324] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a route for parameterizing isotropic (generalized) Langevin [(G)LE] thermostats with the aim to correct the dynamics of coarse-grained (CG) models with pairwise conservative interactions. The approach is based on the Mori-Zwanzig formalism and derives the memory kernels from Q-projected time correlation functions. Bottom-up informed (GLE and LE) thermostats for a CG star-polymer melt are investigated, and it is demonstrated that the inclusion of memory in the CG simulation leads to predictions of polymer diffusion in quantitative agreement with fine-grained simulations. Interestingly, memory effects are observed in the diffusive regime. We demonstrate that previously neglected cross-correlations between the "irrelevant" and the CG degree of freedom are important and lie at the origin of shortcomings in previous CG simulations.
Collapse
Affiliation(s)
- Viktor Klippenstein
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Nico F A van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| |
Collapse
|
18
|
Klippenstein V, Tripathy M, Jung G, Schmid F, van der Vegt NFA. Introducing Memory in Coarse-Grained Molecular Simulations. J Phys Chem B 2021; 125:4931-4954. [PMID: 33982567 PMCID: PMC8154603 DOI: 10.1021/acs.jpcb.1c01120] [Citation(s) in RCA: 40] [Impact Index Per Article: 13.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Preserving the correct dynamics at the coarse-grained (CG) level is a pressing problem in the development of systematic CG models in soft matter simulation. Starting from the seminal idea of simple time-scale mapping, there have been many efforts over the years toward establishing a meticulous connection between the CG and fine-grained (FG) dynamics based on fundamental statistical mechanics approaches. One of the most successful attempts in this context has been the development of CG models based on the Mori-Zwanzig (MZ) theory, where the resulting equation of motion has the form of a generalized Langevin equation (GLE) and closely preserves the underlying FG dynamics. In this Review, we describe some of the recent studies in this regard. We focus on the construction and simulation of dynamically consistent systematic CG models based on the GLE, both in the simple Markovian limit and the non-Markovian case. Some recent studies of physical effects of memory are also discussed. The Review is aimed at summarizing recent developments in the field while highlighting the major challenges and possible future directions.
Collapse
Affiliation(s)
- Viktor Klippenstein
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Madhusmita Tripathy
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| | - Gerhard Jung
- Institut für Theoretische Physik, Universität Innsbruck, Technikerstraße 21 A, A-6020 Innsbruck, Austria
| | - Friederike Schmid
- Institut für Physik, Johannes Gutenberg-Universität Mainz, Staudingerweg 9, 55128 Mainz, Germany
| | - Nico F A van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, 64287 Darmstadt, Germany
| |
Collapse
|
19
|
Ferrer BR, Gomez-Solano JR, Arzola AV. Fluid Viscoelasticity Triggers Fast Transitions of a Brownian Particle in a Double Well Optical Potential. PHYSICAL REVIEW LETTERS 2021; 126:108001. [PMID: 33784172 DOI: 10.1103/physrevlett.126.108001] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2020] [Revised: 02/03/2021] [Accepted: 02/09/2021] [Indexed: 06/12/2023]
Abstract
Thermally activated transitions are ubiquitous in nature, occurring in complex environments which are typically conceived as ideal viscous fluids. We report the first direct observations of a Brownian bead transiting between the wells of a bistable optical potential in a viscoelastic fluid with a single long relaxation time. We precisely characterize both the potential and the fluid, thus enabling a neat comparison between our experimental results and a theoretical model based on the generalized Langevin equation. Our findings reveal a drastic amplification of the transition rates compared to those in a Newtonian fluid, stemming from the relaxation of the fluid during the particle crossing events.
Collapse
Affiliation(s)
- Brandon R Ferrer
- Instituto de Física, Universidad Nacional Autónoma de México, Ciudad de México, Codigo Postal 04510, México
| | - Juan Ruben Gomez-Solano
- Instituto de Física, Universidad Nacional Autónoma de México, Ciudad de México, Codigo Postal 04510, México
| | - Alejandro V Arzola
- Instituto de Física, Universidad Nacional Autónoma de México, Ciudad de México, Codigo Postal 04510, México
| |
Collapse
|
20
|
Yamaguchi T. Molecular dynamics simulation study on the isomerization reaction in a solvent with slow structural relaxation. Chem Phys 2021. [DOI: 10.1016/j.chemphys.2020.111056] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
21
|
Mitterwallner BG, Schreiber C, Daldrop JO, Rädler JO, Netz RR. Non-Markovian data-driven modeling of single-cell motility. Phys Rev E 2021; 101:032408. [PMID: 32289977 DOI: 10.1103/physreve.101.032408] [Citation(s) in RCA: 22] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Accepted: 01/07/2020] [Indexed: 01/23/2023]
Abstract
Trajectories of human breast cancer cells moving on one-dimensional circular tracks are modeled by the non-Markovian version of the Langevin equation that includes an arbitrary memory function. When averaged over cells, the velocity distribution exhibits spurious non-Gaussian behavior, while single cells are characterized by Gaussian velocity distributions. Accordingly, the data are described by a linear memory model which includes different random walk models that were previously used to account for various aspects of cell motility such as migratory persistence, non-Markovian effects, colored noise, and anomalous diffusion. The memory function is extracted from the trajectory data without restrictions or assumptions, thus making our approach truly data driven, and is used for unbiased single-cell comparison. The cell memory displays time-delayed single-exponential negative friction, which clearly distinguishes cell motion from the simple persistent random walk model and suggests a regulatory feedback mechanism that controls cell migration. Based on the extracted memory function we formulate a generalized exactly solvable cell migration model which indicates that negative friction generates cell persistence over long timescales. The nonequilibrium character of cell motion is investigated by mapping the non-Markovian Langevin equation with memory onto a Markovian model that involves a hidden degree of freedom and is equivalent to the underdamped active Ornstein-Uhlenbeck process.
Collapse
Affiliation(s)
- Bernhard G Mitterwallner
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany and Physik Fakultät, Ludwig Maximilians Universität, 80539 München, Germany
| | - Christoph Schreiber
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany and Physik Fakultät, Ludwig Maximilians Universität, 80539 München, Germany
| | - Jan O Daldrop
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany and Physik Fakultät, Ludwig Maximilians Universität, 80539 München, Germany
| | - Joachim O Rädler
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany and Physik Fakultät, Ludwig Maximilians Universität, 80539 München, Germany
| | - Roland R Netz
- Fachbereich Physik, Freie Universität Berlin, 14195 Berlin, Germany and Physik Fakultät, Ludwig Maximilians Universität, 80539 München, Germany
| |
Collapse
|
22
|
Lapolla A, Godec A. Single-file diffusion in a bi-stable potential: Signatures of memory in the barrier-crossing of a tagged-particle. J Chem Phys 2020; 153:194104. [PMID: 33218229 DOI: 10.1063/5.0025785] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We investigate memory effects in barrier-crossing in the overdamped setting. We focus on the scenario where the hidden degrees of freedom relax on exactly the same time scale as the observable. As a prototypical model, we analyze tagged-particle diffusion in a single file confined to a bi-stable potential. We identify the signatures of memory and explain their origin. The emerging memory is a result of the projection of collective many-body eigenmodes onto the motion of a tagged-particle. We are interested in the "confining" (all background particles in front of the tagged-particle) and "pushing" (all background particles behind the tagged-particle) scenarios for which we find non-trivial and qualitatively different relaxation behaviors. Notably and somewhat unexpectedly, at a fixed particle number, we find that the higher the barrier, the stronger the memory effects are. The fact that the external potential alters the memory is important more generally and should be taken into account in applications of generalized Langevin equations. Our results can readily be tested experimentally and may be relevant for understanding transport in biological ion-channels.
Collapse
Affiliation(s)
- Alessio Lapolla
- Mathematical bioPhysics Group, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Göttingen, Germany
| | - Aljaž Godec
- Mathematical bioPhysics Group, Max Planck Institute for Biophysical Chemistry, Am Fassberg 11, 37077 Göttingen, Germany
| |
Collapse
|
23
|
Mitterwallner BG, Lavacchi L, Netz RR. Negative friction memory induces persistent motion. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2020; 43:67. [PMID: 33099707 DOI: 10.1140/epje/i2020-11992-5] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2020] [Accepted: 09/30/2020] [Indexed: 06/11/2023]
Abstract
We investigate the mean-square displacement (MSD) for random motion governed by the generalized Langevin equation for memory functions that contain two different time scales: In the first model, the memory kernel consists of a delta peak and a single-exponential and in the second model of the sum of two exponentials. In particular, we investigate the scenario where the long-time exponential kernel contribution is negative. The competition between positive and negative friction memory contributions produces an enhanced transient persistent regime in the MSD, which is relevant for biological motility and active matter systems.
Collapse
Affiliation(s)
| | - Laura Lavacchi
- Fachbereich Physik, Freie Universität Berlin, 14195, Berlin, Germany
| | - Roland R Netz
- Fachbereich Physik, Freie Universität Berlin, 14195, Berlin, Germany
| |
Collapse
|
24
|
Jack RL, Harris RJ. Giant leaps and long excursions: Fluctuation mechanisms in systems with long-range memory. Phys Rev E 2020; 102:012154. [PMID: 32795045 DOI: 10.1103/physreve.102.012154] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2020] [Accepted: 07/01/2020] [Indexed: 01/15/2023]
Abstract
We analyze large deviations of time-averaged quantities in stochastic processes with long-range memory, where the dynamics at time t depends itself on the value q_{t} of the time-averaged quantity. First we consider the elephant random walk and a Gaussian variant of this model, identifying two mechanisms for unusual fluctuation behavior, which differ from the Markovian case. In particular, the memory can lead to large-deviation principles with reduced speeds and to nonanalytic rate functions. We then explain how the mechanisms operating in these two models are generic for memory-dependent dynamics and show other examples including a non-Markovian simple exclusion process.
Collapse
Affiliation(s)
- Robert L Jack
- Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, United Kingdom.,Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Rosemary J Harris
- School of Mathematical Sciences, Queen Mary University of London, Mile End Road, London E1 4NS, United Kingdom
| |
Collapse
|
25
|
Meyer H, Pelagejcev P, Schilling T. Non-Markovian out-of-equilibrium dynamics: A general numerical procedure to construct time-dependent memory kernels for coarse-grained observables. ACTA ACUST UNITED AC 2020. [DOI: 10.1209/0295-5075/128/40001] [Citation(s) in RCA: 22] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022]
|
26
|
Pyo AGT, Woodside MT. Memory effects in single-molecule force spectroscopy measurements of biomolecular folding. Phys Chem Chem Phys 2019; 21:24527-24534. [PMID: 31663550 DOI: 10.1039/c9cp04197d] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Abstract
Folding is generally assumed to be a Markov process, without memory. When the molecular motion is coupled to that of a probe as in single-molecule force spectroscopy (SMFS) experiments, however, theory predicts that the coupling to a second Markov process should induce memory when monitoring a projection of the full multi-dimensional motion onto a reduced coordinate. We developed a method to evaluate the time constant of the induced memory from its effects on the autocorrelation function, which can be readily determined from experimental data. Applying this method to both simulated SMFS measurements and experimental trajectories of DNA hairpin folding measured by optical tweezers as a model system, we validated the prediction that the linker induces memory. For these measurements, the timescale of the induced memory was found to be similar to the time required for the force probe to respond to changes in the molecule, and in the regime where the experimentally observed dynamics were not significantly perturbed by probe-molecule coupling artifacts. Memory effects are thus a general feature of SMFS measurements induced by the mechanical connection between the molecule and force probe that should be considered when interpreting experimental data.
Collapse
Affiliation(s)
- Andrew G T Pyo
- Department of Physics, University of Alberta, Edmonton, AB T6G 2E1, Canada.
| | | |
Collapse
|
27
|
Kappler J, Hinrichsen VB, Netz RR. Non-Markovian barrier crossing with two-time-scale memory is dominated by the faster memory component. THE EUROPEAN PHYSICAL JOURNAL. E, SOFT MATTER 2019; 42:119. [PMID: 31494784 DOI: 10.1140/epje/i2019-11886-7] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/30/2019] [Accepted: 08/06/2019] [Indexed: 06/10/2023]
Abstract
We investigate non-Markovian barrier-crossing kinetics of a massive particle in one dimension in the presence of a memory function that is the sum of two exponentials with different memory times, [Formula: see text] and [Formula: see text] . Our Langevin simulations for the special case where both exponentials contribute equally to the total friction show that the barrier-crossing time becomes independent of the longer memory time if at least one of the two memory times is larger than the intrinsic diffusion time. When we associate memory effects with coupled degrees of freedom that are orthogonal to a one-dimensional reaction coordinate, this counterintuitive result shows that the faster orthogonal degrees of freedom dominate barrier-crossing kinetics in the non-Markovian limit and that the slower orthogonal degrees become negligible, quite contrary to the standard time-scale separation assumption and with important consequences for the proper setup of coarse-graining procedures in the non-Markovian case. By asymptotic matching and symmetry arguments, we construct a crossover formula for the barrier crossing time that is valid for general multi-exponential memory kernels. This formula can be used to estimate barrier-crossing times for general memory functions for high friction, i.e. in the overdamped regime, as well as for low friction, i.e. in the inertial regime. Typical examples where our results are important include protein folding in the high-friction limit and chemical reactions such as proton-transfer reactions in the low-friction limit.
Collapse
Affiliation(s)
- Julian Kappler
- Freie Universität Berlin, Fachbereich Physik, Berlin, Germany
| | | | - Roland R Netz
- Freie Universität Berlin, Fachbereich Physik, Berlin, Germany.
| |
Collapse
|
28
|
Lee HS, Ahn SH, Darve EF. The multi-dimensional generalized Langevin equation for conformational motion of proteins. J Chem Phys 2019; 150:174113. [PMID: 31067888 DOI: 10.1063/1.5055573] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022] Open
Abstract
Using the generalized Langevin equation (GLE) is a promising approach to build coarse-grained (CG) models of molecular systems since the GLE model often leads to more accurate thermodynamic and kinetic predictions than Brownian dynamics or Langevin models by including a more sophisticated friction with memory. The GLE approach has been used for CG coordinates such as the center of mass of a group of atoms with pairwise decomposition and for a single CG coordinate. We present a GLE approach when CG coordinates are multiple generalized coordinates, defined, in general, as nonlinear functions of microscopic atomic coordinates. The CG model for multiple generalized coordinates is described by the multidimensional GLE from the Mori-Zwanzig formalism, which includes an exact memory matrix. We first present a method to compute the memory matrix in a multidimensional GLE using trajectories of a full system. Then, in order to reduce the computational cost of computing the multidimensional friction with memory, we introduce a method that maps the GLE to an extended Markovian system. In addition, we study the effect of using a nonconstant mass matrix in the CG model. In particular, we include mass-dependent terms in the mean force. We used the proposed CG model to describe the conformational motion of a solvated alanine dipeptide system, with two dihedral angles as the CG coordinates. We showed that the CG model can accurately reproduce two important kinetic quantities: the velocity autocorrelation function and the distribution of first passage times.
Collapse
Affiliation(s)
- Hee Sun Lee
- Mechanical Engineering Department, Stanford University, Stanford, California 94305, USA
| | - Surl-Hee Ahn
- Chemistry Department, Stanford University, Stanford, California 94305, USA
| | - Eric F Darve
- Mechanical Engineering Department, Stanford University, Stanford, California 94305, USA
| |
Collapse
|
29
|
Deichmann G, van der Vegt NFA. Bottom-up approach to represent dynamic properties in coarse-grained molecular simulations. J Chem Phys 2019; 149:244114. [PMID: 30599732 DOI: 10.1063/1.5064369] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Several molecular coarse-graining methods have been proposed in recent years to derive chemical- and state-point transferable force fields. While these force fields describe structural and thermodynamic properties in good agreement with fine-grained models and experiments, dynamic properties are usually overestimated. Herein, we examine if the long-time dynamic properties of molecular coarse-grained (CG) systems can be correctly represented by employing a dissipative particle dynamics (DPD) thermostat, which is "bottom-up informed" by means of a variant of the Markovian Mori-Zwanzig (MZ) DPD coarse-graining method. We report single-site and multiple-site CG models for a monomer, dimer, and 24mer based on 2,2-dimethyl propane as a chemical repeat unit and report data obtained from MZ-DPD simulations of liquids, polymer solutions, and polymer melts. We find that despite incomplete time scale separation of the molecular CG model, MZ-DPD achieves quantitative accuracy in predicting diffusive dynamics in single-component liquids and polymer solutions (24mers in a dimer solvent). We also find that MZ-DPD simulations of molecular penetrant diffusion in polymer networks do not reach quantitative agreement with the fine-grained model. Modeling diffusion governed by the activated barrier crossing of small molecular penetrants in these dense systems requires an accurate description of energy barriers, presumably combined with the treatment of memory effects. The use of a MZ-DPD thermostat extends the scope and applicability of molecular CG models for multicomponent systems where a correct description of the relative diffusion rates of the different components is important.
Collapse
Affiliation(s)
- Gregor Deichmann
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, 64287 Darmstadt, Germany
| | - Nico F A van der Vegt
- Eduard-Zintl-Institut für Anorganische und Physikalische Chemie, Technische Universität Darmstadt, Alarich-Weiss-Straße 10, 64287 Darmstadt, Germany
| |
Collapse
|
30
|
Hridya VM, Mukherjee A. Probing the Viscosity Dependence of Rate: Internal Friction or the Lack of Friction? J Phys Chem B 2018; 122:9081-9086. [DOI: 10.1021/acs.jpcb.8b05585] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- V. M. Hridya
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| | - Arnab Mukherjee
- Department of Chemistry, Indian Institute of Science Education and Research, Pune 411008, India
| |
Collapse
|
31
|
Abstract
The dihedral dynamics of butane in water is known to be rather insensitive to the water viscosity; possible explanations for this involve inertial effects or Kramers' turnover, the finite memory time of friction, and the presence of so-called internal friction. To disentangle these factors, we introduce a method to directly extract the friction memory function from unconstrained simulations in the presence of an arbitrary free-energy landscape. By analysis of the dihedral friction in butane for varying water viscosity, we demonstrate the existence of an internal friction contribution that does not scale linearly with water viscosity. At normal water viscosity, the internal friction turns out to be eight times larger than the solvent friction and thus completely dominates the effective friction. By comparison with simulations of a constrained butane molecule that has the dihedral as the only degree of freedom, we show that internal friction comes from the six additional degrees of freedom in unconstrained butane that are orthogonal to the dihedral angle reaction coordinate. While the insensitivity of butane's dihedral dynamics to water viscosity is solely due to the presence of internal friction, inertial effects nevertheless crucially influence the resultant transition rates. In contrast, non-Markovian effects due to the finite memory time are present but do not significantly influence the dihedral barrier-crossing rate of butane. These results not only settle the character of dihedral dynamics in small solvated molecular systems such as butane, they also have important implications for the folding of polymers and proteins.
Collapse
|