A new topological descriptor for water network structure

Abstract

Bulk water molecular dynamics simulations based on a series of atomistic water potentials (TIP3P, TIP4P/Ew, SPC/E and OPC) are compared using new techniques from the field of topological data analysis. The topological invariants (the different degrees of homology) derived from each simulation frame are used to create a series of persistence diagrams from the atomic positions. These are averaged over the simulation time using the persistence image formalism, before being normalised by their total magnitude (the L1 norm) to ensure a size independent descriptor (L1NPI). We demonstrate that the L1NPI formalism is suitable for the analysis of systems where the number of molecules varies by at least a factor of 10. Using standard machine learning techniques, a basic linear SVM, it is shown that differences in water models are able to be isolated to different degrees of homology. In particular, whereas first degree homology is able to distinguish between all atomistic potentials studied, OPC is the only potential that differs in its second degree homology. The L1 normalised persistence images are then used in the comparison of a series of Stillinger–Weber potential simulations to the atomistic potentials and the effects of changing the strength of three-body interactions on the structures is easily evident in L1NPI space, with a reduction in variance of structures as interaction strength increases being the most obvious result. Furthermore, there is a clear tracking in L1NPI space of the λ parameter. The L1NPI formalism presents a useful new technique for the analysis of water and other materials. It is approximately size-independent, and has been shown to contain information as to real structures in the system. We finally present a perspective on the use of L1NPIs and other persistent homology techniques as a descriptor for water solubility. Electronic supplementary material The online version of this article (10.1186/s13321-019-0369-0) contains supplementary material, which is available to authorized users.

A Topological Data Analysis perspective on noncovalent interactions in relativistic calculations

“…In practice, it allows users to discriminate important features from nonsignificant configurations and provides the necessary basis for multiscale data analysis . Recent years have witnessed some successful applications of TDA to chemistry, where it helped to understand hydrogen‐bonding networks in ion and molecule aggregates, aqueous solubility of molecules, stability of fullerenes, molecular transition pathways, conformational spaces of molecules, or the bonding patterns in molecular systems …”

Topological nature of the liquid–liquid phase transition in tetrahedral liquids

“…The anomalous behaviour of liquid water, along with water polyamorphism, was then suggested to be a manifestation of the presence of this LLPT. Recent numerical studies 8,9 , as well as experiments [10][11][12] (despite the difficulties induced by the rapid formation of ice around the predicted critical temperature and pressure conditions), have provided strong support to this fascinating hypothesis.Several order parameters, based on geometric or energetic criteria, have been proposed to characterize the transition [13][14][15][16][17][18][19][20] . In particular, the local structure index 13 and the r 5 parameter 14 -two order parameters that effectively measure the distance between the first and second coordination shells-have been routinely applied.…”

Topological Analysis of Molecular Dynamics Simulations using the Euler Characteristic

“…20,25 Persistent homology has been shown to provide powerful characterizations of the topology and geometry of data. [25][26][27][28][29][30] However, the outputs of these methods (e.g., persistence diagrams) can be difficult to directly integrate into data analysis and ML tasks without further transformation (e.g., vectorization and smoothing) and hyperparameter optimization. 21 The main goal of this work is to demonstrate that topology can be applied to a broad range of molecular simulation settings.…”