151
|
Affiliation(s)
- Yuhong Liu
- Department of Chemistry, University of the Pacific, Stockton, California, USA
| | - Anthony D. Dutoi
- Department of Chemistry, University of the Pacific, Stockton, California, USA
| |
Collapse
|
152
|
Abstract
Two polymorphs of the drug compound metergoline (C25H29N3O2) were investigated in detail by solid-state NMR measurements. The results have been analysed by an advanced procedure, which uses experimental input together with the results of quantum chemical calculations that were performed for molecular crystals. In this way, it was possible to assign the total of 40 1H–13C correlation pairs in a highly complex system, namely, in the dynamically disordered polymorph with two independent molecules in the unit cell of a large volume of 4234 Å3. For the simpler polymorph, which exhibits only small-amplitude motions and has just one molecule in the unit cell with a volume of 529.0 Å3, the values of the principal elements of the 13C chemical shift tensors were measured. Additionally, for this polymorph, a set of crystal structure predictions were generated, and the {13C, 1H} isotropic and 13C anisotropic chemical shielding data were computed while using the gauge-including projector augmented-wave approach combined with the “revised Perdew-Burke-Ernzerhof“ exchange-correlation functional (GIPAW-RPBE). The experimental and theoretical results were combined in an application of the newly developed strategy to polymorph discrimination. This research thus opens up new routes towards more accurate characterization of the polymorphism of drug formulations.
Collapse
|
153
|
Košata J, Merkl P, Teeratchanan P, Hermann A. Stability of Hydrogen Hydrates from Second-Order Møller-Plesset Perturbation Theory. J Phys Chem Lett 2018; 9:5624-5629. [PMID: 30179509 DOI: 10.1021/acs.jpclett.8b02274] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/08/2023]
Abstract
The formation of gas hydrates and clathrates critically depends on the interaction between the host water network and the guest gas species. Density functional calculations can struggle to quantitatively capture these dispersion-type interactions. Here, we report wave function-based calculations on hydrogen hydrates that combine periodic Hartree-Fock with a localized treatment of electronic correlation. We show that local second-order Møller-Plesset perturbation theory (LMP2) reproduces the stability of the different filled-ice-like hydrates in excellent agreement with experimental data. In contrast to various dispersion-corrected density functional theory implementations, LMP2 correctly identifies the pressures needed to stabilize the C0, C1, and C2 hydrates and does not find a spurious region of stability for an ice-Ih-based dihydrate. Our results suggest that LMP2 or similar approaches can provide quantitative insights into the mechanisms of formation and eventual decomposition of molecular host-guest compounds.
Collapse
Affiliation(s)
- Jan Košata
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| | - Padryk Merkl
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| | - Pattanasak Teeratchanan
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| | - Andreas Hermann
- Centre for Science at Extreme Conditions and SUPA, School of Physics and Astronomy , The University of Edinburgh , Edinburgh EH9 3FD , United Kingdom
| |
Collapse
|
154
|
Price SL. Control and prediction of the organic solid state: a challenge to theory and experiment †. Proc Math Phys Eng Sci 2018; 474:20180351. [PMID: 30333710 PMCID: PMC6189584 DOI: 10.1098/rspa.2018.0351] [Citation(s) in RCA: 31] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2018] [Accepted: 08/15/2018] [Indexed: 11/12/2022] Open
Abstract
The ability of theoretical chemists to quantitatively model the weak forces between organic molecules is being exploited to predict their crystal structures and estimate their physical properties. Evolving crystal structure prediction methods are increasingly being used to aid the design of organic functional materials and provide information about thermodynamically plausible polymorphs of speciality organic materials to aid, for example, pharmaceutical development. However, the increasingly sophisticated experimental studies for detecting the range of organic solid-state behaviours provide many challenges for improving quantitative theories that form the basis for the computer modelling. It is challenging to calculate the relative thermodynamic stability of different organic crystal structures, let alone understand the kinetic effects that determine which polymorphs can be observed and are practically important. However, collaborations between experiment and theory are reaching the stage of devising experiments to target the first crystallization of new polymorphs or create novel organic molecular materials.
Collapse
Affiliation(s)
- Sarah L. Price
- Department of Chemistry, University College London, 20 Gordon St, London WC1H 0AJ, UK
| |
Collapse
|
155
|
Pérez-Jiménez AJ, Brémond E, Adamo C, Sancho-García JC. Communication: Accurate description of interaction energies and three-body effects in weakly bound molecular complexes by PBE-QIDH models. J Chem Phys 2018; 149:041101. [PMID: 30068200 DOI: 10.1063/1.5042153] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We apply a recently developed parameter-free double-hybrid density functional belonging to the quadratic-integrand double-hybrid model to calculate association energies (ΔE) and three-body effects (Δ3E) arising from intermolecular interactions in weakly bound supramolecular complexes (i.e., the dataset 3B-69). The model behaves very accurately for trimer association energies and is found to outperform widely used density functional approximations while approaching the accuracy of more costly ab initio methods for three-body effects. The results are further improved when we add some specific corrections for the remaining dispersion interactions, D3(BJ) or VV10 for two-body effects and Axilrod-Teller-Muto for three-body effects, leading to marginal deviations (less than 1 kcal/mol for ΔE and around 0.03-0.04 kcal/mol for Δ3E) with respect to benchmark results.
Collapse
Affiliation(s)
- A J Pérez-Jiménez
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| | - E Brémond
- Université Paris Diderot, Sorbonne Paris Cité, ITODYS, UMR CNRS 7086, F-75013 Paris, France
| | - C Adamo
- Chimie ParisTech, PSL Research University, CNRS, Institut de Recherche de Chimie Paris, F-75005 Paris, France
| | - J C Sancho-García
- Department of Physical Chemistry, University of Alicante, E-03080 Alicante, Spain
| |
Collapse
|
156
|
Li X, Curtis FS, Rose T, Schober C, Vazquez-Mayagoitia A, Reuter K, Oberhofer H, Marom N. Genarris: Random generation of molecular crystal structures and fast screening with a Harris approximation. J Chem Phys 2018; 148:241701. [PMID: 29960303 DOI: 10.1063/1.5014038] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Abstract
We present Genarris, a Python package that performs configuration space screening for molecular crystals of rigid molecules by random sampling with physical constraints. For fast energy evaluations, Genarris employs a Harris approximation, whereby the total density of a molecular crystal is constructed via superposition of single molecule densities. Dispersion-inclusive density functional theory is then used for the Harris density without performing a self-consistency cycle. Genarris uses machine learning for clustering, based on a relative coordinate descriptor developed specifically for molecular crystals, which is shown to be robust in identifying packing motif similarity. In addition to random structure generation, Genarris offers three workflows based on different sequences of successive clustering and selection steps: the "Rigorous" workflow is an exhaustive exploration of the potential energy landscape, the "Energy" workflow produces a set of low energy structures, and the "Diverse" workflow produces a maximally diverse set of structures. The latter is recommended for generating initial populations for genetic algorithms. Here, the implementation of Genarris is reported and its application is demonstrated for three test cases.
Collapse
Affiliation(s)
- Xiayue Li
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Farren S Curtis
- Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Timothy Rose
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| | - Christoph Schober
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universiät München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | | | - Karsten Reuter
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universiät München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Harald Oberhofer
- Chair for Theoretical Chemistry and Catalysis Research Center, Technische Universiät München, Lichtenbergstr. 4, D-85747 Garching, Germany
| | - Noa Marom
- Department of Materials Science and Engineering, Carnegie Mellon University, Pittsburgh, Pennsylvania 15213, USA
| |
Collapse
|
157
|
Caldeweyher E, Brandenburg JG. Simplified DFT methods for consistent structures and energies of large systems. JOURNAL OF PHYSICS. CONDENSED MATTER : AN INSTITUTE OF PHYSICS JOURNAL 2018; 30:213001. [PMID: 29633964 DOI: 10.1088/1361-648x/aabcfb] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/27/2023]
Abstract
Kohn-Sham density functional theory (DFT) is routinely used for the fast electronic structure computation of large systems and will most likely continue to be the method of choice for the generation of reliable geometries in the foreseeable future. Here, we present a hierarchy of simplified DFT methods designed for consistent structures and non-covalent interactions of large systems with particular focus on molecular crystals. The covered methods are a minimal basis set Hartree-Fock (HF-3c), a small basis set screened exchange hybrid functional (HSE-3c), and a generalized gradient approximated functional evaluated in a medium-sized basis set (B97-3c), all augmented with semi-classical correction potentials. We give an overview on the methods design, a comprehensive evaluation on established benchmark sets for geometries and lattice energies of molecular crystals, and highlight some realistic applications on large organic crystals with several hundreds of atoms in the primitive unit cell.
Collapse
Affiliation(s)
- Eike Caldeweyher
- Mulliken Center for Theoretical Chemistry, Institute for Physical and Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115 Bonn, Germany
| | | |
Collapse
|
158
|
Červinka C, Beran GJO. Ab initio prediction of the polymorph phase diagram for crystalline methanol. Chem Sci 2018; 9:4622-4629. [PMID: 29899955 PMCID: PMC5969506 DOI: 10.1039/c8sc01237g] [Citation(s) in RCA: 41] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Accepted: 04/13/2018] [Indexed: 11/21/2022] Open
Abstract
Organic crystals frequently adopt multiple distinct polymorphs exhibiting different properties. The ability to predict not only what crystal forms might occur, but under what experimental thermodynamic conditions those polymorphs are stable would be immensely valuable to the pharmaceutical industry and others. Starting only from knowledge of the experimental crystal structures, this study successfully predicts the methanol crystal polymorph phase diagram from first-principles quantum chemistry, mapping out the thermodynamic regions of stability for three polymorphs over the range 0-400 K and 0-6 GPa. The agreement between the predicted and experimental phase diagrams corresponds to predicting the relative polymorph free energies to within ∼0.5 kJ mol-1 accuracy, which is achieved by employing fragment-based second-order Møller-Plesset perturbation theory and coupled cluster theory plus a quasi-harmonic treatment of the phonons.
Collapse
Affiliation(s)
- Ctirad Červinka
- Department of Physical Chemistry , University of Chemistry and Technology Prague , Technická 5 , CZ-16628 Prague 6 , Czech Republic .
| | - Gregory J O Beran
- Department of Chemistry , University of California , Riverside , California 92521 , USA .
| |
Collapse
|
159
|
Liu C, Kremer K, Bereau T. Polymorphism of Syndiotactic Polystyrene Crystals from Multiscale Simulations. ADVANCED THEORY AND SIMULATIONS 2018. [DOI: 10.1002/adts.201800024] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/25/2023]
Affiliation(s)
- Chan Liu
- Max Planck Institute for Polymer Research; Ackermannweg 10 55128 Mainz Germany
| | - Kurt Kremer
- Max Planck Institute for Polymer Research; Ackermannweg 10 55128 Mainz Germany
| | - Tristan Bereau
- Max Planck Institute for Polymer Research; Ackermannweg 10 55128 Mainz Germany
| |
Collapse
|
160
|
Ruggiero MT, Zhang W, Bond AD, Mittleman DM, Zeitler JA. Uncovering the Connection Between Low-Frequency Dynamics and Phase Transformation Phenomena in Molecular Solids. PHYSICAL REVIEW LETTERS 2018; 120:196002. [PMID: 29799217 DOI: 10.1103/physrevlett.120.196002] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/02/2017] [Indexed: 06/08/2023]
Abstract
The low-frequency motions of molecules in the condensed phase have been shown to be vital to a large number of physical properties and processes. However, in the case of disordered systems, it is often difficult to elucidate the atomic-level details surrounding these phenomena. In this work, we have performed an extensive experimental and computational study on the molecular solid camphor, which exhibits a rich and complex structure-dynamics relationship, and undergoes an order-disorder transition near ambient conditions. The combination of x-ray diffraction, variable temperature and pressure terahertz time-domain spectroscopy, ab initio molecular dynamics, and periodic density functional theory calculations enables a complete picture of the phase transition to be obtained, inclusive of mechanistic, structural, and thermodynamic phenomena. Additionally, the low-frequency vibrations of a disordered solid are characterized for the first time with atomic-level precision, uncovering a clear link between such motions and the phase transformation. Overall, this combination of methods allows for significant details to be obtained for disordered solids and the associated transformations, providing a framework that can be directly applied for a wide range of similar systems.
Collapse
Affiliation(s)
- Michael T Ruggiero
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
- Present Address: Department of Chemistry, University of Vermont, 82 University Place, Burlington, Vermont 05405, USA
| | - Wei Zhang
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - Andrew D Bond
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Daniel M Mittleman
- School of Engineering, Brown University, Providence, Rhode Island 02912, USA
| | - J Axel Zeitler
- Department of Chemical Engineering and Biotechnology, University of Cambridge, Philippa Fawcett Drive, Cambridge CB3 0AS, United Kingdom
| |
Collapse
|
161
|
Romanova A, Lyssenko K, Ananyev I. Estimations of energy of noncovalent bonding from integrals over interatomic zero-flux surfaces: Correlation trends and beyond. J Comput Chem 2018; 39:1607-1616. [PMID: 29701303 DOI: 10.1002/jcc.25235] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2018] [Revised: 03/28/2018] [Accepted: 04/03/2018] [Indexed: 12/24/2022]
Abstract
Bonding energies of 50 associates composed by neutral molecules (atoms) and bounded by various weak noncovalent interactions are calculated within the DFT framework using the PBE0/aug-cc-pVTZ combination. The electronic virial and electron density values at bond critical points together with their integrals over interatomic surfaces are tested to check their ability to estimate bonding energies. Two correlations schemes dealing with integrals over interatomic surface are suggested to estimate bonding energy of any noncovalent interaction. The physical meaning of explored and several known correlations is discussed. Methods to estimate interatomic surface integrals of electronic virial and electron density are proposed. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Anna Romanova
- A.N. Nesmeyanov Institute of Organoelement Compounds RAS, Vavilova 28, Moscow, 119991, Russia.,D.I. Mendeleev University of Chemical Technology of Russia, Miusskaya sq. 9, Moscow, 125047, Russia
| | - Konstantin Lyssenko
- A.N. Nesmeyanov Institute of Organoelement Compounds RAS, Vavilova 28, Moscow, 119991, Russia
| | - Ivan Ananyev
- A.N. Nesmeyanov Institute of Organoelement Compounds RAS, Vavilova 28, Moscow, 119991, Russia.,D.I. Mendeleev University of Chemical Technology of Russia, Miusskaya sq. 9, Moscow, 125047, Russia
| |
Collapse
|
162
|
Lao KU, Herbert JM. Atomic Orbital Implementation of Extended Symmetry-Adapted Perturbation Theory (XSAPT) and Benchmark Calculations for Large Supramolecular Complexes. J Chem Theory Comput 2018; 14:2955-2978. [DOI: 10.1021/acs.jctc.8b00058] [Citation(s) in RCA: 34] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Ka Un Lao
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - John M. Herbert
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
163
|
Tao XD, Chai WX, Song L, Wei QH, Shi HS, Qin LS. Two luminescent pseudo-polymorphic cuprous complexes with different optical properties: Synthesis, characterization and TD-DFT calculations. Polyhedron 2018. [DOI: 10.1016/j.poly.2017.12.039] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
|
164
|
LeBlanc LM, Otero-de-la-Roza A, Johnson ER. Composite and Low-Cost Approaches for Molecular Crystal Structure Prediction. J Chem Theory Comput 2018; 14:2265-2276. [PMID: 29498837 DOI: 10.1021/acs.jctc.7b01179] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Molecular crystal structure prediction (CSP) requires evaluating differences in lattice energy between candidate crystal structures accurately and efficiently. In this work, we explore and compare several low-cost alternatives to dispersion-corrected density-functional theory (DFT) in the plane-waves/pseudopotential approximation, the most accurate and general approach used for CSP at present. Three types of low-cost methods are considered: DFT with a small basis set of finite-support numerical orbitals (the SIESTA method), dispersion-corrected Gaussian small or minimal-basis-set Hartree-Fock and DFT with additional empirical corrections (HF-3c and PBEh-3c), and self-consistent-charge dispersion-corrected density-functional tight binding (SCC-DFTB3-D3). In addition, we study the performance of composite methods that comprise a geometry optimization using a low-cost approach followed by a single-point calculation using the accurate but comparatively expensive B86bPBE-XDM functional. All methods were tested for their abilities to produce absolute lattice energies, relative lattice energies, and crystal geometries. We show that assessing various methods by their ability to produce absolute lattice energies can be misleading and that relative lattice energies are a much better indicator of performance in CSP. The EE14 set of relative solubilities of homochiral and heterochiral chiral crystals is proposed for relative lattice-energy benchmarking. Our results show that PBE-D2 plus a DZP basis set of numerical orbitals coupled with a final B86bPBE-XDM single-point calculation gives excellent performance at a fraction of the cost of a full B86bPBE-XDM calculation, although the results are sensitive to the particular details of the computational protocol. The B86bPBE-XDM//PBE-D2/DZP method was subsequently tested in a practical CSP application from our recent work on the crystal structure of the enantiopure and racemate forms of 1-aza[6]helicene, a chiral organic semiconductor. Our results show that this multilevel method is able to correctly reproduce the energy ranking of both crystal forms.
Collapse
Affiliation(s)
- Luc M LeBlanc
- Department of Chemistry , Dalhousie University , 6274 Coburg Road , P.O. Box 15000, Halifax , Nova Scotia , Canada B3H 4R2
| | - Alberto Otero-de-la-Roza
- Department of Chemistry , University of British Columbia, Okanagan , 3247 University Way , Kelowna , British Columbia , Canada V1V 1V7
| | - Erin R Johnson
- Department of Chemistry , Dalhousie University , 6274 Coburg Road , P.O. Box 15000, Halifax , Nova Scotia , Canada B3H 4R2
| |
Collapse
|
165
|
Xu P, Guidez EB, Bertoni C, Gordon MS. Perspective:Ab initioforce field methods derived from quantum mechanics. J Chem Phys 2018. [DOI: 10.1063/1.5009551] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Affiliation(s)
- Peng Xu
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Emilie B. Guidez
- Department of Chemistry, University of Colorado Denver, Denver, Colorado 80217, USA
| | - Colleen Bertoni
- Argonne Leadership Computing Facility, Argonne, Illinois 60439, USA
| | - Mark S. Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| |
Collapse
|
166
|
Azuri I, Hirsch A, Reilly AM, Tkatchenko A, Kendler S, Hod O, Kronik L. Terahertz spectroscopy of 2,4,6-trinitrotoluene molecular solids from first principles. Beilstein J Org Chem 2018; 14:381-388. [PMID: 29507643 PMCID: PMC5815306 DOI: 10.3762/bjoc.14.26] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/29/2017] [Accepted: 01/22/2018] [Indexed: 11/25/2022] Open
Abstract
We present a computational analysis of the terahertz spectra of the monoclinic and the orthorhombic polymorphs of 2,4,6-trinitrotoluene. Very good agreement with experimental data is found when using density functional theory that includes Tkatchenko–Scheffler pair-wise dispersion interactions. Furthermore, we show that for these polymorphs the theoretical results are only weakly affected by many-body dispersion contributions. The absence of dispersion interactions, however, causes sizable shifts in vibrational frequencies and directly affects the spatial character of the vibrational modes. Mode assignment allows for a distinction between the contributions of the monoclinic and orthorhombic polymorphs and shows that modes in the range from 0 to ca. 3.3 THz comprise both inter- and intramolecular vibrations, with the former dominating below ca. 1.5 THz. We also find that intramolecular contributions primarily involve the nitro and methyl groups. Finally, we present a prediction for the terahertz spectrum of 1,3,5-trinitrobenzene, showing that a modest chemical change leads to a markedly different terahertz spectrum.
Collapse
Affiliation(s)
- Ido Azuri
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Anna Hirsch
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel
| | - Anthony M Reilly
- School of Chemical Sciences, Dublin City University, Glasnevin, Dublin 9, Ireland
| | - Alexandre Tkatchenko
- Physics and Materials Research Unit, University of Luxembourg, L-1511 Luxembourg
| | - Shai Kendler
- Israel Institute for Biological Research, Ness Ziona 74100, Israel
| | - Oded Hod
- Department of Physical Chemistry, School of Chemistry, The Raymond and Beverly Sackler Faculty of Exact Sciences and The Sackler Center for Computational Molecular and Materials Science, Tel Aviv University, Tel Aviv 6997801, Israel
| | - Leeor Kronik
- Department of Materials and Interfaces, Weizmann Institute of Science, Rehovoth 76100, Israel
| |
Collapse
|
167
|
Curtis F, Li X, Rose T, Vázquez-Mayagoitia Á, Bhattacharya S, Ghiringhelli LM, Marom N. GAtor: A First-Principles Genetic Algorithm for Molecular Crystal Structure Prediction. J Chem Theory Comput 2018; 14:2246-2264. [PMID: 29481740 DOI: 10.1021/acs.jctc.7b01152] [Citation(s) in RCA: 48] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We present the implementation of GAtor, a massively parallel, first-principles genetic algorithm (GA) for molecular crystal structure prediction. GAtor is written in Python and currently interfaces with the FHI-aims code to perform local optimizations and energy evaluations using dispersion-inclusive density functional theory (DFT). GAtor offers a variety of fitness evaluation, selection, crossover, and mutation schemes. Breeding operators designed specifically for molecular crystals provide a balance between exploration and exploitation. Evolutionary niching is implemented in GAtor by using machine learning to cluster the dynamically updated population by structural similarity and then employing a cluster-based fitness function. Evolutionary niching promotes uniform sampling of the potential energy surface by evolving several subpopulations, which helps overcome initial pool biases and selection biases (genetic drift). The various settings offered by GAtor increase the likelihood of locating numerous low-energy minima, including those located in disconnected, hard to reach regions of the potential energy landscape. The best structures generated are re-relaxed and re-ranked using a hierarchy of increasingly accurate DFT functionals and dispersion methods. GAtor is applied to a chemically diverse set of four past blind test targets, characterized by different types of intermolecular interactions. The experimentally observed structures and other low-energy structures are found for all four targets. In particular, for Target II, 5-cyano-3-hydroxythiophene, the top ranked putative crystal structure is a Z' = 2 structure with P1̅ symmetry and a scaffold packing motif, which has not been reported previously.
Collapse
Affiliation(s)
- Farren Curtis
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Xiayue Li
- Google , Mountain View , California 94030 , United States.,Department of Materials Science and Engineering , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Timothy Rose
- Department of Materials Science and Engineering , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| | - Álvaro Vázquez-Mayagoitia
- Argonne Leadership Computing Facility , Argonne National Laboratory , Lemont , Illinois 60439 , United States
| | - Saswata Bhattacharya
- Department of Physics , Indian Institute of Technology Delhi , Hauz Khas , New Delhi 110016 , India
| | - Luca M Ghiringhelli
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6 , 14195 , Berlin , Germany
| | - Noa Marom
- Department of Physics , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States.,Department of Materials Science and Engineering , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States.,Department of Chemistry , Carnegie Mellon University , Pittsburgh , Pennsylvania 15213 , United States
| |
Collapse
|
168
|
Brandenburg JG, Bannwarth C, Hansen A, Grimme S. B97-3c: A revised low-cost variant of the B97-D density functional method. J Chem Phys 2018; 148:064104. [DOI: 10.1063/1.5012601] [Citation(s) in RCA: 255] [Impact Index Per Article: 42.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/22/2023] Open
Affiliation(s)
- Jan Gerit Brandenburg
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AH, United Kingdom
- Thomas Young Centre, University College London, Gower Street, London WC1E 6BT, United Kingdom
| | - Christoph Bannwarth
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Rheinische Friedrich-Wilhelms Universität Bonn, Beringstraße 4, 53115 Bonn, Germany
| |
Collapse
|
169
|
Thomas SP, Spackman PR, Jayatilaka D, Spackman MA. Accurate Lattice Energies for Molecular Crystals from Experimental Crystal Structures. J Chem Theory Comput 2018; 14:1614-1623. [DOI: 10.1021/acs.jctc.7b01200] [Citation(s) in RCA: 112] [Impact Index Per Article: 18.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Sajesh P. Thomas
- School of Molecular Sciences, University of Western Australia, Perth 6009, Australia
| | - Peter R. Spackman
- School of Molecular Sciences, University of Western Australia, Perth 6009, Australia
| | - Dylan Jayatilaka
- School of Molecular Sciences, University of Western Australia, Perth 6009, Australia
| | - Mark A. Spackman
- School of Molecular Sciences, University of Western Australia, Perth 6009, Australia
| |
Collapse
|
170
|
Abstract
Computational approaches based on the fundamental laws of quantum mechanics are now integral to almost all materials design initiatives in academia and industry. If computational materials science is genuinely going to deliver on its promises, then an electronic structure method with consistently high accuracy is urgently needed. We show that, thanks to recent algorithmic advances and the strategy developed in our manuscript, quantum Monte Carlo yields extremely accurate predictions for the lattice energies of materials at a surprisingly modest computational cost. It is thus no longer a technique that requires a world-leading computational facility to obtain meaningful results. While we focus on molecular crystals, the significance of our findings extends to all classes of materials. Computer simulation plays a central role in modern-day materials science. The utility of a given computational approach depends largely on the balance it provides between accuracy and computational cost. Molecular crystals are a class of materials of great technological importance which are challenging for even the most sophisticated ab initio electronic structure theories to accurately describe. This is partly because they are held together by a balance of weak intermolecular forces but also because the primitive cells of molecular crystals are often substantially larger than those of atomic solids. Here, we demonstrate that diffusion quantum Monte Carlo (DMC) delivers subchemical accuracy for a diverse set of molecular crystals at a surprisingly moderate computational cost. As such, we anticipate that DMC can play an important role in understanding and predicting the properties of a large number of molecular crystals, including those built from relatively large molecules which are far beyond reach of other high-accuracy methods.
Collapse
|
171
|
Solid-State Testing of a Van-Der-Waals-Corrected Exchange-Correlation Functional Based on the Semiclassical Atom Theory. COMPUTATION 2018. [DOI: 10.3390/computation6010007] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
172
|
Mortazavi M, Brandenburg JG, Maurer RJ, Tkatchenko A. Structure and Stability of Molecular Crystals with Many-Body Dispersion-Inclusive Density Functional Tight Binding. J Phys Chem Lett 2018; 9:399-405. [PMID: 29298075 DOI: 10.1021/acs.jpclett.7b03234] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/07/2023]
Abstract
Accurate prediction of structure and stability of molecular crystals is crucial in materials science and requires reliable modeling of long-range dispersion interactions. Semiempirical electronic structure methods are computationally more efficient than their ab initio counterparts, allowing structure sampling with significant speedups. We combine the Tkatchenko-Scheffler van der Waals method (TS) and the many-body dispersion method (MBD) with third-order density functional tight-binding (DFTB3) via a charge population-based method. We find an overall good performance for the X23 benchmark database of molecular crystals, despite an underestimation of crystal volume that can be traced to the DFTB parametrization. We achieve accurate lattice energy predictions with DFT+MBD energetics on top of vdW-inclusive DFTB3 structures, resulting in a speedup of up to 3000 times compared with a full DFT treatment. This suggests that vdW-inclusive DFTB3 can serve as a viable structural prescreening tool in crystal structure prediction.
Collapse
Affiliation(s)
- Majid Mortazavi
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
| | - Jan Gerit Brandenburg
- Department of Chemistry, University College London , 20 Gordon Street, WC1H 0AJ London, United Kingdom
- London Centre for Nanotechnology, University College London , 17-19 Gordon Street, WC1H 0AJ London, United Kingdom
- Thomas Young Centre, University College London , Gower Street, WC1E 6BT London, United Kingdom
| | - Reinhard J Maurer
- Department of Chemistry and Centre for Scientific Computing, University of Warwick , Gibbet Hill Road, Coventry CV4 7AL, United Kingdom
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
- Physics and Materials Science Research Unit, University of Luxembourg , L-1511, Luxembourg
| |
Collapse
|
173
|
Usvyat D, Maschio L, Schütz M. Periodic and fragment models based on the local correlation approach. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2018. [DOI: 10.1002/wcms.1357] [Citation(s) in RCA: 25] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/15/2022]
Affiliation(s)
- Denis Usvyat
- Institut für ChemieHumboldt‐Universität zu BerlinBerlinGermany
| | - Lorenzo Maschio
- Dipartimento di Chimica and NIS (Nanostructured Interfaces and Surfaces) CentreUniversità di TorinoTorinoItaly
| | - Martin Schütz
- Institut für ChemieHumboldt‐Universität zu BerlinBerlinGermany
| |
Collapse
|
174
|
Dolgonos GA, Loboda OA, Boese AD. Development of Embedded and Performance of Density Functional Methods for Molecular Crystals. J Phys Chem A 2018; 122:708-713. [PMID: 29265819 DOI: 10.1021/acs.jpca.7b12467] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We report an alternative quantum mechanical:quantum mechanical (QM:QM) method to the currently used periodic density functional calculations including dispersion and investigate its performance with respect to main structural and energetic properties of the X23 set of molecular crystals. By setting the goal of reproducing reference periodic BLYP+D3 values and by embedding BLYP+D3 into DFTB, we obtain results similar to those of periodic BLYP+D3-typically within 1-2% in lattice energies and ∼0.4% in cell volumes. The accuracy of this QM:QM method in comparison to DFTB+D and DFT+D for the X23 set of molecular crystals is discussed.
Collapse
Affiliation(s)
- Grygoriy A Dolgonos
- Institute of Chemistry, University of Graz , Heinrichstrasse 28/IV, 8010 Graz, Austria
| | - Oleksandr A Loboda
- Institute of Chemistry, University of Graz , Heinrichstrasse 28/IV, 8010 Graz, Austria
| | - A Daniel Boese
- Institute of Chemistry, University of Graz , Heinrichstrasse 28/IV, 8010 Graz, Austria
| |
Collapse
|
175
|
Iuzzolino L, McCabe P, Price SL, Brandenburg JG. Crystal structure prediction of flexible pharmaceutical-like molecules: density functional tight-binding as an intermediate optimisation method and for free energy estimation. Faraday Discuss 2018; 211:275-296. [PMID: 30035288 DOI: 10.1039/c8fd00010g] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Successful methodologies for theoretical crystal structure prediction (CSP) on flexible pharmaceutical-like organic molecules explore the lattice energy surface to find a set of plausible crystal structures. The initial search stages of CSP studies use relatively simple lattice energy approximations as hundreds of thousands of minima have to be considered. These generated crystal structures often have poor molecular geometries, as well as inaccurate lattice energy rankings, and performing reasonably accurate but computationally affordable optimisations of the crystal structures generated in a search would be highly desirable. Here, we seek to explore whether semi-empirical quantum-mechanical methods can perform this task. We employed the dispersion-corrected tight-binding Hamiltonian (DFTB3-D3) to relax all the inter- and intra-molecular degrees of freedom of several thousands of generated crystal structures of five pharmaceutical-like molecules, saving a large amount of computational effort compared to earlier studies. The computational cost scales better with molecular size and flexibility than other CSP methods, suggesting that it could be extended to even larger and more flexible molecules. On average, this optimisation improved the average reproduction of the eight experimental crystal structures (RMSD15) and experimental conformers (RMSD1) by 4% and 23%, respectively. The intermolecular interactions were then further optimised using distributed multipoles, derived from the molecular wave-functions, to accurately describe the electrostatic components of the intermolecular energies. In all cases, the experimental crystal structures are close to the top of the lattice energy ranking. Phonon calculations on some of the lowest energy structures were also performed with DFTB3-D3 methods to calculate the vibrational component of the Helmholtz free energy, providing further insights into the solid-state behaviour of the target molecules. We conclude that DFTB3-D3 is a cost-effective method for optimising flexible molecules, bridging the gap between the approximate methods used in CSP searches for generating crystal structures and more accurate methods required in the final energy ranking.
Collapse
Affiliation(s)
- Luca Iuzzolino
- Department of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, UK.
| | | | | | | |
Collapse
|
176
|
Custodio JMF, Michelini LJ, de Castro MRC, Vaz WF, Neves BJ, Cravo PVL, Barreto FS, Filho MOM, Perez CN, Napolitano HB. Structural insights into a novel anticancer sulfonamide chalcone. NEW J CHEM 2018. [DOI: 10.1039/c7nj03523c] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Natural products have stood out due to their wide range of biological activities.
Collapse
Affiliation(s)
- Jean M. F. Custodio
- Instituto de Química
- Universidade Federal de Goiás
- Goiânia
- Brazil
- Ciências Exatas e Tecnológicas
| | | | | | - Wesley F. Vaz
- Ciências Exatas e Tecnológicas
- Universidade Estadual de Goiás
- Anápolis
- Brazil
| | - Bruno J. Neves
- Programa de Pós-Graduação em Sociedade
- Tecnologia e Meio Ambiente
- Centro Universitário de Anápolis
- Anápolis
- Brazil
| | - Pedro V. L. Cravo
- Programa de Pós-Graduação em Sociedade
- Tecnologia e Meio Ambiente
- Centro Universitário de Anápolis
- Anápolis
- Brazil
| | | | | | | | - Hamilton B. Napolitano
- Ciências Exatas e Tecnológicas
- Universidade Estadual de Goiás
- Anápolis
- Brazil
- Programa de Pós-Graduação em Sociedade
| |
Collapse
|
177
|
Curtis F, Rose T, Marom N. Evolutionary niching in the GAtor genetic algorithm for molecular crystal structure prediction. Faraday Discuss 2018; 211:61-77. [DOI: 10.1039/c8fd00067k] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
The effects of evolutionary niching are investigated for the crystal structure prediction of 1,3-dibromo-2-chloro-5-fluorobenzene.
Collapse
Affiliation(s)
- Farren Curtis
- Department of Materials Science and Engineering
- Carnegie Mellon University
- Pittsburgh
- USA
- Department of Physics
| | - Timothy Rose
- Department of Materials Science and Engineering
- Carnegie Mellon University
- Pittsburgh
- USA
| | - Noa Marom
- Department of Materials Science and Engineering
- Carnegie Mellon University
- Pittsburgh
- USA
- Department of Physics
| |
Collapse
|
178
|
McKinley JL, Beran GJO. Identifying pragmatic quasi-harmonic electronic structure approaches for modeling molecular crystal thermal expansion. Faraday Discuss 2018; 211:181-207. [PMID: 30027972 DOI: 10.1039/c8fd00048d] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Quasi-harmonic approaches provide an economical route to modeling the temperature dependence of molecular crystal structures and properties. Several studies have demonstrated good performance of these models, at least for rigid molecules, when using fragment-based approaches with correlated wavefunction techniques. Many others have found success employing dispersion-corrected density functional theory (DFT). Here, a hierarchy of models in which the energies, geometries, and phonons are computed either with correlated methods or DFT are examined to identify which combinations produce useful predictions for properties such as the molar volume, enthalpy, and entropy as a function of temperature. The results demonstrate that refining DFT geometries and phonons with single-point energies based on dispersion-corrected second-order Møller-Plesset perturbation theory can provide clear improvements in the molar volumes and enthalpies compared to those obtained from DFT alone. Predicted entropies, which are governed by vibrational contributions, benefit less clearly from the hybrid schemes. Using these hybrid techniques, the room-temperature thermochemistry of acetaminophen (paracetamol) is predicted to address the discrepancy between two experimental sublimation enthalpy measurements.
Collapse
Affiliation(s)
- Jessica L McKinley
- Department of Chemistry, University of California, Riverside, California 92521, USA.
| | | |
Collapse
|
179
|
Li Q, Zaczek AJ, Korter TM, Zeitler JA, Ruggiero MT. Methyl-rotation dynamics in metal–organic frameworks probed with terahertz spectroscopy. Chem Commun (Camb) 2018; 54:5776-5779. [DOI: 10.1039/c8cc02650e] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
In ZIF-8 and its cobalt analogue ZIF-67, the imidazolate methyl-groups, which point directly into the void space, have been shown to freely rotate – even down to cryogenic temperatures.
Collapse
Affiliation(s)
- Qi Li
- Department of Chemical Engineering and Biotechnology, University of Cambridge
- Philippa Fawcett Drive
- Cambridge
- UK
| | - Adam J. Zaczek
- Department of Chemistry, Syracuse University
- 1-014 Center for Science and Technology
- Syracuse
- USA
| | - Timothy M. Korter
- Department of Chemistry, Syracuse University
- 1-014 Center for Science and Technology
- Syracuse
- USA
| | - J. Axel Zeitler
- Department of Chemical Engineering and Biotechnology, University of Cambridge
- Philippa Fawcett Drive
- Cambridge
- UK
| | - Michael T. Ruggiero
- Department of Chemical Engineering and Biotechnology, University of Cambridge
- Philippa Fawcett Drive
- Cambridge
- UK
- Department of Chemistry
| |
Collapse
|
180
|
Tan M, Shtukenberg AG, Zhu S, Xu W, Dooryhee E, Nichols S, Ward MD, Kahr B, Zhu Q. ROY revisited, again: the eighth solved structure. Faraday Discuss 2018; 211:477-491. [DOI: 10.1039/c8fd00039e] [Citation(s) in RCA: 36] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
X-ray powder diffraction and crystal structure prediction algorithms are used in synergy to establish the crystal structure of the eighth polymorph of ROY, form R05.
Collapse
Affiliation(s)
- Melissa Tan
- Department of Chemistry
- Molecular Design Institute
- New York University
- New York City
- USA
| | | | - Shengcai Zhu
- Department of Physics and Astronomy
- High Pressure Science and Engineering Center
- University of Nevada Las Vegas
- USA
| | - Wenqian Xu
- X-ray Science Division
- Advanced Photon Source
- Argonne National Laboratory
- Argonne
- USA
| | - Eric Dooryhee
- Photon Sciences Division
- National Synchrotron Light Source II
- Brookhaven National Laboratory
- Upton
- USA
| | - Shane M. Nichols
- Department of Chemistry
- Molecular Design Institute
- New York University
- New York City
- USA
| | - Michael D. Ward
- Department of Chemistry
- Molecular Design Institute
- New York University
- New York City
- USA
| | - Bart Kahr
- Department of Chemistry
- Molecular Design Institute
- New York University
- New York City
- USA
| | - Qiang Zhu
- Department of Physics and Astronomy
- High Pressure Science and Engineering Center
- University of Nevada Las Vegas
- USA
| |
Collapse
|
181
|
Ruggiero MT, Kölbel J, Li Q, Zeitler JA. Predicting the structures and associated phase transition mechanisms in disordered crystals via a combination of experimental and theoretical methods. Faraday Discuss 2018; 211:425-439. [DOI: 10.1039/c8fd00042e] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
Experimental terahertz time-domain spectroscopy and theoretical solid-state ab initio density functional theory and molecular dynamics simulations are used to elucidate the structures, dynamics, and phase transformation processes of molecular crystals undergoing a solid-state order–disorder transition.
Collapse
Affiliation(s)
- Michael T. Ruggiero
- Department of Chemical Engineering and Biotechnology
- University of Cambridge
- Cambridge
- UK
| | - Johanna Kölbel
- Department of Chemical Engineering and Biotechnology
- University of Cambridge
- Cambridge
- UK
| | - Qi Li
- Department of Chemical Engineering and Biotechnology
- University of Cambridge
- Cambridge
- UK
| | - J. Axel Zeitler
- Department of Chemical Engineering and Biotechnology
- University of Cambridge
- Cambridge
- UK
| |
Collapse
|
182
|
Hoja J, Tkatchenko A. First-principles stability ranking of molecular crystal polymorphs with the DFT+MBD approach. Faraday Discuss 2018; 211:253-274. [PMID: 30042995 DOI: 10.1039/c8fd00066b] [Citation(s) in RCA: 32] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022]
Abstract
The ability to accurately calculate the relative stabilities of numerous polymorphs of a given molecular crystal is crucial for the success of any molecular crystal structure prediction (CSP) approach. We have recently presented a hierarchical CSP procedure based on van-der-Waals-inclusive density functional theory [Hoja et al., 2018, arXiv:1803.07503], which yields excellent stability rankings for molecular crystals involving rigid molecules, salts, co-crystals, and highly polymorphic drug-like molecules. This approach includes many-body dispersion effects, exact exchange, as well as vibrational free energies. Here, we discuss in detail the impact of these effects on the obtained stability rankings. In addition, we assess the impact of the approximations used in our hierarchical procedure. We show that our procedure is generally robust to 1-2 kJ mol-1 for the systems in the latest CSP blind test but vibrational free energies for crystals involving flexible molecules would benefit from directly including many-body dispersion interactions. In addition, we also discuss the effect of temperature on the structure of molecular crystals and a simple but effective method for estimating anharmonic effects.
Collapse
Affiliation(s)
- Johannes Hoja
- Physics and Materials Science Research Unit, University of Luxembourg, L-1511 Luxembourg, Luxembourg.
| | | |
Collapse
|
183
|
Yuan X, Luo K, Liu N, Ji X, Liu C, He J, Tian G, Zhao Y, Yu D. Cluster-model DFT simulations of the infrared spectra of triazine-based molecular crystals. Phys Chem Chem Phys 2018; 20:20779-20784. [DOI: 10.1039/c8cp01550c] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
A cluster-model approach is developed to simulate the IR spectra of triazine-based molecular crystals, and the distinct hydrogen-bonding environments of the crystallographically independent molecules can be clearly revealed.
Collapse
Affiliation(s)
- Xiaohong Yuan
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
| | - Kun Luo
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
- Hebei Key Laboratory of Microstructural Material Physics
| | - Nan Liu
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
| | - Xueqiang Ji
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
| | - Chao Liu
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
| | - Julong He
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
| | - Guangjun Tian
- Hebei Key Laboratory of Microstructural Material Physics
- School of Science
- Yanshan University
- Qinhuangdao 066004
- China
| | - Yuanchun Zhao
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
| | - Dongli Yu
- State Key Laboratory of Metastable Materials Science and Technology
- Yanshan University
- Qinhuangdao 066004
- China
| |
Collapse
|
184
|
Spackman MA. Towards the use of experimental electron densities to estimate reliable lattice energies. CrystEngComm 2018. [DOI: 10.1039/c8ce01108g] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
Lattice energies derived from experimental charge densities are critically assessed, with a view to encouraging further research of this nature.
Collapse
Affiliation(s)
- Mark A. Spackman
- School of Molecular Sciences
- University of Western Australia
- Perth
- Australia
| |
Collapse
|
185
|
Mohebifar M, Johnson ER, Rowley CN. Evaluating Force-Field London Dispersion Coefficients Using the Exchange-Hole Dipole Moment Model. J Chem Theory Comput 2017; 13:6146-6157. [DOI: 10.1021/acs.jctc.7b00522] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Mohamad Mohebifar
- Department
of Chemistry, Memorial University of Newfoundland, St. John’s, Newfoundland
and Labrador, Canada
| | - Erin R. Johnson
- Department
of Chemistry, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Christopher N. Rowley
- Department
of Chemistry, Memorial University of Newfoundland, St. John’s, Newfoundland
and Labrador, Canada
| |
Collapse
|
186
|
Hartman JD, Balaji A, Beran GJO. Improved Electrostatic Embedding for Fragment-Based Chemical Shift Calculations in Molecular Crystals. J Chem Theory Comput 2017; 13:6043-6051. [PMID: 29139294 DOI: 10.1021/acs.jctc.7b00677] [Citation(s) in RCA: 30] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/01/2023]
Abstract
Fragment-based methods predict nuclear magnetic resonance (NMR) chemical shielding tensors in molecular crystals with high accuracy and computational efficiency. Such methods typically employ electrostatic embedding to mimic the crystalline environment, and the quality of the results can be sensitive to the embedding treatment. To improve the quality of this embedding environment for fragment-based molecular crystal property calculations, we borrow ideas from the embedded ion method to incorporate self-consistently polarized Madelung field effects. The self-consistent reproduction of the Madelung potential (SCRMP) model developed here constructs an array of point charges that incorporates self-consistent lattice polarization and which reproduces the Madelung potential at all atomic sites involved in the quantum mechanical region of the system. The performance of fragment- and cluster-based 1H, 13C, 14N, and 17O chemical shift predictions using SCRMP and density functionals like PBE and PBE0 are assessed. The improved embedding model results in substantial improvements in the predicted 17O chemical shifts and modest improvements in the 15N ones. Finally, the performance of the model is demonstrated by examining the assignment of the two oxygen chemical shifts in the challenging γ-polymorph of glycine. Overall, the SCRMP-embedded NMR chemical shift predictions are on par with or more accurate than those obtained with the widely used gauge-including projector augmented wave (GIPAW) model.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California, Riverside , Riverside, California 92521, United States
| | - Ashwin Balaji
- Department of Chemistry, University of California, Riverside , Riverside, California 92521, United States
| | - Gregory J O Beran
- Department of Chemistry, University of California, Riverside , Riverside, California 92521, United States
| |
Collapse
|
187
|
Aina AA, Misquitta AJ, Price SL. From dimers to the solid-state: Distributed intermolecular force-fields for pyridine. J Chem Phys 2017; 147:161722. [DOI: 10.1063/1.4999789] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Alexander A. Aina
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| | - Alston J. Misquitta
- School of Physics and Astronomy, Queen Mary, University of London, London E1 4NS, United Kingdom
| | - Sarah L. Price
- Department of Chemistry, University College London, London WC1H 0AJ, United Kingdom
| |
Collapse
|
188
|
Meded V, Knorr N, Neumann T, Nelles G, Wenzel W, von Wrochem F. Structural origins of the cohesive energy in metal-terpyridine oligomer thin-films. Phys Chem Chem Phys 2017; 19:27952-27959. [PMID: 28949351 DOI: 10.1039/c7cp05488b] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
FeII-terpyridine based oligomers have attracted considerable interest as key constituents for the realization of highly robust, ultra-thin ordered layers of metal center oligomers (MCOs) for organic electronics applications. By using molecular simulations and nanotribology investigations, we report on the origins of the surprisingly high mechanical and thermal stability in this type of MCO layers, which finds its expression in nanowear resistance values of up to 1.5 μN for the MCO films, as well as in a thermal stability of two-terminal MCO junctions to temperatures up to ∼100 °C under electrical load. A theoretical analysis of the fundamental cohesive forces among the constituents within the context of an electrostatic model reveal that the cohesive energy is essentially based on Coulomb interactions among the ionic constituents of the oligomers, leading to an estimated cohesive energy per molar mass of 0.0132 eV mol g-1 for MCO layers that advantageously compare to the 0.0061 eV mol g-1 reported for pentacene crystals.
Collapse
Affiliation(s)
- Velimir Meded
- Institute of Nanotechnology, Karlsruhe Institute of Technology, 76021 Karlsruhe, Germany.
| | | | | | | | | | | |
Collapse
|
189
|
Zilka M, Sturniolo S, Brown SP, Yates JR. Visualising crystal packing interactions in solid-state NMR: Concepts and applications. J Chem Phys 2017; 147:144203. [DOI: 10.1063/1.4996750] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Affiliation(s)
- Miri Zilka
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Simone Sturniolo
- Scientific Computing Department, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX,
United Kingdom
| | - Steven P. Brown
- Department of Physics, University of Warwick, Coventry CV4 7AL, United Kingdom
| | - Jonathan R. Yates
- Department of Materials, University of Oxford, Oxford OX1 3PH, United Kingdom
| |
Collapse
|
190
|
Zilka M, Dudenko DV, Hughes CE, Williams PA, Sturniolo S, Franks WT, Pickard CJ, Yates JR, Harris KDM, Brown SP. Ab initio random structure searching of organic molecular solids: assessment and validation against experimental data. Phys Chem Chem Phys 2017; 19:25949-25960. [PMID: 28944393 PMCID: PMC5779078 DOI: 10.1039/c7cp04186a] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Accepted: 08/30/2017] [Indexed: 02/03/2023]
Abstract
This paper explores the capability of using the DFT-D ab initio random structure searching (AIRSS) method to generate crystal structures of organic molecular materials, focusing on a system (m-aminobenzoic acid; m-ABA) that is known from experimental studies to exhibit abundant polymorphism. Within the structural constraints selected for the AIRSS calculations (specifically, centrosymmetric structures with Z = 4 for zwitterionic m-ABA molecules), the method is shown to successfully generate the two known polymorphs of m-ABA (form III and form IV) that have these structural features. We highlight various issues that are encountered in comparing crystal structures generated by AIRSS to experimental powder X-ray diffraction (XRD) data and solid-state magic-angle spinning (MAS) NMR data, demonstrating successful fitting for some of the lowest energy structures from the AIRSS calculations against experimental low-temperature powder XRD data for known polymorphs of m-ABA, and showing that comparison of computed and experimental solid-state NMR parameters allows different hydrogen-bonding motifs to be discriminated.
Collapse
Affiliation(s)
- Miri Zilka
- Department of Physics, University of Warwick, Coventry CV4 7AL, UK.
| | - Dmytro V Dudenko
- Department of Physics, University of Warwick, Coventry CV4 7AL, UK. and School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK.
| | - Colan E Hughes
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK.
| | - P Andrew Williams
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK.
| | - Simone Sturniolo
- Scientific Computing Department, Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire OX11 0QX, UK
| | - W Trent Franks
- Department of Physics, University of Warwick, Coventry CV4 7AL, UK.
| | - Chris J Pickard
- Department of Materials Science & Metallurgy, University of Cambridge, 27 Charles Babbage Road, Cambridge CB3 0FS, UK
| | - Jonathan R Yates
- Department of Materials, University of Oxford, Oxford OX1 3PH, UK.
| | - Kenneth D M Harris
- School of Chemistry, Cardiff University, Park Place, Cardiff CF10 3AT, UK.
| | - Steven P Brown
- Department of Physics, University of Warwick, Coventry CV4 7AL, UK.
| |
Collapse
|
191
|
Vener MV, Chernyshov IY, Rykounov AA, Filarowski A. Structural and spectroscopic features of proton hydrates in the crystalline state. Solid-state DFT study on HCl and triflic acid hydrates. Mol Phys 2017. [DOI: 10.1080/00268976.2017.1380860] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- M. V. Vener
- Quantum Chemistry Department, Mendeleev University of Chemical Technology, Moscow, Russia
- Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, Moscow, Russia
| | - I. Yu. Chernyshov
- Quantum Chemistry Department, Mendeleev University of Chemical Technology, Moscow, Russia
- Kurnakov Institute of General and Inorganic Chemistry, Russian Academy of Sciences, Moscow, Russia
| | - A. A. Rykounov
- Theoretical Department, Russian Federal Nuclear Center – All-Russian Research Institute of Technical Physics (RFNC-VNIITF), Snezhinsk, Russia
| | - A. Filarowski
- Faculty of Chemistry, University of Wrocław, Wrocław, Poland
- Department of Physics, Industrial University of Tyumen, Tyumen, Russia
| |
Collapse
|
192
|
Sontising W, Heit YN, McKinley JL, Beran GJO. Theoretical predictions suggest carbon dioxide phases III and VII are identical. Chem Sci 2017; 8:7374-7382. [PMID: 29163888 PMCID: PMC5672840 DOI: 10.1039/c7sc03267f] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2017] [Accepted: 09/03/2017] [Indexed: 11/21/2022] Open
Abstract
Accurate electronic structure calculations for the structures and simulated Raman spectra of high-pressure carbon dioxide suggest phases III and VII are identical, and the phase diagram should be revised.
Solid carbon dioxide exhibits a rich phase diagram at high pressures. Metastable phase III is formed by compressing dry ice above ∼10–12 GPa. Phase VII occurs at similar pressures but higher temperatures, and its stability region is disconnected from III on the phase diagram. Comparison of large-basis-set quasi-harmonic second-order Møller–Plesset perturbation theory calculations and experiment suggests that the long-accepted structure of phase III is problematic. The experimental phase III and VII structures both relax to the same phase VII structure. Furthermore, Raman spectra predicted for phase VII are in good agreement with those observed experimentally for both phase III and VII, while those for the purported phase III structure agree poorly with experimental observations. Crystal structure prediction is employed to search for other potential structures which might account for phase III, but none are found. Together, these results suggest that phases III and VII are likely identical.
Collapse
Affiliation(s)
- Watit Sontising
- Department of Chemistry , University of California , Riverside , California 92521 , USA .
| | - Yonaton N Heit
- Department of Chemistry , University of California , Riverside , California 92521 , USA .
| | - Jessica L McKinley
- Department of Chemistry , University of California , Riverside , California 92521 , USA .
| | - Gregory J O Beran
- Department of Chemistry , University of California , Riverside , California 92521 , USA .
| |
Collapse
|
193
|
Liu L, Brandenburg JG, Grimme S. On the hydrogen activation by frustrated Lewis pairs in the solid state: benchmark studies and theoretical insights. PHILOSOPHICAL TRANSACTIONS. SERIES A, MATHEMATICAL, PHYSICAL, AND ENGINEERING SCIENCES 2017; 375:20170006. [PMID: 28739964 PMCID: PMC5540839 DOI: 10.1098/rsta.2017.0006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Accepted: 02/15/2017] [Indexed: 05/26/2023]
Abstract
Recently, the concept of small molecule activation by frustrated Lewis pairs (FLPs) has been expanded to the solid state showing a variety of interesting reactivities. Therefore, there is a need to establish a computational protocol to investigate such systems theoretically. In the present study, we selected several FLPs and applied multiple levels of theory, ranging from a semi-empirical tight-binding Hamiltonian to dispersion corrected hybrid density functionals. Their performance is benchmarked for the computation of crystal geometries, thermostatistical contributions, and reaction energies. We show that the computationally efficient HF-3c method gives accurate crystal structures and is numerically stable and sufficiently fast for routine applications. This method also gives reliable values for the thermostatistical contributions to Gibbs free energies. The meta-generalized gradient approximated TPSS-D3 evaluated in a projector augmented plane wave basis set is able to produce sufficiently accurate reaction electronic energies. The established protocol is intended to support experimental studies and to predict new reactions in the emerging field of solid-state FLPs.This article is part of the themed issue 'Frustrated Lewis pair chemistry'.
Collapse
Affiliation(s)
- Lei Liu
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstrasse 4, 53115 Bonn, Germany
| | - Jan Gerit Brandenburg
- London Centre for Nanotechnology, University College London, 17-19 Gordon Street, WC1H 0AH London, UK
- Thomas Young Centre, University College London, Gower Street, WC1E 6BT London, UK
| | - Stefan Grimme
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstrasse 4, 53115 Bonn, Germany
| |
Collapse
|
194
|
Hapka M, Rajchel Ł, Modrzejewski M, Schäffer R, Chałasiński G, Szczęśniak MM. The nature of three-body interactions in DFT: Exchange and polarization effects. J Chem Phys 2017; 147:084106. [PMID: 28863509 DOI: 10.1063/1.4986291] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We propose a physically motivated decomposition of density functional theory (DFT) 3-body nonadditive interaction energies into the exchange and density-deformation (polarization) components. The exchange component represents the effect of the Pauli exclusion in the wave function of the trimer and is found to be challenging for density functional approximations (DFAs). The remaining density-deformation nonadditivity is less dependent upon the DFAs. Numerical demonstration is carried out for rare gas atom trimers, Ar2-HX (X = F, Cl) complexes, and small hydrogen-bonded and van der Waals molecular systems. None of the tested semilocal, hybrid, and range-separated DFAs properly accounts for the nonadditive exchange in dispersion-bonded trimers. By contrast, for hydrogen-bonded systems, range-separated DFAs achieve a qualitative agreement to within 20% of the reference exchange energy. A reliable performance for all systems is obtained only when the monomers interact through the Hartree-Fock potential in the dispersion-free Pauli blockade scheme. Additionally, we identify the nonadditive second-order exchange-dispersion energy as an important but overlooked contribution in force-field-like dispersion corrections. Our results suggest that range-separated functionals do not include this component, although semilocal and global hybrid DFAs appear to imitate it in the short range.
Collapse
Affiliation(s)
- Michał Hapka
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Łukasz Rajchel
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
| | - Marcin Modrzejewski
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | - Rainer Schäffer
- Faculty of Chemistry, University of Duisburg-Essen, Universitätsstraße 5, 45117 Essen, Germany
| | - Grzegorz Chałasiński
- Faculty of Chemistry, University of Warsaw, ul. L. Pasteura 1, 02-093 Warsaw, Poland
| | | |
Collapse
|
195
|
Shtukenberg AG, Zhu Q, Carter DJ, Vogt L, Hoja J, Schneider E, Song H, Pokroy B, Polishchuk I, Tkatchenko A, Oganov AR, Rohl AL, Tuckerman ME, Kahr B. Powder diffraction and crystal structure prediction identify four new coumarin polymorphs. Chem Sci 2017; 8:4926-4940. [PMID: 28959416 PMCID: PMC5607859 DOI: 10.1039/c7sc00168a] [Citation(s) in RCA: 64] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2017] [Accepted: 04/30/2017] [Indexed: 11/21/2022] Open
Abstract
Coumarin, a simple, commodity chemical isolated from beans in 1820, has, to date, only yielded one solid state structure. Here, we report a rich polymorphism of coumarin grown from the melt. Four new metastable forms were identified and their crystal structures were solved using a combination of computational crystal structure prediction algorithms and X-ray powder diffraction. With five crystal structures, coumarin has become one of the few rigid molecules showing extensive polymorphism at ambient conditions. We demonstrate the crucial role of advanced electronic structure calculations including many-body dispersion effects for accurate ranking of the stability of coumarin polymorphs and the need to account for anharmonic vibrational contributions to their free energy. As such, coumarin is a model system for studying weak intermolecular interactions, crystallization mechanisms, and kinetic effects.
Collapse
Affiliation(s)
- Alexander G Shtukenberg
- Department of Chemistry , Molecular Design Institute , New York University , New York City , NY 10003 , USA .
| | - Qiang Zhu
- Department of Physics and Astronomy , High Pressure Science and Engineering Center , University of Nevada Las Vegas , Nevada 89154 , USA .
- Department of Geosciences , Stony Brook University , Stony Brook , NY 11794 , USA
| | - Damien J Carter
- Curtin Institute for Computation and Department of Chemistry , Curtin University , P.O. Box U1987 , Perth , 6845 , Western Australia , Australia
| | - Leslie Vogt
- Department of Chemistry , New York University , New York City , NY 10003 , USA
| | - Johannes Hoja
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6 , Berlin , 14195 , Germany
- Physics and Materials Science Research Unit , University of Luxembourg , 1511 Luxembourg , Luxembourg
| | - Elia Schneider
- Department of Chemistry , New York University , New York City , NY 10003 , USA
| | - Hongxing Song
- Department of Chemistry , New York University , New York City , NY 10003 , USA
| | - Boaz Pokroy
- Department of Materials Science and Engineering , Russell Berrie Nanotechnology Institute , Technion Israel Institute of Technology , Haifa 32000 , Israel
| | - Iryna Polishchuk
- Department of Materials Science and Engineering , Russell Berrie Nanotechnology Institute , Technion Israel Institute of Technology , Haifa 32000 , Israel
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6 , Berlin , 14195 , Germany
- Physics and Materials Science Research Unit , University of Luxembourg , 1511 Luxembourg , Luxembourg
| | - Artem R Oganov
- Skolkovo Institute of Science and Technology , Skolkovo Innovation Center , 3 Nobel St. , Moscow 143026 , Russia
- Department of Geosciences , Stony Brook University , Stony Brook , NY 11794 , USA
| | - Andrew L Rohl
- Curtin Institute for Computation and Department of Chemistry , Curtin University , P.O. Box U1987 , Perth , 6845 , Western Australia , Australia
| | - Mark E Tuckerman
- Department of Chemistry , New York University , New York City , NY 10003 , USA
- Courant Institute of Mathematical Sciences , New York University , New York City , NY 10003 , USA
- New York University-East China Normal University Center for Computational Chemistry at NYU Shanghai , 3663 Zhongshan Road North , Shanghai 200062 , China
| | - Bart Kahr
- Department of Chemistry , Molecular Design Institute , New York University , New York City , NY 10003 , USA .
- Department of Advanced Science and Engineering (TWIns) , Waseda University , Wakamatsucho, 3-2 , Shinjuku , 162-0056 Tokyo , Japan
| |
Collapse
|
196
|
Soss SE, Flynn PF, Iuliucci RJ, Young RP, Mueller LJ, Hartman J, Beran GJO, Harper JK. Measuring and Modeling Highly Accurate
15
N Chemical Shift Tensors in a Peptide. Chemphyschem 2017; 18:2225-2232. [DOI: 10.1002/cphc.201700357] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2017] [Revised: 06/08/2017] [Indexed: 11/06/2022]
Affiliation(s)
- Sarah E. Soss
- Department of Chemistry University of Utah Salt Lake City UT 84112 USA
| | - Peter F. Flynn
- Department of Chemistry University of Utah Salt Lake City UT 84112 USA
| | - Robbie J. Iuliucci
- Department of Chemistry Washington and Jefferson College 60 Lincoln Street Washington PA 15301 USA
| | - Robert P. Young
- Department of Chemistry University of California Riverside CA 92521 USA
| | | | - Joshua Hartman
- Department of Chemistry University of California Riverside CA 92521 USA
| | | | - James K. Harper
- Department of Chemistry University of Central Florida 4111 Libra Drive Orlando FL 32816 USA
| |
Collapse
|
197
|
Deringer VL, George J, Dronskowski R, Englert U. Plane-Wave Density Functional Theory Meets Molecular Crystals: Thermal Ellipsoids and Intermolecular Interactions. Acc Chem Res 2017; 50:1231-1239. [PMID: 28467707 DOI: 10.1021/acs.accounts.7b00067] [Citation(s) in RCA: 41] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/28/2023]
Abstract
Molecular compounds, organic and inorganic, crystallize in diverse and complex structures. They continue to inspire synthetic efforts and "crystal engineering", with implications ranging from fundamental questions to pharmaceutical research. The structural complexity of molecular solids is linked with diverse intermolecular interactions: hydrogen bonding with all its facets, halogen bonding, and other secondary bonding mechanisms of recent interest (and debate). Today, high-resolution diffraction experiments allow unprecedented insight into the structures of molecular crystals. Despite their usefulness, however, these experiments also face problems: hydrogen atoms are challenging to locate, and thermal effects may complicate matters. Moreover, even if the structure of a crystal is precisely known, this does not yet reveal the nature and strength of the intermolecular forces that hold it together. In this Account, we show that periodic plane-wave-based density functional theory (DFT) can be a useful, and sometimes unexpected, complement to molecular crystallography. Initially developed in the solid-state physics communities to treat inorganic solids, periodic DFT can be applied to molecular crystals just as well: theoretical structural optimizations "help out" by accurately localizing the elusive hydrogen atoms, reaching neutron-diffraction quality with much less expensive measurement equipment. In addition, phonon computations, again developed by physicists, can quantify the thermal motion of atoms and thus predict anisotropic displacement parameters and ORTEP ellipsoids "from scratch". But the synergy between experiment and theory goes much further than that. Once a structure has been accurately determined, computations give new and detailed insights into the aforementioned intermolecular interactions. For example, it has been debated whether short hydrogen bonds in solids have covalent character, and we have added a new twist to this discussion using an orbital-based theory that once more had been developed for inorganic solids. However, there is more to a crystal structure than a handful of short contacts between neighboring residues. We hence have used dimensionally resolved analyses to dissect crystalline networks in a systematic fashion, one spatial direction at a time. Initially applied to hydrogen bonding, these techniques can be seamlessly extended to halogen, chalcogen, and pnictogen bonding, quantifying bond strength and cooperativity in truly infinite networks. Finally, these methods promise to be useful for (bio)polymers, as we have recently exemplified for α-chitin. At the interface of increasingly accurate and popular DFT methods, ever-improving crystallographic expertise, and new challenging, chemical questions, we believe that combined experimental and theoretical studies of molecular crystals are just beginning to pick up speed.
Collapse
Affiliation(s)
- Volker L. Deringer
- Institute
of Inorganic Chemistry and ‡Jülich−Aachen Research
Alliance (JARA-HPC), RWTH Aachen University, Landoltweg 1, 52056 Aachen, Germany
| | - Janine George
- Institute
of Inorganic Chemistry and ‡Jülich−Aachen Research
Alliance (JARA-HPC), RWTH Aachen University, Landoltweg 1, 52056 Aachen, Germany
| | - Richard Dronskowski
- Institute
of Inorganic Chemistry and ‡Jülich−Aachen Research
Alliance (JARA-HPC), RWTH Aachen University, Landoltweg 1, 52056 Aachen, Germany
| | - Ulli Englert
- Institute
of Inorganic Chemistry and ‡Jülich−Aachen Research
Alliance (JARA-HPC), RWTH Aachen University, Landoltweg 1, 52056 Aachen, Germany
| |
Collapse
|
198
|
Hermann J, DiStasio RA, Tkatchenko A. First-Principles Models for van der Waals Interactions in Molecules and Materials: Concepts, Theory, and Applications. Chem Rev 2017; 117:4714-4758. [PMID: 28272886 DOI: 10.1021/acs.chemrev.6b00446] [Citation(s) in RCA: 255] [Impact Index Per Article: 36.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/29/2023]
Abstract
Noncovalent van der Waals (vdW) or dispersion forces are ubiquitous in nature and influence the structure, stability, dynamics, and function of molecules and materials throughout chemistry, biology, physics, and materials science. These forces are quantum mechanical in origin and arise from electrostatic interactions between fluctuations in the electronic charge density. Here, we explore the conceptual and mathematical ingredients required for an exact treatment of vdW interactions, and present a systematic and unified framework for classifying the current first-principles vdW methods based on the adiabatic-connection fluctuation-dissipation (ACFD) theorem (namely the Rutgers-Chalmers vdW-DF, Vydrov-Van Voorhis (VV), exchange-hole dipole moment (XDM), Tkatchenko-Scheffler (TS), many-body dispersion (MBD), and random-phase approximation (RPA) approaches). Particular attention is paid to the intriguing nature of many-body vdW interactions, whose fundamental relevance has recently been highlighted in several landmark experiments. The performance of these models in predicting binding energetics as well as structural, electronic, and thermodynamic properties is connected with the theoretical concepts and provides a numerical summary of the state-of-the-art in the field. We conclude with a roadmap of the conceptual, methodological, practical, and numerical challenges that remain in obtaining a universally applicable and truly predictive vdW method for realistic molecular systems and materials.
Collapse
Affiliation(s)
- Jan Hermann
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany
| | - Robert A DiStasio
- Department of Chemistry and Chemical Biology, Cornell University , Ithaca, New York 14853, United States
| | - Alexandre Tkatchenko
- Fritz-Haber-Institut der Max-Planck-Gesellschaft , Faradayweg 4-6, 14195 Berlin, Germany.,Physics and Materials Science Research Unit, University of Luxembourg , L-1511 Luxembourg, Luxembourg
| |
Collapse
|
199
|
Li X, Neumann MA, van de Streek J. The application of tailor-made force fields and molecular dynamics for NMR crystallography: a case study of free base cocaine. IUCRJ 2017; 4:175-184. [PMID: 28250956 PMCID: PMC5330528 DOI: 10.1107/s2052252517001415] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/19/2016] [Accepted: 01/27/2017] [Indexed: 05/04/2023]
Abstract
Motional averaging has been proven to be significant in predicting the chemical shifts in ab initio solid-state NMR calculations, and the applicability of motional averaging with molecular dynamics has been shown to depend on the accuracy of the molecular mechanical force field. The performance of a fully automatically generated tailor-made force field (TMFF) for the dynamic aspects of NMR crystallography is evaluated and compared with existing benchmarks, including static dispersion-corrected density functional theory calculations and the COMPASS force field. The crystal structure of free base cocaine is used as an example. The results reveal that, even though the TMFF outperforms the COMPASS force field for representing the energies and conformations of predicted structures, it does not give significant improvement in the accuracy of NMR calculations. Further studies should direct more attention to anisotropic chemical shifts and development of the method of solid-state NMR calculations.
Collapse
Affiliation(s)
- Xiaozhou Li
- Department of Pharmacy, University of Copenhagen, Universitetsparken 2, Copenhagen DK-2100, Denmark
| | - Marcus A. Neumann
- Avant-garde Materials Simulation Deutschland GmbH, Rosa-Luxemberg-Strasse 14, Freiburg D-79100, Germany
| | - Jacco van de Streek
- Department of Pharmacy, University of Copenhagen, Universitetsparken 2, Copenhagen DK-2100, Denmark
| |
Collapse
|
200
|
Fang T, Li Y, Li S. Generalized energy‐based fragmentation approach for modeling condensed phase systems. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2017. [DOI: 10.1002/wcms.1297] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
Affiliation(s)
- Tao Fang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational ChemistryNanjing University Nanjing P. R. China
| | - Yunzhi Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational ChemistryNanjing University Nanjing P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of Education, Institute of Theoretical and Computational ChemistryNanjing University Nanjing P. R. China
| |
Collapse
|