The ReaxFF reactive force-field approach has significantly extended the applicability of reactive molecular dynamics simulations to a wide range of material properties and processes.ReaxFF parameters are commonly trai...The ReaxFF reactive force-field approach has significantly extended the applicability of reactive molecular dynamics simulations to a wide range of material properties and processes.ReaxFF parameters are commonly trained to fit a predefined set of quantummechanical data,but it remains uncertain how accurately the quantities of interest are described when applied to complex chemical reactions.Here,we present a dynamic approach based on multiobjective genetic algorithm for the training of ReaxFF parameters and uncertainty quantification of simulated quantities of interest.ReaxFF parameters are trained by directly fitting reactive molecular dynamics trajectories against quantum molecular dynamics trajectories on the fly,where the Pareto optimal front for the multiple quantities of interest provides an ensemble of ReaxFF models for uncertainty quantification.Our in situ multiobjective genetic algorithm workflow achieves scalability by eliminating the file I/O bottleneck using interprocess communications.The in situ multiobjective genetic algorithm workflow has been applied to high-temperature sulfidation of MoO_(3) by H_(2)S precursor,which is an essential reaction step for chemical vapor deposition synthesis of MoS_(2) layers.Our work suggests a new reactive molecular dynamics simulation approach for far-from-equilibrium chemical processes,which quantitatively reproduces quantum molecular dynamics simulations while providing error bars.展开更多
Hetero-structures made from vertically stacked monolayers of transition metal dichalcogenides hold great potential for optoelectronic and thermoelectric devices.Discovery of the optimal layered material for specific a...Hetero-structures made from vertically stacked monolayers of transition metal dichalcogenides hold great potential for optoelectronic and thermoelectric devices.Discovery of the optimal layered material for specific applications necessitates the estimation of key material properties,such as electronic band structure and thermal transport coefficients.However,screening of material properties via brute force ab initio calculations of the entire material structure space exceeds the limits of current computing resources.Moreover,the functional dependence of material properties on the structures is often complicated,making simplistic statistical procedures for prediction difficult to employ without large amounts of data collection.Here,we present a Gaussian process regression model,which predicts material properties of an input hetero-structure,as well as an active learning model based on Bayesian optimization,which can efficiently discover the optimal hetero-structure using a minimal number of ab initio calculations.The electronic band gap,conduction/valence band dispersions,and thermoelectric performance are used as representative material properties for prediction and optimization.The Materials Project platform is used for electronic structure computation,while the BoltzTraP code is used to compute thermoelectric properties.Bayesian optimization is shown to significantly reduce the computational cost of discovering the optimal structure when compared with finding an optimal structure by building a regression model to predict material properties.The models can be used for predictions with respect to any material property and our software,including data preparation code based on the Python Materials Genomics(PyMatGen)library as well as python-based machine learning code,is available open source.展开更多
Mechanical behavior of 2D materials such as MoS_(2) can be tuned by the ancient art of kirigami.Experiments and atomistic simulations show that 2D materials can be stretched more than 50%by strategic insertion of cuts...Mechanical behavior of 2D materials such as MoS_(2) can be tuned by the ancient art of kirigami.Experiments and atomistic simulations show that 2D materials can be stretched more than 50%by strategic insertion of cuts.However,designing kirigami structures with desired mechanical properties is highly sensitive to the pattern and location of kirigami cuts.We use reinforcement learning(RL)to generate a wide range of highly stretchable MoS_(2) kirigami structures.The RL agent is trained by a small fraction(1.45%)of molecular dynamics simulation data,randomly sampled from a search space of over 4 million candidates for MoS_(2)kirigami structures with 6 cuts.After training,the RL agent not only proposes 6-cut kirigami structures that have stretchability above 45%,but also gains mechanistic insight to propose highly stretchable(above 40%)kirigami structures consisting of 8 and 10 cuts from a search space of billion candidates as zero-shot predictions.展开更多
Predictive materials synthesis is the primary bottleneck in realizing functional and quantum materials.Strategies for synthesis of promising materials are currently identified by time-consuming trial and error and the...Predictive materials synthesis is the primary bottleneck in realizing functional and quantum materials.Strategies for synthesis of promising materials are currently identified by time-consuming trial and error and there are no known predictive schemes to design synthesis parameters for materials.We use offline reinforcement learning(RL)to predict optimal synthesis schedules,i.e.,a time-sequence of reaction conditions like temperatures and concentrations,for the synthesis of semiconducting monolayer MoS2 using chemical vapor deposition.The RL agent,trained on 10,000 computational synthesis simulations,learned threshold temperatures and chemical potentials for onset of chemical reactions and predicted previously unknown synthesis schedules that produce well-sulfidized crystalline,phase-pure MoS2.The model can be extended to multi-task objectives such as predicting profiles for synthesis of complex structures including multi-phase heterostructures and can predict long-time behavior of reacting systems,far beyond the domain of molecular dynamics simulations,making these predictions directly relevant to experimental synthesis.展开更多
The dielectric constant(ϵ)is a critical parameter utilized in the design of polymeric dielectrics for energy storage capacitors,microelectronic devices,and high-voltage insulations.However,agile discovery of polymer d...The dielectric constant(ϵ)is a critical parameter utilized in the design of polymeric dielectrics for energy storage capacitors,microelectronic devices,and high-voltage insulations.However,agile discovery of polymer dielectrics with desirableϵremains a challenge,especially for high-energy,high-temperature applications.To aid accelerated polymer dielectrics discovery,we have developed a machine-learning(ML)-based model to instantly and accurately predict the frequency-dependentϵof polymers with the frequency range spanning 15 orders of magnitude.Our model is trained using a dataset of 1210 experimentally measuredϵvalues at different frequencies,an advanced polymer fingerprinting scheme and the Gaussian process regression algorithm.展开更多
基金This work was supported as part of the Computational Materials Sciences Program funded by the U.S.Department of Energy,Office of Science,Basic Energy Sciences,under Award Number DE-SC0014607.
文摘The ReaxFF reactive force-field approach has significantly extended the applicability of reactive molecular dynamics simulations to a wide range of material properties and processes.ReaxFF parameters are commonly trained to fit a predefined set of quantummechanical data,but it remains uncertain how accurately the quantities of interest are described when applied to complex chemical reactions.Here,we present a dynamic approach based on multiobjective genetic algorithm for the training of ReaxFF parameters and uncertainty quantification of simulated quantities of interest.ReaxFF parameters are trained by directly fitting reactive molecular dynamics trajectories against quantum molecular dynamics trajectories on the fly,where the Pareto optimal front for the multiple quantities of interest provides an ensemble of ReaxFF models for uncertainty quantification.Our in situ multiobjective genetic algorithm workflow achieves scalability by eliminating the file I/O bottleneck using interprocess communications.The in situ multiobjective genetic algorithm workflow has been applied to high-temperature sulfidation of MoO_(3) by H_(2)S precursor,which is an essential reaction step for chemical vapor deposition synthesis of MoS_(2) layers.Our work suggests a new reactive molecular dynamics simulation approach for far-from-equilibrium chemical processes,which quantitatively reproduces quantum molecular dynamics simulations while providing error bars.
基金This work was supported as part of the Computational Materials Sciences Program funded by the US Department of Energy,Office of Science,Basic Energy Sciences,under Award Number DE-SC0014607Calculations were performed at the Center for High Performance Computing of the University of Southern California,as well as the National Energy Research Scientific Computing Center,a DOE Office of Science User Facility supported by the Office of Science of the US Department of Energy under Contract No.DE-AC02-05CH11231.
文摘Hetero-structures made from vertically stacked monolayers of transition metal dichalcogenides hold great potential for optoelectronic and thermoelectric devices.Discovery of the optimal layered material for specific applications necessitates the estimation of key material properties,such as electronic band structure and thermal transport coefficients.However,screening of material properties via brute force ab initio calculations of the entire material structure space exceeds the limits of current computing resources.Moreover,the functional dependence of material properties on the structures is often complicated,making simplistic statistical procedures for prediction difficult to employ without large amounts of data collection.Here,we present a Gaussian process regression model,which predicts material properties of an input hetero-structure,as well as an active learning model based on Bayesian optimization,which can efficiently discover the optimal hetero-structure using a minimal number of ab initio calculations.The electronic band gap,conduction/valence band dispersions,and thermoelectric performance are used as representative material properties for prediction and optimization.The Materials Project platform is used for electronic structure computation,while the BoltzTraP code is used to compute thermoelectric properties.Bayesian optimization is shown to significantly reduce the computational cost of discovering the optimal structure when compared with finding an optimal structure by building a regression model to predict material properties.The models can be used for predictions with respect to any material property and our software,including data preparation code based on the Python Materials Genomics(PyMatGen)library as well as python-based machine learning code,is available open source.
基金This work was supported by National Science Foundation,Future Manufacturing Program,Award 2036359This research was partly supported by Aurora Early Science programs and used resources of the Argonne Leadership Computing Facility,which is a DOE Office of Science User Facility supported under Contract DE-AC02-06CH11357Computations were performed at the Argonne Leadership Computing Facility under the DOE INCITE and Aurora Early Science programs and at the Center for Advanced Research Computing of the University of Southern California.We would like to thank one of the reviewers for asking us to examine if RL can be used to propose a highstretchability kirigami structure with 10 cuts,which led to zero-shot predictions for 8-and 10-cut structures that have stretchability above 40%.
文摘Mechanical behavior of 2D materials such as MoS_(2) can be tuned by the ancient art of kirigami.Experiments and atomistic simulations show that 2D materials can be stretched more than 50%by strategic insertion of cuts.However,designing kirigami structures with desired mechanical properties is highly sensitive to the pattern and location of kirigami cuts.We use reinforcement learning(RL)to generate a wide range of highly stretchable MoS_(2) kirigami structures.The RL agent is trained by a small fraction(1.45%)of molecular dynamics simulation data,randomly sampled from a search space of over 4 million candidates for MoS_(2)kirigami structures with 6 cuts.After training,the RL agent not only proposes 6-cut kirigami structures that have stretchability above 45%,but also gains mechanistic insight to propose highly stretchable(above 40%)kirigami structures consisting of 8 and 10 cuts from a search space of billion candidates as zero-shot predictions.
基金This work was supported as part of the Computational Materials Sciences Program funded by the U.S.Department of Energy,Office of Science,Basic Energy Sciences,under Award Number DE-SC0014607This research was partly supported by Aurora Early Science programs and used resources of the Argonne Leadership Computing Facility,which is a DOE Office of Science User Facility supported under Contract DEAC02-06CH11357Computations were performed at the Argonne Leadership Computing Facility under the DOE INCITE and Aurora Early Science programs and at the Center for Advanced Research Computing of the University of Southern California.
文摘Predictive materials synthesis is the primary bottleneck in realizing functional and quantum materials.Strategies for synthesis of promising materials are currently identified by time-consuming trial and error and there are no known predictive schemes to design synthesis parameters for materials.We use offline reinforcement learning(RL)to predict optimal synthesis schedules,i.e.,a time-sequence of reaction conditions like temperatures and concentrations,for the synthesis of semiconducting monolayer MoS2 using chemical vapor deposition.The RL agent,trained on 10,000 computational synthesis simulations,learned threshold temperatures and chemical potentials for onset of chemical reactions and predicted previously unknown synthesis schedules that produce well-sulfidized crystalline,phase-pure MoS2.The model can be extended to multi-task objectives such as predicting profiles for synthesis of complex structures including multi-phase heterostructures and can predict long-time behavior of reacting systems,far beyond the domain of molecular dynamics simulations,making these predictions directly relevant to experimental synthesis.
基金This work is supported by the Office of Naval Research through N0014-17-1-2656,a Multi-University Research Initiative(MURI)grant.
文摘The dielectric constant(ϵ)is a critical parameter utilized in the design of polymeric dielectrics for energy storage capacitors,microelectronic devices,and high-voltage insulations.However,agile discovery of polymer dielectrics with desirableϵremains a challenge,especially for high-energy,high-temperature applications.To aid accelerated polymer dielectrics discovery,we have developed a machine-learning(ML)-based model to instantly and accurately predict the frequency-dependentϵof polymers with the frequency range spanning 15 orders of magnitude.Our model is trained using a dataset of 1210 experimentally measuredϵvalues at different frequencies,an advanced polymer fingerprinting scheme and the Gaussian process regression algorithm.