An official website of the United States government
Transition path theory computes statistics from ensembles of reactive trajectories. A common strategy for sampling reactive trajectories is to control the branching and pruning of trajectories so as to enhance the sampling of low probability segments. However, it can be challenging to apply transition path theory to data from such methods because determining whether configurations and trajectory segments are part of reactive trajectories requires looking backward and forward in time. Here, we show how this issue can be overcome efficiently by introducing simple data structures. We illustrate the approach in the context of nonequilibrium umbrella sampling, but the strategy is general and can be used to obtain transition path theory statistics from other methods that sample segments of unbiased trajectories.
more »
« less
Augmented transition path theory for sequences of events
Transition path theory provides a statistical description of the dynamics of a reaction in terms of local spatial quantities. In its original formulation, it is limited to reactions that consist of trajectories flowing from a reactant set A to a product set B. We extend the basic concepts and principles of transition path theory to reactions in which trajectories exhibit a specified sequence of events and illustrate the utility of this generalization on examples.
more »
« less
- Award ID(s):
- 2054306
- PAR ID:
- 10444872
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 157
- Issue:
- 9
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract Dynamic motion often controls selectivity in reactions featuring two consecutive potential‐energy transition states. Here we report density functional theory (DFT)‐based direct dynamics trajectories and machine learning classification analysis for cyclopentadienone dimerization and a N 2 extrusion reaction leading to semibullvalene. These reactions have consecutive transition states, and there is dynamic selectivity that determines which of two possible C‐C bonds is formed after the first transition state. For cyclopentadienone dimerization with a bispericyclic first transition state, machine learning analysis using transition‐state based features provided >90% trajectory classification accuracy, but only using AdaBoost and random forest algorithms. Many other relatively sophisticated machine learning algorithms showed poor accuracy despite the obvious motion responsible for selectivity. Feature importance analysis confirmed that the sigmatropic rearrangement vibrational motion in the bispericyclic transition state provides prediction of which of the second C‐C bonds is dynamically formed. For the reaction leading to semibullvalene, machine learning analysis provides solid accuracy for classifying trajectories and predicting which C‐C bond is formed and which C‐C bond is broken immediately after N 2 ejection. Like the cyclopentadienone dimerization reaction, machine learning feature importance analysis showed that the sigmatropic rearrangement vibrational motion in the N 2 extrusion transition state determines which C‐C bond is formed and which is broken. Surprisingly, machine learning struggles to predict which trajectories undergo a subsequent [3,3] sigmatropic rearrangement process, which isomerizes equivalent forms of semibullvalene.more » « less
-
Ulam’s method is a popular discretization scheme for stochastic operators that involves the construction of a transition probability matrix controlling a Markov chain on a set of cells covering some domain. We consider an application to satellite-tracked undrogued surface-ocean drifting buoy trajectories obtained from the National Oceanic and Atmospheric Administration Global Drifter Program dataset. Motivated by the motion of Sargassum in the tropical Atlantic, we apply Transition Path Theory (TPT) to drifters originating off the west coast of Africa to the Gulf of Mexico. We find that the most common case of a regular covering by equal longitude–latitude side cells can lead to a large instability in the computed transition times as a function of the number of cells used. We propose a different covering based on a clustering of the trajectory data that is stable against the number of cells in the covering. We also propose a generalization of the standard transition time statistic of TPT that can be used to construct a partition of the domain of interest into weakly dynamically connected regions.more » « less
-
null (Ed.)An understanding of the role that spin states play in semiconductor surface chemical reactions is currently limited. Herein, we provide evidence of a nonadiabatic reaction involving a localized singlet to triplet thermal excitation of the Si(100) surface dimer dangling bond. By comparing the β-hydrogen elimination kinetics of ethyl adsorbates probed by thermal desorption experiments to electronic structure calculation results, we determined that a coverage-dependent change in mechanism occurs. At low coverage, a nonadiabatic, inter-dimer mechanism is dominant, while adiabatic mechanisms become dominant at higher coverage. Computational results indicate that the spin crossover is rapid near room temperature and the nonadiabatic path is accelerated by a barrier that is 40 kJ mol −1 less than the adiabatic path. Simulated thermal desorption reactions using nonadiabatic transition state theory (NA-TST) for the surface dimer intersystem crossing are in close agreement with experimental observations.more » « less
-
Abstract Advances in single-cell technologies allow scrutinizing of heterogeneous cell states, however, detecting cell-state transitions from snap-shot single-cell transcriptome data remains challenging. To investigate cells with transient properties or mixed identities, we present MuTrans, a method based on multiscale reduction technique to identify the underlying stochastic dynamics that prescribes cell-fate transitions. By iteratively unifying transition dynamics across multiple scales, MuTrans constructs the cell-fate dynamical manifold that depicts progression of cell-state transitions, and distinguishes stable and transition cells. In addition, MuTrans quantifies the likelihood of all possible transition trajectories between cell states using coarse-grained transition path theory. Downstream analysis identifies distinct genes that mark the transient states or drive the transitions. The method is consistent with the well-established Langevin equation and transition rate theory. Applying MuTrans to datasets collected from five different single-cell experimental platforms, we show its capability and scalability to robustly unravel complex cell fate dynamics induced by transition cells in systems such as tumor EMT, iPSC differentiation and blood cell differentiation. Overall, our method bridges data-driven and model-based approaches on cell-fate transitions at single-cell resolution.more » « less
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- https://doi.org/10.1063/5.0098587
