1
|
Rehn DR, Fink A, Dempwolff AL, Dreuw A. Analytical Gradients for Electron-Attached and Ionized States for the Algebraic-Diagrammatic Construction Scheme for the Electron Propagator up to Third Order. J Phys Chem A 2024; 128:8795-8802. [PMID: 39320963 DOI: 10.1021/acs.jpca.4c04435] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/27/2024]
Abstract
The derivation and implementation of analytical gradients for methods based on the non-Dyson algebraic diagrammatic construction for the electron propagator, IP-ADC and EA-ADC, up to the third order is presented. Using nuclear gradients, ground-state equilibrium structures for small open-shell systems are calculated. In addition, we investigated the performance of IP/EA-ADC methods for the calculation of adiabatic ionization potentials and electron affinities for medium-sized organic molecules.
Collapse
Affiliation(s)
- Dirk R Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Fink
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
2
|
Stocks R, Palethorpe E, Barca GMJ. Multi-GPU RI-HF Energies and Analytic Gradients─Toward High-Throughput Ab Initio Molecular Dynamics. J Chem Theory Comput 2024; 20:7503-7515. [PMID: 39192710 DOI: 10.1021/acs.jctc.4c00877] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/29/2024]
Abstract
This article presents an optimized algorithm and implementation for calculating resolution-of-the-identity Hartree-Fock (RI-HF) energies and analytic gradients using multiple graphics processing units (GPUs). The algorithm is especially designed for high throughput ab initio molecular dynamics simulations of small and medium size molecules (10-100 atoms). Key innovations of this work include the exploitation of multi-GPU parallelism and a workload balancing scheme that efficiently distributes computational tasks among GPUs. Our implementation also employs techniques for symmetry utilization, integral screening, and leveraging sparsity to optimize memory usage. Computational results show that the implementation achieves significant performance improvements, including over 3 × speedups in single GPU AIMD throughput compared to previous GPU-accelerated RI-HF and traditional HF methods. Furthermore, utilizing multiple GPUs can provide superlinear speedup when the additional aggregate GPU memory allows for the storage of decompressed three-center integrals.
Collapse
Affiliation(s)
- Ryan Stocks
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| | - Elise Palethorpe
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| | - Giuseppe M J Barca
- School of Computing and Information Systems, Melbourne University, Melbourne, VIC 3052, Australia
- QDX Technologies, Dickson, ACT 2602, Australia
| |
Collapse
|
3
|
Iyer GR, Whelpley N, Tiihonen J, Kent PRC, Krogel JT, Rubenstein BM. Force-Free Identification of Minimum-Energy Pathways and Transition States for Stochastic Electronic Structure Theories. J Chem Theory Comput 2024; 20:7416-7429. [PMID: 39172163 DOI: 10.1021/acs.jctc.4c00214] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 08/23/2024]
Abstract
The accurate mapping of potential energy surfaces (PESs) is crucial to our understanding of the numerous physical and chemical processes mediated by atomic rearrangements, such as conformational changes and chemical reactions, and the thermodynamic and kinetic feasibility of these processes. Stochastic electronic structure theories, e.g., Quantum Monte Carlo (QMC) methods, enable highly accurate total energy calculations that in principle can be used to construct the PES. However, their stochastic nature poses a challenge to the computation and use of forces and Hessians, which are typically required in algorithms for minimum-energy pathway (MEP) and transition state (TS) identification, such as the nudged elastic band (NEB) algorithm and its climbing image formulation. Here, we present strategies that utilize the surrogate Hessian line-search method, previously developed for QMC structural optimization, to efficiently identify MEP and TS structures without requiring force calculations at the level of the stochastic electronic structure theory. By modifying the surrogate Hessian algorithm to operate in path-orthogonal subspaces and at saddle points, we show that it is possible to identify MEPs and TSs by using a force-free QMC approach. We demonstrate these strategies via two examples, the inversion of the ammonia (NH3) molecule and the nucleophilic substitution (SN2) reaction F- + CH3F → FCH3 + F-. We validate our results using Density Functional Theory (DFT)- and Coupled Cluster (CCSD, CCSD(T))-based NEB calculations. We then introduce a hybrid DFT-QMC approach to compute thermodynamic and kinetic quantities, free energy differences, rate constants, and equilibrium constants that incorporates stochastically optimized structures and their energies, and show that this scheme improves upon DFT accuracy. Our methods generalize straightforwardly to other systems and other high-accuracy theories that similarly face challenges computing energy gradients, paving the way for highly accurate PES mapping, transition state determination, and thermodynamic and kinetic calculations at significantly reduced computational expense.
Collapse
Affiliation(s)
- Gopal R Iyer
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Noah Whelpley
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| | - Juha Tiihonen
- Department of Physics, Nanoscience Center, University of Jyväskylä, Jyväskylä 40014, Finland
| | - Paul R C Kent
- Computational Sciences and Engineering Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Jaron T Krogel
- Materials Science and Technology Division, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Brenda M Rubenstein
- Department of Chemistry, Brown University, Providence, Rhode Island 02912, United States
| |
Collapse
|
4
|
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
|
5
|
Sanz García J, Maskri R, Mitrushchenkov A, Joubert-Doriol L. Optimizing Conical Intersections without Explicit Use of Non-Adiabatic Couplings. J Chem Theory Comput 2024; 20:5643-5654. [PMID: 38888629 DOI: 10.1021/acs.jctc.4c00326] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/20/2024]
Abstract
We present two alternative methods for optimizing minimum energy conical intersection (MECI) molecular geometries without knowledge of the derivative coupling (DC). These methods are based on the utilization of Lagrange multipliers: (i) one method uses an approximate calculation of the DC, while the other (ii) do not require the DC. Both methods use the fact that information on the DC is contained in the Hessian of the squared energy difference. Tests done on a set of small molecular systems, in comparison with other methods, show the ability of the proposed methods to optimize MECIs. Finally, we apply the methods to the furimamide molecule, to optimize and characterize its S1/S2 MECI, and to optimizing the S0/S1 MECI of the silver trimer.
Collapse
Affiliation(s)
- Juan Sanz García
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Rosa Maskri
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Alexander Mitrushchenkov
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| | - Loïc Joubert-Doriol
- Univ Gustave Eiffel, Univ Paris Est Creteil, CNRS, UMR 8208, MSME, F-77454 Marne-la-Vallée, France
| |
Collapse
|
6
|
Fantasia A, Rovaris F, Abou El Kheir O, Marzegalli A, Lanzoni D, Pessina L, Xiao P, Zhou C, Li L, Henkelman G, Scalise E, Montalenti F. Development of a machine learning interatomic potential for exploring pressure-dependent kinetics of phase transitions in germanium. J Chem Phys 2024; 161:014110. [PMID: 38953439 DOI: 10.1063/5.0214588] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2024] [Accepted: 06/15/2024] [Indexed: 07/04/2024] Open
Abstract
We introduce a data-driven potential aimed at the investigation of pressure-dependent phase transitions in bulk germanium, including the estimate of kinetic barriers. This is achieved by suitably building a database including several configurations along minimum energy paths, as computed using the solid-state nudged elastic band method. After training the model based on density functional theory (DFT)-computed energies, forces, and stresses, we provide validation and rigorously test the potential on unexplored paths. The resulting agreement with the DFT calculations is remarkable in a wide range of pressures. The potential is exploited in large-scale isothermal-isobaric simulations, displaying local nucleation in the R8 to β-Sn pressure-induced phase transformation, taken here as an illustrative example.
Collapse
Affiliation(s)
- A Fantasia
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| | - F Rovaris
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| | - O Abou El Kheir
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| | - A Marzegalli
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| | - D Lanzoni
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| | - L Pessina
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| | - P Xiao
- Department of Physics and Atmospheric Science, Dalhousie University, 1453 Lord Dalhousie Drive, Halifax, Nova Scotia B3H 4R2, Canada
| | - C Zhou
- Department of Materials Science and Engineering, Southern University of Science and Technology, 1088 Xueyuan Avenue, 518055 Shenzhen, China
| | - L Li
- Department of Materials Science and Engineering, Southern University of Science and Technology, 1088 Xueyuan Avenue, 518055 Shenzhen, China
| | - G Henkelman
- Department of Chemistry, The University of Texas at Austin, 105 East 24th Street STOP A5300 Austin, Texas 78712, USA
| | - E Scalise
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| | - F Montalenti
- Department of Materials Science, University of Milano-Bicocca, 20125 Milano, Italy
| |
Collapse
|
7
|
Barhoumi A, Ryachi K, Belghiti ME, Chafi M, Tounsi A, Syed A, Idrissi ME, Wong LS, Zeroual A. Chromatography Scrutiny, Molecular Docking, Clarifying the Selectivities and the Mechanism of [3 + 2] Cycloloaddition Reaction between Linallol and Chlorobenzene-Nitrile-oxide. J Fluoresc 2024; 34:1913-1929. [PMID: 37668770 DOI: 10.1007/s10895-023-03411-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2023] [Accepted: 08/23/2023] [Indexed: 09/06/2023]
Abstract
Employing the Molecular Electron Density Theory, [3 + 2] cycloaddition processes between 4-chlorobenzenenitrileoxide and linalool, have been applied using the DFT/B3LYP/6-311(d,p) method, activation, reaction energies and the reactivity indices are calculated. In an investigation of conceptual DFT indices, LIL-1 will contribute to this reaction as a nucleophile, whilst NOX-2 will participate as an electrophile. This cyclization is regio, chemo and stereospecific, as demonstrated by the reaction and activation energies, in clear agreement with the experiment's results, in addition, ELF analysis revealed that the mechanism for this cycloaddition occurs in two steps. Furthermore, a docking study was conducted on the products studied, and the interaction with the protein protease COVID-19 (PDB ID: 6LU7), our results indicate that the presence of the -OH group increases the affinity of these products, moreover, adsorption study by chromatography was made on silica gel as support; our outcome reveals that the -OH group creates an intramolecular hydrogen bond in the product P2, while in the product P3 will create a hydrogen bond with the silica gel which makes the two products P2 and P3 are very easy to separate by chromatography, this result is in excellent agreement with the Rf retention value. The study might provide a fundamental for developing natural anti-viral compound in promoting human health.
Collapse
Affiliation(s)
- Ali Barhoumi
- Molecular Modelling and Spectroscopy Research Team, Faculty of Science, Chouaïb Doukkali University, P.O. Box 20, 24000, El Jadida, Morocco
| | - Kamal Ryachi
- Agro-Industrial, Environmental and Ecological Processes Team, Faculty of Science and Techniques of Beni Mellal, Sultan Moulay Slimane University, Beni Mellal, Morocco
| | - Mohammed Elalaoui Belghiti
- Laboratory of Physical Chemistry of Materials, Ben M'Sick Faculty of Sciences, Hassan II University, Casablanca, Morocco
- Laboratory of Nernest Technology, 163 Willington Street, Sherbrooke, QC J1H5C7, Canada
| | - Mohammed Chafi
- LIPE, Higher School of Technology, Hassan II University, Casablanca, Morocco
| | - Abdessamad Tounsi
- Agro-Industrial, Environmental and Ecological Processes Team, Faculty of Science and Techniques of Beni Mellal, Sultan Moulay Slimane University, Beni Mellal, Morocco
| | - Asad Syed
- Department of Botany and Microbiology, College of Science, King Saud University, P.O. Box 2455, 11451, Riyadh, Saudi Arabia
| | - Mohammed El Idrissi
- Team of Chemical Processes and Applied Materials, Faculty Polydisciplinary, Sultan Moulay Slimane University, Beni-Mellal, Morocco.
| | - Ling Shing Wong
- Faculty of Health and Life Sciences, INTI International University, Putra Nilai, 71800, Nilai, Negeri Sembilan, Malaysia
| | - Abdellah Zeroual
- Molecular Modelling and Spectroscopy Research Team, Faculty of Science, Chouaïb Doukkali University, P.O. Box 20, 24000, El Jadida, Morocco
| |
Collapse
|
8
|
Weymuth T, Unsleber JP, Türtscher PL, Steiner M, Sobez JG, Müller CH, Mörchen M, Klasovita V, Grimmel SA, Eckhoff M, Csizi KS, Bosia F, Bensberg M, Reiher M. SCINE-Software for chemical interaction networks. J Chem Phys 2024; 160:222501. [PMID: 38857173 DOI: 10.1063/5.0206974] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2024] [Accepted: 05/09/2024] [Indexed: 06/12/2024] Open
Abstract
The software for chemical interaction networks (SCINE) project aims at pushing the frontier of quantum chemical calculations on molecular structures to a new level. While calculations on individual structures as well as on simple relations between them have become routine in chemistry, new developments have pushed the frontier in the field to high-throughput calculations. Chemical relations may be created by a search for specific molecular properties in a molecular design attempt, or they can be defined by a set of elementary reaction steps that form a chemical reaction network. The software modules of SCINE have been designed to facilitate such studies. The features of the modules are (i) general applicability of the applied methodologies ranging from electronic structure (no restriction to specific elements of the periodic table) to microkinetic modeling (with little restrictions on molecularity), full modularity so that SCINE modules can also be applied as stand-alone programs or be exchanged for external software packages that fulfill a similar purpose (to increase options for computational campaigns and to provide alternatives in case of tasks that are hard or impossible to accomplish with certain programs), (ii) high stability and autonomous operations so that control and steering by an operator are as easy as possible, and (iii) easy embedding into complex heterogeneous environments for molecular structures taken individually or in the context of a reaction network. A graphical user interface unites all modules and ensures interoperability. All components of the software have been made available as open source and free of charge.
Collapse
Affiliation(s)
- Thomas Weymuth
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan P Unsleber
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Paul L Türtscher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Miguel Steiner
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jan-Grimo Sobez
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Charlotte H Müller
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Maximilian Mörchen
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Veronika Klasovita
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Stephanie A Grimmel
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Marco Eckhoff
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Katja-Sophia Csizi
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Francesco Bosia
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Moritz Bensberg
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Markus Reiher
- ETH Zurich, Department of Chemistry and Applied Biosciences, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
9
|
Panadés-Barrueta RL, Duflot D, Soto J, Martínez-Núñez E, Peláez D. Automatic Determination of the Non-Covalent Stable Conformations of the NO 2-Pyrene Cluster in Full Dimensionality (81D) Using the vdW-TSSCDS Approach. Chemphyschem 2024; 25:e202301001. [PMID: 38662437 DOI: 10.1002/cphc.202301001] [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: 12/28/2023] [Revised: 02/27/2024] [Indexed: 05/24/2024]
Abstract
We present the detailed topographical characterisation (stationary points and minimum energy paths connecting them) of the full dimensional (81D) intermolecular potential energy surface associated with the non-covalent interactions between the NO2 radical and the pyrene (C16H10) molecule. The whole procedure is (quasi) fully automated. We have used our recent algorithm vdW-TSSCDS as implemented on the freely-available AutoMekin software package. To this end, a series of inexpensive classical trajectories using forces from a low-level (semi-empirical) theory are used to sample the configuration space of the system in the search for candidates to first order saddle points. These guess structures are determined by means of a graph-theory based algorithm using the concept of adjacency matrix. Low-level optimizations are followed by re-optimizations at a final high-level of theory (DFT and CCSD(T)-F12 in our case.). The resulting set of stationary points and paths connecting them constitutes the so-called reaction network. In the case of NO2-pyrene, this network exhibits four major basins which can be characterized by their point-group symmetry. A central one, of global C2 symmetry, comprises the global minimum (as well as all other permutationally related conformers) together with the corresponding C2v saddle points connecting them. This central basin is connected to three others of lower C1 symmetry. The latter can be distinguished by the projection of the position of the NO2 nitrogen atom on the pyrene plane in combination with the relative orientation of the oxygen pair pointing either inwards, outwards, upwards or downwards.
Collapse
Affiliation(s)
- Ramón L Panadés-Barrueta
- Faculty of Chemistry and Food Chemistry, Technische Universität Dresden, 01069, Dresden, Germany
| | - Denis Duflot
- Univ. Lille, CNRS, UMR 8523, PhLAM - Physique des Lasers Atomes et Molécules, F-59000, Lille, France
| | - Juan Soto
- Departamento de Química Física, Facultad de Ciencias, Universidad de Málaga, Málaga, Spain
| | - Emilio Martínez-Núñez
- Departamento de Química Física, Facultade de Química, Universidade de Santiago de Compostela, 15782, Santiago de Compostela, Spain
| | - Daniel Peláez
- Université Paris-Saclay, CNRS, Institut des Sciences Moléculaires d'Orsay, 91405, Orsay, France
| |
Collapse
|
10
|
Koda SI, Saito S. Locating Transition States by Variational Reaction Path Optimization with an Energy-Derivative-Free Objective Function. J Chem Theory Comput 2024; 20:2798-2811. [PMID: 38513192 DOI: 10.1021/acs.jctc.3c01246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2024]
Abstract
Locating transition states is essential for understanding molecular reactions. We propose a double-ended transition state search method by revisiting a variational reaction path optimization method known as the MaxFlux method. Although its original purpose is to add temperature effects to reaction paths, we conversely let the temperature approach zero to obtain an asymptotically exact minimum energy path and its corresponding transition state in variational formalism with an energy-derivative-free objective function. Using several numerical techniques to directly optimize the objective function, the present method reliably finds transition states with low computational cost. In particular, only three force evaluations per iteration are sufficient. This is confirmed on a variety of molecular reactions where the nudged elastic band method often fails. The present method is implemented in Python using the Atomic Simulation Environment and is 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
|
11
|
Gong Q, Man Q, Zhao J, Li Y, Dou M, Wang Q, Wu YC, Guo GP. Simulating chemical reaction dynamics on quantum computer. J Chem Phys 2024; 160:124103. [PMID: 38526102 DOI: 10.1063/5.0192036] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2023] [Accepted: 03/06/2024] [Indexed: 03/26/2024] Open
Abstract
The electronic energies of molecules have been successfully evaluated on quantum computers. However, more attention is paid to the dynamics simulation of molecules in practical applications. Based on the variational quantum eigensolver (VQE) algorithm, Fedorov et al. proposed a correlated sampling (CS) method and demonstrated the vibrational dynamics of H2 molecules [J. Chem. Phys. 154, 164103 (2021)]. In this study, we have developed a quantum approach by extending the CS method based on the VQE algorithm (labeled eCS-VQE) for simulating chemical reaction dynamics. First, the CS method is extended to the three-dimensional cases for calculation of first-order energy gradients, and then, it is further generalized to calculate the second-order gradients of energies. By calculating atomic forces and vibrational frequencies for H2, LiH, H+ + H2, and Cl- + CH3Cl systems, we have seen that the approach has achieved the CCSD level of accuracy. Thus, we have simulated dynamics processes for two typical chemical reactions, hydrogen exchange and chlorine substitution, and obtained high-precision reaction dynamics trajectories consistent with the classical methods. Our eCS-VQE approach, as measurement expectations and ground-state wave functions can be reused, is less demanding in quantum computing resources and is, therefore, a feasible means for the dynamics simulation of chemical reactions on the current noisy intermediate-scale quantum-era quantum devices.
Collapse
Affiliation(s)
- Qiankun Gong
- Origin Quantum Computing Company Limited, Hefei, Anhui 230026, China
| | - Qingmin Man
- Origin Quantum Computing Company Limited, Hefei, Anhui 230026, China
| | - Jianyu Zhao
- Origin Quantum Computing Company Limited, Hefei, Anhui 230026, China
| | - Ye Li
- Origin Quantum Computing Company Limited, Hefei, Anhui 230026, China
| | - Menghan Dou
- Origin Quantum Computing Company Limited, Hefei, Anhui 230026, China
| | - Qingchun Wang
- Institute of Artificial Intelligence, Hefei Comprehensive National Science Center, Hefei, Anhui 230088, China
| | - Yu-Chun Wu
- Institute of Artificial Intelligence, Hefei Comprehensive National Science Center, Hefei, Anhui 230088, China
- CAS Key Laboratory of Quantum Information, School of Physics, University of Science and Technology of China, Hefei 230026, China
| | - Guo-Ping Guo
- Origin Quantum Computing Company Limited, Hefei, Anhui 230026, China
- Institute of Artificial Intelligence, Hefei Comprehensive National Science Center, Hefei, Anhui 230088, China
- CAS Key Laboratory of Quantum Information, School of Physics, University of Science and Technology of China, Hefei 230026, China
| |
Collapse
|
12
|
Stocks R, Palethorpe E, Barca GMJ. High-Performance Multi-GPU Analytic RI-MP2 Energy Gradients. J Chem Theory Comput 2024; 20:2505-2519. [PMID: 38456899 DOI: 10.1021/acs.jctc.3c01424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2024]
Abstract
This article presents a novel algorithm for the calculation of analytic energy gradients from second-order Møller-Plesset perturbation theory within the Resolution-of-the-Identity approximation (RI-MP2), which is designed to achieve high performance on clusters with multiple graphical processing units (GPUs). The algorithm uses GPUs for all major steps of the calculation, including integral generation, formation of all required intermediate tensors, solution of the Z-vector equation and gradient accumulation. The implementation in the EXtreme Scale Electronic Structure System (EXESS) software package includes a tailored, highly efficient, multistream scheduling system to hide CPU-GPU data transfer latencies and allows nodes with 8 A100 GPUs to operate at over 80% of theoretical peak floating-point performance. Comparative performance analysis shows a significant reduction in computational time relative to traditional multicore CPU-based methods, with our approach achieving up to a 95-fold speedup over the single-node performance of established software such as Q-Chem and ORCA. Additionally, we demonstrate that pairing our implementation with the molecular fragmentation framework in EXESS can drastically lower the computational scaling of RI-MP2 gradient calculations from quintic to subquadratic, enabling further substantial savings in runtime while retaining high numerical accuracy in the resulting gradients.
Collapse
Affiliation(s)
- Ryan Stocks
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| | - Elise Palethorpe
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| | - Giuseppe M J Barca
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
13
|
Bach RD, Schlegel HB. Mechanism of the Sharpless Epoxidation Reaction: A DFT Study. J Phys Chem A 2024; 128:2072-2091. [PMID: 38452484 DOI: 10.1021/acs.jpca.3c08476] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2024]
Abstract
The Sharpless reaction is an enantioselective epoxidation of prochiral allylic alcohols that employs a Ti(IV) catalyst formed from titanium tetra(isopropoxide), Ti(O-i-Pr)4, diethyl tartrate (DET), and the oxidizing agent tert-butyl hydroperoxide. The M06-2X DFT functional with the 6-311+G(d,p) basis set has been employed to model the structures and energetics of the Sharpless epoxidation reaction. The monomeric tetracoordinate titanium(IV) diethyltartrate is thermodynamically strongly favored to dimerize, producing a pentacoordinate catalyst, [Ti(DET)(O-i-Pr)2]2, that is a more reactive chiral epoxidation catalyst. The rapid ligand exchange reactions needed to generate the "loaded" catalyst and to repeat the overall catalytic cycle have been examined and are found to have activation energies that are much lower than the epoxidation barriers. The transition structures and activation energies for the enantioselective epoxidation of allyl alcohol, trans-methyl-allyl alcohol and trans-tert-butyl-allyl alcohol with the "loaded" Sharpless catalyst, [Ti2(DET)2(O-i-Pr)2-(OAllyl)-(OOt-Bu)], are presented. The effect of the C═O···Ti interactions on the activation energies and the significance of the O-C-C═C dihedral angle on the enantioselectivity of the epoxidation reaction are discussed.
Collapse
Affiliation(s)
- Robert D Bach
- Department of Chemistry and Biochemistry, University of Delaware, Newark, Delaware 19716, United States
| | - H Bernhard Schlegel
- Department of Chemistry, Wayne State University, Detroit, Michigan 48202, United States
| |
Collapse
|
14
|
Pracht P, Grimme S, Bannwarth C, Bohle F, Ehlert S, Feldmann G, Gorges J, Müller M, Neudecker T, Plett C, Spicher S, Steinbach P, Wesołowski PA, Zeller F. CREST-A program for the exploration of low-energy molecular chemical space. J Chem Phys 2024; 160:114110. [PMID: 38511658 DOI: 10.1063/5.0197592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2024] [Accepted: 02/29/2024] [Indexed: 03/22/2024] Open
Abstract
Conformer-rotamer sampling tool (CREST) is an open-source program for the efficient and automated exploration of molecular chemical space. Originally developed in Pracht et al. [Phys. Chem. Chem. Phys. 22, 7169 (2020)] as an automated driver for calculations at the extended tight-binding level (xTB), it offers a variety of molecular- and metadynamics simulations, geometry optimization, and molecular structure analysis capabilities. Implemented algorithms include automated procedures for conformational sampling, explicit solvation studies, the calculation of absolute molecular entropy, and the identification of molecular protonation and deprotonation sites. Calculations are set up to run concurrently, providing efficient single-node parallelization. CREST is designed to require minimal user input and comes with an implementation of the GFNn-xTB Hamiltonians and the GFN-FF force-field. Furthermore, interfaces to any quantum chemistry and force-field software can easily be created. In this article, we present recent developments in the CREST code and show a selection of applications for the most important features of the program. An important novelty is the refactored calculation backend, which provides significant speed-up for sampling of small or medium-sized drug molecules and allows for more sophisticated setups, for example, quantum mechanics/molecular mechanics and minimum energy crossing point calculations.
Collapse
Affiliation(s)
- Philipp Pracht
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Christoph Bannwarth
- Institute for Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056 Aachen, Germany
| | - Fabian Bohle
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Sebastian Ehlert
- AI4Science, Microsoft Research, Evert van de Beekstraat 354, 1118 CZ Schiphol, The Netherlands
| | - Gereon Feldmann
- Institute for Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056 Aachen, Germany
| | - Johannes Gorges
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Marcel Müller
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | - Tim Neudecker
- Institute for Physical and Theoretical Chemistry, University of Bremen, 28359 Bremen, Germany
| | - Christoph Plett
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | | | - Pit Steinbach
- Institute for Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056 Aachen, Germany
| | - Patryk A Wesołowski
- Yusuf Hamied Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Felix Zeller
- Institute for Physical and Theoretical Chemistry, University of Bremen, 28359 Bremen, Germany
| |
Collapse
|
15
|
McFarlane NR, Harvey JN. Exploration of biochemical reactivity with a QM/MM growing string method. Phys Chem Chem Phys 2024; 26:5999-6007. [PMID: 38293892 DOI: 10.1039/d3cp05772k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
In this work, we have implemented the single-ended growing string method using a hybrid internal/Cartesian coordinate scheme within our in-house QM/MM package, QoMMMa, representing the first implementation of the growing string method in the QM/MM framework. The goal of the implementation was to facilitate generation of QM/MM reaction pathways with minimal user input, and also to improve the quality of the pathways generated as compared to the widely used adiabatic mapping approach. We have validated the algorithm against a reaction which has been studied extensively in previous computational investigations - the Claisen rearrangement catalysed by chorismate mutase. The nature of the transition state and the height of the barrier was predicted well using our algorithm, where more than 88% of the pathways generated were deemed to be of production quality. Directly compared to using adiabatic mapping, we found that while our QM/MM single-ended growing string method is slightly less efficient, it readily produces reaction pathways with fewer discontinuites and thus minimises the need for involved remapping of unsatisfactory energy profiles.
Collapse
Affiliation(s)
- Neil R McFarlane
- Department of Chemistry, KU Leuven, B-3001 Leuven, Celestijnenlaan 200f, 2404, Belgium.
| | - Jeremy N Harvey
- Department of Chemistry, KU Leuven, B-3001 Leuven, Celestijnenlaan 200f, 2404, Belgium.
| |
Collapse
|
16
|
Jitaru SC, Enache AC, Cojocaru C, Drochioiu G, Petre BA, Gradinaru VR. Self-Assembly of a Novel Pentapeptide into Hydrogelated Dendritic Architecture: Synthesis, Properties, Molecular Docking and Prospective Applications. Gels 2024; 10:86. [PMID: 38391416 PMCID: PMC10887771 DOI: 10.3390/gels10020086] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2023] [Revised: 01/16/2024] [Accepted: 01/19/2024] [Indexed: 02/24/2024] Open
Abstract
Currently, ultrashort oligopeptides consisting of fewer than eight amino acids represent a cutting-edge frontier in materials science, particularly in the realm of hydrogel formation. By employing solid-phase synthesis with the Fmoc/tBu approach, a novel pentapeptide, FEYNF-NH2, was designed, inspired by a previously studied sequence chosen from hen egg-white lysozyme (FESNF-NH2). Qualitative peptide analysis was based on reverse-phase high performance liquid chromatography (RP-HPLC), while further purification was accomplished using solid-phase extraction (SPE). Exact molecular ion confirmation was achieved by matrix-assisted laser desorption-ionization mass spectrometry (MALDI-ToF MS) using two different matrices (HCCA and DHB). Additionally, the molecular ion of interest was subjected to tandem mass spectrometry (MS/MS) employing collision-induced dissociation (CID) to confirm the synthesized peptide structure. A combination of research techniques, including Fourier-transform infrared spectroscopy (FTIR), fluorescence analysis, transmission electron microscopy, polarized light microscopy, and Congo red staining assay, were carefully employed to glean valuable insights into the self-assembly phenomena and gelation process of the modified FEYNF-NH2 peptide. Furthermore, molecular docking simulations were conducted to deepen our understanding of the mechanisms underlying the pentapeptide's supramolecular assembly formation and intermolecular interactions. Our study provides potential insights into amyloid research and proposes a novel peptide for advancements in materials science. In this regard, in silico studies were performed to explore the FEYNF peptide's ability to form polyplexes.
Collapse
Affiliation(s)
- Stefania-Claudia Jitaru
- Faculty of Chemistry, "Alexandru Ioan Cuza" University, 11 Carol I Bd., 700506 Iasi, Romania
| | - Andra-Cristina Enache
- "Petru Poni" Institute of Macromolecular Chemistry, 41-A Grigore Ghica Voda Alley, 700487 Iasi, Romania
| | - Corneliu Cojocaru
- "Petru Poni" Institute of Macromolecular Chemistry, 41-A Grigore Ghica Voda Alley, 700487 Iasi, Romania
| | - Gabi Drochioiu
- Faculty of Chemistry, "Alexandru Ioan Cuza" University, 11 Carol I Bd., 700506 Iasi, Romania
| | - Brindusa-Alina Petre
- Faculty of Chemistry, "Alexandru Ioan Cuza" University, 11 Carol I Bd., 700506 Iasi, Romania
- TRANSCEND-Regional Institute of Oncology, 700483 Iasi, Romania
| | - Vasile-Robert Gradinaru
- Faculty of Chemistry, "Alexandru Ioan Cuza" University, 11 Carol I Bd., 700506 Iasi, Romania
| |
Collapse
|
17
|
Kim S, Woo J, Kim WY. Diffusion-based generative AI for exploring transition states from 2D molecular graphs. Nat Commun 2024; 15:341. [PMID: 38184661 PMCID: PMC10771475 DOI: 10.1038/s41467-023-44629-6] [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: 05/10/2023] [Accepted: 12/21/2023] [Indexed: 01/08/2024] Open
Abstract
The exploration of transition state (TS) geometries is crucial for elucidating chemical reaction mechanisms and modeling their kinetics. Recently, machine learning (ML) models have shown remarkable performance for prediction of TS geometries. However, they require 3D conformations of reactants and products often with their appropriate orientations as input, which demands substantial efforts and computational cost. Here, we propose a generative approach based on the stochastic diffusion method, namely TSDiff, for prediction of TS geometries just from 2D molecular graphs. TSDiff outperforms the existing ML models with 3D geometries in terms of both accuracy and efficiency. Moreover, it enables to sample various TS conformations, because it learns the distribution of TS geometries for diverse reactions in training. Thus, TSDiff finds more favorable reaction pathways with lower barrier heights than those in the reference database. These results demonstrate that TSDiff shows promising potential for an efficient and reliable TS exploration.
Collapse
Affiliation(s)
- Seonghwan Kim
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, 34141, Daejeon, Republic of Korea
| | - Jeheon Woo
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, 34141, Daejeon, Republic of Korea
| | - Woo Youn Kim
- Department of Chemistry, KAIST, 291 Daehak-ro, Yuseong-gu, 34141, Daejeon, Republic of Korea.
- AI Institute, KAIST, 291 Daehak-ro, Yuseong-gu, 34141, Daejeon, Republic of Korea.
| |
Collapse
|
18
|
Modee R, Mehta S, Laghuvarapu S, Priyakumar UD. MolOpt: Autonomous Molecular Geometry Optimization Using Multiagent Reinforcement Learning. J Phys Chem B 2023; 127:10295-10303. [PMID: 38013420 DOI: 10.1021/acs.jpcb.3c04771] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Most optimization problems require the user to select an algorithm and, to some extent, also tune it for better performance. Although intuition and knowledge about the problem can speed up these selection and fine-tuning processes, users often use trial-and-error methodologies, which can be time-consuming and inefficient. With all of that in mind and much more, the concept of "learned optimizers", "learning to learn", and "meta-learning" has been gathering attention in recent years. In this article, we propose MolOpt that uses multiagent reinforcement learning (MARL) for autonomous molecular geometry optimization (MGO). Typically MGO algorithms are hand-designed, but MolOpt uses MARL to learn a learned optimizer (policy) that can perform MGO without the need for other hand-designed optimizers. We cast MGO as a MARL problem, where each agent corresponds to a single atom in the molecule. MolOpt performs MGO by minimizing the forces on each atom of the molecule. Our experiments demonstrate the generalizing ability of MolOpt for the MGO of propane, pentane, heptane, hexane, and octane when trained on ethane, butane, and isobutane. In terms of performance, MolOpt outperforms the MDMin optimizer and demonstrates performance similar to that of the FIRE optimizer. However, it does not surpass the BFGS optimizer. The results demonstrate that MolOpt has the potential to introduce innovative advancements in MGO by providing a novel approach using reinforcement learning (RL), which may open up new research directions for MGO. Overall, this work serves as a proof-of-concept for the potential of MARL in MGO.
Collapse
Affiliation(s)
- Rohit Modee
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500032, India
| | - Sarvesh Mehta
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500032, India
| | - Siddhartha Laghuvarapu
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500032, India
| | - U Deva Priyakumar
- Center for Computational Natural Sciences and Bioinformatics, International Institute of Information Technology, Hyderabad 500032, India
| |
Collapse
|
19
|
Chang YC, Li YP. Integrating Chemical Information into Reinforcement Learning for Enhanced Molecular Geometry Optimization. J Chem Theory Comput 2023. [PMID: 38012608 DOI: 10.1021/acs.jctc.3c00696] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2023]
Abstract
Geometry optimization is a crucial step in computational chemistry, and the efficiency of optimization algorithms plays a pivotal role in reducing computational costs. In this study, we introduce a novel reinforcement-learning-based optimizer that surpasses traditional methods in terms of efficiency. What sets our model apart is its ability to incorporate chemical information into the optimization process. By exploring different state representations that integrate gradients, displacements, primitive type labels, and additional chemical information from the SchNet model, our reinforcement learning optimizer achieves exceptional results. It demonstrates an average reduction of about 50% or more in optimization steps compared to the conventional optimization algorithms that we examined when dealing with challenging initial geometries. Moreover, the reinforcement learning optimizer exhibits promising transferability across various levels of theory, emphasizing its versatility and potential for enhancing molecular geometry optimization. This research highlights the significance of leveraging reinforcement learning algorithms to harness chemical knowledge, paving the way for future advancements in computational chemistry.
Collapse
Affiliation(s)
- Yu-Cheng Chang
- Department of Chemical Engineering, National Taiwan University, No. 1, Sect. 4, Roosevelt Road, Taipei 10617, Taiwan
| | - Yi-Pei Li
- Department of Chemical Engineering, National Taiwan University, No. 1, Sect. 4, Roosevelt Road, Taipei 10617, Taiwan
- Taiwan International Graduate Program on Sustainable Chemical Science and Technology (TIGP-SCST), Academia Sinica, No. 128, Sec. 2, Academia Road, Taipei 11529, Taiwan
| |
Collapse
|
20
|
Domenichini G, Dellago C. Molecular Hessian matrices from a machine learning random forest regression algorithm. J Chem Phys 2023; 159:194111. [PMID: 37982481 DOI: 10.1063/5.0169384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2023] [Accepted: 10/27/2023] [Indexed: 11/21/2023] Open
Abstract
In this article, we present a machine learning model to obtain fast and accurate estimates of the molecular Hessian matrix. In this model, based on a random forest, the second derivatives of the energy with respect to redundant internal coordinates are learned individually. The internal coordinates together with their specific representation guarantee rotational and translational invariance. The model is trained on a subset of the QM7 dataset but is shown to be applicable to larger molecules picked from the QM9 dataset. From the predicted Hessian, it is also possible to obtain reasonable estimates of the vibrational frequencies, normal modes, and zero point energies of the molecules.
Collapse
Affiliation(s)
- Giorgio Domenichini
- Faculty of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| | - Christoph Dellago
- Faculty of Physics, University of Vienna, Kolingasse 14-16, 1090 Vienna, Austria
| |
Collapse
|
21
|
López-Sosa L, Calaminici P, Köster AM. Cartesian constraints in QM/MM optimizations. J Comput Chem 2023; 44:2358-2368. [PMID: 37635671 DOI: 10.1002/jcc.27202] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/12/2023] [Revised: 07/17/2023] [Accepted: 07/18/2023] [Indexed: 08/29/2023]
Abstract
With the rise of quantum mechanical/molecular mechanical (QM/MM) methods, the interest in the calculation of molecular assemblies has increased considerably. The structures and dynamics of such assemblies are usually governed to a large extend by intermolecular interactions. As a result, the corresponding potential energy surfaces are topological rich and possess many shallow minima. Therefore, local structure optimizations of QM/MM molecular assemblies can be challenging, in particular if optimization constraints are imposed. To overcome this problem, structure optimization in normal coordinate space is advocated. To do so, the external degrees of freedom of a molecule are separated from the internal ones by a projector matrix in the space of the Cartesian coordinates. Here we extend this approach to Cartesian constraints. To this end, we devise an algorithm that adds the Cartesian constraints directly to the projector matrix and in this way eliminates them from the reduced coordinate space in which the molecule is optimized. To analyze the performance and stability of the constrained optimization algorithm in normal coordinate space, we present constrained minimizations of small molecular systems and amino acids in gas phase as well as water employing QM/MM constrained optimizations. All calculations are performed in the framework of auxiliary density functional theory as implemented in the program deMon2k.
Collapse
Affiliation(s)
- L López-Sosa
- Departamento de Química, CINVESTAV, Mexico, Mexico
| | - P Calaminici
- Departamento de Química, CINVESTAV, Mexico, Mexico
| | - A M Köster
- Departamento de Química, CINVESTAV, Mexico, Mexico
| |
Collapse
|
22
|
Shajan A, Manathunga M, Götz AW, Merz KM. Geometry Optimization: A Comparison of Different Open-Source Geometry Optimizers. J Chem Theory Comput 2023; 19:7533-7541. [PMID: 37870541 DOI: 10.1021/acs.jctc.3c00188] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/24/2023]
Abstract
Based on a series of energy minimizations with starting structures obtained from the Baker test set of 30 organic molecules, a comparison is made between various open-source geometry optimization codes that are interfaced with the open-source QUantum Interaction Computational Kernel (QUICK) program for gradient and energy calculations. The findings demonstrate how the choice of the coordinate system influences the optimization process to reach an equilibrium structure. With fewer steps, internal coordinates outperform Cartesian coordinates, while the choice of the initial Hessian and Hessian update method in quasi-Newton approaches made by different optimization algorithms also contributes to the rate of convergence. Furthermore, an available open-source machine learning method based on Gaussian process regression (GPR) was evaluated for energy minimizations over surrogate potential energy surfaces with both Cartesian and internal coordinates with internal coordinates outperforming Cartesian. Overall, geomeTRIC and DL-FIND with their default optimization method as well as with the GPR-based model using Hartree-Fock theory with the 6-31G** basis set needed a comparable number of geometry optimization steps to the approach of Baker using a unit matrix as the initial Hessian to reach the optimized geometry. On the other hand, the Berny and Sella offerings in ASE outperformed the other algorithms. Based on this, we recommend using the file-based approaches, ASE/Berny and ASE/Sella, for large-scale optimization efforts, while if using a single executable is preferable, we now distribute QUICK integrated with DL-FIND.
Collapse
Affiliation(s)
- Akhil Shajan
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
| | - Madushanka Manathunga
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| | - Andreas W Götz
- San Diego Supercomputer Center, University of California San Diego, La Jolla, California 92093-0505, United States
| | - Kenneth M Merz
- Department of Chemistry, Michigan State University, East Lansing, Michigan 48824, United States
- Department of Biochemistry and Molecular Biology, Michigan State University, East Lansing, Michigan 48824, United States
| |
Collapse
|
23
|
Yang WH, Yu FQ, Huang R, Shao GF, Liu TD, Wen YH. Structural Determination and Hierarchical Evolution of Transition Metal Clusters Based on an Improved Self-Adaptive Differential Evolution with Neighborhood Search Algorithm. J Chem Inf Model 2023; 63:6727-6739. [PMID: 37853630 DOI: 10.1021/acs.jcim.3c01331] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2023]
Abstract
Determining the optimal structures and clarifying the corresponding hierarchical evolution of transition metal clusters are of fundamental importance for their applications. The global optimization of clusters containing a large number of atoms, however, is a vastly challenging task encountered in many fields of physics and chemistry. In this work, a high-efficiency self-adaptive differential evolution with neighborhood search (SaNSDE) algorithm, which introduced an optimized cross-operation and an improved Basin Hopping module, was employed to search the lowest-energy structures of CoN, PtN, and FeN (N = 3-200) clusters. The performance of the SaNSDE algorithm was first evaluated by comparing our results with the parallel results collected in the Cambridge Cluster Database (CCD). Subsequently, different analytical methods were introduced to investigate the structural and energetic properties of these clusters systematically, and special attention was paid to elucidating the structural evolution with cluster size by exploring their overall shape, atomic arrangement, structural similarity, and growth pattern. By comparison with those results listed in the CCD, 13 lower-energy structures of FeN clusters were discovered. Moreover, our results reveal that the clusters of three metals had different magic numbers with superior stable structures, most of which possessed high symmetry. The structural evolution of Co, Pt, and Fe clusters could be, respectively, considered as predominantly closed-shell icosahedral, Marks decahedral, and disordered icosahedral-ring growth. Further, the formation of shell structures was discovered, and the clusters with hcp-, fcc-, and bcc-like configurations were ascertained. Nevertheless, the growth of the clusters was not simply atom-to-atom piling up on a given cluster despite gradual saturation of the coordination number toward its bulk limit. Our work identifies the general growth trends for such a wide region of cluster sizes, which would be unbearably expensive in first-principles calculations, and advances the development of global optimization algorithms for the structural prediction of clusters.
Collapse
Affiliation(s)
- Wei-Hua Yang
- Department of Physics, Xiamen University, Xiamen 361005, China
| | - Fang-Qi Yu
- Department of Physics, Xiamen University, Xiamen 361005, China
| | - Rao Huang
- Department of Physics, Xiamen University, Xiamen 361005, China
| | - Gui-Fang Shao
- Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Tun-Dong Liu
- Pen-Tung Sah Institute of Micro-Nano Science and Technology, Xiamen University, Xiamen 361005, China
| | - Yu-Hua Wen
- Department of Physics, Xiamen University, Xiamen 361005, China
| |
Collapse
|
24
|
Talbot JJ, Arias-Martinez JE, Cotton SJ, Head-Gordon M. Fantastical excited state optimized structures and where to find them. J Chem Phys 2023; 159:171102. [PMID: 37916588 DOI: 10.1063/5.0172015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2023] [Accepted: 10/16/2023] [Indexed: 11/03/2023] Open
Abstract
The quantum chemistry community has developed analytic forces for approximate electronic excited states to enable walking on excited state potential energy surfaces (PES). One can thereby computationally characterize excited state minima and saddle points. Always implicit in using this machinery is the fact that an excited state PES only exists within the realm of the Born-Oppenheimer approximation, where the nuclear and electronic degrees of freedom separate. This work demonstrates through ab initio calculations and simple nonadiabatic dynamics that some excited state minimum structures are fantastical: they appear to exist as stable configurations only as a consequence of the PES construct, rather than being physically observable. Each fantastical structure exhibits an unphysically high predicted harmonic frequency and associated force constant. This fact can serve as a valuable diagnostic of when an optimized excited state structure is non-observable. The origin of this phenomenon can be attributed to the coupling between different electronic states. As PESs approach one another, the upper surface can form a minimum that is very close to a near-touching point. The force constant, evaluated at this minimum, relates to the strength of the electronic coupling rather than to any characteristic excited state vibration. Nonadiabatic dynamics results using a Landau-Zener model illustrate that fantastical excited state structures have extremely short lifetimes on the order of a few femtoseconds. Their appearance in a calculation signals the presence of a nearby conical intersection through which the system will rapidly cross to a lower surface.
Collapse
Affiliation(s)
- Justin J Talbot
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
| | - Juan E Arias-Martinez
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Stephen J Cotton
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California Berkeley, Berkeley, California 94720, USA
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
25
|
Li Manni G, Fdez. Galván I, Alavi A, Aleotti F, Aquilante F, Autschbach J, Avagliano D, Baiardi A, Bao JJ, Battaglia S, Birnoschi L, Blanco-González A, Bokarev SI, Broer R, Cacciari R, Calio PB, Carlson RK, Carvalho Couto R, Cerdán L, Chibotaru LF, Chilton NF, Church JR, Conti I, Coriani S, Cuéllar-Zuquin J, Daoud RE, Dattani N, Decleva P, de Graaf C, Delcey M, De Vico L, Dobrautz W, Dong SS, Feng R, Ferré N, Filatov(Gulak) M, Gagliardi L, Garavelli M, González L, Guan Y, Guo M, Hennefarth MR, Hermes MR, Hoyer CE, Huix-Rotllant M, Jaiswal VK, Kaiser A, Kaliakin DS, Khamesian M, King DS, Kochetov V, Krośnicki M, Kumaar AA, Larsson ED, Lehtola S, Lepetit MB, Lischka H, López Ríos P, Lundberg M, Ma D, Mai S, Marquetand P, Merritt ICD, Montorsi F, Mörchen M, Nenov A, Nguyen VHA, Nishimoto Y, Oakley MS, Olivucci M, Oppel M, Padula D, Pandharkar R, Phung QM, Plasser F, Raggi G, Rebolini E, Reiher M, Rivalta I, Roca-Sanjuán D, Romig T, Safari AA, Sánchez-Mansilla A, Sand AM, Schapiro I, Scott TR, Segarra-Martí J, Segatta F, Sergentu DC, Sharma P, Shepard R, Shu Y, Staab JK, Straatsma TP, Sørensen LK, Tenorio BNC, Truhlar DG, Ungur L, Vacher M, Veryazov V, Voß TA, Weser O, Wu D, Yang X, Yarkony D, Zhou C, Zobel JP, Lindh R. The OpenMolcas Web: A Community-Driven Approach to Advancing Computational Chemistry. J Chem Theory Comput 2023; 19:6933-6991. [PMID: 37216210 PMCID: PMC10601490 DOI: 10.1021/acs.jctc.3c00182] [Citation(s) in RCA: 76] [Impact Index Per Article: 76.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2023] [Indexed: 05/24/2023]
Abstract
The developments of the open-source OpenMolcas chemistry software environment since spring 2020 are described, with a focus on novel functionalities accessible in the stable branch of the package or via interfaces with other packages. These developments span a wide range of topics in computational chemistry and are presented in thematic sections: electronic structure theory, electronic spectroscopy simulations, analytic gradients and molecular structure optimizations, ab initio molecular dynamics, and other new features. This report offers an overview of the chemical phenomena and processes OpenMolcas can address, while showing that OpenMolcas is an attractive platform for state-of-the-art atomistic computer simulations.
Collapse
Affiliation(s)
- Giovanni Li Manni
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Ignacio Fdez. Galván
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Ali Alavi
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Yusuf Hamied
Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Flavia Aleotti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Francesco Aquilante
- Theory and
Simulation of Materials (THEOS) and National Centre for Computational
Design and Discovery of Novel Materials (MARVEL), École Polytechnique Fédérale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland
| | - Jochen Autschbach
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
| | - Davide Avagliano
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Alberto Baiardi
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Jie J. Bao
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Stefano Battaglia
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Letitia Birnoschi
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Alejandro Blanco-González
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Sergey I. Bokarev
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
- Chemistry
Department, School of Natural Sciences, Technical University of Munich, Lichtenbergstr. 4, 85748 Garching, Germany
| | - Ria Broer
- Theoretical
Chemistry, Zernike Institute for Advanced Materials, University of Groningen, Nijenborgh 4, 9747AG Groningen, The Netherlands
| | - Roberto Cacciari
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Paul B. Calio
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Rebecca K. Carlson
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Rafael Carvalho Couto
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luis Cerdán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
- Instituto
de Óptica (IO−CSIC), Consejo
Superior de Investigaciones Científicas, 28006, Madrid, Spain
| | - Liviu F. Chibotaru
- Department
of Chemistry, KU Leuven, Celestijnenlaan 200F, 3001 Leuven, Belgium
| | - Nicholas F. Chilton
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | | | - Irene Conti
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Sonia Coriani
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Juliana Cuéllar-Zuquin
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Razan E. Daoud
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Nike Dattani
- HPQC Labs, Waterloo, N2T 2K9 Ontario Canada
- HPQC College, Waterloo, N2T 2K9 Ontario Canada
| | - Piero Decleva
- Istituto
Officina dei Materiali IOM-CNR and Dipartimento di Scienze Chimiche
e Farmaceutiche, Università degli
Studi di Trieste, I-34121 Trieste, Italy
| | - Coen de Graaf
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
- ICREA, Pg. Lluís
Companys 23, 08010 Barcelona, Spain
| | - Mickaël
G. Delcey
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Luca De Vico
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Werner Dobrautz
- Chalmers
University of Technology, Department of Chemistry
and Chemical Engineering, 41296 Gothenburg, Sweden
| | - Sijia S. Dong
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Chemical Biology, Department of Physics, and Department
of Chemical Engineering, Northeastern University, Boston, Massachusetts 02115, United States
| | - Rulin Feng
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Department
of Chemistry, Fudan University, Shanghai 200433, China
| | - Nicolas Ferré
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | | | - Laura Gagliardi
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Marco Garavelli
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Leticia González
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Yafu Guan
- State Key
Laboratory of Molecular Reaction Dynamics and Center for Theoretical
Computational Chemistry, Dalian Institute
of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, People’s Republic of China
| | - Meiyuan Guo
- SSRL, SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Matthew R. Hennefarth
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Chad E. Hoyer
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Miquel Huix-Rotllant
- Institut
de Chimie Radicalaire (UMR-7273), Aix-Marseille
Univ, CNRS, ICR 13013 Marseille, France
| | - Vishal Kumar Jaiswal
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Andy Kaiser
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Danil S. Kaliakin
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - Marjan Khamesian
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Daniel S. King
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Vladislav Kochetov
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Marek Krośnicki
- Institute
of Theoretical Physics and Astrophysics, Faculty of Mathematics, Physics
and Informatics, University of Gdańsk, ul Wita Stwosza 57, 80-952, Gdańsk, Poland
| | | | - Ernst D. Larsson
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Susi Lehtola
- Molecular
Sciences Software Institute, Blacksburg, Virginia 24061, United States
- Department
of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 University of Helsinki, Finland
| | - Marie-Bernadette Lepetit
- Condensed
Matter Theory Group, Institut Néel, CNRS UPR 2940, 38042 Grenoble, France
- Theory
Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409-1061, United States
| | - Pablo López Ríos
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Marcus Lundberg
- Department
of Chemistry − Ångström Laboratory, Uppsala University, SE-75120 Uppsala, Sweden
| | - Dongxia Ma
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Sebastian Mai
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Philipp Marquetand
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | | | - Francesco Montorsi
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Maximilian Mörchen
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Artur Nenov
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Vu Ha Anh Nguyen
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Yoshio Nishimoto
- Graduate
School of Science, Kyoto University, Kyoto 606-8502, Japan
| | - Meagan S. Oakley
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Massimo Olivucci
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Markus Oppel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Daniele Padula
- Dipartimento
di Biotecnologie, Chimica e Farmacia, Università
di Siena, Via A. Moro 2, 53100 Siena, Italy
| | - Riddhish Pandharkar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
| | - Quan Manh Phung
- Department
of Chemistry, Graduate School of Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8602, Japan
- Institute
of Transformative Bio-Molecules (WPI-ITbM), Nagoya University, Furo-cho, Chikusa-ku, Nagoya, Aichi 464-8601, Japan
| | - Felix Plasser
- Department
of Chemistry, Loughborough University, Loughborough, LE11 3TU, U.K.
| | - Gerardo Raggi
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Quantum
Materials and Software LTD, 128 City Road, London, EC1V 2NX, United Kingdom
| | - Elisa Rebolini
- Scientific
Computing Group, Institut Laue Langevin, 38042 Grenoble, France
| | - Markus Reiher
- ETH Zurich, Laboratory for Physical Chemistry, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Ivan Rivalta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Daniel Roca-Sanjuán
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Thies Romig
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Arta Anushirwan Safari
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Aitor Sánchez-Mansilla
- Department
of Physical and Inorganic Chemistry, Universitat
Rovira i Virgili, Tarragona 43007, Spain
| | - Andrew M. Sand
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry and Biochemistry, Butler University, Indianapolis, Indiana 46208, United States
| | - Igor Schapiro
- Institute
of Chemistry, The Hebrew University of Jerusalem, Jerusalem 91904, Israel
| | - Thais R. Scott
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
- Department
of Chemistry, Pritzker School of Molecular Engineering, James Franck
Institute, Chicago Center for Theoretical Chemistry, The University of Chicago, Chicago, Illinois 60637, United States
- Department
of Chemistry, University of California, Irvine, California 92697, United States
| | - Javier Segarra-Martí
- Instituto
de Ciencia Molecular, Universitat de València, Catedrático José Beltrán
Martínez n. 2, 46980 Paterna, Spain
| | - Francesco Segatta
- Department
of Industrial Chemistry “Toso Montanari”, University of Bologna, 40136 Bologna, Italy
| | - Dumitru-Claudiu Sergentu
- Department
of Chemistry, University at Buffalo, State
University of New York, Buffalo, New York 14260-3000, United States
- Laboratory
RA-03, RECENT AIR, A. I. Cuza University of Iaşi, RA-03 Laboratory (RECENT AIR), Iaşi 700506, Romania
| | - Prachi Sharma
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, USA
| | - Yinan Shu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Jakob K. Staab
- The Department
of Chemistry, The University of Manchester, M13 9PL, Manchester, U.K.
| | - Tjerk P. Straatsma
- National
Center for Computational Sciences, Oak Ridge
National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
- Department
of Chemistry and Biochemistry, University
of Alabama, Tuscaloosa, Alabama 35487-0336, United States
| | | | - Bruno Nunes Cabral Tenorio
- Department
of Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, 2800 Kongens Lyngby, Denmark
| | - Donald G. Truhlar
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Liviu Ungur
- Department
of Chemistry, National University of Singapore, 3 Science Drive 3, 117543 Singapore
| | - Morgane Vacher
- Nantes
Université, CNRS, CEISAM, UMR 6230, F-44000 Nantes, France
| | - Valera Veryazov
- Division
of Theoretical Chemistry, Chemical Centre, Lund University, P.O. Box 124, SE-22100, Lund, Sweden
| | - Torben Arne Voß
- Institut
für Physik, Universität Rostock, Albert-Einstein-Str. 23-24, 18059 Rostock, Germany
| | - Oskar Weser
- Electronic
Structure Theory Department, Max Planck
Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - Dihua Wu
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - Xuchun Yang
- Chemistry
Department, Bowling Green State University, Overmann Hall, Bowling Green, Ohio 43403, United States
| | - David Yarkony
- Department
of Chemistry, Johns Hopkins University, Baltimore, Maryland 21218, United States
| | - Chen Zhou
- Department
of Chemistry, Chemical Theory Center, and Minnesota Supercomputing
Institute, University of Minnesota, Minneapolis, Minnesota 55455-0431, United
States
| | - J. Patrick Zobel
- Institute
of Theoretical Chemistry, Faculty of Chemistry, University of Vienna, Währinger Straße 17, A-1090 Vienna, Austria
| | - Roland Lindh
- Department
of Chemistry − BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala
Center for Computational Chemistry (UC3), Uppsala University, PO Box 576, SE-751 23 Uppsala. Sweden
| |
Collapse
|
26
|
Tripathy V, Raghavachari K. Fragment-based models for dissociation of strong acids in water: Electrostatic embedding minimizes the dependence on the fragmentation schemes. J Chem Phys 2023; 159:124106. [PMID: 38127382 DOI: 10.1063/5.0164089] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Accepted: 08/28/2023] [Indexed: 12/23/2023] Open
Abstract
Fragmentation methods such as MIM (Molecules-in-Molecules) provide a route to accurately model large systems and have been successful in predicting their structures, energies, and spectroscopic properties. However, their use is often limited to systems at equilibrium due to the inherent complications in the choice of fragments in systems away from equilibrium. Furthermore, the presence of charges resulting from any heterolytic bond breaking may increase the fragmentation error. We have previously suggested EE-MIM (Electrostatically Embedded Molecules-In-Molecules) as a method to mitigate the errors resulting from the missing long-range interactions in molecular clusters in equilibrium. Here, we show that the same method can be applied to improve the performance of MIM to solve the longstanding problem of dependency of the fragmentation energy error on the choice of the fragmentation scheme. We chose four widely used acid dissociation reactions (HCl, HClO4, HNO3, and H2SO4) as test cases due to their importance in chemical processes and complex reaction potential energy surfaces. Electrostatic embedding improves the performance at both one and two-layer MIM as shown by lower EE-MIM1 and EE-MIM2 errors. The EE-MIM errors are also demonstrated to be less dependent on the choice of the fragmentation scheme by analyzing the variation in fragmentation energy at the points with more than one possible fragmentation scheme (points where the fragmentation scheme changes). EE-MIM2 with M06-2X as the low-level resulted in a variation of less than 1 kcal/mol for all the cases and 1 kJ/mol for all but three cases, rendering our method fragmentation scheme-independent for acid dissociation processes.
Collapse
Affiliation(s)
- Vikrant Tripathy
- Department of Chemistry, Indiana University, Bloomington, Indiana 47405, USA
| | | |
Collapse
|
27
|
Kaluarachchige Don UI, Palmer Z, Ward CL, Lord RL, Groysman S. Combining [Mo VIO 3] and [M 0(CO) 3] (M = Mo, Cr) Fragments within the Same Complex: Synthesis and Reactivity of the Single Oxo-Bridged Heterobimetallics Supported by Xanthene-Based Heterodinucleating Ligands. Inorg Chem 2023; 62:15063-15075. [PMID: 37677846 DOI: 10.1021/acs.inorgchem.3c01929] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/09/2023]
Abstract
A functional model of Mo-Cu carbon monoxide dehydrogenase (CODH) enzyme requires the presence of an oxidant (metal-oxo) and a metal-bound carbonyl in close proximity. In this work, we report the synthesis, characterization, and reactivity of a heterobimetallic complex combining Mo(VI) trioxo with Mo(0) tricarbonyl. The formation of the heterobimetallic complex is facilitated by the xanthene-bridged heterodinucleating ligand containing a hard catecholate chelate and a soft iminopyridine chelate. A catechol-coordinated square-pyramidal [MoVIO3] fragment interacts directly with the iminopyridine-bound [Mo0(CO)3] fragment via a single (oxo) bridge, with the overall disposition being related to the proposed first step in the CODH mechanism, where square-pyramidal [MoVIO2S] interacts with the [Cu-CO] via a single sulfido bridge. Our attempt to obtain a sulfido-bridged analogue (using [MoO3S]2- precursor) led to a mixture of products possibly containing different (oxo and sulfido) bridges. Despite a direct interaction between Mo(VI) and Mo(0) segments, no internal redox is observed, with the high lying occupied MOs being mostly d-π orbitals at Mo0(CO)3 and the low lying unoccupied MOs being d-π orbitals at MoVIO3. Due to the overall rigid structure, the heterobimetallic complex was found to be stable up to 100 °C in DMF-d7 (based on 1H NMR). The decomposition of the complex above this temperature does not produce CO2 (based on gas chromatography), dissociating stable Mo(CO)3(DMF)3 instead (based on IR). We also synthesized and studied the reactivity of the Mo(VI)/Cr(0) analogue. While this complex demonstrated more facile decomposition, no CO2 production was observed. Density functional theory calculations suggest that the formation of [CO2]2- and its subsequent reductive elimination is endergonic in the present system, likely due to the stability of fac-Mo0(CO)3 and the relative nucleophilic character of the carbonyl carbon engendered by back donation from Mo(0). The calculations also indicate that the replacement of one oxo by sulfido (both terminal and bridging), replacement of catechol with dithiolene, and replacement of Mo(0) with Cr(0) does not affect significantly the energetics of the process, likely requiring the use a less stable and less π-basic CO anchor.
Collapse
Affiliation(s)
| | - Zsolt Palmer
- Department of Chemistry, Grand Valley State University, 1 Campus Drive, Allendale, Michigan 49401, United States
| | - Cassandra L Ward
- Lumigen Instrument Center, Wayne State University, 5101 Cass Avenue, Detroit, Michigan 48202, United States
| | - Richard L Lord
- Department of Chemistry, Grand Valley State University, 1 Campus Drive, Allendale, Michigan 49401, United States
| | - Stanislav Groysman
- Department of Chemistry, Wayne State University, 5101 Cass Ave. Detroit, Michigan 48202, United States
| |
Collapse
|
28
|
Fdez Galván I, Lindh R. Smooth Things Come in Threes: A Diabatic Surrogate Model for Conical Intersection Optimization. J Chem Theory Comput 2023. [PMID: 37192531 DOI: 10.1021/acs.jctc.3c00389] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/18/2023]
Abstract
The optimization of conical intersection structures is complicated by the nondifferentiability of the adiabatic potential energy surfaces. In this work, we build a pseudodiabatic surrogate model, based on Gaussian process regression, formed by three smooth and differentiable surfaces that can adequately reproduce the adiabatic surfaces. Using this model with the restricted variance optimization method results in a notable decrease of the overall computational effort required to obtain minimum energy crossing points.
Collapse
Affiliation(s)
- Ignacio Fdez Galván
- Department of Chemistry-BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| | - Roland Lindh
- Department of Chemistry-BMC, Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
- Uppsala Center for Computational Chemistry (UC3), Uppsala University, P.O. Box 576, SE-75123 Uppsala, Sweden
| |
Collapse
|
29
|
R S, Mahalakshmi S, Vetrivelan V, Irfan A, Muthu S. Absorption wavelength (TD-DFT) and adsorption of metal chalcogen clusters with methyl nicotinate: Structural, electronic, IRI, SERS, pharmacological and antiviral studies (HIV and omicron). Heliyon 2023; 9:e16066. [PMID: 37234664 PMCID: PMC10208831 DOI: 10.1016/j.heliyon.2023.e16066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Revised: 04/27/2023] [Accepted: 05/04/2023] [Indexed: 05/28/2023] Open
Abstract
The DFT B3LYP-LAND2DZ technique is used to examine interactions of Methyl nicotinate with copper selenide and zinc selenide clusters. The existence of reactive sites is determined using ESP maps and Fukui data. The energy variations between HOMO and LUMO are utilised to calculate various energy parameters. The Atoms in Molecules and ELF (Electron Localisation Function) maps are employed to investigate the topology of the molecule. The Interaction Region Indicator is used to determine the existence of non-covalent zones in the molecule. The UV-Vis spectrum using the TD-DFT method and DOS graphs are used to obtain the theoretical determination of electronic transition and properties. Structural analysis of the compound is obtained using theoretical IR spectra. To explore the adsorption of copper selenide and zinc selenide clusters on the Methyl nicotinate, the adsorption energy and theoretical SERS spectra are employed. Furthermore, pharmacological investigations are carried out to confirm the drug's non-toxicity. The compound's antiviral efficacy against HIV and Omicron is demonstrated via protein-ligand docking.
Collapse
Affiliation(s)
- Sravanthi R
- Department of Physics, Ethiraj College for Women, Chennai, 600008, Tamil Nadu, India
- University of Madras, Chennai, 600005, Tamil Nadu, India
| | - S. Mahalakshmi
- Department of Physics, Ethiraj College for Women, Chennai, 600008, Tamil Nadu, India
| | - V. Vetrivelan
- Department of Physics, Government College of Engineering, Srirangam, Trichy 620012, Tamil Nadu, India
| | - Ahmad Irfan
- Department of Chemistry, College of Science, King Khalid University, P. O. Box 9004, Abha, 61413, Saudi Arabia
| | - S. Muthu
- Department of Physics, Arignar Anna Govt.Arts College, Cheyyar, 604407, Tamil Nadu, India
| |
Collapse
|
30
|
Eaby AC, Myburgh DC, Kosimov A, Kwit M, Esterhuysen C, Janiak AM, Barbour LJ. Dehydration of a crystal hydrate at subglacial temperatures. Nature 2023; 616:288-292. [PMID: 37045922 PMCID: PMC10097597 DOI: 10.1038/s41586-023-05749-7] [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: 04/04/2022] [Accepted: 01/23/2023] [Indexed: 04/14/2023]
Abstract
Water is one of the most important substances on our planet1. It is ubiquitous in its solid, liquid and vaporous states and all known biological systems depend on its unique chemical and physical properties. Moreover, many materials exist as water adducts, chief among which are crystal hydrates (a specific class of inclusion compound), which usually retain water indefinitely at subambient temperatures2. We describe a porous organic crystal that readily and reversibly adsorbs water into 1-nm-wide channels at more than 55% relative humidity. The water uptake/release is chromogenic, thus providing a convenient visual indication of the hydration state of the crystal over a wide temperature range. The complementary techniques of X-ray diffraction, optical microscopy, differential scanning calorimetry and molecular simulations were used to establish that the nanoconfined water is in a state of flux above -70 °C, thus allowing low-temperature dehydration to occur. We were able to determine the kinetics of dehydration over a wide temperature range, including well below 0 °C which, owing to the presence of atmospheric moisture, is usually challenging to accomplish. This discovery unlocks opportunities for designing materials that capture/release water over a range of temperatures that extend well below the freezing point of bulk water.
Collapse
Affiliation(s)
- Alan C Eaby
- Department of Chemistry and Polymer Science, Stellenbosch University, Stellenbosch, South Africa
| | - Dirkie C Myburgh
- Department of Chemistry and Polymer Science, Stellenbosch University, Stellenbosch, South Africa
| | - Akmal Kosimov
- Faculty of Chemistry, Adam Mickiewicz University, Poznań, Poland
| | - Marcin Kwit
- Faculty of Chemistry, Adam Mickiewicz University, Poznań, Poland
| | - Catharine Esterhuysen
- Department of Chemistry and Polymer Science, Stellenbosch University, Stellenbosch, South Africa.
| | | | - Leonard J Barbour
- Department of Chemistry and Polymer Science, Stellenbosch University, Stellenbosch, South Africa.
| |
Collapse
|
31
|
Choi S. Prediction of transition state structures of gas-phase chemical reactions via machine learning. Nat Commun 2023; 14:1168. [PMID: 36859495 PMCID: PMC9977841 DOI: 10.1038/s41467-023-36823-3] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2022] [Accepted: 02/15/2023] [Indexed: 03/03/2023] Open
Abstract
The elucidation of transition state (TS) structures is essential for understanding the mechanisms of chemical reactions and exploring reaction networks. Despite significant advances in computational approaches, TS searching remains a challenging problem owing to the difficulty of constructing an initial structure and heavy computational costs. In this paper, a machine learning (ML) model for predicting the TS structures of general organic reactions is proposed. The proposed model derives the interatomic distances of a TS structure from atomic pair features reflecting reactant, product, and linearly interpolated structures. The model exhibits excellent accuracy, particularly for atomic pairs in which bond formation or breakage occurs. The predicted TS structures yield a high success ratio (93.8%) for quantum chemical saddle point optimizations, and 88.8% of the optimization results have energy errors of less than 0.1 kcal mol-1. Additionally, as a proof of concept, the exploration of multiple reaction paths of an organic reaction is demonstrated based on ML inferences. I envision that the proposed approach will aid in the construction of initial geometries for TS optimization and reaction path exploration.
Collapse
Affiliation(s)
- Sunghwan Choi
- Division of National Supercomputing, Korea Institute of Science and Technology Information, 245 Daehak-ro, Yuseong-gu, 34141, Daejeon, Republic of Korea.
| |
Collapse
|
32
|
Xie Z, Song Y, Peng F, Li J, Cheng Y, Zhang L, Ma Y, Tian Y, Luo Z, Ma H. Kylin 1.0: An ab-initio density matrix renormalization group quantum chemistry program. J Comput Chem 2023; 44:1316-1328. [PMID: 36809661 DOI: 10.1002/jcc.27085] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/28/2022] [Revised: 01/15/2023] [Accepted: 01/19/2023] [Indexed: 02/23/2023]
Abstract
The accurate evaluation of electron correlations is highly necessary for the proper descriptions of the electronic structures in strongly correlated molecules, ranging from bond-dissociating molecules, polyradicals, to large conjugated molecules and transition metal complexes. For this purpose, in this paper, a new ab-initio quantum chemistry program Kylin 1.0 for electron correlation calculations at various quantum many-body levels, including configuration interaction (CI), perturbation theory (PT), and density matrix renormalization group (DMRG), is presented. Furthermore, fundamental quantum chemistry methods such as Hartree-Fock self-consistent field (HF-SCF) and the complete active space SCF (CASSCF) are also implemented. The Kylin 1.0 program possesses the following features: (1) a matrix product operator (MPO) formulation-based efficient DMRG implementation for describing static electron correlation within a large active space composed of more than 100 orbitals, supporting both U 1 n × U 1 S z $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{U}{(1)}_{{\mathrm{S}}_{\mathrm{z}}} $$ and U 1 n × SU 2 S $$ \mathrm{U}{(1)}_{\mathrm{n}}\times \mathrm{SU}{(2)}_{\mathrm{S}} $$ symmetries; (2) an efficient second-order DMRG-self-consistent field (SCF) implementation; (3) an externally contracted multi-reference CI (MRCI) and Epstein-Nesbet PT with DMRG reference wave functions for including the remaining dynamic electron correlation outside the large active spaces. In this paper, we introduce the capabilities and numerical benchmark examples of the Kylin 1.0 program.
Collapse
Affiliation(s)
- Zhaoxuan Xie
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yinxuan Song
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Fangwen Peng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Jianhao Li
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yifan Cheng
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Lingzhi Zhang
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Yingjin Ma
- Computer Network Information Center, Chinese Academy of Sciences, Beijing, China
| | - Yingqi Tian
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Zhen Luo
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China
| | - Haibo Ma
- School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, China.,Qingdao Institute for Theoretical and Computational Sciences, Qingdao Institute of Frontier and Interdisciplinary Science, Shandong University, Qingdao, China
| |
Collapse
|
33
|
Weser O, Hein-Janke B, Mata RA. Automated handling of complex chemical structures in Z-matrix coordinates-The chemcoord library. J Comput Chem 2023; 44:710-726. [PMID: 36541725 DOI: 10.1002/jcc.27029] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/11/2022] [Revised: 10/05/2022] [Accepted: 10/09/2022] [Indexed: 12/24/2022]
Abstract
In this work, we present a fully automated method for the construction of chemically meaningful sets of hierarchical nonredundant internal coordinates (ICs; also commonly denoted as Z-matrices) from the Cartesian coordinates of a molecular system. Particular focus is placed on avoiding ill-definitions of angles and dihedrals due to linear arrangements of atoms, to consistently guarantee a well-defined transformation to Cartesian coordinates, even after structural changes. The representations thus obtained are particularly well suited for pathway construction in double-ended methods for transition state search and optimizations with nonlinear constraints. Analytical gradients for the transformation between the coordinate systems were derived for analytical geometry optimizations purely in Z-matrix coordinates. The geometry optimization was coupled with a Symbolic Algebra package to support arbitrary nonlinear constraints in Z-matrix coordinates, while retaining analytical energy gradient conversion. The difference to the commonly used nonhierarchical IC transformations is discussed. Sample applications are provided for a number of common chemical reactions and illustrative examples.
Collapse
Affiliation(s)
- Oskar Weser
- Electronic Structure Theory Department, Max-Planck-Institute for Solid State Research, Stuttgart, Germany.,Institute of Physical Chemistry, University of Goettingen, Goettingen, Germany
| | - Björn Hein-Janke
- Institute of Physical Chemistry, University of Goettingen, Goettingen, Germany
| | - Ricardo A Mata
- Institute of Physical Chemistry, University of Goettingen, Goettingen, Germany
| |
Collapse
|
34
|
Stevens JE, Pefley CM, Piatkowski A, Smith ZR, Ognanovich N. Density functional theory investigation of mechanisms of degradation reactions of sulfonated PEEK membranes with OH radicals in fuel cells: Addition-elimination reactions and acid catalyzed water elimination. RESEARCH SQUARE 2023:rs.3.rs-2565467. [PMID: 36798331 PMCID: PMC9934756 DOI: 10.21203/rs.3.rs-2565467/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/12/2023]
Abstract
Sulfonated polyether (ether) ketone, or sulfonated PEEK (sPEEK) membranes are one possible candidate for proton-transfer membranes in hydrogen fuel cells. Reaction with hydroxy radicals is expected to be a significant source of degradation of these membranes during fuel cell operation. In this work, the reactivity of the sPEEK polymer molecule with OH radicals is studied by M062X hybrid density functional calculations of the energetics of several reaction paths in a water environment as modeled by polarized continuum model (PCM) calculations. Reactants, products, encounter minima and transition states are optimized for a reaction pathway in which OH addition is followed by acid-catalyzed water elimination which cationizes the polymer, degradation is expected to follow this reaction as the unstable cation then undergoes bond-breaking or other reactions. Two pathways for this acid-catalyzed cationization, one in which a water molecule plays the role of an additional co-catalyst, are reported. Further calculations explore reaction pathways in which addition of OH to the polymer is followed by bond breaking reactions which would break the polymer chain or the bond between the polymer and sulfonyl groups. Examination of the free energy barriers to all these reactions, relative to reactants, suggest that these direct bond-breaking reactions may compete somewhat with acid-catalyzed water elimination following OH addition.
Collapse
|
35
|
Stevens JE, Pefley CM, Piatkowski A, Smith ZR, Ognanovich N. Density functional theory investigation of mechanisms of degradation reactions of sulfonated PEEK membranes with OH radicals in fuel cells: addition-elimination reactions and acid catalyzed water elimination. Theor Chem Acc 2023; 142:49. [PMID: 37124478 PMCID: PMC10129967 DOI: 10.1007/s00214-023-02981-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/08/2023] [Accepted: 03/30/2023] [Indexed: 05/02/2023]
Abstract
Sulfonated polyether (ether) ketone or sulfonated PEEK (sPEEK) membranes are one possible candidate for proton-transfer membranes in hydrogen fuel cells. Reaction with hydroxy radicals is expected to be a significant source of degradation of these membranes during fuel cell operation. In this work, the reactivity of the sPEEK polymer molecule with OH radicals is studied by M062X hybrid density functional calculations of the energetics of several reaction paths in a water environment as modeled by polarized continuum model calculations. Reactants, products, encounter minima and transition states are optimized for a reaction pathway in which OH addition is followed by acid-catalyzed water elimination which cationizes the polymer, degradation is expected to follow this reaction as the unstable cation then undergoes bond-breaking or other reactions. Two pathways for this acid-catalyzed cationization, one in which a water molecule plays the role of an additional co-catalyst, are reported. Further calculations explore reaction pathways in which addition of OH to the polymer is followed by bond breaking reactions which would break the polymer chain or the bond between the polymer and sulfonyl groups. Examination of the free energy barriers to all these reactions, relative to reactants, suggests that these direct bond-breaking reactions may compete somewhat with acid-catalyzed water elimination following OH addition. Supplementary Information The online version contains supplementary material available at 10.1007/s00214-023-02981-2.
Collapse
Affiliation(s)
- Jonathan E. Stevens
- Department of Chemistry and Biochemistry, University of Detroit Mercy, Detroit, MI 48221 USA
| | - Courtney M. Pefley
- Department of Chemistry and Biochemistry, University of Detroit Mercy, Detroit, MI 48221 USA
| | - Alice Piatkowski
- School of Dentistry, University of Detroit Mercy, Detroit, MI 48208 USA
| | | | | |
Collapse
|
36
|
Sugisaki K, Wakimoto H, Toyota K, Sato K, Shiomi D, Takui T. Quantum Algorithm for Numerical Energy Gradient Calculations at the Full Configuration Interaction Level of Theory. J Phys Chem Lett 2022; 13:11105-11111. [PMID: 36444985 PMCID: PMC9743205 DOI: 10.1021/acs.jpclett.2c02737] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2022] [Accepted: 11/07/2022] [Indexed: 06/16/2023]
Abstract
A Bayesian phase difference estimation (BPDE) algorithm allows us to compute the energy gap of two electronic states of a given Hamiltonian directly by utilizing the quantum superposition of their wave functions. Here we report an extension of the BPDE algorithm to the direct calculation of the energy difference of two molecular geometries. We apply the BPDE algorithm for the calculation of numerical energy gradients based on the two-point finite-difference method, enabling us to execute geometry optimization of one-dimensional molecules at the full-CI level on a quantum computer. Results of numerical quantum circuit simulations of the geometry optimization of the H2 molecule with the STO-3G and 6-31G basis sets, the LiH and BeH2 molecules at the full-CI/STO-3G level, and the N2 molecule at the CASCI(6e,6o)/6-311G* level are given.
Collapse
Affiliation(s)
- Kenji Sugisaki
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
- JSTPRESTO, 4-1-8 Honcho, Kawaguchi, Saitama, 332-0012, Japan
- Centre
for Quantum Engineering, Research and Education (CQuERE), TCG Centres for Research and Education in Science
and Technology (TCG CREST), Sector V,
Salt Lake, Kolkata700091, India
| | - Hiroyuki Wakimoto
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Kazuo Toyota
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Kazunobu Sato
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Daisuke Shiomi
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
| | - Takeji Takui
- Department
of Chemistry, Graduate School of Science, Osaka Metropolitan University, 3-3-138 Sugimoto, Sumiyoshi-ku, Osaka558-8585, Japan
- Research
Support Department/University Research Administrator Center, University
Administration Division, Osaka Metropolitan
University, 3-3-138 Sugimoto,
Sumiyoshi-ku, Osaka558-8585, Japan
| |
Collapse
|
37
|
Ouahdi Z, Ourhriss N, Aitouna AO, Barhoumi A, Belghiti ME, Moubarik A, El Alaoui El Abdallaoui H, El Idrissi M, Zeroual A. Exploration of the mechanism, chemospecificity, regiospecificity and stereoselectivity of the cycloaddition reaction between 9α-hydroxyparthenolide and nitrilimine: MEDT study. Theor Chem Acc 2022. [DOI: 10.1007/s00214-022-02913-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
|
38
|
Pracht P, Bannwarth C. Fast Screening of Minimum Energy Crossing Points with Semiempirical Tight-Binding Methods. J Chem Theory Comput 2022; 18:6370-6385. [PMID: 36121838 DOI: 10.1021/acs.jctc.2c00578] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The investigation of photochemical processes is a highly active field in computational chemistry. One research direction is the automated exploration and identification of minimum energy conical intersection (MECI) geometries. However, due to the immense technical effort required to calculate nonadiabatic potential energy landscapes, the routine application of such computational protocols is severely limited. In this study, we will discuss the prospect of combining adiabatic potential energy surfaces from semiempirical quantum mechanical calculations with specialized confinement potential and metadynamics simulations to identify S0/T1 minimum energy crossing point (MECP) geometries. It is shown that MECPs calculated at the GFN2-xTB level can provide suitable approximations to high-level S0/S1ab initio conical intersection geometries at a fraction of the computational cost. Reference MECIs of benzene are studied to illustrate the basic concept. An example application of the presented protocol is demonstrated for a set of photoswitch molecules.
Collapse
Affiliation(s)
- Philipp Pracht
- Institute of Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056Aachen, Germany
| | - Christoph Bannwarth
- Institute of Physical Chemistry, RWTH Aachen University, Melatener Str. 20, 52056Aachen, Germany
| |
Collapse
|
39
|
Aziridination Reactivity of a Manganese(II) Complex with a Bulky Chelating Bis(Alkoxide) Ligand. Molecules 2022; 27:molecules27185751. [PMID: 36144492 PMCID: PMC9505844 DOI: 10.3390/molecules27185751] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2022] [Revised: 08/30/2022] [Accepted: 08/30/2022] [Indexed: 11/17/2022] Open
Abstract
Treatment of Mn(N(SiMe3)2)2(THF)2 with bulky chelating bis(alkoxide) ligand [1,1′:4′,1′′-terphenyl]-2,2′′-diylbis(diphenylmethanol) (H2[O-terphenyl-O]Ph) formed a seesaw manganese(II) complex Mn[O-terphenyl-O]Ph(THF)2, characterized by structural, spectroscopic, magnetic, and analytical methods. The reactivity of Mn[O-terphenyl-O]Ph(THF)2 with various nitrene precursors was investigated. No reaction was observed between Mn[O-terphenyl-O]Ph(THF)2 and aryl azides. In contrast, the treatment of Mn[O-terphenyl-O]Ph(THF)2 with iminoiodinane PhINTs (Ts = p-toluenesulfonyl) was consistent with the formation of a metal-nitrene complex. In the presence of styrene, the reaction led to the formation of aziridine. Combining varying ratios of styrene and PhINTs in different solvents with 10 mol% of Mn[O-terphenyl-O]Ph(THF)2 at room temperature produced 2-phenylaziridine in up to a 79% yield. Exploration of the reactivity of Mn[O-terphenyl-O]Ph(THF)2 with various olefins revealed (1) moderate aziridination yields for p-substituted styrenes, irrespective of the electronic nature of the substituent; (2) moderate yield for 1,1′-disubstituted α-methylstyrene; (3) no aziridination for aliphatic α-olefins; (4) complex product mixtures for the β-substituted styrenes. DFT calculations suggest that iminoiodinane is oxidatively added upon binding to Mn, and the resulting formal imido intermediate has a high-spin Mn(III) center antiferromagnetically coupled to an imidyl radical. This imidyl radical reacts with styrene to form a sextet intermediate that readily reductively eliminates the formation of a sextet Mn(II) aziridine complex.
Collapse
|
40
|
Zhao X, Xu W, Chen X, Lin S, Li X, He L, Liao X, Ye G. A comparison of hydrogen abstraction reaction between allyl-type monomers with thioxanthone-based photoinitiators without amine synergists. Front Chem 2022; 10:967836. [PMID: 36118315 PMCID: PMC9478512 DOI: 10.3389/fchem.2022.967836] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2022] [Accepted: 07/25/2022] [Indexed: 11/18/2022] Open
Abstract
The photodriven radical-mediated [3 + 2] cyclization reaction was found to yield polymers efficiently without being hindered by degradative chain transfer. The first reaction is a hydrogen abstraction process in which one hydrogen atom migrates from the α-methylene group of an allyl monomer to the triplet state (or fragments) of the photoinitiator, thus yielding primary allyl radicals as primary radicals and then begins chain propagation via a 3 + 2 cyclization reaction. Allyl ether monomers were found to be significantly higher than other allyl monomers even with the absence of amine-like synergists. In order to clarify the procedure of the hydrogen abstraction mechanism, we used four allyl-type monomers as hydrogen donors and three thioxanthone photoinitiators as hydrogen acceptors by the quantum chemistry method in terms of geometry and energy. The results were interpreted with transition-state theory and the interaction/deformation model. Then, the tunneling factors of hydrogen abstraction reactions were also investigated by Eckart’s correction. The results show allyl ether systems are more reactive than other allyl systems, and it would provide us with new insights into these hydrogen abstractions.
Collapse
Affiliation(s)
- Xiaotian Zhao
- Department of Pharmacy, Chengdu Second Peoples Hospital, Chengdu, China
| | - Wen Xu
- Department of Dermatology, Chengdu Second Peoples Hospital, Chengdu, China
| | - Xi Chen
- Department of Pharmacy, Chengdu Second Peoples Hospital, Chengdu, China
| | - Shibo Lin
- Department of Pharmacy, Chengdu Second Peoples Hospital, Chengdu, China
| | - Xuanhao Li
- Department of Pharmacy, Chengdu Second Peoples Hospital, Chengdu, China
| | - Lihui He
- Department of Pharmacy, Chengdu Second Peoples Hospital, Chengdu, China
| | - Xu Liao
- Department of Pharmacy, Chengdu Second Peoples Hospital, Chengdu, China
| | - Guodong Ye
- Guangzhou Municipal and Guangdong Provincial Key Laboratory of Molecular Target & Clinical Pharmacology, the NMPA and State Key Laboratory of Respiratory Disease, School of Pharmaceutical Sciences and the Fifth Affiliated Hospital, Guangzhou Medical University, Guangzhou, China
- *Correspondence: Guodong Ye,
| |
Collapse
|
41
|
Cammi R, Chen B. Studying and exploring potential energy surfaces of compressed molecules: a fresh theory from the eXtreme Pressure Polarizable Continuum Model. J Chem Phys 2022; 157:114101. [DOI: 10.1063/5.0104269] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
We present a new theory for studying and exploring the potential energy surface of compressed molecular systems as described within the XP-PCM framework. The effective potential energy surface is defined by the sum of the electronic energy of the compressed system and the pressure-volume work that is necessary in order to create the compression cavity at the given condition of pressure. We show that the resulting total energy Gt is related to the electronic energy by a Legendre transform, in which the pressure and volume of the compression cavity are the conjugate variables. We present an analytical expression for the evaluation of the gradient of the total energy ∇Gt to be used for the geometry optimization of equilibrium geometries and transition states of compressed molecular systems. We also show that, as a result of the Legendre transform property, the potential energy surface can be studied explicitly as function of the pressure, leading to an explicit connection with the well-known Hammond postulate. As a proof of concept, we present the application of the theory to studying and determining of the optimized geometry of compressed methane and the transition state of electrocyclic ring-closure of hexatriene and of H-transfer between two methyl radicals.
Collapse
Affiliation(s)
- Roberto Cammi
- Dipartimento di Scienze Chimica della Vita e della Sostenibilità Ambientale, Università degli Studi di Parma, Italy
| | - Bo Chen
- Donostia international physics center, Spain
| |
Collapse
|
42
|
Geiger J, Settels V, Deglmann P, Schäfer A, Bergeler M. Automated input structure generation for single-ended reaction path optimizations. J Comput Chem 2022; 43:1662-1674. [PMID: 35866245 DOI: 10.1002/jcc.26969] [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: 03/15/2022] [Revised: 06/21/2022] [Accepted: 06/22/2022] [Indexed: 11/05/2022]
Abstract
The exploration of a reaction network requires highly automated workflows to avoid error-prone and time-consuming manual steps. In this respect, a major bottleneck is the search for transition-state (TS) structures, which frequently fails and, therefore, makes (manual) revision necessary. In this work, we present a technique for obtaining suitable input structures for automated TS searches based on single-ended reaction path optimization algorithms, which makes subsequent TS searches via this method significantly more robust. First, possible input structures are generated based on the spatial alignment of the reactants. The appropriate orientation of reacting groups is achieved via stepwise rotations along selected torsional degrees of freedom. Second, a ranking of the obtained structures is performed according to selected geometric criteria. The main goals are to properly align the reactive atoms, to avoid hindrance within the reaction channel and to resolve steric clashes between the reactants. The developed procedure has been carefully tested on a variety of examples and provides suitable input structures for TS searches within seconds. The method is in daily use in an industrial setting.
Collapse
Affiliation(s)
- Julian Geiger
- Institute of Chemical Research of Catalonia (ICIQ), The Barcelona Institute of Science and Technology, Tarragona, Spain
| | | | | | | | - Maike Bergeler
- Institute of Chemical Research of Catalonia (ICIQ), The Barcelona Institute of Science and Technology, Tarragona, Spain
| |
Collapse
|
43
|
Liu Y, Qi H, Lei M. Elastic Image Pair Method for Finding Transition States on Potential Energy Surfaces Using Only First Derivatives. J Chem Theory Comput 2022; 18:5108-5115. [PMID: 35771528 DOI: 10.1021/acs.jctc.2c00137] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Herein, an elastic image pair (EIP) method is proposed to search transition states between two potential-energy minima using only first derivatives. In this method, two images are generated, and the spring forces are added to the images to control the distance between the two images. Transition states are reached when the force and the distance of the image pair are both converged. A set of test molecules is optimized using the EIP method, which shows its efficiency in transition state searching compared to other methods. This new method is more stable and reliable in finding transition states with much less computations.
Collapse
Affiliation(s)
- Yangqiu Liu
- State Key Laboratory of Chemical Resource Engineering, Institute of Computational Chemistry, College of Chemistry, Beijing University of Chemical Technology, Beijing 100029, China
| | - Hexiang Qi
- State Key Laboratory of Chemical Resource Engineering, Institute of Computational Chemistry, College of Chemistry, Beijing University of Chemical Technology, Beijing 100029, China
| | - Ming Lei
- State Key Laboratory of Chemical Resource Engineering, Institute of Computational Chemistry, College of Chemistry, Beijing University of Chemical Technology, Beijing 100029, China
| |
Collapse
|
44
|
Schnack-Petersen AK, Koch H, Coriani S, Kjønstad EF. Efficient implementation of molecular CCSD gradients with Cholesky-decomposed electron repulsion integrals. J Chem Phys 2022; 156:244111. [DOI: 10.1063/5.0087261] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023] Open
Abstract
We present an efficient implementation of ground and excited state coupled cluster singles and doubles (CCSD) gradients based on Cholesky-decomposed electron repulsion integrals. Cholesky decomposition and density fitting are both inner projection methods, and, thus, similar implementation schemes can be applied for both methods. One well-known advantage of inner projection methods, which we exploit in our implementation, is that one can avoid storing large V3 O and V4 arrays by instead considering three-index intermediates. Furthermore, our implementation does not require the formation and storage of Cholesky vector derivatives. The new implementation is shown to perform well, with less than 10% of the time spent calculating the gradients in geometry optimizations. Furthermore, the computational time per optimization cycle is significantly lower compared to other implementations based on an inner projection method. We showcase the capabilities of the implementation by optimizing the geometry of the retinal molecule (C20H28O) at the CCSD/aug-cc-pVDZ level of theory.
Collapse
Affiliation(s)
| | - Henrik Koch
- Scuola Normale Superiore, Piazza dei Cavaleri 7, 56126 Pisa, Italy
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Sonia Coriani
- Department of Chemistry, Technical University of Denmark, 2800 Kongens Lyngby, Denmark
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| | - Eirik F. Kjønstad
- Department of Chemistry, Norwegian University of Science and Technology, 7491 Trondheim, Norway
| |
Collapse
|
45
|
Mills AW, Goings JJ, Beck D, Yang C, Li X. Exploring Potential Energy Surfaces Using Reinforcement Machine Learning. J Chem Inf Model 2022; 62:3169-3179. [PMID: 35709515 DOI: 10.1021/acs.jcim.2c00373] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Reinforcement machine learning is implemented to survey a series of model potential energy surfaces and ultimately identify the global minima point. Through sophisticated reward function design, the introduction of an optimizing target, and incorporating physically motivated actions, the reinforcement learning agent is capable of demonstrating advanced decision making. These improved actions allow the agent to successfully converge to an optimal solution more rapidly when compared to an agent trained without the aforementioned modifications. This study showcases the conceptual feasibility of using reinforcement machine learning to solve difficult environments, namely, potential energy surfaces, with multiple, seemingly, correct solutions in the form of local minima regions. Through these results, we hope to encourage extending reinforcement learning to more complicated optimization problems and using these novel techniques to efficiently solve traditionally challenging problems in chemistry.
Collapse
Affiliation(s)
- Alexis W Mills
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Joshua J Goings
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - David Beck
- Department of Chemical Engineering, University of Washington, Seattle, Washington 98195, United States
| | - Chao Yang
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Xiaosong Li
- Department of Chemistry, University of Washington, Seattle, Washington 98195, United States
| |
Collapse
|
46
|
Sumiya Y, Harabuchi Y, Nagata Y, Maeda S. Quantum Chemical Calculations to Trace Back Reaction Paths for the Prediction of Reactants. JACS AU 2022; 2:1181-1188. [PMID: 35647604 PMCID: PMC9131471 DOI: 10.1021/jacsau.2c00157] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/08/2022] [Revised: 03/31/2022] [Accepted: 04/08/2022] [Indexed: 06/12/2023]
Abstract
The long-due development of a computational method for the ab initio prediction of chemical reactants that provide a target compound has been hampered by the combinatorial explosion that occurs when reactions consist of multiple elementary reaction processes. To address this challenge, we have developed a quantum chemical calculation method that can enumerate the reactant candidates from a given target compound by combining an exhaustive automated reaction path search method with a kinetics method for narrowing down the possibilities. Two conventional name reactions were then assessed by tracing back the reaction paths using this new method to determine whether the known reactants could be identified. Our method is expected to be a powerful tool for the prediction of reactants and the discovery of new reactions.
Collapse
Affiliation(s)
- Yosuke Sumiya
- Department
of Chemistry, Faculty of Science, Hokkaido
University, Sapporo, Hokkaido 060-0810, Japan
| | - Yu Harabuchi
- Department
of Chemistry, Faculty of Science, Hokkaido
University, Sapporo, Hokkaido 060-0810, Japan
- Institute
for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido 001-0021, Japan
- ERATO
Maeda Artificial Intelligence for Chemical Reaction Design and Discovery
Project, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
| | - Yuuya Nagata
- Institute
for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido 001-0021, Japan
- ERATO
Maeda Artificial Intelligence for Chemical Reaction Design and Discovery
Project, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
| | - Satoshi Maeda
- Department
of Chemistry, Faculty of Science, Hokkaido
University, Sapporo, Hokkaido 060-0810, Japan
- Institute
for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, Hokkaido 001-0021, Japan
- ERATO
Maeda Artificial Intelligence for Chemical Reaction Design and Discovery
Project, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan
- Research
and Services Division of Materials Data and Integrated System (MaDIS), National Institute for Materials Science (NIMS), Tsukuba, Ibaraki 305-0044, Japan
| |
Collapse
|
47
|
Yang WH, Li YM, Bi JX, Huang R, Shao GF, Fan TE, Liu TD, Wen YH. An Improved Self-Adaptive Differential Evolution with the Neighborhood Search Algorithm for Global Optimization of Bimetallic Clusters. J Chem Inf Model 2022; 62:2398-2408. [PMID: 35533292 DOI: 10.1021/acs.jcim.1c01570] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Global optimization of multicomponent cluster structures is considerably time-consuming due to the existence of a vast number of isomers. In this work, we proposed an improved self-adaptive differential evolution with the neighborhood search (SaNSDE) algorithm and applied it to the global optimization of bimetallic cluster structures. The cross operation was optimized, and an improved basin hopping module was introduced to enhance the searching efficiency of SaNSDE optimization. Taking (PtNi)N (N = 38 or 55) bimetallic clusters as examples, their structures were predicted by using this algorithm. The traditional SaNSDE algorithm was carried out for comparison with the improved SaNSDE algorithm. For all the optimized clusters, the excess energy and the second difference of the energy were calculated to examine their relative stabilities. Meanwhile, the bond order parameters were adopted to quantitatively characterize the cluster structures. The results reveal that the improved SaNSDE algorithm possessed significantly higher searching capability and faster convergence speed than the traditional SaNSDE algorithm. Furthermore, the lowest-energy configurations of (PtNi)38 clusters could be classified as the truncated octahedral and disordered structures. In contrast, all the optimal (PtNi)55 clusters were approximately icosahedral. Our work fully demonstrates the high efficiency of the improved algorithm and advances the development of global optimization algorithms and the structural prediction of multicomponent clusters.
Collapse
Affiliation(s)
- Wei-Hua Yang
- Department of Physics, Xiamen University, Xiamen 361005, China
| | - Ya-Meng Li
- Department of Physics, Xiamen University, Xiamen 361005, China
| | - Jian-Xiang Bi
- Department of Physics, Xiamen University, Xiamen 361005, China
| | - Rao Huang
- Department of Physics, Xiamen University, Xiamen 361005, China
| | - Gui-Fang Shao
- Department of Automation, Xiamen University, Xiamen 361102, China
| | - Tian-E Fan
- College of Automation, Chongqing University of Posts and Telecommunications, Chongqing 400065, China
| | - Tun-Dong Liu
- Department of Automation, Xiamen University, Xiamen 361102, China
| | - Yu-Hua Wen
- Department of Physics, Xiamen University, Xiamen 361005, China
| |
Collapse
|
48
|
Dzib E, Merino G. The hindered rotor theory: A review. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1583] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Affiliation(s)
- Eugenia Dzib
- Departamento de Física Aplicada Centro de Investigación y de Estudios Avanzados Merida Mexico
| | - Gabriel Merino
- Departamento de Física Aplicada Centro de Investigación y de Estudios Avanzados Merida Mexico
| |
Collapse
|
49
|
Tsutsumi T, Ono Y, Taketsugu T. Reaction Space Projector (ReSPer) for Visualizing Dynamic Reaction Routes Based on Reduced-Dimension Space. Top Curr Chem (Cham) 2022; 380:19. [PMID: 35266073 DOI: 10.1007/s41061-022-00377-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/31/2021] [Accepted: 02/21/2022] [Indexed: 11/26/2022]
Abstract
To analyze chemical reaction dynamics based on a reaction path network, we have developed the "Reaction Space Projector" (ReSPer) method with the aid of the dimensionality reduction method. This program has two functions: the construction of a reduced-dimensionality reaction space from a molecular structure dataset, and the projection of dynamic trajectories into the low-dimensional reaction space. In this paper, we apply ReSPer to isomerization and bifurcation reactions of the Au5 cluster and succeed in analyzing dynamic reaction routes involved in multiple elementary reaction processes, constructing complicated networks (called "closed islands") of nuclear permutation-inversion (NPI) isomerization reactions, and elucidating dynamic behaviors in bifurcation reactions with reference to bundles of trajectories. Interestingly, in the second application, we find a correspondence between the contribution ratios in the ability to visualize and the symmetry of the morphology of closed islands. In addition, the third application suggests the existence of boundaries that determine the selectivity in bifurcation reactions, which was discussed in the phase space. The ReSPer program is a versatile and robust tool to clarify dynamic reaction mechanisms based on the reduced-dimensionality reaction space without prior knowledge of target reactions.
Collapse
Affiliation(s)
- Takuro Tsutsumi
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, 060-0810, Japan
| | - Yuriko Ono
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, 001-0021, Japan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of Science, Hokkaido University, Sapporo, 060-0810, Japan.
- Institute for Chemical Reaction Design and Discovery (WPI-ICReDD), Hokkaido University, Sapporo, 001-0021, Japan.
| |
Collapse
|
50
|
Parsaei-Khomami A, Badiei A, Ghavami ZS, Ghasemi JB. A new fluorescence probe for simultaneous determination of Fe2+ and Fe3+ by orthogonal signal correction-principal component regression. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.131978] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
|