QM7-X: A comprehensive dataset of quantum-mechanical properties spanning the chemical space of small organic molecules

We introduce QM7-X, a comprehensive dataset of 42 physicochemical properties for $\approx$ 4.2 M equilibrium and non-equilibrium structures of small organic molecules with up to seven non-hydrogen (C, N, O, S, Cl) atoms. To span this fundamentally important region of chemical compound space (CCS), QM7-X includes an exhaustive sampling of (meta-)stable equilibrium structures - comprised of constitutional/structural isomers and stereoisomers, e.g., enantiomers and diastereomers (including cis-/trans- and conformational isomers) - as well as 100 non-equilibrium structural variations thereof to reach a total of $\approx$ 4.2 M molecular structures. Computed at the tightly converged quantum-mechanical PBE0+MBD level of theory, QM7-X contains global (molecular) and local (atom-in-a-molecule) properties ranging from ground state quantities (such as atomization energies and dipole moments) to response quantities (such as polarizability tensors and dispersion coefficients). By providing a systematic, extensive, and tightly-converged dataset of quantum-mechanically computed physicochemical properties, we expect that QM7-X will play a critical role in the development of next-generation machine-learning based models for exploring greater swaths of CCS and performing in silico design of molecules with targeted properties

Tunable internal quantum well alignment in rationally designed oligomer-based perovskite films deposited by resonant infrared matrix-assisted pulsed laser evaporation

Hybrid perovskites incorporating conjugated organic cations enable unusual charge carrier interactions among organic and inorganic structural components, but are difficult to prepare as films due to disparate component chemical/physical characteristics (e.g., solubility, thermal stability). Here we demonstrate that resonant infrared matrix-assisted pulsed laser evaporation (RIR-MAPLE) mitigates these challenges, enabling facile deposition of lead-halide-based perovskite films incorporating variable-length oligothiophene cations. Density functional theory (DFT) predicts suitable organic and inorganic moieties that form quantum-well-like structures with targeted luminescence or exciton separation/quenching. RIR-MAPLE-deposited films enable confirmation of these predictions by optical measurements, which further display excited state behavior transcending traditional quantum-well models-i.e., with appropriate selection of specially synthesized organic/inorganic moieties, intercomponent carrier transfer efficiently converts excitons from singlet to triplet states in organics for which intersystem crossing cannot ordinarily compete with recombination. These observations demonstrate the uniquely versatile excited-state behavior in hybrid perovskite quantum wells, and the value of integrating DFT, organic synthesis, RIR-MAPLE and spectroscopy for screening/preparing rationally devised complex structures

MADNESS: A Multiresolution, Adaptive Numerical Environment for Scientific Simulation

MADNESS (multiresolution adaptive numerical environment for scientific simulation) is a high-level software environment for solving integral and differential equations in many dimensions that uses adaptive and fast harmonic analysis methods with guaranteed precision based on multiresolution analysis and separated representations. Underpinning the numerical capabilities is a powerful petascale parallel programming environment that aims to increase both programmer productivity and code scalability. This paper describes the features and capabilities of MADNESS and briefly discusses some current applications in chemistry and several areas of physics

Report on the sixth blind test of organic crystal structure prediction methods.

