1
|
Naik S, Wiggins S. Finding normally hyperbolic invariant manifolds in two and three degrees of freedom with Hénon-Heiles-type potential. Phys Rev E 2019; 100:022204. [PMID: 31574621 DOI: 10.1103/physreve.100.022204] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2019] [Indexed: 11/07/2022]
Abstract
We present a method based on a Lagrangian descriptor for revealing the high-dimensional phase space structures that are of interest in nonlinear Hamiltonian systems with index-1 saddle. These phase space structures include a normally hyperbolic invariant manifold and its stable and unstable manifolds, which act as codimension-1 barriers to phase space transport. In this article, finding the invariant manifolds in high-dimensional phase space will constitute identifying coordinates on these invariant manifolds. The method of Lagrangian descriptor is demonstrated by applying to classical two and three degrees of freedom Hamiltonian systems which have implications for myriad applications in chemistry, engineering, and physics.
Collapse
Affiliation(s)
- Shibabrat Naik
- School of Mathematics, University of Bristol, University Walk, Clifton BS8 1TW, Bristol, United Kingdom
| | - Stephen Wiggins
- School of Mathematics, University of Bristol, University Walk, Clifton BS8 1TW, Bristol, United Kingdom
| |
Collapse
|
2
|
Bartsch T, Revuelta F, Benito RM, Borondo F. Finite-barrier corrections for multidimensional barriers in colored noise. Phys Rev E 2019; 99:052211. [PMID: 31212507 DOI: 10.1103/physreve.99.052211] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2018] [Indexed: 06/09/2023]
Abstract
The usual identification of reactive trajectories for the calculation of reaction rates requires very time-consuming simulations, particularly if the environment presents memory effects. In this paper, we develop a method that permits the identification of reactive trajectories in a system under the action of a stochastic colored driving. This method is based on the perturbative computation of the invariant structures that act as separatrices for reactivity. Furthermore, using this perturbative scheme, we have obtained a formally exact expression for the reaction rate in multidimensional systems coupled to colored noisy environments.
Collapse
Affiliation(s)
- Thomas Bartsch
- Centre for Nonlinear Mathematics and Applications, Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, England, United Kingdom
| | - F Revuelta
- Grupo de Sistemas Complejos, E. T. S. Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Avenida Puerta de Hierro 2-4, 28040 Madrid, Spain
- Instituto de Ciencias Matemáticas, Cantoblanco, 28049 Madrid, Spain
| | - R M Benito
- Grupo de Sistemas Complejos, E. T. S. Ingeniería Agronómica, Alimentaria y de Biosistemas, Universidad Politécnica de Madrid, Avenida Puerta de Hierro 2-4, 28040 Madrid, Spain
| | - F Borondo
- Instituto de Ciencias Matemáticas, Cantoblanco, 28049 Madrid, Spain
- Departamento de Química, Universidad Autónoma de Madrid, Cantoblanco, 28049 Madrid, Spain
| |
Collapse
|
3
|
Craven GT, Bartsch T, Hernandez R. Persistence of transition-state structure in chemical reactions driven by fields oscillating in time. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2014; 89:040801. [PMID: 24827174 DOI: 10.1103/physreve.89.040801] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/27/2013] [Indexed: 06/03/2023]
Abstract
Chemical reactions subjected to time-varying external forces cannot generally be described through a fixed bottleneck near the transition-state barrier or dividing surface. A naive dividing surface attached to the instantaneous, but moving, barrier top also fails to be recrossing-free. We construct a moving dividing surface in phase space over a transition-state trajectory. This surface is recrossing-free for both Hamiltonian and dissipative dynamics. This is confirmed even for strongly anharmonic barriers using simulation. The power of transition-state theory is thereby applicable to chemical reactions and other activated processes even when the bottlenecks are time dependent and move across space.
Collapse
Affiliation(s)
- Galen T Craven
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| | - Thomas Bartsch
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, United Kingdom
| | - Rigoberto Hernandez
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400, USA
| |
Collapse
|
4
|
Bartsch T, Revuelta F, Benito RM, Borondo F. Reaction rate calculation with time-dependent invariant manifolds. J Chem Phys 2012; 136:224510. [DOI: 10.1063/1.4726125] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
5
|
Revuelta F, Bartsch T, Benito RM, Borondo F. Communication: Transition state theory for dissipative systems without a dividing surface. J Chem Phys 2012; 136:091102. [DOI: 10.1063/1.3692182] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
6
|
|
7
|
Time-Dependent Transition State Theory. ADVANCES IN CHEMICAL PHYSICS 2008. [DOI: 10.1002/9780470371572.ch4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|
8
|
Bartsch T, Uzer T, Moix JM, Hernandez R. Transition-State Theory Rate Calculations with a Recrossing-Free Moving Dividing Surface. J Phys Chem B 2007; 112:206-12. [DOI: 10.1021/jp0755600] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Thomas Bartsch
- Department of Mathematical Sciences, Loughborough University, Loughborough LE11 3TU, UK
| | - T. Uzer
- Center for Nonlinear Science, Georgia Institute of Technology, Atlanta, Georgia 30332-0430
| | - Jeremy M. Moix
- Chemical Physics Department, The Weizmann Institute of Science, 76100 Rehovot, Israel
| | - Rigoberto Hernandez
- Center for Computational Molecular Sciences & Technology, School of Chemistry & Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332-0400
| |
Collapse
|
9
|
Bartsch T, Uzer T, Moix JM, Hernandez R. Identifying reactive trajectories using a moving transition state. J Chem Phys 2006; 124:244310. [PMID: 16821980 DOI: 10.1063/1.2206587] [Citation(s) in RCA: 42] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A time-dependent no-recrossing dividing surface is shown to lead to a new criterion for identifying reactive trajectories well before they are evolved to infinite time. Numerical dynamics simulations of a dissipative anharmonic two-dimensional system confirm the efficiency of this approach. The results are compared to the standard fixed transition state dividing surface that is well-known to suffer from recrossings and therefore requires trajectories to be evolved over a long time interval before they can reliably be classified as reactive or nonreactive. The moving dividing surface can be used to identify reactive trajectories in harmonic or moderately anharmonic systems with considerably lower numerical effort or even without any simulation at all.
Collapse
Affiliation(s)
- Thomas Bartsch
- Center for Nonlinear Science, Georgia Institute of Technology, Atlanta, Georgia 30332-0430, USA
| | | | | | | |
Collapse
|
10
|
Hansen P, Mitchell KA, Delos JB. Escape of trajectories from a vase-shaped cavity. PHYSICAL REVIEW. E, STATISTICAL, NONLINEAR, AND SOFT MATTER PHYSICS 2006; 73:066226. [PMID: 16906965 DOI: 10.1103/physreve.73.066226] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2006] [Indexed: 05/11/2023]
Abstract
We consider the escape of ballistic trajectories from an open, vase-shaped cavity. Such a system serves as a model for microwaves escaping from a cavity or electrons escaping from a microjunction. Fixing the initial position of a particle and recording its escape time as a function of the initial launch direction, the resulting escape-time plot shows "epistrophic fractal" structure--repeated structure within structure at all levels of resolution with new features appearing in the fractal at longer time scales. By launching trajectories simultaneously in all directions (modeling an outgoing wave), a detector placed outside the cavity would measure a train of escaping pulses. We connect the structure of this chaotic pulse train with the fractal structure of the escape-time plot.
Collapse
Affiliation(s)
- Paul Hansen
- Department of Physics, College of William and Mary, Williamsburg, Virginia 23187-8795, USA.
| | | | | |
Collapse
|
11
|
Waalkens H, Burbanks A, Wiggins S. A formula to compute the microcanonical volume of reactive initial conditions in transition state theory. ACTA ACUST UNITED AC 2005. [DOI: 10.1088/0305-4470/38/45/l03] [Citation(s) in RCA: 29] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
|
12
|
Waalkens H, Burbanks A, Wiggins S. Efficient procedure to compute the microcanonical volume of initial conditions that lead to escape trajectories from a multidimensional potential well. PHYSICAL REVIEW LETTERS 2005; 95:084301. [PMID: 16196860 DOI: 10.1103/physrevlett.95.084301] [Citation(s) in RCA: 21] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/21/2004] [Revised: 03/10/2005] [Indexed: 05/04/2023]
Abstract
A procedure is presented for computing the phase space volume of initial conditions for trajectories that escape or "react" from a multidimensional potential well. The procedure combines a phase space transition state theory, which allows one to construct dividing surfaces that are free of local recrossing and that minimize the directional flux, and a classical spectral theorem. The procedure gives the volume of reactive initial conditions in terms of a sum over each entrance channel of the well of the product of the phase space flux across the dividing surface associated with the channel and the mean residence time in the well of trajectories which enter through the channel. This approach is illustrated for HCN isomerization in three dimensions, for which the method is several orders of magnitude more efficient than standard Monte Carlo sampling.
Collapse
Affiliation(s)
- Holger Waalkens
- School of Mathematics, University Walk, University of Bristol, Bristol BS8 1TW, United Kingdom
| | | | | |
Collapse
|
13
|
Bartsch T, Hernandez R, Uzer T. Transition state in a noisy environment. PHYSICAL REVIEW LETTERS 2005; 95:058301. [PMID: 16090923 DOI: 10.1103/physrevlett.95.058301] [Citation(s) in RCA: 63] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/21/2005] [Indexed: 05/03/2023]
Abstract
Transition state theory overestimates reaction rates in solution because conventional dividing surfaces between reagents and products are crossed many times by the same reactive trajectory. We describe a recipe for constructing a time-dependent dividing surface free of such recrossings in the presence of noise. The no-recrossing limit of transition state theory thus becomes generally available for the description of reactions in a fluctuating environment.
Collapse
Affiliation(s)
- Thomas Bartsch
- Center for Nonlinear Science, Georgia Institute of Technology, Atlanta, Georgia 30332-0430, USA
| | | | | |
Collapse
|
14
|
Waalkens H, Wiggins S. Direct construction of a dividing surface of minimal flux for multi-degree-of-freedom systems that cannot be recrossed. ACTA ACUST UNITED AC 2004. [DOI: 10.1088/0305-4470/37/35/l02] [Citation(s) in RCA: 91] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
15
|
Waalkens H, Burbanks A, Wiggins S. A computational procedure to detect a new type of high-dimensional chaotic saddle and its application to the 3D Hill's problem. ACTA ACUST UNITED AC 2004. [DOI: 10.1088/0305-4470/37/24/l04] [Citation(s) in RCA: 44] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
|
16
|
Wiggins S, Wiesenfeld L, Jaffé C, Uzer T. Impenetrable barriers in phase-space. PHYSICAL REVIEW LETTERS 2001; 86:5478-5481. [PMID: 11415280 DOI: 10.1103/physrevlett.86.5478] [Citation(s) in RCA: 87] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/29/2000] [Indexed: 05/23/2023]
Abstract
Dynamical systems theory is used to construct a general phase-space version of transition state theory. Special multidimensional separatrices are found which act as impenetrable barriers in phase-space between reacting and nonreacting trajectories. The elusive momentum-dependent transition state between reactants and products is thereby characterized. A practical algorithm is presented and applied to a strongly coupled Hamiltonian.
Collapse
Affiliation(s)
- S Wiggins
- School of Mathematics, University of Bristol, Bristol BS8 1TW, United Kingdom
| | | | | | | |
Collapse
|