Method that utilizes variety-based matrix completion (VMC) aided by polynomial interpolation to recover projected Hessian eigenvalues constituting the minimum energy path of a reaction. This utilizes MATLAB code for the PVMC implementation and Python scripts for computing ZCT constants.
Main.m- Runs PVMC and HVMC on one of the reaction ground truth matricescone_method.m- Function to find the most optimal columns to samplevmc_step.m- VMC algorithm (from: https://github.com/gregongie/vmc) described in G. Ongie, R. Willett, R. Nowak, L. Balzano. "Algebraic Variety Models for High-Rank Matrix Completion", in ICML 2017. Available online: https://arxiv.org/abs/1703.09631plot_results.m- Provides visualizations of ground truth and PVMC-recovered matrices, as well as errors and gradient term influences for each iteration.pvmc_step.m- PVMC algorithm optimized for recovering column-sampled matrices.
AnalyzeResults.py- Calculates ZCT coefficients and free energy errors of the PVMC recovered matrix in results.mat as generated by Main.mcalcFreeEnergy.py- Supporting functions to compute zero-point energy and vibrational free energy contributionsZCT.py- Supporting functions to compute ZCTspline.py- Supporting functions to compute create spline interpolations for ZCT integrals.
Figures/- Contains saved figures from running AnalyzeResults.pyMatrixMATs/- Contains the ground truth matrices and 's' vectors in MAT filesMatrixPKLs/- Contains the ground truth matrices and 's' vectors in PKL filesResultMATs/- Contains theresult.matfiles of PVMC runs for all reactions at their 24-column minimum sampling densities
- In
Main.m, modifyK(number of MEP points or columns to sample in addition to stationary points) and thesystemto simulate based on the files in MatrixMATs. Other parameters may be adjusted. RunMain.m, which results inresults.mat - To calculate ZCT, run
AnalyzeResults.pyafter the creation ofresults.matin the previous steps.
Stephen Jon Quiton, Jeongmin Chae, Selin Bac, Kareesa Kron, Urbashi Mitra, Shaama Sharada.