1
|
Petrov K, Csóka J, Kállay M. Analytic Gradients for Density Fitting MP2 Using Natural Auxiliary Functions. J Phys Chem A 2024; 128:6566-6580. [PMID: 39074307 PMCID: PMC11317987 DOI: 10.1021/acs.jpca.4c02822] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/30/2024] [Revised: 07/02/2024] [Accepted: 07/16/2024] [Indexed: 07/31/2024]
Abstract
The natural auxiliary function (NAF) approach is an approximation to decrease the size of the auxiliary basis set required for quantum chemical calculations utilizing the density fitting technique. It has been proven efficient to speed up various correlation models, such as second-order Møller-Plesset (MP2) theory and coupled-cluster methods. Here, for the first time, we discuss the theory of analytic derivatives for correlation methods employing the NAF approximation on the example of MP2. A detailed algorithm for the gradient calculation with the NAF approximation is proposed in the framework of the method of Lagrange multipliers. To assess the effect of the NAF approximation on gradients and optimized geometric parameters, a series of benchmark calculations on small and medium-sized systems was performed. Our results demonstrate that, for MP2, sufficiently accurate gradients and geometries can be achieved with a moderate time reduction of 15-20% for both small and medium-sized molecules.
Collapse
Affiliation(s)
- Klára Petrov
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN−BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA−BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - József Csóka
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN−BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA−BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| | - Mihály Kállay
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- HUN-REN−BME
Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
- MTA−BME
Lendület Quantum Chemistry Research Group, Műegyetem rkp. 3., H-1111 Budapest, Hungary
| |
Collapse
|
2
|
Panchagnula K, Graf D, Johnson ER, Thom AJW. Targeting spectroscopic accuracy for dispersion bound systems from ab initio techniques: Translational eigenstates of Ne@C70. J Chem Phys 2024; 161:054308. [PMID: 39092939 DOI: 10.1063/5.0223298] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2024] [Accepted: 07/17/2024] [Indexed: 08/04/2024] Open
Abstract
We investigate the endofullerene system Ne@C70 by constructing a three-dimensional Potential Energy Surface (PES) describing the translational motion of the Ne atom. This is constructed from electronic structure calculations from a plethora of methods, including MP2, SCS-MP2, SOS-MP2, RPA@PBE, and C(HF)-RPA, which were previously used for He@C60 in Panchagnula et al. [J. Chem. Phys. 160, 104303 (2024)], alongside B86bPBE-25X-XDM and B86bPBE-50X-XDM. The reduction in symmetry moving from C60 to C70 introduces a double well potential along the anisotropic direction, which forms a test of the sensitivity and effectiveness of the electronic structure methods. The nuclear Hamiltonian is diagonalized using a symmetrized double minimum basis set outlined in Panchagnula and Thom [J. Chem. Phys. 159, 164308 (2023)], with translational energies having error bars ±1 and ±2 cm-1. We find no consistency between electronic structure methods as they find a range of barrier heights and minima positions of the double well and different translational eigenspectra, which also differ from the Lennard-Jones (LJ) PES given in Mandziuk and Bačić [J. Chem. Phys. 101, 2126-2140 (1994)]. We find that generating effective LJ parameters for each electronic structure method cannot reproduce the full PES nor recreate the eigenstates, and this suggests that the LJ form of the PES, while simple, may not be best suited to describe these systems. Even though MP2 and RPA@PBE performed best for He@C60, due to the lack of concordance between all electronic structure methods, we require more experimental data in order to properly validate the choice.
Collapse
Affiliation(s)
- K Panchagnula
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - D Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- Department of Chemistry, University of Munich (LMU), Munich, Germany
| | - E R Johnson
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
- Department of Chemistry, Dalhousie University, 6243 Alumni Crescent, Halifax, Nova Scotia B3H 4R2, Canada
| | - A J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
3
|
Zheng Y, Ni Z, Wang Y, Li W, Li S. Analytical Gradient Using Cluster-in-Molecule RI-MP2 Method for the Geometry Optimizations of Large Systems. J Chem Theory Comput 2024; 20:3626-3636. [PMID: 38626287 DOI: 10.1021/acs.jctc.4c00087] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/18/2024]
Abstract
We present an efficient analytical energy gradient algorithm for the cluster-in-molecule resolution-of-identity second-order Møller-Plesset perturbation (CIM-RI-MP2) method based on the Lagrange multiplier method. Our algorithm independently constructs the Lagrangian formalism within each cluster, avoiding the solution of the coupled-perturbed Hartree-Fock (CPHF) equation for the whole system. Due to this feature, the computational cost of the CIM-RI-MP2 gradients is much lower than that of other local MP2 algorithms. Benchmark calculations of several molecules containing up to 312 atoms demonstrate the general applicability of our CIM-RI-MP2 gradient algorithm. The optimized structure of a 244-atom molecule using the CIM-RI-MP2 method with the cc-pVDZ basis set is in good agreement with the corresponding crystal structure. A single-point gradient calculation conducted for a molecular cage containing 972 atoms and 9612 basis functions takes 48 h on 25 nodes, utilizing a total of 600 CPU cores. The present CIM-RI-MP2 gradient program is applicable for obtaining the optimized geometries of large systems with hundreds of atoms.
Collapse
Affiliation(s)
- Yang Zheng
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology of Ministry of Education, Hangzhou Normal University, Hangzhou 311121, P. R. China
| | - Yuqi Wang
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Wei Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| | - Shuhua Li
- Key Laboratory of Mesoscopic Chemistry of Ministry of Education, New Cornerstone Science Laboratory, School of Chemistry and Chemical Engineering, Nanjing University, Nanjing, Jiangsu 210023, P. R. China
| |
Collapse
|
4
|
Stocks R, Palethorpe E, Barca GMJ. High-Performance Multi-GPU Analytic RI-MP2 Energy Gradients. J Chem Theory Comput 2024; 20:2505-2519. [PMID: 38456899 DOI: 10.1021/acs.jctc.3c01424] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/09/2024]
Abstract
This article presents a novel algorithm for the calculation of analytic energy gradients from second-order Møller-Plesset perturbation theory within the Resolution-of-the-Identity approximation (RI-MP2), which is designed to achieve high performance on clusters with multiple graphical processing units (GPUs). The algorithm uses GPUs for all major steps of the calculation, including integral generation, formation of all required intermediate tensors, solution of the Z-vector equation and gradient accumulation. The implementation in the EXtreme Scale Electronic Structure System (EXESS) software package includes a tailored, highly efficient, multistream scheduling system to hide CPU-GPU data transfer latencies and allows nodes with 8 A100 GPUs to operate at over 80% of theoretical peak floating-point performance. Comparative performance analysis shows a significant reduction in computational time relative to traditional multicore CPU-based methods, with our approach achieving up to a 95-fold speedup over the single-node performance of established software such as Q-Chem and ORCA. Additionally, we demonstrate that pairing our implementation with the molecular fragmentation framework in EXESS can drastically lower the computational scaling of RI-MP2 gradient calculations from quintic to subquadratic, enabling further substantial savings in runtime while retaining high numerical accuracy in the resulting gradients.
Collapse
Affiliation(s)
- Ryan Stocks
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| | - Elise Palethorpe
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| | - Giuseppe M J Barca
- School of Computing, Australian National University, Canberra, ACT 2601, Australia
| |
Collapse
|
5
|
Panchagnula K, Graf D, Albertani FEA, Thom AJW. Translational eigenstates of He@C60 from four-dimensional ab initio potential energy surfaces interpolated using Gaussian process regression. J Chem Phys 2024; 160:104303. [PMID: 38465682 DOI: 10.1063/5.0197903] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/15/2024] [Accepted: 02/22/2024] [Indexed: 03/12/2024] Open
Abstract
We investigate the endofullerene system 3He@C60 with a four-dimensional potential energy surface (PES) to include the three He translational degrees of freedom and C60 cage radius. We compare second order Møller-Plesset perturbation theory (MP2), spin component scaled-MP2, scaled opposite spin-MP2, random phase approximation (RPA)@Perdew, Burke, and Ernzerhof (PBE), and corrected Hartree-Fock-RPA to calibrate and gain confidence in the choice of electronic structure method. Due to the high cost of these calculations, the PES is interpolated using Gaussian Process Regression (GPR), owing to its effectiveness with sparse training data. The PES is split into a two-dimensional radial surface, to which corrections are applied to achieve an overall four-dimensional surface. The nuclear Hamiltonian is diagonalized to generate the in-cage translational/vibrational eigenstates. The degeneracy of the three-dimensional harmonic oscillator energies with principal quantum number n is lifted due to the anharmonicity in the radial potential. The (2l + 1)-fold degeneracy of the angular momentum states is also weakly lifted, due to the angular dependence in the potential. We calculate the fundamental frequency to range between 96 and 110 cm-1 depending on the electronic structure method used. Error bars of the eigenstate energies were calculated from the GPR and are on the order of ∼±1.5 cm-1. Wavefunctions are also compared by considering their overlap and Hellinger distance to the one-dimensional empirical potential. As with the energies, the two ab initio methods MP2 and RPA@PBE show the best agreement. While MP2 has better agreement than RPA@PBE, due to its higher computational efficiency and comparable performance, we recommend RPA as an alternative electronic structure method of choice to MP2 for these systems.
Collapse
Affiliation(s)
- K Panchagnula
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - D Graf
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - F E A Albertani
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| | - A J W Thom
- Yusuf Hamied Department of Chemistry, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
6
|
Iyengar SS, Ricard TC, Zhu X. Reformulation of All ONIOM-Type Molecular Fragmentation Approaches and Many-Body Theories Using Graph-Theory-Based Projection Operators: Applications to Dynamics, Molecular Potential Surfaces, Machine Learning, and Quantum Computing. J Phys Chem A 2024; 128:466-478. [PMID: 38180503 DOI: 10.1021/acs.jpca.3c05630] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/06/2024]
Abstract
We present a graph-theory-based reformulation of all ONIOM-based molecular fragmentation methods. We discuss applications to (a) accurate post-Hartree-Fock AIMD that can be conducted at DFT cost for medium-sized systems, (b) hybrid DFT condensed-phase studies at the cost of pure density functionals, (c) reduced cost on-the-fly large basis gas-phase AIMD and condensed-phase studies, (d) post-Hartree-Fock-level potential surfaces at DFT cost to obtain quantum nuclear effects, and (e) novel transfer machine learning protocols derived from these measures. Additionally, in previous work, the unifying strategy discussed here has been used to construct new quantum computing algorithms. Thus, we conclude that this reformulation is robust and accurate.
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
| | - 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
| | - Xiao Zhu
- 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
|
7
|
Daboe M, Parlak C, Direm A, Alver Ö, Ramasami P. Interaction between escitalopram and ibuprofen or paracetamol: DFT and molecular docking on the drug-drug interactions. J Biomol Struct Dyn 2024; 42:672-686. [PMID: 37042928 DOI: 10.1080/07391102.2023.2195004] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2022] [Accepted: 03/17/2023] [Indexed: 04/13/2023]
Abstract
A large number of drugs are introduced each year to treat different diseases. Most of the time, patients suffer from more than one health problem which makes it necessary to take multiple drugs. When drugs are combined, the problem of drug-drug interaction becomes relevant. In this work, we studied the drug-drug interaction between escitalopram and ibuprofen or paracetamol using density functional theory and quantum theory of atoms in molecules. The results suggest that following the interactions, the activity of drugs changes according to site of interaction. Most reactive and most stable interactions would be preferable for the purpose of use. The in silico drug-likeness studies show that escitalopram and paracetamol couple is more bioavailable than escitalopram and ibuprofen couple. Moreover, in order to gain additional insights into the mentioned drugs' interactions, the drugs were docked separately and jointly against the potential targets for antidepressants and NSAIDs, namely 6HIS and 2PXX. The molecular docking results showed a potential improvement of the effectiveness of the drugs after combining by forming hydrogen bonds, hydrophobic contacts and π…π stacking.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Musa Daboe
- Department of Physics, Science Faculty, Ege University, Izmir, Turkey
| | - Cemal Parlak
- Department of Physics, Science Faculty, Ege University, Izmir, Turkey
| | - Amani Direm
- Department of Matter Sciences, Faculty of Sciences and Technology, Abbes Laghrour University, Khenchela, Algeria
- Laboratory of Structures, Properties and Interatomic Interactions LASPI2A, Faculty of Sciences and Technology, Abbes Laghrour University, Khenchela, Algeria
| | - Özgür Alver
- Department of Physics, Science Faculty, Eskişehir Technical University, Eskişehir, Turkey
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, Mauritius
- Department of Chemistry College of Science, Engineering and Technology, University of South Africa, Pretoria, South Africa
| |
Collapse
|
8
|
Ricard TC, Zhu X, Iyengar SS. Capturing Weak Interactions in Surface Adsorbate Systems at Coupled Cluster Accuracy: A Graph-Theoretic Molecular Fragmentation Approach Improved through Machine Learning. J Chem Theory Comput 2023. [PMID: 38019639 DOI: 10.1021/acs.jctc.3c00955] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2023]
Abstract
The accurate and efficient study of the interactions of organic matter with the surface of water is critical to a wide range of applications. For example, environmental studies have found that acidic polyfluorinated alkyl substances, especially perfluorooctanoic acid (PFOA), have spread throughout the environment and bioaccumulate into human populations residing near contaminated watersheds, leading to many systemic maladies. Thus, the study of the interactions of PFOA with water surfaces became important for the mitigation of their activity as pollutants and threats to public health. However, theoretical study of the interactions of such organic adsorbates on the surface of water, and their bulk concerted properties, often necessitates the use of ab initio methods to properly incorporate the long-range electronic properties that govern these extended systems. Notable theoretical treatments of "on-water" reactions thus far have employed hybrid DFT and semilocal DFT, but the interactions involved are weak interactions that may be best described using post-Hartree-Fock theory. Here, we aim to demonstrate the utility of a graph-theoretic approach to molecular fragmentation that accurately captures the critical "weak" interactions while maintaining an efficient ab initio treatment of the long-range periodic interactions that underpin the physics of extended systems. We apply this graph-theoretical treatment to study PFOA on the surface of water as a model system for the study of weak interactions seen in the wide range of surface interactions and reactions. The approach divides a system into a set of vertices, that are then connected through edges, faces, and higher order graph theoretic objects known as simplexes, to represent a collection of locally interacting subsystems. These subsystems are then used to construct ab initio molecular dynamics simulations and for computing multidimensional potential energy surfaces. To further improve the computational efficiency of our graph theoretic fragmentation method, we use a recently developed transfer learning protocol to construct the full system potential energy from a family of neural networks each designed to accurately model the behavior of individual simplexes. We use a unique multidimensional clustering algorithm, based on the k-means clustering methodology, to define our training space for each separate simplex. These models are used to extrapolate the energies for molecular dynamics trajectories at PFOA water interfaces, at less than one-tenth the cost as compared to a regular molecular fragmentation-based dynamics calculation with excellent agreement with couple cluster level of full system potential energies.
Collapse
Affiliation(s)
- Timothy C Ricard
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Xiao Zhu
- 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
|
9
|
Wang Y, Guo Y, Neese F, Valeev EF, Li W, Li S. Cluster-in-Molecule Approach with Explicitly Correlated Methods for Large Molecules. J Chem Theory Comput 2023; 19:8076-8089. [PMID: 37920973 DOI: 10.1021/acs.jctc.3c00627] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/04/2023]
Abstract
In this article, we present a series of explicitly correlated local correlation methods developed under the cluster-in-molecule (CIM) framework, including explicitly correlated second-order Møller-Plesset perturbation (MP2), coupled-cluster singles and doubles (CCSD), domain-based local pair natural orbital CCSD (DLPNO-CCSD), and DLPNO-CCSD with perturbative triples (DLPNO-CCSD(T)). In these methods, F12 correction is decomposed into contributions from each occupied local molecular orbital and then evaluated independently in a given cluster, which consists of a subset of localized orbitals. These newly developed methods allow F12 calculations of large molecules (up to 145 atoms for quasi-one-dimensional systems) on a single node. We use these methods to investigate the relative stability between extended and folded alkane C30H62, the relative stability of four secondary structures of a polyglycine Ace(Gly)10NH2, and the binding energies of two host-guest complexes. The results demonstrate that the combination of CIM with F12 methods is a promising way to investigate large molecules with small basis set errors.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong 266237, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, New Cornerstone Science Laboratory, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing 210023, P. R. China
| |
Collapse
|
10
|
Drontschenko V, Bangerter FH, Ochsenfeld C. Analytical Second-Order Properties for the Random Phase Approximation: Nuclear Magnetic Resonance Shieldings. J Chem Theory Comput 2023; 19:7542-7554. [PMID: 37863033 DOI: 10.1021/acs.jctc.3c00542] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/22/2023]
Abstract
A method for the analytical computation of nuclear magnetic resonance (NMR) shieldings within the direct random phase approximation (RPA) is presented. As a starting point, we use the RPA ground-state energy expression within the resolution-of-the-identity approximation in the atomic-orbital formalism. As has been shown in a recent benchmark study using numerical second derivatives [Glasbrenner, M. J. Chem. Theory Comput. 2022, 18, 192], RPA based on a Hartree-Fock reference shows accuracies comparable to coupled cluster singles and doubles (CCSD) for NMR chemical shieldings. Together with the much lower computational cost of RPA, it has emerged as an accurate method for the computation of NMR shieldings. Therefore, we aim to extend the applicability of RPA NMR to larger systems by introducing analytical second-order derivatives, making it a viable method for the accurate and efficient computation of NMR chemical shieldings.
Collapse
Affiliation(s)
- Viktoria Drontschenko
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Felix H Bangerter
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), D-81377 Munich, Germany
- Max Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| |
Collapse
|
11
|
Wang Z, Aldossary A, Shi T, Liu Y, Li XS, Head-Gordon M. Local Second-Order Møller-Plesset Theory with a Single Threshold Using Orthogonal Virtual Orbitals: Theory, Implementation, and Assessment. J Chem Theory Comput 2023; 19:7577-7591. [PMID: 37877899 DOI: 10.1021/acs.jctc.3c00744] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2023]
Abstract
It has long been clear that electron correlation methods exhibit unphysical compute scalings with molecular size, which has motivated the development of local correlation methods to discard effectively zero contributions in a controlled way to yield an approximate correlation energy. The ideal local correlation method should have a single numerical threshold that controls the dropping of terms with the ability to have that threshold set small enough so that the correlation energy is reproduced to enough significant figures such that the result is chemically identical. This work reports such a method for the second-order Møller-Plesset (MP2) theory. The theory, implementation, and testing of this local MP2 theory are reported. Thresholds ranging from 10-5 to 10-8 and basis sets ranging from split valence plus polarization through to quadruple-ζ are assessed for local MP2 calculations on a range of molecules, including linear chains and molecules with two- and three-dimensional character. The implementation is shared memory parallel via OpenMP and yields roughly 50% parallel efficiency with 16 cores for a large job. Considerable efforts were made to minimize memory demands, which increased as thresholds were tightened. A variety of relative energy calculations are presented as a function of threshold to provide some guidance to users on how to obtain adequate precision at a low compute cost. It is particularly clear that derivative properties require tighter thresholds in order to achieve an adequate precision.
Collapse
Affiliation(s)
- Zhenling Wang
- Department of Chemistry, University of California, Berkeley, California 94720, United States
- Chemical Sciences Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Abdulrahman Aldossary
- Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Tianyi Shi
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Yang Liu
- Applied Mathematics and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Xiaoye S Li
- Applied Mathematics and 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
|
12
|
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
|
13
|
Werner HJ, Hansen A. Accurate Calculation of Isomerization and Conformational Energies of Larger Molecules Using Explicitly Correlated Local Coupled Cluster Methods in Molpro and ORCA. J Chem Theory Comput 2023; 19:7007-7030. [PMID: 37486154 DOI: 10.1021/acs.jctc.3c00270] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/25/2023]
Abstract
An overview of the approximations in the explicitly correlated local coupled cluster methods PNO-LCCSD(T)-F12 in Molpro and DLPNO-CCSD(T)F12 in ORCA is given. Options to select the domains of projected atomic orbitals (PAOs), pair natural orbitals (PNOs), and triples natural orbitals (TNOs) in both programs are described and compared in detail. The two programs are applied to compute isomerization and conformational energies of the ISOL24 and ACONFL test sets, where the former is part of the GMTKN55 benchmark suite. Thorough studies of basis set effects are presented for selected systems. These revealed large intramolecular basis set superposition effects that make it practically impossible to reliably determine the complete basis set (CBS) limits without including explicitly correlated terms. The latter strongly reduce the basis set dependence and at the same time also errors caused by the local domain approximations. On the basis of these studies, the PNO-LCCSD(T)-F12 method is applied to determine new reference energies for the above-mentioned benchmark sets. We are confident that our results should agree within a few tenths of a kcal mol-1 with the (unknown) CCSD(T)/CBS values, which therefore allowed us to define computational settings for accurate explicitly correlated local coupled cluster methods with moderate computational effort. With these protocols, especially PNO-LCCSD(T)-F12b/AVTZ', reliable reference values for comprehensive benchmark sets can be generated efficiently. This can significantly advance the development and evaluation of the performance of approximate electronic structure methods, especially improved density functional approximations or machine learning approaches.
Collapse
Affiliation(s)
- Hans-Joachim Werner
- Institut für Theoretische Chemie, Universität Stuttgart, Pfaffenwaldring 55, D-70569 Stuttgart, Germany
| | - Andreas Hansen
- Mulliken Center for Theoretical Chemistry, Universität Bonn, Beringstrasse 4, D-53115 Bonn, Germany
| |
Collapse
|
14
|
Di Felice R, Mayes ML, Richard RM, Williams-Young DB, Chan GKL, de Jong WA, Govind N, Head-Gordon M, Hermes MR, Kowalski K, Li X, Lischka H, Mueller KT, Mutlu E, Niklasson AMN, Pederson MR, Peng B, Shepard R, Valeev EF, van Schilfgaarde M, Vlaisavljevich B, Windus TL, Xantheas SS, Zhang X, Zimmerman PM. A Perspective on Sustainable Computational Chemistry Software Development and Integration. J Chem Theory Comput 2023; 19:7056-7076. [PMID: 37769271 PMCID: PMC10601486 DOI: 10.1021/acs.jctc.3c00419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2023] [Indexed: 09/30/2023]
Abstract
The power of quantum chemistry to predict the ground and excited state properties of complex chemical systems has driven the development of computational quantum chemistry software, integrating advances in theory, applied mathematics, and computer science. The emergence of new computational paradigms associated with exascale technologies also poses significant challenges that require a flexible forward strategy to take full advantage of existing and forthcoming computational resources. In this context, the sustainability and interoperability of computational chemistry software development are among the most pressing issues. In this perspective, we discuss software infrastructure needs and investments with an eye to fully utilize exascale resources and provide unique computational tools for next-generation science problems and scientific discoveries.
Collapse
Affiliation(s)
- Rosa Di Felice
- Departments
of Physics and Astronomy and Quantitative and Computational Biology, University of Southern California, Los Angeles, California 90089, United States
- CNR-NANO
Modena, Modena 41125, Italy
| | - Maricris L. Mayes
- Department
of Chemistry and Biochemistry, University
of Massachusetts Dartmouth, North Dartmouth, Massachusetts 02747, United States
| | | | | | - Garnet Kin-Lic Chan
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Wibe A. de Jong
- Lawrence
Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Niranjan Govind
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Martin Head-Gordon
- Pitzer Center
for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, United States
| | - Matthew R. Hermes
- Department
of Chemistry, Chicago Center for Theoretical Chemistry, University of Chicago, Chicago, Illinois 60637, United States
| | - Karol Kowalski
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Xiaosong Li
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
| | - Hans Lischka
- Department
of Chemistry and Biochemistry, Texas Tech
University, Lubbock, Texas 79409, United States
| | - Karl T. Mueller
- Physical
and Computational Sciences Directorate, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Erdal Mutlu
- Advanced
Computing, Mathematics, and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Anders M. N. Niklasson
- Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, United States
| | - Mark R. Pederson
- Department
of Physics, The University of Texas at El
Paso, El Paso, Texas 79968, United States
| | - Bo Peng
- Physical
Sciences Division, Pacific Northwest National
Laboratory, Richland, Washington 99354, United States
| | - Ron Shepard
- Chemical
Sciences and Engineering Division, Argonne
National Laboratory, Lemont, Illinois 60439, United States
| | - Edward F. Valeev
- Department
of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | | | - Bess Vlaisavljevich
- Department
of Chemistry, University of South Dakota, Vermillion, South Dakota 57069, United States
| | - Theresa L. Windus
- Department
of Chemistry, Iowa State University and
Ames Laboratory, Ames, Iowa 50011, United States
| | - Sotiris S. Xantheas
- Department
of Chemistry, University of Washington, Seattle, Washington 98195, United States
- Advanced
Computing, Mathematics and Data Division, Pacific Northwest National Laboratory, Richland, Washington 99354, United States
| | - Xing Zhang
- Division
of Chemistry and Chemical Engineering, California
Institute of Technology, Pasadena, California 91125, United States
| | - Paul M. Zimmerman
- Department
of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| |
Collapse
|
15
|
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
|
16
|
Demel O, Lecours MJ, Nooijen M. Further investigations into a Laplace MP2 method using range separated Coulomb potential and orbital selective virtuals: Multipole correction, OSV extrapolation, and critical assessment. J Chem Phys 2023; 158:114120. [PMID: 36948803 DOI: 10.1063/5.0135113] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/03/2023] Open
Abstract
We report further investigations to aid the development of a Laplace MP2 (second-order Møller Plesset) method with a range separated Coulomb potential partitioned into short- and long-range parts. The implementation of the method extensively uses sparse matrix algebra, density fitting techniques for the short-range part, and a Fourier transformation in spherical coordinates for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas virtual space is described by orbital specific virtual orbitals (OSVs) associated with localized molecular orbitals. The Fourier transform is deficient for very large distances between localized occupied orbitals, and a multipole expansion for widely separated pairs is introduced for the direct MP2 contribution, which is applicable also to non-Coulombic potentials that do not satisfy the Laplace equation. For the exchange contribution, an efficient screening of contributing localized occupied pairs is employed, which is discussed more completely here. To mitigate errors due to the truncation of OSVs, a simple and efficient extrapolation procedure is used to obtain results close to MP2 for the full basis set of atomic orbitals Using a suitable set of default parameters, the accuracy of the approach is demonstrated. The current implementation of the approach is not very efficient, and the aim of this paper is to introduce and critically discuss ideas that can have more general applicability beyond MP2 calculations for large molecules.
Collapse
Affiliation(s)
- Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Michael J Lecours
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Marcel Nooijen
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
17
|
Spadetto E, Philipsen PHT, Förster A, Visscher L. Toward Pair Atomic Density Fitting for Correlation Energies with Benchmark Accuracy. J Chem Theory Comput 2023; 19:1499-1516. [PMID: 36787494 PMCID: PMC10018742 DOI: 10.1021/acs.jctc.2c01201] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/16/2023]
Abstract
Pair atomic density fitting (PADF) has been identified as a promising strategy to reduce the scaling with system size of quantum chemical methods for the calculation of the correlation energy like the direct random-phase approximation (RPA) or second-order Møller-Plesset perturbation theory (MP2). PADF can however introduce large errors in correlation energies as the two-electron interaction energy is not guaranteed to be bounded from below. This issue can be partially alleviated by using very large fit sets, but this comes at the price of reduced efficiency and having to deal with near-linear dependencies in the fit set. One posibility is to use global density fitting (DF), but in this work, we introduce an alternative methodology to overcome this problem that preserves the intrinsically favorable scaling of PADF. We first regularize the Fock matrix by projecting out parts of the basis set which gives rise to orbital products that are hard to describe by PADF. After having thus obtained a reliable self-consistent field solution, we then also apply this projector to the orbital coefficient matrix to improve the precision of PADF-MP2 and PADF-RPA. We systematically assess the accuracy of this new approach in a numerical atomic orbital framework using Slater type orbitals (STO) and correlation consistent Gaussian type basis sets up to quintuple-ζ quality for systems with more than 200 atoms. For the small and medium systems in the S66 database we show the maximum deviation of PADF-MP2 and PADF-RPA relative correlation energies to DF-MP2 and DF-RPA reference results to be 0.07 and 0.14 kcal/mol, respectively. When the new projector method is used, the errors only slightly increase for large molecules and also when moderately sized fit sets are used the resulting errors are well under control. Finally, we demonstrate the computational efficiency of our algorithm by calculating the interaction energies of large, non-covalently bound complexes with more than 1000 atoms and 20000 atomic orbitals at the RPA@PBE/CC-pVTZ level of theory.
Collapse
Affiliation(s)
- Edoardo Spadetto
- Software for Chemistry and Materials NV, NL-1081HV Amsterdam, The Netherlands
| | | | - Arno Förster
- Software for Chemistry and Materials NV, NL-1081HV Amsterdam, The Netherlands.,Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| | - Lucas Visscher
- Theoretical Chemistry, Vrije Universiteit, De Boelelaan 1083, NL-1081 HV Amsterdam, The Netherlands
| |
Collapse
|
18
|
Kaya MO, Demirci T, Ozdemir O, Calisir U, Sonmez F, Arslan M. Synthesis, inhibition effects, molecular docking and theoretical studies as Paraoxonase 1 (PON1) inhibitors of novel 1,4-dihydropyridine substituted sulfonamide derivatives. Med Chem Res 2023. [DOI: 10.1007/s00044-023-03029-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/05/2023]
|
19
|
Wang Z, Aldossary A, Head-Gordon M. Sparsity of the electron repulsion integral tensor using different localized virtual orbital representations in local second-order Møller-Plesset theory. J Chem Phys 2023; 158:064105. [PMID: 36792513 DOI: 10.1063/5.0134764] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/20/2023] Open
Abstract
Utilizing localized orbitals, local correlation theory can reduce the unphysically high system-size scaling of post-Hartree-Fock (post-HF) methods to linear scaling in insulating molecules. The sparsity of the four-index electron repulsion integral (ERI) tensor is central to achieving this reduction. For second-order Møller-Plesset theory (MP2), one of the simplest post-HF methods, only the (ia|jb) ERIs are needed, coupling occupied orbitals i, j and virtuals a, b. In this paper, we compare the numerical sparsity (called the "ragged list") and two other approaches revealing the low-rank sparsity of the ERI. The ragged list requires only one set of (localized) virtual orbitals, and we find that the orthogonal valence virtual-hard virtual set of virtuals originally proposed by Subotnik et al. gives the sparsest ERI tensor. To further compress the ERI tensor, the pair natural orbital (PNO) type representation uses different sets of virtual orbitals for different occupied orbital pairs, while the occupied-specific virtual (OSV) approach uses different virtuals for each occupied orbital. Our results indicate that while the low-rank PNO representation achieves significant rank reduction, it also requires more memory than the ragged list. The OSV approach requires similar memory to that of the ragged list, but it involves greater algorithmic complexity. An approximation (called the "fixed sparsity pattern") for solving the local MP2 equations using the numerically sparse ERI tensor is proposed and tested to be sufficiently accurate and to have highly controllable error. A low-scaling local MP2 algorithm based on the ragged list and the fixed sparsity pattern is therefore promising.
Collapse
Affiliation(s)
- Zhenling Wang
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Abdulrahman Aldossary
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| | - Martin Head-Gordon
- Pitzer Center for Theoretical Chemistry, Department of Chemistry, University of California, Berkeley, California 94720, USA
| |
Collapse
|
20
|
Convergence of Møller–Plesset perturbation theory for excited reference states. ADVANCES IN QUANTUM CHEMISTRY 2023. [DOI: 10.1016/bs.aiq.2023.01.012] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
|
21
|
Drontschenko V, Graf D, Laqua H, Ochsenfeld C. Efficient Method for the Computation of Frozen-Core Nuclear Gradients within the Random Phase Approximation. J Chem Theory Comput 2022; 18:7359-7372. [PMID: 36331398 DOI: 10.1021/acs.jctc.2c00774] [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/06/2022]
Abstract
A method for the evaluation of analytical frozen-core gradients within the random phase approximation is presented. We outline an efficient way to evaluate the response of the density of active electrons arising only when introducing the frozen-core approximation and constituting the main difficulty, together with the response of the standard Kohn-Sham density. The general framework allows to extend the outlined procedure to related electron correlation methods in the atomic orbital basis that require the evaluation of density responses, such as second-order Møller-Plesset perturbation theory or coupled cluster variants. By using Cholesky decomposed densities─which reintroduce the occupied index in the time-determining steps─we are able to achieve speedups of 20-30% (depending on the size of the basis set) by using the frozen-core approximation, which is of similar magnitude as for molecular orbital formulations. We further show that the errors introduced by the frozen-core approximation are practically insignificant for molecular geometries.
Collapse
Affiliation(s)
- Viktoria Drontschenko
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany
| | - Daniel Graf
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany
| | - Henryk Laqua
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), 81377 Munich, Germany.,Max Planck Institute for Solid State Research, D-70569 Stuttgart, Germany
| |
Collapse
|
22
|
İzzettinYILMAZER M. Theoretical Study of the Adsorption Mechanism of Ioversol Contrast Agent on the Single-Walled Carbon Nanotube Surface. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.134504] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
|
23
|
Lesiuk M. When Gold Is Not Enough: Platinum Standard of Quantum Chemistry with N7 Cost. J Chem Theory Comput 2022; 18:6537-6556. [PMID: 36314739 PMCID: PMC9648181 DOI: 10.1021/acs.jctc.2c00460] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
![]()
In this paper, we extend the rank-reduced coupled-cluster
formalism
to the calculation of non-iterative energy corrections due to quadruple
excitations. There are two major components of the proposed formalism.
The first is an approximate compression of the quadruple excitation
amplitudes using the Tucker format. The second is a modified functional
used for the evaluation of the corrections which gives exactly the
same results for the exact amplitudes, but is less susceptible to
errors resulting from the aforementioned compression. We show, both
theoretically and numerically, that the computational cost of the
proposed method scales as the seventh power of the system size. Using
reference results for a set of small molecules, the method is calibrated
to deliver relative accuracy of a few percent in energy corrections.
To illustrate the potential of the theory, we calculate the isomerization
energy of ortho/meta benzyne (C6H4) and the barrier height for the Cope rearrangement
in bullvalene (C10H10). The method retains a
near-black-box nature of the conventional coupled-cluster formalism
and depends on only one additional parameter that controls the accuracy.
Collapse
Affiliation(s)
- Michał Lesiuk
- Faculty of Chemistry, University of Warsaw, Pasteura 1, Warsaw 02-093, Poland
| |
Collapse
|
24
|
Wang Y, Ni Z, Neese F, Li W, Guo Y, Li S. Cluster-in-Molecule Method Combined with the Domain-Based Local Pair Natural Orbital Approach for Electron Correlation Calculations of Periodic Systems. J Chem Theory Comput 2022; 18:6510-6521. [PMID: 36240189 DOI: 10.1021/acs.jctc.2c00412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The cluster-in-molecule (CIM) method was extended to systems with periodic boundary conditions (PBCs) in a previous work (PBC-CIM) [J. Chem. Theory Comput.2019, 15, 2933], which is able to compute the electronic structures of periodic systems at second-order Møller-Plesset perturbation theory (MP2) and coupled cluster singles and doubles (CCSD) levels. However, the high computational costs of CCSD with respect to the size of clusters limit the usage of PBC-CIM to crystals with small or medium unit cells. In this work, we further develop the PBC-CIM method by employing the domain-based local pair natural orbital (DLPNO) methods for the electron correlation calculations of clusters to reduce the computational costs. The combined approach allows CCSD with perturbative triples, denoted as CCSD(T), to be computationally available for accurate descriptions of periodic systems. The distant-pair correction is also implemented to improve the accuracy of PBC-CIM. As in the molecular cases, the distant pair correction significantly improves the accuracy of various PBC-CIM methods with few additional costs. The PBC-CIM-DLPNO-CCSD(T) approach has been applied to investigate the optimized lattice parameter of the cubic LiCl crystal and two adsorption problems (CO on the NaCl(100) surface and H2O on the h-BN surface). The results show that the CIM-DLPNO-CCSD(T) method offers accurate and efficient descriptions for the studied systems. Another application to the cohesive energy of the acetic acid crystal reveals that large basis sets are necessary for reliable calculations on the cohesive energies of molecular crystals.
Collapse
Affiliation(s)
- Yuqi Wang
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Zhigang Ni
- College of Material, Chemistry and Chemical Engineering, Key Laboratory of Organosilicon Chemistry and Material Technology of Ministry of Education, Hangzhou Normal University, Hangzhou311121, P. R. China
| | - Frank Neese
- Max Planck Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, Mülheim an der RuhrD-45470, Germany
| | - Wei Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| | - Yang Guo
- Qingdao Institute for Theoretical and Computational Sciences, Shandong University, Qingdao, Shandong266237, P. R. China
| | - Shuhua Li
- School of Chemistry and Chemical Engineering, Key Laboratory of Mesoscopic Chemistry of MOE, Institute of Theoretical and Computational Chemistry, Nanjing University, Nanjing210023, P. R. China
| |
Collapse
|
25
|
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
|
26
|
Kadri S, Direm A, Athmani H, El Bali B, Parlak C, Hebbachi R. Dual inhibition of S. aureus TyrRS and S. aureus gyrase by two 4-amino-4′-acetyldiphenyl sulfide-based Schiff bases: Structural features, DFT study, Hirshfeld surface analysis and molecular docking. INORG CHEM COMMUN 2022. [DOI: 10.1016/j.inoche.2022.109779] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
|
27
|
Whitten JL. Correction of Residual Errors in Configuration Interaction Electronic Structure Calculations. J Chem Phys 2022; 157:124101. [DOI: 10.1063/5.0098793] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Methods for correcting residual energy errors of configuration interaction (CI) calculations of molecules and other electronic systems are discussed based on the assumption that the energy defect can be mapped onto atomic regions. The methods do not consider the detailed nature of excitations, but instead define a defect energy per electron that that is unique to a specific atom. Defect energy contributions are determined from calculations on diatomic and hydride molecules and then applied to other systems. Calculated energies are compared with experimental thermodynamic and spectroscopic data for a set of forty-one mainly organic molecules representing a wide range of bonding environments. The most stringent test is based on a severely truncated virtual space in which higher spherical harmonic basis functions are removed. The errors of the initial CI calculations are large, but in each case, including defect corrections brings calculated CI energies into agreement with experimental values. The method is also applied to a NIST compilation of coupled-cluster calculations that employ a larger basis set and no truncation of the virtual space. The corrections show excellent consistency with total energies in very good agreement with experimental values. An extension of the method is applied to dmsn states of Sc, Ti, V, Mn, Cr, Fe, Co, Ni and Cu, significantly improving the agreement of calculated transition energies with spectroscopic values.
Collapse
Affiliation(s)
- Jerry L Whitten
- North Carolina State University, North Carolina State University, United States of America
| |
Collapse
|
28
|
Alkan L, Alver Ö, Asley ST, Kulan E, Parlak C. A convenient one-pot synthesis, theoretical studies, and NMR-based conformational analysis of nitroarylidene dimalonates. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2022.133848] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/16/2022]
|
29
|
Parlak C, Alver Ö, Bağlayan Ö, Ramasami P. Theoretical insights of the drug-drug interaction between favipiravir and ibuprofen: a DFT, QTAIM and drug-likeness investigation. J Biomol Struct Dyn 2022:1-8. [DOI: 10.1080/07391102.2022.2066022] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/18/2022]
Affiliation(s)
- Cemal Parlak
- Department of Physics, Science Faculty, Ege University, Izmir, Turkey
| | - Özgür Alver
- Department of Physics, Science Faculty, Eskisehir Technical University, Eskisehir, Turkey
| | - Özge Bağlayan
- Department of Physics, Science Faculty, Eskisehir Technical University, Eskisehir, Turkey
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, Mauritius
- Department of Chemistry, College of Science, Engineering and Technology, University of South Africa, Pretoria, South Africa
| |
Collapse
|
30
|
Yan W, Xu X. Accurate Prediction of Nuclear Magnetic Resonance Parameters via the XYG3 Type of Doubly Hybrid Density Functionals. J Chem Theory Comput 2022; 18:2931-2946. [PMID: 35467852 DOI: 10.1021/acs.jctc.2c00055] [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
Nuclear magnetic resonance (NMR) spectroscopy is one of the most powerful and versatile tools in elucidating molecular structures. To eliminate ambiguities of experimental assignments, accurate calculations of NMR spectra are of great importance. Here, a method for theoretical evaluation of the NMR shielding constants by analytic derivatives using gauge including atomic orbitals (GIAO) has been implemented for the XYG3 type of doubly hybrid density functionals (xDH), namely, the GIAO-xDH method. Benchmark calculations on shielding constants and chemical shifts demonstrate the remarkable accuracy of the GIAO-xDH method, compared to the accurate CCSD(T) references. It is shown here that the XYGJ-OS functional is able to give a mean absolute deviation (MAD) of ∼3.0 ppm in the calculated shielding constants for 13C, 15N, 17O, 19F, while both XYGJ-OS and xDH-PBE0 functionals are able to provide a satisfactory estimation of chemical shifts with MADs of ∼0.03 and 1.0 ppm for 1H and 13C, respectively. The basis set influence upon the method has been examined and a computational scheme considering both accuracy and efficiency has been proposed and tested to predict the experimental 13C chemical shifts of five medium-sized natural product molecules, yielding a MAD of ∼1.0 ppm, which demonstrates the practical feasibility of the GIAO-xDH method.
Collapse
Affiliation(s)
- Wenjie Yan
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, 200433, China
| | - Xin Xu
- Collaborative Innovation Center of Chemistry for Energy Materials, Shanghai Key Laboratory of Molecular Catalysis and Innovative Materials, Ministry of Education Key Laboratory of Computational Physical Sciences, Department of Chemistry, Fudan University, Shanghai, 200433, China
| |
Collapse
|
31
|
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
|
32
|
Toktas U, Sarikahya NB, Parlak C, Ozturk I, Kayalar H. A new iridoid skeleton from Galium asparagifolium and biological activity studies. J Mol Struct 2022. [DOI: 10.1016/j.molstruc.2021.131693] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/20/2022]
|
33
|
Lesiuk M. Quintic-scaling rank-reduced coupled cluster theory with single and double excitations. J Chem Phys 2022; 156:064103. [DOI: 10.1063/5.0071916] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
34
|
Glasbrenner M, Vogler S, Ochsenfeld C. Efficient low-scaling computation of NMR shieldings at the second-order Møller-Plesset perturbation theory level with Cholesky-decomposed densities and an attenuated Coulomb metric. J Chem Phys 2021; 155:224107. [PMID: 34911319 DOI: 10.1063/5.0069956] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A method for the computation of nuclear magnetic resonance (NMR) shieldings with second-order Møller-Plesset perturbation theory (MP2) is presented which allows to efficiently compute the entire set of shieldings for a given molecular structure. The equations are derived using Laplace-transformed atomic orbital second-order Møller-Plesset perturbation theory as a starting point. The Z-vector approach is employed for minimizing the number of coupled-perturbed self-consistent-field equations that need to be solved. In addition, the method uses the resolution-of-the-identity approximation with an attenuated Coulomb metric and Cholesky decomposition of pseudo-density matrices. The sparsity in the three-center integrals is exploited with sparse linear algebra approaches, leading to reduced computational cost and memory demands. Test calculations show that the deviations from NMR shifts obtained with canonical MP2 are small if appropriate thresholds are used. The performance of the method is illustrated in calculations on DNA strands and on glycine chains with up to 283 atoms and 2864 basis functions.
Collapse
Affiliation(s)
- Michael Glasbrenner
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| | - Sigurd Vogler
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| | - Christian Ochsenfeld
- Chair of Theoretical Chemistry, Department of Chemistry, University of Munich (LMU), Butenandtstr. 7, D-81377 Munich, Germany
| |
Collapse
|
35
|
Altun A, Ghosh S, Riplinger C, Neese F, Bistoni G. Addressing the System-Size Dependence of the Local Approximation Error in Coupled-Cluster Calculations. J Phys Chem A 2021; 125:9932-9939. [PMID: 34730360 PMCID: PMC8607505 DOI: 10.1021/acs.jpca.1c09106] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
Over the last two decades, the local approximation has been successfully used to extend the range of applicability of the "gold standard" singles and doubles coupled-cluster method with perturbative triples CCSD(T) to systems with hundreds of atoms. The local approximation error grows in absolute value with the increasing system size, i.e., by increasing the number of electron pairs in the system. In this study, we demonstrate that the recently introduced two-point extrapolation scheme for approaching the complete pair natural orbital (PNOs) space limit in domain-based pair natural orbital CCSD(T) calculations drastically reduces the dependence of the error on the system size, thus opening up unprecedented opportunities for the calculation of benchmark quality relative energies for large systems.
Collapse
Affiliation(s)
- Ahmet Altun
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Soumen Ghosh
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | | | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| | - Giovanni Bistoni
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, D-45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
36
|
Liang Q, Yang J. Third-Order Many-Body Expansion of OSV-MP2 Wave Function for Low-Order Scaling Analytical Gradient Computation. J Chem Theory Comput 2021; 17:6841-6860. [PMID: 34704757 DOI: 10.1021/acs.jctc.1c00581] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
Abstract
We present a many-body expansion (MBE) formulation and implementation for efficient computation of analytical energy gradients from the orbital-specific-virtual second-order Møllet-Plesset perturbation theory (OSV-MP2) based on our earlier work (Zhou et al. J. Chem. Theory Comput. 2020, 16, 196-210). The third-order MBE(3) expansion of OSV-MP2 amplitudes and density matrices was developed to adopt the orbital-specific clustering and long-range termination schemes, which avoids term-by-term differentiations of the MBE energy bodies. We achieve better efficiency by exploiting the algorithmic sparsity that allows us to prune out insignificant fitting integrals and OSV relaxations. With these approximations, the present implementation is benchmarked on a range of molecules that show an economic scaling in the linear and quadratic regimes for computing MBE(3)-OSV-MP2 amplitude and gradient equations, respectively, and yields normal accuracy comparable to the original OSV-MP2 results. The MPI-3-based parallelism through shared memory one-sided communication is further developed for improving parallel scalability and memory accessibility by sorting the MBE(3) orbital clusters into independent tasks that are distributed on multiple processes across many nodes, supporting both global and local data locations in which selected MBE(3)-OSV-MP2 intermediates of different sizes are distinguished and accordingly placed. The accuracy and efficiency level of our MBE(3)-OSV-MP2 analytical gradient implementation is finally illustrated in two applications: we show that the subtle coordination structure differences of mechanically interlocked Cu-catenane complexes can be distinguished when tuning ligand lengths; and the porphycene molecular dynamics reveals the emergence of the vibrational signature arising from softened N-H stretching associated with hydrogen transfer, using an MP2 level of electron correlation and classical nuclei for the first time.
Collapse
Affiliation(s)
- Qiujiang Liang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| | - Jun Yang
- Department of Chemistry, The University of Hong Kong, Hong Kong SAR, P. R. China
| |
Collapse
|
37
|
Kumar A, DeGregorio N, Iyengar SS. Graph-Theory-Based Molecular Fragmentation for Efficient and Accurate Potential Surface Calculations in Multiple Dimensions. J Chem Theory Comput 2021; 17:6671-6690. [PMID: 34623129 DOI: 10.1021/acs.jctc.1c00065] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
We present a multitopology molecular fragmentation approach, based on graph theory, to calculate multidimensional potential energy surfaces in agreement with post-Hartree-Fock levels of theory but at the density functional theory cost. A molecular assembly is coarse-grained into a set of graph-theoretic nodes that are then connected with edges to represent a collection of locally interacting subsystems up to an arbitrary order. Each of the subsystems is treated at two levels of electronic structure theory, the result being used to construct many-body expansions that are embedded within an ONIOM scheme. These expansions converge rapidly with the many-body order (or graphical rank) of subsystems and capture many-body interactions accurately and efficiently. However, multiple graphs, and hence multiple fragmentation topologies, may be defined in molecular configuration space that may arise during conformational sampling or from reactive, bond breaking and bond formation, events. Obtaining the resultant potential surfaces is an exponential scaling proposition, given the number of electronic structure computations needed. We utilize a family of graph-theoretic representations within a variational scheme to obtain multidimensional potential surfaces at a reduced cost. The fast convergence of the graph-theoretic expansion with increasing order of many-body interactions alleviates the exponential scaling cost for computing potential surfaces, with the need to only use molecular fragments that contain a fewer number of quantum nuclear degrees of freedom compared to the full system. This is because the dimensionality of the conformational space sampled by the fragment subsystems is much smaller than the full molecular configurational space. Additionally, we also introduce a multidimensional clustering algorithm, based on physically defined criteria, to reduce the number of energy calculations by orders of magnitude. The molecular systems benchmarked include coupled proton motion in protonated water wires. The potential energy surfaces and multidimensional nuclear eigenstates obtained are shown to be in very good agreement with those from explicit post-Hartree-Fock calculations that become prohibitive as the number of quantum nuclear dimensions grows. The developments here provide a rigorous and efficient alternative to this important chemical physics problem.
Collapse
Affiliation(s)
- Anup Kumar
- Department of Chemistry and Department of Physics, Indiana University, 800 E. Kirkwood Avenue, Bloomington, Indiana 47405, United States
| | - Nicole DeGregorio
- 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
|
38
|
Parlak C, Alver Ö, Ouma CNM, Rhyman L, Ramasami P. Interaction between favipiravir and hydroxychloroquine and their combined drug assessment: in silico investigations. ACTA ACUST UNITED AC 2021; 76:1471-1478. [PMID: 34744292 PMCID: PMC8562770 DOI: 10.1007/s11696-021-01946-8] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/19/2021] [Accepted: 10/22/2021] [Indexed: 12/23/2022]
Abstract
Hydroxychloroquine (HCQ) and favipiravir (FPV) are known to be effective antivirals, and there are reports about their use to fight the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) despite that these are not conclusive. The use of combined drugs is common in drug discovery, and thus, we investigated HCQ and FPV as a combined drug. The density functional theory method was used for the optimization of geometries, spectroscopic analysis and calculation of reactivity parameters. The quantum theory of atoms in molecules was applied to explain the nature of the hydrogen bonds and confirm the higher stability of the combined drug. We also evaluated the absorption, distribution, metabolism and excretion (ADME) parameters to assess their drug actions jointly using SwissADME. The preliminary findings of our theoretical study are promising for further investigations of more potent and selective antiviral drugs.
Collapse
Affiliation(s)
- Cemal Parlak
- Department of Physics, Science Faculty, Ege University, Izmir, 35100 Turkey
| | - Özgür Alver
- Department of Physics, Science Faculty, Eskisehir Technical University, Eskisehir, Turkey
| | - Cecil Naphtaly Moro Ouma
- HySA-Infrastructure CoC, Faculty of Engineering, North-West University, Private Bag X6001, Potchefstroom, 2531 South Africa
| | - Lydia Rhyman
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837 Mauritius.,Centre for Natural Product Research, Department of Chemical Sciences, University of Johannesburg, Doornfontein Campus, Johannesburg, 2028 South Africa
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, 80837 Mauritius.,Centre for Natural Product Research, Department of Chemical Sciences, University of Johannesburg, Doornfontein Campus, Johannesburg, 2028 South Africa
| |
Collapse
|
39
|
Demel O, Lecours MJ, Habrovský R, Nooijen M. Toward Laplace MP2 method using range separated Coulomb potential and orbital selective virtuals. J Chem Phys 2021; 155:154104. [PMID: 34686052 DOI: 10.1063/5.0060099] [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/14/2022] Open
Abstract
We report the development of a new Laplace MP2 (second-order Møller-Plesset) implementation using a range separated Coulomb potential, partitioned into short- and long-range parts. The implementation heavily relies on the use of sparse matrix algebra, density fitting techniques for the short-range Coulomb interactions, while a Fourier transformation in spherical coordinates is used for the long-range part of the potential. Localized molecular orbitals are employed for the occupied space, whereas orbital specific virtual orbitals associated with localized molecular orbitals are obtained from the exchange matrix associated with specific localized occupied orbitals. The range separated potential is crucial to achieve efficient treatment of the direct term in the MP2, while extensive screening is employed to reduce the expense of the exchange contribution in MP2. The focus of this paper is on controllable accuracy and linear scaling of the data entering the algorithm.
Collapse
Affiliation(s)
- Ondřej Demel
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Michael J Lecours
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| | - Richard Habrovský
- J. Heyrovský Institute of Physical Chemistry, v.v.i., Academy of Sciences of the Czech Republic, Dolejškova 3, 18223 Prague 8, Czech Republic
| | - Marcel Nooijen
- University of Waterloo, 200 University Avenue West, Waterloo, Ontario N2L 3G1, Canada
| |
Collapse
|
40
|
Parlak C, Alver Ö, Ouma CNM, Rhyman L, Ramasami P. Can the Antivirals Remdesivir and Favipiravir Work Better Jointly? In Silico Insights. Drug Res (Stuttg) 2021; 72:34-40. [PMID: 34535038 DOI: 10.1055/a-1585-1323] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023]
Abstract
The proverb "Old is gold" is applicable in drug discovery and the proverb "All that Glitters is not Gold" is also appropriate. In the COVID-19 era, there has been a race for drugs to be effective against SARS-CoV-2. There are reports about the uses of Remdesivir and Favipiravir as existing antivirals against virus but none have been conclusive so far. In the attempts for innovations, the combination of drugs is also under trials. Therefore, we used the density functional theory method and quantum theory of atoms in molecules to investigate drug-drug interactions involving Remdesivir and Favipiravir. The computed parameters were related to the antiviral actions of both drugs together. The results indicate enhanced antiviral activity and it will be worthy to consider additional investigations with the combination of these two drugs.
Collapse
Affiliation(s)
- Cemal Parlak
- Department of Physics, Science Faculty, Ege University, Izmir, Turkey
| | - Özgür Alver
- Department of Physics, Science Faculty, Eskisehir Technical University, Eskisehir, Turkey
| | - Cecil N M Ouma
- HySA-Infrastructure CoC, North-West University, Faculty of Engineering, Potchefstroom, South Africa
| | - Lydia Rhyman
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, Mauritius.,Centre for Natural Product Research, Department of Chemical Sciences, University of Johannesburg, Doornfontein Campus, Johannesburg, South Africa
| | - Ponnadurai Ramasami
- Computational Chemistry Group, Department of Chemistry, Faculty of Science, University of Mauritius, Réduit, Mauritius.,Centre for Natural Product Research, Department of Chemical Sciences, University of Johannesburg, Doornfontein Campus, Johannesburg, South Africa
| |
Collapse
|
41
|
Parlak C, Alver Ö, Bağlayan Ö. Quantum mechanical simulation of Molnupiravir drug interaction with Si-doped C60 fullerene. COMPUT THEOR CHEM 2021. [DOI: 10.1016/j.comptc.2021.113336] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/16/2022]
|
42
|
Stoychev GL, Auer AA, Gauss J, Neese F. DLPNO-MP2 second derivatives for the computation of polarizabilities and NMR shieldings. J Chem Phys 2021; 154:164110. [PMID: 33940835 DOI: 10.1063/5.0047125] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
We present a derivation and efficient implementation of the formally complete analytic second derivatives for the domain-based local pair natural orbital second order Møller-Plesset perturbation theory (MP2) method, applicable to electric or magnetic field-response properties but not yet to harmonic frequencies. We also discuss the occurrence and avoidance of numerical instability issues related to singular linear equation systems and near linear dependences in the projected atomic orbital domains. A series of benchmark calculations on medium-sized systems is performed to assess the effect of the local approximation on calculated nuclear magnetic resonance shieldings and the static dipole polarizabilities. Relative deviations from the resolution of the identity-based MP2 (RI-MP2) reference for both properties are below 0.5% with the default truncation thresholds. For large systems, our implementation achieves quadratic effective scaling, is more efficient than RI-MP2 starting at 280 correlated electrons, and is never more than 5-20 times slower than the equivalent Hartree-Fock property calculation. The largest calculation performed here was on the vancomycin molecule with 176 atoms, 542 correlated electrons, and 4700 basis functions and took 3.3 days on 12 central processing unit cores.
Collapse
Affiliation(s)
- Georgi L Stoychev
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Alexander A Auer
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Jürgen Gauss
- Department Chemie, Johannes Gutenberg-Universität Mainz, Duesbergweg 10-14, 55128 Mainz, Germany
| | - Frank Neese
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| |
Collapse
|
43
|
Szabó PB, Csóka J, Kállay M, Nagy PR. Linear-Scaling Open-Shell MP2 Approach: Algorithm, Benchmarks, and Large-Scale Applications. J Chem Theory Comput 2021; 17:2886-2905. [PMID: 33819030 PMCID: PMC8154337 DOI: 10.1021/acs.jctc.1c00093] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
![]()
A linear-scaling
local second-order Møller–Plesset
(MP2) method is presented for high-spin open-shell molecules based
on restricted open-shell (RO) reference functions. The open-shell
local MP2 (LMP2) approach inherits the iteration- and redundancy-free
formulation and the completely integral-direct, OpenMP-parallel, and
memory and disk use economic algorithms of our closed-shell LMP2 implementation.
By utilizing restricted local molecular orbitals for the demanding
integral transformation step and by introducing a novel long-range
spin-polarization approximation, the computational cost of RO-LMP2
approaches that of closed-shell LMP2. Extensive benchmarks were performed
for reactions of radicals, ionization potentials, as well as spin-state
splittings of carbenes and transition-metal complexes. Compared to
the conventional MP2 reference for systems of up to 175 atoms, local
errors of at most 0.1 kcal/mol were found, which are well below the
intrinsic accuracy of MP2. RO-LMP2 computations are presented for
challenging protein models of up to 601 atoms and 11 000 basis
functions, which involve either spin states of a complexed iron ion
or a highly delocalized singly occupied orbital. The corresponding
runtimes of 9–15 h obtained with a single, many-core CPU demonstrate
that MP2, as well as spin-scaled MP2 and double-hybrid density functional
methods, become widely accessible for open-shell systems of unprecedented
size and complexity.
Collapse
Affiliation(s)
- P Bernát Szabó
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - József Csóka
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Mihály Kállay
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| | - Péter R Nagy
- Department of Physical Chemistry and Materials Science, Budapest University of Technology and Economics, P.O. Box 91, H-1521 Budapest, Hungary
| |
Collapse
|
44
|
Fujimori T, Kobayashi M, Taketsugu T. Energy-based automatic determination of buffer region in the divide-and-conquer second-order Møller-Plesset perturbation theory. J Comput Chem 2021; 42:620-629. [PMID: 33534916 PMCID: PMC7986104 DOI: 10.1002/jcc.26486] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2020] [Revised: 12/19/2020] [Accepted: 01/15/2021] [Indexed: 11/21/2022]
Abstract
In the linear‐scaling divide‐and‐conquer (DC) electronic structure method, each subsystem is calculated together with the neighboring buffer region, the size of which affects the energy error introduced by the fragmentation in the DC method. The DC self‐consistent field calculation utilizes a scheme to automatically determine the appropriate buffer region that is as compact as possible for reducing the computational time while maintaining acceptable accuracy (J. Comput. Chem. 2018, 39, 909). To extend the automatic determination scheme of the buffer region to the DC second‐order Møller–Plesset perturbation (MP2) calculation, a scheme for estimating the subsystem MP2 correlation energy contribution from each atom in the buffer region is proposed. The estimation is based on the atomic orbital Laplace MP2 formalism. Based on this, an automatic buffer determination scheme for the DC‐MP2 calculation is constructed and its performance for several types of systems is assessed.
Collapse
Affiliation(s)
- Toshikazu Fujimori
- Graduate School of Chemical Sciences and EngineeringHokkaido UniversitySapporoJapan
| | - Masato Kobayashi
- Department of Chemistry, Faculty of ScienceHokkaido UniversitySapporoJapan
- WPI‐ICReDDHokkaido UniversitySapporoJapan
- ESICB, Kyoto UniversityKyotoJapan
| | - Tetsuya Taketsugu
- Department of Chemistry, Faculty of ScienceHokkaido UniversitySapporoJapan
- WPI‐ICReDDHokkaido UniversitySapporoJapan
- ESICB, Kyoto UniversityKyotoJapan
| |
Collapse
|
45
|
Kowalski K, Bair R, Bauman NP, Boschen JS, Bylaska EJ, Daily J, de Jong WA, Dunning T, Govind N, Harrison RJ, Keçeli M, Keipert K, Krishnamoorthy S, Kumar S, Mutlu E, Palmer B, Panyala A, Peng B, Richard RM, Straatsma TP, Sushko P, Valeev EF, Valiev M, van Dam HJJ, Waldrop JM, Williams-Young DB, Yang C, Zalewski M, Windus TL. From NWChem to NWChemEx: Evolving with the Computational Chemistry Landscape. Chem Rev 2021; 121:4962-4998. [PMID: 33788546 DOI: 10.1021/acs.chemrev.0c00998] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
Since the advent of the first computers, chemists have been at the forefront of using computers to understand and solve complex chemical problems. As the hardware and software have evolved, so have the theoretical and computational chemistry methods and algorithms. Parallel computers clearly changed the common computing paradigm in the late 1970s and 80s, and the field has again seen a paradigm shift with the advent of graphical processing units. This review explores the challenges and some of the solutions in transforming software from the terascale to the petascale and now to the upcoming exascale computers. While discussing the field in general, NWChem and its redesign, NWChemEx, will be highlighted as one of the early codesign projects to take advantage of massively parallel computers and emerging software standards to enable large scientific challenges to be tackled.
Collapse
Affiliation(s)
- Karol Kowalski
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Raymond Bair
- Argonne National Laboratory, Lemont, Illinois 60439, United States
| | - Nicholas P Bauman
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | - Eric J Bylaska
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Jeff Daily
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Wibe A de Jong
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Thom Dunning
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Niranjan Govind
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Robert J Harrison
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794, United States
| | - Murat Keçeli
- Argonne National Laboratory, Lemont, Illinois 60439, United States
| | | | | | - Suraj Kumar
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Erdal Mutlu
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bruce Palmer
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Ajay Panyala
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Bo Peng
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | - T P Straatsma
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831-6373, United States
| | - Peter Sushko
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Edward F Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Marat Valiev
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | | | | | | | - Chao Yang
- Lawrence Berkeley National Laboratory, Berkeley, California 94720, United States
| | - Marcin Zalewski
- Pacific Northwest National Laboratory, Richland, Washington 99352, United States
| | - Theresa L Windus
- Department of Chemistry, Iowa State University and Ames Laboratory, Ames, Iowa 50011, United States
| |
Collapse
|
46
|
Ma Q, Werner HJ. Scalable Electron Correlation Methods. 8. Explicitly Correlated Open-Shell Coupled-Cluster with Pair Natural Orbitals PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. J Chem Theory Comput 2021; 17:902-926. [PMID: 33405921 DOI: 10.1021/acs.jctc.0c01129] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We present explicitly correlated open-shell pair natural orbital local coupled-cluster methods, PNO-RCCSD(T)-F12 and PNO-UCCSD(T)-F12. The methods are extensions of our previously reported PNO-R/UCCSD methods (J. Chem. Theory Comput., 2020, 16, 3135-3151, https://pubs.acs.org/doi/10.1021/acs.jctc.0c00192) with additions of explicit correlation and perturbative triples corrections. The explicit correlation treatment follows the spin-orbital CCSD-F12b theory using Ansatz 3*A, which is found to yield comparable or better basis set convergence than the more rigorous Ansatz 3C in computed ionization potentials and reaction energies using double- to quaduple-ζ basis sets. The perturbative triples correction is adapted from the spin-orbital (T) theory to use triples natural orbitals (TNOs). To address the coupling due to off-diagonal Fock matrix elements, the local triples amplitudes are iteratively solved using small domains of TNOs, and a semicanonical (T0) domain correction with larger domains is applied to reduce the domain errors. The performance of the methods is demonstrated through benchmark calculations on ionization potentials, radical stabilization energies, reaction energies of fragmentations and rearrangements in radical cations, and spin-state energy differences of iron complexes. For a few test sets where canonical calculations are feasible, PNO-RCCSD(T)-F12 results agree with the canonical ones to within 0.4 kcal mol-1, and this maximum error is reduced to below 0.2 kcal mol-1 when large local domains are used. For larger systems, results using different thresholds for the local approximations are compared to demonstrate that 1 kcal mol-1 level of accuracy can be achieved using our default settings. For a couple of difficult cases, it is demonstrated that the errors from individual approximations are only a fraction of 1 kcal mol-1, and the overall accuracy of the method does not rely on error compensations. In contrast to canonical calculations, the use of spin-orbitals does not lead to a significant increase of computational time and memory usage in the most expensive steps of PNO-R/UCCSD(T)-F12 calculations. The only exception is the iterative solution of the (T) amplitudes, which can be avoided without significant errors by using a perturbative treatment of the off-diagonal coupling, known as (T1) approximation. For most systems, even the semicanonical approximation (T0) leads only to small errors in relative energies. Our program is well parallelized and capable of computing accurate correlation energies for molecules with 100-200 atoms using augmented triple-ζ basis sets in less than a day of elapsed time on a small computer cluster.
Collapse
Affiliation(s)
- Qianli Ma
- 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
| |
Collapse
|
47
|
Low-Scaling Tensor Hypercontraction in the Cholesky Molecular Orbital Basis Applied to Second-Order Møller-Plesset Perturbation Theory. J Chem Theory Comput 2020; 17:211-221. [PMID: 33375790 DOI: 10.1021/acs.jctc.0c00934] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
We employ various reduced scaling techniques to accelerate the recently developed least-squares tensor hypercontraction (LS-THC) approximation [Parrish, R. M., Hohenstein, E. G., Martínez, T. J., Sherrill, C. D. J. Chem. Phys. 137, 224106 (2012)] for electron repulsion integrals (ERIs) and apply it to second-order Møller-Plesset perturbation theory (MP2). The grid-projected ERI tensors are efficiently constructed using a localized Cholesky molecular orbital basis from density-fitted integrals with an attenuated Coulomb metric. Additionally, rigorous integral screening and the natural blocking matrix format are applied to reduce the complexity of this step. By recasting the equations to form the quantized representation of the 1/r operator Z into the form of a system of linear equations, the bottleneck of inverting the grid metric via pseudoinversion is removed. This leads to a reduced scaling THC algorithm and application to MP2 yields the (sub-)quadratically scaling THC-ω-RI-CDD-SOS-MP2 method. The efficiency of this method is assessed for various systems including DNA fragments with over 8000 basis functions and the subquadratic scaling is illustrated.
Collapse
|
48
|
Whitten JL. Estimates of electron correlation based on density expansions. J Chem Phys 2020; 153:244103. [PMID: 33380085 DOI: 10.1063/5.0031279] [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/14/2022] Open
Abstract
Two methods for estimating the correlation energy of molecules and other electronic systems are discussed based on the assumption that the correlation energy can be partitioned between atomic regions. In the first method, the electron density is expanded in terms of atomic contributions using rigorous electron repulsion bounds, and in the second method, correlation contributions are associated with basis function pairs. These methods do not consider the detailed nature of localized excitations but instead define a correlation energy per electron factor that is unique to a specific atom. The correlation factors are basis function dependent and are determined by configuration interaction (CI) calculations on diatomic and hydride molecules. The correlation energy estimates are compared with the results of high-level CI calculations for a test set of 27 molecules representing a wide range of bonding environments (average error of 2.6%). An extension based on truncated CI calculations in which d-type and hydrogen p-type functions are eliminated from the virtual space combined with estimates of dynamical correlation contributions using atomic correlation factors is discussed and applied to the dissociation of several molecules.
Collapse
Affiliation(s)
- Jerry L Whitten
- Department of Chemistry, North Carolina State University, Raleigh, North Carolina 27695, USA
| |
Collapse
|
49
|
Head-Marsden K, Flick J, Ciccarino CJ, Narang P. Quantum Information and Algorithms for Correlated Quantum Matter. Chem Rev 2020; 121:3061-3120. [PMID: 33326218 DOI: 10.1021/acs.chemrev.0c00620] [Citation(s) in RCA: 37] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
Discoveries in quantum materials, which are characterized by the strongly quantum-mechanical nature of electrons and atoms, have revealed exotic properties that arise from correlations. It is the promise of quantum materials for quantum information science superimposed with the potential of new computational quantum algorithms to discover new quantum materials that inspires this Review. We anticipate that quantum materials to be discovered and developed in the next years will transform the areas of quantum information processing including communication, storage, and computing. Simultaneously, efforts toward developing new quantum algorithmic approaches for quantum simulation and advanced calculation methods for many-body quantum systems enable major advances toward functional quantum materials and their deployment. The advent of quantum computing brings new possibilities for eliminating the exponential complexity that has stymied simulation of correlated quantum systems on high-performance classical computers. Here, we review new algorithms and computational approaches to predict and understand the behavior of correlated quantum matter. The strongly interdisciplinary nature of the topics covered necessitates a common language to integrate ideas from these fields. We aim to provide this common language while weaving together fields across electronic structure theory, quantum electrodynamics, algorithm design, and open quantum systems. Our Review is timely in presenting the state-of-the-art in the field toward algorithms with nonexponential complexity for correlated quantum matter with applications in grand-challenge problems. Looking to the future, at the intersection of quantum information science and algorithms for correlated quantum matter, we envision seminal advances in predicting many-body quantum states and describing excitonic quantum matter and large-scale entangled states, a better understanding of high-temperature superconductivity, and quantifying open quantum system dynamics.
Collapse
Affiliation(s)
- Kade Head-Marsden
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Johannes Flick
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, United States
| | - Christopher J Ciccarino
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States.,Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts 02138, United States
| | - Prineha Narang
- John A. Paulson School of Engineering and Applied Sciences, Harvard University, Cambridge, Massachusetts 02138, United States
| |
Collapse
|
50
|
Calvin JA, Peng C, Rishi V, Kumar A, Valeev EF. Many-Body Quantum Chemistry on Massively Parallel Computers. Chem Rev 2020; 121:1203-1231. [DOI: 10.1021/acs.chemrev.0c00006] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023]
Affiliation(s)
- Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, United States
| |
Collapse
|