1
|
Jia L, Brémond É, Zaida L, Gaüzère B, Tognetti V, Joubert L. Predicting redox potentials by graph-based machine learning methods. J Comput Chem 2024; 45:2383-2396. [PMID: 38923574 DOI: 10.1002/jcc.27380] [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: 11/08/2023] [Revised: 03/25/2024] [Accepted: 04/19/2024] [Indexed: 06/28/2024]
Abstract
The evaluation of oxidation and reduction potentials is a pivotal task in various chemical fields. However, their accurate prediction by theoretical computations, which is a complementary task and sometimes the only alternative to experimental measurement, may be often resource-intensive and time-consuming. This paper addresses this challenge through the application of machine learning techniques, with a particular focus on graph-based methods (such as graph edit distances, graph kernels, and graph neural networks) that are reviewed to enlighten their deep links with theoretical chemistry. To this aim, we establish the ORedOx159 database, a comprehensive, homogeneous (with reference values stemming from density functional theory calculations), and reliable resource containing 318 one-electron reduction and oxidation reactions and featuring 159 large organic compounds. Subsequently, we provide an instructive overview of the good practice in machine learning and of commonly utilized machine learning models. We then assess their predictive performances on the ORedOx159 dataset through extensive analyses. Our simulations using descriptors that are computed in an almost instantaneous way result in a notable improvement in prediction accuracy, with mean absolute error (MAE) values equal to 5.6 kcal mol- 1 for reduction and 7.2 kcal mol- 1 for oxidation potentials, which paves a way toward efficient in silico design of new electrochemical systems.
Collapse
Affiliation(s)
- Linlin Jia
- The PRG Group, Institute of Computer Science, University of Bern, Bern, Switzerland
| | - Éric Brémond
- Université Paris Cité, ITODYS, CNRS, Paris, France
| | | | - Benoit Gaüzère
- LITIS, Univ Rouen Normandie, INSA Rouen Normandie, Université Le Havre Normandie, Normandie Univ, Rouen, France
| | - Vincent Tognetti
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, Mont St Aignan Cedex, France
| | - Laurent Joubert
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, Mont St Aignan Cedex, France
| |
Collapse
|
2
|
Pakzad F, Eskandari K. Exploring the influence of metal cations on individual hydrogen bonds in Watson-Crick guanine-cytosine DNA base pair: An interacting quantum atoms analysis. J Comput Chem 2024; 45:2397-2408. [PMID: 38922952 DOI: 10.1002/jcc.27441] [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: 01/30/2024] [Revised: 04/25/2024] [Accepted: 05/16/2024] [Indexed: 06/28/2024]
Abstract
This study delves into the nature of individual hydrogen bonds and the relationship between metal cations and hydrogen bonding in the Watson-Crick guanine-cytosine (GC) base pair and its alkali and alkaline earth cation-containing complexes (Mn+-GC). The findings reveal how metal cations affect the nature and strength of individual hydrogen bonds. The study employs interacting quantum atoms (IQA) analysis to comprehensively understand three individual hydrogen bonds within the GC base pair and its cationic derivatives. These analyses unveil the nature and strength of hydrogen bonds and serve as a valuable reference for exploring the impact of cations (and other factors) on each hydrogen bond. All the H ⋯ D interactions (H is hydrogen and D is oxygen or nitrogen) in the GC base pair are primarily electrostatic in nature, with the charge transfer component playing a substantial role. Introducing a metal cation perturbs all H ⋯ D interatomic interactions in the system, weakening the nearest hydrogen bond to the cation (indicated by a) and reinforcing the other (b and c) interactions. Notably, the interaction a, the strongest H ⋯ D interaction in the GC base pair, becomes the weakest in the Mn+-GC complexes. A broader perspective on the stability of GC and Mn+-GC complexes is provided through interacting quantum fragments (IQF) analysis. This approach considers all pairwise interactions between fragments and intra-fragment components, offering a complete view of the factors that stabilize and destabilize GC and Mn+-GC complexes. The IQF analysis underscores the importance of electron sharing, with the dominant contribution arising from the inter-fragment exchange-correlation term, in shaping and sustaining GC and Mn+-GC complexes. From this point of view, alkaline and alkaline earth cations have distinct effects, with alkaline cations generally weakening inter-fragment interactions and alkaline earth cations strengthening them. In addition, IQA and IQF calculations demonstrate that the hydration of cations led to small changes in the hydrogen bonding network. Finally, the IQA interatomic energies associated with the hydrogen bonds and also inter-fragment interaction energies provide robust indicators for characterizing hydrogen bonds and complex stability, showing a strong correlation with total interaction energies.
Collapse
Affiliation(s)
- F Pakzad
- Department of Chemistry, Isfahan University of Technology, Isfahan, Iran
| | - K Eskandari
- Department of Chemistry, Isfahan University of Technology, Isfahan, Iran
| |
Collapse
|
3
|
Tognetti V, Joubert L. Exchange-correlation effects in interatomic energies for pure density functionals and their application to the molecular energy prediction. J Comput Chem 2024; 45:2270-2283. [PMID: 38847367 DOI: 10.1002/jcc.27431] [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/28/2023] [Revised: 04/02/2024] [Accepted: 04/28/2024] [Indexed: 08/15/2024]
Abstract
In this proof-of-concept paper, we show how exchange-correlation effects can be simply recovered for interatomic energies within the interacting quantum atoms decomposition when local, gradient generalized, or meta-gradient generalized approximations are used in density functional theory (DFT) calculations. We also demonstrate how inhomogeneity and non-local effects can be introduced even from a pure local scheme, without resorting to any orbital information. Finally, we provide numerical evidence on a database of selected energetic molecules that this decomposition scheme can be efficiently used to build accurate models for the prediction of molecular energies from an initial "cheap" DFT calculation.
Collapse
Affiliation(s)
- Vincent Tognetti
- University of Rouen Normandie, INSA Rouen Normandie, CNRS, Normandie University, Rouen, France
| | - Laurent Joubert
- University of Rouen Normandie, INSA Rouen Normandie, CNRS, Normandie University, Rouen, France
| |
Collapse
|
4
|
Marie N, Ma JA, Tognetti V, Cahard D. Photocatalyzed Cascade Hydrogen Atom Transfers for Assembly of Multi-Substituted α-SCF 3 and α-SCF 2H Cyclopentanones. Angew Chem Int Ed Engl 2024; 63:e202407689. [PMID: 38845586 DOI: 10.1002/anie.202407689] [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/23/2024] [Indexed: 07/23/2024]
Abstract
A photocatalyzed formal (3+2) cycloaddition has been developed to construct original polysubstituted α-SCF3 cyclopentanones in a regio- and diastereoselective manner. This building block approach leverages trifluoromethylthio alkynes and branched/linear aldehydes, as readily available reaction partners, in consecutive hydrogen atom transfers and C-C bond formations. Difluoromethylthio alkynes are also compatible substrates. Furthermore, the potential for telescoped reaction starting from alcohols instead of aldehydes was demonstrated, as well as process automatization and scale-up under continuous microflow conditions. This prompted density functional theory (DFT) calculations to support a radical-mediated cascade process.
Collapse
Affiliation(s)
- Nicolas Marie
- CNRS, UMR 6014 COBRA, Univ Rouen Normandie, INSA Rouen Normandie, Normandie Univ, INC3M FR 3038, F-76000, Rouen, France
| | - Jun-An Ma
- Department of Chemistry, Tianjin Key Laboratory of Molecular Optoelectronic Sciences, and Frontiers Science Center for Synthetic Biology (Ministry of Education), Tianjin University, Tianjin, 300072, China
| | - Vincent Tognetti
- CNRS, UMR 6014 COBRA, Univ Rouen Normandie, INSA Rouen Normandie, Normandie Univ, INC3M FR 3038, F-76000, Rouen, France
| | - Dominique Cahard
- CNRS, UMR 6014 COBRA, Univ Rouen Normandie, INSA Rouen Normandie, Normandie Univ, INC3M FR 3038, F-76000, Rouen, France
| |
Collapse
|
5
|
Cador A, Morell C, Tognetti V, Joubert L, Popelier PLA. Determining the Factors Accounting for Reaction Selectivity: A Relative Energy Gradient - Interacting Quantum Atoms and Natural Bonding Orbitals Study. Chemphyschem 2024; 25:e202400163. [PMID: 38747261 DOI: 10.1002/cphc.202400163] [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: 02/13/2024] [Revised: 05/14/2024] [Indexed: 07/14/2024]
Abstract
Identifying the main physicochemical properties accounting for the course of a reaction is of utmost importance to rationalize chemical syntheses. To this aim, the relative energy gradient (REG) method is an appealing approach because it is an unbiased and automatic process to extract the most relevant pieces of energy information. Initially formulated within the interacting quantum atoms (IQA) framework for a single reaction, here we extend the REG method to natural bond orbitals (NBO) analysis and to the case of two competitive processes. This development enables the determination of the driving forces of any chemical selectivity. We illustrate the extended REG method on the case study of ring opening in cyclobutenes, which is an important instance of the so-called torquoselectivity.
Collapse
Affiliation(s)
- Aël Cador
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, 76821, Mont St Aignan Cedex, France
| | - Christophe Morell
- University of Lyon, Université Lyon 1 et CNRS UMR 5280, Institut des Sciences Analytiques, 5 rue de la Doua, 69100, Villeurbanne, France
| | - Vincent Tognetti
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, 76821, Mont St Aignan Cedex, France
| | - Laurent Joubert
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, 76821, Mont St Aignan Cedex, France
| | - Paul L A Popelier
- Department of Chemistry, University of Manchester, Manchester, M13 9PL, Great Britain
| |
Collapse
|
6
|
Ahirwar MB, Gadre SR, Deshmukh MM. Molecular Tailoring Approach for the Direct Estimation of Individual Noncovalent Interaction Energies in Molecular Systems. J Phys Chem A 2024; 128:6099-6115. [PMID: 39037864 DOI: 10.1021/acs.jpca.4c01176] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/24/2024]
Abstract
The noncovalent interactions (NCIs) are omnipresent in chemistry, physics, and biology. The study of such interactions offers insights into various physicochemical phenomena. Some indirect approaches proposed in the literature for exploring the NCIs are briefly reviewed in Section 1 of this Perspective. These include: (i) Shift in the stretching frequency of an X-Y bond involved in X-Y···Z interaction. (ii) Topological analysis of molecular electron density. (iii) Empirical equations derived employing experimental and theoretical quantities. However, a direct method for estimating individual intramolecular/intermolecular interaction energies has been conspicuous by its absence from the literature. We have developed a molecular tailoring approach (MTA)-based method enabling a direct and reliable estimation of the energy of intra- as well as intermolecular interactions. This method offers a direct and reliable estimation of these interactions, in particular of the hydrogen bonds (HB) in molecules/weakly bound clusters along with the respective cooperativity contribution. In Section 2, the basis of our method is discussed, along with some illustrative examples. The application of this method to a variety of molecules and clusters, with a special emphasis on estimating the HB energy along with the energy of other NCIs is presented in Section 3. Section 4 discusses some computational strategies for applying our method to large molecular clusters. The last Section provides a summary and a discussion on future developments.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling, & Simulation, Savitribai Phule Pune University, Pune 411007, India
- Department of Chemistry, Savitribai Phule Pune University, Pune 411007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| |
Collapse
|
7
|
Cui LJ, Liu YQ, Pan S, Cui ZH, Frenking G. Unusual quadruple bonds featuring collective interaction-type σ bonds between first octal-row atoms in the alkaline-earth compounds AeOLi 2 (Ae = Be-Ba). Chem Sci 2024:d4sc01979b. [PMID: 39176245 PMCID: PMC11337084 DOI: 10.1039/d4sc01979b] [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/25/2024] [Accepted: 07/31/2024] [Indexed: 08/24/2024] Open
Abstract
Quantum chemical calculations are reported for the complexes of alkaline earth metals AeOLi2 (Ae = Be-Ba) at the BP86-D3(BJ)/def2-QZVPP and CCSD(T)/def2-QZVPPQZVPP levels. The nature of the Ae-OLi2 bond has been analyzed with a variety of methods. The AeOLi2 molecules exhibit an unprecedented σ donor bond Ae→OLi2 where the (n)s2 lone-pair electrons of the Ae atom are donated to vacant O-Li2 antibonding orbitals having the largest coefficient at lithium. This is a covalent bond where the accumulation of the associated electronic charge is located at two positions above and below the Ae-OLi2 axis. The bifurcated component of orbital interactions is structurally related to the recently proposed collective bonding model, but exhibits a completely different type of bonding. The most stable isomer of AeOLi2 has a C 2v geometry and a singlet (1A1) electronic ground state. The bond dissociation energy (BDE) of the Ae-OLi2 bonds exhibits a zig-zag trend from BeOLi2 to BaOLi2, with BeOLi2 having the largest BDE (D e = 73.0 kcal mol-1) and MgOLi2 possessing the lowest BDE (D e = 42.3 kcal mol-1) at the CCSD(T) level. The calculation of the atomic partial charges by the Hirshfeld and Voronoi methods suggests that Be and Mg carry small negative charges in the lighter molecules whereas the heavier atoms Ca-Ba have small positive charges. In contrast, the NBO and QTAIM methods give positive charges for all Ae atoms that are larger for Ca-Ba than that calculated by the Hirshfeld and Voronoi approaches. The molecules AeOLi2 have large dipole moments where the negative end is at the Ae atom with the polarity Ae→OLi2. The largest dipole moments are predicted for the lighter species BeOLi2 and MgOLi2 and the smallest value is calculated for BaOLi2. The calculation of the vibrational spectra shows a significant red-shift toward lower wave numbers for the Ae-OLi2 stretching mode in comparison to diatomic AeO. Besides the Ae→OLi2 σ-donor bonds there are also three dative bonds due to Ae←OLi2 backdonation which consist of one σ bond and two π bonds. The appearance of strong Ae→OLi2 σ donation leads to quadruple bonds AeOLi2 in all systems AeOLi2, even for the lightest species with Ae = Be, Mg. The valence orbitals of Ca, Sr, and Ba, which are involved in the dative interactions, are the (n)s and (n-1)d AOs whereas Be and Mg use their (n)s and (n)p AOs. The EDA-NOCV results are supported by the AdNDP calculations which give four 2c-2e bonding orbitals. Three bonding orbitals have occupation numbers ∼2. One σ orbital has smaller occupation numbers between 1.32 and 1.73 due to the delocalization to the lithium atoms. The analysis of the electronic structure with the ELF method suggests multicenter bonds with mainly trisynaptic and tetrasynaptic basins, which also support the results of the EDA-NOCV calculations.
Collapse
Affiliation(s)
- Li-Juan Cui
- Institute of Atomic and Molecular Physics, Jilin University Changchun 130023 China
| | - Yu-Qian Liu
- Institute of Atomic and Molecular Physics, Jilin University Changchun 130023 China
| | - Sudip Pan
- Institute of Atomic and Molecular Physics, Jilin University Changchun 130023 China
| | - Zhong-Hua Cui
- Institute of Atomic and Molecular Physics, Jilin University Changchun 130023 China
- Key Laboratory of Physics and Technology for Advanced Batteries (Ministry of Education), Jilin University Changchun 130023 China
| | - Gernot Frenking
- Philipps-Universität Marburg Hans-Meerwein-Strasse 4 D-35043 Marburg Germany
- Institute of Advanced Synthesis, School of Chemistry and Molecular Engineering, Nanjing Tech University Nanjing 211816 China
| |
Collapse
|
8
|
Genoni A, Martín Pendás Á. Critical assessment of the x-ray restrained wave function approach: Advantages, drawbacks, and perspectives for density functional theory and periodic ab initio calculations. J Chem Phys 2024; 160:234108. [PMID: 38899684 DOI: 10.1063/5.0213247] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2024] [Accepted: 06/03/2024] [Indexed: 06/21/2024] Open
Abstract
The x-ray restrained wave function (XRW) method is a quantum crystallographic technique to extract wave functions compatible with experimental x-ray diffraction data. The approach looks for wave functions that minimize the energies of the investigated systems and also reproduce sets of x-ray structure factors. Given the strict relationship between x-ray structure factors and electron distributions, the strategy practically allows determining wave functions that correspond to given (usually experimental) electron densities. In this work, the capabilities of the XRW approach were further tested. The aim was to evaluate whether the XRW technique could serve as a tool for suggesting new exchange-correlation functionals for density functional theory or refining existing ones. Additionally, the ability of the method to address the influences of the crystalline environment was also assessed. The outcomes of XRW computations were thus compared to those of traditional gas-phase, embedding quantum mechanics/molecular mechanics, and fully periodic calculations. The results revealed that, irrespective of the initial conditions, the XRW computations practically yield a consensus electron density, in contrast to the currently employed density functional approximations (DFAs), which tend to give a too large range of electron distributions. This is encouraging in view of exploiting the XRW technique to develop improved functionals. Conversely, the calculations also emphasized that the XRW method is limited in its ability to effectively address the influences of the crystalline environment. This underscores the need for a periodic XRW technique, which would allow further untangling the shortcomings of DFAs from those inherent to the XRW approach.
Collapse
Affiliation(s)
- Alessandro Genoni
- Université de Lorraine & CNRS, Laboratoire de Physique et Chimie Théoriques (LPCT), UMR CNRS 7019, 1 Boulevard Arago, 57078 Metz, France
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, Avenida Julian Clavería 8, 33006 Oviedo, Spain
| |
Collapse
|
9
|
Mikkelsen JES, Jensen F. Ambiguities in Decomposing Molecular Polarizability into Atomic Charge Flow and Induced Dipole Contributions. J Phys Chem A 2024; 128:4168-4175. [PMID: 38743593 DOI: 10.1021/acs.jpca.4c01890] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/16/2024]
Abstract
The molecular dipole polarizability can be decomposed into components corresponding to the charge flow between atoms and changes in atomic dipole moments. Such decompositions are recognized to depend on how atoms are defined within a molecule, as, for example, by Hirshfeld, iterative Stockholder, or quantum topology partitioning of the electron density. For some of these, however, there are significant differences between the numerical results obtained by analytical response methods and finite field calculations. We show that this difference is due to analytical response methods accounting for (only) the change in electron density by a perturbation, while finite field methods may also include a component corresponding to a perturbation-dependent change in the definition of an atom within a molecule. For some atom-in-molecule definitions, such as the iterative Hirshfeld, iterative Stockholder, and quantum topology methods, the latter effect significantly increases the charge flow component. The decomposition of molecular polarizability into atomic charge flow and induced dipole components thus depends on whether the atom-in-molecule definition is taken to be perturbation-dependent.
Collapse
Affiliation(s)
- Jonas E S Mikkelsen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus, Denmark
| | - Frank Jensen
- Department of Chemistry, Aarhus University, Langelandsgade 140, DK-8000 Aarhus, Denmark
| |
Collapse
|
10
|
Gallegos M, Vassilev-Galindo V, Poltavsky I, Martín Pendás Á, Tkatchenko A. Explainable chemical artificial intelligence from accurate machine learning of real-space chemical descriptors. Nat Commun 2024; 15:4345. [PMID: 38773090 DOI: 10.1038/s41467-024-48567-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/11/2023] [Accepted: 04/24/2024] [Indexed: 05/23/2024] Open
Abstract
Machine-learned computational chemistry has led to a paradoxical situation in which molecular properties can be accurately predicted, but they are difficult to interpret. Explainable AI (XAI) tools can be used to analyze complex models, but they are highly dependent on the AI technique and the origin of the reference data. Alternatively, interpretable real-space tools can be employed directly, but they are often expensive to compute. To address this dilemma between explainability and accuracy, we developed SchNet4AIM, a SchNet-based architecture capable of dealing with local one-body (atomic) and two-body (interatomic) descriptors. The performance of SchNet4AIM is tested by predicting a wide collection of real-space quantities ranging from atomic charges and delocalization indices to pairwise interaction energies. The accuracy and speed of SchNet4AIM breaks the bottleneck that has prevented the use of real-space chemical descriptors in complex systems. We show that the group delocalization indices, arising from our physically rigorous atomistic predictions, provide reliable indicators of supramolecular binding events, thus contributing to the development of Explainable Chemical Artificial Intelligence (XCAI) models.
Collapse
Affiliation(s)
- Miguel Gallegos
- Department of Analytical and Physical Chemistry, University of Oviedo, E-33006, Oviedo, Spain
| | | | - Igor Poltavsky
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg
| | - Ángel Martín Pendás
- Department of Analytical and Physical Chemistry, University of Oviedo, E-33006, Oviedo, Spain.
| | - Alexandre Tkatchenko
- Department of Physics and Materials Science, University of Luxembourg, L-1511, Luxembourg City, Luxembourg.
| |
Collapse
|
11
|
Pérez-Bitrián A, Munárriz J, Krause KB, Schlögl J, Hoffmann KF, Sturm JS, Hadi AN, Teutloff C, Wiesner A, Limberg C, Riedel S. Questing for homoleptic mononuclear manganese complexes with monodentate O-donor ligands. Chem Sci 2024; 15:5564-5572. [PMID: 38638238 PMCID: PMC11023055 DOI: 10.1039/d4sc00543k] [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: 01/23/2024] [Accepted: 02/26/2024] [Indexed: 04/20/2024] Open
Abstract
Compounds containing Mn-O bonds are of utmost importance in biological systems and catalytic processes. Nevertheless, mononuclear manganese complexes containing all O-donor ligands are still rare. Taking advantage of the low tendency of the pentafluoroorthotellurate ligand (teflate, OTeF5) to bridge metal centers, we have synthesized two homoleptic manganese complexes with monomeric structures and an all O-donor coordination sphere. The tetrahedrally distorted MnII anion, [Mn(OTeF5)4]2-, can be described as a high spin d5 complex (S = 5/2), as found experimentally (magnetic susceptibility measurements and EPR spectroscopy) and using theoretical calculations (DFT and CASSCF/NEVPT2). The high spin d4 electronic configuration (S = 2) of the MnIII anion, [Mn(OTeF5)5]2-, was also determined experimentally and theoretically, and a square pyramidal geometry was found to be the most stable one for this complex. Finally, the bonding situation in both complexes was investigated by means of the Interacting Quantum Atoms (IQA) methodology and compared to that of hypothetical mononuclear fluoromanganates. Within each pair of [MnXn]2- (n = 4, 5) species (X = OTeF5, F), the Mn-X interaction is found to be comparable, therefore proving that the similar electronic properties of the teflate and the fluoride are also responsible for the stabilization of these unique species.
Collapse
Affiliation(s)
- Alberto Pérez-Bitrián
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin Fabeckstraße 34/36 Berlin 14195 Germany
- Institut für Chemie, Humboldt-Universität zu Berlin Brook-Taylor-Straße 2 Berlin 12489 Germany
| | - Julen Munárriz
- Departamento de Química Física and Instituto de Biocomputación y Física de Sistemas Complejos (BIFI), Universidad de Zaragoza Pedro Cerbuna 12 Zaragoza 50009 Spain
| | - Konstantin B Krause
- Institut für Chemie, Humboldt-Universität zu Berlin Brook-Taylor-Straße 2 Berlin 12489 Germany
| | - Johanna Schlögl
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin Fabeckstraße 34/36 Berlin 14195 Germany
| | - Kurt F Hoffmann
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin Fabeckstraße 34/36 Berlin 14195 Germany
| | - Johanna S Sturm
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin Fabeckstraße 34/36 Berlin 14195 Germany
| | - Amiera N Hadi
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin Fabeckstraße 34/36 Berlin 14195 Germany
| | - Christian Teutloff
- Fachbereich Physik, Freie Universität Berlin Arnimallee 14 Berlin 14195 Germany
| | - Anja Wiesner
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin Fabeckstraße 34/36 Berlin 14195 Germany
| | - Christian Limberg
- Institut für Chemie, Humboldt-Universität zu Berlin Brook-Taylor-Straße 2 Berlin 12489 Germany
| | - Sebastian Riedel
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin Fabeckstraße 34/36 Berlin 14195 Germany
| |
Collapse
|
12
|
Falcioni F, Bennett S, Stroer-Jarvis P, Popelier PLA. Probing Non-Covalent Interactions through Molecular Balances: A REG-IQA Study. Molecules 2024; 29:1043. [PMID: 38474554 DOI: 10.3390/molecules29051043] [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: 02/01/2024] [Revised: 02/19/2024] [Accepted: 02/21/2024] [Indexed: 03/14/2024] Open
Abstract
The interaction energies of two series of molecular balances (1-X with X = H, Me, OMe, NMe2 and 2-Y with Y = H, CN, NO2, OMe, NMe2) designed to probe carbonyl…carbonyl interactions were analysed at the B3LYP/6-311++G(d,p)-D3 level of theory using the energy partitioning method of Interacting Quantum Atoms/Fragments (IQA/IQF). The partitioned energies are analysed by the Relative Energy Gradient (REG) method, which calculates the correlation between these energies and the total energy of a system, thereby explaining the role atoms have in the energetic behaviour of the total system. The traditional "back-of-the-envelope" open and closed conformations of molecular balances do not correspond to those of the lowest energy. Hence, more care needs to be taken when considering which geometries to use for comparison with the experiment. The REG-IQA method shows that the 1-H and 1-OMe balances behave differently to the 1-Me and 1-NMe2 balances because the latter show more prominent electrostatics between carbonyl groups and undergoes a larger dihedral rotation due to the bulkiness of the functional groups. For the 2-Y balance, REG-IQA shows the same behaviour across the series as the 1-H and 1-OMe balances. From an atomistic point of view, the formation of the closed conformer is favoured by polarisation and charge-transfer effects on the amide bond across all balances and is counterbalanced by a de-pyramidalisation of the amide nitrogen. Moreover, focusing on the oxygen of the amide carbonyl and the α-carbon of the remaining carbonyl group, electrostatics have a major role in the formation of the closed conformer, which goes against the well-known n-π* interaction orbital overlap concept. However, REG-IQF shows that exchange-correlation energies overtake electrostatics for all the 2-Y balances when working with fragments around the carbonyl groups, while they act on par with electrostatics for the 1-OMe and 1-NMe2. REG-IQF also shows that exchange-correlation energies in the 2-Y balance are correlated to the inductive electron-donating and -withdrawing trends on aromatic groups. We demonstrate that methods such as REG-IQA/IQF can help with the fine-tuning of molecular balances prior to the experiment and that the energies that govern the probed interactions are highly dependent on the atoms and functional groups involved.
Collapse
Affiliation(s)
- Fabio Falcioni
- Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Sophie Bennett
- Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Pallas Stroer-Jarvis
- Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Paul L A Popelier
- Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, UK
| |
Collapse
|
13
|
Ferrer M, Alkorta I, Elguero J, Oliva-Enrich JM. Capture of CO 2 by Melamine Derivatives: A DFT Study Combining the Relative Energy Gradient Method with an Interaction Energy Partitioning Scheme. J Phys Chem A 2024; 128:1288-1296. [PMID: 38351470 PMCID: PMC10895662 DOI: 10.1021/acs.jpca.3c08412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/26/2023] [Revised: 01/26/2024] [Accepted: 01/31/2024] [Indexed: 02/23/2024]
Abstract
A theoretical study of the interaction between melamine and CO2 was carried out using density functional theory (DFT) with the B3LYP-D3(BJ)/aug-cc-pVTZ level of theory. The presence of anions interacting with melamine transforms the weakly bonded tetrel complexes into adducts. Thus, melamine acts as an FLP (frustrated Lewis pair) with acid groups (NHs as hydrogen bond donors) and a base group (N of the triazine ring). The application of the relative energy gradient formalism (REG) along the reaction coordinate has demonstrated that the ability of the melamine-anion systems to capture CO2 is linked to its capacity to polarize the CO2 molecule. These results have been confirmed by placing the melamine:CO2 complex in a uniform electric field with different strengths.
Collapse
Affiliation(s)
- Maxime Ferrer
- Instituto
de Química Médica (CSIC), Juan de la Cierva, 3, E-28006 Madrid, Spain
- PhD
Program in Theoretical Chemistry and Computational Modeling, Doctoral
School, Universidad Autónoma de Madrid, 28049 Madrid, Spain
| | - Ibon Alkorta
- Instituto
de Química Médica (CSIC), Juan de la Cierva, 3, E-28006 Madrid, Spain
| | - Jose Elguero
- Instituto
de Química Médica (CSIC), Juan de la Cierva, 3, E-28006 Madrid, Spain
| | - Josep M. Oliva-Enrich
- Instituto
de Química-Física “Blas-Cabrera” (CSIC), Serrano, 119, E-28006 Madrid, Spain
| |
Collapse
|
14
|
Ahirwar MB, Khire SS, Gadre SR, Deshmukh MM. Hydrogen bond energy estimation (H-BEE) in large molecular clusters: A Python program for quantum chemical investigations. J Comput Chem 2024; 45:274-283. [PMID: 37792345 DOI: 10.1002/jcc.27237] [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: 09/13/2023] [Accepted: 09/18/2023] [Indexed: 10/05/2023]
Abstract
A procedure, derived from the fragmentation-based molecular tailoring approach (MTA), has been proposed and extensively applied by Deshmukh and Gadre for directly estimating the individual hydrogen bond (HB) energies and cooperativity contributions in molecular clusters. However, the manual fragmentation and high computational cost of correlated quantum chemical methods make the application of this method to large molecular clusters quite formidable. In this article, we report an in-house developed software for automated hydrogen bond energy estimation (H-BEE) in large molecular clusters. This user-friendly software is essentially written in Python and executed on a Linux platform with the Gaussian package at the backend. Two approximations to the MTA-based procedure, viz. the first spherical shell (SS1) and the Fragments-in-Fragments (Frags-in-Frags), enabling cost-effective, automated evaluation of HB energies and cooperativity contributions, are also implemented in this software. The software has been extensively tested on a variety of molecular clusters and is expected to be of immense use, especially in conjunction with correlated methods such as MP2, CCSD(T), and so forth.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| | | | - Shridhar R Gadre
- Department of Scientific Computing, Modelling & Simulation, Savitribai Phule Pune University, Pune, India
- Department of Chemistry, Savitribai Phule Pune University, Pune, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| |
Collapse
|
15
|
Sahre MJ, von Rudorff GF, Marquetand P, von Lilienfeld OA. Transferability of atomic energies from alchemical decomposition. J Chem Phys 2024; 160:054106. [PMID: 38341696 DOI: 10.1063/5.0187298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2023] [Accepted: 01/09/2024] [Indexed: 02/13/2024] Open
Abstract
We study alchemical atomic energy partitioning as a method to estimate atomization energies from atomic contributions, which are defined in physically rigorous and general ways through the use of the uniform electron gas as a joint reference. We analyze quantitatively the relation between atomic energies and their local environment using a dataset of 1325 organic molecules. The atomic energies are transferable across various molecules, enabling the prediction of atomization energies with a mean absolute error of 23 kcal/mol, comparable to simple statistical estimates but potentially more robust given their grounding in the physics-based decomposition scheme. A comparative analysis with other decomposition methods highlights its sensitivity to electrostatic variations, underlining its potential as a representation of the environment as well as in studying processes like diffusion in solids characterized by significant electrostatic shifts.
Collapse
Affiliation(s)
- Michael J Sahre
- Vienna Doctoral School in Chemistry (DoSChem) and Institute of Theoretical Chemistry and Faculty of Physics, University of Vienna, 1090 Vienna, Austria
| | - Guido Falk von Rudorff
- Department of Chemistry, University Kassel, Heinrich-Plett-Str.40, 34132 Kassel, Germany
- Center for Interdisciplinary Nanostructure Science and Technology (CINSaT), Heinrich-Plett-Straße 40, 34132 Kassel, Germany
| | - Philipp Marquetand
- Faculty of Chemistry, Institute of Theoretical Chemistry, University of Vienna, Währinger Str. 17, 1090 Vienna, Austria
| | - O Anatole von Lilienfeld
- Vienna Doctoral School in Chemistry (DoSChem) and Institute of Theoretical Chemistry and Faculty of Physics, University of Vienna, 1090 Vienna, Austria
- Chemical Physics Theory Group, Department of Chemistry, University of Toronto, St. George Campus, Toronto, M5S 3H6 Ontario, Canada
- Department of Materials Science and Engineering, University of Toronto, St. George Campus, Toronto, M5S 3E4 Ontario, Canada
- Vector Institute for Artificial Intelligence, Toronto, M5S 1M1 Ontario, Canada
- ML Group, Technische Universität Berlin and Institute for the Foundations of Learning and Data, 10587 Berlin, Germany
- Berlin Institute for the Foundations of Learning and Data, 10587 Berlin, Germany
- Department of Physics, University of Toronto, St. George Campus, Toronto, M5S 1A7 Ontario, Canada
| |
Collapse
|
16
|
Deltsov ID, Ananyev IV, Meerov DB, Fershtat LL. Expanding the Limits of Organic Energetic Materials: High-Performance Alliance of 1,3,4-Thiadiazole and Furazan Scaffolds. J Org Chem 2024; 89:174-182. [PMID: 38079563 DOI: 10.1021/acs.joc.3c01858] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2024]
Abstract
A majority of known and newly synthesized energetic materials comprise polynitrogen or nitrogen-oxygen heterocycles with various explosophores. However, available structural combinations of these organic scaffolds are finite and are about to reach their limits. Herein, we present the design and synthesis of a series of sulfur-containing polyazole structures comprising 1,3,4-thiadiazole and furazan rings linked by C-C bonds and enriched with energetic nitro and azo functionalities. In terms of detonation performance, all synthesized 1,3,4-thiadiazole-furazan assemblies (D = 7.7-7.9 km s-1; P = 26-28 GPa) lie between the powerful explosive TATB (D = 8.0 km s-1; P = 31 GPa) and melt-cast material TNT (D = 6.9 km s-1; P = 23 GPa). In the synthesized series, azo-bridged derivative 5 seems to be most practically interesting, as it combines a relatively high energetic performance (D = 7.9 km s-1; P = 28 GPa), a very high thermal stability (271 °C), and insensitivity to friction. By these functional properties, 5 outperforms the benchmark heat-resistant explosive hexanitrostilbene (HNS). To the best of our knowledge, this is the first example of an energetic alliance of furazan and 1,3,4-thiadiazole scaffolds and a rare case of sulfur-containing high-energy materials, which can certainly be considered as an evolutionary step in energetic materials science.
Collapse
Affiliation(s)
- Ilya D Deltsov
- N. D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, Moscow 119991, Russian Federation
| | - Ivan V Ananyev
- N. S. Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, GSP-1, Leninsky Prospect, 31, Moscow 119991, Russian Federation
| | - Dmitry B Meerov
- N. N. Semenov Federal Research Centre for Chemical Physics, Russian Academy of Sciences, 4 Kosygin Strasse, Moscow 119991, Russian Federation
| | - Leonid L Fershtat
- N. D. Zelinsky Institute of Organic Chemistry, Russian Academy of Sciences, Leninsky Prospect 47, Moscow 119991, Russian Federation
- National Research University Higher School of Economics, Myasnitskaya strasse, 20, Moscow 101000, Russian Federation
| |
Collapse
|
17
|
Goli M, Bressanini D, Shahbazian S. On the nature of the two-positron bond: evidence for a novel bond type. Phys Chem Chem Phys 2023; 25:29531-29547. [PMID: 37905569 DOI: 10.1039/d3cp03003b] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The nature of the newly proposed two-positron bond in (PsH)2, which is composed of two protons, four electrons and two positrons, is considered in this contribution. The study is done at the multi-component-Hartree-Fock (MC-HF) and the Diffusion Monte Carlo (DMC) levels of theory by comparing ab initio data, analyzing the spatial structure of the DMC wavefunction, and applying the multi-component quantum theory of atoms in molecules and the two-component interacting quantum atoms energy partitioning schemes to the MC-HF wavefunction. The analysis demonstrates that (PsH)2 to a good approximation may be conceived of as two slightly perturbed PsH atoms, bonded through a two-positron bond. In contrast to the usual two-electron bonds, the positron exchange phenomenon is quite marginal in the considered two-positron bond. The dominant stabilizing mechanism of bonding is a novel type of classical electrostatic interaction between the positrons, which are mainly localized between nuclei, and the surrounding electrons. To emphasize its uniqueness, this mechanism of bonding is proposed to be called gluonic which has also been previously identified as the main driving mechanism behind formation of the one-positron bond in [H-,e+,H-]. We conclude that the studied two-positron bond should not be classified as a covalent bond and it must be seen as a brand-new type of bond, foreign to the electronic bonding modes discovered so far in the purely electronic systems.
Collapse
Affiliation(s)
- Mohammad Goli
- School of Nano Science, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran.
| | - Dario Bressanini
- Dipartimento di Scienza e Alta Tecnologia, Università dell'Insubria, Como, Italy.
| | - Shant Shahbazian
- Department of Physics, Shahid Beheshti University, Evin, Tehran 19839-69411, Iran.
| |
Collapse
|
18
|
Vásquez-Espinal A, Pino-Rios R. Strong carbon - noble gas covalent bond and fluxionality in hypercoordinate compounds. Phys Chem Chem Phys 2023; 25:27468-27474. [PMID: 37800185 DOI: 10.1039/d3cp03576j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/07/2023]
Abstract
Thermodynamic, kinetic, and chemical bonding analysis at the coupled cluster level has been carried out for a series of hypercoordinated carbon compounds with formula CH4Ng2+ (Ng = He-Rn). Results show that these compounds could be stable at room temperature and Born-Oppenheimer molecular dynamics simulations (BOMD) in conjunction with activation energies indicate high kinetic stability. In addition, all chemical bonding descriptors agree with a strong C-Ng covalent bond and a bonding pattern similar to that of CH5+. Finally, BOMD simulations showed that these compounds are fluxional, with a continuous formation/breaking of H-H and C-H bonds. To the best of the authors' knowledge, these results represent the first series of fluxional compounds possessing a covalent bond between a main group element and a noble gas atom.
Collapse
Affiliation(s)
- Alejandro Vásquez-Espinal
- Química y Farmacia, Facultad de Ciencias de la Salud, Universidad Arturo Prat. Casilla 121, Iquique 1100000, Chile.
| | - Ricardo Pino-Rios
- Química y Farmacia, Facultad de Ciencias de la Salud, Universidad Arturo Prat. Casilla 121, Iquique 1100000, Chile.
- Instituto de Estudios de la Salud, Universidad Arturo Prat, Iquique, 1100000, Chile
| |
Collapse
|
19
|
Anisimov AA, Ananyev IV. Electron density-based protocol to recover the interacting quantum atoms components of intermolecular binding energy. J Chem Phys 2023; 159:124113. [PMID: 38127385 DOI: 10.1063/5.0167874] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2023] [Accepted: 09/11/2023] [Indexed: 12/23/2023] Open
Abstract
A new approach for obtaining interacting quantum atoms-defined components of binding energy of intermolecular interactions, which bypasses the use of standard six-dimensional integrals and two-particle reduced density matrix (2-RDM) reconstruction, is proposed. To examine this approach, three datasets calculated within the density functional theory framework using the def2-TZVP basis have been explored. The first two, containing 53 weakly bound bimolecular associates and 13 molecular clusters taken from the crystal, were used in protocol refinement, and the third one containing other 20 bimolecular and three cluster systems served as a validation reference. In addition, to verify the performance of the proposed approach on an exact 2-RDM, calculations within the coupled cluster formalism were performed for part of the first set systems using the cc-pVTZ basis set. The process of optimization of the proposed parametric model is considered, and the role of various energy contributions in the formation of non-covalent interactions is discussed with regard to the obtained trends.
Collapse
Affiliation(s)
- Aleksei A Anisimov
- A.N. Nesmeyanov Institute of Organoelement Compounds, Russian Academy of Sciences, Vavilova Str. 28, GSP-1, Moscow 119334, Russian Federation
- National Research University Higher School of Economics, Myasnitskaya Str. 20, Moscow 101000, Russian Federation
| | - Ivan V Ananyev
- N.S. Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, GSP-1, Leninsky prospect 31, Moscow 119991, Russian Federation
| |
Collapse
|
20
|
Triestram L, Falcioni F, Popelier PLA. Interacting Quantum Atoms and Multipolar Electrostatic Study of XH···π Interactions. ACS OMEGA 2023; 8:34844-34851. [PMID: 37779962 PMCID: PMC10535255 DOI: 10.1021/acsomega.3c04149] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 06/12/2023] [Accepted: 08/31/2023] [Indexed: 10/03/2023]
Abstract
The interaction energies of nine XH···π (X = C, N, and O) benzene-containing van der Waals complexes were analyzed, at the atomic and fragment levels, using QTAIM multipolar electrostatics and the energy partitioning method interacting quantum atoms/fragment (IQA/IQF). These descriptors were paired with the relative energy gradient method, which solidifies the connection between quantum mechanical properties and chemical interpretation. This combination provides a precise understanding, both qualitative and quantitative, of the nature of these interactions, which are ubiquitous in biochemical systems. The formation of the OH···π and NH···π systems is electrostatically driven, with the Qzz component of the quadrupole moment of the benzene carbons interacting with the charges of X and H in XH. There is the unexpectedly intramonomeric role of X-H (X = O, N) where its electrostatic energy helps the formation of the complex and its covalent energy thwarts it. However, the CH···π interaction is governed by exchange-correlation energies, thereby establishing a covalent character, as opposed to the literature's designation as a noncovalent interaction. Moreover, dispersion energy is relevant, statically and in absolute terms, but less relevant compared to other energy components in terms of the formation of the complex. Multipolar electrostatics are similar across all systems.
Collapse
Affiliation(s)
- Lena Triestram
- Department of Chemistry, University
of Manchester, Manchester M13 9PL, Great
Britain
| | - Fabio Falcioni
- Department of Chemistry, University
of Manchester, Manchester M13 9PL, Great
Britain
| | - Paul L. A. Popelier
- Department of Chemistry, University
of Manchester, Manchester M13 9PL, Great
Britain
| |
Collapse
|
21
|
Cador A, Tognetti V, Joubert L, Popelier PLA. Aza-Michael Addition in Explicit Solvent: A Relative Energy Gradient-Interacting Quantum Atoms Study. Chemphyschem 2023:e202300529. [PMID: 37728125 DOI: 10.1002/cphc.202300529] [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: 07/27/2023] [Revised: 09/09/2023] [Accepted: 09/17/2023] [Indexed: 09/21/2023]
Abstract
Aza-Michael additions are key reactions in organic synthesis. We investigate, from a theoretical and computational point of view, several examples ranging from weak to strong electrophiles in dimethylsulfoxide treated as explicit solvent. We use the REG-IQA method, which is a quantum topological energy decomposition (Interacting Quantum Atoms, IQA) coupled to a chemical-interpretation calculator (Relative Energy Gradient, REG). We focus on the rate-limiting addition step in order to unravel the different events taking place in this step, and understand the influence of solvent on the reaction, with an eye on predicting the Mayr electrophilicity. For the first time, a link is established between an REG-IQA analysis and experimental values.
Collapse
Affiliation(s)
- Aël Cador
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, 76821, Mont St, Aignan Cedex, France
| | - Vincent Tognetti
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, 76821, Mont St, Aignan Cedex, France
| | - Laurent Joubert
- Normandy Univ., COBRA UMR 6014 & FR 3038, Université de Rouen, INSA Rouen, CNRS, 1 rue Tesnière, 76821, Mont St, Aignan Cedex, France
| | - Paul L A Popelier
- Department of Chemistry, The University of Manchester, Manchester, M13 9PL, Great Britain
| |
Collapse
|
22
|
Pérez-Bitrián A, Munárriz J, Sturm JS, Wegener D, Krause KB, Wiesner A, Limberg C, Riedel S. Further Perspectives on the Teflate versus Fluoride Analogy: The Case of a Co(II) Pentafluoroorthotellurate Complex. Inorg Chem 2023; 62:12947-12953. [PMID: 37505485 DOI: 10.1021/acs.inorgchem.3c01730] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/29/2023]
Abstract
The pentafluoroorthotellurate group (teflate, OTeF5) is considered as a bulky analogue of fluoride, yet its coordination behavior in transition metal complexes is not fully understood. By reaction of [CoCl4]2- and neat ClOTeF5, we synthesized the first cobalt teflate complex, [Co(OTeF5)4]2-, which exhibits moisture-resistant Co-OTeF5 bonds. Through a combined experimental and theoretical (DFT and NEVPT2) study, the properties and electronic structure of this species have been investigated. It exhibits a distorted tetrahedral structure around the cobalt center and can be described as a d7 system with a quartet (S = 3/2) ground state. A comparative bonding analysis of the (pseudo)tetrahedral [CoX4]2- anions (X = OTeF5, F, Cl) revealed that the strength of the Co-X interaction is similar in the three cases, being the strongest in [Co(OTeF5)4]2-. In addition, an analysis of the charge of the Co center reinforced the similar electron-withdrawing properties of the teflate and fluoride ligands. Therefore, the [Co(OTeF5)4]2- anion constitutes an analogue of the polymeric [CoF4]2- in terms of electronic properties, but with a monomeric structure.
Collapse
Affiliation(s)
- Alberto Pérez-Bitrián
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/36, 14195 Berlin, Germany
| | - Julen Munárriz
- Departamento de Química Física y Analítica, Universidad de Oviedo, Campus Universitario de El Cristo, Julián Clavería no. 8, 33006 Oviedo, Spain
| | - Johanna S Sturm
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/36, 14195 Berlin, Germany
| | - Daniel Wegener
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/36, 14195 Berlin, Germany
| | - Konstantin B Krause
- Institut für Chemie, Humboldt-Universität zu Berlin, Brook-Taylor-Straße 2, 12489 Berlin, Germany
| | - Anja Wiesner
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/36, 14195 Berlin, Germany
| | - Christian Limberg
- Institut für Chemie, Humboldt-Universität zu Berlin, Brook-Taylor-Straße 2, 12489 Berlin, Germany
| | - Sebastian Riedel
- Fachbereich Biologie, Chemie, Pharmazie, Institut für Chemie und Biochemie - Anorganische Chemie, Freie Universität Berlin, Fabeckstraße 34/36, 14195 Berlin, Germany
| |
Collapse
|
23
|
Gimferrer M, Salvador P. Exact decompositions of the total KS-DFT exchange-correlation energy into one- and two-center terms. J Chem Phys 2023; 158:234105. [PMID: 37326158 DOI: 10.1063/5.0142778] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2023] [Accepted: 05/25/2023] [Indexed: 06/17/2023] Open
Abstract
In the so-called Interacting Quantum Atoms (IQA) approach, the molecular energy is numerically decomposed as a sum of atomic and diatomic contributions. While proper formulations have been put forward for both Hartree-Fock and post-Hartree-Fock wavefunctions, this is not the case for the Kohn-Sham density functional theory (KS-DFT). In this work, we critically analyze the performance of two fully additive approaches for the IQA decomposition of the KS-DFT energy, namely, the one from Francisco et al., which uses atomic scaling factors, and that from Salvador and Mayer based upon the bond order density (SM-IQA). Atomic and diatomic exchange-correlation (xc) energy components are obtained for a molecular test set comprising different bond types and multiplicities and along the reaction coordinate of a Diels-Alder reaction. Both methodologies behave similarly for all systems considered. In general, the SM-IQA diatomic xc components are less negative than the Hartree-Fock ones, which is in good agreement with the known effect of electron correlation upon (most) covalent bonds. In addition, a new general scheme to minimize the numerical error of the sum of two-electron energy contributions (i.e., Coulomb and exact exchange) in the framework of overlapping atoms is described in detail.
Collapse
Affiliation(s)
- Martí Gimferrer
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| | - Pedro Salvador
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, Maria Aurèlia Capmany 69, 17003 Girona, Catalonia, Spain
| |
Collapse
|
24
|
Ahirwar MB, Deshmukh MM. Fragments-in-fragments method for efficient and reliable estimates of individual hydrogen bond energies in large molecular clusters. J Comput Chem 2023. [PMID: 37191018 DOI: 10.1002/jcc.27133] [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/2023] [Revised: 05/01/2023] [Accepted: 05/03/2023] [Indexed: 05/17/2023]
Abstract
The knowledge of individual hydrogen bond (HB) strength in molecular clusters is indispensable to get insights into the bulk properties of condensed systems. Recently, we have developed the molecular tailoring approach based (MTA-based) method for the estimation of individual HB energy in molecular clusters. However, the direct use of this MTA-based method to large molecular clusters becomes progressively difficult with the increase in the size of a cluster. To overcome this caveat, herein, we propose the use of linear scaling method (such as the original MTA method) for the estimation of single-point (SP) energies of large-sized parent molecular cluster and their respective fragments. Because the fragments of the MTA-based method, for the estimation of HB energy, are further fragmented, this proposed strategy is called as Fragments-in-Fragments (Frags-in-Frags) method. The SP energies of fragments and parent cluster calculated by the Frags-in-Frags approach were utilized to estimate the individual HB energy. The estimated individual HB energies, in various molecular clusters, by Frags-in-Frags method are found to be in excellent linear agreement with their MTA-based counterparts (R2 = 0.9975 of 348 data points). The difference being less than 0.5 kcal/mol in most of the cases. Furthermore, RMSD is 0.43 kcal/mol, MAE is 0.33 kcal/mol, and the standard deviation is 0.44 kcal/mol. Importantly, the Frags-in-Frags method not only enables the reliable estimation of HB energy in large molecular clusters but also requires less computational time and can be possible even with off-the-shelf hardware.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar, India
| |
Collapse
|
25
|
Ahirwar MB, Gadre SR, Deshmukh MM. On the Short-Range Nature of Cooperativity in Hydrogen-Bonded Large Molecular Clusters. J Phys Chem A 2023; 127:4394-4406. [PMID: 37186960 DOI: 10.1021/acs.jpca.3c00359] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/17/2023]
Abstract
The variation in the hydrogen bond (HB) strength has considerable consequences on the physicochemical properties of molecular clusters. Such a variation mainly arises due to the cooperative/anti-cooperative networking effect of neighboring molecules connected by HBs. In the present work, we systematically study the effect of neighboring molecules on the strength of an individual HB and the respective cooperativity contribution toward each of them in a variety of molecular clusters. For this purpose, we propose a use of a small model of a large molecular cluster called the spherical shell-1 (SS1) model. This SS1 model is constructed by placingg the spheres of an appropriate radius centered on X and Y atoms of the X-H···Y HB under consideration. The molecules falling within these spheres constitute the SS1 model. Utilizing this SS1 model, the individual HB energies are calculated within the molecular tailoring approach-based framework and the results are compared with their actual counterparts. It is found that the SS1 is a reasonably good model of large molecular clusters, providing 81-99% of the total HB energy estimated using the actual molecular clusters. This in turn suggests that the maximum cooperativity contribution toward a particular HB is due to the fewer number of molecules (in the SS1 model) directly interacting with two molecules involved in its formation. We further demonstrate that the remaining part of the energy or cooperativity (∼1 to 19%) is captured by the molecules falling in the second spherical shell (SS2) centered on the hetero-atom of the molecules in the SS1 model. The effect of increasing size of a cluster on the strength of a particular HB, calculated by the SS1 model, is also investigated. The calculated value of the HB energy remains unchanged with the increase in the size of a cluster, emphasizing the short-ranged nature of the HB cooperativity in neutral molecular clusters.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| | - Shridhar R Gadre
- Department of Scientific Computing, Modelling, and Simulation, Savitribai Phule Pune University, Pune 411007, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar 470003, India
| |
Collapse
|
26
|
Hercigonja M, Milovanović B, Etinski M, Petković M. Decorated crown ethers as selective ion traps: Solvent’s role in crown’s preference towards a specific ion. J Mol Liq 2023. [DOI: 10.1016/j.molliq.2023.121791] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/08/2023]
|
27
|
Sahre MJ, von Rudorff GF, von Lilienfeld OA. Quantum Alchemy Based Bonding Trends and Their Link to Hammett's Equation and Pauling's Electronegativity Model. J Am Chem Soc 2023; 145:5899-5908. [PMID: 36862462 DOI: 10.1021/jacs.2c13393] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023]
Abstract
We present an intuitive and general analytical approximation estimating the energy of covalent single and double bonds between participating atoms in terms of their respective nuclear charges with just three parameters, [EAB ≈ a - bZAZB + c(ZA7/3 + ZB7/3) ]. The functional form of our expression models an alchemical atomic energy decomposition between participating atoms A and B. After calibration, reasonably accurate bond dissociation energy estimates are obtained for hydrogen-saturated diatomics composed of p-block elements coming from the same row 2 ≤ n ≤ 4 in the periodic table. Corresponding changes in bond dissociation energies due to substitution of atom B by C can be obtained via simple formulas. While being of different functional form and origin, our model is as simple and accurate as Pauling's well-known electronegativity model. Analysis indicates that the model's response in covalent bonding to variation in nuclear charge is near-linear, which is consistent with Hammett's equation.
Collapse
Affiliation(s)
- Michael J Sahre
- Faculty of Physics, University of Vienna, Vienna, 1090, Austria.,Vienna Doctoral School in Chemistry (DoSChem), University of Vienna, Vienna, 1090, Austria
| | | | - O Anatole von Lilienfeld
- Vector Institute for Artificial Intelligence, Toronto, M5S 1M1, Canada.,Departments of Chemistry, Materials Science and Engineering, and Physics, University of Toronto, St. George Campus, Toronto, M5R 0A3, Canada.,Machine Learning Group, Technische Universität Berlin and Institute for the Foundations of Learning and Data, Berlin, 10587, Germany
| |
Collapse
|
28
|
The Ultrashort Spike-Ring Interaction in Substituted Iron Maiden Molecules. Molecules 2023; 28:molecules28052244. [PMID: 36903489 PMCID: PMC10004599 DOI: 10.3390/molecules28052244] [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: 02/06/2023] [Revised: 02/23/2023] [Accepted: 02/24/2023] [Indexed: 03/06/2023] Open
Abstract
The in forms of molecular iron maidens are known for their unique ultrashort interaction between the apical hydrogen atom or its small substituent and the surface of the benzene ring. It is generally believed that this forced ultrashort X⋯π contact is associated with high steric hindrance, which is responsible for specific properties of iron maiden molecules. The main aim of this article is to investigate the influence of significant charge enrichment or depletion of the benzene ring on the characteristics of the ultrashort C-X⋯π contact in iron maiden molecules. For this purpose, three strongly electron-donating (-NH2) or strongly electron-withdrawing (-CN) groups were inserted into the benzene ring of in-[34,10][7]metacyclophane and its halogenated (X = F, Cl, Br) derivatives. It is shown that, despite such extremely electron-donating or electron-accepting properties, the considered iron maiden molecules surprisingly reveal quite high resistance to changes in electronic properties.
Collapse
|
29
|
Goli M, Shahbazian S. MC-QTAIM analysis reveals an exotic bond in coherently quantum superposed malonaldehyde. Phys Chem Chem Phys 2023; 25:5718-5730. [PMID: 36744327 DOI: 10.1039/d2cp05499j] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
The proton between the two oxygen atoms of the malonaldehyde molecule experiences an effective double-well potential in which the proton's wavefunction is delocalized between the two wells. Herein we employ a state-of-the-art multi-component quantum theory of atoms in molecules partitioning scheme to obtain the molecular structure, i.e. atoms in molecules and bonding network, from the superposed ab initio wavefunctions of malonaldehyde. In contrast to the familiar clamped-proton portrayal of malonaldehyde, in which the proton forms a hydrogen basin, for the superposed states the hydrogen basin disappears and two novel hybrid oxygen-hydrogen basins appear instead, with an even distribution of the proton population between the two basins. The interaction between the hybrid basins is stabilizing thanks to an unprecedented mechanism. This involves the stabilizing classical Coulomb interaction of the one-proton density in one of the basins with one-electron density in the other basin. This stabilizing mechanism yields a bond foreign to the known bonding modes in chemistry.
Collapse
Affiliation(s)
- Mohammad Goli
- School of Nano Science, Institute for Research in Fundamental Sciences (IPM), Tehran 19395-5531, Iran.
| | - Shant Shahbazian
- Department of Physics, Shahid Beheshti University, Evin, Tehran, Iran.
| |
Collapse
|
30
|
Ahirwar MB, Deshmukh MM. Two-Step ONIOM Method for the Accurate Estimation of Individual Hydrogen Bond Energy in Large Molecular Clusters. J Phys Chem A 2023; 127:1219-1232. [PMID: 36705264 DOI: 10.1021/acs.jpca.2c08087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
The study of molecular clusters to understand the properties of condensed systems has been the subject of immense interest. To get insight into these properties, the knowledge of various noncovalent interactions present in these molecular clusters is indispensable. Our recently developed molecular tailoring approach-based (MTA-based) method for the estimation of the individual hydrogen bond (HB) energy in molecular clusters is useful for this purpose. However, the direct application of this MTA-based method becomes progressively difficult with the increase in the size of the cluster. This is because of the difficulty in the evaluation of single-point energy at the correlated level of theory. To overcome this caveat, herein, we propose a two-step method within the our own N-layer integrated molecular orbital molecular mechanics (ONIOM) framework. In this method, the HB energy evaluated by the MTA-based method employing the actual molecular cluster at a low Hartree-Fock (HF) level of theory is added to the difference in the HB energies evaluated by the MTA-based method, employing an appropriate small model system, called the shell-1 model, calculated at high (MP2) and low (HF) levels of theory. The shell-1 model of a large molecular cluster is made up of only a few molecules that are in direct contact (by a single HB) with the two molecules involved in the formation of an HB under consideration. We tested this proposed two-step ONIOM method to estimate the individual HB energies in various molecular clusters, viz., water (Wn, n = 10-16, 18 and 20), (H2O2)12, (H2O3)8, (NH3)n and strongly interacting (HF)15 and (HF)m(W)n clusters. Furthermore, these estimated individual HB energies by the ONIOM method are compared with those calculated by the MTA-based method using actual molecular clusters. The estimated individual HB energies by the ONIOM method, in all these clusters, are in excellent linear one-to-one agreement (R2 = 0.9996) with those calculated by the MTA-based method using actual molecular clusters. Furthermore, the small values of root-mean-square deviation (0.06), mean absolute error (0.04), |ΔEmax| (0.21) and Sε (0.06) suggest that this two-step ONIOM method is a pragmatic approach to provide accurate estimates of individual HB energies in large molecular clusters.
Collapse
Affiliation(s)
- Mini Bharati Ahirwar
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar470003, India
| | - Milind M Deshmukh
- Department of Chemistry, Dr. Harisingh Gour Vishwavidyalaya (A Central University), Sagar470003, India
| |
Collapse
|
31
|
Feitosa LF, Campos RB, Richter WE. Energetics and electronics of polar Diels-Alder reactions at the atomic level: QTAIM and IQA analyses of complete IRC paths. J Mol Graph Model 2023; 118:108326. [PMID: 36166996 DOI: 10.1016/j.jmgm.2022.108326] [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: 08/05/2022] [Revised: 09/05/2022] [Accepted: 09/06/2022] [Indexed: 11/28/2022]
Abstract
The mechanism of Diels-Alder reactions between cyclopentadiene and several cyanoethylenes was studied by means of Density Functional Theory calculations using QTAIM and IQA (Interacting Quantum Atoms) analyses along complete IRC paths. Each geometry from the IRC had its wavefunction computed and the topology of the electronic density for it was then evaluated. By means of IQA, the global energetic profile was partitioned among the various atoms in the molecule, providing insight into what atoms are the main ones responsible for the magnitude of the energy barriers. The (a)synchronicity of the reaction mechanisms featuring non-symmetrically substituted dienophiles was characterized, from QTAIM, by the electron densities and Laplacians over the LCP's as well as by the different atomic energy barriers obtained from IQA. The magnitude of the atomic barrier nicely explains the (a)synchronicity of the reaction mechanisms, and the degree of (a)synchronicity is nicely revealed by the difference between the earlier and later bond breaking and bond formations. The main conclusion is that important energetic and electronic changes are occurring before and after the position of the transition state structure, mainly for those asynchronous mechanisms, and although these provide essential insight into the reaction mechanism, most studies cannot assess this kind of information because they are focusing solely on reactants, transition states, and products. We advocate that the additional computational effort required for such analyses is more than compensated by the great amount of useful information it provides.
Collapse
Affiliation(s)
- Lucas Freitas Feitosa
- Department of Chemistry and Biology, Federal University of Technology - Paraná [UTFPR], 81.280-340, Curitiba PR, Brazil
| | - Renan Borsoi Campos
- Department of Chemistry and Biology, Federal University of Technology - Paraná [UTFPR], 81.280-340, Curitiba PR, Brazil
| | - Wagner Eduardo Richter
- Department of Chemistry, Federal University of Technology - Paraná [UTFPR], 84.017-220, Ponta Grossa PR, Brazil.
| |
Collapse
|
32
|
Job N, Thirumoorthy K. Chemical Bonding Perspective on Low-Lying SiC 4H 2 Isomers: Conceptual Quantum Chemical Views. J Phys Chem A 2022; 126:9366-9374. [PMID: 36512698 DOI: 10.1021/acs.jpca.2c06230] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Abstract
The nature of the chemical bonding in seven low-lying isomers of SiC4H2 is analyzed through quantum chemical concepts. Out of the seven, four isomers, 1-ethynyl-3-silacycloprop-1(2)-en-3-ylidene (1), diethynylsilylidene (2), 1-sila-1,2,3,4-pentatetraenylidene (4), and 1,3-butadiynylsilylidene (5), have already been identified in the laboratory. The other three isomers, 2-methylenesilabicyclo[1.1.0]but-1(3)-en-4-ylidene (3), 4-sila-2-methylenebicyclo[1.1.0]but-1(3)-en-4-ylidene (6), and 3-ethynyl-1-silapropadienylidene (7) remain elusive in the laboratory to date (J. Phys. Chem. A, 2020, 124, 987-1002). Deep insight into the characteristics of chemical bonding is explored with different bonding analysis tools. Quantum theory of atoms in molecules (QTAIM), interaction quantum atoms analysis, natural bond orbital analysis, adaptive natural density partitioning, electron localization function (ELF), Laplacian of electron density, energy decomposition analysis, atomic charge analysis, bond order analysis, and frontier molecular orbital analysis are employed in the present work to gain a better understanding of the chemical bonding perspective in SiC4H2 isomers. Different quantum chemical topology approaches (QTAIM, ELF, and Laplacian of electron density) are employed to complement each other. The obtained results dictate that the lone pair of the silicon atom participate in delocalization and influences the structural stability of isomers.
Collapse
Affiliation(s)
- Nisha Job
- Department of Chemistry, School of Advanced Sciences, Vellore Institute of Technology, Vellore 632014, Tamil Nadu, India
| | - Krishnan Thirumoorthy
- Department of Chemistry, School of Advanced Sciences, Vellore Institute of Technology, Vellore 632014, Tamil Nadu, India
| |
Collapse
|
33
|
Nature of Beryllium, Magnesium, and Zinc Bonds in Carbene⋯MX 2 (M = Be, Mg, Zn; X = H, Br) Dimers Revealed by the IQA, ETS-NOCV and LED Methods. Int J Mol Sci 2022; 23:ijms232314668. [PMID: 36498996 PMCID: PMC9738500 DOI: 10.3390/ijms232314668] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/24/2022] [Revised: 11/21/2022] [Accepted: 11/22/2022] [Indexed: 11/25/2022] Open
Abstract
The nature of beryllium−, magnesium− and zinc−carbene bonds in the cyclopropenylidene⋯MX2 (M = Be, Mg, Zn; X = H, Br) and imidazol-2-ylidene⋯MBr2 dimers is investigated by the joint use of the topological QTAIM-based IQA decomposition scheme, the molecular orbital-based ETS-NOCV charge and energy decomposition method, and the LED energy decomposition approach based on the state-of-the-art DLPNO-CCSD(T) method. All these methods show that the C⋯M bond strengthens according to the following order: Zn < Mg << Be. Electrostatics is proved to be the dominant bond component, whereas the orbital component is far less important. It is shown that QTAIM/IQA underestimates electrostatic contribution for zinc bonds with respect to both ETS-NOCV and LED schemes. The σ carbene→MX2 donation appears to be much more important than the MX2→ carbene back-donation of π symmetry. The substitution of hydrogen atoms by bromine (X in MX2) strengthens the metal−carbene bond in all cases. The physical origin of rotational barriers has been unveiled by the ETS-NOCV approach.
Collapse
|
34
|
The IQA Energy Partition in a Drug Design Setting: A Hepatitis C Virus RNA-Dependent RNA Polymerase (NS5B) Case Study. Pharmaceuticals (Basel) 2022; 15:ph15101237. [PMID: 36297349 PMCID: PMC9609620 DOI: 10.3390/ph15101237] [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: 09/09/2022] [Revised: 09/30/2022] [Accepted: 10/01/2022] [Indexed: 11/05/2022] Open
Abstract
The interaction of the thumb site II of the NS5B protein of hepatitis C virus and a pair of drug candidates was studied using a topological energy decomposition method called interacting quantum atoms (IQA). The atomic energies were then processed by the relative energy gradient (REG) method, which extracts chemical insight by computation based on minimal assumptions. REG reveals the most important IQA energy contributions, by atom and energy type (electrostatics, sterics, and exchange–correlation), that are responsible for the behaviour of the whole system, systematically from a short-range ligand–pocket interaction until a distance of approximately 22 Å. The degree of covalency in various key interatomic interactions can be quantified. No exchange–correlation contribution is responsible for the changes in the energy profile of both pocket–ligand systems investigated in the ligand–pocket distances equal to or greater than that of the global minimum. Regarding the hydrogen bonds in the system, a “neighbour effect” was observed thanks to the REG method, which states that a carbon atom would rather not have its covalent neighbour oxygen form a hydrogen bond. The combination of IQA and REG enables the automatic identification of the pharmacophore in the ligands. The coarser Interacting Quantum Fragments (IQF) enables the determination of which amino acids of the pocket contribute most to the binding and the type of energy of said binding. This work is an example of the contribution topological energy decomposition methods can make to fragment-based drug design.
Collapse
|
35
|
Zapata-Acevedo CA, Guevara-Vela JM, Popelier PLA, Rocha Rinza T. Binding Energy Partition of Promising IRAK-4 Inhibitor (Zimlovisertib) for the Treatment of COVID-19 Pneumonia. Chemphyschem 2022; 23:e202200455. [PMID: 36044560 PMCID: PMC9538207 DOI: 10.1002/cphc.202200455] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Revised: 08/19/2022] [Indexed: 01/05/2023]
Abstract
The technique of Fragment-Based Drug Design (FBDD) considers the interactions of different moieties of molecules with biological targets for the rational construction of potential drugs. One basic assumption of FBDD is that the different functional groups of a ligand interact with a biological target in an approximately additive, that is, independent manner. We investigated the interactions of different fragments of ligands and Interleukin-1 Receptor-Associated Kinase 4 (IRAK-4) throughout the FBDD design of Zimlovisertib, a promising anti-inflammatory, currently in trials to be used for the treatment of COVID-19 pneumonia. We utilised state-of-the-art methods of wave function analyses mainly the Interacting Quantum Atoms (IQA) energy partition for this purpose. By means of IQA, we assessed the suitability of every change to the ligand in the five stages of FBDD which led to Zimlovisertib on a quantitative basis. We determined the energetics of the interaction of different functional groups in the ligands with the IRAK-4 protein target and thereby demonstrated the adequacy (or lack thereof) of the changes made across the design of this drug. This analysis permits to verify whether a given alteration of a prospective drug leads to the intended tuning of non-covalent interactions with its protein objective. Overall, we expect that the methods exploited in this paper will prove valuable in the understanding and control of chemical modifications across FBDD processes.
Collapse
Affiliation(s)
- César Arturo Zapata-Acevedo
- Tecnologico de Monterrey: Instituto Tecnologico y de Estudios Superiores de MonterreyChemistryAv. Carlos Lazo 100Santa Fe, La Loma01389Álvaro ObregónMEXICO
| | | | - Paul L. A. Popelier
- UoM: The University of ManchesterChemistryOxford RoadM13 9PLManchesterUNITED KINGDOM
| | - Tomás Rocha Rinza
- Institute Of Chemistry, National Autonomous University of MexicoDepartment of Physical ChemistryCircuito Exterior, Ciudad Universitaria04510Mexico CityMEXICO
| |
Collapse
|
36
|
Non-covalent interactions from a Quantum Chemical Topology perspective. J Mol Model 2022; 28:276. [PMID: 36006513 PMCID: PMC9411098 DOI: 10.1007/s00894-022-05188-7] [Citation(s) in RCA: 18] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Accepted: 02/07/2022] [Indexed: 11/12/2022]
Abstract
About half a century after its little-known beginnings, the quantum topological approach called QTAIM has grown into a widespread, but still not mainstream, methodology of interpretational quantum chemistry. Although often confused in textbooks with yet another population analysis, be it perhaps an elegant but somewhat esoteric one, QTAIM has been enriched with about a dozen other research areas sharing its main mathematical language, such as Interacting Quantum Atoms (IQA) or Electron Localisation Function (ELF), to form an overarching approach called Quantum Chemical Topology (QCT). Instead of reviewing the latter’s role in understanding non-covalent interactions, we propose a number of ideas emerging from the full consequences of the space-filling nature of topological atoms, and discuss how they (will) impact on interatomic interactions, including non-covalent ones. The architecture of a force field called FFLUX, which is based on these ideas, is outlined. A new method called Relative Energy Gradient (REG) is put forward, which is able, by computation, to detect which fragments of a given molecular assembly govern the energetic behaviour of this whole assembly. This method can offer insight into the typical balance of competing atomic energies both in covalent and non-covalent case studies. A brief discussion on so-called bond critical points is given, highlighting concerns about their meaning, mainly in the arena of non-covalent interactions.
Collapse
|
37
|
Salazar-Lozas H, Guevara-Vela JM, Pendás ÁM, Francisco E, Rocha-Rinza T. Partition of the electronic energy of the PM7 method via the interacting quantum atoms approach. Phys Chem Chem Phys 2022; 24:19521-19530. [PMID: 35938407 DOI: 10.1039/d2cp02013k] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Partitions of the electronic energy such as that provided by the Interacting Quantum Atoms (IQA) approach have given valuable insights for numerous chemical systems and processes. Unfortunately, this kind of analysis may involve the integration of scalar fields over very irregular volumes, a condition which leads to a large and often prohibitive computational effort. These circumstances have limited the use of these energy partitions to systems comprising a few tens of atoms at most. On the other hand, semiempirical methods have proved useful in the study of systems of several thousands of atoms. Therefore, the goal of this work is to carry out partitions of the semiempirical method PM7 in compliance with the IQA approach. For this purpose, we computed one- and two-atomic energetic contributions whose sum equals the PM7 electronic energy. We illustrate how one might exploit the partition of electronic energies computed via the PM7 method by considering small organic and inorganic molecules and the energetics of individual hydrogen bond interactions within several water clusters which include (H2O)30, (H2O)50 and (H2O)100. We also considered the solvation of the amphiphilic caprylate anion to exemplify how to exploit the energy partition proposed in this paper. Overall, this investigation shows how the approach put forward herein might give further insights of the interactions occurring within complex systems in physical and biological chemistry.
Collapse
Affiliation(s)
- Hugo Salazar-Lozas
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Alcaldía Coyoacán C.P. 04510, Ciudad de México, Mexico.
| | | | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Universidad de Oviedo, Julián Claveria 8, 33006, Oviedo, Spain
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Universidad de Oviedo, Julián Claveria 8, 33006, Oviedo, Spain
| | - Tomás Rocha-Rinza
- Instituto de Química, Universidad Nacional Autónoma de México, Circuito Exterior, Ciudad Universitaria, Alcaldía Coyoacán C.P. 04510, Ciudad de México, Mexico.
| |
Collapse
|
38
|
Jabłoński M. The physical nature of the ultrashort spike-ring interaction in iron maiden molecules. J Comput Chem 2022; 43:1206-1220. [PMID: 35593685 DOI: 10.1002/jcc.26879] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/25/2022] [Revised: 04/11/2022] [Accepted: 04/13/2022] [Indexed: 11/08/2022]
Abstract
The so-called 'iron maiden' molecules belong to one of the most interesting subgroups of cyclophanes due to the presence of the ultrashort interaction between the CX apical bond and the benzene ring. This article presents an in-depth theoretical study of 16 'iron maiden' molecules, in which X = H, F, Cl or Br and the side chains are of various lengths and types: CSC, CSCC, CCC, and CCCC. It is shown that the H → F → Cl → Br substitution leads to a significant expansion of the 'iron maiden' molecule. Shorter chains lead to more pronounced effects, while insertion of sulfur atoms into the side chains lowers them. Structural changes are associated with an increase in energetic destabilization of X. Moreover, unlike for H, in the case of X = halogen, the out → in isomerization is energetically disadvantageous. The 'iron maiden' molecules are characterized by the presence of only three X⋯CAr bond paths. Particularly noteworthy are unusually large (even up to 32) values of the X⋯CAr bond ellipticity, which results from flat electron density distribution. The X⋯π interaction in each of the investigated 'iron maiden' molecule turned out to be multi-center, stabilizing and almost purely covalent in nature as indicated by the definitely dominant percentage (94.8%-101.6%) of the exchange-correlation energy. The spatial hindrance within the 'iron maiden' molecules appears to be not so much due to the X⋯π repulsion, but due to unfavorable steric interactions between X and the CC side bonds. It is also confirmed that some CH⋯HC interactions in aliphatic chains can be very weakly stabilizing.
Collapse
Affiliation(s)
- Mirosław Jabłoński
- Faculty of Chemistry, Nicolaus Copernicus University in Toruń, Toruń, Poland
| |
Collapse
|
39
|
Martín Pendás Á, Francisco E. The role of references and the elusive nature of the chemical bond. Nat Commun 2022; 13:3327. [PMID: 35680893 PMCID: PMC9184482 DOI: 10.1038/s41467-022-31036-6] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2022] [Accepted: 05/30/2022] [Indexed: 11/16/2022] Open
Abstract
Chemical bonding theory is of utmost importance to chemistry, and a standard paradigm in which quantum mechanical interference drives the kinetic energy lowering of two approaching fragments has emerged. Here we report that both internal and external reference biases remain in this model, leaving plenty of unexplored territory. We show how the former biases affect the notion of wavefunction interference, which is purportedly recognized as the most basic bonding mechanism. The latter influence how bonding models are chosen. We demonstrate that the use of real space analyses are as reference-less as possible, advocating for their use. Delocalisation emerges as the reference-less equivalent to interference and the ultimate root of bonding. Atoms (or fragments) in molecules should be understood as a statistical mixture of components differing in electron number, spin, etc. The theory of chemical bonding relies on arbitrary references. Here the authors report a fundamental study on the chemical bond showing that considering the binding fragments as objects in real space enables to eliminate inherent biases.
Collapse
Affiliation(s)
- Ángel Martín Pendás
- Departamento de Química Física y Analítica, Universidad de Oviedo, 33006, Oviedo, Spain.
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Universidad de Oviedo, 33006, Oviedo, Spain
| |
Collapse
|
40
|
Kalita AJ, Sarmah K, Borah RR, Yashmin F, Mazumder LJ, Purkayastha SK, Das K, Dutta T, Guha AK. Missing Recipe in the Na‐B Bond in NaBH
3
–
Cluster. ChemistrySelect 2022. [DOI: 10.1002/slct.202201536] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Amlan J. Kalita
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| | - Kangkan Sarmah
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| | - Ritam R. Borah
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| | - Farnaz Yashmin
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| | - Lakhya J. Mazumder
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| | | | - Kanwaki Das
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| | - Trisha Dutta
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| | - Ankur K. Guha
- Advanced Computational Chemistry Centre Cotton University Panbazar Guwahati Assam India- 781001
| |
Collapse
|
41
|
Jara-Cortés J, Matta CF, Hernández-Trujillo J. A fast approximate extension of the interacting quantum atoms energy decomposition to excited states. J Comput Chem 2022; 43:1068-1078. [PMID: 35470908 DOI: 10.1002/jcc.26863] [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: 12/21/2021] [Revised: 02/25/2022] [Accepted: 03/23/2022] [Indexed: 11/11/2022]
Abstract
An approach is developed for the fast calculation of the interacting quantum atoms energy decomposition (IQA) from the information contained in the first order reduced density matrix only. The proposed methodology utilizes an approximate exchange-correlation density from Density Matrix Functional Theory without the need to evaluate the correlation-exchange contribution directly. Instead, weight factors are estimated to decompose the exact Vxc into atomic and pairwise contributions. In this way, the sum of the IQA contributions recovers the energy obtained from the electronic structure calculation. This method can, hence, be applied to obtain atomic contributions in excited states on the same footing as in their ground states using any method that delivers the reduced first-order density matrix. In this way, one can locate chromophores from first principles quantum chemical calculations. Test calculations on the ground and excited states of a set of small molecules indicate that the scaled atomic contributions reproduce vertical electronic transition energies calculated exactly. This approach may be useful to extend the applicability of the IQA approach in the study of large photochemical systems especially when the calculations of the second order reduced density matrices is prohibitive or not possible.
Collapse
Affiliation(s)
- Jesús Jara-Cortés
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic, Mexico
| | - Chérif F Matta
- Department of Chemistry and Physics, Mount Saint Vincent University, Halifax, Nova Scotia, Canada
| | - Jesús Hernández-Trujillo
- Departamento de Física y Química Teórica, Facultad de Química, UNAM. Circuito Escolar, Ciudad Universitaria, Mexico City, Mexico
| |
Collapse
|
42
|
Ochoa-Resendiz D, Gamboa-Suárez A, Hernández-Lamoneda R. Halogen bonding and rotational disorder in chlorine clathrate hydrate cages. J Chem Phys 2022; 156:124302. [DOI: 10.1063/5.0082604] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a detailed theoretical characterization of the structure and interactions in dichlorine clathrate hydrate cages. In the case of the dodecahedral cage there is clear evidence of the presence of halogen bonding whereas in the tetrakaidecahedral the expected signatures are there but in a weaker form. Comparison is made with the available structural data from X-ray experiments, where the rotational motion of the dichlorine has been taken into account through Monte Carlo simulations illustrating delocalization effects associated with sampling multiple minima, specifically for the larger cage. Finally, the intermolecular potentials have been calculated with local correlation methods and energy decomposition analysis has been applied to shed light on the nature of the interactions.
Collapse
|
43
|
Müller F, Wickemeyer L, Schwabedissen J, Ertl M, Neumann B, Stammler HG, Monkowius U, Mitzel NW. Synthesis, structural and photophysical properties of dimethylphosphino(perfluoro-)phenylene-based gold(I) dimers. Dalton Trans 2022; 51:1955-1967. [PMID: 35023528 DOI: 10.1039/d1dt03658k] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Starting with 1,2-dibromobenzene and 1,2,3,4-tetrafluorobenzene, dimethyl(2-(trimethylstannyl)phenyl)phosphane, Me2P(o-C6H4)SnMe3, and dimethyl-[2,3,4,5-tetrafluoro-6-(trimethylstannyl)phenyl]phosphane, Me2P(o-C6F4)SnMe3, were synthesized and used in tin-gold exchange reactions to prepare two gold(I) dimers, bis[(2-dimethylphosphino)phenyl]di-gold(I), [Au2(μ-2-C6H4PMe2)2], and bis[(2-dimethylphosphino)-3,4,5,6-tetrafluorophenyl]di-gold(I), [Au2(μ-2-C6F4PMe2)2], respectively. Both tin precursor molecules, as well as the gold(I) complexes, were characterized by multinuclear NMR spectroscopy, CHN analysis and X-ray diffraction experiments. Both gold(I) dimers were further investigated by the means of computational as well as photophysical methods.
Collapse
Affiliation(s)
- Felix Müller
- Lehrstuhl für Anorganische Chemie und Strukturchemie; Centrum für Molekulare Materialien CM2, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany.
| | - Lucas Wickemeyer
- Lehrstuhl für Anorganische Chemie und Strukturchemie; Centrum für Molekulare Materialien CM2, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany.
| | - Jan Schwabedissen
- Lehrstuhl für Anorganische Chemie und Strukturchemie; Centrum für Molekulare Materialien CM2, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany.
| | - Martin Ertl
- School of Education, Chemistry, Johannes Kepler University Linz, Altenbergerstr. 69, A-4040 Linz, Austria
| | - Beate Neumann
- Lehrstuhl für Anorganische Chemie und Strukturchemie; Centrum für Molekulare Materialien CM2, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany.
| | - Hans-Georg Stammler
- Lehrstuhl für Anorganische Chemie und Strukturchemie; Centrum für Molekulare Materialien CM2, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany.
| | - Uwe Monkowius
- School of Education, Chemistry, Johannes Kepler University Linz, Altenbergerstr. 69, A-4040 Linz, Austria
| | - Norbert W Mitzel
- Lehrstuhl für Anorganische Chemie und Strukturchemie; Centrum für Molekulare Materialien CM2, Fakultät für Chemie, Universität Bielefeld, Universitätsstraße 25, D-33615 Bielefeld, Germany.
| |
Collapse
|
44
|
Jara-Cortés J, Leal-Sánchez E, Francisco E, Pérez-Pimienta JA, Martín Pendás Á, Hernández-Trujillo J. Implementation of the interacting quantum atom energy decomposition using the CASPT2 method. Phys Chem Chem Phys 2021; 23:27508-27519. [PMID: 34874377 DOI: 10.1039/d1cp02837e] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
We present an implementation of the interacting quantum atom (IQA) energy decomposition scheme using the complete active space second-order perturbation theory (CASPT2). This combination yields a real-space interpretation tool with a proper account of the static and dynamic correlation that is particularly relevant for the description of processes in electronic excited states. The IQA/CASPT2 approach allows determination of the energy redistribution that takes place along a photophysical/photochemical deactivation path in terms of self- and interatomic contributions. The applicability of the method is illustrated by the description of representative processes spanning different bonding regimes: noble gas excimer and exciplex formation, the reaction of ozone with a chlorine atom, and the photodissociations of formaldehyde and cyclobutane. These examples show the versatility of using CASPT2 with the significant information provided by the IQA partition to describe chemical processes with a large multiconfigurational character.
Collapse
Affiliation(s)
- Jesús Jara-Cortés
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Edith Leal-Sánchez
- Departamento de Física y Química Teórica, Facultad de Química, UNAM, México City 04510, Mexico
| | - Evelio Francisco
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | - José A Pérez-Pimienta
- Unidad Académica de Ciencias Básicas e Ingenierías, Universidad Autónoma de Nayarit, Tepic 63155, Mexico.
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Faculta de Química, Universidad de Oviedo, Oviedo 33006, Spain
| | | |
Collapse
|
45
|
Karachi SS, Eskandari K. Boron Triel Bonds: A Quantum Chemical Topology Perspective. ChemistrySelect 2021. [DOI: 10.1002/slct.202103002] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Sara S. Karachi
- Department of Chemistry Isfahan University of Technology Isfahan 84156-83111 Iran
| | - Kiamars Eskandari
- Department of Chemistry Isfahan University of Technology Isfahan 84156-83111 Iran
| |
Collapse
|
46
|
Symons BCB, Bane MK, Popelier PLA. DL_FFLUX: A Parallel, Quantum Chemical Topology Force Field. J Chem Theory Comput 2021; 17:7043-7055. [PMID: 34617748 PMCID: PMC8582247 DOI: 10.1021/acs.jctc.1c00595] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
![]()
DL_FFLUX is a force
field based on quantum chemical topology that
can perform molecular dynamics for flexible molecules endowed with
polarizable atomic multipole moments (up to hexadecapole). Using the
machine learning method kriging (aka Gaussian process regression),
DL_FFLUX has access to atomic properties (energy, charge, dipole moment,
etc.) with quantum mechanical accuracy. Newly optimized and parallelized
using domain decomposition Message Passing Interface (MPI), DL_FFLUX
is now able to deliver this rigorous methodology at scale while still
in reasonable time frames. DL_FFLUX is delivered as an add-on to the
widely distributed molecular dynamics code DL_POLY 4.08. For the systems
studied here (103–105 atoms), DL_FFLUX
is shown to add minimal computational cost to the standard DL_POLY
package. In fact, the optimization of the electrostatics in DL_FFLUX
means that, when high-rank multipole moments are enabled, DL_FFLUX
is up to 1.25× faster than standard DL_POLY. The parallel DL_FFLUX
preserves the quality of the scaling of MPI implementation in standard
DL_POLY. For the first time, it is feasible to use the full capability
of DL_FFLUX to study systems that are large enough to be of real-world
interest. For example, a fully flexible, high-rank polarized (up to
and including quadrupole moments) 1 ns simulation of a system of 10 125
atoms (3375 water molecules) takes 30 h (wall time) on 18 cores.
Collapse
Affiliation(s)
- Benjamin C B Symons
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain.,Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| | - Michael K Bane
- High End Compute LTD, 23 Welby Street, Manchester M13 0EL, Great Britainhttps://highendcompute.co.uk.,Department of Computing and Mathematics, Manchester Metropolitan University, Manchester M15 6BH, Great Britain
| | - Paul L A Popelier
- Manchester Institute of Biotechnology (MIB), 131 Princess Street, Manchester M1 7DN, Great Britain.,Department of Chemistry, University of Manchester, Oxford Road, Manchester M13 9PL, Great Britain
| |
Collapse
|
47
|
Cornaton Y, Djukic JP. Noncovalent Interactions in Organometallic Chemistry: From Cohesion to Reactivity, a New Chapter. Acc Chem Res 2021; 54:3828-3840. [PMID: 34617728 DOI: 10.1021/acs.accounts.1c00393] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Noncovalent interactions (NCIs) have long interested a vast community of chemists who investigated their "canonical categories" derived from descriptive crystallography, e.g., H-bonds, π-π interactions, halogen/chalcogen/tetrel bonds, cation-π and C-H-π interactions, metallophilic interactions in the broad sense, etc. Recent developments in theoretical chemistry have enabled the treatment of noncovalent interactions under new auspices: dispersion-force-inclusive density functionals have emerged, which are reliable for modeling small to large molecular systems. It is possible to perform the full analysis of the contributions of London, Debye, and Keesom forces, i.e., the main components of van der Waals forces, by the DFT-D and ab initio methods at a reasonable computational cost. Our research has been focusing for now 15 years on the role of NCIs in the cohesion of organometallic complexes. NCIs are not only effective in Werner's secondary coordination sphere but also in the metal's primary one. The stabilization of electron-unsaturated transition metal complexes by hemichelation, metal-metal donor-acceptor complexes, and self-aggregation of cationic Rh(I) chromophores have indeed outlined the significance of the London dispersion force as an attractive force operating throughout the whole molecule or molecular assembly. The recent outburst of interest in C-H bond functionalization led us to address the broader question of reaction and catalyst engineering: although one can now satisfactorily analyze bonding and molecular cohesion in transition-metal-based organometallic systems, can modern theoretical methods guide reactivity exploration and the engineering of novel catalytic systems? We addressed this question by investigating the ambiphilic metal-ligand activation/concerted metalation-deprotonation mechanism involved in transition-metal-catalyzed directed C-H bond functionalization. This endeavor was initiated having in scope the construction of a rationale for the transposition of 4-5d metal chemistry to earth-abundant 3d metals. In this base-assisted mechanism of C-H bond metalation, agostic interactions are necessary but not sufficient because C-H bond breaking actually relies on the attractive NCI coding of a proton-transfer step and the minimization of metal-H repulsion. This Account introduces the recent shift of our research toward the construction of an NCI-inclusive paradigm of chemical reactivity engineering based on experimental efforts propped up by state-of-the-art theoretical tools.
Collapse
Affiliation(s)
- Yann Cornaton
- Laboratoire de Chimie et Systémique Organométallique, Institut de Chimie de Strasbourg UMR7177, CNRS - Université de Strasbourg, 4 rue Blaise Pascal, F-67070 Strasbourg, France
| | - Jean-Pierre Djukic
- Laboratoire de Chimie et Systémique Organométallique, Institut de Chimie de Strasbourg UMR7177, CNRS - Université de Strasbourg, 4 rue Blaise Pascal, F-67070 Strasbourg, France
| |
Collapse
|
48
|
Barrales-Martínez C, Gutiérrez-Oliva S, Toro-Labbé A, Pendás ÁM. Interacting Quantum Atoms Analysis of the Reaction Force: A Tool to Analyze Driving and Retarding Forces in Chemical Reactions. Chemphyschem 2021; 22:1976-1988. [PMID: 34293240 DOI: 10.1002/cphc.202100428] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2021] [Revised: 07/21/2021] [Indexed: 12/22/2022]
Abstract
The analysis of the reaction force and its topology has provided a wide range of fruitful concepts in the theory of chemical reactivity over the years, allowing to identify chemically relevant regions along a reaction profile. The reaction force (RF), a projection of the Hellmann-Feynman forces acting on the nuclei of a molecular system onto a suitable reaction coordinate, is partitioned using the interacting quantum atoms approach (IQA). The exact IQA molecular energy decomposition is now shown to open a unique window to identify and quantify the chemical entities that drive or retard a chemical reaction. The RF/IQA coupling offers an extraordinarily detailed view of the type and number of elementary processes that take reactants into products, as tested on two sets of simple reactions.
Collapse
Affiliation(s)
- César Barrales-Martínez
- Departamento de Química Orgánica y Fisicoquímica, Facultad de Ciencias Químicas y Farmacéuticas, Universidad de Chile, Sergio Livingstone 1007, Independencia, Santiago, Chile
| | - Soledad Gutiérrez-Oliva
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago, Chile
| | - Alejandro Toro-Labbé
- Laboratorio de Química Teórica Computacional (QTC), Departamento de Química-Física, Facultad de Química y de Farmacia, Pontificia Universidad Católica de Chile, Av. Vicuña Mackenna 4860, Macul, Santiago, Chile
| | - Ángel Martín Pendás
- Departamento de Química Física y Analítica, Facultad de Química, Universidad de Oviedo, 33006, Oviedo, Spain
| |
Collapse
|
49
|
Parida R, Inostroza Rivera R, Sinha S, Giri S. Can Superhalogen Ligand Make More Reactive Frustrated Lewis Pairs? ChemistrySelect 2021. [DOI: 10.1002/slct.202102189] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Affiliation(s)
- Rakesh Parida
- School of Applied Science and Humanities Haldia Institute of Technology Haldia 721657 India
- Department of Chemistry National Institute of Technology Rourkela Odisha 769008 India
| | - Ricardo Inostroza Rivera
- Departmento de Quimica Organica Facultad de Química y de Farmacia Pontificia Universidad Católica de Chile Macul 7820436 Santiago Chile
| | - Swapan Sinha
- School of Applied Science and Humanities Haldia Institute of Technology Haldia 721657 India
| | - Santanab Giri
- School of Applied Science and Humanities Haldia Institute of Technology Haldia 721657 India
| |
Collapse
|
50
|
Pino‐Rios R, Inostroza D, Tiznado W. Neither too Classic nor too Exotic: One‐Electron Na⋅B Bond in NaBH
3
−
Cluster. Angew Chem Int Ed Engl 2021. [DOI: 10.1002/ange.202101403] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/11/2022]
Affiliation(s)
- Ricardo Pino‐Rios
- Laboratorio de Química teórica Facultad de Química y Biología Universidad de Santiago de Chile (USACH) Av. Libertador Bernardo O'Higgins 3363 Santiago, Estación Central, Región Metropolitana Chile
| | - Diego Inostroza
- Universidad Andres Bello Programa de Doctorado en Fisicoquímica Molecular Facultad de Ciencias Exactas Santiago Chile
- Computational and Theoretical Chemistry Group Departamento de Ciencias Químicas Facultad de Ciencias Exactas Universidad Andres Bello República 498 Santiago Chile
| | - William Tiznado
- Universidad Andres Bello Programa de Doctorado en Fisicoquímica Molecular Facultad de Ciencias Exactas Santiago Chile
- Computational and Theoretical Chemistry Group Departamento de Ciencias Químicas Facultad de Ciencias Exactas Universidad Andres Bello República 498 Santiago Chile
| |
Collapse
|