Journal
JOURNAL OF CHEMICAL PHYSICS
Volume 142, Issue 4, Pages -Publisher
AIP Publishing
DOI: 10.1063/1.4905957
Keywords
-
Funding
- U.S. Department of Energy Office of Science, Office of Advanced Scientific Computing Research, and Applied Mathematics program [DE-SC0010549]
- U.S. Department of Energy (DOE) [DE-SC0010549] Funding Source: U.S. Department of Energy (DOE)
Ask authors/readers for more resources
Kinetic Monte Carlo simulation is an integral tool in the study of complex physical phenomena present in applications ranging from heterogeneous catalysis to biological systems to crystal growth and atmospheric sciences. Sensitivity analysis is useful for identifying important parameters and rate-determining steps, but the finite-difference application of sensitivity analysis is computationally demanding. Techniques based on the likelihood ratio method reduce the computational cost of sensitivity analysis by obtaining all gradient information in a single run. However, we show that disparity in time scales of microscopic events, which is ubiquitous in real systems, introduces drastic statistical noise into derivative estimates for parameters affecting the fast events. In this work, the steady-state likelihood ratio sensitivity analysis is extended to singularly perturbed systems by invoking partial equilibration for fast reactions, that is, by working on the fast and slow manifolds of the chemistry. Derivatives on each time scale are computed independently and combined to the desired sensitivity coefficients to considerably reduce the noise in derivative estimates for stiff systems. The approach is demonstrated in an analytically solvable linear system. (C) 2015 AIP Publishing LLC.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available