1
|
Iyengar SS, Zhang JH, Saha D, Ricard TC. Graph-| Q⟩⟨ C|: A Quantum Algorithm with Reduced Quantum Circuit Depth for Electronic Structure. J Phys Chem A 2023; 127:9334-9345. [PMID: 37906738 DOI: 10.1021/acs.jpca.3c04261] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2023]
Abstract
The accurate determination of chemical properties is known to have a critical impact on multiple fundamental chemical problems but is deeply hindered by the steep algebraic scaling of electron correlation calculations and the exponential scaling of quantum nuclear dynamics. With the advent of new quantum computing hardware and associated developments in creating new paradigms for quantum software, this avenue has been recognized as perhaps one way to address exponentially complex challenges in quantum chemistry and molecular dynamics. In this paper, we discuss a new approach to drastically reduce the quantum circuit depth (by several orders of magnitude) and help improve the accuracy in the quantum computation of electron correlation energies for large molecular systems. The method is derived from a graph-theoretic approach to molecular fragmentation and enables us to create a family of projection operators that decompose quantum circuits into separate unitary processes. Some of these processes can be treated on quantum hardware and others on classical hardware in a completely asynchronous and parallel fashion. Numerical benchmarks are provided through the computation of unitary coupled-cluster singles and doubles (UCCSD) energies for medium-sized protonated and neutral water clusters using the new quantum algorithms presented here.
Collapse
Affiliation(s)
- Srinivasan S Iyengar
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Juncheng Harry Zhang
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Debadrita Saha
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Timothy C Ricard
- Department of Chemistry, Department of Physics, and the Indiana University Quantum Science and Engineering Center (IU-QSEC), Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
2
|
Zhang JH, Iyengar SS. Graph-| Q⟩⟨ C|, a Graph-Based Quantum/Classical Algorithm for Efficient Electronic Structure on Hybrid Quantum/Classical Hardware Systems: Improved Quantum Circuit Depth Performance. J Chem Theory Comput 2022; 18:2885-2899. [PMID: 35412836 DOI: 10.1021/acs.jctc.1c01303] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a procedure to reduce the depth of quantum circuits and improve the accuracy of results in computing post-Hartree-Fock electronic structure energies in large molecular systems. The method is based on molecular fragmentation where a molecular system is divided into overlapping fragments through a graph-theoretic procedure. This allows us to create a set of projection operators that decompose the unitary evolution of the full system into separate sets of processes, some of which can be treated on quantum hardware and others on classical hardware. Thus, we develop a procedure for an electronic structure that can be asynchronously spawned onto a potentially large ensemble of classical and quantum hardware systems. We demonstrate this method by computing Unitary Coupled Cluster Singles and Doubles (UCCSD) energies for a set of [H2]n clusters, with n ranging from 4 to 128. We implement our methodology using quantum circuits, and when these quantum circuits are processed on a quantum simulator, we obtain energies in agreement with the UCCSD energies in the milli-hartree energy range. We also show that our circuit decomposition approach yields up to 9 orders of magnitude reduction in the number of CNOT gates and quantum circuit depth for the large-sized clusters when compared to a standard quantum circuit implementation available on IBM's Quantum Information Science kit, known as Qiskit.
Collapse
Affiliation(s)
- Juncheng Harry Zhang
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
3
|
Han Y, Luo H, Lu Q, Liu Z, Liu J, Zhang J, Wei Z, Li J. Quantum Mechanical-Based Stability Evaluation of Crystal Structures for HIV-Targeted Drug Cabotegravir. Molecules 2021; 26:molecules26237178. [PMID: 34885762 PMCID: PMC8659202 DOI: 10.3390/molecules26237178] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Revised: 11/21/2021] [Accepted: 11/23/2021] [Indexed: 11/23/2022] Open
Abstract
The long-acting parenteral formulation of the HIV integrase inhibitor cabotegravir (GSK744) is currently being developed to prevent HIV infections, benefiting from infrequent dosing and high efficacy. The crystal structure can affect the bioavailability and efficacy of cabotegravir. However, the stability determination of crystal structures of GSK744 have remained a challenge. Here, we introduced an ab initio protocol to determine the stability of the crystal structures of pharmaceutical molecules, which were obtained from crystal structure prediction process starting from the molecular diagram. Using GSK744 as a case study, the ab initio predicted that Gibbs free energy provides reliable further refinement of the predicted crystal structures and presents its capability for becoming a crystal stability determination approach in the future. The proposed work can assist in the comprehensive screening of pharmaceutical design and can provide structural predictions and stability evaluation for pharmaceutical crystals.
Collapse
Affiliation(s)
- Yanqiang Han
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
| | - Hongyuan Luo
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
| | - Qianqian Lu
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
| | - Zeying Liu
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
| | - Jinyun Liu
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Provincial Engineering Laboratory for New-Energy Vehicle Battery Energy-Storage Materials, School of Chemistry and Materials Science, Anhui Normal University, Wuhu 241002, China
- Correspondence: (J.L.); (Z.W.); (J.L.)
| | - Jiarui Zhang
- Division of Computational Biomedicine, Boston University School of Medicine, Boston, MA 02118, USA;
| | - Zhiyun Wei
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
- Correspondence: (J.L.); (Z.W.); (J.L.)
| | - Jinjin Li
- Shanghai Key Laboratory of Maternal Fetal Medicine, Shanghai First Maternity and Infant Hospital, School of Medicine, Tongji University, Shanghai 200092, China; (Y.H.); (Z.L.)
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai 200240, China; (H.L.); (Q.L.)
- Correspondence: (J.L.); (Z.W.); (J.L.)
| |
Collapse
|
4
|
Tang J, Han Y, Ali I, Luo H, Nowak A, Li J. Stability and phase transition investigation of olanzapine polymorphs. Chem Phys Lett 2021. [DOI: 10.1016/j.cplett.2021.138384] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2022]
|
5
|
Han Y, Wang Z, Li J. Neural Networks Accelerate the Ab Initio Prediction of Solid-Solid Phase Transitions at High Pressures. J Phys Chem Lett 2021; 12:132-137. [PMID: 33314933 DOI: 10.1021/acs.jpclett.0c03101] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/12/2023]
Abstract
High-level ab initio chemical calculations, such as second-order Møller-Plesset perturbation (MP2), are highly accurate but time-consuming, making it inefficient to apply to macromolecular systems. Here, we propose a newly efficient approach based on the neural network and fragment method to predict the Gibbs free energy, structural characteristics, and thus phase transition of solid crystal structures. The proposed approach has the same prediction accuracy as the MP2 calculation but is hundreds of times faster than the MP2. The predicted structures and phase transitions of two selected ice phases (IX and XV) under extreme conditions are in excellent agreement with the MP2 calculations and experimental results but with an extremely low computational cost. It not only predicts the high-pressure structures and phase diagrams of solid systems accurately and efficiently but also solves the problem of extreme calculation cost during a high-precision theoretical study on high-pressure molecular crystals with potentially essential applications.
Collapse
|
6
|
Qin X, Hirata S. Anharmonic Phonon Dispersion in Polyethylene. J Phys Chem B 2020; 124:10477-10485. [PMID: 33169996 DOI: 10.1021/acs.jpcb.0c08493] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The second-order Green's function method for anharmonic crystals has been applied to an infinite, periodic chain of polyethylene taking into account up to quartic force constants. The frequency-independent approximation to the Dyson self-energy gives rise to numerous divergent resonances, which are fortuitous. Instead, solving the Dyson equation self-consistently with a frequency-dependent self-energy resists divergences from resonances or zero-frequency acoustic vibrations. The calculated anharmonic phonon dispersion, which nonetheless displays many true resonances, and anharmonic phonon density of states furnish hitherto unknown details that explain smaller features of observed vibrational spectra.
Collapse
Affiliation(s)
- Xiuyi Qin
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, Urbana, Illinois 61801, United States
| |
Collapse
|
7
|
Huang L, Han Y, Liu J, He X, Li J. Ab Initio Prediction of the Phase Transition for Solid Ammonia at High Pressures. Sci Rep 2020; 10:7546. [PMID: 32372007 PMCID: PMC7200730 DOI: 10.1038/s41598-020-64030-3] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Accepted: 04/08/2020] [Indexed: 11/11/2022] Open
Abstract
Ammonia is one of the most basic components on the planet and its high-pressure characteristics play an important role in planetary science. Solid ammonia crystals frequently adopt multiple distinct polymorphs exhibiting different properties. Predicting the crystal structure of these polymorphs and under what thermodynamic conditions these polymorphs are stable would be of great value to environmental industry and other fields. Theoretical calculations based on the classical force fields and density-functional theory (DFT) are versatile methods but lack of accurate description of weak intermolecular interactions for molecular crystals. In this study, we employ an ab initio computational study on the solid ammonia at high pressures, using the second-order Møller-Plesset perturbation (MP2) theory and the coupled cluster singles, doubles, and perturbative triples (CCSD(T)) theory along with the embedded fragmentation method. The proposed algorithm is capable of performing large-scale calculations using high-level wavefunction theories, and accurately describing covalent, ionic, hydrogen bonding, and dispersion interactions within molecular crystals, and therefore can predict the crystal structures, Raman spectra and phase transition of solid ammonia phases I and IV accurately. We confirm the crystal structures of solid ammonia phases I and IV that have been controversial for a long time and predict their phase transition that occurs at 1.17 GPa and 210 K with small temperature dependence, which is in line with experiment.
Collapse
Affiliation(s)
- Lei Huang
- Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China
| | - Yanqiang Han
- Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai, 200240, China
| | - Jinyun Liu
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu, Anhui, 241000, China.
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development, School of Chemistry and Molecular Engineering, East China Normal University, Shanghai, 200062, China.
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, Shanghai, 200062, China.
| | - Jinjin Li
- Key laboratory for Thin Film and Microfabrication of the Ministry of Education, Department of Micro/Nano-electronics, Shanghai Jiao Tong University, Shanghai, 200240, China.
- Key Laboratory of Functional Molecular Solids of the Ministry of Education, Anhui Laboratory of Molecule-Based Materials, College of Chemistry and Materials Science, Anhui Normal University, Wuhu, Anhui, 241000, China.
| |
Collapse
|
8
|
Huang L, Han Y, He X, Li J. Ab initio-enabled phase transition prediction of solid carbon dioxide at ultra-high temperatures. RSC Adv 2020; 10:236-243. [PMID: 35492555 PMCID: PMC9049158 DOI: 10.1039/c9ra06478h] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2019] [Accepted: 11/12/2019] [Indexed: 12/27/2022] Open
Abstract
Carbon dioxide is one of the fundamental chemical species on Earth, but its solid-phase behavior at high pressures is still far from well understood and some phases remain uncertain or unknown, which increases the challenge to predict its structures.
Collapse
Affiliation(s)
- Lei Huang
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education
- Department of Micro/Nano-Electronics
- Shanghai Jiao Tong University
- Shanghai 200240
- China
| | - Yanqiang Han
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education
- Department of Micro/Nano-Electronics
- Shanghai Jiao Tong University
- Shanghai 200240
- China
| | - Xiao He
- Shanghai Engineering Research Center of Molecular Therapeutics and New Drug Development
- School of Chemistry and Molecular Engineering
- East China Normal University
- Shanghai
- China
| | - Jinjin Li
- Key Laboratory for Thin Film and Microfabrication of Ministry of Education
- Department of Micro/Nano-Electronics
- Shanghai Jiao Tong University
- Shanghai 200240
- China
| |
Collapse
|
9
|
Crystal Structure Optimization and Gibbs Free Energy Comparison of Five Sulfathiazole Polymorphs by the Embedded Fragment QM Method at the DFT Level. CRYSTALS 2019. [DOI: 10.3390/cryst9050256] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
Abstract
Molecular crystal plays an important role in many fields of science and technology, but it often crystallizes in different polymorphs with different physical properties. To guide the experimental synthesis of candidate materials, the atomic-scale model is frequently used to predict the most stable polymorph and its structural properties. Here, we show how an ab initio method can be used to achieve a rapid and accurate prediction of sulfathiazole crystal polymorphs (an antibiotic drug), based on the Gibbs free energy calculation and Raman spectra analysis. At the atmospheric pressure and the temperature of 300 K, we demonstrate that form III (FIII) is the most stable structure of sulfathiazole. The agreement between the predicted and experimental crystal structures corresponds to the order of stability for five sulfathiazole polymorphs as FI < FV < FIV < FII < FIII, which is achieved by employing the density functional theory (DFT) calculations.
Collapse
|
10
|
Ricard TC, Haycraft C, Iyengar SS. Adaptive, Geometric Networks for Efficient Coarse-Grained Ab Initio Molecular Dynamics with Post-Hartree–Fock Accuracy. J Chem Theory Comput 2018; 14:2852-2866. [DOI: 10.1021/acs.jctc.8b00186] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Affiliation(s)
- Timothy C. Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Cody Haycraft
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Srinivasan S. Iyengar
- Department of Chemistry and Department of Physics, Indiana University, 800 East Kirkwood Avenue, Bloomington, Indiana 47405, United States
| |
Collapse
|
11
|
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
|
12
|
Salim MA, Willow SY, Hirata S. Ice Ih anomalies: Thermal contraction, anomalous volume isotope effect, and pressure-induced amorphization. J Chem Phys 2016; 144:204503. [DOI: 10.1063/1.4951687] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Michael A. Salim
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - Soohaeng Yoo Willow
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| | - So Hirata
- Department of Chemistry, University of Illinois at Urbana-Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, USA
| |
Collapse
|
13
|
Fang T, Jia J, Li S. Vibrational Spectra of Molecular Crystals with the Generalized Energy-Based Fragmentation Approach. J Phys Chem A 2016; 120:2700-11. [DOI: 10.1021/acs.jpca.5b10927] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/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 Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Junteng Jia
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| | - Shuhua Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210093, People’s Republic of China
| |
Collapse
|
14
|
Abstract
Interest in molecular crystals has grown thanks to their relevance to pharmaceuticals, organic semiconductor materials, foods, and many other applications. Electronic structure methods have become an increasingly important tool for modeling molecular crystals and polymorphism. This article reviews electronic structure techniques used to model molecular crystals, including periodic density functional theory, periodic second-order Møller-Plesset perturbation theory, fragment-based electronic structure methods, and diffusion Monte Carlo. It also discusses the use of these models for predicting a variety of crystal properties that are relevant to the study of polymorphism, including lattice energies, structures, crystal structure prediction, polymorphism, phase diagrams, vibrational spectroscopies, and nuclear magnetic resonance spectroscopy. Finally, tools for analyzing crystal structures and intermolecular interactions are briefly discussed.
Collapse
Affiliation(s)
- Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
15
|
[FHF]−—The Strongest Hydrogen Bond under the Influence of External Interactions. CRYSTALS 2015. [DOI: 10.3390/cryst6010003] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
16
|
Hartman JD, Beran GJO. Fragment-Based Electronic Structure Approach for Computing Nuclear Magnetic Resonance Chemical Shifts in Molecular Crystals. J Chem Theory Comput 2015; 10:4862-72. [PMID: 26584373 DOI: 10.1021/ct500749h] [Citation(s) in RCA: 43] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Abstract
First-principles chemical shielding tensor predictions play a critical role in studying molecular crystal structures using nuclear magnetic resonance. Fragment-based electronic structure methods have dramatically improved the ability to model molecular crystal structures and energetics using high-level electronic structure methods. Here, a many-body expansion fragment approach is applied to the calculation of chemical shielding tensors in molecular crystals. First, the impact of truncating the many-body expansion at different orders and the role of electrostatic embedding are examined on a series of molecular clusters extracted from molecular crystals. Second, the ability of these techniques to assign three polymorphic forms of the drug sulfanilamide to the corresponding experimental (13)C spectra is assessed. This challenging example requires discriminating among spectra whose (13)C chemical shifts differ by only a few parts per million (ppm) across the different polymorphs. Fragment-based PBE0/6-311+G(2d,p) level chemical shielding predictions correctly assign these three polymorphs and reproduce the sulfanilamide experimental (13)C chemical shifts with 1 ppm accuracy. The results demonstrate that fragment approaches are competitive with the widely used gauge-invariant projector augmented wave (GIPAW) periodic density functional theory calculations.
Collapse
Affiliation(s)
- Joshua D Hartman
- Department of Chemistry, University of California , Riverside, California 92521, United States
| | - Gregory J O Beran
- Department of Chemistry, University of California , Riverside, California 92521, United States
| |
Collapse
|
17
|
Byrd JN, Jindal N, Molt RW, Bartlett RJ, Sanders BA, Lotrich VF. Molecular cluster perturbation theory. I. Formalism. Mol Phys 2015. [DOI: 10.1080/00268976.2015.1036145] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
18
|
Li J, Sode O, Hirata S. Second-Order Many-Body Perturbation Study on Thermal Expansion of Solid Carbon Dioxide. J Chem Theory Comput 2014; 11:224-9. [DOI: 10.1021/ct500983k] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022]
Affiliation(s)
- Jinjin Li
- Department
of Chemistry, University of Illinois at Urbana-Champaign, 600
South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Olaseni Sode
- Department
of Chemistry, University of Chicago, 5735 South Ellis Avenue, Chicago, Illinois 60637, United States
| | - So Hirata
- Department
of Chemistry, University of Illinois at Urbana-Champaign, 600
South Mathews Avenue, Urbana, Illinois 61801, United States
- CREST, Japan Science and Technology Agency, 4-1-8 Honcho, Kawaguchi, Saitama 332-0012, Japan
| |
Collapse
|
19
|
Fang T, Li W, Gu F, Li S. Accurate Prediction of Lattice Energies and Structures of Molecular Crystals with Molecular Quantum Chemistry Methods. J Chem Theory Comput 2014; 11:91-8. [DOI: 10.1021/ct500833k] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Tao Fang
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| | - Wei Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| | - Fangwei Gu
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| | - Shuhua Li
- School of Chemistry and Chemical
Engineering, Key Laboratory of Mesoscopic Chemistry of Ministry of
Education, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing, 210093, People’s Republic of China
| |
Collapse
|
20
|
Heit Y, Beran GJO. Exploiting space-group symmetry in fragment-based molecular crystal calculations. J Comput Chem 2014; 35:2205-14. [DOI: 10.1002/jcc.23737] [Citation(s) in RCA: 17] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2014] [Revised: 07/23/2014] [Accepted: 08/31/2014] [Indexed: 02/03/2023]
Affiliation(s)
- Yonaton Heit
- Department of Chemistry; University of California; Riverside California 92521
| | - Gregory J. O. Beran
- Department of Chemistry; University of California; Riverside California 92521
| |
Collapse
|
21
|
Hirata S, Gilliard K, He X, Li J, Sode O. Ab initio molecular crystal structures, spectra, and phase diagrams. Acc Chem Res 2014; 47:2721-30. [PMID: 24754304 DOI: 10.1021/ar500041m] [Citation(s) in RCA: 68] [Impact Index Per Article: 6.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Conspectus Molecular crystals are chemists' solids in the sense that their structures and properties can be understood in terms of those of the constituent molecules merely perturbed by a crystalline environment. They form a large and important class of solids including ices of atmospheric species, drugs, explosives, and even some organic optoelectronic materials and supramolecular assemblies. Recently, surprisingly simple yet extremely efficient, versatile, easily implemented, and systematically accurate electronic structure methods for molecular crystals have been developed. The methods, collectively referred to as the embedded-fragment scheme, divide a crystal into monomers and overlapping dimers and apply modern molecular electronic structure methods and software to these fragments of the crystal that are embedded in a self-consistently determined crystalline electrostatic field. They enable facile applications of accurate but otherwise prohibitively expensive ab initio molecular orbital theories such as Møller-Plesset perturbation and coupled-cluster theories to a broad range of properties of solids such as internal energies, enthalpies, structures, equation of state, phonon dispersion curves and density of states, infrared and Raman spectra (including band intensities and sometimes anharmonic effects), inelastic neutron scattering spectra, heat capacities, Gibbs energies, and phase diagrams, while accounting for many-body electrostatic (namely, induction or polarization) effects as well as two-body exchange and dispersion interactions from first principles. They can fundamentally alter the role of computing in the studies of molecular crystals in the same way ab initio molecular orbital theories have transformed research practices in gas-phase physical chemistry and synthetic chemistry in the last half century. In this Account, after a brief summary of formalisms and algorithms, we discuss applications of these methods performed in our group as compelling illustrations of their unprecedented power in addressing some of the outstanding problems of solid-state chemistry, high-pressure chemistry, or geochemistry. They are the structure and spectra of ice Ih, in particular, the origin of two peaks in the hydrogen-bond-stretching region of its inelastic neutron scattering spectra, a solid-solid phase transition from CO2-I to elusive, metastable CO2-III, pressure tuning of Fermi resonance in solid CO2, and the structure and spectra of solid formic acid, all at the level of second-order Møller-Plesset perturbation theory or higher.
Collapse
Affiliation(s)
- So Hirata
- Department of Chemistry, University of Illinois at Urbana−Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Kandis Gilliard
- Department of Chemistry, University of Illinois at Urbana−Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Xiao He
- Department of Chemistry, University of Illinois at Urbana−Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Jinjin Li
- Department of Chemistry, University of Illinois at Urbana−Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| | - Olaseni Sode
- Department of Chemistry, University of Illinois at Urbana−Champaign, 600 South Mathews Avenue, Urbana, Illinois 61801, United States
| |
Collapse
|
22
|
Gilliard K, Sode O, Hirata S. Second-order many-body perturbation and coupled-cluster singles and doubles study of ice VIII. J Chem Phys 2014; 140:174507. [DOI: 10.1063/1.4873919] [Citation(s) in RCA: 19] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
23
|
Li J, Sode O, Voth GA, Hirata S. A solid–solid phase transition in carbon dioxide at high pressures and intermediate temperatures. Nat Commun 2013; 4:2647. [DOI: 10.1038/ncomms3647] [Citation(s) in RCA: 47] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/23/2013] [Accepted: 09/19/2013] [Indexed: 11/09/2022] Open
|
24
|
Liu H, Wang Y, Bowman JM. Vibrational analysis of an ice Ih model from 0 to 4000 cm(-1) using the ab initio WHBB potential energy surface. J Phys Chem B 2013; 117:10046-52. [PMID: 23924359 DOI: 10.1021/jp405865c] [Citation(s) in RCA: 18] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/27/2022]
Abstract
We present an analysis of the vibrational modes of a model of hexagonal ice, ice Ih, comprised of 192 monomers with a core region of 105 monomers, using the ab initio WHBB potential energy surface [Wang, Y.; Shepler, B.; Braams, B.; Bowman, J. M. J. Chem. Phys. 2011, 134, 094509]. A standard normal-mode analysis and a local-monomer normal-mode analysis of 105 core monomers are performed to obtain harmonic frequencies and state densities of the "pseudo-translation" (0-400 cm(-1)), "libration" (500-1100 cm(-1)), monomer bend fundamental (~1600 cm), and O-H stretch (~3000-3700 cm(-1)) bands. In addition, the coupled local-monomer model is used to obtain the vibrational density of states in the bend fundamental and O-H stretch regions. The harmonic and local-monomer vibrational density of states obtained from core monomers are in good agreement with those of inelastic neutron scattering spectra, especially the latter, which accounts for anharmonic coupling of monomer modes. Full deuteration is also considered, and the vibrational density of states is again compared to experiment, where good agreement is found.
Collapse
Affiliation(s)
- Hanchao Liu
- Cherry L. Emerson Center for Scientific Computation and Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | | | | |
Collapse
|
25
|
Sode O, Keçeli M, Yagi K, Hirata S. Fermi resonance in solid CO2 under pressure. J Chem Phys 2013; 138:074501. [DOI: 10.1063/1.4790537] [Citation(s) in RCA: 45] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
26
|
He X, Sode O, Xantheas SS, Hirata S. Second-order many-body perturbation study of ice Ih. J Chem Phys 2012. [DOI: 10.1063/1.4767898] [Citation(s) in RCA: 66] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
27
|
|