1
|
Gamboni G, Belpassi L, Belanzoni P. The Chemical Bond at the Bottom of the Periodic Table: The Case of the Heavy Astatine and the Super-Heavy Tennessine. Chemphyschem 2024; 25:e202400310. [PMID: 38708605 DOI: 10.1002/cphc.202400310] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2024] [Accepted: 05/03/2024] [Indexed: 05/07/2024]
Abstract
In this work, we study the chemical bond in molecules containing heavy and super-heavy elements according to the current state-of-the-art bonding models. An Energy Decomposition Analysis in combination with Natural Orbital for Chemical Valence (EDA-NOCV) within the relativistic four-component Dirac-Kohn-Sham (DKS) framework is employed, which allows to successfully include the spin-orbit coupling (SOC) effects on the chemical bond description. Simple halogen-bonded adducts ClX⋯L (X=At, Ts; L=NH3, Br-, H2O, CO) of astatine and tennessine have been selected to assess a trend on descending along a group, while modulating the ClX⋯L bond features through the different electronic nature of the ligand L. Interesting effects caused by SOC have been revealed: i) a huge increase of the ClTs dipole moment (which is almost twice as that of ClAt), ii) a lowering of the ClX⋯L bonding energy arising from different contributions to the ClX…L interaction energy strongly depending on the nature of L, iii) a quenching of one of the π back-donation components to the bond. In the ClTs(CO) adduct, the back-donation from ClTs to CO becomes the most important component. The analysis of the electronic structure of the ClX dimers allows for a clear interpretation of the SOC effects in these systems.
Collapse
Affiliation(s)
- Giulia Gamboni
- Department of Chemistry, Biology and Biotechnology, University of Perugia, via Elce di Sotto 8, 06123, Perugia, Italy
| | - Leonardo Belpassi
- CNR Institute of Chemical Science and Technologies "Giulio Natta" (CNR-SCITEC), via Elce di Sotto 8, 06123, Perugia, Italy
| | - Paola Belanzoni
- Department of Chemistry, Biology and Biotechnology, University of Perugia, via Elce di Sotto 8, 06123, Perugia, Italy
- CNR Institute of Chemical Science and Technologies "Giulio Natta" (CNR-SCITEC), via Elce di Sotto 8, 06123, Perugia, Italy
| |
Collapse
|
2
|
Solov’yov AV, Verkhovtsev AV, Mason NJ, Amos RA, Bald I, Baldacchino G, Dromey B, Falk M, Fedor J, Gerhards L, Hausmann M, Hildenbrand G, Hrabovský M, Kadlec S, Kočišek J, Lépine F, Ming S, Nisbet A, Ricketts K, Sala L, Schlathölter T, Wheatley AEH, Solov’yov IA. Condensed Matter Systems Exposed to Radiation: Multiscale Theory, Simulations, and Experiment. Chem Rev 2024; 124:8014-8129. [PMID: 38842266 PMCID: PMC11240271 DOI: 10.1021/acs.chemrev.3c00902] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/06/2023] [Revised: 05/02/2024] [Accepted: 05/10/2024] [Indexed: 06/07/2024]
Abstract
This roadmap reviews the new, highly interdisciplinary research field studying the behavior of condensed matter systems exposed to radiation. The Review highlights several recent advances in the field and provides a roadmap for the development of the field over the next decade. Condensed matter systems exposed to radiation can be inorganic, organic, or biological, finite or infinite, composed of different molecular species or materials, exist in different phases, and operate under different thermodynamic conditions. Many of the key phenomena related to the behavior of irradiated systems are very similar and can be understood based on the same fundamental theoretical principles and computational approaches. The multiscale nature of such phenomena requires the quantitative description of the radiation-induced effects occurring at different spatial and temporal scales, ranging from the atomic to the macroscopic, and the interlinks between such descriptions. The multiscale nature of the effects and the similarity of their manifestation in systems of different origins necessarily bring together different disciplines, such as physics, chemistry, biology, materials science, nanoscience, and biomedical research, demonstrating the numerous interlinks and commonalities between them. This research field is highly relevant to many novel and emerging technologies and medical applications.
Collapse
Affiliation(s)
| | | | - Nigel J. Mason
- School
of Physics and Astronomy, University of
Kent, Canterbury CT2 7NH, United
Kingdom
| | - Richard A. Amos
- Department
of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K.
| | - Ilko Bald
- Institute
of Chemistry, University of Potsdam, Karl-Liebknecht-Str. 24-25, 14476 Potsdam, Germany
| | - Gérard Baldacchino
- Université
Paris-Saclay, CEA, LIDYL, 91191 Gif-sur-Yvette, France
- CY Cergy Paris Université,
CEA, LIDYL, 91191 Gif-sur-Yvette, France
| | - Brendan Dromey
- Centre
for Light Matter Interactions, School of Mathematics and Physics, Queen’s University Belfast, Belfast BT7 1NN, United Kingdom
| | - Martin Falk
- Institute
of Biophysics of the Czech Academy of Sciences, Královopolská 135, 61200 Brno, Czech Republic
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Juraj Fedor
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Luca Gerhards
- Institute
of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| | - Michael Hausmann
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
| | - Georg Hildenbrand
- Kirchhoff-Institute
for Physics, Heidelberg University, Im Neuenheimer Feld 227, 69120 Heidelberg, Germany
- Faculty
of Engineering, University of Applied Sciences
Aschaffenburg, Würzburger
Str. 45, 63743 Aschaffenburg, Germany
| | | | - Stanislav Kadlec
- Eaton European
Innovation Center, Bořivojova
2380, 25263 Roztoky, Czech Republic
| | - Jaroslav Kočišek
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Franck Lépine
- Université
Claude Bernard Lyon 1, CNRS, Institut Lumière
Matière, F-69622, Villeurbanne, France
| | - Siyi Ming
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Andrew Nisbet
- Department
of Medical Physics and Biomedical Engineering, University College London, London WC1E 6BT, U.K.
| | - Kate Ricketts
- Department
of Targeted Intervention, University College
London, Gower Street, London WC1E 6BT, United Kingdom
| | - Leo Sala
- J.
Heyrovský Institute of Physical Chemistry, Czech Academy of Sciences, Dolejškova 3, 18223 Prague, Czech Republic
| | - Thomas Schlathölter
- Zernike
Institute for Advanced Materials, University
of Groningen, Nijenborgh
4, 9747 AG Groningen, The Netherlands
- University
College Groningen, University of Groningen, Hoendiepskade 23/24, 9718 BG Groningen, The Netherlands
| | - Andrew E. H. Wheatley
- Yusuf
Hamied Department of Chemistry, University
of Cambridge, Lensfield
Road, Cambridge CB2 1EW, United Kingdom
| | - Ilia A. Solov’yov
- Institute
of Physics, Carl von Ossietzky University, Carl-von-Ossietzky-Str. 9-11, 26129 Oldenburg, Germany
| |
Collapse
|
3
|
Focke K, De Santis M, Wolter M, Martinez B JA, Vallet V, Pereira Gomes AS, Olejniczak M, Jacob CR. Interoperable workflows by exchanging grid-based data between quantum-chemical program packages. J Chem Phys 2024; 160:162503. [PMID: 38686818 DOI: 10.1063/5.0201701] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/31/2024] [Accepted: 04/02/2024] [Indexed: 05/02/2024] Open
Abstract
Quantum-chemical subsystem and embedding methods require complex workflows that may involve multiple quantum-chemical program packages. Moreover, such workflows require the exchange of voluminous data that go beyond simple quantities, such as molecular structures and energies. Here, we describe our approach for addressing this interoperability challenge by exchanging electron densities and embedding potentials as grid-based data. We describe the approach that we have implemented to this end in a dedicated code, PyEmbed, currently part of a Python scripting framework. We discuss how it has facilitated the development of quantum-chemical subsystem and embedding methods and highlight several applications that have been enabled by PyEmbed, including wave-function theory (WFT) in density-functional theory (DFT) embedding schemes mixing non-relativistic and relativistic electronic structure methods, real-time time-dependent DFT-in-DFT approaches, the density-based many-body expansion, and workflows including real-space data analysis and visualization. Our approach demonstrates, in particular, the merits of exchanging (complex) grid-based data and, in general, the potential of modular software development in quantum chemistry, which hinges upon libraries that facilitate interoperability.
Collapse
Affiliation(s)
- Kevin Focke
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Matteo De Santis
- CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, Univ. Lille, F-59000 Lille, France
| | - Mario Wolter
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| | - Jessica A Martinez B
- CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, Univ. Lille, F-59000 Lille, France
- Department of Chemistry, Rutgers University, Newark, New Jersey 07102, USA
| | - Valérie Vallet
- CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, Univ. Lille, F-59000 Lille, France
| | | | - Małgorzata Olejniczak
- Centre of New Technologies, University of Warsaw, S. Banacha 2c, 02-097 Warsaw, Poland
| | - Christoph R Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstraße 17, 38106 Braunschweig, Germany
| |
Collapse
|
4
|
Brakestad A, Jensen SR, Tantardini C, Pitteloud Q, Wind P, Užulis J, Gulans A, Hopmann KH, Frediani L. Scalar Relativistic Effects with Multiwavelets: Implementation and Benchmark. J Chem Theory Comput 2024; 20:728-737. [PMID: 38181377 PMCID: PMC10809714 DOI: 10.1021/acs.jctc.3c01095] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2023] [Revised: 12/18/2023] [Accepted: 12/19/2023] [Indexed: 01/07/2024]
Abstract
The importance of relativistic effects in quantum chemistry is widely recognized, not only for heavier elements but throughout the periodic table. At the same time, relativistic effects are strongest in the nuclear region, where the description of electrons through a linear combination of atomic orbitals becomes more challenging. Furthermore, the choice of basis sets for heavier elements is limited compared with lighter elements where precise basis sets are available. Thanks to the framework of multiresolution analysis, multiwavelets provide an appealing alternative to overcoming this challenge: they lead to robust error control and adaptive algorithms that automatically refine the basis set description until the desired precision is reached. This allows one to achieve a proper description of the nuclear region. In this work, we extended the multiwavelet-based code MRChem to the scalar zero-order regular approximation framework. We validated our implementation by comparing the total energies for a small set of elements and molecules. To confirm the validity of our implementation, we compared both against a radial numerical code for atoms and the plane-wave-based code EXCITING.
Collapse
Affiliation(s)
- Anders Brakestad
- Hylleraas
Centre for Quantum Molecular Sciences, UiT
The Arctic University of Norway, Tromsø 9037, Norway
- Department
of Chemistry, UiT The Arctic University
of Norway, Tromsø 9037, Norway
| | - Stig Rune Jensen
- Hylleraas
Centre for Quantum Molecular Sciences, UiT
The Arctic University of Norway, Tromsø 9037, Norway
- Department
of Chemistry, UiT The Arctic University
of Norway, Tromsø 9037, Norway
| | - Christian Tantardini
- Hylleraas
Centre for Quantum Molecular Sciences, UiT
The Arctic University of Norway, Tromsø 9037, Norway
- Department
of Chemistry, UiT The Arctic University
of Norway, Tromsø 9037, Norway
- Department
of Materials Science and NanoEngineering, Rice University, Houston, Texas 77005, United States
| | - Quentin Pitteloud
- Hylleraas
Centre for Quantum Molecular Sciences, UiT
The Arctic University of Norway, Tromsø 9037, Norway
- Department
of Chemistry, UiT The Arctic University
of Norway, Tromsø 9037, Norway
| | - Peter Wind
- Hylleraas
Centre for Quantum Molecular Sciences, UiT
The Arctic University of Norway, Tromsø 9037, Norway
- Department
of Chemistry, UiT The Arctic University
of Norway, Tromsø 9037, Norway
| | - Jānis Užulis
- Department
of Physics, University of Latvia, Jelgavas iela 3, Riga, Latvia 1004, Latvia
| | - Andris Gulans
- Department
of Physics, University of Latvia, Jelgavas iela 3, Riga, Latvia 1004, Latvia
| | | | - Luca Frediani
- Hylleraas
Centre for Quantum Molecular Sciences, UiT
The Arctic University of Norway, Tromsø 9037, Norway
- Department
of Chemistry, UiT The Arctic University
of Norway, Tromsø 9037, Norway
| |
Collapse
|
5
|
Rueda Espinosa KJ, Kananenka AA, Rusakov AA. Novel Computational Chemistry Infrastructure for Simulating Astatide in Water: From Basis Sets to Force Fields Using Particle Swarm Optimization. J Chem Theory Comput 2023; 19:7998-8012. [PMID: 38014419 DOI: 10.1021/acs.jctc.3c00826] [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/2023]
Abstract
Using the example of astatine, the heaviest naturally occurring halogen whose isotope At-211 has promising medical applications, we propose a new infrastructure for large-scale computational models of heavy elements with strong relativistic effects. In particular, we focus on developing an accurate force field for At- in water based on reliable relativistic density functional theory (DFT) calculations. To ensure the reliability of such calculations, we design novel basis sets for relativistic DFT, via the particle swarm optimization algorithm to optimize the coefficients of the new basis sets and the polarization-consistent basis set idea's extension to heavy elements to eliminate the basis set error from DFT calculations. The resulting basis sets enable the well-grounded evaluation of relativistic DFT against "gold-standard" CCSD(T) results. Accounting for strong relativistic effects, including spin-orbit interaction, via our redesigned infrastructure, we elucidate a noticeable dissimilarity between At- and I- in halide-water force field parameters, radial distribution functions, diffusion coefficients, and hydration energies. This work establishes the framework for the systematic development of polarization-consistent basis sets for relativistic DFT and accurate force fields for molecular dynamics simulations to be used in large-scale models of complex molecular systems with elements from the bottom of the periodic table, including actinides and even superheavy elements.
Collapse
Affiliation(s)
- Kennet J Rueda Espinosa
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Alexei A Kananenka
- Department of Physics and Astronomy, University of Delaware, Newark, Delaware 19716, United States
| | - Alexander A Rusakov
- Department of Chemistry, Oakland University, Rochester, Michigan 48309, United States
| |
Collapse
|
6
|
Mátyus E, Ferenc D, Jeszenszki P, Margócsy Á. The Bethe-Salpeter QED Wave Equation for Bound-State Computations of Atoms and Molecules. ACS PHYSICAL CHEMISTRY AU 2023; 3:222-240. [PMID: 37249939 PMCID: PMC10214514 DOI: 10.1021/acsphyschemau.2c00062] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 11/04/2022] [Revised: 01/07/2023] [Accepted: 01/09/2023] [Indexed: 05/31/2023]
Abstract
Interactions in atomic and molecular systems are dominated by electromagnetic forces and the theoretical framework must be in the quantum regime. The physical theory for the combination of quantum mechanics and electromagnetism, quantum electrodynamics has been "established" by the mid-twentieth century, primarily as a scattering theory. To describe atoms and molecules, it is important to consider bound states. In the nonrelativistic quantum mechanics framework, bound states can be efficiently computed using robust and general methodologies with systematic approximations developed for solving wave equations. With the sight of the development of a computational quantum electrodynamics framework for atomic and molecular matter, the field theoretic Bethe-Salpeter wave equation expressed in space-time coordinates, its exact equal-time variant, and emergence of a relativistic wave equation, is reviewed. A computational framework, with initial applications and future challenges in relation with precision spectroscopy, is also highlighted.
Collapse
|
7
|
Franzke YJ. Reducing Exact Two-Component Theory for NMR Couplings to a One-Component Approach: Efficiency and Accuracy. J Chem Theory Comput 2023; 19:2010-2028. [PMID: 36939092 DOI: 10.1021/acs.jctc.2c01248] [Citation(s) in RCA: 9] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/21/2023]
Abstract
The self-consistent and complex spin-orbit exact two-component (X2C) formalism for NMR spin-spin coupling constants [ J. Chem. Theory Comput. 17, 2021, 3874-3994] is reduced to a scalar one-component ansatz. This way, the first-order response term can be partitioned into the Fermi-contact (FC) and spin-dipole (SD) interactions as well as the paramagnetic spin-orbit (PSO) contribution. The FC+SD terms are real and symmetric, while the PSO term is purely imaginary and antisymmetric. The relativistic one-component approach is combined with a modern density functional treatment up to local hybrid functionals including the response of the current density. Computational demands are reduced by factors of 8-24 as shown for a large tin compound consisting of 137 atoms. Limitations of the current ansatz are critically assessed for Sn, Pb, Pd, and Pt compounds, i.e. the one-component treatment is not sufficient for tin compounds featuring a few heavy halogen atoms.
Collapse
Affiliation(s)
- Yannick J Franzke
- Fachbereich Chemie, Philipps-Universität Marburg, Hans-Meerwein-Str. 4, 35032 Marburg, Germany
| |
Collapse
|
8
|
Wang X, Sharma S. Relativistic Semistochastic Heat-Bath Configuration Interaction. J Chem Theory Comput 2023; 19:848-855. [PMID: 36700783 DOI: 10.1021/acs.jctc.2c01025] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
In this work we present the extension of semistochastic heat-bath configuration interaction (SHCI) to work with any two-component and four-component Hamiltonian. The vertical detachment energy (VDE) of AuH2- and zero-field splitting (ZFS) of NpO22+ is calculated by correlating more than 100 spinors in both cases. This work demonstrates the capability of SHCI to treat problems where both relativistic effect and electron correlation are important.
Collapse
Affiliation(s)
- Xubo Wang
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado80309, United States
| | - Sandeep Sharma
- Department of Chemistry, University of Colorado Boulder, Boulder, Colorado80309, United States
| |
Collapse
|
9
|
De Santis M, Sorbelli D, Vallet V, Gomes AS, Storchi L, Belpassi L. Frozen-Density Embedding for Including Environmental Effects in the Dirac-Kohn-Sham Theory: An Implementation Based on Density Fitting and Prototyping Techniques. J Chem Theory Comput 2022; 18:5992-6009. [PMID: 36172757 PMCID: PMC9558305 DOI: 10.1021/acs.jctc.2c00499] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2022] [Indexed: 11/28/2022]
Abstract
Frozen density embedding (FDE) represents an embedding scheme in which environmental effects are included from first-principles calculations by considering the surrounding system explicitly by means of its electron density. In the present paper, we extend the full four-component relativistic Dirac-Kohn-Sham (DKS) method, as implemented in the BERTHA code, to include environmental and confinement effects with the FDE scheme (DKS-in-DFT FDE). The implementation, based on the auxiliary density fitting techniques, has been enormously facilitated by BERTHA's python API (PyBERTHA), which facilitates the interoperability with other FDE implementations available through the PyADF framework. The accuracy and numerical stability of this new implementation, also using different auxiliary fitting basis sets, has been demonstrated on the simple NH3-H2O system, in comparison with a reference nonrelativistic implementation. The computational performance has been evaluated on a series of gold clusters (Aun, with n = 2, 4, 8) embedded into an increasing number of water molecules (5, 10, 20, 40, and 80 water molecules). We found that the procedure scales approximately linearly both with the size of the frozen surrounding environment (consistent with the underpinnings of the FDE approach) and with the size of the active system (in line with the use of density fitting). Finally, we applied the code to a series of heavy (Rn) and super-heavy elements (Cn, Fl, Og) embedded in a C60 cage to explore the confinement effect induced by C60 on their electronic structure. We compare the results from our simulations, with respect to more-approximate models employed in the atomic physics literature. Our results indicate that the specific interactions described by FDE are able to improve upon the cruder approximations currently employed, and, thus, they provide a basis from which to generate more-realistic radial potentials for confined atoms.
Collapse
Affiliation(s)
- Matteo De Santis
- Univ.
Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Diego Sorbelli
- Dipartimento
di Chimica, Biologia e Biotecnologie, Università
degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Istituto
di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle
Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Valérie Vallet
- Univ.
Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | | | - Loriano Storchi
- Istituto
di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle
Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Dipartimento
di Farmacia, Università degli Studi
‘G. D’Annunzio’, Via dei Vestini 31, 66100 Chieti, Italy
| | - Leonardo Belpassi
- Istituto
di Scienze e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle
Ricerche c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| |
Collapse
|
10
|
Knecht S, Repisky M, Jensen HJA, Saue T. Exact two-component Hamiltonians for relativistic quantum chemistry: Two-electron picture-change corrections made simple. J Chem Phys 2022; 157:114106. [PMID: 36137811 DOI: 10.1063/5.0095112] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Based on self-consistent field (SCF) atomic mean-field (amf) quantities, we present two simple yet computationally efficient and numerically accurate matrix-algebraic approaches to correct both scalar-relativistic and spin-orbit two-electron picture-change effects (PCEs) arising within an exact two-component (X2C) Hamiltonian framework. Both approaches, dubbed amfX2C and e(xtended)amfX2C, allow us to uniquely tailor PCE corrections to mean-field models, viz. Hartree-Fock or Kohn-Sham DFT, in the latter case also avoiding the need for a point-wise calculation of exchange-correlation PCE corrections. We assess the numerical performance of these PCE correction models on spinor energies of group 18 (closed-shell) and group 16 (open-shell) diatomic molecules, achieving a consistent ≈10-5 Hartree accuracy compared to reference four-component data. Additional tests include SCF calculations of molecular properties such as absolute contact density and contact density shifts in copernicium fluoride compounds (CnFn, n = 2,4,6), as well as equation-of-motion coupled-cluster calculations of x-ray core-ionization energies of 5d- and 6d-containing molecules, where we observe an excellent agreement with reference data. To conclude, we are confident that our (e)amfX2C PCE correction models constitute a fundamental milestone toward a universal and reliable relativistic two-component quantum-chemical approach, maintaining the accuracy of the parent four-component one at a fraction of its computational cost.
Collapse
Affiliation(s)
- Stefan Knecht
- Algorithmiq Ltd, Kanavakatu 3C, FI-00160 Helsinki, Finland
| | - Michal Repisky
- Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, UiT-The Arctic University of Norway, N-9037 Tromsø, Norway
| | - Hans Jørgen Aagaard Jensen
- Department of Physics, Chemistry and Pharmacy, University of Southern Denmark, Campusvej 55, DK-5230 Odense M, Denmark
| | - Trond Saue
- Laboratoire de Chimie et Physique Quantiques (CNRS UMR 5626), Université Toulouse III - Paul Sabatier, 118 Route de Narbonne, F-31062 Toulouse Cedex, France
| |
Collapse
|
11
|
Yuan X, Visscher L, Gomes ASP. Assessing MP2 frozen natural orbitals in relativistic correlated electronic structure calculations. J Chem Phys 2022; 156:224108. [PMID: 35705406 DOI: 10.1063/5.0087243] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
The high computational scaling with the basis set size and the number of correlated electrons is a bottleneck limiting applications of coupled cluster algorithms, in particular for calculations based on two- or four-component relativistic Hamiltonians, which often employ uncontracted basis sets. This problem may be alleviated by replacing canonical Hartree-Fock virtual orbitals by natural orbitals (NOs). In this paper, we describe the implementation of a module for generating NOs for correlated wavefunctions and, in particular, second order Møller-Plesset perturbation frozen natural orbitals (MP2FNOs) as a component of our novel implementation of relativistic coupled cluster theory for massively parallel architectures [Pototschnig et al. J. Chem. Theory Comput. 17, 5509, (2021)]. Our implementation can manipulate complex or quaternion density matrices, thus allowing for the generation of both Kramers-restricted and Kramers-unrestricted MP2FNOs. Furthermore, NOs are re-expressed in the parent atomic orbital (AO) basis, allowing for generating coupled cluster singles and doubles NOs in the AO basis for further analysis. By investigating the truncation errors of MP2FNOs for both the correlation energy and molecular properties-electric field gradients at the nuclei, electric dipole and quadrupole moments for hydrogen halides HX (X = F-Ts), and parity-violating energy differences for H2Z2 (Z = O-Se)-we find MP2FNOs accelerate the convergence of the correlation energy in a roughly uniform manner across the Periodic Table. It is possible to obtain reliable estimates for both energies and the molecular properties considered with virtual molecular orbital spaces truncated to about half the size of the full spaces.
Collapse
Affiliation(s)
- Xiang Yuan
- Univ. Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, F-59000 Lille, France
| | - Lucas Visscher
- Department of Chemistry and Pharmaceutical Sciences, Faculty of Science, Vrije Universiteit Amsterdam, de Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | | |
Collapse
|
12
|
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
|
13
|
Nakai H. Development of Linear-Scaling Relativistic Quantum Chemistry Covering the Periodic Table. BULLETIN OF THE CHEMICAL SOCIETY OF JAPAN 2021. [DOI: 10.1246/bcsj.20210091] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 11/12/2022]
Affiliation(s)
- Hiromi Nakai
- Department of Chemistry and Biochemistry, School of Advanced Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Waseda Research Institute for Science and Engineering (WISE), Waseda University, 3-4-1 Okubo, Shinjuku-ku, Tokyo 169-8555, Japan
- Elements Strategy Initiative for Catalysts and Batteries (ESICB), Kyoto University, Katsura, Kyoto 615-8520, Japan
| |
Collapse
|
14
|
Jeszenszki P, Ferenc D, Mátyus E. All-order explicitly correlated relativistic computations for atoms and molecules. J Chem Phys 2021; 154:224110. [PMID: 34241227 DOI: 10.1063/5.0051237] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
A variational solution procedure is reported for the many-particle no-pair Dirac-Coulomb and Dirac-Coulomb-Breit Hamiltonians aiming at a parts-per-billion (ppb) convergence of the atomic and molecular energies, described within the fixed nuclei approximation. The procedure is tested for nuclear charge numbers from Z = 1 (hydrogen) to 28 (iron). Already for the lowest Z values, a significant difference is observed from leading-order Foldy-Woythusen perturbation theory, but the observed deviations are smaller than the estimated self-energy and vacuum polarization corrections.
Collapse
Affiliation(s)
- Péter Jeszenszki
- Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, Budapest H-1117, Hungary
| | - Dávid Ferenc
- Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, Budapest H-1117, Hungary
| | - Edit Mátyus
- Institute of Chemistry, ELTE, Eötvös Loránd University, Pázmány Péter sétány 1/A, Budapest H-1117, Hungary
| |
Collapse
|
15
|
Halbert L, Vidal ML, Shee A, Coriani S, Severo Pereira Gomes A. Relativistic EOM-CCSD for Core-Excited and Core-Ionized State Energies Based on the Four-Component Dirac-Coulomb(-Gaunt) Hamiltonian. J Chem Theory Comput 2021; 17:3583-3598. [PMID: 33944570 DOI: 10.1021/acs.jctc.0c01203] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
We report an implementation of the core-valence separation approach to the four-component relativistic Hamiltonian-based equation-of-motion coupled-cluster with singles and doubles theory (CVS-EOM-CCSD) for the calculation of relativistic core-ionization potentials and core-excitation energies. With this implementation, which is capable of exploiting double group symmetry, we investigate the effects of the different CVS-EOM-CCSD variants and the use of different Hamiltonians based on the exact two-component (X2C) framework on the energies of different core-ionized and -excited states in halogen- (CH3I, HX, and X-, X = Cl-At) and xenon-containing (Xe, XeF2) species. Our results show that the X2C molecular mean-field approach [Sikkema, J.; J. Chem. Phys. 2009, 131, 124116], based on four-component Dirac-Coulomb mean-field calculations (2DCM), is capable of providing core excitations and ionization energies that are nearly indistinguishable from the reference four-component energies for up to and including fifth-row elements. We observe that two-electron integrals over the small-component basis sets lead to non-negligible contributions to core binding energies for the K and L edges for atoms such as iodine or astatine and that the approach based on Dirac-Coulomb-Gaunt mean-field calculations (2DCGM) are significantly more accurate than X2C calculations for which screened two-electron spin-orbit interactions are included via atomic mean-field integrals.
Collapse
Affiliation(s)
- Loïc Halbert
- CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, Université de Lille, F-59000 Lille, France
| | - Marta L Vidal
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - Avijit Shee
- Department of Chemistry, University of Michigan, Ann Arbor, Michigan 48109, United States
| | - Sonia Coriani
- DTU Chemistry-Department of Chemistry, Technical University of Denmark, DK-2800 Kongens Lyngby, Denmark
| | - André Severo Pereira Gomes
- CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, Université de Lille, F-59000 Lille, France
| |
Collapse
|
16
|
Sorbelli D, De Santis M, Belanzoni P, Belpassi L. Spectroscopic/Bond Property Relationship in Group 11 Dihydrides via Relativistic Four-Component Methods. J Phys Chem A 2020; 124:10565-10579. [PMID: 33327724 PMCID: PMC8016197 DOI: 10.1021/acs.jpca.0c09043] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2020] [Revised: 11/13/2020] [Indexed: 11/30/2022]
Abstract
Group 11 dihydrides MH2- (M = Cu, Ag, Au, Rg) have been much less studied than the corresponding MH compounds, despite having potentially several interesting applications in chemical research. In this work, their main spectroscopic constants (bond lengths, dissociation energies, and force constants) have been evaluated by means of highly accurate relativistic four-component coupled cluster (4c-CCSD(T)) calculations in combination with large basis sets. Periodic trends have been quantitatively explained by the charge-displacement/natural orbitals for chemical valence (CD-NOCV) analysis based on the four-component relativistic Dirac-Kohn-Sham method, which allows a consistent picture of the nature of the M-H bond to be obtained on going down the periodic table in terms of Dewar-Chatt-Duncanson bonding components. A strong ligand-to-metal donation drives the M-H bond and it is responsible for the heterolytic (HM···H-) dissociation energies to increase monotonically from Cu to Rg, with RgH2- showing the strongest and most covalent M-H bond. The "V"-shaped trend observed for the bond lengths, dissociation energies, and stretching frequencies can be explained in terms of relativistic effects and, in particular, of the relativistically enhanced sd hybridization occurring at the metal, which affects the metal-ligand distances in heavy transition-metal complexes. The sd hybridization is very small for Cu and Ag, whereas it becomes increasingly important for Au and Rg, being responsible for the increasing covalent character of the bond, the sizable contraction of the Au-H and Rg-H bonds, and the observed trend. This work rationalizes the spectroscopic/bond property relationship in group 11 dihydrides within highly accurate relativistic quantum chemistry methods, paving the way for their applications in chemical bond investigations involving heavy and superheavy elements.
Collapse
Affiliation(s)
- Diego Sorbelli
- Department
of Chemistry, Biology and Biotechnology, University of Perugia, via Elce di Sotto 8, 06123 Perugia, Italy
| | - Matteo De Santis
- Department
of Chemistry, Biology and Biotechnology, University of Perugia, via Elce di Sotto 8, 06123 Perugia, Italy
- CNR
Institute of Chemical Science and Technologies “Giulio Natta”
(CNR-SCITEC), c/o Department of Chemistry, Biology and Biotechnology, University of Perugia, 06123 Perugia, Italy
| | - Paola Belanzoni
- Department
of Chemistry, Biology and Biotechnology, University of Perugia, via Elce di Sotto 8, 06123 Perugia, Italy
- CNR
Institute of Chemical Science and Technologies “Giulio Natta”
(CNR-SCITEC), c/o Department of Chemistry, Biology and Biotechnology, University of Perugia, 06123 Perugia, Italy
- Consortium
for Computational Molecular and Materials Sciences (CMS)2, via Elce di Sotto 8, 06123 Perugia, Italy
| | - Leonardo Belpassi
- CNR
Institute of Chemical Science and Technologies “Giulio Natta”
(CNR-SCITEC), c/o Department of Chemistry, Biology and Biotechnology, University of Perugia, 06123 Perugia, Italy
- Consortium
for Computational Molecular and Materials Sciences (CMS)2, via Elce di Sotto 8, 06123 Perugia, Italy
| |
Collapse
|
17
|
De Santis M, Belpassi L, Jacob CR, Severo Pereira Gomes A, Tarantelli F, Visscher L, Storchi L. Environmental Effects with Frozen-Density Embedding in Real-Time Time-Dependent Density Functional Theory Using Localized Basis Functions. J Chem Theory Comput 2020; 16:5695-5711. [PMID: 32786918 PMCID: PMC8009524 DOI: 10.1021/acs.jctc.0c00603] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2020] [Indexed: 12/14/2022]
Abstract
Frozen-density embedding (FDE) represents a versatile embedding scheme to describe the environmental effect on electron dynamics in molecular systems. The extension of the general theory of FDE to the real-time time-dependent Kohn-Sham method has previously been presented and implemented in plane waves and periodic boundary conditions [Pavanello, M.; J. Chem. Phys. 2015, 142, 154116]. In the current paper, we extend our recent formulation of the real-time time-dependent Kohn-Sham method based on localized basis set functions and developed within the Psi4NumPy framework to the FDE scheme. The latter has been implemented in its "uncoupled" flavor (in which the time evolution is only carried out for the active subsystem, while the environment subsystems remain at their ground state), using and adapting the FDE implementation already available in the PyEmbed module of the scripting framework PyADF. The implementation was facilitated by the fact that both Psi4NumPy and PyADF, being native Python API, provided an ideal framework of development using the Python advantages in terms of code readability and reusability. We employed this new implementation to investigate the stability of the time-propagation procedure, which is based on an efficient predictor/corrector second-order midpoint Magnus propagator employing an exact diagonalization, in combination with the FDE scheme. We demonstrate that the inclusion of the FDE potential does not introduce any numerical instability in time propagation of the density matrix of the active subsystem, and in the limit of the weak external field, the numerical results for low-lying transition energies are consistent with those obtained using the reference FDE calculations based on the linear-response TDDFT. The method is found to give stable numerical results also in the presence of a strong external field inducing nonlinear effects. Preliminary results are reported for high harmonic generation (HHG) of a water molecule embedded in a small water cluster. The effect of the embedding potential is evident in the HHG spectrum reducing the number of the well-resolved high harmonics at high energy with respect to the free water. This is consistent with a shift toward lower ionization energy passing from an isolated water molecule to a small water cluster. The computational burden for the propagation step increases approximately linearly with the size of the surrounding frozen environment. Furthermore, we have also shown that the updating frequency of the embedding potential may be significantly reduced, much less than one per time step, without jeopardizing the accuracy of the transition energies.
Collapse
Affiliation(s)
- Matteo De Santis
- Dipartimento di
Chimica, Biologia e Biotecnologie, Università
degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Leonardo Belpassi
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Christoph R. Jacob
- Institute of Physical and Theoretical Chemistry, Technische Universität Braunschweig, Gaußstr. 17, 38106 Braunschweig, Germany
| | | | - Francesco Tarantelli
- Dipartimento di
Chimica, Biologia e Biotecnologie, Università
degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
| | - Lucas Visscher
- Theoretical Chemistry, Faculty of Science, Vrije Universiteit Amsterdam, De Boelelaan 1083, 1081 HV Amsterdam, The Netherlands
| | - Loriano Storchi
- Istituto di Scienze
e Tecnologie Chimiche (SCITEC), Consiglio Nazionale delle Ricerche
c/o Dipartimento di Chimica, Biologia e Biotecnologie, Università degli Studi di Perugia, Via Elce di Sotto 8, 06123 Perugia, Italy
- Dipartimento di Farmacia, Università
degli Studi ‘G. D’Annunzio’, Via dei Vestini 31, 66100 Chieti, Italy
| |
Collapse
|
18
|
Sherrill CD, Manolopoulos DE, Martínez TJ, Michaelides A. Electronic structure software. J Chem Phys 2020; 153:070401. [DOI: 10.1063/5.0023185] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/19/2022] Open
Affiliation(s)
- 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-0400, USA
| | - David E. Manolopoulos
- Physical and Theoretical Chemistry Laboratory, Department of Chemistry, Oxford University, South Parks Road, Oxford OX1
3QZ, United Kingdom
| | - Todd J. Martínez
- Department of Chemistry and the PULSE Institute, Stanford University, Stanford, California 94305,
USA
| | - Angelos Michaelides
- Thomas Young Centre, London Centre for Nanotechnology and Department of Physics and Astronomy, University College London, Gower Street, London WC1E 6BT, United Kingdom
| |
Collapse
|