-
Molecular-Resolution Imaging of Ice Crystallized from Liquid Water
Authors:
Jingshan S. Du,
Suvo Banik,
Henry Chan,
Birk Fritsch,
Ying Xia,
Andreas Hutzler,
Subramanian K. R. S. Sankaranarayanan,
James J. De Yoreo
Abstract:
Despite the ubiquity of ice, a molecular-resolution image of ice crystallized from liquid water or the resulting defect structure has never been obtained. Here, we report the stabilization and angstrom-resolution electron imaging of ice Ih crystallized from liquid water. We combine lattice mapping with molecular dynamics simulations to reveal that ice formation is highly tolerant to nanoscale defe…
▽ More
Despite the ubiquity of ice, a molecular-resolution image of ice crystallized from liquid water or the resulting defect structure has never been obtained. Here, we report the stabilization and angstrom-resolution electron imaging of ice Ih crystallized from liquid water. We combine lattice mapping with molecular dynamics simulations to reveal that ice formation is highly tolerant to nanoscale defects such as misoriented subdomains and trapped gas bubbles, which are stabilized by molecular-scale structural motifs. Importantly, bubble surfaces adopt low-energy nanofacets and create negligible strain fields in the surrounding crystal. These bubbles can dynamically nucleate, grow, migrate, dissolve, and coalesce under electron irradiation and be monitored in situ near a steady state. This work opens the door to understanding water crystallization behaviors at an unprecedented spatial resolution.
△ Less
Submitted 4 June, 2024; v1 submitted 2 June, 2024;
originally announced June 2024.
-
Lorentz violating backgrounds from quadratic, shift-symmetric, ultralight dark matter
Authors:
Minyuan Jiang,
Benjamin D. Pecjak,
Gilad Perez,
Somasundaram Sankaranarayanan
Abstract:
We consider an effective theory for a shift-symmetric, quadratically-coupled, ultralight spin-0 field. The leading CP-conserving interactions with Standard Model fields in the effective theory arise at dimension 8. We discuss the renormalization group evolution and positivity bounds on these operators, as well as their possible UV origins. Assuming that the spin-0 field is associated with an ultra…
▽ More
We consider an effective theory for a shift-symmetric, quadratically-coupled, ultralight spin-0 field. The leading CP-conserving interactions with Standard Model fields in the effective theory arise at dimension 8. We discuss the renormalization group evolution and positivity bounds on these operators, as well as their possible UV origins. Assuming that the spin-0 field is associated with an ultralight dark matter candidate, we discuss the effects of the dimension-8 operators on experiments searching for the oscillation of fundamental constants and Lorentz violation. We find that the direct bounds on these two effects are of similar strength but rather weak, corresponding to a UV cutoff scale of keV order, as they are mediated by dimension-8 operators.
△ Less
Submitted 15 May, 2024; v1 submitted 26 April, 2024;
originally announced April 2024.
-
Evolutionary Search and Theoretical Study of Silicene Grain Boundaries' Mechanical Properties
Authors:
Jianan Zhang,
Aditya Koneru,
Subramanian K. R. S. Sankaranarayanan,
Carmen M. Lilley
Abstract:
Defects such as grain boundaries (GBs) are almost inevitable during the synthesis process of 2D materials. To take advantage of the fascinating properties of 2D materials, understanding the nature and impact of various GB structures on the pristine 2D sheet is crucial. In this work, using an evolutionary algorithm search, we predict a wide variety of silicene GB structures with very different atom…
▽ More
Defects such as grain boundaries (GBs) are almost inevitable during the synthesis process of 2D materials. To take advantage of the fascinating properties of 2D materials, understanding the nature and impact of various GB structures on the pristine 2D sheet is crucial. In this work, using an evolutionary algorithm search, we predict a wide variety of silicene GB structures with very different atomic structures compared to those found in graphene or hexagonal boron-nitride. Twenty-one GBs with the lowest energy were validated by density functional theory (DFT) - a majority of which were previously unreported to our best knowledge. Based on the diversity of the GB predictions, we found that the formation energy and mechanical properties can be dramatically altered by adatoms positions within a GB and certain types of atomic structures, such as four-atom rings. To study the mechanical behavior of these GBs, we apply strain to the GB structures stepwise and use DFT calculations to investigate the mechanical properties of 9 representative structures. It is observed that GB structures based on pentagon-heptagon pairs are likely to have similar or higher in-plane stiffness and strength compared with the zigzag orientation of pristine silicene. However, an adatom located at the hollow site of a heptagon ring can significantly deteriorate the mechanical strength. For all the structures, the in-plane stiffness and strength were found to decrease with increasing formation energy. For the failure behavior of GB structures, it was found that GB structures based on pentagon-heptagon pairs have failure behavior similar to graphene. We also found that the GB structures with atoms positioned outside of the 2D plane tend to experience phase transitions before failure. Utilizing the evolutionary algorithm, we locate diverse silicene GBs and obtain useful information for their mechanical properties.
△ Less
Submitted 19 February, 2024; v1 submitted 29 September, 2023;
originally announced September 2023.
-
Infrared Nanoimaging of Hydrogenated Perovskite Nickelate Synaptic Devices
Authors:
Sampath Gamage,
Sukriti Manna,
Marc Zajac,
Steven Hancock,
Qi Wang,
Sarabpreet Singh,
Mahdi Ghafariasl,
Kun Yao,
Tom Tiwald,
Tae Joon Park,
David P. Landau,
Haidan Wen,
Subramanian Sankaranarayanan,
Pierre Darancet,
Shriram Ramanathan,
Yohannes Abate
Abstract:
Solid-state devices made from correlated oxides such as perovskite nickelates are promising for neuromorphic computing by mimicking biological synaptic function. However, comprehending dopant action at the nanoscale poses a formidable challenge to understanding the elementary mechanisms involved. Here, we perform operando infrared nanoimaging of hydrogen-doped correlated perovskite, neodymium nick…
▽ More
Solid-state devices made from correlated oxides such as perovskite nickelates are promising for neuromorphic computing by mimicking biological synaptic function. However, comprehending dopant action at the nanoscale poses a formidable challenge to understanding the elementary mechanisms involved. Here, we perform operando infrared nanoimaging of hydrogen-doped correlated perovskite, neodymium nickel oxide (H-NdNiO3) devices and reveal how an applied field perturbs dopant distribution at the nanoscale. This perturbation leads to stripe phases of varying conductivity perpendicular to the applied field, which define the macroscale electrical characteristics of the devices. Hyperspectral nano-FTIR imaging in conjunction with density functional theory calculations unveil a real-space map of multiple vibrational states of H-NNO associated with OH stretching modes and their dependence on the dopant concentration. Moreover, the localization of excess charges induces an out-of-plane lattice expansion in NNO which was confirmed by in-situ - x-ray diffraction and creates a strain that acts as a barrier against further diffusion. Our results and the techniques presented here hold great potential to the rapidly growing field of memristors and neuromorphic devices wherein nanoscale ion motion is fundamentally responsible for function.
△ Less
Submitted 29 August, 2023;
originally announced September 2023.
-
Accelerated Design of Block Copolymers: An Unbiased Exploration Strategy via Fusion of Molecular Dynamics Simulations and Machine Learning
Authors:
Jan Michael Y. Carrillo,
Vijith P,
Tarak K. Patra,
Zhan Chen,
Thomas P. Russell,
Subramanian KRS Sankaranarayanan,
Bobby G. Sumpter,
Rohit Batra
Abstract:
Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatbilization, chemical transformations and separations. s-BCPs are star-shaped macromolecules comprised of linear chains of different chemical blocks (e.g., solvophilic and solvophobic blocks) that are covalently joined at one junction point. Various parameters of these macromolec…
▽ More
Star block copolymers (s-BCPs) have potential applications as novel surfactants or amphiphiles for emulsification, compatbilization, chemical transformations and separations. s-BCPs are star-shaped macromolecules comprised of linear chains of different chemical blocks (e.g., solvophilic and solvophobic blocks) that are covalently joined at one junction point. Various parameters of these macromolecules can be tuned to obtain desired surface properties, including the number of arms, composition of the arms, and the degree-of-polymerization of the blocks (or the length of the arm). This makes identification of the optimal s-BCP design highly non-trivial as the total number of plausible s-BCPs architectures is experimentally or computationally intractable. In this work, we use molecular dynamics (MD) simulations coupled with reinforcement learning based Monte Carlo tree search (MCTS) to identify s-BCPs designs that minimize the interfacial tension between polar and non-polar solvents. We first validate the MCTS approach for design of small- and medium-sized s-BCPs, and then use it to efficiently identify sequences of copolymer blocks for large-sized s-BCPs. The structural origins of interfacial tension in these systems are also identified using the configurations obtained from MD simulations. Chemical insights on the arrangement of copolymer blocks that promote lower interfacial tension were mined using machine learning (ML) techniques. Overall, this work provides an efficient approach to solve design problems via fusion of simulations and ML and provide important groundwork for future experimental investigation of s-BCPs sequences for various applications.
△ Less
Submitted 16 August, 2023;
originally announced August 2023.
-
AutoPhaseNN: Unsupervised Physics-aware Deep Learning of 3D Nanoscale Bragg Coherent Diffraction Imaging
Authors:
Yudong Yao,
Henry Chan,
Subramanian Sankaranarayanan,
Prasanna Balaprakash,
Ross J. Harder,
Mathew J. Cherukara
Abstract:
The problem of phase retrieval, or the algorithmic recovery of lost phase information from measured intensity alone, underlies various imaging methods from astronomy to nanoscale imaging. Traditional methods of phase retrieval are iterative in nature, and are therefore computationally expensive and time consuming. More recently, deep learning (DL) models have been developed to either provide learn…
▽ More
The problem of phase retrieval, or the algorithmic recovery of lost phase information from measured intensity alone, underlies various imaging methods from astronomy to nanoscale imaging. Traditional methods of phase retrieval are iterative in nature, and are therefore computationally expensive and time consuming. More recently, deep learning (DL) models have been developed to either provide learned priors to iterative phase retrieval or in some cases completely replace phase retrieval with networks that learn to recover the lost phase information from measured intensity alone. However, such models require vast amounts of labeled data, which can only be obtained through simulation or performing computationally prohibitive phase retrieval on hundreds of or even thousands of experimental datasets. Using a 3D nanoscale X-ray imaging modality (Bragg Coherent Diffraction Imaging or BCDI) as a representative technique, we demonstrate AutoPhaseNN, a DL-based approach which learns to solve the phase problem without labeled data. By incorporating the physics of the imaging technique into the DL model during training, AutoPhaseNN learns to invert 3D BCDI data from reciprocal space to real space in a single shot without ever being shown real space images. Once trained, AutoPhaseNN is about one hundred times faster than traditional iterative phase retrieval methods while providing comparable image quality.
△ Less
Submitted 4 April, 2022; v1 submitted 28 September, 2021;
originally announced September 2021.
-
Real-time 3D Nanoscale Coherent Imaging via Physics-aware Deep Learning
Authors:
Henry Chan,
Youssef S. G. Nashed,
Saugat Kandel,
Stephan Hruszkewycz,
Subramanian Sankaranarayanan,
Ross J. Harder,
Mathew J. Cherukara
Abstract:
Phase retrieval, the problem of recovering lost phase information from measured intensity alone, is an inverse problem that is widely faced in various imaging modalities ranging from astronomy to nanoscale imaging. The current process of phase recovery is iterative in nature. As a result, the image formation is time-consuming and computationally expensive, precluding real-time imaging. Here, we us…
▽ More
Phase retrieval, the problem of recovering lost phase information from measured intensity alone, is an inverse problem that is widely faced in various imaging modalities ranging from astronomy to nanoscale imaging. The current process of phase recovery is iterative in nature. As a result, the image formation is time-consuming and computationally expensive, precluding real-time imaging. Here, we use 3D nanoscale X-ray imaging as a representative example to develop a deep learning model to address this phase retrieval problem. We introduce 3D-CDI-NN, a deep convolutional neural network and differential programming framework trained to predict 3D structure and strain solely from input 3D X-ray coherent scattering data. Our networks are designed to be "physics-aware" in multiple aspects; in that the physics of x-ray scattering process is explicitly enforced in the training of the network, and the training data are drawn from atomistic simulations that are representative of the physics of the material. We further refine the neural network prediction through a physics-based optimization procedure to enable maximum accuracy at lowest computational cost. 3D-CDI-NN can invert a 3D coherent diffraction pattern to real-space structure and strain hundreds of times faster than traditional iterative phase retrieval methods, with negligible loss in accuracy. Our integrated machine learning and differential programming solution to the phase retrieval problem is broadly applicable across inverse problems in other application areas.
△ Less
Submitted 16 June, 2020;
originally announced June 2020.
-
Active Learning A Neural Network Model For Gold Clusters \& Bulk From Sparse First Principles Training Data
Authors:
Troy D Loeffler,
Sukriti Manna,
Tarak K Patra,
Henry Chan,
Badri Narayanan,
Subramanian Sankaranarayanan
Abstract:
Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Classical MD typically employ predefined functional forms which limits t…
▽ More
Small metal clusters are of fundamental scientific interest and of tremendous significance in catalysis. These nanoscale clusters display diverse geometries and structural motifs depending on the cluster size; a knowledge of this size-dependent structural motifs and their dynamical evolution has been of longstanding interest. Classical MD typically employ predefined functional forms which limits their ability to capture such complex size-dependent structural and dynamical transformation. Neural Network (NN) based potentials represent flexible alternatives and in principle, well-trained NN potentials can provide high level of flexibility, transferability and accuracy on-par with the reference model used for training. A major challenge, however, is that NN models are interpolative and requires large quantities of training data to ensure that the model adequately samples the energy landscape both near and far-from-equilibrium. Here, we introduce an active learning (AL) scheme that trains a NN model on-the-fly with minimal amount of first-principles based training data. Our AL workflow is initiated with a sparse training dataset (1 to 5 data points) and is updated on-the-fly via a Nested Ensemble Monte Carlo scheme that iteratively queries the energy landscape in regions of failure and updates the training pool to improve the network performance. Using a representative system of gold clusters, we demonstrate that our AL workflow can train a NN with ~500 total reference calculations. Our NN predictions are within 30 meV/atom and 40 meV/Åof the reference DFT calculations. Moreover, our AL-NN model also adequately captures the various size-dependent structural and dynamical properties of gold clusters in excellent agreement with DFT calculations and available experiments.
△ Less
Submitted 5 June, 2020;
originally announced June 2020.
-
Screening of Therapeutic Agents for COVID-19 using Machine Learning and Ensemble Docking Simulations
Authors:
Rohit Batra,
Henry Chan,
Ganesh Kamath,
Rampi Ramprasad,
Mathew J. Cherukara,
Subramanian Sankaranarayanan
Abstract:
The world has witnessed unprecedented human and economic loss from the COVID-19 disease, caused by the novel coronavirus SARS-CoV-2. Extensive research is being conducted across the globe to identify therapeutic agents against the SARS-CoV-2. Here, we use a powerful and efficient computational strategy by combining machine learning (ML) based models and high-fidelity ensemble docking simulations t…
▽ More
The world has witnessed unprecedented human and economic loss from the COVID-19 disease, caused by the novel coronavirus SARS-CoV-2. Extensive research is being conducted across the globe to identify therapeutic agents against the SARS-CoV-2. Here, we use a powerful and efficient computational strategy by combining machine learning (ML) based models and high-fidelity ensemble docking simulations to enable rapid screening of possible therapeutic molecules (or ligands). Our screening is based on the binding affinity to either the isolated SARS-CoV-2 S-protein at its host receptor region or to the Sprotein-human ACE2 interface complex, thereby potentially limiting and/or disrupting the host-virus interactions. We first apply our screening strategy to two drug datasets (CureFFI and DrugCentral) to identify hundreds of ligands that bind strongly to the aforementioned two systems. Candidate ligands were then validated by all atom docking simulations. The validated ML models were subsequently used to screen a large bio-molecule dataset (with nearly a million entries) to provide a rank-ordered list of ~19,000 potentially useful compounds for further validation. Overall, this work not only expands our knowledge of small-molecule treatment against COVID-19, but also provides an efficient pathway to perform high-throughput computational drug screening by combining quick ML surrogate models with expensive high-fidelity simulations, for accelerating the therapeutic cure of diseases.
△ Less
Submitted 7 April, 2020;
originally announced April 2020.
-
Machine Learning for Multi-fidelity Scale Bridging and Dynamical Simulations of Materials
Authors:
Rohit Batra,
Subramanian Sankaranarayanan
Abstract:
Molecular dynamics (MD) is a powerful and popular tool for understanding the dynamical evolution of materials at the nano and mesoscopic scales. There are various flavors of MD ranging from the high fidelity albeit computationally expensive ab-initio MD to relatively lower fidelity but much more efficient classical MD such as atomistic and coarse-grained models. Each of these different flavors of…
▽ More
Molecular dynamics (MD) is a powerful and popular tool for understanding the dynamical evolution of materials at the nano and mesoscopic scales. There are various flavors of MD ranging from the high fidelity albeit computationally expensive ab-initio MD to relatively lower fidelity but much more efficient classical MD such as atomistic and coarse-grained models. Each of these different flavors of MD have been independently used by materials scientists to bring about breakthroughs in materials discovery and design. A significant gulf exists between the various MD flavors, each having varying levels of fidelity. The accuracy of DFT or ab-initio MD is generally much higher than that of classical atomistic simulations which is higher than that of coarse-grained models. Multi-fidelity scale bridging to combine the accuracy and flexibility of ab-initio MD with efficiency classical MD has been a longstanding goal. The advent of big-data analytics has brought to the forefront powerful machine learning methods that can be deployed to achieve this goal. Here, we provide our perspective on the challenges in multi-fidelity scale bridging and trace the developments leading up to the use of machine learning algorithms and data-science towards addressing this grand challenge.
△ Less
Submitted 1 April, 2020;
originally announced April 2020.
-
A coarse-grained deep neural network model for liquid water
Authors:
Tarak K Patra,
Troy D. Loeffler,
Henry Chan,
Mathew J. Cherukara,
Badri Narayanan,
Subramanian K. R. S. Sankaranarayanan
Abstract:
We introduce a coarse-grained deep neural network model (CG-DNN) for liquid water that utilizes 50 rotational and translational invariant coordinates, and is trained exclusively against energies of ~30,000 bulk water configurations. Our CG-DNN potential accurately predicts both the energies and molecular forces of water; within 0.9 meV/molecule and 54 meV/angstrom of a reference (coarse-grained bo…
▽ More
We introduce a coarse-grained deep neural network model (CG-DNN) for liquid water that utilizes 50 rotational and translational invariant coordinates, and is trained exclusively against energies of ~30,000 bulk water configurations. Our CG-DNN potential accurately predicts both the energies and molecular forces of water; within 0.9 meV/molecule and 54 meV/angstrom of a reference (coarse-grained bond-order potential) model. The CG-DNN water model also provides good prediction of several structural, thermodynamic, and temperature dependent properties of liquid water, with values close to that obtained from the reference model. More importantly, CG-DNN captures the well-known density anomaly of liquid water observed in experiments. Our work lays the groundwork for a scheme where existing empirical water models can be utilized to develop fully flexible neural network framework that can subsequently be trained against sparse data from high-fidelity albeit expensive beyond-DFT calculations.
△ Less
Submitted 14 October, 2019; v1 submitted 1 October, 2019;
originally announced October 2019.
-
Active Learning the Coarse-Grained Energy Landscape For Water Clusters From Sparse Training Data
Authors:
Troy D. Loeffler,
Tarak K. Patra,
Henry Chan,
Mathew Cherukara,
Subramanian K. R. S. Sankaranarayanan
Abstract:
ANNs are currently trained by generating large quantities (On the order of $10^{4}$ or greater) of structural data in hopes that the ANN has adequately sampled the energy landscape both near and far-from-equilibrium. This can, however, be a bit prohibitive when it comes to more accurate levels of quantum theory. As such it is desirable to train a model using the absolute minimal data set possible,…
▽ More
ANNs are currently trained by generating large quantities (On the order of $10^{4}$ or greater) of structural data in hopes that the ANN has adequately sampled the energy landscape both near and far-from-equilibrium. This can, however, be a bit prohibitive when it comes to more accurate levels of quantum theory. As such it is desirable to train a model using the absolute minimal data set possible, especially when costs of high-fidelity calculations such as CCSD and QMC are high. Here, we present an Active Learning approach that iteratively trains an ANN model to faithfully replicate the coarse-grained energy surface of water clusters using only 426 total structures in its training data. Our active learning workflow starts with a sparse training dataset which is continually updated via a Monte Carlo scheme that sparsely queries the energy landscape and tests the network performance. Next, the network is retrained with an updated training set that includes failed configurations/energies from previous iteration until convergence is attained. Once trained, we generate an extensive test set of ~100,000 configurations sampled across clusters ranging from 1 to 200 molecules and demonstrate that the trained network adequately reproduces the energies (within mean absolute error (MAE) of ~ 2 meV/molecule) and forces (MAE ~ 40 meV/{Å}) compared to the reference model. More importantly, the trained ANN model also accurately captures both the structure as well as the free energy as a function of the various cluster sizes. Overall, this study reports a new active learning scheme with promising strategy to develop accurate force-fields for molecular simulations using extremely sparse training data sets.
△ Less
Submitted 2 November, 2019; v1 submitted 25 September, 2019;
originally announced September 2019.
-
Graphene Overcoats for Ultra-High Storage Density Magnetic Media
Authors:
N. Dwivedi,
A. K. Ott,
K. Sasikumar,
C. Dou,
R. J. Yeo,
B. Narayanan,
U. Sassi,
D. De Fazio,
G. Soavi,
T. Dutta,
S. K. R. S. Sankaranarayanan,
A. C. Ferrari,
C. S. Bhatia
Abstract:
Hard disk drives (HDDs) are used as secondary storage in a number of digital electronic devices owing to low cost ($<$0.1\$/GB at 2016 prices) and large data storage capacity (10TB with a 3.5 inch HDD). Due to the exponentially increasing amount of data, there is a need to increase areal storage densities beyond$\sim$1Tb/in$^2$. This requires the thickness of carbon overcoats (COCs) to be$<…
▽ More
Hard disk drives (HDDs) are used as secondary storage in a number of digital electronic devices owing to low cost ($<$0.1\$/GB at 2016 prices) and large data storage capacity (10TB with a 3.5 inch HDD). Due to the exponentially increasing amount of data, there is a need to increase areal storage densities beyond$\sim$1Tb/in$^2$. This requires the thickness of carbon overcoats (COCs) to be$<$2nm. Friction, wear, corrosion, and thermal stability are critical concerns$<$2nm, where most of the protective properties of current COCs are lost. This limits current technology and restricts COC integration with heat assisted magnetic recording technology (HAMR), since this also requires laser irradiation stability. Here we show that graphene-based overcoats can overcome all these limitations. 2-4 layers of graphene enable two-fold reduction in friction and provide better corrosion and wear than state-of-the-art COCs. A single graphene layer is enough to reduce corrosion$\sim$2.5 times. We also show that graphene can withstand HAMR conditions. Thus, graphene-based overcoats can enable ultrahigh areal density HDDs$>$10Tb/in$^2$.
△ Less
Submitted 2 June, 2019;
originally announced June 2019.
-
Phase Segmentation in Atom-Probe Tomography Using Deep Learning-Based Edge Detection
Authors:
Sandeep Madireddy,
Ding-Wen Chung,
Troy Loeffler,
Subramanian K. R. S. Sankaranarayanan,
David N. Seidman,
Prasanna Balaprakash,
Olle Heinonen
Abstract:
Atom-probe tomography (APT) facilitates nano- and atomic-scale characterization and analysis of microstructural features. Specifically, APT is well suited to study the interfacial properties of granular or heterophase systems. Traditionally, the identification of the interface between, for precipitate and matrix phases, in APT data has been obtained either by extracting iso-concentration surfaces…
▽ More
Atom-probe tomography (APT) facilitates nano- and atomic-scale characterization and analysis of microstructural features. Specifically, APT is well suited to study the interfacial properties of granular or heterophase systems. Traditionally, the identification of the interface between, for precipitate and matrix phases, in APT data has been obtained either by extracting iso-concentration surfaces based on a user-supplied concentration value or by manually perturbing the concentration value until the iso-concentration surface qualitatively matches the interface. These approaches are subjective, not scalable, and may lead to inconsistencies due to local composition inhomogeneities.
We propose a digital image segmentation approach based on deep neural networks that transfer learned knowledge from natural images to automatically segment the data obtained from APT into different phases. This approach not only provides an efficient way to segment the data and extract interfacial properties but does so without the need for expensive interface labeling for training the segmentation model.
We consider here a system with a precipitate phase in a matrix and with three different interface modalities---layered, isolated, and interconnected---that are obtained for different relative geometries of the precipitate phase. We demonstrate the accuracy of our segmentation approach through qualitative visualization of the interfaces, as well as through quantitative comparisons with proximity histograms obtained by using more traditional approaches.
△ Less
Submitted 10 April, 2019;
originally announced April 2019.
-
Comparing optimization strategies for force field parameterization
Authors:
Fatih G. Sen,
Badri Narayanan,
Jeffrey Larson,
Alper Kinaci,
Kiran Sasikumar,
Michael J. Davis,
Stefan M. Wild,
Stephen K. Gray,
Subramanian K. R. S. Sankaranarayanan,
Maria K. Y. Chan
Abstract:
Classical molecular dynamics (MD) simulations enable modeling of materials and examination of microscopic details that are not accessible experimentally. The predictive capability of MD relies on the force field (FF) used to describe interatomic interactions. FF parameters are typically determined to reproduce selected material properties computed from density functional theory (DFT) and/or measur…
▽ More
Classical molecular dynamics (MD) simulations enable modeling of materials and examination of microscopic details that are not accessible experimentally. The predictive capability of MD relies on the force field (FF) used to describe interatomic interactions. FF parameters are typically determined to reproduce selected material properties computed from density functional theory (DFT) and/or measured experimentally. A common practice in parameterizing FFs is to use least-squares local minimization algorithms. Genetic algorithms (GAs) have also been demonstrated as a viable global optimization approach, even for complex FFs. However, an understanding of the relative effectiveness and efficiency of different optimization techniques for the determination of FF parameters is still lacking. In this work, we evaluate various FF parameter optimization schemes, using as example a training data set calculated from DFT for different polymorphs of Ir$O_2$. The Morse functional form is chosen for the pairwise interactions and the optimization of the parameters against the training data is carried out using (1) multi-start local optimization algorithms: Simplex, Levenberg-Marquardt, and POUNDERS, (2) single-objective GA, and (3) multi-objective GA. Using random search as a baseline, we compare the algorithms in terms of reaching the lowest error, and number of function evaluations. We also compare the effectiveness of different approaches for FF parameterization using a test data set with known ground truth (i.e generated from a specific Morse FF). We find that the performance of optimization approaches differs when using the Test data vs. the DFT data. Overall, this study provides insight for selecting a suitable optimization method for FF parameterization, which in turn can enable more accurate prediction of material properties and chemical phenomena.
△ Less
Submitted 1 December, 2018;
originally announced December 2018.
-
Methodology of Parameterization of Molecular Mechanics Force Field From Quantum Chemistry Calculations using Genetic Algorithm: A case study of methanol
Authors:
Ying Li,
Hui Li,
Maria K. Y. Chan,
Subramanian Sankaranarayanan,
Benoît Rouxb
Abstract:
In molecular dynamics (MD) simulation, force field determines the capability of an individual model in capturing physical and chemistry properties. The method for generating proper parameters of the force field form is the key component for computational research in chemistry, biochemistry, and condensed-phase physics. Our study showed that the feasibility to predict experimental condensed phase p…
▽ More
In molecular dynamics (MD) simulation, force field determines the capability of an individual model in capturing physical and chemistry properties. The method for generating proper parameters of the force field form is the key component for computational research in chemistry, biochemistry, and condensed-phase physics. Our study showed that the feasibility to predict experimental condensed phase properties (i.e., density and heat of vaporization) of methanol through problem specific force field from only quantum chemistry information. To acquire the satisfying parameter sets of the force field, the genetic algorithm (GA) is the main optimization method. For electrostatic potential energy, we optimized both the electrostatic parameters of methanol using the GA method, which leads to low deviations of between the quantum mechanics (QM) calculations and the GA optimized parameters. We optimized the van der Waals (vdW) parameters both using GA and guided GA methods by calibrating interaction energy of various methanol homo-clusters, such as nonamers, undecamers, or tridecamers. Excellent agreement between the training dataset from QM calculations (i.e., MP2) and GA optimized parameters can be achieved. However, only the guided GA method, which eliminates the overestimation of interaction energy from MP2 calculations in the optimization process, provides proper vdW parameters for MD simulation to get the condensed phase properties (i.e., density and heat of vaporization) of methanol. Throughout the whole optimization process, the experimental value were not involved in the objective functions, but were only used for the purpose of justifying models (i.e., nonamers, undecamers, or tridecamers) and validating methods (i.e., GA or guided GA). Our method shows the possibility of developing descriptive polarizable force field using only QM calculations.
△ Less
Submitted 15 July, 2016; v1 submitted 14 July, 2016;
originally announced July 2016.