TY - JOUR
T1 - Physics-guided estimation of mean first-passage times from censored nucleation trajectories
AU - Huang, Liang Yao
AU - Lai, Pin Kuang
AU - Lin, Shiang Tai
N1 - Publisher Copyright:
© 2025 Author(s).
PY - 2025/7/28
Y1 - 2025/7/28
N2 - Mean first-passage time (MFPT) analysis is a powerful tool for extracting thermodynamic and kinetic parameters of nucleation, including the critical nucleus size, nucleation barrier, and nucleation rate, from molecular dynamics (MD) simulations. However, accurate MFPT estimation typically requires extensive sampling and long simulation durations, making it computationally expensive, especially for rare nucleation events. Here, we present an efficient method to estimate MFPT from fixed-length MD simulations by leveraging the single-exponential tail (SET) behavior of first-passage time distributions in monomer-transition-based nucleation processes. Our approach allows for the estimation of MFPT for censored samples using information from uncensored ones via the characteristic decay rate of the SET, eliminating the need to observe complete nucleation trajectories. We validate this method using both artificial nucleation datasets and unbiased MD simulations of CO2 hydrate nucleation. The results demonstrate that, with an appropriate choice of simulation length and number of independent runs, the SET-based approach accurately recovers the full MFPT and nucleation parameters while substantially reducing computational costs. This work provides a practical and broadly applicable framework for accelerating nucleation studies in molecular simulations.
AB - Mean first-passage time (MFPT) analysis is a powerful tool for extracting thermodynamic and kinetic parameters of nucleation, including the critical nucleus size, nucleation barrier, and nucleation rate, from molecular dynamics (MD) simulations. However, accurate MFPT estimation typically requires extensive sampling and long simulation durations, making it computationally expensive, especially for rare nucleation events. Here, we present an efficient method to estimate MFPT from fixed-length MD simulations by leveraging the single-exponential tail (SET) behavior of first-passage time distributions in monomer-transition-based nucleation processes. Our approach allows for the estimation of MFPT for censored samples using information from uncensored ones via the characteristic decay rate of the SET, eliminating the need to observe complete nucleation trajectories. We validate this method using both artificial nucleation datasets and unbiased MD simulations of CO2 hydrate nucleation. The results demonstrate that, with an appropriate choice of simulation length and number of independent runs, the SET-based approach accurately recovers the full MFPT and nucleation parameters while substantially reducing computational costs. This work provides a practical and broadly applicable framework for accelerating nucleation studies in molecular simulations.
UR - https://www.scopus.com/pages/publications/105011588918
UR - https://www.scopus.com/pages/publications/105011588918#tab=citedBy
U2 - 10.1063/5.0280948
DO - 10.1063/5.0280948
M3 - Article
C2 - 40693614
AN - SCOPUS:105011588918
SN - 0021-9606
VL - 163
JO - Journal of Chemical Physics
JF - Journal of Chemical Physics
IS - 4
M1 - 044101
ER -