1
|
Desmarais JK, Maul J, Civalleri B, Erba A, Vignale G, Pittalis S. Spin Currents via the Gauge Principle for Meta-Generalized Gradient Exchange-Correlation Functionals. PHYSICAL REVIEW LETTERS 2024; 132:256401. [PMID: 38996240 DOI: 10.1103/physrevlett.132.256401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/14/2023] [Accepted: 05/15/2024] [Indexed: 07/14/2024]
Abstract
The prominence of density functional theory in the field of electronic structure computation stems from its ability to usefully balance accuracy and computational effort. At the base of this ability is a functional of the electron density: the exchange-correlation energy. This functional satisfies known exact conditions that guide the derivation of approximations. The strongly constrained and appropriately normed (SCAN) approximation stands out as a successful, modern, example. In this Letter, we demonstrate how the SU(2) gauge invariance of the exchange-correlation functional in spin current density functional theory allows us to add an explicit dependence on spin currents in the SCAN functional (here called JSCAN)-and similar meta-generalized-gradient functional approximations-solely invoking first principles. In passing, a spin-current dependent generalization of the electron localization function (here called JELF) is also derived. The extended forms are implemented in a developer's version of the crystal23 program. Applications on molecules and materials confirm the practical relevance of the extensions.
Collapse
Affiliation(s)
| | | | | | | | | | - Stefano Pittalis
- Istituto Nanoscienze, Consiglio Nazionale delle Ricerche, Via Campi 213A, I-41125 Modena, Italy
| |
Collapse
|
2
|
Erba A, Desmarais JK, Casassa S, Civalleri B, Donà L, Bush IJ, Searle B, Maschio L, Edith-Daga L, Cossard A, Ribaldone C, Ascrizzi E, Marana NL, Flament JP, Kirtman B. CRYSTAL23: A Program for Computational Solid State Physics and Chemistry. J Chem Theory Comput 2023; 19:6891-6932. [PMID: 36502394 PMCID: PMC10601489 DOI: 10.1021/acs.jctc.2c00958] [Citation(s) in RCA: 20] [Impact Index Per Article: 20.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2022] [Indexed: 12/14/2022]
Abstract
The Crystal program for quantum-mechanical simulations of materials has been bridging the realm of molecular quantum chemistry to the realm of solid state physics for many years, since its first public version released back in 1988. This peculiarity stems from the use of atom-centered basis functions within a linear combination of atomic orbitals (LCAO) approach and from the corresponding efficiency in the evaluation of the exact Fock exchange series. In particular, this has led to the implementation of a rich variety of hybrid density functional approximations since 1998. Nowadays, it is acknowledged by a broad community of solid state chemists and physicists that the inclusion of a fraction of Fock exchange in the exchange-correlation potential of the density functional theory is key to a better description of many properties of materials (electronic, magnetic, mechanical, spintronic, lattice-dynamical, etc.). Here, the main developments made to the program in the last five years (i.e., since the previous release, Crystal17) are presented and some of their most noteworthy applications reviewed.
Collapse
Affiliation(s)
- Alessandro Erba
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jacques K. Desmarais
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Silvia Casassa
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Bartolomeo Civalleri
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Lorenzo Donà
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Ian J. Bush
- STFC
Rutherford Appleton Laboratory, Chilton Didcot, Oxfordshire OX11 0QX, United Kingdom
| | - Barry Searle
- SFTC
Daresbury Laboratory, Daresbury, Cheshire WA4 4AD, United Kingdom
| | - Lorenzo Maschio
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Loredana Edith-Daga
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Alessandro Cossard
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Chiara Ribaldone
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Eleonora Ascrizzi
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Naiara L. Marana
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université
de Lille, CNRS, UMR 8523 — PhLAM — Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department
of Chemistry and Biochemistry, University
of California, Santa
Barbara, California 93106, United States
| |
Collapse
|
3
|
Desmarais JK, Boccuni A, Flament JP, Kirtman B, Erba A. Perturbation Theory Treatment of Spin-Orbit Coupling. III: Coupled Perturbed Method for Solids. J Chem Theory Comput 2023; 19:1853-1863. [PMID: 36917759 DOI: 10.1021/acs.jctc.3c00088] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/16/2023]
Abstract
A previously proposed noncanonical coupled-perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment for spin-orbit coupling is here generalized to infinite periodic systems. The scalar-relativistic periodic KS-DFT/HF solution, obtained with a relativistic effective core potential, is taken as the zeroth-order approximation. Explicit expressions are given for the total energy through third-order, which satisfy the 2N + 1 rule (i.e., requiring only the first-order perturbed wave function for determining the energy through third-order). Expressions for additional second-order corrections to the perturbed wave function (as well as related one-electron properties) are worked out at the uncoupled-perturbed level of theory. The approach is implemented in the Crystal program and validated with calculations of the total energy, electronic band structure, and density variables of spin-current DFT on the tungsten dichalcogenide hexagonal bilayer series (i.e., WSe2, WTe2, WPo2, WLv2), including 6p and 7p elements as a stress test. The computed properties through second- or third-order match well with those from reference two-component self-consistent field (2c-SCF) calculations. For total energies, E(3) was found to consistently improve the agreement against the 2c-SCF reference values. For electronic band structures, visible differences w.r.t. 2c-SCF remained through second-order in only the single-most difficult case of WLv2. As for density variables of spin-current DFT, the perturbed electron density, being vanishing in first-order, is the most challenging for the perturbation theory approach. The visible differences in the electron densities are, however, largest close to the core region of atoms and smaller in the valence region. Perturbed spin-current densities, on the other hand, are well reproduced in all tested cases.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Alberto Boccuni
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523 ─ PhLAM ─ Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| |
Collapse
|
4
|
Desmarais JK, Erba A, Flament JP, Kirtman B. Perturbation Theory Treatment of Spin-Orbit Coupling, Part I: Double Perturbation Theory Based on a Single-Reference Initial Approximation. J Chem Theory Comput 2021; 17:4697-4711. [PMID: 34288690 DOI: 10.1021/acs.jctc.1c00343] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
We develop a perturbation theory for solving the many-body Dirac equation within a given relativistic effective-core potential approximation. Starting from a scalar-relativistic unrestricted Hartree-Fock (SR UHF) solution, we carry out a double perturbation expansion in terms of spin-orbit coupling (SOC) and the electron fluctuation potential. Computationally convenient energy expressions are derived through fourth order in SOC, second order in the electron fluctuation potential, and a total of third order in the coupling between the two. Illustrative calculations on the halogen series of neutral and singly positive diatomic molecules show that the perturbation expansion is well-converged by taking into account only the leading (nonvanishing) term at each order of the electron fluctuation potential. Our perturbation theory approach provides a computationally attractive alternative to a two-component self-consistent field treatment of SOC. In addition, it includes coupling with the fluctuation potential through third order and can be extended (in principle) to multireference calculations, when necessary for both closed- and open-shell cases, using quasi-degenerate perturbation theory.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy.,IPREM, E2S UPPA, CNRS, Université de Pau et des Pays de l'Adour, 64053 Pau, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, CNRS, Université de Lille, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, Santa Barbara, California 93106, United States
| |
Collapse
|
5
|
Desmarais JK, Erba A, Flament JP, Kirtman B. Perturbation Theory Treatment of Spin-Orbit Coupling II: A Coupled Perturbed Kohn-Sham Method. J Chem Theory Comput 2021; 17:4712-4732. [PMID: 34286577 DOI: 10.1021/acs.jctc.1c00460] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
A noncanonical coupled perturbed Kohn-Sham density functional theory (KS-DFT)/Hartree-Fock (HF) treatment of spin-orbit coupling (SOC) is provided. We take the scalar-relativistic KS-DFT/HF solution, obtained with a relativistic effective core potential, as the zeroth-order approximation. Explicit expressions are given for the total energy through the 4th order, which satisfy the 2n + 1 rule. Second-order expressions are provided for orbital energies and density variables of spin-current DFT. Test calculations are carried out on the halogen homonuclear diatomic and hydride molecules, including 6p and 7p elements, as well as open-shell negative ions. The computed properties through second or third order match well with those from reference two-component self-consistent field calculations for total and orbital energies as well as spin-current densities. In only one case (At2-) did a significant deviation occur for the remaining density variables. Our coupled perturbation theory approach provides an efficient way of adding the effect of SOC to a scalar-relativistic single-reference KS-DFT/HF treatment, in particular because it does not require diagonalization in the two-component spinor basis, leading to saving factors on the number of required floating-point operations that may exceed one order of magnitude.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy.,Université de Pau et des Pays de l'Adour, E2S UPPA, CNRS, IPREM, 64000 Pau, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523-PhLAM-Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Bernard Kirtman
- Department of Chemistry and Biochemistry, University of California, Santa Barbara, California 93106, United States
| |
Collapse
|
6
|
Cossard A, Casassa S, Gatti C, Desmarais JK, Erba A. Topology of the Electron Density and of Its Laplacian from Periodic LCAO Calculations on f-Electron Materials: The Case of Cesium Uranyl Chloride. Molecules 2021; 26:4227. [PMID: 34299502 PMCID: PMC8303866 DOI: 10.3390/molecules26144227] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2021] [Revised: 07/06/2021] [Accepted: 07/08/2021] [Indexed: 11/28/2022] Open
Abstract
The chemistry of f-electrons in lanthanide and actinide materials is yet to be fully rationalized. Quantum-mechanical simulations can provide useful complementary insight to that obtained from experiments. The quantum theory of atoms in molecules and crystals (QTAIMAC), through thorough topological analysis of the electron density (often complemented by that of its Laplacian) constitutes a general and robust theoretical framework to analyze chemical bonding features from a computed wave function. Here, we present the extension of the Topond module (previously limited to work in terms of s-, p- and d-type basis functions only) of the Crystal program to f- and g-type basis functions within the linear combination of atomic orbitals (LCAO) approach. This allows for an effective QTAIMAC analysis of chemical bonding of lanthanide and actinide materials. The new implemented algorithms are applied to the analysis of the spatial distribution of the electron density and its Laplacian of the cesium uranyl chloride, Cs2UO2Cl4, crystal. Discrepancies between the present theoretical description of chemical bonding and that obtained from a previously reconstructed electron density by experimental X-ray diffraction are illustrated and discussed.
Collapse
Affiliation(s)
- Alessandro Cossard
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy; (A.C.); (S.C.)
| | - Silvia Casassa
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy; (A.C.); (S.C.)
| | - Carlo Gatti
- CNR-SCITEC, Istituto di Scienze e Tecnologie Chimiche “Giulio Natta”, Via C. Golgi 19, 20133 Milano, Italy;
| | - Jacques K. Desmarais
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy; (A.C.); (S.C.)
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy; (A.C.); (S.C.)
| |
Collapse
|
7
|
Desmarais JK, Komorovsky S, Flament JP, Erba A. Spin–orbit coupling from a two-component self-consistent approach. II. Non-collinear density functional theories. J Chem Phys 2021; 154:204110. [DOI: 10.1063/5.0051447] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022] Open
Affiliation(s)
- Jacques K. Desmarais
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
- Université de Pau et des Pays de l’Adour, E2S UPPA, CNRS, IPREM, Pau, France
| | - Stanislav Komorovsky
- Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská cesta 9, SK-84536 Bratislava, Slovakia
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523—PhLAM—Physique des Lasers, Atomes et Molécules, 59000 Lille, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| |
Collapse
|
8
|
Desmarais JK, Erba A, Pan Y, Civalleri B, Tse JS. Mechanisms for Pressure-Induced Isostructural Phase Transitions in EuO. PHYSICAL REVIEW LETTERS 2021; 126:196404. [PMID: 34047588 DOI: 10.1103/physrevlett.126.196404] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/05/2021] [Revised: 03/01/2021] [Accepted: 04/22/2021] [Indexed: 06/12/2023]
Abstract
We study pressure-induced isostructural electronic phase transitions in the prototypical mixed valence and strongly correlated material EuO using the global-hybrid density functional theory. The simultaneous presence in the valence of highly localized d- and f-type bands and itinerant s- and p-type states, as well as the half-filled f-type orbital shell with seven unpaired electrons on each Eu atom, have made the description of the electronic features of this system a challenge. The electronic band structure, density of states, and atomic oxidation states of EuO are analyzed in the 0-50 GPa pressure range. An insulator-to-metal transition at about 12 GPa of pressure was identified. The second isostructural transition at approximately 30-35 GPa, previously believed to be driven by an oxidation from Eu(II) to Eu(III), is shown instead to be associated with a change in the occupation of the Eu d orbitals, as can be determined from the analysis of the corresponding atomic orbital populations. The Eu d band is confined by the surrounding oxygens and split by the crystal field, which results in orbitals of e_{g} symmetry (i.e., d_{x^{2}-y^{2}} and d_{2z^{2}-x^{2}-y^{2}}, pointing along the Eu-O direction) being abruptly depopulated at the transition as a means to alleviate electron-electron repulsion in the highly compressed structures.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica e NIS centro interdipartimentale, Università di Torino, Via P. Giuria 7, 10125 Torino, Italy
- Equipe de Chimie Physique, IPREM UMR5254, Université de Pau et des Pays de lAdour, 64053, Pau, CEDEX 9, France
- Department of Geological Sciences, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E2, Canada
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E2, Canada
| | - Alessandro Erba
- Dipartimento di Chimica e NIS centro interdipartimentale, Università di Torino, Via P. Giuria 7, 10125 Torino, Italy
| | - Yuanming Pan
- Department of Geological Sciences, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E2, Canada
| | - Bartolomeo Civalleri
- Dipartimento di Chimica e NIS centro interdipartimentale, Università di Torino, Via P. Giuria 7, 10125 Torino, Italy
| | - John S Tse
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E2, Canada
| |
Collapse
|
9
|
Cossard A, Desmarais JK, Casassa S, Gatti C, Erba A. Charge Density Analysis of Actinide Compounds from the Quantum Theory of Atoms in Molecules and Crystals. J Phys Chem Lett 2021; 12:1862-1868. [PMID: 33577336 PMCID: PMC8028320 DOI: 10.1021/acs.jpclett.1c00100] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Accepted: 02/10/2021] [Indexed: 05/17/2023]
Abstract
The nature of chemical bonding in actinide compounds (molecular complexes and materials) remains elusive in many respects. A thorough analysis of their electron charge distribution can prove decisive in elucidating bonding trends and oxidation states along the series. However, the accurate determination and robust analysis of the charge density of actinide compounds pose several challenges from both experimental and theoretical perspectives. Significant advances have recently been made on the experimental reconstruction and topological analysis of the charge density of actinide materials [Gianopoulos et al. IUCrJ, 2019, 6, 895]. Here, we discuss complementary advances on the theoretical side, which allow for the accurate determination of the charge density of actinide materials from quantum-mechanical simulations in the bulk. In particular, the extension of the Topond software implementing Bader's quantum theory of atoms in molecules and crystals (QTAIMAC) to f- and g-type basis functions is introduced, which allows for an effective study of lanthanides and actinides in the bulk and in vacuo, on the same grounds. Chemical bonding of the tetraphenyl phosphate uranium hexafluoride cocrystal [PPh4+][UF6-] is investigated, whose experimental charge density is available for comparison. Crystal packing effects on the charge density and chemical bonding are quantified and discussed. The methodology presented here allows reproducing all subtle features of the topology of the Laplacian of the experimental charge density. Such a remarkable qualitative and quantitative agreement represents a strong mutual validation of both approaches-experimental and computational-for charge density analysis of actinide compounds.
Collapse
Affiliation(s)
- Alessandro Cossard
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Jacques K. Desmarais
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Silvia Casassa
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| | - Carlo Gatti
- CNR-SCITEC,
Istituto di Scienze e Tecnologie Chimiche “Giulio Natta”, via C. Golgi 19, 20133 Milano, Italy
| | - Alessandro Erba
- Dipartimento
di Chimica, Università di Torino, via Giuria 5, 10125 Torino, Italy
| |
Collapse
|
10
|
Zöllner MS, Saghatchi A, Mujica V, Herrmann C. Influence of Electronic Structure Modeling and Junction Structure on First-Principles Chiral Induced Spin Selectivity. J Chem Theory Comput 2020; 16:7357-7371. [PMID: 33167619 DOI: 10.1021/acs.jctc.0c00621] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
We have carried out a comprehensive study of the influence of electronic structure modeling and junction structure description on the first-principles calculation of the spin polarization in molecular junctions caused by the chiral induced spin selectivity (CISS) effect. We explore the limits and the sensitivity to modeling decisions of a Landauer/Green's function/two-component density functional theory approach to CISS. We find that although the CISS effect is entirely attributed in the literature to molecular spin filtering, spin-orbit coupling being partially inherited from the metal electrodes plays an important role in our calculations on ideal carbon helices, even though this effect cannot explain the experimental conductance results. Its magnitude depends considerably on the shape, size, and material of the metal clusters modeling the electrodes. Also, a pronounced dependence on the specific description of exchange interaction and spin-orbit coupling is manifest in our approach. This is important because the interplay between exchange effects and spin-orbit coupling may play an important role in the description of the junction magnetic response. Our calculations are relevant for the whole field of spin-polarized electron transport and electron transfer, because there is still an open discussion in the literature about the detailed underlying mechanism and the magnitude of physical parameters that need to be included to achieve a consistent description of the CISS effect: seemingly good quantitative agreement between simulation and the experiment can be caused by error compensation, because spin polarization as contained in a Landauer/Green's function/two-component density functional theory approach depends strongly on computational and structural parameters.
Collapse
Affiliation(s)
| | - Aida Saghatchi
- Department of Chemistry, University of Hamburg, 20146 Hamburg, Germany
| | - Vladimiro Mujica
- School of Molecular Sciences, Arizona State University, Tempe, Arizona 85287-1604, United States.,Kimika Fakultatea, Euskal Herriko Unibertsitatea and Donostia International Physics Center (DIPC), Donostia, Euskadi P.K. 1072, 20080, Spain
| | - Carmen Herrmann
- Department of Chemistry, University of Hamburg, 20146 Hamburg, Germany
| |
Collapse
|
11
|
Desmarais JK, Flament JP, Erba A. Spin-orbit coupling from a two-component self-consistent approach. II. Non-collinear density functional theories. J Chem Phys 2019; 151:074108. [DOI: 10.1063/1.5114902] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Affiliation(s)
- Jacques K. Desmarais
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
- Department of Geological Sciences, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E2,
Canada
- Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, Saskatchewan S7N 5E2,
Canada
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molècules, 59000 Lille,
France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
| |
Collapse
|
12
|
Desmarais JK, Flament JP, Erba A. Spin-orbit coupling from a two-component self-consistent approach. I. Generalized Hartree-Fock theory. J Chem Phys 2019; 151:074107. [PMID: 31438694 DOI: 10.1063/1.5114901] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
Abstract
Formal and computational aspects are discussed for a self-consistent treatment of spin-orbit coupling within the two-component generalization of the Hartree-Fock theory. A molecular implementation into the CRYSTAL program is illustrated, where the standard one-component code (typical of Hartree-Fock and Kohn-Sham spin-unrestricted methodologies) is extended to work in terms of two-component spinors. When passing from a one- to a two-component description, the Fock and density matrices become complex. Furthermore, apart from the αα and ββ diagonal spin blocks, one has also to deal with the αβ and βα off-diagonal spin blocks. These latter blocks require special care as, for open-shell electronic configurations, certain constraints of the one-component code have to be relaxed. This formalism intrinsically allows us to treat local magnetic torque as well as noncollinear magnetization and orbital current-density. An original scheme to impose a specified noncollinear magnetization on each atomic center as a starting guess to the self-consistent procedure is presented. This approach turns out to be essential to surpass local minima in the rugged energy landscape and allows possible convergence to the ground-state solution in all of the discussed test cases.
Collapse
Affiliation(s)
- Jacques K Desmarais
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
| | - Jean-Pierre Flament
- Université de Lille, CNRS, UMR 8523 - PhLAM - Physique des Lasers Atomes et Molécules, 59000 Lille, France
| | - Alessandro Erba
- Dipartimento di Chimica, Università di Torino, Via Giuria 5, 10125 Torino, Italy
| |
Collapse
|