1
|
Pathirage PDVS, Phillips JT, Vogiatzis KD. Exploration of the Two-Electron Excitation Space with Data-Driven Coupled Cluster. J Phys Chem A 2024. [PMID: 38422511 DOI: 10.1021/acs.jpca.3c06600] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/02/2024]
Abstract
Computational cost limits the applicability of post-Hartree-Fock methods such as coupled-cluster on larger molecular systems. The data-driven coupled-cluster (DDCC) method applies machine learning to predict the coupled-cluster two-electron amplitudes (t2) using data from second-order perturbation theory (MP2). One major limitation of the DDCC models is the size of training sets that increases exponentially with the system size. Effective sampling of the amplitude space can resolve this issue. Five different amplitude selection techniques that reduce the amount of data used for training were evaluated, an approach that also prevents model overfitting and increases the portability of data-driven coupled-cluster singles and doubles to more complex molecules or larger basis sets. In combination with a localized orbital formalism to predict the CCSD t2 amplitudes, we have achieved a 10-fold error reduction for energy calculations.
Collapse
Affiliation(s)
- P D Varuna S Pathirage
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600, United States
| | - Justin T Phillips
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600, United States
| | - Konstantinos D Vogiatzis
- Department of Chemistry, University of Tennessee, Knoxville, Tennessee 37996-1600, United States
| |
Collapse
|
2
|
Galvez Vallejo JL, Snowdon C, Stocks R, Kazemian F, Yan Yu FC, Seidl C, Seeger Z, Alkan M, Poole D, Westheimer BM, Basha M, De La Pierre M, Rendell A, Izgorodina EI, Gordon MS, Barca GMJ. Toward an extreme-scale electronic structure system. J Chem Phys 2023; 159:044112. [PMID: 37497819 DOI: 10.1063/5.0156399] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2023] [Accepted: 07/03/2023] [Indexed: 07/28/2023] Open
Abstract
Electronic structure calculations have the potential to predict key matter transformations for applications of strategic technological importance, from drug discovery to material science and catalysis. However, a predictive physicochemical characterization of these processes often requires accurate quantum chemical modeling of complex molecular systems with hundreds to thousands of atoms. Due to the computationally demanding nature of electronic structure calculations and the complexity of modern high-performance computing hardware, quantum chemistry software has historically failed to operate at such large molecular scales with accuracy and speed that are useful in practice. In this paper, novel algorithms and software are presented that enable extreme-scale quantum chemistry capabilities with particular emphasis on exascale calculations. This includes the development and application of the multi-Graphics Processing Unit (GPU) library LibCChem 2.0 as part of the General Atomic and Molecular Electronic Structure System package and of the standalone Extreme-scale Electronic Structure System (EXESS), designed from the ground up for scaling on thousands of GPUs to perform high-performance accurate quantum chemistry calculations at unprecedented speed and molecular scales. Among various results, we report that the EXESS implementation enables Hartree-Fock/cc-pVDZ plus RI-MP2/cc-pVDZ/cc-pVDZ-RIFIT calculations on an ionic liquid system with 623 016 electrons and 146 592 atoms in less than 45 min using 27 600 GPUs on the Summit supercomputer with a 94.6% parallel efficiency.
Collapse
Affiliation(s)
| | - Calum Snowdon
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Ryan Stocks
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Fazeleh Kazemian
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Fiona Chuo Yan Yu
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Christopher Seidl
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| | - Zoe Seeger
- School of Chemistry, Monash University, Clayton 3800, VIC, Australia
| | - Melisa Alkan
- Department of Chemistry, Iowa State University, Ames, Iowa 50011-3111, USA
| | - David Poole
- School of Chemistry and Biochemistry, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Bryce M Westheimer
- Department of Chemistry, Iowa State University, Ames, Iowa 50011-3111, USA
| | - Mehaboob Basha
- Pawsey Supercomputing Research Centre, Kensington, WA 6151, Australia
| | | | - Alistair Rendell
- College of Science and Engineering, Flinders University, Adelaide, SA 5042, Australia
| | | | | | - Giuseppe M J Barca
- School of Computing, Australian National University, Canberra 2601, ACT, Australia
| |
Collapse
|
3
|
Naim C, Besalú-Sala P, Zaleśny R, Luis JM, Castet F, Matito E. Are Accelerated and Enhanced Wave Function Methods Accurate to Compute Static Linear and Nonlinear Optical Properties? J Chem Theory Comput 2023; 19:1753-1764. [PMID: 36862983 DOI: 10.1021/acs.jctc.2c01212] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/04/2023]
Abstract
Key components of organic-based electro-optic devices are challenging to design or optimize because they exhibit nonlinear optical responses, which are difficult to model or rationalize. Computational chemistry furnishes the tools to investigate extensive collections of molecules in the quest for target compounds. Among the electronic structure methods that provide static nonlinear optical properties (SNLOPs), density functional approximations (DFAs) are often preferred because of their low cost/accuracy ratio. However, the accuracy of the SNLOPs critically depends on the amount of exact exchange and electron correlation included in the DFA, precluding the reliable calculation of many molecular systems. In this scenario, wave function methods such as MP2, CCSD, and CCSD(T) constitute a reliable alternative to compute SNLOPs. Unfortunately, the computational cost of these methods significantly restricts the size of molecules to study, a limitation that hampers the identification of molecules with significant nonlinear optical responses. This paper analyzes various flavors and alternatives to MP2, CCSD, and CCSD(T) methods that either drastically reduce the computational cost or improve their performance but were scarcely and unsystematically employed to compute SNLOPs. In particular, we have tested RI-MP2, RIJK-MP2, RIJCOSX-MP2 (with GridX2 and GridX4 setups), LMP2, SCS-MP2, SOS-MP2, DLPNO-MP2, LNO-CCSD, LNO-CCSD(T), DLPNO-CCSD, DLPNO-CCSD(T0), and DLPNO-CCSD(T1). Our results indicate that all these methods can be safely employed to calculate the dipole moment and the polarizability with average relative errors below 5% with respect to CCSD(T). On the other hand, the calculation of higher-order properties represents a challenge for LNO and DLPNO methods, which present severe numerical instabilities in computing the single-point field-dependent energies. RI-MP2, RIJK-MP2, or RIJCOSX-MP2 are cost-effective methods to compute first and second hyperpolarizabilities with a marginal average error with respect to canonical MP2 (up to 5% for β and up to 11% for γ). More accurate hyperpolarizabilities can be obtained with DLPNO-CCSD(T1); however, this method cannot be employed to obtain reliable second hyperpolarizabilities. These results open the way to obtain accurate nonlinear optical properties at a computational cost that can compete with current DFAs.
Collapse
Affiliation(s)
- Carmelo Naim
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, 20018 Donostia, Euskadi, Spain.,Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France.,Polimero eta Material Aurreratuak: Fisika, Kimika eta Teknologia, Kimika Fakultatea, Euskal Herriko Unibertsitatea UPV/EHU, 20080 Donostia, Euskadi, Spain
| | - Pau Besalú-Sala
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia, Spain
| | - Robert Zaleśny
- Faculty of Chemistry, Wrocław University of Science and Technology, Wyb. Wyspiańskiego 27, PL-50370 Wrocław, Poland
| | - Josep M Luis
- Institut de Química Computacional i Catàlisi and Departament de Química, Universitat de Girona, 17003 Girona, Catalonia, Spain
| | - Frédéric Castet
- Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France
| | - Eduard Matito
- Donostia International Physics Center (DIPC), Manuel Lardizabal Ibilbidea 4, 20018 Donostia, Euskadi, Spain.,Ikerbasque Foundation for Science, 48011 Bilbao, Euskadi, Spain
| |
Collapse
|
4
|
Bosia F, Zheng P, Vaucher A, Weymuth T, Dral PO, Reiher M. Ultra-fast semi-empirical quantum chemistry for high-throughput computational campaigns with Sparrow. J Chem Phys 2023; 158:054118. [PMID: 36754821 DOI: 10.1063/5.0136404] [Citation(s) in RCA: 11] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023] Open
Abstract
Semi-empirical quantum chemical approaches are known to compromise accuracy for the feasibility of calculations on huge molecules. However, the need for ultrafast calculations in interactive quantum mechanical studies, high-throughput virtual screening, and data-driven machine learning has shifted the emphasis toward calculation runtimes recently. This comes with new constraints for the software implementation as many fast calculations would suffer from a large overhead of the manual setup and other procedures that are comparatively fast when studying a single molecular structure, but which become prohibitively slow for high-throughput demands. In this work, we discuss the effect of various well-established semi-empirical approximations on calculation speed and relate this to data transfer rates from the raw-data source computer to the results of the visualization front end. For the former, we consider desktop computers, local high performance computing, and remote cloud services in order to elucidate the effect on interactive calculations, for web and cloud interfaces in local applications, and in world-wide interactive virtual sessions. The models discussed in this work have been implemented into our open-source software SCINE Sparrow.
Collapse
Affiliation(s)
- Francesco Bosia
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Peikun Zheng
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Alain Vaucher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Thomas Weymuth
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| | - Pavlo O Dral
- State Key Laboratory of Physical Chemistry of Solid Surfaces, Fujian Provincial Key Laboratory of Theoretical and Computational Chemistry, Department of Chemistry, and College of Chemistry and Chemical Engineering, Xiamen University, Xiamen 361005, China
| | - Markus Reiher
- Laboratory of Physical Chemistry, ETH Zurich, Vladimir-Prelog-Weg 2, 8093 Zurich, Switzerland
| |
Collapse
|
5
|
Hesselmann A, Werner HJ, Knowles PJ. Thermochemical evaluation of adaptive and fixed density functional theory quadrature schemes. J Chem Phys 2022; 157:234106. [PMID: 36550055 DOI: 10.1063/5.0119622] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
A systematic study is made of the accuracy and efficiency of a number of existing quadrature schemes for molecular Kohn-Sham Density-Functional Theory (DFT) using 408 molecules and 254 chemical reactions. Included are the fixed SG-x (x = 0-3) grids of Gill et al., Dasgupta, and Herbert, the 3-zone grids of Treutler and Ahlrichs, a fixed five-zone grid implemented in Molpro, and a new adaptive grid scheme. While all methods provide a systematic reduction of errors upon extension of the grid sizes, significant differences are observed in the accuracies for similar grid sizes with various approaches. For the tests in this work, the SG-x fixed grids are less suitable to achieve high accuracies in the DFT integration, while our new adaptive grid performed best among the schemes studied in this work. The extra computational time to generate the adaptive grid scales linearly with molecular size and is negligible compared with the time needed for the self-consistent field iterations for large molecules. A comparison of the grid accuracies using various density functionals shows that meta-GGA functionals need larger integration grids than GGA functionals to reach the same degree of accuracy, confirming previous investigations of the numerical stability of meta-GGA functionals. On the other hand, the grid integration errors are almost independent of the basis set, and the basis set errors are mostly much larger than the errors caused by the numerical integrations, even when using the smallest grids tested in this work.
Collapse
Affiliation(s)
- Andreas Hesselmann
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Peter J Knowles
- School of Chemistry, Cardiff University, Main Building, Park Place, Cardiff CF10 3AT, United Kingdom
| |
Collapse
|
6
|
Dutkiewicz Z. Computational methods for calculation of protein-ligand binding affinities in structure-based drug design. PHYSICAL SCIENCES REVIEWS 2022. [DOI: 10.1515/psr-2020-0034] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023]
Abstract
Abstract
Drug design is an expensive and time-consuming process. Any method that allows reducing the time the costs of the drug development project can have great practical value for the pharmaceutical industry. In structure-based drug design, affinity prediction methods are of great importance. The majority of methods used to predict binding free energy in protein-ligand complexes use molecular mechanics methods. However, many limitations of these methods in describing interactions exist. An attempt to go beyond these limits is the application of quantum-mechanical description for all or only part of the analyzed system. However, the extensive use of quantum mechanical (QM) approaches in drug discovery is still a demanding challenge. This chapter briefly reviews selected methods used to calculate protein-ligand binding affinity applied in virtual screening (VS), rescoring of docked poses, and lead optimization stage, including QM methods based on molecular simulations.
Collapse
Affiliation(s)
- Zbigniew Dutkiewicz
- Department of Chemical Technology of Drugs , Poznan University of Medical Sciences , ul. Grunwaldzka 6 , 60-780 Poznań , Poznan , 60-780, Poland
| |
Collapse
|
7
|
Sacchetta F, Graf D, Laqua H, Ambroise MA, Kussmann J, Dreuw A, Ochsenfeld C. An effective sub-quadratic scaling atomic-orbital reformulation of the scaled opposite-spin RI-CC2 ground-state model using Cholesky-decomposed densities and an attenuated Coulomb metric. J Chem Phys 2022; 157:104104. [DOI: 10.1063/5.0098719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
An atomic-orbital reformulation of the Laplace-transformed scaled opposite-spin (SOS) coupled cluster singles and doubles (CC2) model within the resolution of the identity (RI) approximation (SOS-RI-CC2) is presented that extends its applicability to molecules with several hundreds of atoms and triple-zeta basis sets. We exploit sparse linear algebra and an attenuated Coulomb metric to decrease the disk space demands and the computational efforts. In this way, an effective sub-quadratic computational scaling is achieved with our ω-SOS-CDD-RI-CC2 model. Moreover, Cholesky decomposition of the ground-state one-electron density matrix reduces the prefactor, allowing for an early crossover with the molecular orbital formulation. The accuracy and performance of the presented method are investigated for various molecular systems.
Collapse
Affiliation(s)
- F. Sacchetta
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - D. Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - H. Laqua
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - M. A. Ambroise
- Chair of Theoretical and Computational Chemistry, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - J. Kussmann
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - A. Dreuw
- Chair of Theoretical and Computational Chemistry, Interdisciplinary Center for Scientific Computing, Heidelberg University, Heidelberg, Germany
| | - C. Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Munich, Germany
| |
Collapse
|
8
|
Yang Z, Hajlasz N, Kulik HJ. Computational Modeling of Conformer Stability in Benenodin-1, a Thermally Actuated Lasso Peptide Switch. J Phys Chem B 2022; 126:3398-3406. [PMID: 35481742 DOI: 10.1021/acs.jpcb.2c00762] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Benenodin-1 is a thermally actuated lasso peptide rotaxane switch with two primary translational isomers that differ in the relative position of the residue Gln15. The conversion from one conformer to the other involves substantial enthalpy-entropy compensation: one conformer is energetically favored and the other is entropically favored. Here, we take a multi-scale quantum mechanical (QM) and classical molecular dynamic (MD) approach to reveal residue-specific sources of these differences in stability. QM reveals that the two benenodin-1 conformers involve distinct hydrogen bonding networks, with the enthalpically favored conformer having more intra-peptide hydrogen bonds between the Gln15 side chain and nearby residues. The evaluation of configurational entropy over the MD-sampled geometries reveals that the entropically favored conformer has enhanced conformational flexibility. By computing the by-residue-sum entropies, we identify the role of Gln15 and neighboring Glu14 in mediating the entropic variation during the switching process. These computational insights help explain the effects of Glu14Ala and Gln15Ala mutations on the conformational population of benenodin-1 observed experimentally.
Collapse
Affiliation(s)
- Zhongyue Yang
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Natalia Hajlasz
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
9
|
Shpiro B, Fabian M, Rabani E, Baer R. Forces from Stochastic Density Functional Theory under Nonorthogonal Atom-Centered Basis Sets. J Chem Theory Comput 2022; 18:1458-1466. [PMID: 35099187 PMCID: PMC8908760 DOI: 10.1021/acs.jctc.1c00794] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2021] [Indexed: 11/28/2022]
Abstract
We develop a formalism for calculating forces on the nuclei within the linear-scaling stochastic density functional theory (sDFT) in a nonorthogonal atom-centered basis set representation (Fabian et al. Wiley Interdiscip. Rev.: Comput. Mol. Sci. 2019, 9, e1412, 10.1002/wcms.1412) and apply it to the Tryptophan Zipper 2 (Trp-zip2) peptide solvated in water. We use an embedded-fragment approach to reduce the statistical errors (fluctuation and systematic bias), where the entire peptide is the main fragment and the remaining 425 water molecules are grouped into small fragments. We analyze the magnitude of the statistical errors in the forces and find that the systematic bias is of the order of 0.065 eV/Å (∼1.2 × 10-3Eh/a0) when 120 stochastic orbitals are used, independently of system size. This magnitude of bias is sufficiently small to ensure that the bond lengths estimated by stochastic DFT (within a Langevin molecular dynamics simulation) will deviate by less than 1% from those predicted by a deterministic calculation.
Collapse
Affiliation(s)
- Ben Shpiro
- Fritz
Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Marcel
David Fabian
- Fritz
Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| | - Eran Rabani
- Department
of Chemistry, University of California, Berkeley, California 94720, United States
- Materials
Sciences Division, Lawrence Berkeley National
Laboratory, Berkeley, California 94720, United States
- The
Raymond and Beverly Sackler Center of Computational Molecular and
Materials Science, Tel Aviv University, Tel Aviv 69978, Israel
| | - Roi Baer
- Fritz
Haber Center for Molecular Dynamics and Institute of Chemistry, The Hebrew University of Jerusalem, Jerusalem 9190401, Israel
| |
Collapse
|
10
|
Lafiosca P, Gómez S, Giovannini T, Cappelli C. Absorption Properties of Large Complex Molecular Systems: The DFTB/Fluctuating Charge Approach. J Chem Theory Comput 2022; 18:1765-1779. [PMID: 35184553 PMCID: PMC8908768 DOI: 10.1021/acs.jctc.1c01066] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
![]()
We report on the
first formulation of a novel polarizable QM/MM
approach, where the density functional tight binding (DFTB) is coupled
to the fluctuating charge (FQ) force field. The resulting method (DFTB/FQ)
is then extended to the linear response within the TD-DFTB framework
and challenged to study absorption spectra of large condensed-phase
systems.
Collapse
Affiliation(s)
- Piero Lafiosca
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Sara Gómez
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Tommaso Giovannini
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| | - Chiara Cappelli
- Scuola Normale Superiore, Classe di Scienze, Piazza dei Cavalieri 7, 56126 Pisa, Italy
| |
Collapse
|
11
|
Nandy A, Duan C, Taylor MG, Liu F, Steeves AH, Kulik HJ. Computational Discovery of Transition-metal Complexes: From High-throughput Screening to Machine Learning. Chem Rev 2021; 121:9927-10000. [PMID: 34260198 DOI: 10.1021/acs.chemrev.1c00347] [Citation(s) in RCA: 70] [Impact Index Per Article: 23.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Transition-metal complexes are attractive targets for the design of catalysts and functional materials. The behavior of the metal-organic bond, while very tunable for achieving target properties, is challenging to predict and necessitates searching a wide and complex space to identify needles in haystacks for target applications. This review will focus on the techniques that make high-throughput search of transition-metal chemical space feasible for the discovery of complexes with desirable properties. The review will cover the development, promise, and limitations of "traditional" computational chemistry (i.e., force field, semiempirical, and density functional theory methods) as it pertains to data generation for inorganic molecular discovery. The review will also discuss the opportunities and limitations in leveraging experimental data sources. We will focus on how advances in statistical modeling, artificial intelligence, multiobjective optimization, and automation accelerate discovery of lead compounds and design rules. The overall objective of this review is to showcase how bringing together advances from diverse areas of computational chemistry and computer science have enabled the rapid uncovering of structure-property relationships in transition-metal chemistry. We aim to highlight how unique considerations in motifs of metal-organic bonding (e.g., variable spin and oxidation state, and bonding strength/nature) set them and their discovery apart from more commonly considered organic molecules. We will also highlight how uncertainty and relative data scarcity in transition-metal chemistry motivate specific developments in machine learning representations, model training, and in computational chemistry. Finally, we will conclude with an outlook of areas of opportunity for the accelerated discovery of transition-metal complexes.
Collapse
Affiliation(s)
- Aditya Nandy
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Chenru Duan
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Michael G Taylor
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Adam H Steeves
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
12
|
Desmarais JK, Erba A, Flament JP, Kirtman B. Perturbation Theory Treatment of Spin-Orbit Coupling II: A Coupled Perturbed Kohn-Sham Method. J Chem Theory Comput 2021; 17:4712-4732. [PMID: 34286577 DOI: 10.1021/acs.jctc.1c00460] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A noncanonical coupled perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment of spin-orbit coupling (SOC) is provided. We take the scalar-relativistic KS-DFT/HF solution, obtained with a relativistic effective core potential, as the zeroth-order approximation. Explicit expressions are given for the total energy through the 4th order, which satisfy the 2n + 1 rule. Second-order expressions are provided for orbital energies and density variables of spin-current DFT. Test calculations are carried out on the halogen homonuclear diatomic and hydride molecules, including 6p and 7p elements, as well as open-shell negative ions. The computed properties through second or third order match well with those from reference two-component self-consistent field calculations for total and orbital energies as well as spin-current densities. In only one case (At2-) did a significant deviation occur for the remaining density variables. Our coupled perturbation theory approach provides an efficient way of adding the effect of SOC to a scalar-relativistic single-reference KS-DFT/HF treatment, in particular because it does not require diagonalization in the two-component spinor basis, leading to saving factors on the number of required floating-point operations that may exceed one order of magnitude.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy.,Université de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, 64000 Pau, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
13
|
Liu F, Filatov M, Martínez TJ. Analytical derivatives of the individual state energies in ensemble density functional theory. II. Implementation on graphical processing units (GPUs). J Chem Phys 2021; 154:104108. [DOI: 10.1063/5.0041389] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/24/2023] Open
Affiliation(s)
- Fang Liu
- Department of Chemistry, Emory University, Atlanta, Georgia 30322, USA
| | - Michael Filatov
- Department of Chemistry, Kyungpook National University, Daegu 702-701, South Korea
| | - Todd J. Martínez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
14
|
Wolter M, von Looz M, Meyerhenke H, Jacob CR. Systematic Partitioning of Proteins for Quantum-Chemical Fragmentation Methods Using Graph Algorithms. J Chem Theory Comput 2021; 17:1355-1367. [PMID: 33591754 DOI: 10.1021/acs.jctc.0c01054] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
Quantum-chemical fragmentation methods offer an efficient approach for the treatment of large proteins, in particular if local target quantities such as protein-ligand interaction energies, enzymatic reaction energies, or spectroscopic properties of embedded chromophores are sought. However, the accuracy that is achievable for such local target quantities intricately depends on how the protein is partitioned into smaller fragments. While the commonly employed naı̈ve approach of using fragments with a fixed size is widely used, it can result in large and unpredictable errors when varying the fragment size. Here, we present a systematic partitioning scheme that aims at minimizing the fragmentation error of a local target quantity for a given maximum fragment size. To this end, we construct a weighted graph representation of the protein, in which the amino acids constitute the nodes. These nodes are connected by edges weighted with an estimate for the fragmentation error that is expected when cutting this edge. This allows us to employ graph partitioning algorithms provided by computer science to determine near-optimal partitions of the protein. We apply this scheme to a test set of six proteins representing various prototypical applications of quantum-chemical fragmentation methods using a simplified molecular fractionation with conjugate caps (MFCC) approach with hydrogen caps. We show that our graph-based scheme consistently improves upon the naı̈ve approach.
Collapse
Affiliation(s)
- Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstrasse 17, 38106 Braunschweig, Germany
| | - Moritz von Looz
- Department of Computer Science, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
| | - Henning Meyerhenke
- Department of Computer Science, Humboldt-Universität zu Berlin, Unter den Linden 6, 10099 Berlin, Germany
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstrasse 17, 38106 Braunschweig, Germany
| |
Collapse
|
15
|
Ryley MS, Withnall M, Irons TJP, Helgaker T, Teale AM. Robust All-Electron Optimization in Orbital-Free Density-Functional Theory Using the Trust-Region Image Method. J Phys Chem A 2021; 125:459-475. [PMID: 33356245 DOI: 10.1021/acs.jpca.0c09502] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a Gaussian-basis implementation of orbital-free density-functional theory (OF-DFT) in which the trust-region image method (TRIM) is used for optimization. This second-order optimization scheme has been constructed to provide benchmark all-electron results with very tight convergence of the particle-number constraint, associated chemical potential, and electron density. It is demonstrated that, by preserving the saddle-point nature of the optimization and simultaneously optimizing the density and chemical potential, an order of magnitude reduction in the number of iterations required for convergence is obtained. The approach is compared and contrasted with a new implementation of the nested optimization scheme put forward by Chan, Cohen, and Handy. Our implementation allows for semilocal kinetic-energy (and exchange-correlation) functionals to be handled self-consistently in all-electron calculations. The all-electron Gaussian-basis setting for these calculations will enable direct comparison with a wide range of standard high-accuracy quantum-chemical methods as well as with Kohn-Sham density-functional theory. We expect that the present implementation will provide a useful tool for analyzing the performance of approximate kinetic-energy functionals in finite systems.
Collapse
Affiliation(s)
- Matthew S Ryley
- School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, U.K
| | - Michael Withnall
- School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, U.K
| | - Tom J P Irons
- School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, U.K
| | - Trygve Helgaker
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| | - Andrew M Teale
- School of Chemistry, University of Nottingham, University Park, Nottingham, NG7 2RD, U.K.,Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
| |
Collapse
|
16
|
Mehmood R, Kulik HJ. Both Configuration and QM Region Size Matter: Zinc Stability in QM/MM Models of DNA Methyltransferase. J Chem Theory Comput 2020; 16:3121-3134. [PMID: 32243149 DOI: 10.1021/acs.jctc.0c00153] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/18/2022]
Abstract
Quantum-mechanical/molecular-mechanical (QM/MM) methods are essential to the study of metalloproteins, but the relative importance of sampling and degree of QM treatment in achieving quantitative predictions is poorly understood. We study the relative magnitude of configurational and QM-region sensitivity of energetic and electronic properties in a representative Zn2+ metal binding site of a DNA methyltransferase. To quantify property variations, we analyze snapshots extracted from 250 ns of molecular dynamics simulation. To understand the degree of QM-region sensitivity, we perform analysis using QM regions ranging from a minimal 49-atom region consisting only of the Zn2+ metal and its four coordinating Cys residues up to a 628-atom QM region that includes residues within 12 Å of the metal center. Over the configurations sampled, we observe that illustrative properties (e.g., rigid Zn2+ removal energy) exhibit large fluctuations that are well captured with even minimal QM regions. Nevertheless, for both energetic and electronic properties, we observe a slow approach to asymptotic limits with similarly large changes in absolute values that converge only with larger (ca. 300-atom) QM region sizes. For the smaller QM regions, the electronic description of Zn2+ binding is incomplete: the metal binds too tightly and is too stabilized by the strong electrostatic potential of MM point charges, and the Zn-S bond covalency is overestimated. Overall, this work suggests that efficient sampling with QM/MM in small QM regions is an effective method to explore the influence of enzyme structure on target properties. At the same time, accurate descriptions of electronic and energetic properties require a larger QM region than the minimal metal-coordinating residues in order to converge treatment of both metal-local bonding and the overall electrostatic environment.
Collapse
Affiliation(s)
- Rimsha Mehmood
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States.,Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
17
|
Mackie CJ, Gonthier JF, Head-Gordon M. Compressed intramolecular dispersion interactions. J Chem Phys 2020; 152:024112. [PMID: 31941286 DOI: 10.1063/1.5126716] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
The feasibility of the compression of localized virtual orbitals is explored in the context of intramolecular long-range dispersion interactions. Singular value decomposition (SVD) of coupled cluster doubles amplitudes associated with the dispersion interactions is analyzed for a number of long-chain systems, including saturated and unsaturated hydrocarbons and a silane chain. Further decomposition of the most important amplitudes obtained from these SVDs allows for the analysis of the dispersion-specific virtual orbitals that are naturally localized. Consistent with previous work on intermolecular dispersion interactions in dimers, it is found that three important geminals arise and account for the majority of dispersion interactions at the long range, even in the many body intramolecular case. Furthermore, it is shown that as few as three localized virtual orbitals per occupied orbital can be enough to capture all pairwise long-range dispersion interactions within a molecule.
Collapse
Affiliation(s)
- Cameron J Mackie
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Jérôme F Gonthier
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Kenneth S. Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
18
|
Rinkevicius Z, Li X, Vahtras O, Ahmadzadeh K, Brand M, Ringholm M, List NH, Scheurer M, Scott M, Dreuw A, Norman P. VeloxChem: A Python‐driven density‐functional theory program for spectroscopy simulations in high‐performance computing environments. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1457] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Zilvinas Rinkevicius
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
- Department of Physics Kaunas University of Technology Kaunas Lithuania
| | - Xin Li
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Olav Vahtras
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Karan Ahmadzadeh
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Manuel Brand
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Magnus Ringholm
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| | - Nanna Holmgaard List
- Department of Chemistry and the PULSE Institute Stanford University Stanford California
- SLAC National Accelerator Laboratory Menlo Park California
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing Ruprecht‐Karls University Heidelberg Germany
| | - Mikael Scott
- Interdisciplinary Center for Scientific Computing Ruprecht‐Karls University Heidelberg Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing Ruprecht‐Karls University Heidelberg Germany
| | - Patrick Norman
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health KTH Royal Institute of Technology Stockholm Sweden
| |
Collapse
|
19
|
Lee J, Lin L, Head-Gordon M. Systematically Improvable Tensor Hypercontraction: Interpolative Separable Density-Fitting for Molecules Applied to Exact Exchange, Second- and Third-Order Møller–Plesset Perturbation Theory. J Chem Theory Comput 2019; 16:243-263. [DOI: 10.1021/acs.jctc.9b00820] [Citation(s) in RCA: 25] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Lin Lin
- Department of Mathematics, University of California, Berkeley, California 94720, United States
- Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| |
Collapse
|
20
|
Veccham SP, Lee J, Head-Gordon M. Making many-body interactions nearly pairwise additive: The polarized many-body expansion approach. J Chem Phys 2019; 151:194101. [DOI: 10.1063/1.5125802] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/10/2023] Open
Affiliation(s)
- Srimukh Prasad Veccham
- Department of Chemistry, University of California, Berkeley, California 94720 USA, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Joonho Lee
- Department of Chemistry, University of California, Berkeley, California 94720 USA, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Department of Chemistry, University of California, Berkeley, California 94720 USA, and Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, USA
| |
Collapse
|
21
|
Kulik HJ. Making machine learning a useful tool in the accelerated discovery of transition metal complexes. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1439] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Affiliation(s)
- Heather J. Kulik
- Department of Chemical Engineering Massachusetts Institute of Technology Cambridge Massachusetts
| |
Collapse
|
22
|
Stochastic density functional theory. WILEY INTERDISCIPLINARY REVIEWS-COMPUTATIONAL MOLECULAR SCIENCE 2019. [DOI: 10.1002/wcms.1412] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/07/2022]
|
23
|
Qi HW, Kulik HJ. Evaluating Unexpectedly Short Non-covalent Distances in X-ray Crystal Structures of Proteins with Electronic Structure Analysis. J Chem Inf Model 2019; 59:2199-2211. [DOI: 10.1021/acs.jcim.9b00144] [Citation(s) in RCA: 30] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Helena W. Qi
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
- Department of Chemistry, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
24
|
Janet JP, Liu F, Nandy A, Duan C, Yang T, Lin S, Kulik HJ. Designing in the Face of Uncertainty: Exploiting Electronic Structure and Machine Learning Models for Discovery in Inorganic Chemistry. Inorg Chem 2019; 58:10592-10606. [PMID: 30834738 DOI: 10.1021/acs.inorgchem.9b00109] [Citation(s) in RCA: 59] [Impact Index Per Article: 11.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Recent transformative advances in computing power and algorithms have made computational chemistry central to the discovery and design of new molecules and materials. First-principles simulations are increasingly accurate and applicable to large systems with the speed needed for high-throughput computational screening. Despite these strides, the combinatorial challenges associated with the vastness of chemical space mean that more than just fast and accurate computational tools are needed for accelerated chemical discovery. In transition-metal chemistry and catalysis, unique challenges arise. The variable spin, oxidation state, and coordination environments favored by elements with well-localized d or f electrons provide great opportunity for tailoring properties in catalytic or functional (e.g., magnetic) materials but also add layers of uncertainty to any design strategy. We outline five key mandates for realizing computationally driven accelerated discovery in inorganic chemistry: (i) fully automated simulation of new compounds, (ii) knowledge of prediction sensitivity or accuracy, (iii) faster-than-fast property prediction methods, (iv) maps for rapid chemical space traversal, and (v) a means to reveal design rules on the kilocompound scale. Through case studies in open-shell transition-metal chemistry, we describe how advances in methodology and software in each of these areas bring about new chemical insights. We conclude with our outlook on the next steps in this process toward realizing fully autonomous discovery in inorganic chemistry using computational chemistry.
Collapse
Affiliation(s)
- Jon Paul Janet
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Fang Liu
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Aditya Nandy
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Chenru Duan
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States.,Department of Chemistry , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Tzuhsiung Yang
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Sean Lin
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| | - Heather J Kulik
- Department of Chemical Engineering , Massachusetts Institute of Technology , Cambridge , Massachusetts 02139 , United States
| |
Collapse
|
25
|
Wu X, Dral PO, Koslowski A, Thiel W. Big data analysis of ab Initio molecular integrals in the neglect of diatomic differential overlap approximation. J Comput Chem 2019; 40:638-649. [PMID: 30549072 DOI: 10.1002/jcc.25748] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2018] [Revised: 10/15/2018] [Accepted: 10/16/2018] [Indexed: 12/14/2022]
Abstract
Most modern semiempirical quantum-chemical (SQC) methods are based on the neglect of diatomic differential overlap (NDDO) approximation to ab initio molecular integrals. Here, we check the validity of this approximation by computing all relevant integrals for 32 typical organic molecules using Gaussian-type orbitals and various basis sets (from valence-only minimal to all-electron triple-ζ basis sets) covering in total more than 15.6 million one-electron (1-e) and 10.3 billion two-electron (2-e) integrals. The integrals are calculated in the nonorthogonal atomic basis and then transformed by symmetric orthogonalization to the Löwdin basis. In the case of the 1-e integrals, we find strong orthogonalization effects that need to be included in SQC models, for example, by strategies such as those adopted in the available OMx methods. For the valence-only minimal basis, we confirm that the 2-e Coulomb integrals in the Löwdin basis are quantitatively close to their counterparts in the atomic basis and that the 2-e exchange integrals can be safely neglected in line with the NDDO approximation. For larger all-electron basis sets, there are strong multishell orthogonalization effects that lead to more irregular patterns in the transformed 2-e integrals and thus cast doubt on the validity of the NDDO approximation for extended basis sets. Focusing on the valence-only minimal basis, we find that some of the NDDO-neglected integrals are reduced but remain sizable after the transformation to the Löwdin basis; this is true for the two-center 2-e hybrid integrals, the three-center 1-e nuclear attraction integrals, and the corresponding three-center 2-e hybrid integrals. We consider a scheme with a valence-only minimal basis that includes such terms as a possible strategy to go beyond the NDDO integral approximation in attempts to improve SQC methods. © 2018 Wiley Periodicals, Inc.
Collapse
Affiliation(s)
- Xin Wu
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| | - Pavlo O Dral
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| | - Axel Koslowski
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| | - Walter Thiel
- Max-Planck-Institut für Kohlenforschung, 45470, Mülheim an der Ruhr, Germany
| |
Collapse
|
26
|
Grimme S, Schreiner PR. Computerchemie: das Schicksal aktueller Methoden und zukünftige Herausforderungen. Angew Chem Int Ed Engl 2017. [DOI: 10.1002/ange.201709943] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry; Universität Bonn; Beringstraße 4 53115 Bonn Deutschland
| | - Peter R. Schreiner
- Institut für Organische Chemie; Justus-Liebig-Universität; Heinrich-Buff-Ring 17 35392 Gießen Deutschland
| |
Collapse
|
27
|
Grimme S, Schreiner PR. Computational Chemistry: The Fate of Current Methods and Future Challenges. Angew Chem Int Ed Engl 2017; 57:4170-4176. [PMID: 29105929 DOI: 10.1002/anie.201709943] [Citation(s) in RCA: 106] [Impact Index Per Article: 15.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Indexed: 11/12/2022]
Abstract
"Where do we go from here?" is the underlying question regarding the future (perhaps foreseeable) developments in computational chemistry. Although this young discipline has already permeated practically all of chemistry, it is likely to become even more powerful with the rapid development of computational hard- and software.
Collapse
Affiliation(s)
- Stefan Grimme
- Mulliken Center for Theoretical Chemistry, University of Bonn, Beringstr. 4, 53115, Bonn, Germany
| | - Peter R Schreiner
- Institute of Organic Chemistry, Justus-Liebig University, Heinrich-Buff-Ring 17, 35392, Giessen, Germany
| |
Collapse
|
28
|
Ramos-Cordoba E, Matito E. Local Descriptors of Dynamic and Nondynamic Correlation. J Chem Theory Comput 2017; 13:2705-2711. [DOI: 10.1021/acs.jctc.7b00293] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Affiliation(s)
- Eloy Ramos-Cordoba
- Kimika
Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, and Donostia International Physics Center (DIPC), P.K. 1072, 20080 Donostia, Euskadi Spain
| | - Eduard Matito
- Kimika
Fakultatea, Euskal Herriko Unibertsitatea, UPV/EHU, and Donostia International Physics Center (DIPC), P.K. 1072, 20080 Donostia, Euskadi Spain
- IKERBASQUE, Basque Foundation for Science, 48011 Bilbao, Spain
| |
Collapse
|
29
|
Karelina M, Kulik HJ. Systematic Quantum Mechanical Region Determination in QM/MM Simulation. J Chem Theory Comput 2017; 13:563-576. [DOI: 10.1021/acs.jctc.6b01049] [Citation(s) in RCA: 57] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Affiliation(s)
- Maria Karelina
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| | - Heather J. Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, United States
| |
Collapse
|
30
|
Kumar C, Kjærgaard T, Helgaker T, Fliegl H. Nuclei-selected atomic-orbital response-theory formulation for the calculation of NMR shielding tensors using density-fitting. J Chem Phys 2016; 145:234108. [DOI: 10.1063/1.4972212] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023] Open
Affiliation(s)
- Chandan Kumar
- Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Thomas Kjærgaard
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000 Aarhus C, Denmark
| | - Trygve Helgaker
- Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| | - Heike Fliegl
- Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, P.O. Box 1033, Blindern, 0315 Oslo, Norway
| |
Collapse
|
31
|
Akimov AV. Nonadiabatic Molecular Dynamics with Tight-Binding Fragment Molecular Orbitals. J Chem Theory Comput 2016; 12:5719-5736. [DOI: 10.1021/acs.jctc.6b00955] [Citation(s) in RCA: 33] [Impact Index Per Article: 4.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023]
Affiliation(s)
- Alexey V. Akimov
- Department of Chemistry, University at Buffalo, The State University of New York, Buffalo, New York 14260-3000, United States
| |
Collapse
|
32
|
Zen A, Roch LM, Cox SJ, Hu XL, Sorella S, Alfè D, Michaelides A. Toward Accurate Adsorption Energetics on Clay Surfaces. THE JOURNAL OF PHYSICAL CHEMISTRY. C, NANOMATERIALS AND INTERFACES 2016; 120:26402-26413. [PMID: 27917256 PMCID: PMC5126707 DOI: 10.1021/acs.jpcc.6b09559] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2016] [Revised: 10/31/2016] [Indexed: 05/31/2023]
Abstract
Clay minerals are ubiquitous in nature, and the manner in which they interact with their surroundings has important industrial and environmental implications. Consequently, a molecular-level understanding of the adsorption of molecules on clay surfaces is crucial. In this regard computer simulations play an important role, yet the accuracy of widely used empirical force fields (FF) and density functional theory (DFT) exchange-correlation functionals is often unclear in adsorption systems dominated by weak interactions. Herein we present results from quantum Monte Carlo (QMC) for water and methanol adsorption on the prototypical clay kaolinite. To the best of our knowledge, this is the first time QMC has been used to investigate adsorption at a complex, natural surface such as a clay. As well as being valuable in their own right, the QMC benchmarks obtained provide reference data against which the performance of cheaper DFT methods can be tested. Indeed using various DFT exchange-correlation functionals yields a very broad range of adsorption energies, and it is unclear a priori which evaluation is better. QMC reveals that in the systems considered here it is essential to account for van der Waals (vdW) dispersion forces since this alters both the absolute and relative adsorption energies of water and methanol. We show, via FF simulations, that incorrect relative energies can lead to significant changes in the interfacial densities of water and methanol solutions at the kaolinite interface. Despite the clear improvements offered by the vdW-corrected and the vdW-inclusive functionals, absolute adsorption energies are often overestimated, suggesting that the treatment of vdW forces in DFT is not yet a solved problem.
Collapse
Affiliation(s)
- Andrea Zen
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Loïc M. Roch
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Stephen J. Cox
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Chemistry, University College London, 20 Gordon Street, London WC1H 0AJ, United Kingdom
| | - Xiao Liang Hu
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Sandro Sorella
- SISSA−International
School for Advanced Studies, Via Bonomea 26, 34136 Trieste, Italy
- INFM
Democritos National Simulation Center, 34151 Trieste, Italy
| | - Dario Alfè
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Earth Sciences, University College London, Gower Street, London WC1E 6BT, United
Kingdom
| | - Angelos Michaelides
- Thomas
Young Centre and London Centre for Nanotechnology, 17−19 Gordon Street, London WC1H 0AH, United
Kingdom
- Department
of Physics and Astronomy, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|
33
|
Lewis CA, Calvin JA, Valeev EF. Clustered Low-Rank Tensor Format: Introduction and Application to Fast Construction of Hartree–Fock Exchange. J Chem Theory Comput 2016; 12:5868-5880. [DOI: 10.1021/acs.jctc.6b00884] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Affiliation(s)
- Cannada A. Lewis
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Justus A. Calvin
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|
34
|
Kulik H, Zhang J, Klinman J, Martínez TJ. How Large Should the QM Region Be in QM/MM Calculations? The Case of Catechol O-Methyltransferase. J Phys Chem B 2016; 120:11381-11394. [PMID: 27704827 PMCID: PMC5108028 DOI: 10.1021/acs.jpcb.6b07814] [Citation(s) in RCA: 139] [Impact Index Per Article: 17.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/02/2016] [Revised: 09/09/2016] [Indexed: 01/29/2023]
Abstract
Hybrid quantum mechanical-molecular mechanical (QM/MM) simulations are widely used in studies of enzymatic catalysis. Until recently, it has been cost prohibitive to determine the asymptotic limit of key energetic and structural properties with respect to increasingly large QM regions. Leveraging recent advances in electronic structure efficiency and accuracy, we investigate catalytic properties in catechol O-methyltransferase, a prototypical methyltransferase critical to human health. Using QM regions ranging in size from reactants-only (64 atoms) to nearly one-third of the entire protein (940 atoms), we show that properties such as the activation energy approach within chemical accuracy of the large-QM asymptotic limits rather slowly, requiring approximately 500-600 atoms if the QM residues are chosen simply by distance from the substrate. This slow approach to asymptotic limit is due to charge transfer from protein residues to the reacting substrates. Our large QM/MM calculations enable identification of charge separation for fragments in the transition state as a key component of enzymatic methyl transfer rate enhancement. We introduce charge shift analysis that reveals the minimum number of protein residues (approximately 11-16 residues or 200-300 atoms for COMT) needed for quantitative agreement with large-QM simulations. The identified residues are not those that would be typically selected using criteria such as chemical intuition or proximity. These results provide a recipe for a more careful determination of QM region sizes in future QM/MM studies of enzymes.
Collapse
Affiliation(s)
- Heather
J. Kulik
- Department
of Chemistry and PULSE Institute, Stanford
University, Stanford, California 94305, United States
- SLAC
National Accelerator Laboratory, Menlo Park, California 94025, United States
| | - Jianyu Zhang
- Departments
of Chemistry and of Molecular and Cell Biology, and California Institute
for Quantitative Biosciences, University
of California, Berkeley, California 94720, United States
| | - Judith
P. Klinman
- Departments
of Chemistry and of Molecular and Cell Biology, and California Institute
for Quantitative Biosciences, University
of California, Berkeley, California 94720, United States
| | - Todd J. Martínez
- Department
of Chemistry and PULSE Institute, Stanford
University, Stanford, California 94305, United States
- SLAC
National Accelerator Laboratory, Menlo Park, California 94025, United States
| |
Collapse
|
35
|
von Looz M, Wolter M, Jacob CR, Meyerhenke H. Better Partitions of Protein Graphs for Subsystem Quantum Chemistry. EXPERIMENTAL ALGORITHMS 2016. [DOI: 10.1007/978-3-319-38851-9_24] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
36
|
Hollman DS, Schaefer HF, Valeev EF. A tight distance-dependent estimator for screening three-center Coulomb integrals over Gaussian basis functions. J Chem Phys 2015; 142:154106. [DOI: 10.1063/1.4917519] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- David S. Hollman
- Center for Computational Quantum Chemistry, University of Georgia, 1004 Cedar St., Athens, Georgia 30602, USA
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, 1004 Cedar St., Athens, Georgia 30602, USA
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
37
|
Vereecken L, Glowacki DR, Pilling MJ. Theoretical Chemical Kinetics in Tropospheric Chemistry: Methodologies and Applications. Chem Rev 2015; 115:4063-114. [DOI: 10.1021/cr500488p] [Citation(s) in RCA: 146] [Impact Index Per Article: 16.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
Affiliation(s)
- Luc Vereecken
- Max Planck Institute for Chemistry, 55128 Mainz, Germany
| | - David R. Glowacki
- PULSE
Institute and Department of Chemistry, Stanford University, Stanford, California 94305, United States
- SLAC National Accelerator Laboratory, Menlo Park, California 94025, United States
- School
of Chemistry, University of Bristol, Bristol BS8 1TS, United Kingdom
- Department
of Computer Science, University of Bristol, Bristol BS8 1UB, United Kingdom
| | | |
Collapse
|
38
|
McClean JR, Babbush R, Love PJ, Aspuru-Guzik A. Exploiting Locality in Quantum Computation for Quantum Chemistry. J Phys Chem Lett 2014; 5:4368-80. [PMID: 26273989 DOI: 10.1021/jz501649m] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 05/26/2023]
Abstract
Accurate prediction of chemical and material properties from first-principles quantum chemistry is a challenging task on traditional computers. Recent developments in quantum computation offer a route toward highly accurate solutions with polynomial cost; however, this solution still carries a large overhead. In this Perspective, we aim to bring together known results about the locality of physical interactions from quantum chemistry with ideas from quantum computation. We show that the utilization of spatial locality combined with the Bravyi-Kitaev transformation offers an improvement in the scaling of known quantum algorithms for quantum chemistry and provides numerical examples to help illustrate this point. We combine these developments to improve the outlook for the future of quantum chemistry on quantum computers.
Collapse
Affiliation(s)
- Jarrod R McClean
- †Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Ryan Babbush
- †Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Peter J Love
- ‡Department of Physics, Haverford College, Haverford, Pennsylvania 19041, United States
| | - Alán Aspuru-Guzik
- †Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
39
|
Interactive Chemical Reactivity Exploration. Chemphyschem 2014; 15:3301-19. [DOI: 10.1002/cphc.201402342] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2014] [Revised: 07/07/2014] [Indexed: 11/07/2022]
|
40
|
Hedegård ED, List NH, Jensen HJA, Kongsted J. The multi-configuration self-consistent field method within a polarizable embedded framework. J Chem Phys 2013; 139:044101. [DOI: 10.1063/1.4811835] [Citation(s) in RCA: 43] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
41
|
Kiewisch K, Jacob CR, Visscher L. Quantum-Chemical Electron Densities of Proteins and of Selected Protein Sites from Subsystem Density Functional Theory. J Chem Theory Comput 2013; 9:2425-40. [DOI: 10.1021/ct3008759] [Citation(s) in RCA: 38] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Karin Kiewisch
- Amsterdam Center for Multiscale
Modeling, VU University Amsterdam, De Boelelaan
1083, 1081 HV Amsterdam, The Netherlands
| | - Christoph R. Jacob
- Center for Functional Nanostructures
and Institute of Physical Chemistry, Karlsruhe Institute of Technology (KIT), Wolfgang-Gaede-Str. 1a, 76131 Karlsruhe,
Germany
| | - Lucas Visscher
- Amsterdam Center for Multiscale
Modeling, VU University Amsterdam, De Boelelaan
1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
42
|
Maurer SA, Lambrecht DS, Kussmann J, Ochsenfeld C. Efficient distance-including integral screening in linear-scaling Møller-Plesset perturbation theory. J Chem Phys 2013; 138:014101. [DOI: 10.1063/1.4770502] [Citation(s) in RCA: 60] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
43
|
Kruse H, Grimme S. A geometrical correction for the inter- and intra-molecular basis set superposition error in Hartree-Fock and density functional theory calculations for large systems. J Chem Phys 2012; 136:154101. [DOI: 10.1063/1.3700154] [Citation(s) in RCA: 414] [Impact Index Per Article: 34.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
44
|
Maurer SA, Lambrecht DS, Flaig D, Ochsenfeld C. Distance-dependent Schwarz-based integral estimates for two-electron integrals: Reliable tightness vs. rigorous upper bounds. J Chem Phys 2012; 136:144107. [DOI: 10.1063/1.3693908] [Citation(s) in RCA: 64] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
45
|
Mikulskis P, Genheden S, Wichmann K, Ryde U. A semiempirical approach to ligand-binding affinities: Dependence on the Hamiltonian and corrections. J Comput Chem 2012; 33:1179-89. [DOI: 10.1002/jcc.22949] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2011] [Revised: 01/10/2012] [Accepted: 01/16/2012] [Indexed: 12/30/2022]
|
46
|
Bowler DR, Miyazaki T. O(N) methods in electronic structure calculations. REPORTS ON PROGRESS IN PHYSICS. PHYSICAL SOCIETY (GREAT BRITAIN) 2012; 75:036503. [PMID: 22790422 DOI: 10.1088/0034-4885/75/3/036503] [Citation(s) in RCA: 210] [Impact Index Per Article: 17.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/21/2023]
Abstract
Linear-scaling methods, or O(N) methods, have computational and memory requirements which scale linearly with the number of atoms in the system, N, in contrast to standard approaches which scale with the cube of the number of atoms. These methods, which rely on the short-ranged nature of electronic structure, will allow accurate, ab initio simulations of systems of unprecedented size. The theory behind the locality of electronic structure is described and related to physical properties of systems to be modelled, along with a survey of recent developments in real-space methods which are important for efficient use of high-performance computers. The linear-scaling methods proposed to date can be divided into seven different areas, and the applicability, efficiency and advantages of the methods proposed in these areas are then discussed. The applications of linear-scaling methods, as well as the implementations available as computer programs, are considered. Finally, the prospects for and the challenges facing linear-scaling methods are discussed.
Collapse
Affiliation(s)
- D R Bowler
- London Centre for Nanotechnology, UCL, 17-19 Gordon St, London WC1H 0AH, UK.
| | | |
Collapse
|
47
|
Spencer JS, Blunt NS, Foulkes WM. The sign problem and population dynamics in the full configuration interaction quantum Monte Carlo method. J Chem Phys 2012; 136:054110. [DOI: 10.1063/1.3681396] [Citation(s) in RCA: 83] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
48
|
Kristensen K, Ziółkowski M, Jansík B, Kjærgaard T, Jørgensen P. A Locality Analysis of the Divide–Expand–Consolidate Coupled Cluster Amplitude Equations. J Chem Theory Comput 2011; 7:1677-94. [DOI: 10.1021/ct200114k] [Citation(s) in RCA: 85] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Affiliation(s)
- Kasper Kristensen
- Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, DK-8000 Århus C, Denmark
| | - Marcin Ziółkowski
- Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, DK-8000 Århus C, Denmark
| | - Branislav Jansík
- Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, DK-8000 Århus C, Denmark
| | - Thomas Kjærgaard
- Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, DK-8000 Århus C, Denmark
| | - Poul Jørgensen
- Lundbeck Foundation Center for Theoretical Chemistry, Department of Chemistry, University of Aarhus, DK-8000 Århus C, Denmark
| |
Collapse
|
49
|
Jacobson LD, Herbert JM. An efficient, fragment-based electronic structure method for molecular systems: Self-consistent polarization with perturbative two-body exchange and dispersion. J Chem Phys 2011; 134:094118. [DOI: 10.1063/1.3560026] [Citation(s) in RCA: 77] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
50
|
Beer M, Kussmann J, Ochsenfeld C. Nuclei-selected NMR shielding calculations: A sublinear-scaling quantum-chemical method. J Chem Phys 2011; 134:074102. [DOI: 10.1063/1.3526315] [Citation(s) in RCA: 42] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|