Logo image
Transition-Path Theory Calculations on Non-Uniform Meshes in Two and Three Dimensions using Finite Elements
Journal article   Open access   Peer reviewed

Transition-Path Theory Calculations on Non-Uniform Meshes in Two and Three Dimensions using Finite Elements

Mauro Lapelosa and Cameron F Abrams
Computer physics communications, v 184(10), pp 2310-2315
Oct 2013
PMID: 24014889
url
https://doi.org/10.1016/j.cpc.2013.05.017View
Published, Version of Record (VoR) Open

Abstract

rare events committor function reaction kinetics
Rare events between states in complex systems are fundamental in many scientific fields and can be studied by building reaction pathways. A theoretical framework to analyze reaction pathways is provided by transition-path theory (TPT). The central object in TPT is the committor function, which is found by solution of the backward-Kolmogorov equation on a given potential. Once determined, the committor can be used to calculate reactive fluxes and rates, among other important quantities. We demonstrate here that the committor can be calculated using the method of finite elements on non-uniform meshes. We show that this approach makes it feasible to perform TPT calculations on 3D potentials because it requires many fewer degrees of freedom than a regular-mesh finite-difference approach. In various illustrative 2D and 3D problems, we calculate the committor function and reaction rates at different temperatures, and we discuss effects of temperatures and simple entropic barriers on the structure of the committor and the reaction rate constants.

Metrics

5 Record Views
2 citations in Scopus

Details

UN Sustainable Development Goals (SDGs)

This publication has contributed to the advancement of the following goals:

#3 Good Health and Well-Being

InCites Highlights

Data related to this publication, from InCites Benchmarking & Analytics tool:

Web of Science research areas
Computer Science, Interdisciplinary Applications
Physics, Mathematical
Logo image