Searching for local order parameters to classify water structures at triple points

Molecular cluster analysis using local order parameters selected by machine learning

Accurately extracting local molecular structures is essential for understanding the mechanisms of phase and structural transitions. A promising method to characterize the local molecular structure is defining the value of the local order parameter (LOP) for each particle. This work develops the Molecular Assembly structure Learning package for Identifying Order parameters (MALIO), a machine learning package that can propose an optimal (set of) LOP(s) quickly and automatically for a huge number of LOP species and various methods of selecting neighboring particles for the calculation. We applied this package to distinguish between the nematic and smectic phases of uniaxial liquid crystal molecules, and selected candidate LOPs that could be used to precisely observe the nematic-smectic phase transition. The LOP candidates were used to observe the nucleation and subsequent percolation transition, and the effect of the choice of LOP species and neighboring particles on the statistics of local molecular structures (clusters) was examined. The procedure revealed the time evolution of the number of clusters and the dependence of the percolation curve on the number of neighboring particles for each LOP species. The LOP species with the lowest dependence on the number of neighboring particles was the best-performing LOP species in the MALIO screening strategy. These results not only show that machine learning can powerfully screen a huge number of LOP species and suggest only a few promising candidates, but also indicate that MALIO can select the best LOP species.

Screening toward the Development of Fingerprints of Atomic Environments Using Bond-Orientational Order Parameters

A combination of atomic numbers and bond-orientational order parameters is considered a candidate for a simple representation that involves information on both the atomic species and their positional relation. The 504 candidates are applied as the fingerprint of the molecules stored in QM9, a data set of computed geometric, energetic, electronic, and thermodynamic properties for 133 885 stable small organic molecules made up of carbon, hydrogen, oxygen, nitrogen, and fluorine atoms. To screen the fingerprints, a regression analysis of the atomic charges given by Open Babel was performed by supervised machine learning. The regression results indicate that the 60 fingerprints successfully estimate Open Babel charges. The results of the dipole moments, an example of a property expressed by charge and position, also had a high accuracy in comparison with the values computed from Open Babel charges. Therefore, the screened 60 fingerprints have the potential to precisely describe the chemical and structural information on the atomic environment of molecules.

Mining of Effective Local Order Parameters to Classify Ice Polymorphs

Order parameters make it possible to quantify the degree of structural ordering in a material and thus to apply as the reaction coordinates during the free-energy analysis of phase or structure transitions. Furthermore, order parameters are useful in determining the local structures of molecular groups during transition stages. However, identifying or developing local order parameters (LOPs) that are sensitive for specific materials and phases is a non-trivial task. In this study, the ability of LOPs to classify the solid and liquid structures of water at coexistence or triple points is investigated with the aid of supervised machine learning. The classification accuracy of a total of 179,738,433 combinations of 493 LOPs is automatically and systematically compared for water structures at the ice Ih-Ic-liquid coexistence point and the ice III-V-liquid and ice V-VI-liquid triple points. The optimal sets of two LOPs are found for each point, and sets of three LOPs are suggested for better accuracy.

Multistep nucleation of anisotropic molecules

Phase transition of anisotropic materials is ubiquitously observed in physics, biology, materials science, and engineering. Nevertheless, how anisotropy of constituent molecules affects the phase transition dynamics is still poorly understood. Here we investigate numerically the phase transition of a simple model system composed of anisotropic molecules, and report on our discovery of multistep nucleation of nuclei with layered positional ordering (smectic ordering), from a fluid-like nematic phase with orientational order only (no positional order). A trinity of molecular dynamics simulation, machine learning, and molecular cluster analysis yielding free energy landscapes unambiguously demonstrates the dynamics of multistep nucleation process involving characteristic metastable clusters that precede supercritical smectic nuclei and cannot be accounted for by the classical nucleation theory. Our work suggests that molecules of simple shape can exhibit rich and complex nucleation processes, and our numerical approach will provide deeper understanding of phase transitions and resulting structures in anisotropic materials such as biological systems and functional materials.