1
|
de Moura CEV, Sokolov AY. Efficient Spin-Adapted Implementation of Multireference Algebraic Diagrammatic Construction Theory. I. Core-Ionized States and X-ray Photoelectron Spectra. J Phys Chem A 2024; 128:5816-5831. [PMID: 38962857 DOI: 10.1021/acs.jpca.4c03161] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/05/2024]
Abstract
We present an efficient implementation of multireference algebraic diagrammatic construction theory (MR-ADC) for simulating core-ionized states and X-ray photoelectron spectra (XPS). Taking advantage of spin adaptation, automatic code generation, and density fitting, our implementation can perform calculations for molecules with more than 1500 molecular orbitals, incorporating static and dynamic correlation in the ground and excited electronic states. We demonstrate the capabilities of MR-ADC methods by simulating the XPS spectra of substituted ferrocene complexes and azobenzene isomers. For the ground electronic states of these molecules, the XPS spectra computed using the extended second-order MR-ADC method (MR-ADC(2)-X) are in a very good agreement with available experimental results. We further show that MR-ADC can be used as a tool for interpreting or predicting the results of time-resolved XPS measurements by simulating the core ionization spectra of azobenzene along its photoisomerization, including the XPS signatures of excited states and the minimum energy conical intersection. This work is the first in a series of publications reporting the efficient implementations of MR-ADC methods.
Collapse
Affiliation(s)
- Carlos E V de Moura
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, United States
| |
Collapse
|
2
|
Fransson T, Pettersson LGM. TDDFT and the x-ray absorption spectrum of liquid water: Finding the "best" functional. J Chem Phys 2024; 160:234105. [PMID: 38884399 DOI: 10.1063/5.0209719] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2024] [Accepted: 06/03/2024] [Indexed: 06/18/2024] Open
Abstract
We investigate the performance of time-dependent density functional theory (TDDFT) for reproducing high-level reference x-ray absorption spectra of liquid water and water clusters. For this, we apply the integrated absolute difference (IAD) metric, previously used for x-ray emission spectra of liquid water [T. Fransson and L. G. M. Pettersson, J. Chem. Theory Comput. 19, 7333-7342 (2023)], in order to investigate which exchange-correlation (xc) functionals yield TDDFT spectra in best agreement to reference, as well as to investigate the suitability of IAD for x-ray absorption spectroscopy spectrum calculations. Considering highly asymmetric and symmetric six-molecule clusters, it is seen that long-range corrected xc-functionals are required to yield good agreement with the reference coupled cluster (CC) and algebraic-diagrammatic construction spectra, with 100% asymptotic Hartree-Fock exchange resulting in the lowest IADs. The xc-functionals with best agreement to reference have been adopted for larger water clusters, yielding results in line with recently published CC theory, but which still show some discrepancies in the relative intensity of the features compared to experiment.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| | - Lars G M Pettersson
- Department of Physics, AlbaNova University Center, Stockholm University, 10961 Stockholm, Sweden
| |
Collapse
|
3
|
Hicks CB, Martinez TJ. Massively scalable workflows for quantum chemistry: BigChem and ChemCloud. J Chem Phys 2024; 160:142501. [PMID: 38591672 DOI: 10.1063/5.0190834] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/08/2023] [Accepted: 03/14/2024] [Indexed: 04/10/2024] Open
Abstract
Electronic structure theory, i.e., quantum chemistry, is the fundamental building block for many problems in computational chemistry. We present a new distributed computing framework (BigChem), which allows for an efficient solution of many quantum chemistry problems in parallel. BigChem is designed to be easily composable and leverages industry-standard middleware (e.g., Celery, RabbitMQ, and Redis) for distributed approaches to large scale problems. BigChem can harness any collection of worker nodes, including ones on cloud providers (such as AWS or Azure), local clusters, or supercomputer centers (and any mixture of these). BigChem builds upon MolSSI packages, such as QCEngine to standardize the operation of numerous computational chemistry programs, demonstrated here with Psi4, xtb, geomeTRIC, and TeraChem. BigChem delivers full utilization of compute resources at scale, offers a programable canvas for designing sophisticated quantum chemistry workflows, and is fault tolerant to node failures and network disruptions. We demonstrate linear scalability of BigChem running computational chemistry workloads on up to 125 GPUs. Finally, we present ChemCloud, a web API to BigChem and successor to TeraChem Cloud. ChemCloud delivers scalable and secure access to BigChem over the Internet.
Collapse
Affiliation(s)
- Colton B Hicks
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| | - Todd J Martinez
- Department of Chemistry and The PULSE Institute, Stanford University, Stanford, California 94305, USA and SLAC National Accelerator Laboratory, Menlo Park, California 94025, USA
| |
Collapse
|
4
|
Brand M, Dreuw A, Norman P, Li X. Efficient and Parallel Implementation of Real and Complex Response Functions Employing the Second-Order Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator. J Chem Theory Comput 2024; 20:103-113. [PMID: 38117937 PMCID: PMC10782457 DOI: 10.1021/acs.jctc.3c01065] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Revised: 11/29/2023] [Accepted: 12/01/2023] [Indexed: 12/22/2023]
Abstract
We present the implementation of an efficient matrix-folded formalism for the evaluation of complex response functions and the calculation of transition properties at the level of the second-order algebraic-diagrammatic construction (ADC(2)) scheme. The underlying algorithms, in combination with the adopted hybrid MPI/OpenMP parallelization strategy, enabled calculations of the UV/vis spectra of a guanine oligomer series ranging up to 1032 contracted basis functions, thereby utilizing vast computational resources from up to 32,768 CPU cores. Further analysis of the convergence behavior of the involved iterative subspace algorithms revealed the superiority of a frequency-separated treatment of response equations even for a large spectral window, including 101 frequencies. We demonstrate the applicability to general quantum mechanical operators by the first reported electronic circular dichroism spectrum calculated with a complex polarization propagator approach at the ADC(2) level of theory.
Collapse
Affiliation(s)
- Manuel Brand
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, Stockholm SE-100 44, Sweden
| | - Andreas Dreuw
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer Feld 205, Heidelberg 69120, Germany
| | - Patrick Norman
- Division
of Theoretical Chemistry and Biology, School of Engineering Sciences
in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, Stockholm SE-100 44, Sweden
| | - Xin Li
- PDC
Center for High Performance Computing, KTH
Royal Institute of Technology, Stockholm SE-100 44, Sweden
| |
Collapse
|
5
|
Lehtola S. A call to arms: Making the case for more reusable libraries. J Chem Phys 2023; 159:180901. [PMID: 37947507 DOI: 10.1063/5.0175165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/05/2023] [Accepted: 10/23/2023] [Indexed: 11/12/2023] Open
Abstract
The traditional foundation of science lies on the cornerstones of theory and experiment. Theory is used to explain experiment, which in turn guides the development of theory. Since the advent of computers and the development of computational algorithms, computation has risen as the third cornerstone of science, joining theory and experiment on an equal footing. Computation has become an essential part of modern science, amending experiment by enabling accurate comparison of complicated theories to sophisticated experiments, as well as guiding by triage both the design and targets of experiments and the development of novel theories and computational methods. Like experiment, computation relies on continued investment in infrastructure: it requires both hardware (the physical computer on which the calculation is run) as well as software (the source code of the programs that performs the wanted simulations). In this Perspective, I discuss present-day challenges on the software side in computational chemistry, which arise from the fast-paced development of algorithms, programming models, as well as hardware. I argue that many of these challenges could be solved with reusable open source libraries, which are a public good, enhance the reproducibility of science, and accelerate the development and availability of state-of-the-art methods and improved software.
Collapse
Affiliation(s)
- Susi Lehtola
- Department of Chemistry, University of Helsinki, P.O. Box 55, FI-00014 Helsinki, Finland
| |
Collapse
|
6
|
Fransson T, Pettersson LGM. Calibrating TDDFT Calculations of the X-ray Emission Spectrum of Liquid Water: The Effects of Hartree-Fock Exchange. J Chem Theory Comput 2023; 19:7333-7342. [PMID: 37787584 PMCID: PMC10601479 DOI: 10.1021/acs.jctc.3c00728] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2023] [Indexed: 10/04/2023]
Abstract
The structure and dynamics of liquid water continue to be debated, with insight provided by, among others, X-ray emission spectroscopy (XES), which shows a split in the high-energy 1b1 feature. This split is yet to be reproduced by theory, and it remains unclear if these difficulties are related to inaccuracies in dynamics simulations, spectrum calculations, or both. We investigate the performance of different methods for calculating XES of liquid water, focusing on the ability of time-dependent density functional theory (TDDFT) to reproduce reference spectra obtained by high-level coupled cluster and algebraic-diagrammatic construction scheme calculations. A metric for evaluating the agreement between theoretical spectra termed the integrated absolute difference (IAD), which considers the integral of shifted difference spectra, is introduced and used to investigate the performance of different exchange-correlation functionals. We find that computed spectra of symmetric and asymmetric model water structures are strongly and differently influenced by the amount of Hartree-Fock exchange, with best agreement to reference spectra for ∼40-50%. Lower percentages tend to yield high density of contributing states, resulting in too broad features. The method introduced here is useful also for other spectrum calculations, in particular where the performance for ensembles of structures are evaluated.
Collapse
Affiliation(s)
- Thomas Fransson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| | - Lars G. M. Pettersson
- Department of Physics, AlbaNova
University Center, Stockholm University, 109 61 Stockholm, Sweden
| |
Collapse
|
7
|
Schneider F, Papapostolou A, Leitner J, Rehn DR, Dreuw A. Cotton-Mouton Effect Using Algebraic Diagrammatic Construction Schemes in the Intermediate State Representation Formalism. J Phys Chem A 2023; 127:8723-8733. [PMID: 37816160 DOI: 10.1021/acs.jpca.3c04963] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/12/2023]
Abstract
The Cotton-Mouton effect is theoretically investigated for a selected set of molecules by using a novel computational methodology based on algebraic diagrammatic construction (ADC) schemes in the intermediate state representation (ISR) formulation. Therefore, the electronic contributions to the frequency-dependent polarizabilities and, for the first time, to the magnetizabilities as well as mixed electric and magnetic hypermagnetizabilities have been computed in the ADC/ISR framework. In addition to calculation of the Cotton-Mouton constant and the birefringence, the gauge origin dependence of the computed tensors and the applied methodology are thoroughly investigated. The new ADC/ISR methodology, employing the recently presented responsefun package, is applied to a test set of Ne and small molecules (H2, HF, O2, CO2, and benzene) and compared to data from the experiment as well as other ab initio methods. The presented theoretical ab initio ADC/ISR approach is a substantial extension of the available computational methods for the investigation of complex nonlinear properties, however, with a gauge origin dependence inherent to the method that decreases with increasing perturbation order.
Collapse
Affiliation(s)
- Friederike Schneider
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| | - Antonia Papapostolou
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| | - Jonas Leitner
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| | - Dirk R Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| |
Collapse
|
8
|
Papapostolou A, Scheurer M, Dreuw A, Rehn DR. responsefun: Fun with Response Functions in the Algebraic Diagrammatic Construction Framework. J Chem Theory Comput 2023; 19:6375-6391. [PMID: 37676497 DOI: 10.1021/acs.jctc.3c00456] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 09/08/2023]
Abstract
We present the open-source responsefun package, which implements a universally applicable procedure for computing molecular response properties within the algebraic diagrammatic construction (ADC) framework, exploiting the intermediate state representation (ISR) approach. With symbolic mathematics, the user can simply enter textbook sum-over-states (SOS) expressions from time-dependent perturbation theory, which are then automatically translated into the corresponding symbolic ADC/ISR formulations. Using the data structures provided by the hybrid Python/C++ module adcc for calculating excited states with ADC, the specified response property is directly evaluated, and the result is returned to the user. Employing the novel responsefun package, we present the first ADC/ISR calculations of second-order hyperpolarizability tensors and three-photon-absorption matrix elements.
Collapse
Affiliation(s)
- Antonia Papapostolou
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Dirk R Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
9
|
Dreuw A, Papapostolou A, Dempwolff AL. Algebraic Diagrammatic Construction Schemes Employing the Intermediate State Formalism: Theory, Capabilities, and Interpretation. J Phys Chem A 2023; 127:6635-6646. [PMID: 37498297 DOI: 10.1021/acs.jpca.3c02761] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/28/2023]
Abstract
Algebraic diagrammatic construction (ADC) schemes represent a family of ab initio methods for the calculation of excited electronic states and electron-detached and -attached states. All ADC methods have been demonstrated to possess great potential for molecular applications, e.g., for the calculation of absorption or photoelectron spectra or electron attachment processes. ADC originates from Green's function or propagator theory; however, most recent ADC developments heavily rely on the intermediate state representation or effective Liouvillian formalisms, which comprise new ADC methods and computational schemes for high-order properties. The different approaches for the calculation of excitation energies, ionization potentials, and electron affinities are intimately related, and they provide a coherent description of these quantities at equivalent levels of theory and with comparable errors. Most quantum chemical program packages contain ADC methods; however, the most complete ADC suite of methods can be found in the recent release of Q-Chem.
Collapse
Affiliation(s)
- Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Antonia Papapostolou
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| | - Adrian L Dempwolff
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, 69120 Heidelberg, Germany
| |
Collapse
|
10
|
Kulahlioglu AH, Dreuw A. The Multistate Quantum Monte Carlo Algebraic Diagrammatic Construction Method. J Phys Chem A 2023; 127:2161-2175. [PMID: 36847774 DOI: 10.1021/acs.jpca.2c08391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/01/2023]
Abstract
A multistate formulation of the recently developed quantum Monte Carlo (QMC) algebraic diagrammatic construction (ADC) method, QMCADC, is presented. QMCADC solves the Hermitian eigenvalue problem of the second-order ADC scheme for the polarization propagator stochastically by combining ADC schemes with projector quantum Monte Carlo (PQMC). It allows for massively parallel distributed computing and exploits the sparsity of the effective ADC matrix, thereby relaxing memory and processing requirements of ADC methods significantly. Here, the theory and implementation of the multistate variant of QMCADC are described, and our first proof-of-principle calculations for various molecular systems are shown. Indeed, multistate QMCADC enables sampling of an arbitrary number of low-lying excited states and can reproduce their vertical excitation energies with a marginal controllable error. The performance of multistate QMCADC is examined in terms of state-wise and overall accuracy as well as with respect to the balance in the treatments of excited states relatively to each other. The results are very promising as they show bias and imbalances among excited states to diminish as the number of sampling points increases. Furthermore, the impact of the quality of trial wave functions on the vertical excitation energies is investigated. A black-box approach for the generation of high quality trial wave functions internally is given.
Collapse
Affiliation(s)
- Adem Halil Kulahlioglu
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
11
|
Scheurer M, Papapostolou A, Fransson T, Norman P, Dreuw A, Rehn DR. Solving response expressions in the ADC/ISR framework. J Chem Phys 2023; 158:084105. [PMID: 36859074 DOI: 10.1063/5.0139721] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023] Open
Abstract
We present an implementation for the calculation of molecular response properties using the algebraic-diagrammatic construction (ADC)/intermediate state representation approach. For the second-order ADC model [ADC(2)], a memory-efficient ansatz avoiding the storage of double excitation amplitudes is investigated. We compare the performance of different numerical algorithms for the solution of the underlying response equations for ADC(2) and show that our approach also strongly improves the convergence behavior for the investigated algorithms compared with the standard implementation. All routines are implemented in an open-source Python library.
Collapse
Affiliation(s)
- Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Antonia Papapostolou
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Thomas Fransson
- Division of Theoretical Chemistry and Biology, School of Biotechnology, Royal Institute of Technology, SE-10044 Stockholm, Sweden
| | - Patrick Norman
- Division of Theoretical Chemistry and Biology, Royal Institute of Technology, SE-10044 Stockholm, Sweden
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Dirk R Rehn
- Interdisciplinary Center for Scientific Computing, University of Heidelberg, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
12
|
Mester D, Kállay M. Double-Hybrid Density Functional Theory for Core Excitations: Theory and Benchmark Calculations. J Chem Theory Comput 2023; 19:1310-1321. [PMID: 36721871 PMCID: PMC9979613 DOI: 10.1021/acs.jctc.2c01222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023]
Abstract
The double-hybrid (DH) time-dependent density functional theory is extended to core excitations. Two different DH formalisms are presented utilizing the core-valence separation (CVS) approximation. First, a CVS-DH variant is introduced relying on the genuine perturbative second-order correction, while an iterative analogue is also presented using our second-order algebraic-diagrammatic construction [ADC(2)]-based DH ansatz. The performance of the new approaches is tested for the most popular DH functionals using the recently proposed XABOOM [J. Chem. Theory Comput.2021, 17, 1618] benchmark set. In order to make a careful comparison, the accuracy and precision of the methods are also inspected. Our results show that the genuine approaches are highly competitive with the more advanced CVS-ADC(2)-based methods if only excitation energies are required. In contrast, as expected, significant differences are observed in oscillator strengths; however, the precision is acceptable for the genuine functionals as well. Concerning the performance of the CVS-DH approaches, the PBE0-2/CVS-ADC(2) functional is superior, while its spin-opposite-scaled variant is also recommended as a cost-effective alternative. For these approaches, significant improvements are realized in the error measures compared with the popular CVS-ADC(2) method.
Collapse
Affiliation(s)
- Dávid Mester
- Department
of Physical Chemistry and Materials Science, Faculty of Chemical Technology
and Biotechnology, Budapest University of
Technology and Economics, Müegyetem rkp. 3, H-1111Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, 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-1111Budapest, Hungary,ELKH-BME
Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,MTA-BME
Lendület Quantum Chemistry Research Group, Müegyetem rkp. 3, H-1111Budapest, Hungary,
| |
Collapse
|
13
|
Lehtola S, Karttunen AJ. Free and open source software for computational chemistry education. WIRES COMPUTATIONAL MOLECULAR SCIENCE 2022. [DOI: 10.1002/wcms.1610] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Susi Lehtola
- Molecular Sciences Software Institute Blacksburg Virginia USA
| | - Antti J. Karttunen
- Department of Chemistry and Materials Science Aalto University Espoo Finland
| |
Collapse
|
14
|
Brumboiu IE, Fransson T. Core-hole delocalization for modeling X-ray spectroscopies: A cautionary tale. J Chem Phys 2022; 156:214109. [DOI: 10.1063/5.0088195] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
The influence of core-hole delocalization for X-ray photoelectron, X-ray absorption, and X-ray emission spectrum calculations is investigated in detail, using approaches including response theory, transition-potential methods, and ground state schemes. The question of a localized/delocalized vacancy is relevant for systems with symmetrically equivalent atoms, as well as near-degeneracies which can distribute the core-orbitals over several atoms. We show that issues relating to core-hole delocalization are present for calculations considering explicit core-hole states, e.g. when using a core-excited or core-ionized reference state, or for fractional occupation numbers. Including electron correlation eventually alleviates the issues, but even using CCSD(T) there is a noticable discrepancy between core-ionization energies obtained with a localized and delocalized core-hole (0.5 eV for the carbon K-edge). Within density functional theory, the discrepancy correlates to the exchange interaction involving the core orbitals of the same spin symmetry as the delocalized core-hole. The use of a localized core-hole allows for a reasonably good inclusion of relaxation at lower level of theory, whereas the proper symmetry solution involving a delocalized core-hole requires higher levels of theory to account for the correlation effects involved in orbital relaxation. For linear response methods, we further show that if X-ray absorption spectra are modelled by considering symmetry-unique sets of atoms, care has to be taken such that there are no delocalizations of the core orbitals, which would otherwise introduce shifts in absolute energies and relative features.
Collapse
Affiliation(s)
- Iulia Emilia Brumboiu
- Faculty of Physics, Astronomy and Informatics, Nicolaus Copernicus University in Torun Institute of Physics, Poland
| | | |
Collapse
|
15
|
Kulahlioglu AH, Rehn D, Dreuw A. Quantum Monte Carlo formulation of the second order algebraic diagrammatic construction: Toward a massively parallel correlated excited state method. J Chem Phys 2022; 156:044105. [DOI: 10.1063/5.0071091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Adem Halil Kulahlioglu
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Dirk Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
16
|
González-Espinoza CE, Rumble CA, Borgis D, Wesolowski TA. Quantifying Fluctuations of Average Solvent Environments for Embedding Calculations. J Chem Theory Comput 2022; 18:1072-1088. [PMID: 35044168 DOI: 10.1021/acs.jctc.1c01108] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The viability and effectiveness of replacing an ensemble of embedded solute calculations by a single calculation using an average description of the solvent environment are evaluated. This work explores the fluctuations of the average description of the system obtained in two ways: from calculations on an ensemble of geometries and from an average environment constructed from the same ensemble. To this end, classical molecular dynamics simulations of a rigid acetone solute in SPCE water are performed in order to generate an ensemble of solvent environments. From this ensemble of solvent configurations, a number of different approaches for constructing an average solvent environment are employed. We perform a thorough numerical analysis of the fluctuations of the electrostatic potential experienced by the solute, as well as the resulting fluctuations of the solute's electronic density, measured through its dipole moment and fitted atomic point charges. At the same time, we inspect the accuracy of the methods used to construct average environments. Finally, the proposed method for generating the embedding potential from an average environment density is applied to estimate the solvatochromic shift of the first excitation of acetone. In order to account for quantum confinement effects, which may be important in certain cases, the fluctuations in the shift due to the interaction with the solvent are evaluated using frozen-density-embedding theory. Our results demonstrate that, for normally distributed environments, the constructed average environment is a reasonably good representation of a fluctuating molecular solvent environment. We then provide guidance for future comparisons between these theoretical treatments of solute/solvent systems to experimental measurements.
Collapse
Affiliation(s)
| | - Christopher A Rumble
- Université de Genève, Départment de Chimie Physique 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| | - Daniel Borgis
- Ecole Normale Supérieure, Départment de Chimie, UMR 8640, PSL University, Sorbonne University, CNRS, 75005 Paris, France
| | - Tomasz A Wesolowski
- Université de Genève, Départment de Chimie Physique 30, Quai Ernest-Ansermet, CH-1211 Genève 4, Switzerland
| |
Collapse
|
17
|
Huang M, Li C, Evangelista FA. Theoretical Calculation of Core-Excited States along Dissociative Pathways beyond Second-Order Perturbation Theory. J Chem Theory Comput 2021; 18:219-233. [PMID: 34964628 DOI: 10.1021/acs.jctc.1c00884] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022]
Abstract
We extend the multireference driven similarity renormalization (MR-DSRG) method to compute core-excited states by combining it with a GASSCF treatment of orbital relaxation and static electron correlation effects. We consider MR-DSRG treatments of dynamical correlation truncated at the level of perturbation theory (DSRG-MRPT2/3) and iterative linearized approximations with one- and two-body operators [MR-LDSRG(2)] in combination with a spin-free exact-two-component (X2C) one-electron treatment of scalar relativistic effects. This approach is calibrated and tested on a series of 16 core-excited states of five closed- and open-shell diatomic molecules containing first-row elements (C, N, and O). All GASSCF-MR-DSRG theories show excellent agreement with experimental adiabatic transitions energies, with mean absolute errors ranging between 0.17 and 0.35 eV, even for the challenging partially doubly excited states of the N2+ molecule. The vibrational structure of all these transitions, obtained from using a full potential energy scan, shows a mean absolute error as low as 25 meV for DSRG-MRPT2 and 12/13 meV for DSRG-MRPT3 and MR-LDSRG(2). We generally find that a treatment of dynamical correlation that goes beyond the second-order level in perturbation theory improves the accuracy of the potential energy surface, especially in the bond-dissociation region.
Collapse
Affiliation(s)
- Meng Huang
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| | - Chenyang Li
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States.,Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Francesco A Evangelista
- Department of Chemistry and Cherry Emerson Center for Scientific Computation, Emory University, Atlanta, Georgia 30322, United States
| |
Collapse
|
18
|
Smith DGA, Lolinco AT, Glick ZL, Lee J, Alenaizan A, Barnes TA, Borca CH, Di Remigio R, Dotson DL, Ehlert S, Heide AG, Herbst MF, Hermann J, Hicks CB, Horton JT, Hurtado AG, Kraus P, Kruse H, Lee SJR, Misiewicz JP, Naden LN, Ramezanghorbani F, Scheurer M, Schriber JB, Simmonett AC, Steinmetzer J, Wagner JR, Ward L, Welborn M, Altarawy D, Anwar J, Chodera JD, Dreuw A, Kulik HJ, Liu F, Martínez TJ, Matthews DA, Schaefer HF, Šponer J, Turney JM, Wang LP, De Silva N, King RA, Stanton JF, Gordon MS, Windus TL, Sherrill CD, Burns LA. Quantum Chemistry Common Driver and Databases (QCDB) and Quantum Chemistry Engine (QCEngine): Automation and interoperability among computational chemistry programs. J Chem Phys 2021; 155:204801. [PMID: 34852489 DOI: 10.1063/5.0059356] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Community efforts in the computational molecular sciences (CMS) are evolving toward modular, open, and interoperable interfaces that work with existing community codes to provide more functionality and composability than could be achieved with a single program. The Quantum Chemistry Common Driver and Databases (QCDB) project provides such capability through an application programming interface (API) that facilitates interoperability across multiple quantum chemistry software packages. In tandem with the Molecular Sciences Software Institute and their Quantum Chemistry Archive ecosystem, the unique functionalities of several CMS programs are integrated, including CFOUR, GAMESS, NWChem, OpenMM, Psi4, Qcore, TeraChem, and Turbomole, to provide common computational functions, i.e., energy, gradient, and Hessian computations as well as molecular properties such as atomic charges and vibrational frequency analysis. Both standard users and power users benefit from adopting these APIs as they lower the language barrier of input styles and enable a standard layout of variables and data. These designs allow end-to-end interoperable programming of complex computations and provide best practices options by default.
Collapse
Affiliation(s)
- Daniel G A Smith
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | | | - Zachary L Glick
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Jiyoung Lee
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - Asem Alenaizan
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Taylor A Barnes
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Carlos H Borca
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical and Computational Chemistry, UiT, The Arctic University of Norway, N-9037 Tromsø, Norway
| | - David L Dotson
- Open Force Field Initiative, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Sebastian Ehlert
- Mulliken Center for Theoretical Chemistry, Institut für Physikalische und Theoretische Chemie, Universität Bonn, Beringstraße 4, D-53115 Bonn, Germany
| | - Alexander G Heide
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Michael F Herbst
- Applied and Computational Mathematics, RWTH Aachen University, Schinkelstr. 2, 52062 Aachen, Germany
| | - Jan Hermann
- FU Berlin, Department of Mathematics and Computer Science, 14195 Berlin, Germany
| | - Colton B Hicks
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Joshua T Horton
- Department of Chemistry, Lancaster University, Lancaster LA1 4YW, United Kingdom
| | - Adrian G Hurtado
- Institute for Advanced Computational Science, Stony Brook University, Stony Brook, New York 11794-5250, USA
| | - Peter Kraus
- School of Molecular and Life Sciences, Curtin University, GPO Box U1987, Perth 6845, WA, Australia
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | | | - Jonathon P Misiewicz
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Levi N Naden
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | | | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Jeffrey B Schriber
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Andrew C Simmonett
- Laboratory of Computational Biology, National Institutes of Health-National Heart, Lung and Blood Institute, Bethesda, Maryland 20892, USA
| | - Johannes Steinmetzer
- Institute of Physical Chemistry, Friedrich Schiller University Jena, Jena, Germany
| | - Jeffrey R Wagner
- Open Force Field Initiative, University of Colorado Boulder, Boulder, Colorado 80309, USA
| | - Logan Ward
- Data Science and Learning Division, Argonne National Laboratory, Lemont, Illinois 60439, USA
| | - Matthew Welborn
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Doaa Altarawy
- Molecular Sciences Software Institute, Blacksburg, Virginia 24060, USA
| | - Jamshed Anwar
- Department of Chemistry, Lancaster University, Lancaster LA1 4YW, United Kingdom
| | - John D Chodera
- Computational and Systems Biology Program, Sloan Kettering Institute, Memorial Sloan Kettering Cancer Center, New York, New York 10065, USA
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, Im Neuenheimer Feld 205, 69120 Heidelberg, Germany
| | - Heather J Kulik
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Fang Liu
- Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, Massachusetts 02139, USA
| | - Todd J Martínez
- Department of Chemistry, Stanford University, Stanford, California 94305, USA
| | - Devin A Matthews
- The Institute for Computational Engineering and Sciences, The University of Texas at Austin, Austin, Texas 78712, USA
| | - Henry F Schaefer
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Jiří Šponer
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 65 Brno, Czech Republic
| | - Justin M Turney
- Center for Computational Quantum Chemistry, University of Georgia, Athens, Georgia 30602, USA
| | - Lee-Ping Wang
- Department of Chemistry, University of California Davis, Davis, California 95616, USA
| | - Nuwan De Silva
- Department of Chemistry, Iowa State University, Ames, Iowa 50011, USA
| | - Rollin A King
- Department of Chemistry, Bethel University, St. Paul, Minnesota 55112, USA
| | - John F Stanton
- Quantum Theory Project, The University of Florida, 2328 New Physics Building, Gainesville, Florida 32611-8435, USA
| | - Mark S Gordon
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - Theresa L Windus
- Department of Chemistry and Ames Laboratory, Iowa State University, Ames, Iowa 50011, USA
| | - C David Sherrill
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| | - Lori A Burns
- Center for Computational Molecular Science and Technology, School of Chemistry and Biochemistry, and School of Computational Science and Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332, USA
| |
Collapse
|
19
|
Lemmens L, De Vriendt X, Tolstykh D, Huysentruyt T, Bultinck P, Acke G. GQCP: The Ghent Quantum Chemistry Package. J Chem Phys 2021; 155:084802. [PMID: 34470369 DOI: 10.1063/5.0057515] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/14/2023] Open
Abstract
The Ghent Quantum Chemistry Package (GQCP) is an open-source electronic structure software package that aims to provide an intuitive and expressive software framework for electronic structure software development. Its high-level interfaces (accessible through C++ and Python) have been specifically designed to correspond to theoretical concepts, while retaining access to lower-level intermediates and allowing structural run-time modifications of quantum chemical solvers. GQCP focuses on providing quantum chemical method developers with the computational "building blocks" that allow them to flexibly develop proof of principle implementations for new methods and applications up to the level of two-component spinor bases.
Collapse
Affiliation(s)
- Laurent Lemmens
- Ghent Quantum Chemistry Group, Department of Chemistry, Ghent University, Krijgslaan 281 (S3), B-9000 Gent, Belgium
| | - Xeno De Vriendt
- Ghent Quantum Chemistry Group, Department of Chemistry, Ghent University, Krijgslaan 281 (S3), B-9000 Gent, Belgium
| | - Daria Tolstykh
- Ghent Quantum Chemistry Group, Department of Chemistry, Ghent University, Krijgslaan 281 (S3), B-9000 Gent, Belgium
| | - Tobias Huysentruyt
- Ghent Quantum Chemistry Group, Department of Chemistry, Ghent University, Krijgslaan 281 (S3), B-9000 Gent, Belgium
| | - Patrick Bultinck
- Ghent Quantum Chemistry Group, Department of Chemistry, Ghent University, Krijgslaan 281 (S3), B-9000 Gent, Belgium
| | - Guillaume Acke
- Ghent Quantum Chemistry Group, Department of Chemistry, Ghent University, Krijgslaan 281 (S3), B-9000 Gent, Belgium
| |
Collapse
|
20
|
Backhouse OJ, Santana-Bonilla A, Booth GH. Scalable and Predictive Spectra of Correlated Molecules with Moment Truncated Iterated Perturbation Theory. J Phys Chem Lett 2021; 12:7650-7658. [PMID: 34351782 DOI: 10.1021/acs.jpclett.1c02383] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/13/2023]
Abstract
A reliable and efficient computation of the entire single-particle spectrum of correlated molecules is an outstanding challenge in the field of quantum chemistry, with standard density functional theory approaches often giving an inadequate description of excitation energies and gaps. In this work, we expand upon a recently introduced approach that relies on a fully self-consistent many-body perturbation theory coupled to a nonperturbative truncation of the effective dynamics at each step. We show that this yields a low-scaling and accurate method across a diverse benchmark test set that is capable of treating moderate levels of strong correlation effects, and we detail an efficient implementation for applications involving up to ∼1000 orbitals on parallel resources. We then use this method to characterize the spectral properties of the antimalarial drug molecule artemisinin, resolving discrepancies in previous works concerning the active sites of the lowest-energy fundamental excitations of the system.
Collapse
Affiliation(s)
- Oliver J Backhouse
- Department of Physics, King's College London, Strand, London WC2R 2LS, U.K
| | | | - George H Booth
- Department of Physics, King's College London, Strand, London WC2R 2LS, U.K
| |
Collapse
|
21
|
Brumboiu IE, Rehn DR, Dreuw A, Rhee YM, Norman P. Analytical gradients for core-excited states in the algebraic diagrammatic construction (ADC) framework. J Chem Phys 2021; 155:044106. [PMID: 34340367 DOI: 10.1063/5.0058221] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
Expressions for analytical molecular gradients of core-excited states have been derived and implemented for the hierarchy of algebraic diagrammatic construction (ADC) methods up to extended second-order within the core-valence separation (CVS) approximation. We illustrate the use of CVS-ADC gradients by determining relaxed core-excited state potential energy surfaces and optimized geometries for water, formic acid, and benzene. For water, our results show that in the dissociative lowest core-excited state, a linear configuration is preferred. For formic acid, we find that the O K-edge lowest core-excited state is non-planar, a fact that is not captured by the equivalent core approximation where the core-excited atom with its hole is replaced by the "Z + 1" neighboring atom in the periodic table. For benzene, the core-excited state gradients are presented along the Jahn-Teller distorted geometry of the 1s → π* excited state. Our development may pave a new path to studying the dynamics of molecules in their core-excited states.
Collapse
Affiliation(s)
- Iulia Emilia Brumboiu
- Department of Chemistry, Pohang University of Science and Technology (POSTECH), 37673 Pohang, Republic of Korea
| | - Dirk R Rehn
- Interdisciplinary Center for Scientific Computing, Heidelberg University, 69120 Heidelberg, Germany
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, 69120 Heidelberg, Germany
| | - Young Min Rhee
- Department of Chemistry, Korea Advanced Institute of Science and Technology, 34141 Daejeon, Republic of Korea
| | - Patrick Norman
- Department of Theoretical Chemistry and Biology, KTH Royal Institute of Technology, 10691 Stockholm, Sweden
| |
Collapse
|
22
|
Scheurer M, Reinholdt P, Olsen JMH, Dreuw A, Kongsted J. Efficient Open-Source Implementations of Linear-Scaling Polarizable Embedding: Use Octrees to Save the Trees. J Chem Theory Comput 2021; 17:3445-3454. [PMID: 33949862 DOI: 10.1021/acs.jctc.1c00225] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
Abstract
We present open-source implementations of the linear-scaling fast multipole method (FMM) within the polarizable embedding (PE) model for efficient treatment of large polarizable environments in computational spectroscopy simulations. The implementations are tested for accuracy, efficiency, and usability on model systems as well as more realistic biomolecular systems. We explain how FMM parameters affect the calculation of molecular properties and show that PE calculations employing FMM can be carried out in a black-box manner. The efficiency of the linear-scaling approach is demonstrated by simulating the UV/vis spectrum of a chromophore in an environment of more than 1 million polarizable sites. Our implementations are interfaced to several open-source quantum chemistry programs, making computational spectroscopy simulations within the PE model and FMM available to a large variety of methods and a broad user base.
Collapse
Affiliation(s)
- Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Heidelberg University, D-69120 Heidelberg, Germany
| | - Peter Reinholdt
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense, Denmark
| | | | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Heidelberg University, D-69120 Heidelberg, Germany
| | - Jacob Kongsted
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense, Denmark
| |
Collapse
|
23
|
Bilalbegović G, Maksimović A, Valencic LA, Lehtola S. Sulfur Molecules in Space by X-rays: A Computational Study. ACS EARTH & SPACE CHEMISTRY 2021; 5:436-448. [PMID: 33842801 PMCID: PMC8028330 DOI: 10.1021/acsearthspacechem.0c00238] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/29/2020] [Revised: 02/12/2021] [Accepted: 02/15/2021] [Indexed: 06/12/2023]
Abstract
X-ray astronomy lacks high resolution spectra of interstellar dust analogues and molecules, severely hampering interstellar medium studies based on upcoming X-ray missions. Various theoretical approaches may be used to address this problem, but they must first be shown to reproduce reliable spectra compared to the experiment. In this work, we calculate the sulfur K edge X-ray absorption spectra of H2S, SO2, and OCS, whose spectra are already known from X-ray experiments and predict the X-ray spectrum of CS, which as far as we are aware has not been measured, thereby hampering its detection by X-ray telescopes. We chose these four molecules as the astrochemistry of sulfur is an unsolved problem and as the four molecules are already known to exist in space. We consider three types of methods for modeling the X-ray spectra: more accurate calculations with the algebraic-diagrammatic construction (ADC) and the CC2, CCSD, and CC3 coupled cluster (CC) approaches as well as more affordable ones with transition potential density functional theory (TP-DFT). A comparison of our computational results to previously reported experimental spectra shows that the core-valence separation (CVS) approaches CVS-ADC(2)-x and CVS-CC3 generally yield a good qualitative level of agreement with the experiment, suggesting that they can be used for interpreting measured spectra, while the TP-DFT method is not reliable for these molecules. However, quantitative agreement with the experiment is still outside the reach of the computational methods studied in this work.
Collapse
Affiliation(s)
- Goranka Bilalbegović
- Department
of Physics, Faculty of Science, University
of Zagreb, Bijenička
cesta 32, 10000 Zagreb, Croatia
| | - Aleksandar Maksimović
- Center
of Excellence for Advanced Materials and Sensing Devices, Rudjer Bošković Institute, Bijenička cesta 54, 10000 Zagreb, Croatia
| | - Lynne A. Valencic
- NASA
Goddard Space Flight Center, Greenbelt, 20771 Maryland, United States
- Department
of Physics & Astronomy, The Johns Hopkins
University, 366 Bloomberg
Center, 3400 N. Charles Street, Baltimore, 21218 Maryland, United States
| | - Susi Lehtola
- Department
of Chemistry, University of Helsinki, P.O. Box 55, A. I. Virtasen aukio
1, FI-00014 Helsinki, Finland
| |
Collapse
|
24
|
Fransson T, Brumboiu IE, Vidal ML, Norman P, Coriani S, Dreuw A. XABOOM: An X-ray Absorption Benchmark of Organic Molecules Based on Carbon, Nitrogen, and Oxygen 1s → π* Transitions. J Chem Theory Comput 2021; 17:1618-1637. [PMID: 33544612 PMCID: PMC8023667 DOI: 10.1021/acs.jctc.0c01082] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2020] [Indexed: 01/05/2023]
Abstract
The performance of several standard and popular approaches for calculating X-ray absorption spectra at the carbon, nitrogen, and oxygen K-edges of 40 primarily organic molecules up to the size of guanine has been evaluated, focusing on the low-energy and intense 1s → π* transitions. Using results obtained with CVS-ADC(2)-x and fc-CVS-EOM-CCSD as benchmark references, we investigate the performance of CC2, ADC(2), ADC(3/2), and commonly adopted density functional theory (DFT)-based approaches. Here, focus is on precision rather than on accuracy of transition energies and intensities-in other words, we target relative energies and intensities and the spread thereof, rather than absolute values. The use of exchange-correlation functionals tailored for time-dependent DFT calculations of core excitations leads to error spreads similar to those seen for more standard functionals, despite yielding superior absolute energies. Long-range corrected functionals are shown to perform particularly well compared to our reference data, showing error spreads in energy and intensity of 0.2-0.3 eV and ∼10%, respectively, as compared to 0.3-0.6 eV and ∼20% for a typical pure hybrid. In comparing intensities, state mixing can complicate matters, and techniques to avoid this issue are discussed. Furthermore, the influence of basis sets in high-level ab initio calculations is investigated, showing that reasonably accurate results are obtained with the use of 6-311++G**. We name this benchmark suite as XABOOM (X-ray absorption benchmark of organic molecules) and provide molecular structures and ground-state self-consistent field energies and spectroscopic data. We believe that it provides a good assessment of electronic structure theory methods for calculating X-ray absorption spectra and will become useful for future developments in this field.
Collapse
Affiliation(s)
- Thomas Fransson
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
- Fysikum, Stockholm University, Albanova, 10691 Stockholm, Sweden
| | - Iulia E. Brumboiu
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
- Department
of Chemistry, Korea Advanced Institute of
Science and Technology, 34141 Daejeon, Korea
| | - Marta L. Vidal
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
| | - Patrick Norman
- Department
of Theoretical Chemistry and Biology, KTH
Royal Institute of Technology, 10691 Stockholm, Sweden
| | - Sonia Coriani
- DTU
Chemistry, Technical University of Denmark, Kemitorvet Bldg 207, DK-2800 Kongens Lyngby, Denmark
- Department
of Chemistry, NTNU-Norwegian University
of Science and Technology, N-7991 Trondheim, Norway
| | - Andreas Dreuw
- Interdisciplinary
Center for Scientific Computing, Ruprecht-Karls
University, Im Neuenheimer
Feld 205, 69120 Heidelberg, Germany
| |
Collapse
|
25
|
Banerjee S, Sokolov AY. Efficient implementation of the single-reference algebraic diagrammatic construction theory for charged excitations: Applications to the TEMPO radical and DNA base pairs. J Chem Phys 2021; 154:074105. [PMID: 33607870 DOI: 10.1063/5.0040317] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/16/2022] Open
Abstract
We present an efficient implementation of the second- and third-order single-reference algebraic diagrammatic construction (ADC) theory for electron attachment and ionization energies and spectra [EA/IP-ADC(n), n = 2, 3]. Our new EA/IP-ADC program features spin adaptation for closed-shell systems, density fitting for efficient handling of the two-electron integral tensors, and vectorized and parallel implementation of tensor contractions. We demonstrate capabilities of our efficient implementation by applying the EA/IP-ADC(n) (n = 2, 3) methods to compute the photoelectron spectrum of the (2,2,6,6-tetramethylpiperidin-1-yl)oxyl (TEMPO) radical, as well as the vertical and adiabatic electron affinities of TEMPO and two DNA base pairs (guanine-cytosine and adenine-thymine). The spectra and electron affinities computed using large diffuse basis sets with up to 1028 molecular orbitals are found to be in good agreement with the best available results from the experiment and theoretical simulations.
Collapse
Affiliation(s)
- Samragni Banerjee
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Alexander Yu Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| |
Collapse
|
26
|
Backhouse OJ, Booth GH. Efficient Excitations and Spectra within a Perturbative Renormalization Approach. J Chem Theory Comput 2020; 16:6294-6304. [PMID: 32886508 DOI: 10.1021/acs.jctc.0c00701] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We present a self-consistent approach for computing the correlated quasiparticle spectrum of charged excitations in iterative O[N5] computational time. This is based on the auxiliary second-order Green's function approach [Backhouse, O. J. Chem. Theory Comput., 2000], in which a self-consistent effective Hamiltonian is constructed by systematically renormalizing the dynamical effects of the self-energy at second-order perturbation theory. From extensive benchmarking across the W4-11 molecular test set, we show that the iterative renormalization and truncation of the effective dynamical resolution arising from the 2h1p and 1h2p spaces can substantially improve the quality of the resulting ionization potential and electron affinity predictions compared to benchmark values. The resulting method is shown to be superior in accuracy to similarly scaling quantum chemical methods for charged excitations in EOM-CC2 and ADC(2), across this test set, while the self-consistency also removes the dependence on the underlying mean-field reference. The approach also allows for single-shot computation of the entire quasiparticle spectrum, which is applied to the benzoquinone molecule and demonstrates the reduction in the single-particle gap due to the correlated physics, and gives direct access to the localization of the Dyson orbitals.
Collapse
Affiliation(s)
- Oliver J Backhouse
- Department of Physics, King's College London, Strand, London WC2R 2LS, U.K
| | - George H Booth
- Department of Physics, King's College London, Strand, London WC2R 2LS, U.K
| |
Collapse
|
27
|
Schwalbe S, Fiedler L, Kraus J, Kortus J, Trepte K, Lehtola S. PyFLOSIC: Python-based Fermi–Löwdin orbital self-interaction correction. J Chem Phys 2020; 153:084104. [DOI: 10.1063/5.0012519] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Sebastian Schwalbe
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Lenz Fiedler
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Jakob Kraus
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Jens Kortus
- Institute of Theoretical Physics, TU Bergakademie Freiberg, Leipziger Str. 23, D-09599 Freiberg, Germany
| | - Kai Trepte
- Department of Physics, Central Michigan University, Mount Pleasant, Michigan 48859, USA
| | - Susi Lehtola
- Department of Chemistry, University of Helsinki, P.O. Box 55 (A. I. Virtasen Aukio 1), FI-00014 University of Helsinki, Finland
| |
Collapse
|
28
|
Scheurer M, Fransson T, Norman P, Dreuw A, Rehn DR. Complex excited state polarizabilities in the ADC/ISR framework. J Chem Phys 2020; 153:074112. [DOI: 10.1063/5.0012120] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Affiliation(s)
- Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| | - Thomas Fransson
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
- Department of Physics, AlbaNova University Center, Stockholm University, Stockholm SE-106 91, Sweden
| | - Patrick Norman
- Department of Theoretical Chemistry and Biology, School of Engineering Sciences in Chemistry, Biotechnology and Health, KTH Royal Institute of Technology, SE-106 91 Stockholm, Sweden
| | - Andreas Dreuw
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| | - Dirk R. Rehn
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University, Im Neuenheimer Feld 205, D-69120 Heidelberg, Germany
| |
Collapse
|
29
|
Herbst MF, Fransson T. Quantifying the error of the core-valence separation approximation. J Chem Phys 2020; 153:054114. [PMID: 32770930 DOI: 10.1063/5.0013538] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
For the calculation of core-excited states probed through X-ray absorption spectroscopy, the core-valence separation (CVS) scheme has become a vital tool. This approach allows us to target such states with high specificity, albeit introducing an error. We report the implementation of a post-processing step for CVS excitations obtained within the algebraic-diagrammatic construction scheme for the polarization propagator, which removes this error. Based on this, we provide a detailed analysis of the CVS scheme, identifying its accuracy to be dominated by an error balance between two neglected couplings, one between core and valence single excitations and the other between single and double core excitations. The selection of the basis set is shown to be vital for a proper description of both couplings, with tight polarizing functions being necessary for a good balance of errors. The CVS error is confirmed to be stable across multiple systems, with an element-specific spread for K-edge spectrum calculations of only about ±0.02 eV. A systematic lowering of the CVS error by 0.02 eV-0.03 eV is noted when considering excitations to extremely diffuse states, emulating ionization.
Collapse
Affiliation(s)
- Michael F Herbst
- CERMICS, École des Ponts ParisTech, 6-8 Avenue Blaise Pascal, 77455 Marne-la-Vallée, France; Inria Paris, 75589 Paris Cedex 12, France; and Sorbonne Universitée, Institut des Sciences du Calcul et des Données, ISCD, 75005 Paris, France
| | - Thomas Fransson
- Interdisciplinary Center for Scientific Computing, Heidelberg University, 69120 Heidelberg, Germany and Fysikum, Stockholm University, Albanova, 10691 Stockholm, Sweden
| |
Collapse
|
30
|
Peng C, Lewis CA, Wang X, Clement MC, Pierce K, Rishi V, Pavošević F, Slattery S, Zhang J, Teke N, Kumar A, Masteran C, Asadchev A, Calvin JA, Valeev EF. Massively Parallel Quantum Chemistry: A high-performance research platform for electronic structure. J Chem Phys 2020; 153:044120. [DOI: 10.1063/5.0005889] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Chong Peng
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Cannada A. Lewis
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Xiao Wang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | | | - Karl Pierce
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Varun Rishi
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Fabijan Pavošević
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Samuel Slattery
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Jinmei Zhang
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Nakul Teke
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Ashutosh Kumar
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Conner Masteran
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Andrey Asadchev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Justus A. Calvin
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| | - Edward F. Valeev
- Department of Chemistry, Virginia Tech, Blacksburg, Virginia 24061, USA
| |
Collapse
|
31
|
Sun Q, Zhang X, Banerjee S, Bao P, Barbry M, Blunt NS, Bogdanov NA, Booth GH, Chen J, Cui ZH, Eriksen JJ, Gao Y, Guo S, Hermann J, Hermes MR, Koh K, Koval P, Lehtola S, Li Z, Liu J, Mardirossian N, McClain JD, Motta M, Mussard B, Pham HQ, Pulkin A, Purwanto W, Robinson PJ, Ronca E, Sayfutyarova ER, Scheurer M, Schurkus HF, Smith JET, Sun C, Sun SN, Upadhyay S, Wagner LK, Wang X, White A, Whitfield JD, Williamson MJ, Wouters S, Yang J, Yu JM, Zhu T, Berkelbach TC, Sharma S, Sokolov AY, Chan GKL. Recent developments in the PySCF program package. J Chem Phys 2020; 153:024109. [DOI: 10.1063/5.0006074] [Citation(s) in RCA: 151] [Impact Index Per Article: 37.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023] Open
Affiliation(s)
- Qiming Sun
- AxiomQuant Investment Management LLC, Shanghai 200120, China
| | - Xing Zhang
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Samragni Banerjee
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Peng Bao
- Beijing National Laboratory for Molecular Sciences, State Key Laboratory for Structural Chemistry of Unstable and Stable Species, Institute of Chemistry, Chinese Academy of Sciences, Beijing 100190, China
| | - Marc Barbry
- Simbeyond B.V., P.O. Box 513, NL-5600 MB, Eindhoven, The Netherlands
| | - Nick S. Blunt
- Department of Chemistry, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | - Nikolay A. Bogdanov
- Max Planck Institute for Solid State Research, Heisenbergstraße 1, 70569 Stuttgart, Germany
| | - George H. Booth
- Department of Physics, King’s College London, Strand, London WC2R 2LS, United Kingdom
| | - Jia Chen
- Department of Physics, University of Florida, Gainesville, Florida 32611, USA
- Quantum Theory Project, University of Florida, Gainesville, Florida 32611, USA
| | - Zhi-Hao Cui
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Janus J. Eriksen
- School of Chemistry, University of Bristol, Cantock’s Close, Bristol BS8 1TS, United Kingdom
| | - Yang Gao
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, California 91125, USA
| | - Sheng Guo
- Google Inc., Mountain View, California 94043, USA
| | - Jan Hermann
- FU Berlin, Department of Mathematics and Computer Science, Arnimallee 6, 14195 Berlin, Germany
- TU Berlin, Machine Learning Group, Marchstr. 23, 10587 Berlin, Germany
| | - Matthew R. Hermes
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455, USA
| | - Kevin Koh
- Department of Chemistry and Biochemistry, The University of Notre Dame du Lac, 251 Nieuwland Science Hall, Notre Dame, Indiana 46556, USA
| | - Peter Koval
- Simune Atomistics S.L., Avenida Tolosa 76, Donostia-San Sebastian, Spain
| | - Susi Lehtola
- Department of Chemistry, University of Helsinki, P.O. Box 55 (A. I. Virtasen aukio 1), FI-00014 Helsinki, Finland
| | - Zhendong Li
- Key Laboratory of Theoretical and Computational Photochemistry, Ministry of Education, College of Chemistry, Beijing Normal University, Beijing 100875, China
| | - Junzi Liu
- Department of Chemistry, The Johns Hopkins University, Baltimore, Maryland 21218, USA
| | - Narbe Mardirossian
- AMGEN Research, One Amgen Center Drive, Thousand Oaks, California 91320, USA
| | | | - Mario Motta
- IBM Almaden Research Center, San Jose, California 95120, USA
| | - Bastien Mussard
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Hung Q. Pham
- Department of Chemistry, Chemical Theory Center, and Supercomputing Institute, University of Minnesota, 207 Pleasant Street SE, Minneapolis, Minnesota 55455, USA
| | - Artem Pulkin
- QuTech and Kavli Institute of Nanoscience, Delft University of Technology, The Netherlands
| | - Wirawan Purwanto
- Information Technology Services, Old Dominion University, Norfolk, Virginia 23529, USA
| | - Paul J. Robinson
- Department of Chemistry, Columbia University, New York, New York 10027, USA
| | - Enrico Ronca
- Istituto per i Processi Chimico Fisici del CNR (IPCF-CNR), Via G. Moruzzi, 1, 56124 Pisa, Italy
| | - Elvira R. Sayfutyarova
- Department of Chemistry, Yale University, 225 Prospect Street, New Haven, Connecticut 06520, USA
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific Computing, Ruprecht-Karls University of Heidelberg, 205 Im Neuenheimer Feld, 69120 Heidelberg, Germany
| | - Henry F. Schurkus
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - James E. T. Smith
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Chong Sun
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Shi-Ning Sun
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, California 91125, USA
| | - Shiv Upadhyay
- Department of Chemistry, University of Pittsburgh, Pittsburgh, Pennsylvania 15260, USA
| | - Lucas K. Wagner
- Department of Physics and Institute for Condensed Matter Theory, University of Illinois at Urbana-Champaign, Illinois 61801, USA
| | - Xiao Wang
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| | - Alec White
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - James Daniel Whitfield
- Department of Physics and Astronomy, Dartmouth College, Hanover, New Hampshire 03755, USA
| | - Mark J. Williamson
- Department of Chemistry, University of Cambridge, Lensfield Road, Cambridge CB2 1EW, United Kingdom
| | | | - Jun Yang
- Department of Chemistry, The University of Hong Kong, Pokfulam Road, Hong Kong SAR, China
| | - Jason M. Yu
- Department of Chemistry, University of California, Irvine, 1102 Natural Sciences II, Irvine, California 92697-2025, USA
| | - Tianyu Zhu
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| | - Timothy C. Berkelbach
- Department of Chemistry, Columbia University, New York, New York 10027, USA
- Center for Computational Quantum Physics, Flatiron Institute, New York, New York 10010, USA
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado, Boulder, Colorado 80302, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The Ohio State University, Columbus, Ohio 43210, USA
| | - Garnet Kin-Lic Chan
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, California 91125, USA
| |
Collapse
|
32
|
Smith DGA, Burns LA, Simmonett AC, Parrish RM, Schieber MC, Galvelis R, Kraus P, Kruse H, Di Remigio R, Alenaizan A, James AM, Lehtola S, Misiewicz JP, Scheurer M, Shaw RA, Schriber JB, Xie Y, Glick ZL, Sirianni DA, O’Brien JS, Waldrop JM, Kumar A, Hohenstein EG, Pritchard BP, Brooks BR, Schaefer HF, Sokolov AY, Patkowski K, DePrince AE, Bozkaya U, King RA, Evangelista FA, Turney JM, Crawford TD, Sherrill CD. Psi4 1.4: Open-source software for high-throughput quantum chemistry. J Chem Phys 2020; 152:184108. [PMID: 32414239 PMCID: PMC7228781 DOI: 10.1063/5.0006002] [Citation(s) in RCA: 337] [Impact Index Per Article: 84.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2020] [Accepted: 04/12/2020] [Indexed: 12/13/2022] Open
Abstract
PSI4 is a free and open-source ab initio electronic structure program providing implementations of Hartree-Fock, density functional theory, many-body perturbation theory, configuration interaction, density cumulant theory, symmetry-adapted perturbation theory, and coupled-cluster theory. Most of the methods are quite efficient, thanks to density fitting and multi-core parallelism. The program is a hybrid of C++ and Python, and calculations may be run with very simple text files or using the Python API, facilitating post-processing and complex workflows; method developers also have access to most of PSI4's core functionalities via Python. Job specification may be passed using The Molecular Sciences Software Institute (MolSSI) QCSCHEMA data format, facilitating interoperability. A rewrite of our top-level computation driver, and concomitant adoption of the MolSSI QCARCHIVE INFRASTRUCTURE project, makes the latest version of PSI4 well suited to distributed computation of large numbers of independent tasks. The project has fostered the development of independent software components that may be reused in other quantum chemistry programs.
Collapse
Affiliation(s)
| | - Lori A. Burns
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew C. Simmonett
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Robert M. Parrish
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Matthew C. Schieber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | | | - Peter Kraus
- School of Molecular and Life Sciences, Curtin
University, Kent St., Bentley, Perth, Western Australia 6102,
Australia
| | - Holger Kruse
- Institute of Biophysics of the Czech Academy of
Sciences, Královopolská 135, 612 65 Brno, Czech
Republic
| | - Roberto Di Remigio
- Department of Chemistry, Centre for Theoretical
and Computational Chemistry, UiT, The Arctic University of Norway, N-9037
Tromsø, Norway
| | - Asem Alenaizan
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Andrew M. James
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Susi Lehtola
- Department of Chemistry, University of
Helsinki, P.O. Box 55 (A. I. Virtasen aukio 1), FI-00014 Helsinki,
Finland
| | - Jonathon P. Misiewicz
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Maximilian Scheurer
- Interdisciplinary Center for Scientific
Computing, Heidelberg University, D-69120 Heidelberg,
Germany
| | - Robert A. Shaw
- ARC Centre of Excellence in Exciton Science,
School of Science, RMIT University, Melbourne, VIC 3000,
Australia
| | - Jeffrey B. Schriber
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Yi Xie
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Zachary L. Glick
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Dominic A. Sirianni
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Joseph Senan O’Brien
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| | - Jonathan M. Waldrop
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - Ashutosh Kumar
- Department of Chemistry, Virginia
Tech, Blacksburg, Virginia 24061, USA
| | - Edward G. Hohenstein
- SLAC National Accelerator Laboratory, Stanford
PULSE Institute, Menlo Park, California 94025,
USA
| | | | - Bernard R. Brooks
- National Institutes of Health – National Heart,
Lung and Blood Institute, Laboratory of Computational Biology, Bethesda,
Maryland 20892, USA
| | - Henry F. Schaefer
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | - Alexander Yu. Sokolov
- Department of Chemistry and Biochemistry, The
Ohio State University, Columbus, Ohio 43210, USA
| | - Konrad Patkowski
- Department of Chemistry and Biochemistry, Auburn
University, Auburn, Alabama 36849, USA
| | - A. Eugene DePrince
- Department of Chemistry and Biochemistry,
Florida State University, Tallahassee, Florida 32306-4390,
USA
| | - Uğur Bozkaya
- Department of Chemistry, Hacettepe
University, Ankara 06800, Turkey
| | - Rollin A. King
- Department of Chemistry, Bethel
University, St. Paul, Minnesota 55112, USA
| | | | - Justin M. Turney
- Center for Computational Quantum Chemistry,
University of Georgia, Athens, Georgia 30602, USA
| | | | - C. David Sherrill
- Center for Computational Molecular Science and
Technology, School of Chemistry and Biochemistry, School of Computational Science and
Engineering, Georgia Institute of Technology, Atlanta, Georgia 30332-0400,
USA
| |
Collapse
|