1
|
Skogh M, Dobrautz W, Lolur P, Warren C, Biznárová J, Osman A, Tancredi G, Bylander J, Rahm M. The electron density: a fidelity witness for quantum computation. Chem Sci 2024; 15:2257-2265. [PMID: 38332826 PMCID: PMC10848700 DOI: 10.1039/d3sc05269a] [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: 10/05/2023] [Accepted: 12/19/2023] [Indexed: 02/10/2024] Open
Abstract
There is currently no combination of quantum hardware and algorithms that can provide an advantage over conventional calculations of molecules or materials. However, if or when such a point is reached, new strategies will be needed to verify predictions made using quantum devices. We propose that the electron density, obtained through experimental or computational means, can serve as a robust benchmark for validating the accuracy of quantum computation of chemistry. An initial exploration into topological features of electron densities, facilitated by quantum computation, is presented here as a proof of concept. Additionally, we examine the effects of constraining and symmetrizing measured one-particle reduced density matrices on noise-driven errors in the electron density distribution. We emphasize the potential benefits and future need for high-quality electron densities derived from diffraction experiments for validating classically intractable quantum computations of materials.
Collapse
Affiliation(s)
- Mårten Skogh
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology Gothenburg Sweden
- Data Science & Modelling, Pharmaceutical Science, R&D, AstraZeneca Gothenburg Sweden
| | - Werner Dobrautz
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology Gothenburg Sweden
| | - Phalgun Lolur
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology Gothenburg Sweden
| | - Christopher Warren
- Department of Microtechnology and Nanoscience MC2, Chalmers University of Technology Gothenburg Sweden
| | - Janka Biznárová
- Department of Microtechnology and Nanoscience MC2, Chalmers University of Technology Gothenburg Sweden
| | - Amr Osman
- Department of Microtechnology and Nanoscience MC2, Chalmers University of Technology Gothenburg Sweden
| | - Giovanna Tancredi
- Department of Microtechnology and Nanoscience MC2, Chalmers University of Technology Gothenburg Sweden
| | - Jonas Bylander
- Department of Microtechnology and Nanoscience MC2, Chalmers University of Technology Gothenburg Sweden
| | - Martin Rahm
- Department of Chemistry and Chemical Engineering, Chalmers University of Technology Gothenburg Sweden
| |
Collapse
|
2
|
Stanton SA, Du JJ, Lai F, Stanton G, Hawkins BA, Ong JA, Groundwater PW, Platts JA, Hibbs DE. Understanding Hygroscopicity of Theophylline via a Novel Cocrystal Polymorph: A Charge Density Study. J Phys Chem A 2021; 125:9736-9756. [PMID: 34731566 DOI: 10.1021/acs.jpca.0c09536] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The charge density distribution in a novel cocrystal (1) complex of 1,3-dimethylxanthine (theophylline) and propanedioic acid (malonic acid) has been determined. The molecules crystallize in the triclinic, centrosymmetric space group P1̅, with four independent molecules (Z = 4) in the asymmetric unit (two molecules each of theophylline and malonic acid). Theophylline has a notably high hygroscopic nature, and numerous cocrystals have shown a significant improvement in stability to humidity. A charge density study of the novel polymorph has identified interesting theoretical results correlating the stability enhancement of theophylline via cocrystallization. Topological analysis of the electron density highlighted key differences (up to 17.8) in Laplacian (∇2ρ) between the experimental (EXP) and single-point (SP) models, mainly around intermolecular-bonded carbonyls. Further investigation via molecular electrostatic potential maps reaffirmed that the charge redistribution enhanced intramolecular hydrogen bonding, predominantly for N(2') and N(2) (61.2 and 61.8 kJ mol-1, respectively). An overall weaker lattice energy of the triclinic form (-126.1 kJ mol-1) compared to that of the monoclinic form (-133.8 kJ mol-1) suggests a lower energy threshold to overcome to initiate dissociation. Future work via physical testing of the novel cocrystal in both dissolution and solubility will further solidify the correlation between theoretical and experimental results.
Collapse
Affiliation(s)
- Stephen A Stanton
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| | - Jonathan J Du
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| | - Felcia Lai
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| | - Gyte Stanton
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| | - Bryson A Hawkins
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| | - Jennifer A Ong
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| | - Paul W Groundwater
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| | - James A Platts
- School of Chemistry, Cardiff University, Cardiff CF10 3AT, U.K
| | - David E Hibbs
- Sydney School of Pharmacy, Faculty of Medicine and Health, The University of Sydney, Sydney, NSW 2006, Australia
| |
Collapse
|
3
|
Evarestov RA, Kuzmin A. Topological analysis of chemical bonding in the layered FePSe 3 upon pressure-induced phase transitions. J Comput Chem 2020; 41:2610-2623. [PMID: 32905632 DOI: 10.1002/jcc.26416] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2020] [Revised: 08/13/2020] [Accepted: 08/20/2020] [Indexed: 12/31/2022]
Abstract
Two pressure-induced phase transitions have been theoretically studied in the layered iron phosphorus triselenide (FePSe3 ). Topological analysis of chemical bonding in FePSe3 has been performed based on the results of first-principles calculations within the periodic linear combination of atomic orbitals (LCAO) method with hybrid Hartree-Fock-DFT B3LYP functional. The first transition at about 6 GPa is accompanied by the symmetry change from R 3 ¯ to C2/m, whereas the semiconductor-to-metal transition (SMT) occurs at about 13 GPa leading to the symmetry change from C2/m to P 3 ¯ 1 m . We found that the collapse of the band gap at about 13 GPa occurs due to changes in the electronic structure of FePSe3 induced by relative displacements of phosphorus or selenium atoms along the c-axis direction under pressure. The results of the topological analysis of the electron density and its Laplacian demonstrate that the pressure changes not only the interatomic distances but also the bond nature between the intralayer and interlayer phosphorus atoms. The interlayer P-P interactions are absent in two non-metallic FePSe3 phases while after SMT the intralayer P-P interactions weaken and the interlayer P-P interactions appear.
Collapse
Affiliation(s)
- Robert A Evarestov
- Department of Quantum Chemistry, Saint Petersburg State University, St. Petersburg, Russian Federation
| | - Alexei Kuzmin
- Institute of Solid State Physics, University of Latvia, Riga, Latvia
| |
Collapse
|
4
|
Korlyukov AA, Nelyubina YV. Quantum chemical methods in charge density studies from X-ray diffraction data. RUSSIAN CHEMICAL REVIEWS 2019. [DOI: 10.1070/rcr4866] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/08/2022]
|
5
|
Shi MW, Thomas SP, Hathwar VR, Edwards AJ, Piltz RO, Jayatilaka D, Koutsantonis GA, Overgaard J, Nishibori E, Iversen BB, Spackman MA. Measurement of Electric Fields Experienced by Urea Guest Molecules in the 18-Crown-6/Urea (1:5) Host-Guest Complex: An Experimental Reference Point for Electric-Field-Assisted Catalysis. J Am Chem Soc 2019; 141:3965-3976. [PMID: 30761898 DOI: 10.1021/jacs.8b12927] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
High-resolution synchrotron and neutron single-crystal diffraction data of 18-crown-6/(pentakis)urea measured at 30 K are combined, with the aim of better appreciating the electrostatics associated with intermolecular interactions in condensed matter. With two 18-crown-6 molecules and five different urea molecules in the crystal, this represents the most ambitious combined X-ray/synchrotron and neutron experimental charge density analysis to date on a cocrystal or host-guest system incorporating such a large number of unique molecules. The dipole moments of the five urea guest molecules in the crystal are enhanced considerably compared to values determined for isolated molecules, and 2D maps of the electrostatic potential and electric field show clearly how the urea molecules are oriented with dipole moments aligned along the electric field exerted by their molecular neighbors. Experimental electric fields in the range of 10-19 GV m-1, obtained for the five different urea environments, corroborate independent measurements of electric fields in the active sites of enzymes and provide an important experimental reference point for recent discussions focused on electric-field-assisted catalysis.
Collapse
Affiliation(s)
- Ming W Shi
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - Sajesh P Thomas
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia.,Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Venkatesha R Hathwar
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark.,Division of Physics, Faculty of Pure and Applied Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8571 , Japan
| | - Alison J Edwards
- Australian Nuclear Science and Technology Organization , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , New South Wales 2234 , Australia
| | - Ross O Piltz
- Australian Nuclear Science and Technology Organization , Australian Centre for Neutron Scattering , New Illawarra Road , Lucas Heights , New South Wales 2234 , Australia
| | - Dylan Jayatilaka
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - George A Koutsantonis
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| | - Jacob Overgaard
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Eiji Nishibori
- Division of Physics, Faculty of Pure and Applied Sciences , University of Tsukuba , 1-1-1 Tennodai , Tsukuba , Ibaraki 305-8571 , Japan
| | - Bo B Iversen
- Center for Materials Crystallography and Department of Chemistry , Aarhus University , Langelandsgade 140 , DK-8000 Aarhus C , Denmark
| | - Mark A Spackman
- School of Molecular Sciences , University of Western Australia , 35 Stirling Highway , Crawley , WA 6009 , Australia
| |
Collapse
|
6
|
Nikolova V, Cheshmedzhieva D, Ilieva S, Galabov B. Atomic Charges in Describing Properties of Aromatic Molecules. J Org Chem 2019; 84:1908-1915. [PMID: 30620875 DOI: 10.1021/acs.joc.8b02908] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The performance of four frequently employed population analysis methods is assessed by comparisons with experimentally derived properties of monosubstituted benzene derivatives. The analysis is based on the expected dependence between site reactivities and electron densities at the respective ring carbon atoms. The correspondence between charges obtained from Mulliken, NPA, Hirshfeld, and QTAIM approaches and the σ0m and σ0p aromatic substituent constants is examined. The series of molecules investigated includes benzene and 18 monosubstituted derivatives. The atomic charges are derived using the B3LYP, ωB97X-D density functional, and MP2 MO methods combined with the 6-311++G(3df,2pd) basis set. A quantitative correspondence between Hirshfeld charges and σ0 constants is established. Application of Møller-Plesset second-order perturbation theory (MP2) wave functions appears to be essential in obtaining a more realistic electron density distribution. NPA and QTAIM charges provide in most cases a satisfactory description of the substituent effects. The net transfer of charges between substituents and the aromatic ring is assessed.
Collapse
Affiliation(s)
- Valia Nikolova
- Department of Chemistry and Pharmacy , University of Sofia , Sofia 1164 , Bulgaria
| | - Diana Cheshmedzhieva
- Department of Chemistry and Pharmacy , University of Sofia , Sofia 1164 , Bulgaria
| | - Sonia Ilieva
- Department of Chemistry and Pharmacy , University of Sofia , Sofia 1164 , Bulgaria
| | - Boris Galabov
- Department of Chemistry and Pharmacy , University of Sofia , Sofia 1164 , Bulgaria
| |
Collapse
|
7
|
Varadwaj PR, Varadwaj A, Marques HM, Yamashita K. Significance of hydrogen bonding and other noncovalent interactions in determining octahedral tilting in the CH 3NH 3PbI 3 hybrid organic-inorganic halide perovskite solar cell semiconductor. Sci Rep 2019; 9:50. [PMID: 30631082 PMCID: PMC6328624 DOI: 10.1038/s41598-018-36218-1] [Citation(s) in RCA: 53] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2018] [Accepted: 11/16/2018] [Indexed: 12/19/2022] Open
Abstract
The CH3NH3PbI3 (methylammonium lead triiodide) perovskite semiconductor system has been viewed as a blockbuster research material during the last five years. Because of its complicated architecture, several of its technological, physical and geometrical issues have been examined many times. Yet this has not assisted in overcoming a number of problems in the field nor in enabling the material to be marketed. For instance, these studies have not clarified the nature and type of hydrogen bonding and other noncovalent interactions involved; the origin of hysteresis; the actual role of the methylammonium cation; the nature of polarity associated with the tetragonal geometry; the unusual origin of various frontier orbital contributions to the conduction band minimum; the underlying phenomena of spin-orbit coupling that causes significant bandgap reduction; and the nature of direct-to-indirect bandgap transition features. Arising from many recent reports, it is now a common belief that the I···H–N interaction formed between the inorganic framework and the ammonium group of CH3NH3+ is the only hydrogen bonded interaction responsible for all temperature-dependent geometrical polymorphs of the system, including the most stable one that persists at low-temperatures, and the significance of all other noncovalent interactions has been overlooked. This study focussed only on the low temperature orthorhombic polymorph of CH3NH3PbI3 and CD3ND3PbI3, where D refers deuterium. Together with QTAIM, DORI and RDG based charge density analyses, the results of density functional theory calculations with PBE with and without van der Waals corrections demonstrate that the prevailing view of hydrogen bonding in CH3NH3PbI3 is misleading as it does not alone determine the a−b+a− tilting pattern of the PbI64− octahedra. This study suggests that it is not only the I···H/D–N, but also the I···H/D–C hydrogen/deuterium bonding and other noncovalent interactions (viz. tetrel-, pnictogen- and lump-hole bonding interactions) that are ubiquitous in the orthorhombic CH3NH3PbI3/CD3ND3PbI3 perovskite geometry. Their interplay determines the overall geometry of the polymorph, and are therefore responsible in part for the emergence of the functional optical properties of this material. This study also suggests that these interactions should not be regarded as the sole determinants of octahedral tilting since lattice dynamics is known to play a critical role as well, a common feature in many inorganic perovskites both in the presence and the absence of the encaged cation, as in CsPbI3/WO3 perovskites, for example.
Collapse
Affiliation(s)
- Pradeep R Varadwaj
- Department of Chemical System Engineering, School of Engineering, The University of Tokyo 7-3-1, Hongo, Bunkyo-ku, 113-8656, Japan. .,CREST-JST, 7 Gobancho, Chiyoda-ku, Tokyo, 102-0076, Japan. .,The National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, 305-8560, Japan.
| | - Arpita Varadwaj
- Department of Chemical System Engineering, School of Engineering, The University of Tokyo 7-3-1, Hongo, Bunkyo-ku, 113-8656, Japan.,CREST-JST, 7 Gobancho, Chiyoda-ku, Tokyo, 102-0076, Japan.,The National Institute of Advanced Industrial Science and Technology (AIST), Tsukuba, Ibaraki, 305-8560, Japan
| | - Helder M Marques
- Molecular Sciences Institute, School of Chemistry, University of the Witwatersrand, Johannesburg, 2050, South Africa
| | - Koichi Yamashita
- Department of Chemical System Engineering, School of Engineering, The University of Tokyo 7-3-1, Hongo, Bunkyo-ku, 113-8656, Japan.,CREST-JST, 7 Gobancho, Chiyoda-ku, Tokyo, 102-0076, Japan
| |
Collapse
|
8
|
Lopes TO, Machado DFS, Risko C, Brédas JL, de Oliveira HCB. Bond Ellipticity Alternation: An Accurate Descriptor of the Nonlinear Optical Properties of π-Conjugated Chromophores. J Phys Chem Lett 2018; 9:1377-1383. [PMID: 29498530 DOI: 10.1021/acs.jpclett.8b00478] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
Well-defined structure-property relationships offer a conceptual basis to afford a priori design principles to develop novel π-conjugated molecular and polymer materials for nonlinear optical (NLO) applications. Here, we introduce the bond ellipticity alternation (BEA) as a robust parameter to assess the NLO characteristics of organic chromophores and illustrate its effectiveness in the case of streptocyanines. BEA is based on the symmetry of the electron density, a physical observable that can be determined from experimental X-ray electron densities or from quantum-chemical calculations. Through comparisons to the well-established bond-length alternation and π-bond order alternation parameters, we demonstrate the generality of BEA to foreshadow NLO characteristics and underline that, in the case of large electric fields, BEA is a more reliable descriptor. Hence, this study introduces BEA as a prominent descriptor of organic chromophores of interest for NLO applications.
Collapse
Affiliation(s)
- Thiago O Lopes
- Laboratório de Estrutura Eletrônica e Dinâmica Molecular (LEEDMOL), Institute of Chemistry , University of Brasília , Campus Darcy Ribeiro , Brasília , Brazil
| | - Daniel F Scalabrini Machado
- Laboratório de Estrutura Eletrônica e Dinâmica Molecular (LEEDMOL), Institute of Chemistry , University of Brasília , Campus Darcy Ribeiro , Brasília , Brazil
| | - Chad Risko
- Department of Chemistry and Center for Applied Energy Research , University of Kentucky , Lexington , Kentucky 40506-0055 , United States
| | - Jean-Luc Brédas
- School of Chemistry and Biochemistry and Center for Organic Photonics and Electronics , Georgia Institute of Technology , Atlanta , Georgia 30332-0400 , United States
| | - Heibbe C B de Oliveira
- Laboratório de Estrutura Eletrônica e Dinâmica Molecular (LEEDMOL), Institute of Chemistry , University of Brasília , Campus Darcy Ribeiro , Brasília , Brazil
| |
Collapse
|
9
|
Merli M, Pavese A. Electron-density critical points analysis and catastrophe theory to forecast structure instability in periodic solids. ACTA CRYSTALLOGRAPHICA A-FOUNDATION AND ADVANCES 2018; 74:102-111. [PMID: 29493539 DOI: 10.1107/s2053273317018381] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Accepted: 12/22/2017] [Indexed: 11/10/2022]
Abstract
The critical points analysis of electron density, i.e. ρ(x), from ab initio calculations is used in combination with the catastrophe theory to show a correlation between ρ(x) topology and the appearance of instability that may lead to transformations of crystal structures, as a function of pressure/temperature. In particular, this study focuses on the evolution of coalescing non-degenerate critical points, i.e. such that ∇ρ(xc) = 0 and λ1, λ2, λ3 ≠ 0 [λ being the eigenvalues of the Hessian of ρ(x) at xc], towards degenerate critical points, i.e. ∇ρ(xc) = 0 and at least one λ equal to zero. The catastrophe theory formalism provides a mathematical tool to model ρ(x) in the neighbourhood of xc and allows one to rationalize the occurrence of instability in terms of electron-density topology and Gibbs energy. The phase/state transitions that TiO2 (rutile structure), MgO (periclase structure) and Al2O3 (corundum structure) undergo because of pressure and/or temperature are here discussed. An agreement of 3-5% is observed between the theoretical model and experimental pressure/temperature of transformation.
Collapse
Affiliation(s)
- Marcello Merli
- Earth and Sea Sciences Department (DiSTeM), University of Palermo, Via Archirafi 36, Palermo, 90126, Italy
| | - Alessandro Pavese
- Earth Sciences Department, University of Turin, Via Valperga Caluso 35, Turin, 10100, Italy
| |
Collapse
|