[Paper Review] Analysis of molecular dynamics simulation data via statistical distances between covariance matrices

Yusuke Ono, Takumi Sato+2 more|arXiv (Cornell University)|Mar 18, 2026

Material Dynamics and Properties0 citations

TL;DR

The paper develops a framework that analyzes MD trajectories by constructing Toeplitz-structured covariance matrices from time windows, measures dissimilarities between states with Euclidean distance, and reduces to low-dimensional embeddings via PCA, linking embeddings to diffusion and phase differences.

ABSTRACT

Molecular dynamics (MD) simulations are powerful tools for elucidating the macroscopic physical properties of materials from microscopic atomic behaviors. However, the massive, high-dimensional datasets generated by MD simulations pose a significant challenge for analysis, necessitating efficient dimensionality reduction and feature extraction techniques. While existing methods such as principal component analysis and unsupervised learning have been utilized, issues regarding data efficiency and computational cost remain. In this study, we propose a statistical analysis framework focusing on the analysis of the particle data distributions through their covariance matrices, corresponding to the second-order moments of MD trajectory data. Discrepancies between system states are quantified using statistical distances between these covariance matrices. By applying dimensionality reduction to the resulting distance matrix, we extract lower-dimensional features that characterize the systems' dynamics. We validate the proposed method using Lennard-Jones (LJ) particle systems under different temperature conditions, as well as separate bulk systems of ice and liquid water. The results of LJ particles demonstrate an approximately linear correlation between the first principal component obtained through dimensionality reduction of the distance matrix and the diffusion coefficient. This suggests that global physical properties can be effectively inferred from local statistical information, such as covariance matrices, offering a data-efficient alternative for analyzing complex molecular systems. Furthermore, in the case of separate bulk systems of ice and liquid water, the method successfully distinguishes between the two phases, highlighting its potential for characterizing phase transitions and structural differences in molecular systems.

Motivation & Objective

• Motivate the need for data-efficient analysis of massive MD datasets.
• Propose a covariance-matrix based descriptor for MD states derived from time-windowed particle data.
• Quantify state dissimilarities with a statistical distance between covariance matrices.
• Demonstrate the approach on Lennard-Jones systems and bulk ice/liquid water to connect embeddings with physical properties.

Proposed method

• Partition MD time-series into sub-windows of length N to form data matrices.
• Construct 3N x 3N block covariance matrices with nine N x N blocks and Toeplitz structure for R_alpha_beta to ensure symmetry.
• Estimate correlation blocks r_k^{alpha beta} from data and assemble R_m accordingly.
• Compute Euclidean distance between covariance matrices using Frobenius norm as the dissimilarity measure.
• Compute the ensemble mean of covariance matrices across K segments and apply PCA on the resulting distance matrix to obtain a 2D embedding.
• Correlate the leading principal component with macroscopic properties such as the diffusion coefficient.

Figure 1 : Schematic illustration of the proposed method.

Experimental results

Research questions

• RQ1 Can covariance-matrix based statistics capture thermodynamic and transport differences between MD states?
• RQ2 Does a distance-based embedding reflect changes in temperature or phase similar to conventional metrics?
• RQ3 Is it possible to infer macroscopic properties (e.g., diffusion) from local short-time fluctuations?

Key findings

• For Lennard-Jones systems, the first principal component (PC1) of the embedding correlates linearly with the diffusion coefficient across temperatures.
• Distance matrices reveal distinct distributions corresponding to temperature changes, enabling state discrimination.
• In bulk ice vs. liquid water, the method distinguishes phases via the distance histograms, indicating sensitivity to structural/dynamical differences.
• The approach demonstrates data-efficient reconstruction of meaningful physical variation from short-time velocity data (eight steps in the LJ case).
• The framework is extensible to other SPD-matrix distances and higher-order statistics for richer dynamical descriptions.

Figure 2 : Computation for the covariance matrices.