1
|
Pototschnig JV, Papadopoulos A, Lyakh DI, Repisky M, Halbert L, Severo Pereira Gomes A, Jensen HJA, Visscher L. Implementation of Relativistic Coupled Cluster Theory for Massively Parallel GPU-Accelerated Computing Architectures. J Chem Theory Comput 2021; 17:5509-5529. [PMID: 34370471 PMCID: PMC8444343 DOI: 10.1021/acs.jctc.1c00260] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
In this paper, we report reimplementation of the core algorithms of relativistic coupled cluster theory aimed at modern heterogeneous high-performance computational infrastructures. The code is designed for parallel execution on many compute nodes with optional GPU coprocessing, accomplished via the new ExaTENSOR back end. The resulting ExaCorr module is primarily intended for calculations of molecules with one or more heavy elements, as relativistic effects on the electronic structure are included from the outset. In the current work, we thereby focus on exact two-component methods and demonstrate the accuracy and performance of the software. The module can be used as a stand-alone program requiring a set of molecular orbital coefficients as the starting point, but it is also interfaced to the DIRAC program that can be used to generate these. We therefore also briefly discuss an improvement of the parallel computing aspects of the relativistic self-consistent field algorithm of the DIRAC program.
Collapse
Affiliation(s)
- Johann V Pototschnig
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, de Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Anastasios Papadopoulos
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, de Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Dmitry I Lyakh
- National Center for Computational Sciences, Oak Ridge National Laboratory, Oak Ridge, Tennessee 37831, United States
| | - Michal Repisky
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Loïc Halbert
- Universite de Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers, Atomes et Molecules, F-59000 Lille, France
| | - André Severo Pereira Gomes
- Universite de Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers, Atomes et Molecules, F-59000 Lille, France
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, de Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| |
Collapse
|
2
|
Gaul K, Berger R. Quasi-relativistic study of nuclear electric quadrupole coupling constants in chiral molecules containing heavy elements. Mol Phys 2020. [DOI: 10.1080/00268976.2020.1797199] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Affiliation(s)
- Konstantin Gaul
- Fachbereich Chemie, Philipps-Universität Marburg, Marburg, Germany
| | - Robert Berger
- Fachbereich Chemie, Philipps-Universität Marburg, Marburg, Germany
| |
Collapse
|
3
|
Saue T, Bast R, Gomes ASP, Jensen HJA, Visscher L, Aucar IA, Di Remigio R, Dyall KG, Eliav E, Fasshauer E, Fleig T, Halbert L, Hedegård ED, Helmich-Paris B, Iliaš M, Jacob CR, Knecht S, Laerdahl JK, Vidal ML, Nayak MK, Olejniczak M, Olsen JMH, Pernpointner M, Senjean B, Shee A, Sunaga A, van Stralen JNP. The DIRAC code for relativistic molecular calculations. J Chem Phys 2020; 152:204104. [PMID: 32486677 DOI: 10.1063/5.0004844] [Citation(s) in RCA: 122] [Impact Index Per Article: 30.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
DIRAC is a freely distributed general-purpose program system for one-, two-, and four-component relativistic molecular calculations at the level of Hartree-Fock, Kohn-Sham (including range-separated theory), multiconfigurational self-consistent-field, multireference configuration interaction, electron propagator, and various flavors of coupled cluster theory. At the self-consistent-field level, a highly original scheme, based on quaternion algebra, is implemented for the treatment of both spatial and time reversal symmetry. DIRAC features a very general module for the calculation of molecular properties that to a large extent may be defined by the user and further analyzed through a powerful visualization module. It allows for the inclusion of environmental effects through three different classes of increasingly sophisticated embedding approaches: the implicit solvation polarizable continuum model, the explicit polarizable embedding model, and the frozen density embedding model.
Collapse
Affiliation(s)
- Trond Saue
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Radovan Bast
- Department of Information Technology, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - André Severo Pereira Gomes
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Hans Jørgen Aa Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, DK-5230 Odense M, Denmark
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Ignacio Agustín Aucar
- Instituto de Modelado e Innovación Tecnológica, CONICET, and Departamento de Física-Facultad de Ciencias Exactas y Naturales, UNNE, Avda. Libertad 5460, W3404AAS Corrientes, Argentina
| | - Roberto Di Remigio
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Kenneth G Dyall
- Dirac Solutions, 10527 NW Lost Park Drive, Portland, Oregon 97229, USA
| | - Ephraim Eliav
- School of Chemistry, Tel Aviv University, Ramat Aviv, Tel Aviv 69978, Israel
| | - Elke Fasshauer
- Department of Physics and Astronomy, Aarhus University, Ny Munkegade 120, 8000 Aarhus, Denmark
| | - Timo Fleig
- Laboratoire de Chimie et Physique Quantique, UMR 5626 CNRS-Université Toulouse III-Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse, France
| | - Loïc Halbert
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, F-59000 Lille, France
| | - Erik Donovan Hedegård
- Division of Theoretical Chemistry, Lund University, Chemical Centre, P.O. Box 124, SE-221 00 Lund, Sweden
| | - Benjamin Helmich-Paris
- Max-Planck-Institut für Kohlenforschung, Kaiser-Wilhelm-Platz 1, 45470 Mülheim an der Ruhr, Germany
| | - Miroslav Iliaš
- Department of Chemistry, Faculty of Natural Sciences, Matej Bel University, Tajovského 40, 974 01 Banská Bystrica, Slovakia
| | - Christoph R Jacob
- Technische Universität Braunschweig, Institute of Physical and Theoretical Chemistry, Gaußstr. 17, 38106 Braunschweig, Germany
| | - Stefan Knecht
- ETH Zürich, Laboratorium für Physikalische Chemie, Vladimir-Prelog-Weg 2, 8093 Zürich, Switzerland
| | - Jon K Laerdahl
- Department of Microbiology, Oslo University Hospital, Oslo, Norway
| | - Marta L Vidal
- Department of Chemistry, Technical University of Denmark, 2800 Kgs. Lyngby, Denmark
| | - Malaya K Nayak
- Theoretical Chemistry Section, Bhabha Atomic Research Centre, Trombay, Mumbai 400085, India
| | - Małgorzata Olejniczak
- Centre of New Technologies, University of Warsaw, S. Banacha 2c, 02-097 Warsaw, Poland
| | - Jógvan Magnus Haugaard Olsen
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT The Arctic University of Norway, N-9037 Tromsø, Norway
| | | | - Bruno Senjean
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, USA
| | - Ayaki Sunaga
- Department of Chemistry, Tokyo Metropolitan University, 1-1 Minami-Osawa, Hachioji-city, Tokyo 192-0397, Japan
| | - Joost N P van Stralen
- Department of Chemistry and Pharmaceutical Sciences, Vrije Universiteit Amsterdam, NL-1081HV Amsterdam, The Netherlands
| |
Collapse
|
4
|
Affiliation(s)
- Pekka Pyykkö
- Department of Chemistry, University of Helsinki, Helsinki, Finland
| |
Collapse
|
5
|
Santiago RT, Haiduke RLA. New density functional parameterizations to accurate calculations of electric field gradient variations among compounds. J Comput Chem 2015; 36:2125-30. [PMID: 26284820 DOI: 10.1002/jcc.24052] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2015] [Revised: 06/15/2015] [Accepted: 07/30/2015] [Indexed: 11/09/2022]
Abstract
This research provides a performance investigation of density functional theory and also proposes new functional parameterizations to deal with electric field gradient (EFG) calculations at nuclear positions. The entire procedure is conducted within the four-component formalism. First, we noticed that traditional hybrid and long-range corrected functionals are more efficient in the description of EFG variations for a set of elements (indium, antimony, iodine, lutetium, and hafnium) among linear molecules. Thus, we selected the PBE0, B3LYP, and CAM-B3LYP functionals and promoted a reoptimization of their parameters for a better description of these EFG changes. The PBE0q variant developed here showed an overall promising performance in a validation test conducted with potassium, iodine, copper, and gold. In general, the correlation coefficients found in linear regressions between experimental nuclear quadrupole coupling constants and calculated EFGs are improved while the systematic EFG errors also decrease as a result of this reparameterization.
Collapse
Affiliation(s)
- Régis Tadeu Santiago
- Departamento De Química E Física Molecular, Instituto De Química De São Carlos, Universidade De São Paulo, Avenida Trabalhador São-Carlense, 400, CP 780, São Carlos, SP, 13560-970, Brazil
| | - Roberto Luiz Andrade Haiduke
- Departamento De Química E Física Molecular, Instituto De Química De São Carlos, Universidade De São Paulo, Avenida Trabalhador São-Carlense, 400, CP 780, São Carlos, SP, 13560-970, Brazil
| |
Collapse
|
6
|
Santiago RT, Teodoro TQ, Haiduke RLA. The nuclear electric quadrupole moment of copper. Phys Chem Chem Phys 2014; 16:11590-6. [PMID: 24806277 DOI: 10.1039/c4cp00706a] [Citation(s) in RCA: 9] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
Abstract
The nuclear electric quadrupole moment (NQM) of the (63)Cu nucleus was determined from an indirect approach by combining accurate experimental nuclear quadrupole coupling constants (NQCCs) with relativistic Dirac-Coulomb coupled cluster calculations of the electric field gradient (EFG). The data obtained at the highest level of calculation, DC-CCSD-T, from 14 linear molecules containing the copper atom give rise to an indicated NQM of -198(10) mbarn. Such result slightly deviates from the previously accepted standard value given by the muonic method, -220(15) mbarn, although the error bars are superimposed.
Collapse
Affiliation(s)
- Régis Tadeu Santiago
- Departamento de Química e Física Molecular, Instituto de Química de São Carlos, Universidade de São Paulo, Av. Trabalhador São-carlense, 400-CP 780, 13560-970-São Carlos, SP, Brazil.
| | | | | |
Collapse
|
7
|
Stopkowicz S, Cheng L, Harding ME, Puzzarini C, Gauss J. The bromine nuclear quadrupole moment revisited. Mol Phys 2013. [DOI: 10.1080/00268976.2013.796072] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
Affiliation(s)
- Stella Stopkowicz
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
- b Centre for Theoretical and Computational Chemistry, Department of Chemistry , University of Oslo , Oslo , N-0315 , Norway
| | - Lan Cheng
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
- c Institute for Theoretical Chemistry, Department of Chemistry and Biochemistry , The University of Texas at Austin , Austin , TX , 78712 , USA
| | - Michael E. Harding
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
- d Institut für Nanotechnologie, Karlsruher Institut für Technologie (KIT) , Karlsruhe , D-76021 , Germany
| | - Cristina Puzzarini
- e Dipartimento di Chimica ‘Giacomo Ciamician’ , Università di Bologna , Via Selmi 2, Bologna , I-40126 , Italy
| | - Jürgen Gauss
- a Institut für Physikalische Chemie , Universität Mainz , Mainz , D-55099 , Germany
| |
Collapse
|
8
|
Filatov M, Zou W, Cremer D. Relativistically corrected electric field gradients calculated with the normalized elimination of the small component formalism. J Chem Phys 2012; 137:054113. [DOI: 10.1063/1.4742175] [Citation(s) in RCA: 20] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
|
9
|
|
10
|
Arcisauskaite V, Knecht S, Sauer SPA, Hemmingsen L. Fully relativistic coupled cluster and DFT study of electric field gradients at Hg in 199Hg compounds. Phys Chem Chem Phys 2012; 14:2651-7. [DOI: 10.1039/c2cp23080a] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022]
|
11
|
Styszyński J. Why do we Need Relativistic Computational Methods? CHALLENGES AND ADVANCES IN COMPUTATIONAL CHEMISTRY AND PHYSICS 2010. [DOI: 10.1007/978-1-4020-9975-5_3] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/30/2023]
|
12
|
Abstract
Electric field gradients at the nuclei of gallim and indium are determined by finite field calculations of the atomic energies as functions of the nuclear quadrupole moments. The four-component Dirac–Coulomb–Gaunt Hamiltonian serves as framework, and all electrons are correlated by Fock-space coupled cluster with single and double excitations or by single reference coupled cluster with approximate triples. Large, converged basis sets (e.g., 28s24p20d13f5g4h for In) and virtual spaces are used. Together with experimental nuclear quadrupole coupling constants, known with high precision, the calculated electric field gradients yield the nuclear quadrupole moments. For 69Ga, we get Q = 174(3) mb, in agreement with the earlier 171(2) mb obtained from molecular calculations. The 115In moment is Q = 772(5) mb, considerably lower than the previously accepted 810 mb, and in good agreement with the recent molecular value of 770(8) mb.
Collapse
Affiliation(s)
- Hana Yakobi
- School of Chemistry, Tel Aviv University, 69978 Tel Aviv, Israel
| | - Ephraim Eliav
- School of Chemistry, Tel Aviv University, 69978 Tel Aviv, Israel
| | - Uzi Kaldor
- School of Chemistry, Tel Aviv University, 69978 Tel Aviv, Israel
| |
Collapse
|
13
|
|
14
|
|
15
|
Kellö V, Sadlej AJ. The Nuclear Quadrupole Moment of 14N from Accurate Electric Field Gradient Calculations and Microwave Spectra of NP Molecule. ACTA ACUST UNITED AC 2007. [DOI: 10.1135/cccc20070064] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/05/2022]
Abstract
Extensive series of relativistic coupled cluster calculations of the electric field gradient at N in NP has been carried out. The accurate value of the calculated electric field gradient, combined with the highly accurate experimental value of the nuclear quadrupole coupling constant for the 14N nucleus, gives the 'molecular' value of the nuclear quadrupole moment Q(14N) = 20.46 mb. This result perfectly agrees with the value (20.44 ± 0.03 mb) determined from atomic calculations and atomic spectra. The present study involves also extensive investigations of basis sets which must be used in highly accurate calculations of electric field gradients.
Collapse
|
16
|
Haiduke RLA, da Silva ABF, Visscher L. The nuclear electric quadrupole moment of antimony from the molecular method. J Chem Phys 2006; 125:64301. [PMID: 16942280 DOI: 10.1063/1.2234369] [Citation(s) in RCA: 25] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Relativistic Dirac-Coulomb (DC) Hartree-Fock calculations are employed to obtain the analytic electric field gradient (EFG) on the antimony nucleus in the SbN, SbP, SbF, and SbCl molecules. The electronic correlation contribution to the EFGs is included with the DC-CCSD(T) and DC-CCSD-T approaches, also in the four-component framework, using a finite-difference method. The total EFG results, along with the experimental nuclear quadrupole coupling constants from microwave spectroscopy, allow to derive the nuclear quadrupole moments of (121)Sb and (123)Sb, respectively, as -543(11) and -692(14) mb.
Collapse
Affiliation(s)
- Roberto L A Haiduke
- Departamento de Química e Física Molecular, Instituto de Química de São Carlos, C.P. 780, São Carlos, São Paulo 13560-970, Brazil
| | | | | |
Collapse
|
17
|
Verbraak H, van Stralen JNP, Bouwman J, de Klerk JS, Verdes D, Linnartz H, Bickelhaupt FM. High-resolution infrared spectroscopy of the charge-transfer complex [Ar–N2]+∙: A combined experimental/theoretical study. J Chem Phys 2005; 123:144305. [PMID: 16238389 DOI: 10.1063/1.2039083] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
A combined experimental and theoretical study of the charge-transfer complex [Ar-N(2)](+) is presented. Nearly 50 transitions split by spin-rotation interaction have been observed by means of infrared diode laser absorption spectroscopy in a supersonic planar plasma expansion. The band origin is at 2272.2563(18) cm(-1) and rotational constants in the ground and vibrationally (NN-stretch) excited state amount to 0.128701(8) cm(-1) and 0.128203(8) cm(-1), respectively. The interpretation of the data in terms of a charge switch upon complexation is supported by new ab initio calculations. The best estimate for a linear equilibrium structure yields R(e)(NN)=1.102 A and R(e)(Ar-N)=2.190 A. Predictions for molecular parameters not directly available from the experimental results are presented as well. Furthermore, the electronic structure and Ar-N bonding mechanism of [Ar-N(2)](+) have been analyzed in detail. The Ar-N bond is a textbook example of a classical 2-center-3-electron bond.
Collapse
Affiliation(s)
- H Verbraak
- Laser Centre Amsterdam, Department of Physical Chemistry, Vrije Universiteit, Amsterdam, The Netherlands
| | | | | | | | | | | | | |
Collapse
|
18
|
van Stralen JN, Visscher L, Larsen CV, Jensen HJA. First-order MP2 molecular properties in a relativistic framework. Chem Phys 2005. [DOI: 10.1016/j.chemphys.2004.10.018] [Citation(s) in RCA: 28] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/26/2022]
|
19
|
Some Trends in Relativistic and Electron Correlation Effects in Electric Properties of Small Molecules. ADVANCES IN QUANTUM CHEMISTRY 2005. [DOI: 10.1016/s0065-3276(05)50012-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register]
|