Hostname: page-component-848d4c4894-5nwft Total loading time: 0 Render date: 2024-05-23T09:41:48.564Z Has data issue: false hasContentIssue false

L1 Regularization-Based Model Reduction of Complex Chemistry Molecular Dynamics for Statistical Learning of Kinetic Monte Carlo Models

Published online by Cambridge University Press:  11 February 2016

Qian Yang
Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305, U.S.A.
Carlos A. Sing-Long
Institute for Computational and Mathematical Engineering, Stanford University, Stanford, CA 94305, U.S.A.
Evan J. Reed*
Department of Materials Science and Engineering, Stanford University, Stanford, CA 94305, U.S.A.
Get access


Kinetic Monte Carlo (KMC) methods have been a successful technique for accelerating time scales and increasing system sizes beyond those achievable with fully atomistic simulations. However, a requirement for its success is a priori knowledge of all relevant reaction pathways and their rate coefficients. This can be difficult for systems with complex chemistry, such as shock-compressed materials at high temperatures and pressures or phenolic spacecraft heat shields undergoing pyrolysis, which can consist of hundreds of molecular species and thousands of distinct reactions. In this work, we develop a method for first estimating a KMC model composed of elementary reactions and rate coefficients by using large datasets derived from a few molecular dynamics (MD) simulations of shock compressed liquid methane, and then using L1 regularization to reduce the estimated chemical reaction network. We find that the full network of 2613 reactions can be reduced by 89% while incurring approximately 9% error in the dominant species (CH4) population. We find that the degree of sparsity achievable decreases when similar accuracy is required for additional populations of species.

Copyright © Materials Research Society 2016 

Access options

Get access to the full version of this content by using one of the access options below. (Log in options will check for institutional or personal access. Content may require purchase if you do not have access.)



Wasserman, L., All of Statistics: A Concise Course in Statistical Inference. (Springer, New York, 2004) p. 122128.Google Scholar
Petzold, L. and Zhu, W., AIChE J. 45: 869886 (1999).Google Scholar
Hannemann-Tamás, R., Gábor, A., Szederkényi, G., Hangos, K., Computers & Mathematics with Applications 65, 10 (2013).Google Scholar
Sikalo, N., Hasemann, O., Schulz, C., Kempf, A., Wlokas, I., Intl J. of Chem. Kinetics 46, 1 (2014).Google Scholar
Nagy, T. and Turányi, T., Combustion and Flame 156, 2 (2009).Google Scholar
Chenoweth, K., van Duin, A. C. T., Goddard, W. A. III, J. Phys. Chem. A 112, 5 (2008).Google Scholar
Qi, T., Reed, E. J., J. Phys. Chem. A 116, 42 (2012).Google Scholar
Gillespie, D., J. Comput. Phys. 22, 4 (1976).CrossRefGoogle Scholar
Gillespie, D., J. Chem. Phys. 115, 4 (2001).CrossRefGoogle Scholar
Gillespie, D., J. Chem. Phys. 113, 1 (2000).Google Scholar
Tibshirani, R., J. R. Stat. Soc. B 58, 1 (1996).Google Scholar
Osborne, M. R., Presnell, B., Turlach, B. A., J. Comp. Graph. Stat. 9, 2 (2000).Google Scholar
Grant, M. and Boyd, S., CVX: Matlab software for disciplined convex programming, version 2.0 beta., September 2013.Google Scholar
Grant, M. and Boyd, S., in Recent Advances in Learning and Control, edited by Blondel, V., Boyd, S., Kimura, H. (Springer, New York 2008), p. 95110.CrossRefGoogle Scholar