1
|
Stulajter MM, Rappoport D. Reaction Networks Resemble Low-Dimensional Regular Lattices. J Chem Theory Comput 2024. [PMID: 39236261 DOI: 10.1021/acs.jctc.4c00810] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/07/2024]
Abstract
The computational exploration, manipulation, and design of complex chemical reactions face fundamental challenges related to the high-dimensional nature of potential energy surfaces (PESs) that govern reactivity. Accurately modeling complex reactions is crucial for understanding the chemical processes involved in, for example, organocatalysis, autocatalytic cycles, and one-pot molecular assembly. Our prior research demonstrated that discretizing PESs using heuristics based on bond breaking and bond formation produces a reaction network representation with a low-dimensional structure (metric space). We now find that these stoichiometry-preserving reaction networks possess additional, though approximate, structure and resemble low-dimensional regular lattices with a small amount of random edge rewiring. The heuristics-based discretization thus generates a nonlinear dimensionality reduction by a factor of 10 with an a posteriori error measure (probability of random rewiring). The structure becomes evident through a comparative analysis of CHNO reaction networks of varying stoichiometries against a panel of size-matched generative network models, taking into account their local, metric, and global properties. The generative models include random networks (Erdős-Rényi and bipartite random networks), regular lattices (periodic and nonperiodic), and network models with a tunable level of "randomness" (Watts-Strogatz graphs and regular lattices with random rewiring). The CHNO networks are simultaneously closely matched in all these properties by 3-4-dimensional regular lattices with 10% or less of edges randomly rewired. The effective dimensionality reduction is found to be independent of the system size, stoichiometry, and ruleset, suggesting that search and sampling algorithms for PESs of complex chemical reactions can be effectively leveraged.
Collapse
Affiliation(s)
- Miko M Stulajter
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
- Computational Science Research Center, San Diego State University, San Diego, California 92182, United States
| | - Dmitrij Rappoport
- Department of Chemistry, University of California Irvine, Irvine, California 92697, United States
| |
Collapse
|
2
|
Fakhoury Z, Sosso GC, Habershon S. Contact-Map-Driven Exploration of Heterogeneous Protein-Folding Paths. J Chem Theory Comput 2024. [PMID: 39228261 DOI: 10.1021/acs.jctc.4c00878] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/05/2024]
Abstract
We have recently shown how physically realizable protein-folding pathways can be generated using directed walks in the space of inter-residue contact-maps; combined with a back-transformation to move from protein contact-maps to Cartesian coordinates, we have demonstrated how this approach can generate protein-folding trajectory ensembles without recourse to molecular dynamics. In this article, we demonstrate that this framework can be used to study a challenging protein-folding problem that is known to exhibit two different folding paths which were previously identified through molecular dynamics simulation at several different temperatures. From the viewpoint of protein-folding mechanism prediction, this particular problem is extremely challenging to address, specifically involving folding to an identical nontrivial compact native structure along distinct pathways defined by heterogeneous secondary structural elements. Here, we show how our previously reported contact-map-based protein-folding strategy can be significantly enhanced to enable accurate and robust prediction of heterogeneous folding paths by (i) introducing a novel topologically informed metric for comparing two protein contact maps, (ii) reformulating our graph-represented folding path generation, and (iii) introducing a new and more reliable structural back-mapping algorithm. These changes improve the reliability of generating structurally sound folding intermediates and dramatically decrease the number of physically irrelevant folding intermediates generated by our previous simulation strategy. Most importantly, we demonstrate how our enhanced folding algorithm can successfully identify the alternative folding mechanisms of a multifolding-pathway protein, in line with direct molecular dynamics simulations.
Collapse
Affiliation(s)
- Ziad Fakhoury
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, U.K
| | - Gabriele C Sosso
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, U.K
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, U.K
| |
Collapse
|
3
|
Gilkes J, Storr MT, Maurer RJ, Habershon S. Predicting Long-Time-Scale Kinetics under Variable Experimental Conditions with Kinetica.jl. J Chem Theory Comput 2024; 20:5196-5214. [PMID: 38829777 PMCID: PMC11209948 DOI: 10.1021/acs.jctc.4c00333] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2024] [Revised: 05/13/2024] [Accepted: 05/13/2024] [Indexed: 06/05/2024]
Abstract
Predicting the degradation processes of molecules over long time scales is a key aspect of industrial materials design. However, it is made computationally challenging by the need to construct large networks of chemical reactions that are relevant to the experimental conditions that kinetic models must mirror, with every reaction requiring accurate kinetic data. Here, we showcase Kinetica.jl, a new software package for constructing large-scale chemical reaction networks in a fully automated fashion by exploring chemical reaction space with a kinetics-driven algorithm; coupled to efficient machine-learning models of activation energies for sampled elementary reactions, we show how this approach readily enables generation and kinetic characterization of networks containing ∼103 chemical species and ≃104-105 reactions. Symbolic-numeric modeling of the generated reaction networks is used to allow for flexible, efficient computation of kinetic profiles under experimentally realizable conditions such as continuously variable temperature regimes, enabling direct connection between bottom-up reaction networks and experimental observations. Highly efficient propagation of long-time-scale kinetic profiles is required for automated reaction network refinement and is enabled here by a new discrete kinetic approximation. The resulting Kinetica.jl simulation package therefore enables automated generation, characterization, and long-time-scale modeling of complex chemical reaction systems. We demonstrate this for hydrocarbon pyrolysis simulated over time scales of seconds, using transient temperature profiles representing those of tubular flow reactor experiments.
Collapse
Affiliation(s)
- Joe Gilkes
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, CV4 7AL Coventry, U.K.
- EPSRC
HetSys Centre for Doctoral Training, University
of Warwick, Gibbet Hill
Rd, CV4 7AL Coventry, U.K.
| | | | - Reinhard J. Maurer
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, CV4 7AL Coventry, U.K.
- Department
of Physics, University of Warwick, Gibbet Hill Road, CV4 7AL Coventry, U.K.
| | - Scott Habershon
- Department
of Chemistry, University of Warwick, Gibbet Hill Road, CV4 7AL Coventry, U.K.
| |
Collapse
|
4
|
Ong A, Wong ZC, Chin KLO, Loh WW, Chua MH, Ang SJ, Lim JYC. Enhancing the photocatalytic upcycling of polystyrene to benzoic acid: a combined computational-experimental approach for acridinium catalyst design. Chem Sci 2024; 15:1061-1067. [PMID: 38239702 PMCID: PMC10793207 DOI: 10.1039/d3sc06388g] [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: 11/29/2023] [Accepted: 12/07/2023] [Indexed: 01/22/2024] Open
Abstract
Converting polystyrene into value-added oxygenated aromatic compounds is an attractive end-of-life upcycling strategy. However, identification of appropriate catalysts often involves laborious and time-consuming empirical screening. Herein, after demonstrating the feasibility of using acridinium salts for upcycling polystyrene into benzoic acid by photoredox catalysis for the first time, we applied low-cost descriptor-based combinatorial in silico screening to predict the photocatalytic performance of a family of potential candidates. Through this approach, we identified a non-intuitive fluorinated acridinium catalyst that outperforms other candidates for converting polystyrene to benzoic acid in useful yields at low catalyst loadings (≤5 mol%). In addition, this catalyst also proved effective with real-life polystyrene waste containing dyes and additives. Our study underscores the potential of computer-aided catalyst design for valorizing polymeric waste into essential chemical feedstock for a more sustainable future.
Collapse
Affiliation(s)
- Albert Ong
- Institute of Materials Research and Engineering (IMRE), Agency for Science, Technology and Research (A*STAR) 2 Fusionopolis Way, Innovis #08-03 Singapore 138634 Republic of Singapore
| | - Zi Cheng Wong
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, Connexis, #16-16 Singapore 138632 Republic of Singapore
| | - Kang Le Osmund Chin
- Institute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR) 1 Pesek Road, Jurong Island Singapore 627833 Republic of Singapore
| | - Wei Wei Loh
- Institute of Materials Research and Engineering (IMRE), Agency for Science, Technology and Research (A*STAR) 2 Fusionopolis Way, Innovis #08-03 Singapore 138634 Republic of Singapore
| | - Ming Hui Chua
- Institute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR) 1 Pesek Road, Jurong Island Singapore 627833 Republic of Singapore
| | - Shi Jun Ang
- Institute of High Performance Computing (IHPC), Agency for Science, Technology and Research (A*STAR) 1 Fusionopolis Way, Connexis, #16-16 Singapore 138632 Republic of Singapore
- Institute of Sustainability for Chemicals, Energy and Environment (ISCE2), Agency for Science, Technology and Research (A*STAR) 1 Pesek Road, Jurong Island Singapore 627833 Republic of Singapore
| | - Jason Y C Lim
- Institute of Materials Research and Engineering (IMRE), Agency for Science, Technology and Research (A*STAR) 2 Fusionopolis Way, Innovis #08-03 Singapore 138634 Republic of Singapore
- Department of Materials Science and Engineering, National University of Singapore (NUS) 9 Engineering Drive 1 Singapore 117576 Republic of Singapore
| |
Collapse
|
5
|
Krep L, Schmalz F, Solbach F, Komissarov L, Nevolianis T, Kopp WA, Verstraelen T, Leonhard K. A Reactive Molecular Dynamics Study of Chlorinated Organic Compounds. Part II: A ChemTraYzer Study of Chlorinated Dibenzofuran Formation and Decomposition Processes. Chemphyschem 2022; 24:e202200783. [PMID: 36511423 DOI: 10.1002/cphc.202200783] [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/18/2022] [Revised: 12/09/2022] [Accepted: 12/12/2022] [Indexed: 12/14/2022]
Abstract
In our two-paper series, we first present the development of ReaxFF CHOCl parameters using the recently published ParAMS parametrization tool. In this second part, we update the reactive Molecular Dynamics - Quantum Mechanics coupling scheme ChemTraYzer and combine it with our new ReaxFF parameters from Part I to study formation and decomposition processes of chlorinated dibenzofurans. We introduce a self-learning method for recovering failed transition-state searches that improves the overall ChemTraYzer transition-state search success rate by 10 percentage points to a total of 48 %. With ChemTraYzer, we automatically find and quantify more than 500 reactions using transition state theory and DFT. Among the discovered chlorinated dibenzofuran reactions are numerous reactions that are new to the literature. In three case studies, we discuss the set of reactions that are most relevant to the dibenzofuran literature: (i) bimolecular reactions of the chlorinated-dibenzofuran precursors phenoxy radical and 1,3,5-trichlorobenzene, (ii) dibenzofuran chlorination and pyrolysis, and (iii) oxidation of chlorinated dibenzofurans.
Collapse
Affiliation(s)
- Lukas Krep
- Institute of Technical Thermodynamics, RWTH Aachen University, North Rhine-Westphalia, 52062, Aachen, Germany
| | - Felix Schmalz
- Institute of Technical Thermodynamics, RWTH Aachen University, North Rhine-Westphalia, 52062, Aachen, Germany
| | - Florian Solbach
- Institute of Technical Thermodynamics, RWTH Aachen University, North Rhine-Westphalia, 52062, Aachen, Germany
| | - Leonid Komissarov
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark-Zwijnaarde 46, B-9052, Ghent, Belgium
| | - Thomas Nevolianis
- Institute of Technical Thermodynamics, RWTH Aachen University, North Rhine-Westphalia, 52062, Aachen, Germany
| | - Wassja A Kopp
- Institute of Technical Thermodynamics, RWTH Aachen University, North Rhine-Westphalia, 52062, Aachen, Germany
| | - Toon Verstraelen
- Center for Molecular Modeling (CMM), Ghent University, Technologiepark-Zwijnaarde 46, B-9052, Ghent, Belgium
| | - Kai Leonhard
- Institute of Technical Thermodynamics, RWTH Aachen University, North Rhine-Westphalia, 52062, Aachen, Germany
| |
Collapse
|
6
|
Ismail I, Chantreau Majerus R, Habershon S. Graph-Driven Reaction Discovery: Progress, Challenges, and Future Opportunities. J Phys Chem A 2022; 126:7051-7069. [PMID: 36190262 PMCID: PMC9574932 DOI: 10.1021/acs.jpca.2c06408] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/07/2022] [Revised: 09/22/2022] [Indexed: 11/29/2022]
Abstract
Graph-based descriptors, such as bond-order matrices and adjacency matrices, offer a simple and compact way of categorizing molecular structures; furthermore, such descriptors can be readily used to catalog chemical reactions (i.e., bond-making and -breaking). As such, a number of graph-based methodologies have been developed with the goal of automating the process of generating chemical reaction network models describing the possible mechanistic chemistry in a given set of reactant species. Here, we outline the evolution of these graph-based reaction discovery schemes, with particular emphasis on more recent methods incorporating graph-based methods with semiempirical and ab initio electronic structure calculations, minimum-energy path refinements, and transition state searches. Using representative examples from homogeneous catalysis and interstellar chemistry, we highlight how these schemes increasingly act as "virtual reaction vessels" for interrogating mechanistic questions. Finally, we highlight where challenges remain, including issues of chemical accuracy and calculation speeds, as well as the inherent challenge of dealing with the vast size of accessible chemical reaction space.
Collapse
Affiliation(s)
- Idil Ismail
- Department of Chemistry, University
of Warwick, CoventryCV4 7AL, United Kingdom
| | | | - Scott Habershon
- Department of Chemistry, University
of Warwick, CoventryCV4 7AL, United Kingdom
| |
Collapse
|
7
|
Ismail I, Robertson C, Habershon S. Successes and challenges in using machine-learned activation energies in kinetic simulations. J Chem Phys 2022; 157:014109. [DOI: 10.1063/5.0096027] [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/14/2022] Open
Abstract
The prediction of the thermodynamic and kinetic properties of chemical reactions is increasingly being addressed by machine-learning (ML) methods such as artificial neural networks (ANNs). While a number of recent studies have reported success in predicting chemical reaction activation energies, less attention has focused on how the accuracy of ML predictions filter through to predictions of macroscopic observables. Here, we consider the impact of the uncertainty associated with ML prediction of activation energies on observable properties of chemical reaction networks, as given by microkinetics simulations based on ML-predicted reaction rates. After training an ANN to predict activation energies given standard molecular descriptors for reactants and products alone, we performed microkinetics simulations of three different prototypical reaction networks: formamide decomposition, aldol reactions and decomposition of 3-hydroperoxypropanal. We find that the kinetic modelling predictions can be in excellent agreement with corresponding simulations performed with ab initio calculations, but this is dependent on the inherent energetic landscape of the networks. We use these simulations to suggest some guidelines for when ML-based activation energies can be reliable, and when one should take more care in applications to kinetics modelling.
Collapse
Affiliation(s)
| | | | - Scott Habershon
- Department of Chemistry, University of Warwick, United Kingdom
| |
Collapse
|
8
|
Rasmussen MH, Jensen JH. Fast and automated identification of reactions with low barriers using meta-MD simulations. PEERJ PHYSICAL CHEMISTRY 2022. [DOI: 10.7717/peerj-pchem.22] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
We test our meta-molecular dynamics (MD) based approach for finding low-barrier (<30 kcal/mol) reactions on uni- and bimolecular reactions extracted from the barrier dataset developed by Grambow, Pattanaik & Green (2020). For unimolecular reactions the meta-MD simulations identify 25 of the 26 products found by Grambow, Pattanaik & Green (2020), while the subsequent semiempirical screening eliminates an additional four reactions due to an overestimation of the reaction energies or estimated barrier heights relative to DFT. In addition, our approach identifies 36 reactions not found by Grambow, Pattanaik & Green (2020), 10 of which are <30 kcal/mol. For bimolecular reactions the meta-MD simulations identify 19 of the 20 reactions found by Grambow, Pattanaik & Green (2020), while the subsequent semiempirical screening eliminates an additional reaction. In addition, we find 34 new low-barrier reactions. For bimolecular reactions we found that it is necessary to “encourage” the reactants to go to previously undiscovered products, by including products found by other MD simulations when computing the biasing potential as well as decreasing the size of the molecular cavity in which the MD occurs, until a reaction is observed. We also show that our methodology can find the correct products for two reactions that are more representative of those encountered in synthetic organic chemistry. The meta-MD hyperparameters used in this study thus appear to be generally applicable to finding low-barrier reactions.
Collapse
|
9
|
Krep L, Roy IS, Kopp W, Schmalz F, Huang C, Leonhard K. Efficient Reaction Space Exploration with ChemTraYzer-TAD. J Chem Inf Model 2022; 62:890-902. [PMID: 35142513 DOI: 10.1021/acs.jcim.1c01197] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The development of a reaction model is often a time-consuming process, especially if unknown reactions have to be found and quantified. To alleviate the reaction modeling process, automated procedures for reaction space exploration are highly desired. We present ChemTraYzer-TAD, a new reactive molecular dynamics acceleration technique aimed at efficient reaction space exploration. The new method is based on the basin confinement strategy known from the temperature-accelerated dynamics (TAD) acceleration method. Our method features integrated ChemTraYzer bond-order processing steps for the automatic and on-the-fly determination of the positions of virtual walls in configuration space that confine the system in a potential energy basin. We use the example of 1,3-dioxolane-4-hydroperoxide-2-yl radical oxidation to show that ChemTraYzer-TAD finds more than 100 different parallel reactions for the given set of reactants in less than 2 ns of simulation time. Among the many observed reactions, ChemTraYzer-TAD finds the expected typical low-temperature reactions despite the use of extremely high simulation temperatures up to 5000 K. Our method also finds a new concerted β-scission plus O2 addition with a lower reaction barrier than the literature-known and so-far dominant β-scission.
Collapse
Affiliation(s)
- Lukas Krep
- Institute of Technical Thermodynamics, RWTH Aachen University, Aachen 52062, Germany
| | - Indu Sekhar Roy
- Institute of Technical Thermodynamics, RWTH Aachen University, Aachen 52062, Germany
| | - Wassja Kopp
- Institute of Technical Thermodynamics, RWTH Aachen University, Aachen 52062, Germany
| | - Felix Schmalz
- Institute of Technical Thermodynamics, RWTH Aachen University, Aachen 52062, Germany
| | - Can Huang
- Institute of Technical Thermodynamics, RWTH Aachen University, Aachen 52062, Germany
| | - Kai Leonhard
- Institute of Technical Thermodynamics, RWTH Aachen University, Aachen 52062, Germany
| |
Collapse
|
10
|
Garay-Ruiz D, Álvarez-Moreno M, Bo C, Martínez-Núñez E. New Tools for Taming Complex Reaction Networks: The Unimolecular Decomposition of Indole Revisited. ACS PHYSICAL CHEMISTRY AU 2022; 2:225-236. [PMID: 36855573 PMCID: PMC9718323 DOI: 10.1021/acsphyschemau.1c00051] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The level of detail attained in the computational description of reaction mechanisms can be vastly improved through tools for automated chemical space exploration, particularly for systems of small to medium size. Under this approach, the unimolecular decomposition landscape for indole was explored through the automated reaction mechanism discovery program AutoMeKin. Nevertheless, the sheer complexity of the obtained mechanisms might be a hindrance regarding their chemical interpretation. In this spirit, the new Python library amk-tools has been designed to read and manipulate complex reaction networks, greatly simplifying their overall analysis. The package provides interactive dashboards featuring visualizations of the network, the three-dimensional (3D) molecular structures and vibrational normal modes of all chemical species, and the corresponding energy profiles for selected pathways. The combination of the joined mechanism generation and postprocessing workflow with the rich chemistry of indole decomposition enabled us to find new details of the reaction (obtained at the CCSD(T)/aug-cc-pVTZ//M06-2X/MG3S level of theory) that were not reported before: (i) 16 pathways leading to the formation of HCN and NH3 (via amino radical); (ii) a barrierless reaction between methylene radical and phenyl isocyanide, which might be an operative mechanism under the conditions of the interstellar medium; and (iii) reaction channels leading to both hydrogen cyanide and hydrogen isocyanide, of potential astrochemical interest as the computed HNC/HCN ratios greatly exceed the calculated equilibrium value at very low temperatures. The reported reaction networks can be very valuable to supplement databases of kinetic data, which is of remarkable interest for pyrolysis and astrochemical studies.
Collapse
Affiliation(s)
- Diego Garay-Ruiz
- Institute
of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Avinguda Països Catalans,
16, 43007 Tarragona, Spain,Departament
de Química Física i Inorgànica, Universitat Rovira i Virgili (URV), Marcel·lí Domingo s/n, 43007 Tarragona, Spain
| | - Moises Álvarez-Moreno
- Institute
of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Avinguda Països Catalans,
16, 43007 Tarragona, Spain
| | - Carles Bo
- Institute
of Chemical Research of Catalonia (ICIQ), Barcelona Institute of Science & Technology (BIST), Avinguda Països Catalans,
16, 43007 Tarragona, Spain,Departament
de Química Física i Inorgànica, Universitat Rovira i Virgili (URV), Marcel·lí Domingo s/n, 43007 Tarragona, Spain,
| | - Emilio Martínez-Núñez
- Departmento
de Química Física, Facultade de Química, Universidade de Santiago de Compostela, 15782 Santiago
de Compostela, Spain,
| |
Collapse
|
11
|
Chantreau Majerus R, Robertson C, Habershon S. Assessing and rationalizing the performance of Hessian update schemes for reaction path Hamiltonian rate calculations. J Chem Phys 2021; 155:204112. [PMID: 34852478 DOI: 10.1063/5.0064685] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The reaction path Hamiltonian (RPH) can be used to calculate chemical reaction rate constants, going beyond transition-state theory in taking account of recrossing by providing an approximation to the dynamic transmission coefficient. However, the RPH necessitates the calculation of the Hessian matrix at a number of points along the minimum energy path; the associated computational cost stands as a bottleneck in RPH calculations, especially if one is interested in using high-accuracy electronic structure methods. In this work, four different Hessian update schemes (symmetric rank-1, Powell-symmetric Broyden, Bofill, and TS-BFGS updates) are assessed to see whether or not they reliably reproduce calculated transmission coefficients for three different chemical reactions. Based on the reactions investigated, the symmetric rank-1 Hessian update was the least appropriate for RPH construction, giving different transmission coefficients from the standard analytical Hessian approach, as well as inconsistent frequencies and coupling properties. The Bofill scheme, the Powell-symmetric Broyden scheme, and the TS-BFGS scheme were the most reliable Hessian update methods, with transmission coefficients that were in good agreement with those calculated by the standard RPH calculations. The relative accuracy of the different Hessian update schemes is further rationalized by investigating the approximated Coriolis and curvature coupling terms along the reaction-path, providing insight into when these schemes would be expected to work well. Furthermore, the associated computational cost associated with the RPH calculations was substantially reduced by the tested update schemes. Together, these results provide useful rules-of-thumb for using Hessian update schemes in RPH simulations.
Collapse
Affiliation(s)
- R Chantreau Majerus
- Molecular Analytical Science Centre for Doctoral Training, Senate House, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - C Robertson
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - S Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
12
|
Shannon RJ, Deeks HM, Burfoot E, Clark E, Jones AJ, Mulholland AJ, Glowacki DR. Exploring human-guided strategies for reaction network exploration: Interactive molecular dynamics in virtual reality as a tool for citizen scientists. J Chem Phys 2021; 155:154106. [PMID: 34686059 DOI: 10.1063/5.0062517] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The emerging fields of citizen science and gamification reformulate scientific problems as games or puzzles to be solved. Through engaging the wider non-scientific community, significant breakthroughs may be made by analyzing citizen-gathered data. In parallel, recent advances in virtual reality (VR) technology are increasingly being used within a scientific context and the burgeoning field of interactive molecular dynamics in VR (iMD-VR) allows users to interact with dynamical chemistry simulations in real time. Here, we demonstrate the utility of iMD-VR as a medium for gamification of chemistry research tasks. An iMD-VR "game" was designed to encourage users to explore the reactivity of a particular chemical system, and a cohort of 18 participants was recruited to playtest this game as part of a user study. The reaction game encouraged users to experiment with making chemical reactions between a propyne molecule and an OH radical, and "molecular snapshots" from each game session were then compiled and used to map out reaction pathways. The reaction network generated by users was compared to existing literature networks demonstrating that users in VR capture almost all the important reaction pathways. Further comparisons between humans and an algorithmic method for guiding molecular dynamics show that through using citizen science to explore these kinds of chemical problems, new approaches and strategies start to emerge.
Collapse
Affiliation(s)
- Robin J Shannon
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Helen M Deeks
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Eleanor Burfoot
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Edward Clark
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | - Alex J Jones
- School of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
| | | | - David R Glowacki
- ArtSci Foundation International, 5th floor Mariner House, Bristol, BS1 4QD, United Kingdom
| |
Collapse
|
13
|
Martínez-Núñez E, Barnes GL, Glowacki DR, Kopec S, Peláez D, Rodríguez A, Rodríguez-Fernández R, Shannon RJ, Stewart JJP, Tahoces PG, Vazquez SA. AutoMeKin2021: An open-source program for automated reaction discovery. J Comput Chem 2021; 42:2036-2048. [PMID: 34387374 DOI: 10.1002/jcc.26734] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2021] [Revised: 07/16/2021] [Accepted: 07/27/2021] [Indexed: 01/10/2023]
Abstract
AutoMeKin2021 is an updated version of tsscds2018, a program for the automated discovery of reaction mechanisms (J. Comput. Chem. 2018, 39, 1922). This release features a number of new capabilities: rare-event molecular dynamics simulations to enhance reaction discovery, extension of the original search algorithm to study van der Waals complexes, use of chemical knowledge, a new search algorithm based on bond-order time series analysis, statistics of the chemical reaction networks, a web application to submit jobs, and other features. The source code, manual, installation instructions and the website link are available at: https://rxnkin.usc.es/index.php/AutoMeKin.
Collapse
Affiliation(s)
- Emilio Martínez-Núñez
- Department of Physical Chemistry, University of Santiago de Compostela, Santiago de Compostela, Spain
| | - George L Barnes
- Department of Chemistry and Biochemistry, Siena College, Loudonville, New York, USA
| | - David R Glowacki
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| | - Sabine Kopec
- Institut de Sciences Moléculaires d'Orsay, UMR 8214, Université Paris-Sud - Université Paris-Saclay, Orsay, France
| | - Daniel Peláez
- Institut de Sciences Moléculaires d'Orsay, UMR 8214, Université Paris-Sud - Université Paris-Saclay, Orsay, France
| | - Aurelio Rodríguez
- Galicia Supercomputing Center (CESGA), Santiago de Compostela, Spain
| | | | - Robin J Shannon
- Centre for Computational Chemistry, School of Chemistry, University of Bristol, Bristol, UK
| | | | - Pablo G Tahoces
- Department of Electronics and Computer Science, University of Santiago de Compostela, Santiago de Compostela, Spain
| | - Saulo A Vazquez
- Department of Physical Chemistry, University of Santiago de Compostela, Santiago de Compostela, Spain
| |
Collapse
|
14
|
Abstract
Computational methods have emerged as a powerful tool to augment traditional experimental molecular catalyst design by providing useful predictions of catalyst performance and decreasing the time needed for catalyst screening. In this perspective, we discuss three approaches for computational molecular catalyst design: (i) the reaction mechanism-based approach that calculates all relevant elementary steps, finds the rate and selectivity determining steps, and ultimately makes predictions on catalyst performance based on kinetic analysis, (ii) the descriptor-based approach where physical/chemical considerations are used to find molecular properties as predictors of catalyst performance, and (iii) the data-driven approach where statistical analysis as well as machine learning (ML) methods are used to obtain relationships between available data/features and catalyst performance. Following an introduction to these approaches, we cover their strengths and weaknesses and highlight some recent key applications. Furthermore, we present an outlook on how the currently applied approaches may evolve in the near future by addressing how recent developments in building automated computational workflows and implementing advanced ML models hold promise for reducing human workload, eliminating human bias, and speeding up computational catalyst design at the same time. Finally, we provide our viewpoint on how some of the challenges associated with the up-and-coming approaches driven by automation and ML may be resolved.
Collapse
Affiliation(s)
- Ademola Soyemi
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, AL 35487, USA.
| | - Tibor Szilvási
- Department of Chemical and Biological Engineering, The University of Alabama, Tuscaloosa, AL 35487, USA.
| |
Collapse
|
15
|
Shannon RJ, Martínez-Núñez E, Shalashilin DV, Glowacki DR. ChemDyME: Kinetically Steered, Automated Mechanism Generation through Combined Molecular Dynamics and Master Equation Calculations. J Chem Theory Comput 2021; 17:4901-4912. [PMID: 34283599 DOI: 10.1021/acs.jctc.1c00335] [Citation(s) in RCA: 18] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
In many scientific fields, there is an interest in understanding the way in which chemical networks evolve. The chemical networks which researchers focus upon have become increasingly complex, and this has motivated the development of automated methods for exploring chemical reactivity or conformational change in a "black-box" manner, harnessing modern computing resources to automate mechanism discovery. In this work, we present a new approach to automated mechanism generation which couples molecular dynamics and statistical rate theory to automatically find kinetically important reactions and then solve the time evolution of the species in the evolving network. The key to this chemical network mapping through combined dynamics and ME simulation approach is the concept of "kinetic convergence", whereby the search for new reactions is constrained to those species which are kinetically favorable at the conditions of interest. We demonstrate the capability of the new approach for two systems, a well-studied combustion system and a multiple oxygen addition system relevant to atmospheric aerosol formation.
Collapse
Affiliation(s)
- Robin J Shannon
- School of Chemistry, University of Bristol, Bristol BS8 1TS, U.K
| | - Emilio Martínez-Núñez
- Department of Physical Chemistry, University of Santiago de Compostela, Santiago de Compostela 15705, Spain
| | | | - David R Glowacki
- ArtSci International Foundation, 5th floor Mariner House, Bristol BS1 4QD, U.K
| |
Collapse
|
16
|
Robertson C, Hyland R, Lacey AJD, Havens S, Habershon S. Identifying Barrierless Mechanisms for Benzene Formation in the Interstellar Medium Using Permutationally Invariant Reaction Discovery. J Chem Theory Comput 2021; 17:2307-2322. [DOI: 10.1021/acs.jctc.1c00046] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
| | - Ross Hyland
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Andrew J. D. Lacey
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Sebastian Havens
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Scott Habershon
- Department of Chemistry, University of Warwick, Coventry CV4 7AL, United Kingdom
| |
Collapse
|
17
|
Blau SM, Patel HD, Spotte-Smith EWC, Xie X, Dwaraknath S, Persson KA. A chemically consistent graph architecture for massive reaction networks applied to solid-electrolyte interphase formation. Chem Sci 2021; 12:4931-4939. [PMID: 34163740 PMCID: PMC8179555 DOI: 10.1039/d0sc05647b] [Citation(s) in RCA: 19] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/13/2020] [Accepted: 02/23/2021] [Indexed: 01/09/2023] Open
Abstract
Modeling reactivity with chemical reaction networks could yield fundamental mechanistic understanding that would expedite the development of processes and technologies for energy storage, medicine, catalysis, and more. Thus far, reaction networks have been limited in size by chemically inconsistent graph representations of multi-reactant reactions (e.g. A + B → C) that cannot enforce stoichiometric constraints, precluding the use of optimized shortest-path algorithms. Here, we report a chemically consistent graph architecture that overcomes these limitations using a novel multi-reactant representation and iterative cost-solving procedure. Our approach enables the identification of all low-cost pathways to desired products in massive reaction networks containing reactions of any stoichiometry, allowing for the investigation of vastly more complex systems than previously possible. Leveraging our architecture, we construct the first ever electrochemical reaction network from first-principles thermodynamic calculations to describe the formation of the Li-ion solid electrolyte interphase (SEI), which is critical for passivation of the negative electrode. Using this network comprised of nearly 6000 species and 4.5 million reactions, we interrogate the formation of a key SEI component, lithium ethylene dicarbonate. We automatically identify previously proposed mechanisms as well as multiple novel pathways containing counter-intuitive reactions that have not, to our knowledge, been reported in the literature. We envision that our framework and data-driven methodology will facilitate efforts to engineer the composition-related properties of the SEI - or of any complex chemical process - through selective control of reactivity.
Collapse
Affiliation(s)
- Samuel M Blau
- Energy Technologies Area, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| | - Hetal D Patel
- Department of Materials Science and Engineering, University of California Berkeley CA 94720 USA
- Materials Science Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| | - Evan Walter Clark Spotte-Smith
- Department of Materials Science and Engineering, University of California Berkeley CA 94720 USA
- Materials Science Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| | - Xiaowei Xie
- Materials Science Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
- College of Chemistry, University of California Berkeley CA 94720 USA
| | - Shyam Dwaraknath
- Materials Science Division, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| | - Kristin A Persson
- Department of Materials Science and Engineering, University of California Berkeley CA 94720 USA
- Molecular Foundry, Lawrence Berkeley National Laboratory Berkeley CA 94720 USA
| |
Collapse
|
18
|
Robertson C, Habershon S. Simple position and orientation preconditioning scheme for minimum energy path calculations. J Comput Chem 2021; 42:761-770. [PMID: 33617652 DOI: 10.1002/jcc.26495] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Revised: 01/16/2021] [Accepted: 01/22/2021] [Indexed: 11/08/2022]
Abstract
Minimum-energy path (MEP) calculations, such as those typified by the nudged elastic band method, require input of reactant and product molecular configurations at initialization. In the case of reactions involving more than one molecule, generating initial reactant and product configurations requires careful consideration of the relative position and orientations of the reactive molecules in order to ensure that the resulting MEP calculation proceeds without converging on an alternative reaction-path, and without requiring excessive numbers of optimization iterations; as such, this initial system set-up is most commonly performed "by hand," with an expert user arranging reactive molecules in space to ensure that the following MEP calculation runs smoothly. In this Article, we introduce a simple preconditioning scheme which replaces this labor-intensive, human-knowledge-based step with an automated deterministic computational scheme. In our approach, initial reactant and product configurations are generated such that steric hindrance between reactive molecules is minimized in the reactant and product configurations, while also simultaneously requiring minimal structural differences between the reactants and products. The method is demonstrated using a benchmark test-set of >3400 organic molecular reactions, where comparison of the reactant/product configurations generated using our approach compare very well to initial configurations which were generated on an ad hoc basis.
Collapse
Affiliation(s)
- Christopher Robertson
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry, UK
| | - Scott Habershon
- Department of Chemistry and Centre for Scientific Computing, University of Warwick, Coventry, UK
| |
Collapse
|
19
|
Bryden MA, Zysman-Colman E. Organic thermally activated delayed fluorescence (TADF) compounds used in photocatalysis. Chem Soc Rev 2021; 50:7587-7680. [PMID: 34002736 DOI: 10.1039/d1cs00198a] [Citation(s) in RCA: 135] [Impact Index Per Article: 45.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
Organic compounds that show Thermally Activated Delayed Fluorescence (TADF) have become wildly popular as next-generation emitters in organic light emitting diodes (OLEDs). Since 2016, a subset of these have found increasing use as photocatalysts. This review comprehensively highlights their potential by documenting the diversity of the reactions where an organic TADF photocatalyst can be used in lieu of a noble metal complex photocatalyst. Beyond the small number of TADF photocatalysts that have been used to date, the analysis conducted within this review reveals the wider potential of organic donor-acceptor TADF compounds as photocatalysts. A discussion of the benefits of compounds showing TADF for photocatalysis is presented, which paints a picture of a very promising future for organic photocatalyst development.
Collapse
Affiliation(s)
- Megan Amy Bryden
- Organic Semiconductor Centre, EaStCHEM School of Chemistry, University of St Andrews, St Andrews, KY16 9ST, UK.
| | - Eli Zysman-Colman
- Organic Semiconductor Centre, EaStCHEM School of Chemistry, University of St Andrews, St Andrews, KY16 9ST, UK.
| |
Collapse
|
20
|
DFT-Based Microkinetic Simulations: A Bridge Between Experiment and Theory in Synthetic Chemistry. TOP ORGANOMETAL CHEM 2020. [DOI: 10.1007/3418_2020_44] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|