The sixth blind test of organic crystal structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal and a bulky flexible molecule. This blind test has seen substantial growth in the number of participants, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and `best practices' for performing CSP calculations. All of the targets, apart from a single potentially disordered Z' = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms.The organisers and participants are very grateful to the crystallographers who supplied the candidate structures: Dr. Peter Horton (XXII), Dr. Brian Samas (XXIII), Prof. Bruce Foxman (XXIV), and Prof. Kraig Wheeler (XXV and XXVI). We are also grateful to Dr. Emma Sharp and colleagues at Johnson Matthey (Pharmorphix) for the polymorph screening of XXVI, as well as numerous colleagues at the CCDC for assistance in organising the blind test. Submission 2: We acknowledge Dr. Oliver Korb for numerous useful discussions. Submission 3: The Day group acknowledge the use of the IRIDIS High Performance Computing Facility, and associated support services at the University of Southampton, in the completion of this work. We acknowledge funding from the EPSRC (grants EP/J01110X/1 and EP/K018132/1) and the European Research Council under the European Union’s Seventh Framework Programme (FP/2007-2013)/ERC through grant agreements n. 307358 (ERC-stG- 2012-ANGLE) and n. 321156 (ERC-AG-PE5-ROBOT). Submission 4: I am grateful to Mikhail Kuzminskii for calculations of molecular structures on Gaussian 98 program in the Institute of Organic Chemistry RAS. The Russian Foundation for Basic Research is acknowledged for financial support (14-03-01091). Submission 5: Toine Schreurs provided computer facilities and assistance. I am grateful to Matthew Habgood at AWE company for providing a travel grant. Submission 6: We would like to acknowledge support of this work by GlaxoSmithKline, Merck, and Vertex. Submission 7: The research was financially supported by the VIDI Research Program 700.10.427, which is financed by The Netherlands Organisation for Scientific Research (NWO), and the European Research Council (ERC-2010-StG, grant agreement n. 259510-KISMOL). We acknowledge the support of the Foundation for Fundamental Research on Matter (FOM). Supercomputer facilities were provided by the National Computing Facilities Foundation (NCF). Submission 8: Computer resources were provided by the Center for High Performance Computing at the University of Utah and the Extreme Science and Engineering Discovery Environment (XSEDE), supported by NSF grant number ACI-1053575. MBF and GIP acknowledge the support from the University of Buenos Aires and the Argentinian Research Council. Submission 9: We thank Dr. Bouke van Eijck for his valuable advice on our predicted structure of XXV. We thank the promotion office for TUT programs on advanced simulation engineering (ADSIM), the leading program for training brain information architects (BRAIN), and the information and media center (IMC) at Toyohashi University of Technology for the use of the TUT supercomputer systems and application software. We also thank the ACCMS at Kyoto University for the use of their supercomputer. In addition, we wish to thank financial supports from Conflex Corp. and Ministry of Education, Culture, Sports, Science and Technology. Submission 12: We thank Leslie Leiserowitz from the Weizmann Institute of Science and Geoffrey Hutchinson from the University of Pittsburgh for helpful discussions. We thank Adam Scovel at the Argonne Leadership Computing Facility (ALCF) for technical support. Work at Tulane University was funded by the Louisiana Board of Regents Award # LEQSF(2014-17)-RD-A-10 “Toward Crystal Engineering from First Principles”, by the NSF award # EPS-1003897 “The Louisiana Alliance for Simulation-Guided Materials Applications (LA-SiGMA)”, and by the Tulane Committee on Research Summer Fellowship. Work at the Technical University of Munich was supported by the Solar Technologies Go Hybrid initiative of the State of Bavaria, Germany. Computer time was provided by the Argonne Leadership Computing Facility (ALCF), which is supported by the Office of Science of the U.S. Department of Energy under contract DE-AC02-06CH11357. Submission 13: This work would not have been possible without funding from Khalifa University’s College of Engineering. I would like to acknowledge Prof. Robert Bennell and Prof. Bayan Sharif for supporting me in acquiring the resources needed to carry out this research. Dr. Louise Price is thanked for her guidance on the use of DMACRYS and NEIGHCRYS during the course of this research. She is also thanked for useful discussions and numerous e-mail exchanges concerning the blind test. Prof. Sarah Price is acknowledged for her support and guidance over many years and for providing access to DMACRYS and NEIGHCRYS. Submission 15: The work was supported by the United Kingdom’s Engineering and Physical Sciences Research Council (EPSRC) (EP/J003840/1, EP/J014958/1) and was made possible through access to computational resources and support from the High Performance Computing Cluster at Imperial College London. We are grateful to Professor Sarah L. Price for supplying the DMACRYS code for use within CrystalOptimizer, and to her and her research group for support with DMACRYS and feedback on CrystalPredictor and CrystalOptimizer. Submission 16: R. J. N. acknowledges financial support from the Engineering and Physical Sciences Research Council (EPSRC) of the U.K. [EP/J017639/1]. R. J. N. and C. J. P. acknowledge use of the Archer facilities of the U.K.’s national high-performance computing service (for which access was obtained via the UKCP consortium [EP/K014560/1]). C. J. P. also acknowledges a Leadership Fellowship Grant [EP/K013688/1]. B. M. acknowledges Robinson College, Cambridge, and the Cambridge Philosophical Society for a Henslow Research Fellowship. Submission 17: The work at the University of Delaware was supported by the Army Research Office under Grant W911NF-13-1- 0387 and by the National Science Foundation Grant CHE-1152899. The work at the University of Silesia was supported by the Polish National Science Centre Grant No. DEC-2012/05/B/ST4/00086. Submission 18: We would like to thank Constantinos Pantelides, Claire Adjiman and Isaac Sugden of Imperial College for their support of our use of CrystalPredictor and CrystalOptimizer in this and Submission 19. The CSP work of the group is supported by EPSRC, though grant ESPRC EP/K039229/1, and Eli Lilly. The PhD students support: RKH by a joint UCL Max-Planck Society Magdeburg Impact studentship, REW by a UCL Impact studentship; LI by the Cambridge Crystallographic Data Centre and the M3S Centre for Doctoral Training (EPSRC EP/G036675/1). Submission 19: The potential generation work at the University of Delaware was supported by the Army Research Office under Grant W911NF-13-1-0387 and by the National Science Foundation Grant CHE-1152899. Submission 20: The work at New York University was supported, in part, by the U.S. Army Research Laboratory and the U.S. Army Research Office under contract/grant number W911NF-13-1-0387 (MET and LV) and, in part, by the Materials Research Science and Engineering Center (MRSEC) program of the National Science Foundation under Award Number DMR-1420073 (MET and ES). The work at the University of Delaware was supported by the U.S. Army Research Laboratory and the U.S. Army Research Office under contract/grant number W911NF-13-1- 0387 and by the National Science Foundation Grant CHE-1152899. Submission 21: We thank the National Science Foundation (DMR-1231586), the Government of Russian Federation (Grant No. 14.A12.31.0003), the Foreign Talents Introduction and Academic Exchange Program (No. B08040) and the Russian Science Foundation, project no. 14-43-00052, base organization Photochemistry Center of the Russian Academy of Sciences. Calculations were performed on the Rurik supercomputer at Moscow Institute of Physics and Technology. Submission 22: The computational results presented have been achieved in part using the Vienna Scientific Cluster (VSC). Submission 24: The potential generation work at the University of Delaware was supported by the Army Research Office under Grant W911NF-13-1-0387 and by the National Science Foundation Grant CHE-1152899. Submission 25: J.H. and A.T. acknowledge the support from the Deutsche Forschungsgemeinschaft under the program DFG-SPP 1807. H-Y.K., R.A.D., and R.C. acknowledge support from the Department of Energy (DOE) under Grant Nos. DE-SC0008626. This research used resources of the Argonne Leadership Computing Facility at Argonne National Laboratory, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-06CH11357. This research used resources of the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DEAC02-05CH11231. Additional computational resources were provided by the Terascale Infrastructure for Groundbreaking Research in Science and Engineering (TIGRESS) High Performance Computing Center and Visualization Laboratory at Princeton University.This is the final version of the article. It first appeared from Wiley via http://dx.doi.org/10.1107/S2052520616007447

Report on the sixth blind test of organic crystal-structure prediction methods

The sixth blind test of organic crystal-structure prediction (CSP) methods has been held, with five target systems: a small nearly rigid molecule, a polymorphic former drug candidate, a chloride salt hydrate, a co-crystal, and a bulky flexible molecule. This blind test has seen substantial growth in the number of submissions, with the broad range of prediction methods giving a unique insight into the state of the art in the field. Significant progress has been seen in treating flexible molecules, usage of hierarchical approaches to ranking structures, the application of density-functional approximations, and the establishment of new workflows and "best practices" for performing CSP calculations. All of the targets, apart from a single potentially disordered Z` = 2 polymorph of the drug candidate, were predicted by at least one submission. Despite many remaining challenges, it is clear that CSP methods are becoming more applicable to a wider range of real systems, including salts, hydrates and larger flexible molecules. The results also highlight the potential for CSP calculations to complement and augment experimental studies of organic solid forms

