1
|
Koda SI, Saito S. Flat-Bottom Elastic Network Model for Generating Improved Plausible Reaction Paths. J Chem Theory Comput 2024. [PMID: 39150850 DOI: 10.1021/acs.jctc.4c00792] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/18/2024]
Abstract
Rapid generation of a plausible reaction path connecting a given reactant and product in advance is crucial for the efficient computation of precise reaction paths or transition states. We propose a computationally efficient potential energy based on the molecular structure to generate such paths. This potential energy has a flat bottom consisting of structures without atomic collisions while preserving nonreactive chemical bonds, bond angles, and partial planar structures. By combining this potential energy with the direct MaxFlux method, a recently developed reaction-path/transition-state search method, we can find the shortest plausible path passing within the bottom. Numerical results show that this combination yields lower energy paths compared to the paths obtained by the well-known image-dependent pair potential. We also theoretically investigate the differences between these two potential energies. The proposed potential energy and path generation routine are implemented in our Python version of the direct MaxFlux method, available on GitHub.
Collapse
Affiliation(s)
- Shin-Ichi Koda
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
- School of Physical Sciences, The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| | - Shinji Saito
- Department of Theoretical and Computational Molecular Science, Institute for Molecular Science, Myodaiji, Okazaki, Aichi 444-8585, Japan
- School of Physical Sciences, The Graduate University for Advanced Studies, Myodaiji, Okazaki, Aichi 444-8585, Japan
| |
Collapse
|
2
|
Teng C, Wang Y, Bao JL. Physical Prior Mean Function-Driven Gaussian Processes Search for Minimum-Energy Reaction Paths with a Climbing-Image Nudged Elastic Band: A General Method for Gas-Phase, Interfacial, and Bulk-Phase Reactions. J Chem Theory Comput 2024; 20:4308-4324. [PMID: 38720441 DOI: 10.1021/acs.jctc.4c00291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
The climbing-image nudged elastic band (CI-NEB) method serves as an indispensable tool for computational chemists, offering insight into minimum-energy reaction paths (MEPs) by delineating both transition states (TSs) and intermediate nonstationary structures along reaction coordinates. However, executing CI-NEB calculations for reactions with extensive reaction coordinate spans necessitates a large number of images to ensure a reliable convergence of the MEPs and TS structures, presenting a computationally demanding optimization challenge, even with mildly costly electronic-structure methods. In this study, we advocate for the utilization of physically inspired prior mean function-based Gaussian processes (GPs) to expedite MEP exploration and TS optimization via the CI-NEB method. By incorporating reliable prior physical approximations into potential energy surface (PES) modeling, we demonstrate enhanced efficiency in multidimensional CI-NEB optimization with surrogate-based optimizers. Our physically informed GP approach not only outperforms traditional nonsurrogate-based optimizers in optimization efficiency but also on-the-fly learns the reaction path valley during optimization, culminating in significant advancements. The surrogate PES derived from our optimization exhibits high accuracy compared to true PES references, aligning with our emphasis on leveraging reliable physical priors for robust and efficient posterior mean learning in GPs. Through a systematic benchmark study encompassing various reaction pathways, including gas-phase, bulk-phase, and interfacial/surface reactions, our physical GPs consistently demonstrate superior efficiency and reliability. For instance, they outperform the popular fast inertial relaxation engine optimizer by approximately a factor of 10, showcasing their versatility and efficacy in exploring reaction mechanisms and surface reaction PESs.
Collapse
Affiliation(s)
- Chong Teng
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Yang Wang
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, United States
| |
Collapse
|
3
|
Iwai S, Hayashi Y, Baba T, Kitagawa Y. Acceleration of hydrolytic ring opening of N7-alkylguanine by the terminal carbamoyl group of glycidamide. Chem Commun (Camb) 2024; 60:5014-5017. [PMID: 38577847 DOI: 10.1039/d3cc04997c] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/06/2024]
Abstract
Hydrolytic ring opening of guanine N7-adducts with compounds containing an oxacyclopropane ring, namely glycidamide, glycidol and 1,2-epoxybutane, was analyzed, and the reaction of the glycidamide adduct was the fastest. The differences in the reaction rates were confirmed by theoretical calculations.
Collapse
Affiliation(s)
- Shigenori Iwai
- Division of Chemistry, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan.
| | - Yuta Hayashi
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| | - Tomohiro Baba
- Division of Chemistry, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan.
| | - Yasutaka Kitagawa
- Division of Chemical Engineering, Graduate School of Engineering Science, Osaka University, 1-3 Machikaneyama, Toyonaka, Osaka 560-8531, Japan
| |
Collapse
|
4
|
Bhuiyan FH, Li YS, Kim SH, Martini A. Shear-activation of mechanochemical reactions through molecular deformation. Sci Rep 2024; 14:2992. [PMID: 38316829 PMCID: PMC10844542 DOI: 10.1038/s41598-024-53254-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2023] [Accepted: 01/30/2024] [Indexed: 02/07/2024] Open
Abstract
Mechanical stress can directly activate chemical reactions by reducing the reaction energy barrier. A possible mechanism of such mechanochemical activation is structural deformation of the reactant species. However, the effect of deformation on the reaction energetics is unclear, especially, for shear stress-driven reactions. Here, we investigated shear stress-driven oligomerization reactions of cyclohexene on silica using a combination of reactive molecular dynamics simulations and ball-on-flat tribometer experiments. Both simulations and experiments captured an exponential increase in reaction yield with shear stress. Elemental analysis of ball-on-flat reaction products revealed the presence of oxygen in the polymers, a trend corroborated by the simulations, highlighting the critical role of surface oxygen atoms in oligomerization reactions. Structural analysis of the reacting molecules in simulations indicated the reactants were deformed just before a reaction occurred. Quantitative evidence of shear-induced deformation was established by comparing bond lengths in cyclohexene molecules in equilibrium and prior to reactions. Nudged elastic band calculations showed that the deformation had a small effect on the transition state energy but notably increased the reactant state energy, ultimately leading to a reduction in the energy barrier. Finally, a quantitative relationship was developed between molecular deformation and energy barrier reduction by mechanical stress.
Collapse
Affiliation(s)
- Fakhrul H Bhuiyan
- Department of Mechanical Engineering, University of California Merced, 5200 N. Lake Road, Merced, CA, 95343, USA
| | - Yu-Sheng Li
- Department of Chemical Engineering and Materials Research Institute, Pennsylvania State University, University Park, PA, 16802, USA
| | - Seong H Kim
- Department of Chemical Engineering and Materials Research Institute, Pennsylvania State University, University Park, PA, 16802, USA
| | - Ashlie Martini
- Department of Mechanical Engineering, University of California Merced, 5200 N. Lake Road, Merced, CA, 95343, USA.
| |
Collapse
|
5
|
Mitsuta Y, Asada T. Curvature-weighted nudged elastic band method using the Riemann curvature. J Comput Chem 2023; 44:662-669. [PMID: 36380703 DOI: 10.1002/jcc.27030] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2022] [Revised: 10/03/2022] [Accepted: 10/08/2022] [Indexed: 11/17/2022]
Abstract
The nudged elastic band (NEB) method is utilized to find reaction paths (RPs) using discretized intermediate structures called "images" between a reactant and a product. In fact, NEB calculations do not always converge because of the bent of RPs. Mathematically, more images are needed for complex curves, and here, we focused on the curvature. In this study, we propose a new method for calculating the curvature of the RPs as well as a method for weighting the spring constant of the NEB with the curvature, which we named the curvature weighted NEB (CW-NEB) method. In addition, we will propose the CW-NEB method with the climbing image (CI) method (CW-CI-NEB). To show the efficiency of our method, calculations for the ene-reaction of the CW-CI-NEB method were performed and compared with those of the CI-NEB methods. The CW-CI-NEB methods converged in shorter iteration times than the CI-NEB calculation, which was found by gathering the bent of RPs.
Collapse
Affiliation(s)
- Yuki Mitsuta
- Graduate School of Science, Osaka Metropolitan University, Osaka, Japan.,RIMED, Osaka Prefecture University, Osaka, Japan
| | - Toshio Asada
- Graduate School of Science, Osaka Metropolitan University, Osaka, Japan.,RIMED, Osaka Prefecture University, Osaka, Japan
| |
Collapse
|
6
|
Li R, Barel N, Subramaniyan V, Tibika F, Hoffman R, Tulchinsky Y. Sulfonium-Pincer Ligands Flexibility in Pt(II) Complexes. Organometallics 2023. [DOI: 10.1021/acs.organomet.2c00578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/05/2023]
Affiliation(s)
- Ruiping Li
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem9190401, Israel
| | - Nitsan Barel
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem9190401, Israel
| | | | - Françoise Tibika
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem9190401, Israel
| | - Roy Hoffman
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem9190401, Israel
| | - Yuri Tulchinsky
- Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem9190401, Israel
| |
Collapse
|
7
|
Holst LH, Madsen NG, Toftgård FT, Rønne F, Moise IM, Petersen EI, Fojan P. De novo design of a polycarbonate hydrolase. Protein Eng Des Sel 2023; 36:gzad022. [PMID: 38035789 DOI: 10.1093/protein/gzad022] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Revised: 10/31/2023] [Accepted: 11/14/2023] [Indexed: 12/02/2023] Open
Abstract
Enzymatic degradation of plastics is currently limited to the use of engineered natural enzymes. As of yet, all engineering approaches applied to plastic degrading enzymes retain the natural $\alpha /\beta $-fold. While mutations can be used to increase thermostability, an inherent maximum likely exists for the $\alpha /\beta $-fold. It is thus of interest to introduce catalytic activity toward plastics in a different protein fold to escape the sequence space of plastic degrading enzymes. Here, a method for designing highly thermostable enzymes that can degrade plastics is described. With the help of Rosetta an active site catalysing the hydrolysis of polycarbonate is introduced into a set of thermostable scaffolds. Through computational evaluation, a potential PCase was selected and produced recombinantly in Escherichia coli. Thermal analysis suggests that the design has a melting temperature of >95$^{\circ }$C. Activity toward polycarbonate was confirmed using atomic force spectroscopy (AFM), proving the successful design of a PCase.
Collapse
Affiliation(s)
- Laura H Holst
- Material Science and Engineering Group, Department of Materials and Production, Aalborg University, 9000 Aalborg, Denmark
| | - Niklas G Madsen
- Material Science and Engineering Group, Department of Materials and Production, Aalborg University, 9000 Aalborg, Denmark
| | - Freja T Toftgård
- Material Science and Engineering Group, Department of Materials and Production, Aalborg University, 9000 Aalborg, Denmark
| | - Freja Rønne
- Material Science and Engineering Group, Department of Materials and Production, Aalborg University, 9000 Aalborg, Denmark
| | - Ioana-Malina Moise
- Material Science and Engineering Group, Department of Materials and Production, Aalborg University, 9000 Aalborg, Denmark
| | - Evamaria I Petersen
- Material Science and Engineering Group, Department of Materials and Production, Aalborg University, 9000 Aalborg, Denmark
| | - Peter Fojan
- Material Science and Engineering Group, Department of Materials and Production, Aalborg University, 9000 Aalborg, Denmark
| |
Collapse
|
8
|
Ma T, Wang R, Wang W, Gu W, Yuan Y, Zhang A, Wei J. Studies on the thermal degradation mechanism of polyethylene terephthalate and its 2-carboxy ethyl (phenyl) phosphinic acid copolymers. Polym Degrad Stab 2022. [DOI: 10.1016/j.polymdegradstab.2022.110185] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
9
|
Lavigne C, Gomes G, Pollice R, Aspuru-Guzik A. Guided discovery of chemical reaction pathways with imposed activation. Chem Sci 2022; 13:13857-13871. [PMID: 36544742 PMCID: PMC9710306 DOI: 10.1039/d2sc05135d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Accepted: 11/09/2022] [Indexed: 11/12/2022] Open
Abstract
Computational power and quantum chemical methods have improved immensely since computers were first applied to the study of reactivity, but the de novo prediction of chemical reactions has remained challenging. We show that complex reaction pathways can be efficiently predicted in a guided manner using chemical activation imposed by geometrical constraints of specific reactive modes, which we term imposed activation (IACTA). Our approach is demonstrated on realistic and challenging chemistry, such as a triple cyclization cascade involved in the total synthesis of a natural product, a water-mediated Michael addition, and several oxidative addition reactions of complex drug-like molecules. Notably and in contrast with traditional hand-guided computational chemistry calculations, our method requires minimal human involvement and no prior knowledge of the products or the associated mechanisms. We believe that IACTA will be a transformational tool to screen for chemical reactivity and to study both by-product formation and decomposition pathways in a guided way.
Collapse
Affiliation(s)
- Cyrille Lavigne
- Department of Computer Science, University of Toronto214 College St.TorontoOntarioM5T 3A1Canada
| | - Gabe Gomes
- Department of Computer Science, University of Toronto214 College St.TorontoOntarioM5T 3A1Canada,Chemical Physics Theory Group, Department of Chemistry, University of Toronto80 St George StTorontoOntarioM5S 3H6Canada
| | - Robert Pollice
- Department of Computer Science, University of Toronto214 College St.TorontoOntarioM5T 3A1Canada,Chemical Physics Theory Group, Department of Chemistry, University of Toronto80 St George StTorontoOntarioM5S 3H6Canada
| | - Alán Aspuru-Guzik
- Department of Computer Science, University of Toronto214 College St.TorontoOntarioM5T 3A1Canada,Chemical Physics Theory Group, Department of Chemistry, University of Toronto80 St George StTorontoOntarioM5S 3H6Canada,Department of Chemical Engineering & Applied Chemistry, University of Toronto200 College St.OntarioM5S 3E5Canada,Department of Materials Science & Engineering, University of Toronto184 College St.OntarioM5S 3E4Canada,Vector Institute for Artificial Intelligence661 University Ave Suite 710TorontoOntarioM5G 1M1Canada,Lebovic Fellow, Canadian Institute for Advanced Research (CIFAR)661 University AveTorontoOntarioM5GCanada
| |
Collapse
|
10
|
Abstract
This Perspective reviews the use of Transition Path Sampling methods to study enzymatically catalyzed chemical reactions. First applied by our group to an enzymatic reaction over 15 years ago, the method has uncovered basic principles in enzymatic catalysis such as the protein promoting vibration, and it has also helped harmonize such ideas as electrostatic preorganization with dynamic views of enzyme function. It is now being used to help uncover principles of protein design necessary to artificial enzyme creation.
Collapse
Affiliation(s)
- Steven D Schwartz
- Department of Chemistry and Biochemistry University of Arizona Tucson, Arizona 85721, United States
| |
Collapse
|
11
|
Snyder R, Kim B, Pan X, Shao Y, Pu J. Facilitating ab initio QM/MM free energy simulations by Gaussian process regression with derivative observations. Phys Chem Chem Phys 2022; 24:25134-25143. [PMID: 36222412 PMCID: PMC11095978 DOI: 10.1039/d2cp02820d] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023]
Abstract
In combined quantum mechanical and molecular mechanical (QM/MM) free energy simulations, how to synthesize the accuracy of ab initio (AI) methods with the speed of semiempirical (SE) methods for a cost-effective QM treatment remains a long-standing challenge. In this work, we present a machine-learning-facilitated method for obtaining AI/MM-quality free energy profiles through efficient SE/MM simulations. In particular, we use Gaussian process regression (GPR) to learn the energy and force corrections needed for SE/MM to match with AI/MM results during molecular dynamics simulations. Force matching is enabled in our model by including energy derivatives into the observational targets through the extended-kernel formalism. We demonstrate the effectiveness of this method on the solution-phase SN2 Menshutkin reaction using AM1/MM and B3LYP/6-31+G(d,p)/MM as the base and target levels, respectively. Trained on only 80 configurations sampled along the minimum free energy path (MFEP), the resulting GPR model reduces the average energy error in AM1/MM from 18.2 to 5.8 kcal mol-1 for the 4000-sample testing set with the average force error on the QM atoms decreased from 14.6 to 3.7 kcal mol-1 Å-1. Free energy sampling with the GPR corrections applied (AM1-GPR/MM) produces a free energy barrier of 14.4 kcal mol-1 and a reaction free energy of -34.1 kcal mol-1, in closer agreement with the AI/MM benchmarks and experimental results.
Collapse
Affiliation(s)
- Ryan Snyder
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN 46202, USA.
| | - Bryant Kim
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN 46202, USA.
| | - Xiaoliang Pan
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Pkwy, Norman, OK 73019, USA.
| | - Yihan Shao
- Department of Chemistry and Biochemistry, University of Oklahoma, 101 Stephenson Pkwy, Norman, OK 73019, USA.
| | - Jingzhi Pu
- Department of Chemistry and Chemical Biology, Indiana University-Purdue University Indianapolis, 402 N. Blackford St., Indianapolis, IN 46202, USA.
| |
Collapse
|
12
|
Rutt A, Shen JX, Horton M, Kim J, Lin J, Persson KA. Expanding the Material Search Space for Multivalent Cathodes. ACS APPLIED MATERIALS & INTERFACES 2022; 14:44367-44376. [PMID: 36137562 PMCID: PMC9542693 DOI: 10.1021/acsami.2c11733] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 07/01/2022] [Accepted: 09/05/2022] [Indexed: 06/16/2023]
Abstract
Multivalent batteries are an energy storage technology with the potential to surpass lithium-ion batteries; however, their performance have been limited by the low voltages and poor solid-state ionic mobility of available cathodes. A computational screening approach to identify high-performance multivalent intercalation cathodes among materials that do not contain the working ion of interest has been developed, which greatly expands the search space that can be considered for material discovery. This approach has been applied to magnesium cathodes as a proof of concept, and four resulting candidate materials [NASICON V2(PO4)3, birnessite NaMn4O8, tavorite MnPO4F, and spinel MnO2] are discussed in further detail. In examining the ion migration environment and associated Mg2+ migration energy in these materials, local energy maxima are found to correspond with pathway positions where Mg2+ passes through a plane of anion atoms. While previous studies have established the influence of local coordination on multivalent ion mobility, these results suggest that considering both the type of the local bonding environment and available free volume for the mobile ion along its migration pathway can be significant for improving solid-state mobility.
Collapse
Affiliation(s)
- Ann Rutt
- Department
of Materials Science and Engineering, University
of California, Berkeley California 94720, United States
| | - Jimmy-Xuan Shen
- Department
of Materials Science and Engineering, University
of California, Berkeley California 94720, United States
| | - Matthew Horton
- Materials
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley
California 94720, United
States
| | - Jiyoon Kim
- Department
of Materials Science and Engineering, University
of California, Berkeley California 94720, United States
| | - Jerry Lin
- Department
of Materials Science and Engineering, University
of California, Berkeley California 94720, United States
| | - Kristin A. Persson
- Department
of Materials Science and Engineering, University
of California, Berkeley California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley
California 94720, United
States
| |
Collapse
|
13
|
Hajibabaei A, Umer M, Anand R, Ha M, Kim KS. Fast atomic structure optimization with on-the-fly sparse Gaussian process potentials . JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2022; 34:344007. [PMID: 35675808 DOI: 10.1088/1361-648x/ac76ff] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/31/2022] [Accepted: 06/08/2022] [Indexed: 06/15/2023]
Abstract
We apply on-the-fly machine learning potentials (MLPs) using the sparse Gaussian process regression (SGPR) algorithm for fast optimization of atomic structures. Great acceleration is achieved even in the context of a single local optimization. Although for finding the exact local minimum, due to limited accuracy of MLPs, switching to another algorithm may be needed. For random gold clusters, the forces are reduced to ∼0.1 eV Å-1within less than ten first-principles (FP) calculations. Because of highly transferable MLPs, this algorithm is specially suitable for global optimization methods such as random or evolutionary structure searching or basin hopping. This is demonstrated by sequential optimization of random gold clusters for which, after only a few optimizations, FP calculations were rarely needed.
Collapse
Affiliation(s)
- Amir Hajibabaei
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Muhammad Umer
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Rohit Anand
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Miran Ha
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Republic of Korea
| | - Kwang S Kim
- Center for Superfunctional Materials, Department of Chemistry, Ulsan National Institute of Science and Technology, 50 UNIST-gil, Ulsan 44919, Republic of Korea
| |
Collapse
|
14
|
Abstract
Recent work has demonstrated the promise of using machine-learned surrogates, in particular, Gaussian process (GP) surrogates, in reducing the number of electronic structure calculations (ESCs) needed to perform surrogate model based (SMB) geometry optimization. In this paper, we study geometry meta-optimization with GP surrogates where a SMB optimizer additionally learns from its past "experience" performing geometry optimization. To validate this idea, we start with the simplest setting where a geometry meta-optimizer learns from previous optimizations of the same molecule with different initial-guess geometries. We give empirical evidence that geometry meta-optimization with GP surrogates is effective and requires less tuning compared to SMB optimization with GP surrogates on the ANI-1 dataset of off-equilibrium initial structures of small organic molecules. Unlike SMB optimization where a surrogate should be immediately useful for optimizing a given geometry, a surrogate in geometry meta-optimization has more flexibility because it can distribute its ESC savings across a set of geometries. Indeed, we find that GP surrogates that preserve rotational invariance provide increased marginal ESC savings across geometries. As a more stringent test, we also apply geometry meta-optimization to conformational search on a hand-constructed dataset of hydrocarbons and alcohols. We observe that while SMB optimization and geometry meta-optimization do save on ESCs, they also tend to miss higher energy conformers compared to standard geometry optimization. We believe that further research into characterizing the divergence between GP surrogates and potential energy surfaces is critical not only for advancing geometry meta-optimization but also for exploring the potential of machine-learned surrogates in geometry optimization in general.
Collapse
Affiliation(s)
- Daniel Huang
- Department of Computer Science, San Francisco State University, San Francisco, California 94132, USA
| | - Junwei Lucas Bao
- Department of Chemistry, Boston College, Chestnut Hill, Massachusetts 02467, USA
| | - Jean-Baptiste Tristan
- Department of Computer Science, Boston College, Chestnut Hill, Massachusetts 02467, USA
| |
Collapse
|
15
|
|
16
|
Zhan N, Kitchin JR. Model-Specific to Model-General Uncertainty for Physical Properties. Ind Eng Chem Res 2022. [DOI: 10.1021/acs.iecr.1c04706] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Ni Zhan
- Department of Chemical Engineering, Carnegie Mellon University, 5000 Forbes, Ave, Pittsburgh, Pennsylvania 15213, United States
| | - John R. Kitchin
- Department of Chemical Engineering, Carnegie Mellon University, 5000 Forbes, Ave, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
17
|
Chen L, Zhang X, Chen A, Yao S, Hu X, Zhou Z. Targeted design of advanced electrocatalysts by machine learning. CHINESE JOURNAL OF CATALYSIS 2022. [DOI: 10.1016/s1872-2067(21)63852-4] [Citation(s) in RCA: 12] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/19/2022]
|
18
|
Komp E, Janulaitis N, Valleau S. Progress towards machine learning reaction rate constants. Phys Chem Chem Phys 2021; 24:2692-2705. [PMID: 34935798 DOI: 10.1039/d1cp04422b] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/12/2023]
Abstract
Quantum and classical reaction rate constant calculations come at the cost of exploring potential energy surfaces. Due to the "curse of dimensionality", their evaluation quickly becomes unfeasible as the system size grows. Machine learning algorithms can accelerate the calculation of reaction rate constants by predicting them using low cost input features. In this perspective, we briefly introduce supervised machine learning algorithms in the context of reaction rate constant prediction. We discuss existing and recently created kinetic datasets and input feature representations as well as the use and design of machine learning algorithms to predict reaction rate constants or quantities required for their computation. Amongst these, we first describe the use of machine learning to predict activation, reaction, solvation and dissociation energies. We then look at the use of machine learning to predict reactive force field parameters, reaction rate constants as well as to help accelerate the search for minimum energy paths. Lastly, we provide an outlook on areas which have yet to be explored so as to improve and evaluate the use of machine learning algorithms for chemical reaction rate constants.
Collapse
Affiliation(s)
- Evan Komp
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, USA.
| | - Nida Janulaitis
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, USA.
| | - Stéphanie Valleau
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, USA.
| |
Collapse
|
19
|
Gerz I, Jannuzzi SAV, Hylland KT, Negri C, Wragg DS, Øien‐Ødegaard S, Tilset M, Olsbye U, DeBeer S, Amedjkouh M. Structural Elucidation, Aggregation, and Dynamic Behaviour of N,N,N,N-Copper(I) Schiff Base Complexes in Solid and in Solution: A Combined NMR, X-ray Spectroscopic and Crystallographic Investigation. Eur J Inorg Chem 2021; 2021:4762-4775. [PMID: 35874966 PMCID: PMC9298233 DOI: 10.1002/ejic.202100722] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2021] [Revised: 09/22/2021] [Indexed: 12/30/2022]
Abstract
A series of Cu(I) complexes of bidentate or tetradentate Schiff base ligands bearing either 1-H-imidazole or pyridine moieties were synthesized. The complexes were studied by a combination of NMR and X-ray spectroscopic techniques. The differences between the imidazole- and pyridine-based ligands were examined by 1H, 13C and 15N NMR spectroscopy. The magnitude of the 15Nimine coordination shifts was found to be strongly affected by the nature of the heterocycle in the complexes. These trends showed good correlation with the obtained Cu-Nimine bond lengths from single-crystal X-ray diffraction measurements. Variable-temperature NMR experiments, in combination with diffusion ordered spectroscopy (DOSY) revealed that one of the complexes underwent a temperature-dependent interconversion between a monomer, a dimer and a higher aggregate. The complexes bearing tetradentate imidazole ligands were further studied using Cu K-edge XAS and VtC XES, where DFT-assisted assignment of spectral features suggested that these complexes may form polynuclear oligomers in solid state. Additionally, the Cu(II) analogue of one of the complexes was incorporated into a metal-organic framework (MOF) as a way to obtain discrete, mononuclear complexes in the solid state.
Collapse
Affiliation(s)
- Isabelle Gerz
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| | | | - Knut T. Hylland
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| | - Chiara Negri
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| | - David S. Wragg
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| | - Sigurd Øien‐Ødegaard
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| | - Mats Tilset
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| | - Unni Olsbye
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| | - Serena DeBeer
- Department of Inorganic SpectroscopyMax Planck Institute for Chemical Energy ConversionStiftstraße 34–3645470Mülheim an der RuhrGermany
| | - Mohamed Amedjkouh
- Department of ChemistryUniversity of OsloP. O. Box 1033 Blindern0315OsloNorway
- Centre for Materials Science and NanotechnologyUniversity of OsloP.O. Box 1126 Blindern0316OsloNorway
| |
Collapse
|
20
|
Jin L, Ji Y, Wang H, Ding L, Li Y. First-principles materials simulation and design for alkali and alkaline metal ion batteries accelerated by machine learning. Phys Chem Chem Phys 2021; 23:21470-21483. [PMID: 34570138 DOI: 10.1039/d1cp02963k] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The challenge of regeneration of batteries requires a performance improvement in the alkali/alkaline metal ion battery (AMIB) materials, whereas the traditional research paradigm fully based on experiments and theoretical simulations needs massive research and development investment. During the last decade, machine learning (ML) has made breakthroughs in many complex disciplines, which testifies to their high processing speed and ability to capture relationships. Inspired by these achievements, ML has also been introduced to bring a new paradigm for shortening the development of AMIB materials. In this Perspective, the focus will be on how this new ML technology solves the key problems of redox potentials, ionic conductivity and stability parameters in first-principles materials' simulation and design for AMIBs. It is found that ML not only accelerates the property prediction, but also gives physicochemical insights into AMIB materials' design. In addition, the final part of this paper summarizes current achievements and looks forward to the progress of a novel paradigm in direct/inverse design with the increasing number of databases, skills, and ML technologies for AMIBs.
Collapse
Affiliation(s)
- Lujie Jin
- Institute of Functional Nano & Soft Materials (FUNSOM), Jiangsu Key Laboratory for Carbon-Based Functional Materials & Devices, Soochow University, Suzhou, Jiangsu 215123, China.
| | - Yujin Ji
- Institute of Functional Nano & Soft Materials (FUNSOM), Jiangsu Key Laboratory for Carbon-Based Functional Materials & Devices, Soochow University, Suzhou, Jiangsu 215123, China.
| | - Hongshuai Wang
- Institute of Functional Nano & Soft Materials (FUNSOM), Jiangsu Key Laboratory for Carbon-Based Functional Materials & Devices, Soochow University, Suzhou, Jiangsu 215123, China.
| | - Lifeng Ding
- Department of Chemistry, Xi'an JiaoTong-Liverpool University, 111 Ren'ai Road, Suzhou Dushu Lake, Higher Education Town, Jiangsu Province 215123, China
| | - Youyong Li
- Institute of Functional Nano & Soft Materials (FUNSOM), Jiangsu Key Laboratory for Carbon-Based Functional Materials & Devices, Soochow University, Suzhou, Jiangsu 215123, China. .,Macao Institute of Materials Science and Engineering, Macau University of Science and Technology, Taipa, Macau SAR 999078, China
| |
Collapse
|
21
|
Deringer VL, Bartók AP, Bernstein N, Wilkins DM, Ceriotti M, Csányi G. Gaussian Process Regression for Materials and Molecules. Chem Rev 2021; 121:10073-10141. [PMID: 34398616 PMCID: PMC8391963 DOI: 10.1021/acs.chemrev.1c00022] [Citation(s) in RCA: 232] [Impact Index Per Article: 77.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2021] [Indexed: 12/18/2022]
Abstract
We provide an introduction to Gaussian process regression (GPR) machine-learning methods in computational materials science and chemistry. The focus of the present review is on the regression of atomistic properties: in particular, on the construction of interatomic potentials, or force fields, in the Gaussian Approximation Potential (GAP) framework; beyond this, we also discuss the fitting of arbitrary scalar, vectorial, and tensorial quantities. Methodological aspects of reference data generation, representation, and regression, as well as the question of how a data-driven model may be validated, are reviewed and critically discussed. A survey of applications to a variety of research questions in chemistry and materials science illustrates the rapid growth in the field. A vision is outlined for the development of the methodology in the years to come.
Collapse
Affiliation(s)
- Volker L. Deringer
- Department
of Chemistry, Inorganic Chemistry Laboratory, University of Oxford, Oxford OX1 3QR, United Kingdom
| | - Albert P. Bartók
- Department
of Physics and Warwick Centre for Predictive Modelling, School of
Engineering, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Noam Bernstein
- Center
for Computational Materials Science, U.S.
Naval Research Laboratory, Washington D.C. 20375, United States
| | - David M. Wilkins
- Atomistic
Simulation Centre, School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, Northern Ireland, United Kingdom
| | - Michele Ceriotti
- Laboratory
of Computational Science and Modeling, IMX, École Polytechnique Fédérale de Lausanne, Lausanne 1015, Switzerland
- National
Centre for Computational Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale
de Lausanne, Lausanne, Switzerland
| | - Gábor Csányi
- Engineering
Laboratory, University of Cambridge, Cambridge CB2 1PZ, United Kingdom
| |
Collapse
|
22
|
Westermayr J, Marquetand P. Machine Learning for Electronically Excited States of Molecules. Chem Rev 2021; 121:9873-9926. [PMID: 33211478 PMCID: PMC8391943 DOI: 10.1021/acs.chemrev.0c00749] [Citation(s) in RCA: 167] [Impact Index Per Article: 55.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2020] [Indexed: 12/11/2022]
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna
Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data
Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
23
|
Keith JA, Vassilev-Galindo V, Cheng B, Chmiela S, Gastegger M, Müller KR, Tkatchenko A. Combining Machine Learning and Computational Chemistry for Predictive Insights Into Chemical Systems. Chem Rev 2021; 121:9816-9872. [PMID: 34232033 PMCID: PMC8391798 DOI: 10.1021/acs.chemrev.1c00107] [Citation(s) in RCA: 190] [Impact Index Per Article: 63.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/04/2021] [Indexed: 12/23/2022]
Abstract
Machine learning models are poised to make a transformative impact on chemical sciences by dramatically accelerating computational algorithms and amplifying insights available from computational chemistry methods. However, achieving this requires a confluence and coaction of expertise in computer science and physical sciences. This Review is written for new and experienced researchers working at the intersection of both fields. We first provide concise tutorials of computational chemistry and machine learning methods, showing how insights involving both can be achieved. We follow with a critical review of noteworthy applications that demonstrate how computational chemistry and machine learning can be used together to provide insightful (and useful) predictions in molecular and materials modeling, retrosyntheses, catalysis, and drug design.
Collapse
Affiliation(s)
- John A. Keith
- Department
of Chemical and Petroleum Engineering Swanson School of Engineering, University of Pittsburgh, Pittsburgh, Pennsylvania 15261, United States
| | - Valentin Vassilev-Galindo
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| | - Bingqing Cheng
- Accelerate
Programme for Scientific Discovery, Department
of Computer Science and Technology, 15 J. J. Thomson Avenue, Cambridge CB3 0FD, United Kingdom
| | - Stefan Chmiela
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Michael Gastegger
- Department
of Software Engineering and Theoretical Computer Science, Technische Universität Berlin, 10587, Berlin, Germany
| | - Klaus-Robert Müller
- Machine
Learning Group, Technische Universität
Berlin, 10587, Berlin, Germany
- Department
of Artificial Intelligence, Korea University, Anam-dong, Seongbuk-gu, Seoul, 02841, Korea
- Max-Planck-Institut für Informatik, 66123 Saarbrücken, Germany
- Google Research, Brain Team, 10117 Berlin, Germany
| | - Alexandre Tkatchenko
- Department
of Physics and Materials Science, University
of Luxembourg, L-1511 Luxembourg City, Luxembourg
| |
Collapse
|
24
|
Nandy A, Duan C, Taylor MG, Liu F, Steeves AH, Kulik HJ. Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning. Chem Rev 2021; 121:9927-10000. [PMID: 34260198 DOI: 10.1021/acs.chemrev.1c00347] [Citation(s) in RCA: 70] [Impact Index Per Article: 23.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.
Collapse
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adam H Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
25
|
Abstract
Electronically excited states of molecules are at the heart of photochemistry, photophysics, as well as photobiology and also play a role in material science. Their theoretical description requires highly accurate quantum chemical calculations, which are computationally expensive. In this review, we focus on not only how machine learning is employed to speed up such excited-state simulations but also how this branch of artificial intelligence can be used to advance this exciting research field in all its aspects. Discussed applications of machine learning for excited states include excited-state dynamics simulations, static calculations of absorption spectra, as well as many others. In order to put these studies into context, we discuss the promises and pitfalls of the involved machine learning techniques. Since the latter are mostly based on quantum chemistry calculations, we also provide a short introduction into excited-state electronic structure methods and approaches for nonadiabatic dynamics simulations and describe tricks and problems when using them in machine learning for excited states of molecules.
Collapse
Affiliation(s)
- Julia Westermayr
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
| | - Philipp Marquetand
- Institute of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Vienna Research Platform on Accelerating Photoreaction Discovery, University of Vienna, Währinger Strasse 17, 1090 Vienna, Austria
- Data Science @ Uni Vienna, University of Vienna, Währinger Strasse 29, 1090 Vienna, Austria
| |
Collapse
|
26
|
Born D, Kästner J. Geometry Optimization in Internal Coordinates Based on Gaussian Process Regression: Comparison of Two Approaches. J Chem Theory Comput 2021; 17:5955-5967. [PMID: 34378918 DOI: 10.1021/acs.jctc.1c00517] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Geometry optimization based on Gaussian process regression (GPR) was extended to internal coordinates. We used delocalized internal coordinates composed of distances and several types of angles and compared two methods of including them. In both cases, the GPR surrogate surface is trained on geometries in internal coordinates. In one case, it predicts the gradient in Cartesian coordinates and in the other, in internal coordinates. We tested both methods on a set of 30 small molecules and one larger Rh complex taken from the study of a catalytic mechanism. The former method is slightly more efficient, while the latter method is somewhat more robust. Both methods reduce the number of required optimization steps compared to GPR in Cartesian coordinates or the standard L-BFGS optimizer. We found it advantageous to use automatically adjusted hyperparameters to optimize them.
Collapse
Affiliation(s)
- Daniel Born
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
27
|
Xu J, Cao XM, Hu P. Perspective on computational reaction prediction using machine learning methods in heterogeneous catalysis. Phys Chem Chem Phys 2021; 23:11155-11179. [PMID: 33972971 DOI: 10.1039/d1cp01349a] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022]
Abstract
Heterogeneous catalysis plays a significant role in the modern chemical industry. Towards the rational design of novel catalysts, understanding reactions over surfaces is the most essential aspect. Typical industrial catalytic processes such as syngas conversion and methane utilisation can generate a large reaction network comprising thousands of intermediates and reaction pairs. This complexity not only arises from the permutation of transformations between species but also from the extra reaction channels offered by distinct surface sites. Despite the success in investigating surface reactions at the atomic scale, the huge computational expense of ab initio methods hinders the exploration of such complicated reaction networks. With the proliferation of catalysis studies, machine learning as an emerging tool can take advantage of the accumulated reaction data to emulate the output of ab initio methods towards swift reaction prediction. Here, we briefly summarise the conventional workflow of reaction prediction, including reaction network generation, ab initio thermodynamics and microkinetic modelling. An overview of the frequently used regression models in machine learning is presented. As a promising alternative to full ab initio calculations, machine learning interatomic potentials are highlighted. Furthermore, we survey applications assisted by these methods for accelerating reaction prediction, exploring reaction networks, and computational catalyst design. Finally, we envisage future directions in computationally investigating reactions and implementing machine learning algorithms in heterogeneous catalysis.
Collapse
Affiliation(s)
- Jiayan Xu
- Key Laboratory for Advanced Materials and Joint International Research Laboratory of Precision Chemistry and Molecular Engineering, Feringa Nobel Prize Scientist Joint Research Center, Frontiers Science Center for Materiobiology and Dynamic Chemistry, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 130 Meilong Road, Shanghai, 200237, P. R. China. and School of Chemistry and Chemical Engineering, Queen's University Belfast, Belfast BT9 5AG, UK
| | - Xiao-Ming Cao
- Key Laboratory for Advanced Materials and Joint International Research Laboratory of Precision Chemistry and Molecular Engineering, Feringa Nobel Prize Scientist Joint Research Center, Frontiers Science Center for Materiobiology and Dynamic Chemistry, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 130 Meilong Road, Shanghai, 200237, P. R. China.
| | - P Hu
- Key Laboratory for Advanced Materials and Joint International Research Laboratory of Precision Chemistry and Molecular Engineering, Feringa Nobel Prize Scientist Joint Research Center, Frontiers Science Center for Materiobiology and Dynamic Chemistry, Centre for Computational Chemistry and Research Institute of Industrial Catalysis, School of Chemistry and Molecular Engineering, East China University of Science and Technology, 130 Meilong Road, Shanghai, 200237, P. R. China. and School of Chemistry and Chemical Engineering, Queen's University Belfast, Belfast BT9 5AG, UK
| |
Collapse
|
28
|
Ozkan A, Sitharam M, Flores-Canales JC, Prabhu R, Kurnikova M. Baseline Comparisons of Complementary Sampling Methods for Assembly Driven by Short-Ranged Pair Potentials toward Fast and Flexible Hybridization. J Chem Theory Comput 2021; 17:1967-1987. [PMID: 33576635 DOI: 10.1021/acs.jctc.0c00945] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
This work measures baseline sampling characteristics that highlight fundamental differences between sampling methods for assembly driven by short-ranged pair potentials. Such granular comparison is essential for fast, flexible, and accurate hybridization of complementary methods. Besides sampling speed, efficiency, and accuracy of uniform grid coverage, other sampling characteristics measured are (i) accuracy of covering narrow low energy regions that have low effective dimension (ii) ability to localize sampling to specific basins, and (iii) flexibility in sampling distributions. As a proof of concept, we compare a recently developed geometric methodology EASAL (Efficient Atlasing and Search of Assembly Landscapes) and the traditional Monte Carlo (MC) method for sampling the energy landscape of two assembling trans-membrane helices, driven by short-range pair potentials. By measuring the above-mentioned sampling characteristics, we demonstrate that EASAL provides localized and accurate coverage of crucial regions of the energy landscape of low effective dimension, under flexible sampling distributions, with much fewer samples and computational resources than MC sampling. EASAL's empirically validated theoretical guarantees permit credible extrapolation of these measurements and comparisons to arbitrary number and size of assembling units. Promising avenues for hybridizing the complementary advantages of the two methods are discussed.
Collapse
Affiliation(s)
- Aysegul Ozkan
- CISE Department, University of Florida, Gainesville, Florida 32611-6120, United States
| | - Meera Sitharam
- CISE Department, University of Florida, Gainesville, Florida 32611-6120, United States
| | | | - Rahul Prabhu
- CISE Department, University of Florida, Gainesville, Florida 32611-6120, United States
| | - Maria Kurnikova
- Chemistry Department, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, United States
| |
Collapse
|
29
|
|
30
|
Fdez Galván I, Raggi G, Lindh R. Restricted-Variance Constrained, Reaction Path, and Transition State Molecular Optimizations Using Gradient-Enhanced Kriging. J Chem Theory Comput 2020; 17:571-582. [PMID: 33382621 PMCID: PMC7871327 DOI: 10.1021/acs.jctc.0c01163] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Gaussian process
regression has recently been explored as an alternative
to standard surrogate models in molecular equilibrium geometry optimization.
In particular, the gradient-enhanced Kriging approach in association
with internal coordinates, restricted-variance optimization, and an
efficient and fast estimate of hyperparameters has demonstrated performance
on par or better than standard methods. In this report, we extend
the approach to constrained optimizations and transition states and
benchmark it for a set of reactions. We compare the performance of
the newly developed method with the standard techniques in the location
of transition states and in constrained optimizations, both isolated
and in the context of reaction path computation. The results show
that the method outperforms the current standard in efficiency as
well as in robustness.
Collapse
Affiliation(s)
| | - Gerardo Raggi
- Department of Chemistry - BMC, Uppsala University, Uppsala 75123, Sweden
| | - Roland Lindh
- Department of Chemistry - BMC, Uppsala University, Uppsala 75123, Sweden
| |
Collapse
|
31
|
Garijo del Río E, Kaappa S, Garrido Torres JA, Bligaard T, Jacobsen KW. Machine learning with bond information for local structure optimizations in surface science. J Chem Phys 2020; 153:234116. [DOI: 10.1063/5.0033778] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Affiliation(s)
| | - Sami Kaappa
- Department of Physics, Technical University of Denmark, Kgs. Lyngby, Denmark
| | - José A. Garrido Torres
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, California 94025, USA
- Columbia Electrochemical Energy Center, Department of Chemical Engineering, Columbia University, New York, New York 10027, USA
| | - Thomas Bligaard
- SUNCAT Center for Interface Science and Catalysis, Department of Chemical Engineering, Stanford University, Stanford, California 94305, USA
- Department of Energy Conversion and Storage, Technical University of Denmark, Kgs. Lyngby, Denmark
| | | |
Collapse
|
32
|
Tran A, Tranchida J, Wildey T, Thompson AP. Multi-fidelity machine-learning with uncertainty quantification and Bayesian optimization for materials design: Application to ternary random alloys. J Chem Phys 2020; 153:074705. [PMID: 32828092 DOI: 10.1063/5.0015672] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a scale-bridging approach based on a multi-fidelity (MF) machine-learning (ML) framework leveraging Gaussian processes (GP) to fuse atomistic computational model predictions across multiple levels of fidelity. Through the posterior variance of the MFGP, our framework naturally enables uncertainty quantification, providing estimates of confidence in the predictions. We used density functional theory as high-fidelity prediction, while a ML interatomic potential is used as low-fidelity prediction. Practical materials' design efficiency is demonstrated by reproducing the ternary composition dependence of a quantity of interest (bulk modulus) across the full aluminum-niobium-titanium ternary random alloy composition space. The MFGP is then coupled to a Bayesian optimization procedure, and the computational efficiency of this approach is demonstrated by performing an on-the-fly search for the global optimum of bulk modulus in the ternary composition space. The framework presented in this manuscript is the first application of MFGP to atomistic materials simulations fusing predictions between density functional theory and classical interatomic potential calculations.
Collapse
Affiliation(s)
- Anh Tran
- Optimization and Uncertainty Quantification, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Julien Tranchida
- Computational Multiscale, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Tim Wildey
- Optimization and Uncertainty Quantification, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| | - Aidan P Thompson
- Computational Multiscale, Center for Computing Research, Sandia National Laboratories, Albuquerque, New Mexico 87123, USA
| |
Collapse
|
33
|
Denzel A, Kästner J. Hessian Matrix Update Scheme for Transition State Search Based on Gaussian Process Regression. J Chem Theory Comput 2020; 16:5083-5089. [DOI: 10.1021/acs.jctc.0c00348] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Affiliation(s)
- Alexander Denzel
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| | - Johannes Kästner
- Institute for Theoretical Chemistry, University of Stuttgart, Pfaffenwaldring 55, 70569 Stuttgart, Germany
| |
Collapse
|
34
|
Raggi G, Galván IF, Ritterhoff CL, Vacher M, Lindh R. Restricted-Variance Molecular Geometry Optimization Based on Gradient-Enhanced Kriging. J Chem Theory Comput 2020; 16:3989-4001. [PMID: 32374164 PMCID: PMC7304864 DOI: 10.1021/acs.jctc.0c00257] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
Machine learning techniques, specifically
gradient-enhanced Kriging
(GEK), have been implemented for molecular geometry optimization.
GEK-based optimization has many advantages compared to conventional—step-restricted
second-order truncated expansion—molecular optimization methods.
In particular, the surrogate model given by GEK can have multiple
stationary points, will smoothly converge to the exact model as the
number of sample points increases, and contains an explicit expression
for the expected error of the model function at an arbitrary point.
Machine learning is, however, associated with abundance of data, contrary
to the situation desired for efficient geometry optimizations. In
this paper, we demonstrate how the GEK procedure can be utilized in
a fashion such that in the presence of few data points, the surrogate
surface will in a robust way guide the optimization to a minimum of
a potential energy surface. In this respect, the GEK procedure will
be used to mimic the behavior of a conventional second-order scheme
but retaining the flexibility of the superior machine learning approach.
Moreover, the expected error will be used in the optimizations to
facilitate restricted-variance optimizations. A procedure which relates
the eigenvalues of the approximate guessed Hessian with the individual
characteristic lengths, used in the GEK model, reduces the number
of empirical parameters to optimize to two: the value of the trend
function and the maximum allowed variance. These parameters are determined
using the extended Baker (e-Baker) and part of the Baker transition-state
(Baker-TS) test suites as a training set. The so-created optimization
procedure is tested using the e-Baker, full Baker-TS, and S22 test
suites, at the density functional theory and second-order Møller–Plesset
levels of approximation. The results show that the new method is generally
of similar or better performance than a state-of-the-art conventional
method, even for cases where no significant improvement was expected.
Collapse
Affiliation(s)
- Gerardo Raggi
- Department of Chemistry-BMC, Uppsala University, 751 23 Uppsala, Sweden
| | | | - Christian L Ritterhoff
- Department of Chemistry-BMC, Uppsala University, 751 23 Uppsala, Sweden.,Faculty of Science, Universität Erlangen-Nürnberg, 91054 Erlangen, Germany
| | - Morgane Vacher
- Department of Chemistry-Ångström Laboratory, Uppsala University, 751 21 Uppsala, Sweden.,Université de Nantes, CNRS, CEISAM UMR 6230, F-44000 Nantes, France
| | - Roland Lindh
- Department of Chemistry-BMC, Uppsala University, 751 23 Uppsala, Sweden
| |
Collapse
|
35
|
Abstract
As the quantum chemistry (QC) community embraces machine learning (ML), the number of new methods and applications based on the combination of QC and ML is surging. In this Perspective, a view of the current state of affairs in this new and exciting research field is offered, challenges of using machine learning in quantum chemistry applications are described, and potential future developments are outlined. Specifically, examples of how machine learning is used to improve the accuracy and accelerate quantum chemical research are shown. Generalization and classification of existing techniques are provided to ease the navigation in the sea of literature and to guide researchers entering the field. The emphasis of this Perspective is on supervised machine learning.
Collapse
Affiliation(s)
- Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| |
Collapse
|
36
|
Meyer R, Hauser AW. Geometry optimization using Gaussian process regression in internal coordinate systems. J Chem Phys 2020; 152:084112. [PMID: 32113346 DOI: 10.1063/1.5144603] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Locating the minimum energy structure of molecules, typically referred to as geometry optimization, is one of the first steps of any computational chemistry calculation. Earlier research was mostly dedicated to finding convenient sets of molecule-specific coordinates for a suitable representation of the potential energy surface, where a faster convergence toward the minimum structure can be achieved. More recent approaches, on the other hand, are based on various machine learning techniques and seem to revert to Cartesian coordinates instead for practical reasons. We show that the combination of Gaussian process regression with those coordinate systems employed by state-of-the-art geometry optimizers can significantly improve the performance of this powerful machine learning technique. This is demonstrated on a benchmark set of 30 small covalently bonded molecules.
Collapse
Affiliation(s)
- Ralf Meyer
- Institute of Experimental Physics, Graz University of Technology, Petersgasse 16, 8010 Graz, Austria
| | - Andreas W Hauser
- Institute of Experimental Physics, Graz University of Technology, Petersgasse 16, 8010 Graz, Austria
| |
Collapse
|
37
|
Koistinen OP, Ásgeirsson V, Vehtari A, Jónsson H. Minimum Mode Saddle Point Searches Using Gaussian Process Regression with Inverse-Distance Covariance Function. J Chem Theory Comput 2019; 16:499-509. [DOI: 10.1021/acs.jctc.9b01038] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Olli-Pekka Koistinen
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík 107, Iceland
| | - Vilhjálmur Ásgeirsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík 107, Iceland
| | | | - Hannes Jónsson
- Science Institute and Faculty of Physical Sciences, University of Iceland, Reykjavík 107, Iceland
| |
Collapse
|