Learning reaction coordinates via cross-entropy minimization: Application to alanine dipeptide

Abstract

We propose a cross-entropy minimization method for finding the reaction coordinate from a large number of collective variables in complex molecular systems. This method is an extension of the likelihood maximization approach describing the committor function with a sigmoid. By design, the reaction coordinate as a function of various collective variables is optimized such that the distribution of the committor $p_{\mathrm{B}}^{\ast }$ values generated from molecular dynamics simulations can be described in a sigmoidal manner. We also introduce the $L_2$-norm regularization used in the machine learning field to prevent overfitting when the number of considered collective variables is large. The current method is applied to study the isomerization of alanine dipeptide in vacuum, where 45 dihedral angles are used as candidate variables. The regularization parameter is determined by cross-validation using training and test datasets. It is demonstrated that the optimal reaction coordinate involves important dihedral angles, which are consistent with the previously reported results. Furthermore, the points with $p_{\mathrm{B}}^{\ast }\sim 0.5$ clearly indicate a separatrix distinguishing reactant and product states on the potential of mean force using the extracted dihedral angles.

8 pages, 5 figures, 1 table for main text, 5 pages for supplementary materials, to appear in J. Chem. Phys. for "Classical Molecular Dynamics (MD) Simulations: Codes, Algorithms, Force Fields, and Applications"

References in corpus (5)

• Escaping free-energy minima
• Machine learning for molecular simulation
• Recent Developments in Methods for Identifying Reaction Coordinates
• Artificial Intelligence Assists Discovery of Reaction Coordinates and Mechanisms from Molecular Dynamics Simulations
• Maximum Likelihood Analysis of Reaction Coordinates during Solidification in Ni

Cited by in corpus (9)

• State Predictive Information Bottleneck
• Explaining reaction coordinates of alanine dipeptide isomerization obtained from deep neural networks using Explainable Artificial Intelligence (XAI)
• Driving and characterizing nucleation of urea and glycine polymorphs in water
• Variational deep learning of equilibrium transition path ensembles
• Unveiling interatomic distances influencing the reaction coordinates in alanine dipeptide isomerization: An explainable deep learning approach
• Investigating the hyperparameter space of deep neural network models for reaction coordinates
• Deep learning of committor for ion dissociation and interpretable analysis of solvent effects using atom-centered symmetry functions
• Understanding Reaction Mechanisms from Start to Finish
• Deep learning of committor and explainable artificial intelligence analysis for identifying reaction coordinates