Density-functional approaches to noncovalent interactions: A comparison of dispersion corrections (DFT-D), exchange-hole dipole moment (XDM) theory, and specialized functionals

© 2011 American Institute of Physics. The electronic version of this article is the complete one and can be found at: http://dx.doi.org/10.1063/1.3545971DOI: 10.1063/1.3545971A systematic study of techniques for treating noncovalent interactions within the computationally efficient density functional theory (DFT) framework is presented through comparison to benchmarkquality evaluations of binding strength compiled for molecular complexes of diverse size and nature. In particular, the efficacy of functionals deliberately crafted to encompass long-range forces, a posteriori DFT+dispersion corrections (DFT-D2 and DFT-D3), and exchange-hole dipole moment (XDM) theory is assessed against a large collection (469 energy points) of reference interaction energies at the CCSD(T) level of theory extrapolated to the estimated complete basis set limit. The established S22 [revised in J. Chem. Phys. 132, 144104 (2010)] and JSCH test sets of minimum-energy structures, as well as collections of dispersion-bound (NBC10) and hydrogenbonded (HBC6) dissociation curves and a pairwise decomposition of a protein–ligand reaction site (HSG), comprise the chemical systems for this work. From evaluations of accuracy, consistency, and efficiency for PBE-D, BP86-D, B97-D, PBE0-D, B3LYP-D, B970-D, M05-2X,M06-2X, ωB97X-D, B2PLYP-D, XYG3, and B3LYP-XDM methodologies, it is concluded that distinct, often contrasting, groups of these elicit the best performance within the accessible double-ζ or robust triple-ζ basis set regimes and among hydrogen-bonded or dispersion-dominated complexes. For overall results, M05-2X, B97-D3, and B970-D2 yield superior values in conjunction with aug-cc-pVDZ, for a mean absolute deviation of 0.41 – 0.49 kcal/mol, and B3LYP-D3, B97-D3, ωB97X-D, and B2PLYP-D3 dominate with aug-cc-pVTZ, affording, together with XYG3/6-311+G(3df,2p), a mean absolute deviation of 0.33 – 0.38 kcal/mol

