1
|
Lyu N, Bergold P, Soley MB, Wang C, Batista VS. Holographic Gaussian Boson Sampling with Matrix Product States on 3D cQED Processors. J Chem Theory Comput 2024; 20:6402-6413. [PMID: 38968605 DOI: 10.1021/acs.jctc.4c00384] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/07/2024]
Abstract
We introduce quantum circuits for simulations of multimode state vectors on 3D circuit quantum electrodynamics (cQED) processors using matrix product state representations. The circuits are demonstrated as applied to simulations of molecular docking based on holographic Gaussian boson sampling (GBS), as illustrated for the binding of a thiol-containing aryl sulfonamide ligand to the tumor necrosis factor-α converting enzyme receptor. We show that cQED devices with a modest number of modes could be employed to simulate multimode systems by repurposing working modes through measurement and reinitialization. We anticipate that a wide range of GBS applications could be implemented on compact 3D cQED processors analogously using the holographic approach. Simulations on qubit-based quantum computers could be implemented analogously using circuits that represent continuous variables in terms of truncated expansions of Fock states.
Collapse
Affiliation(s)
- Ningyi Lyu
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
| | - Paul Bergold
- Department of Mathematics, University of Surrey, Guildford GU2 7XH, U.K
| | - Micheline B Soley
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
- Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, Wisconsin 53706, United States
| | - Chen Wang
- Department of Physics, University of Massachusetts-Amherst, Amherst, Massachusetts 01003, United States
| | - Victor S Batista
- Department of Chemistry, Yale University, New Haven, Connecticut 06520, United States
- Yale Quantum Institute, Yale University, New Haven, Connecticut 06511, United States
| |
Collapse
|
2
|
Kalikadien AV, Mirza A, Hossaini AN, Sreenithya A, Pidko EA. Paving the road towards automated homogeneous catalyst design. Chempluschem 2024; 89:e202300702. [PMID: 38279609 DOI: 10.1002/cplu.202300702] [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: 11/29/2023] [Revised: 12/20/2023] [Indexed: 01/28/2024]
Abstract
In the past decade, computational tools have become integral to catalyst design. They continue to offer significant support to experimental organic synthesis and catalysis researchers aiming for optimal reaction outcomes. More recently, data-driven approaches utilizing machine learning have garnered considerable attention for their expansive capabilities. This Perspective provides an overview of diverse initiatives in the realm of computational catalyst design and introduces our automated tools tailored for high-throughput in silico exploration of the chemical space. While valuable insights are gained through methods for high-throughput in silico exploration and analysis of chemical space, their degree of automation and modularity are key. We argue that the integration of data-driven, automated and modular workflows is key to enhancing homogeneous catalyst design on an unprecedented scale, contributing to the advancement of catalysis research.
Collapse
Affiliation(s)
- Adarsh V Kalikadien
- Inorganic Systems Engineering, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, 2629 HZ, Delft, The Netherlands
| | - Adrian Mirza
- Inorganic Systems Engineering, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, 2629 HZ, Delft, The Netherlands
| | - Aydin Najl Hossaini
- Inorganic Systems Engineering, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, 2629 HZ, Delft, The Netherlands
| | - Avadakkam Sreenithya
- Inorganic Systems Engineering, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, 2629 HZ, Delft, The Netherlands
| | - Evgeny A Pidko
- Inorganic Systems Engineering, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, 2629 HZ, Delft, The Netherlands
| |
Collapse
|
3
|
Stan-Bernhardt A, Glinkina L, Hulm A, Ochsenfeld C. Exploring Chemical Space Using Ab Initio Hyperreactor Dynamics. ACS CENTRAL SCIENCE 2024; 10:302-314. [PMID: 38435517 PMCID: PMC10906254 DOI: 10.1021/acscentsci.3c01403] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Revised: 12/20/2023] [Accepted: 12/21/2023] [Indexed: 03/05/2024]
Abstract
In recent years, first-principles exploration of chemical reaction space has provided valuable insights into intricate reaction networks. Here, we introduce ab initio hyperreactor dynamics, which enables rapid screening of the accessible chemical space from a given set of initial molecular species, predicting new synthetic routes that can potentially guide subsequent experimental studies. For this purpose, different hyperdynamics derived bias potentials are applied along with pressure-inducing spherical confinement of the molecular system in ab initio molecular dynamics simulations to efficiently enhance reactivity under mild conditions. To showcase the advantages and flexibility of the hyperreactor approach, we present a systematic study of the method's parameters on a HCN toy model and apply it to a recently introduced experimental model for the prebiotic formation of glycinal and acetamide in interstellar ices, which yields results in line with experimental findings. In addition, we show how the developed framework enables the study of complicated transitions like the first step of a nonenzymatic DNA nucleoside synthesis in an aqueous environment, where the molecular fragmentation problem of earlier nanoreactor approaches is avoided.
Collapse
Affiliation(s)
- Alexandra Stan-Bernhardt
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5, D-81377 München, Germany
| | - Liubov Glinkina
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5, D-81377 München, Germany
| | - Andreas Hulm
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5, D-81377 München, Germany
| | - Christian Ochsenfeld
- Chair
of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstrasse 5, D-81377 München, Germany
- Max
Planck Institute for Solid State Research, Heisenbergstrasse 1, D-70569 Stuttgart, Germany
| |
Collapse
|
4
|
Aboulfath Y, Bougueroua S, Cimas A, Barth D, Gaigeot MP. Time-Resolved Graphs of Polymorphic Cycles for H-Bonded Network Identification in Flexible Biomolecules. J Chem Theory Comput 2024; 20:1019-1035. [PMID: 38236138 DOI: 10.1021/acs.jctc.3c01031] [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: 01/19/2024]
Abstract
A novel approach based on a coarse-grained representation of topological graphs is proposed for the automatic analysis of molecular dynamics (MD) trajectories of hydrogen-bonded (H-Bonded) flexible biomolecules. Herein, our approach models an H-Bonded biomolecule by its H-Bonded cycles and its graph of cycles in which the vertices and links represent the intersections between these cycles. We propose a methodology in which each identified conformer/isomer from the MD is represented by a well-chosen set of H-Bonded cycles called a minimum cycle basis. The key component is the "polycycles" that distinguish the cycles that play the same polymorphic role in the molecule from the ones that lead to an actual conformational change of the molecule. The relevance of our proposed method is evaluated on MD trajectories of gas-phase biomolecules, for which the covalent bonds are unchanged over time and only the hydrogen bonds change over time. The polygraphs and their time evolution are shown to reveal the dynamicity of the metastructure(s) of the H-Bonded biomolecules while providing polymorphic information on the cycles. Such information on the dynamics and changes in the H-bond network, as some cycles change identity while retaining the same role in the overall structure, is not easily captured at the atomic level of representation. Such information can instead be captured by polymorphic cycles.
Collapse
Affiliation(s)
- Ylène Aboulfath
- Université Paris-Saclay, Univ Versailles Saint Quentin, DAVID, 78035 Versailles, France
| | - Sana Bougueroua
- Université Paris-Saclay, Univ Evry, CY Cergy Paris Université, CNRS, LAMBE, 91025 Evry-Courcouronnes, France
| | - Alvaro Cimas
- Université Paris-Saclay, Univ Evry, CY Cergy Paris Université, CNRS, LAMBE, 91025 Evry-Courcouronnes, France
| | - Dominique Barth
- Université Paris-Saclay, Univ Versailles Saint Quentin, DAVID, 78035 Versailles, France
| | - Marie-Pierre Gaigeot
- Université Paris-Saclay, Univ Evry, CY Cergy Paris Université, CNRS, LAMBE, 91025 Evry-Courcouronnes, France
- Institut Universitaire de France (IUF), 75005 Paris, France
| |
Collapse
|
5
|
Sabadell-Rendón A, Kaźmierczak K, Morandi S, Euzenat F, Curulla-Ferré D, López N. Automated MUltiscale simulation environment. DIGITAL DISCOVERY 2023; 2:1721-1732. [PMID: 38054103 PMCID: PMC10694852 DOI: 10.1039/d3dd00163f] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/23/2023] [Accepted: 11/05/2023] [Indexed: 12/07/2023]
Abstract
Multiscale techniques integrating detailed atomistic information on materials and reactions to predict the performance of heterogeneous catalytic full-scale reactors have been suggested but lack seamless implementation. The largest challenges in the multiscale modeling of reactors can be grouped into two main categories: catalytic complexity and the difference between time and length scales of chemical and transport phenomena. Here we introduce the Automated MUltiscale Simulation Environment AMUSE, a workflow that starts from Density Functional Theory (DFT) data, automates the analysis of the reaction networks through graph theory, prepares it for microkinetic modeling, and subsequently integrates the results into a standard open-source Computational Fluid Dynamics (CFD) code. We demonstrate the capabilities of AMUSE by applying it to the unimolecular iso-propanol dehydrogenation reaction and then, increasing the complexity, to the pre-commercial Pd/In2O3 catalyst employed for the CO2 hydrogenation to methanol. The results show that AMUSE allows the computational investigation of heterogeneous catalytic reactions in a comprehensive way, providing essential information for catalyst design from the atomistic to the reactor scale level.
Collapse
Affiliation(s)
- Albert Sabadell-Rendón
- Institute of Chemical Research of Catalonia (ICIQ-CERCA), The Barcelona Institute of Science and Technology, (BIST) Av. Paisos Catalans 16 Tarragona 43007 Spain
| | - Kamila Kaźmierczak
- TotalEnergies, TotalEnergies One Tech Belgium Zone industrielle C, 7181 Feluy Belgium
| | - Santiago Morandi
- Institute of Chemical Research of Catalonia (ICIQ-CERCA), The Barcelona Institute of Science and Technology, (BIST) Av. Paisos Catalans 16 Tarragona 43007 Spain
- Department of Physical and Inorganic Chemistry, Universitat Rovira i Virgili Campus Sescelades, N4 Block, C. Marcel·lí Domingo 1 Tarragona 43007 Spain
| | - Florian Euzenat
- TotalEnergies Research and Technology Gonfreville, Route Industrielle, Carrefour 4, Port 4864 76700 Rogerville France
| | - Daniel Curulla-Ferré
- TotalEnergies, TotalEnergies One Tech Belgium Zone industrielle C, 7181 Feluy Belgium
| | - Núria López
- Institute of Chemical Research of Catalonia (ICIQ-CERCA), The Barcelona Institute of Science and Technology, (BIST) Av. Paisos Catalans 16 Tarragona 43007 Spain
| |
Collapse
|
6
|
Zhang Y, Xu C, Lan Z. Automated Exploration of Reaction Networks and Mechanisms Based on Metadynamics Nanoreactor Simulations. J Chem Theory Comput 2023. [PMID: 38031422 DOI: 10.1021/acs.jctc.3c00752] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2023]
Abstract
We developed an automated approach to construct a complex reaction network and explore the reaction mechanisms for numerous reactant molecules by integrating several theoretical approaches. Nanoreactor-type molecular dynamics was used to generate possible chemical reactions, in which the metadynamics was used to overcome the reaction barriers, and the semiempirical GFN2-xTB method was used to reduce the computational cost. Reaction events were identified from trajectories using the hidden Markov model based on the evolution of the molecular connectivity. This provided the starting points for further transition-state searches at the electronic structure levels of density functional theory to obtain the reaction mechanism. Finally, the entire reaction network containing multiple pathways was built. The feasibility and efficiency of the automated construction of the reaction network were investigated using the HCHO and NH3 biomolecular reaction and the reaction network for a multispecies system comprising dozens of HCN and H2O molecules. The results indicated that the proposed approach provides a valuable and effective tool for the automated exploration of the reaction networks.
Collapse
Affiliation(s)
- Yutai Zhang
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Chao Xu
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| | - Zhenggang Lan
- Guangdong Provincial Key Laboratory of Chemical Pollution and Environmental Safety and MOE Key Laboratory of Environmental Theoretical Chemistry, SCNU Environmental Research Institute, School of Environment, South China Normal University, Guangzhou 510006, P. R. China
| |
Collapse
|
7
|
Hashemi A, Bougueroua S, Gaigeot MP, Pidko EA. HiREX: High-Throughput Reactivity Exploration for Extended Databases of Transition-Metal Catalysts. J Chem Inf Model 2023; 63:6081-6094. [PMID: 37738303 PMCID: PMC10565810 DOI: 10.1021/acs.jcim.3c00660] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Indexed: 09/24/2023]
Abstract
A method is introduced for the automated analysis of reactivity exploration for extended in silico databases of transition-metal catalysts. The proposed workflow is designed to tackle two key challenges for bias-free mechanistic explorations on large databases of catalysts: (1) automated exploration of the chemical space around each catalyst with unique structural and chemical features and (2) automated analysis of the resulting large chemical data sets. To address these challenges, we have extended the application of our previously developed ReNeGate method for bias-free reactivity exploration and implemented an automated analysis procedure to identify the classes of reactivity patterns within specific catalyst groups. Our procedure applied to an extended series of representative Mn(I) pincer complexes revealed correlations between structural and reactive features, pointing to new channels for catalyst transformation under the reaction conditions. Such an automated high-throughput virtual screening of systematically generated hypothetical catalyst data sets opens new opportunities for the design of high-performance catalysts as well as an accelerated method for expert bias-free high-throughput in silico reactivity exploration.
Collapse
Affiliation(s)
- Ali Hashemi
- Inorganic
Systems Engineering, Department of Chemical Engineering, Faculty of
Applied Sciences, Delft University of Technology, Van der Maasweg 9, Delft 2629 HZ, The Netherlands
| | - Sana Bougueroua
- Laboratoire
Analyse et Modélisation pour la Biologie et l’Environnement
(LAMBE) UMR8587, Paris-Saclay, Univ Evry,
CY Cergy Paris Université, CNRS, LAMBE UMR8587, Evry-Courcouronnes 91025, France
| | - Marie-Pierre Gaigeot
- Laboratoire
Analyse et Modélisation pour la Biologie et l’Environnement
(LAMBE) UMR8587, Paris-Saclay, Univ Evry,
CY Cergy Paris Université, CNRS, LAMBE UMR8587, Evry-Courcouronnes 91025, France
| | - Evgeny A. Pidko
- Inorganic
Systems Engineering, Department of Chemical Engineering, Faculty of
Applied Sciences, Delft University of Technology, Van der Maasweg 9, Delft 2629 HZ, The Netherlands
| |
Collapse
|
8
|
Petrus E, Garay-Ruiz D, Reiher M, Bo C. Multi-Time-Scale Simulation of Complex Reactive Mixtures: How Do Polyoxometalates Form? J Am Chem Soc 2023; 145:18920-18930. [PMID: 37496164 DOI: 10.1021/jacs.3c05514] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Understanding the dynamics of reactive mixtures still challenges both experiments and theory. A relevant example can be found in the chemistry of molecular metal-oxide nanoclusters, also known as polyoxometalates. The high number of species potentially involved, the interconnectivity of the reaction network, and the precise control of the pH and concentrations needed in the synthesis of such species make the theoretical/computational treatment of such processes cumbersome. This work addresses this issue relying on a unique combination of recently developed computational methods that tackle the construction, kinetic simulation, and analysis of complex chemical reaction networks. By using the Bell-Evans-Polanyi approximation for estimating activation energies, and an accurate and robust linear scaling for correcting the computed pKa values, we report herein multi-time-scale kinetic simulations for the self-assembly processes of polyoxotungstates that comprise 22 orders of magnitude, from tens of femtoseconds to months of reaction time. This very large time span was required to reproduce very fast processes such as the acid/base equilibria (at 10-12 s), relatively slow reactions such as the formation of key clusters such as the metatungstate (at 103 s), and the very slow assembly of the decatungstate (at 106 s). Analysis of the kinetic data and of the reaction network topology shed light onto the details of the main reaction mechanisms, which explains the origin of kinetic and thermodynamic control followed by the reaction. Simulations at alkaline pH fully reproduce experimental evidence since clusters do not form under those conditions.
Collapse
Affiliation(s)
- Enric Petrus
- Institute of Chemical Research of Catalonia (ICIQ), The Barcelona Institute of Science and Technology (BIST), Avenida Països Catalans, 16, Tarragona 43007, Spain
| | - Diego Garay-Ruiz
- Institute of Chemical Research of Catalonia (ICIQ), The Barcelona Institute of Science and Technology (BIST), Avenida Països Catalans, 16, Tarragona 43007, Spain
| | - Markus Reiher
- Department of Chemistry and Applied Biosciences, ETH Zürich, Vladimir-Prelog-Weg 2, Zürich 8093, Switzerland
| | - Carles Bo
- Institute of Chemical Research of Catalonia (ICIQ), The Barcelona Institute of Science and Technology (BIST), Avenida Països Catalans, 16, Tarragona 43007, Spain
- Departament de Química Física i Inorgànica, Universitat Rovira i Virgili, Marcel•li Domingo s/n, Tarragona 43007, Spain
| |
Collapse
|
9
|
Yang W, Filonenko GA, Pidko EA. Performance of homogeneous catalysts viewed in dynamics. Chem Commun (Camb) 2023; 59:1757-1768. [PMID: 36683401 PMCID: PMC9910057 DOI: 10.1039/d2cc05625a] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
Effective assessment of catalytic performance is the foundation for the rational design and development of new catalysts with superior performance. The ubiquitous screening/optimization studies use reaction yields as the sole performance metric in an approach that often neglects the complexity of the catalytic system and intrinsic reactivities of the catalysts. Using an example of hydrogenation catalysis, we examine the transient behavior of catalysts that are often encountered in activation, deactivation and catalytic turnover processes. Each of these processes and the reaction environment in which they take place are gradually shown to determine the real-time catalyst speciation and the resulting kinetics of the overall catalytic reaction. As a result, the catalyst performance becomes a complex and time-dependent metric defined by multiple descriptors apart from the reaction yield. This behaviour is not limited to hydrogenation catalysis and affects various catalytic transformations. In this feature article, we discuss these catalytically relevant descriptors in an attempt to arrive at a comprehensive depiction of catalytic performance.
Collapse
Affiliation(s)
- Wenjun Yang
- Inorganic Systems Engineering group, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 9, 2629 HZ, Delft, The Netherlands.
| | - Georgy A. Filonenko
- Inorganic Systems Engineering group, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 92629 HZDelftThe Netherlands
| | - Evgeny A. Pidko
- Inorganic Systems Engineering group, Department of Chemical Engineering, Faculty of Applied Sciences, Delft University of Technology, Van der Maasweg 92629 HZDelftThe Netherlands
| |
Collapse
|
10
|
Bougueroua S, Aboulfath Y, Barth D, Gaigeot MP. Algorithmic graph theory for post-processing molecular dynamics trajectories. Mol Phys 2023. [DOI: 10.1080/00268976.2022.2162456] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/21/2023]
Affiliation(s)
- Sana Bougueroua
- Université Paris-Saclay, Univ Evry, CNRS, LAMBE UMR8587, Evry-Courcouronnes, France
| | - Ylène Aboulfath
- Université Paris-Saclay, Univ Versailles SQ, DAVID, Versailles, France
| | - Dominique Barth
- Université Paris-Saclay, Univ Versailles SQ, DAVID, Versailles, France
| | - Marie-Pierre Gaigeot
- Université Paris-Saclay, Univ Evry, CNRS, LAMBE UMR8587, Evry-Courcouronnes, France
| |
Collapse
|