A Combined Machine Learning and High-Energy X-ray Diffraction Approach to Understanding Liquid and Amorphous Metal Oxides

Determining the structure-property relations of liquid and amorphous metal oxides is challenging, due to their variable short-range order and polyhedral connectivity. To predict chemically realistic structures, we have developed a Machine Learned, Gaussian Approximation Potential (GAP) for HfO2, with a focus on enhanced sampling of the training database and accurate density functional theory calculations. By using training datasets for the GAP model at the level of Density Functional Theory-Strongly Constrained and Appropriately Normed (DFT-SCAN) level of theory, our results show that the topology of both the low viscosity liquid and the amorphous form are dominated by edge-shared chains and small corner-shared rings of polyhedra. This topology is shown to be consistent with the structure of other liquid and amorphous transition metal oxides of variable ion size, such as TiO2 and ZrO2. Current limitations of the ML-GAP modeling method for obtaining glass structures and future perspectives are also discussed

Quantum Mechanics Enables "Freedom of Design" in Molecular Property Space

Rational design of molecules with targeted properties requires understanding quantum-mechanical (QM) structure-property/property-property relationships (SPR/PPR) across chemical compound space. We analyze these relationships using the QM7-X dataset---which includes multiple QM properties for ~4.2 M equilibrium and non-equilibrium structures of small (primarily organic) molecules. Instead of providing simple SPR/PPR that strictly follow physicochemical intuition, our analysis uncovers substantial flexibility in molecular property space (MPS) when searching for a single molecule with a desired pair of QM properties or distinct molecules with a targeted set of QM properties. As proof-of-concept, we used Pareto multi-property optimization to search for the most promising (i.e., highly polarizable and electrically stable) molecules for polymeric battery materials; without prior knowledge of this complex manifold of MPS, Pareto front analysis reflected this intrinsic flexibility and identified small directed structural/compositional changes that simultaneously optimize these properties. Our analysis of such extensive QM property data provides compelling evidence for an intrinsic “freedom of design” in MPS, and indicates that rational design of molecules with a diverse array of targeted QM properties is quite feasible