Density functional theory in magnetic fields
|
|
- Juniper Hampton
- 6 years ago
- Views:
Transcription
1 Density functional theory in magnetic fields U. E. Ekström, T. Helgaker, S. Kvaal, E. Sagvolden, E. Tellgren, A. M. Teale Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway School of Chemistry, University of Nottingham, Nottingham, UK Electronic Structure Theory for Strongly Correlated Systems celebrating Per Åke Malmqvist s 6th birthday and his career in quantum chemistry Splendid Hotel La Torre, Mondello, Palermo, Italy, May 3 June 1, 212 T. Helgaker (CTCC, University of Oslo) Calculation of the universal density functional Malmqvist Celebration Palermo / 23
2 Overview 1 Density functional theory (DFT) 2 The universal density functional and its calculation by Lieb maximization 3 Regularization of DFT 4 The adiabatic connection with examples 5 Molecules in magnetic fields 6 Magnetic-field DFT (BDFT) and current DFT (CDFT) Helgaker et al. (CTCC, University of Oslo) Overview Malmqvist Celebration Palermo / 23
3 The universal density functional Consider the ground-state energy with interaction strength λ: E λ [v] = inf Ψ T + i v(r i ) + λw Ψ the Rayleigh Ritz variation principle Ψ We may summarize density-functional theory in two variation principles: ( E λ [v] = inf Fλ [ρ] + (v ρ) ) the Hohenberg Kohn variation principle (1964) ρ ( F λ [ρ] = sup Eλ [v] (v ρ) ) the Lieb variation principle (1983) v these are alternative attempts at sharpening the same inequality into an equality F λ [ρ] E λ [v] (v ρ) E λ [v] F λ [ρ] + (v ρ) Fenchel s inequality The Hohenberg Kohn variation principle underlies all applications of DFT in chemistry its success hinges on the accurate modelling of F λ [ρ] The Lieb variation principle provides a tool for studying F λ [ρ] numerically we have implemented the Lieb variation principle at various levels ab initio theory accurate numerical studies of F λ [ρ] are a guide to insight our motivation is to benchmark approximate F λ [ρ] and to help develop new ones Previous work: studies by Colonna and Savin (1999) for few-electron atoms work by Wu and Yang (23) on Lieb maximizations Helgaker et al. (CTCC, University of Oslo) Lieb s universal density functional Malmqvist Celebration Palermo / 23
4 Conjugate functionals E[v] F [ρ] The ground-state energy E[v] is concave in v by the Rayleigh Ritz variation principle it can therefore be exactly represented by its convex conjugate F [ρ]: E[v] F [ρ] FΡ sup v EvvΡ FΡ FΡ Ρ min v max Ev Ev Ev inf Ρ FΡvΡ vρ vρ Both variation principles represent a Legendre Fenchel transformation the essential point is the concavity of E[v] rather than the details of the Schrödinger equation roughly speaking, the derivatives of E[v] are F [ρ] are each other s inverse functions the density functional F [ρ] is more complicated that E[v] Helgaker et al. (CTCC, University of Oslo) Lieb s universal density functional Malmqvist Celebration Palermo / 23
5 Lieb maximization For a given density ρ(r) and level of theory E λ [v], we wish to calculate F λ [ρ] = max v ( Eλ [v] (v ρ) ) The potential is parameterized as suggested by Wu and Yang 23: v c(r) = v ext(r) + (1 λ)v ref (r) + ct gt(r) t the physical, external potential v ext(r) the Fermi Amaldi reference potential to ensure correct asymptotic behaviour an expansion in Gaussians g t(r) with coefficients c t The maximization may be carried out by Newton-type methods: because of concavity, all maxima are global maxima 5 1 iterations with the exact Hessian to gradient norm iterations with an approximate noninteracting Hessian Hessian regularized by a singular-value decomposition with cutoff 1 6 possible (unphysical) oscillations in the potential do not affect F λ [ρ] Convergence is mostly unproblematic but difficulties may arise for λ from a (near) piecewise linearity in c t from a poor description of E λ [v c] for large c t All code is implemented in DALTON for HF, MP2, CCD, CCSD, and CCSD(T) Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
6 The bundle method for near piecewise linear functions The goal function contains alternating regions of small and large curvature The cutting-plane method minimizes linearized functions (supporting hyperplanes) The bundle method combines Newton s method (Hessian information) with cutting planes Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
7 Moreau Yosida regularization The exact density functional is complicated: F λ [ρ] is convex, unbounded, everywhere discontinuous and not everywhere (sub)differentiable not all densities are v-representable or even N-representable ( F λ [ρ] = sup Eλ [v] (v ρ) ) ( = sup ψv T + λw ψ v + (v ρ v ρ) ) v v in practice, it may be impossible to satisfy ρ v = ρ in a small orbital basis maximization may then give a too large potential v and a too high F λ [ρ] This may be avoided by adding a penalty function to the ground-state energy: E µ λ [v] = E λ[v] µ 2 v 2 penalty for large potentials The Lieb variation principle now gives a bounded regularized universal density functional: F µ ( λ [ρ] = sup E µ λ [v] (v ρ)) v where the Moreau envelope of the density functional is given by ( ) F µ λ [ρ] = inf F λ [n] + 1 n 2µ ρ n 2 infimal convolution it represents the intrinsic energy of the most favourable nearby density n ρ Moreau Yosida density functional is simple: F µ λ [ρ] with is convex, bounded, everywhere continuous and differentiable for µ > all densities are N-representable and even v-representable Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
8 några minuter till. Moreau envelopes Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
9 Choice of regularization metric The Moreau Yosida regularization norms are related by convex conjugation: v 2 = v M v ρ 2 = ρ M 1 ρ convex conjugation Different M and µ lead to different regularizations: E µ λ [v] = E λ[v] µ 2 v M v extrinsic energy with damped potential ( ) F µ λ [ρ] = inf F λ [ρ + ] + 1 2µ M 1 intrinsic energy of favoured nearby density where contains no electrons and only changes the shape Two choices of metric: M = I (overlap) and M = T (kinetic) The kinetic metric connects two established regularizations the regularization of Heaton Burgess, Bulat and Yang (27) for OEPs: v T v = v(r) 2 dr small for smooth potentials the regularization of Zhao, Morrison and Parr (ZMP) (1994) for KS potentials: T 1 = (r 1 )r 1 12 (r 2) dr 1 dr 2 small for oscillatory densities v becomes smoother as density oscillations are removed Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
10 MY regularization for the neon atom with µ = 1 18 aug-cc-pvtz orbital basis; 35s even-tempered potential basis Μ Μ TSVD cutoff 1^18 1. MY S Metric Μ 1^18 1. MY T Metric Μ 1^ TSVD cutoff 1^ MY S Metric Μ 1^ MY T Metric Μ 1^ Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
11 MY regularization for the neon atom with µ = 1 9 aug-cc-pvtz orbital basis; 35s even-tempered potential basis Μ Μ TSVD cutoff 1^9 1. MY S Metric Μ 1^9 1. MY T Metric Μ 1^ TSVD cutoff 1^9 6.5 MY S Metric Μ 1^9 6.5 MY T Metric Μ 1^ Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
12 MY regularization for the neon atom with µ = 1 5 aug-cc-pvtz orbital basis; 35s even-tempered potential basis Μ Μ TSVD cutoff 1^5 1. MY S Metric Μ 1^5 1. MY T Metric Μ 1^ TSVD cutoff 1^5 6.5 MY S Metric Μ 1^5 6.5 MY T Metric Μ 1^ Helgaker et al. (CTCC, University of Oslo) Lieb maximization Malmqvist Celebration Palermo / 23
13 Kohn Sham theory and the adiabatic connection We introduce Kohn Sham theory by expanding about λ = : F λ [ρ] = F [ρ] + λf [ρ] + E c,λ[ρ] = T s[ρ] + λ(j[ρ] K[ρ]) + E c,λ [ρ] The correlation energy is the only term that depends on λ in a nontrivial manner: E c,λ [ρ] = F λ [ρ] F [ρ] = Ψ λ V ee Ψ λ Ψ V ee Ψ = W λ [ρ] AC integrand The correlation energy obtained by integration over the AC integrand: 1 E c[ρ] = W λ [ρ] dλ CCSD and CCSD(T) AC curves for the water molecule: E c Ρ WΛ a.u T c Ρ CCSD 9.5 CCSDT Λ Helgaker et al. (CTCC, University of Oslo) The adiabatic connection Malmqvist Celebration Palermo / 23
14 H 2 dissociation: from dynamical to static correlation Near equilibrium, correlation is predominantly dynamic nearly linear curve since correlation is dominated by doubles (λ 2 in PT) Towards dissociation, static correlation becomes important increased curvature form higher-order contributions in PT At dissociation, correlation is fully static wave function adjusts immediately for λ (degenerate PT) first-order degenerate PT yields constant curve for λ > bohr bohr.25 1 bohr Helgaker et al. (CTCC, University of Oslo) The adiabatic connection Malmqvist Celebration Palermo / 23
15 A two-parameter model for the H 2 molecule From a two-level CI model, we obtain the ground-state energy E CI (λ) = 1 2 E 1 2 E 2 + 4g 2 λ 2, E = h + gλ Differentiation gives a two-parameter model for the AC integrand good least-square fits (full lines) and fits to initial gradient and end point (dashed lines). Wc,Λ a.u R.7 a.u. R 1.4 a.u. R 3. a.u..2 R 5. a.u..25 R 7. a.u. R 1. a.u Only the initial slope (GL2) and the end point (perhaps at λ = ) are needed Helgaker et al. (CTCC, University of Oslo) The adiabatic connection Malmqvist Celebration Palermo / 23
16 BLYP and FCI correlation curves for H 2 The BLYP functional treats correlation as dynamical at all bond distances it was designed for spin-unrestricted theory but is used here in a spin-restricted manner it hence ignores static correlation BLYP FCI BLYP FCI R 1.4 bohr R 3. bohr BLYP BLYP R 5. bohr FCI R 1. bohr FCI Helgaker et al. (CTCC, University of Oslo) The adiabatic connection Malmqvist Celebration Palermo / 23
17 BLYP and FCI exchange correlation curves for H 2 The improved BLYP performance arises from an overestimation of exchange error cancellation between exchange and correlation reduces total error to about one third R 1.4 bohr R 3. bohr FCI exchange BLYP exchange BLYP correlation FCI correlation R 5. bohr R 1. bohr Helgaker et al. (CTCC, University of Oslo) The adiabatic connection Malmqvist Celebration Palermo / 23
18 Molecules in magnetic fields A magnetic field B modifies the kinetic-energy operator: H(B) = T (B) + W + i v(r i ), T (B) = 1 2 i (σ π i ) 2 where σ are the Pauli spin matrices and π i the kinetic-momentum operator: π i = i i + A(r i ), A = 1 B (r O) 2 We have recently developed the London code for molecular calculations in strong fields complex wave functions and London atomic orbitals Hartree Fock, FCI and Kohn Sham models Our calculations have revealed many interesting features: transition from paramagnetism to diamagnetism for all molecules at a critical B c molecular bonding of triplet H 2 and singlet He 2 in strong fields EêkJmol -1 EêkJmol Σ -55 H2 1Σ u He HF s A HF 1s B FCI Rêpm Σ g FCI Rêpm What can we learn about the exchange correlation functional in magnetic fields? Helgaker et al. (CTCC, University of Oslo) DFT in magnetic fields Malmqvist Celebration Palermo / 23
19 BDFT: magnetic-field DFT The simplest approach is to set up a separate DFT for each B: E [v, B] = inf ρ (F [ρ, B] + (v ρ)) magnetic-field DFT (BDFT) (Grayce & Harris 1994) The density functional now depends on the density and on the field strength F [ρ, B] = min Ψ Tπ(B) + W Ψ = Ts[ρ, B] + J[ρ] + Exc[ρ, B] Ψ ρ the noninteracting magnetic response is exactly taken care of by T s[ρ, B] AC correlation curves of H 2 in a perpendicular magnetic field E xc[ρ, B] differs from E xc[ρ] only at stretched geometries an earlier onset of static correlation in strong magnetic fields.68 W H H R=1.4 au B= B=.25 au.45 W.4.35 H H R=5. au B= B=.25 au dw/db λ λ Helgaker et al. (CTCC, University of Oslo) DFT in magnetic fields Malmqvist Celebration Palermo / 23
20 CDFT: current DFT Alternatively, we can build DFT with (ρ, j) as variables, conjugate to (v, A) The current consists of paramagnetic and diamagnetic parts: j = j p + j d = j p + ρ A, j p = i 1 2 ψ (r) ψ(r) + c.c. Note: different physical systems may share the same ground-state wave function: (ρ, j) (v, A) (ρ, j p) Ψ up to convex combinations and gauge transformations Eschrig (21), Capelle & Vignale (21), Pan & Sahni (21) We may now set up current DFT (CDFT): ( E [v, A] = inf F [ρ, jp] + ( v 1 ρ,j 2 A2 ρ ) + ( A j )) the density functional depends on ρ and j p = j ρa: F [ρ, j p] = min Ψ T + W Ψ Ψ ρ, j p Vignale and Rasolt (1987,1988), Pan and Sahni (21) A Lieb variation principle can be set up to calculate F [ρ, j p] we have a pilot implementation... Helgaker et al. (CTCC, University of Oslo) DFT in magnetic fields Malmqvist Celebration Palermo / 23
21 Kohn Sham CDFT theory A Kohn Sham decomposition of the CDFT functional yields (Vignale and Rasolt): F [ρ, j p] = T s[ρ, j p] + J[ρ] + E xc[ρ, ν], ν(r) = jp(r) ρ(r) all gauge dependence is located in the kinetic energy in the usual manner the exchange and correlation functionals are separately gauge invariant they depend on the vorticity rather than the paramagnetic current vorticity The Vignale Rasolt Geldart (VRG) current density functional (1987,1988) ( ) E xc[ρ, j p] = Exc [ρ] + kf (r s) χ L 24π 2 χ 1 (r s) ν(r) dr, L ( ) 3 1/3 ( ) 9π 1/3 r s =, k F (r s) = rs 1 4πρ(r) 4 several parametrizations of the diamagnetic susceptibility ratio χ L /χ L are available We have completed a KS-CDFT implementation to the london program VRG calculations on shieldings in water (Huz II basis) HF LDA VRG(1) VRG(2) VRG(3) O H Lee, Handy and Colwell, JCP 13, 195 (1995) truncated at r s = 1 Helgaker et al. (CTCC, University of Oslo) DFT in magnetic fields Malmqvist Celebration Palermo / 23
22 6 6 1 Density 4 and vorticity in benzene log 1 () B/ B (in plane) log 1 () 6.B/ B (1 bohr above) log 1 () B/ B (1 bohr above) Helgaker et al. (CTCC, University of Oslo) 4 DFT in magnetic fields 4 Malmqvist Celebration Palermo /
23 Summary and acknowledgments We have discussed Lieb s theory of DFT ground-state energy and universal density functional related by conjugation We have discussed Lieb maximizations Moreau Yosida regularization of energy and density functional We have discussed the adabatic connection ab-initio adiabatic-connection curves for benchmarking and modelling We have discussed DFT in magnetic fields B-dependent (BDFT) and current-dependent (CDFT) density functionals This work was supported by Norwegian Research Council European Research Council Helgaker et al. (CTCC, University of Oslo) Summary and acknowledgments Malmqvist Celebration Palermo / 23
The calculation of the universal density functional by Lieb maximization
The calculation of the universal density functional by Lieb maximization Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
More informationImporting ab-initio theory into DFT: Some applications of the Lieb variation principle
Importing ab-initio theory into DFT: Some applications of the Lieb variation principle Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department
More informationAb-initio studies of the adiabatic connection in density-functional theory
Ab-initio studies of the adiabatic connection in density-functional theory Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
More informationThe adiabatic connection
The adiabatic connection Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Dipartimento di Scienze
More informationThe Role of the Hohenberg Kohn Theorem in Density-Functional Theory
The Role of the Hohenberg Kohn Theorem in Density-Functional Theory T. Helgaker, U. E. Ekström, S. Kvaal, E. Sagvolden, A. M. Teale,, E. Tellgren Centre for Theoretical and Computational Chemistry (CTCC),
More informationDensity-functional theory
Density-functional theory Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 13th Sostrup Summer School Quantum Chemistry and Molecular
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Raman Centre for Atomic, Molecular and Optical
More informationFundamentals of Density-Functional Theory
Fundamentals of Density-Functional Theory Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway GdR CORREL Meeting 015 Amphi Sciences Naturelles,
More informationQuantum Chemistry in Magnetic Fields
Quantum Chemistry in Magnetic Fields Trygve Helgaker Hylleraas Centre of Quantum Molecular Sciences, Department of Chemistry, University of Oslo, Norway 11th Triennial Congress of the World Association
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationarxiv: v1 [physics.chem-ph] 29 Jan 2018
Density-Wavefunction Mapping in Degenerate Current-Density-Functional Theory arxiv:1801.09606v1 [physics.chem-ph] 9 Jan 018 Andre Laestadius and Erik I. Tellgren Hylleraas Centre for Quantum Molecular
More informationMolecules in Magnetic Fields
Molecules in Magnetic Fields Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo, Norway
More informationMolecular electronic structure in strong magnetic fields
Trygve Helgaker 1, Mark Hoffmann 1,2 Kai Molecules Lange in strong 1, Alessandro magnetic fields Soncini 1,3, andcctcc Erik Jackson Tellgren 211 1 (CTCC, 1 / 23 Uni Molecular electronic structure in strong
More informationMolecules in strong magnetic fields
Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway
More informationThe Nottingham eprints service makes this work by researchers of the University of Nottingham available open access under the following conditions.
Furness, James W. and Verbeke, Joachim and Tellgren, Erik I. and Stopkowicz, Stella and Ekström, Ulf and Helgaker, Trygve and Teale, Andrew M. (2015) Current density functional theory using meta-generalized
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationExcitation energies from density-functional theory some failures and successes. Trygve Helgaker
1 Excitation energies from density-functional theory some failures and successes Trygve Helgaker Centre for Theoretical and Computational Chemistry, University of Oslo, Norway Ola Lutnæs, Simen Reine,
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway European Summer School in Quantum Chemistry
More informationMolecular Magnetism. Magnetic Resonance Parameters. Trygve Helgaker
Molecular Magnetism Magnetic Resonance Parameters Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Laboratoire de Chimie Théorique,
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Hylleraas Centre, Department of Chemistry, University of Oslo, Norway and Centre for Advanced Study at the Norwegian Academy of Science and Letters, Oslo,
More informationOrbital dependent correlation potentials in ab initio density functional theory
Orbital dependent correlation potentials in ab initio density functional theory noniterative - one step - calculations Ireneusz Grabowski Institute of Physics Nicolaus Copernicus University Toruń, Poland
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationMolecular bonding in strong magnetic fields
Trygve Helgaker 1, Mark Hoffmann 1,2 Kai Molecules Lange in strong 1, Alessandro magnetic fields Soncini 1,3 AMAP,, andzurich, ErikJune Tellgren 1 4 212 1 (CTCC, 1 / 35 Uni Molecular bonding in strong
More informationChemical bonding in strong magnetic fields
Trygve Helgaker 1, Mark Hoffmann 1,2 Chemical Kai Lange bonding in 1, strong Alessandro magnetic fields Soncini 1,3, CMS212, and Erik JuneTellgren 24 27 212 1 (CTCC, 1 / 32 Uni Chemical bonding in strong
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry and
More informationMolecules in strong magnetic fields
Trygve Helgaker 1, Mark Hoffmann 1,2 Kai Molecules Lange in strong 1, Alessandro magnetic fields Soncini 1,3 Nottingham,, and Erik 3th October Tellgren 213 1 (CTCC, 1 / 27 Uni Molecules in strong magnetic
More informationDensity-functional theory in quantum chemistry. Trygve Helgaker. From Quarks to the Nuclear Many-Body Problem
1 Density-functional theory in quantum chemistry Trygve Helgaker Centre for Theoretical and Computational Chemistry, University of Oslo, Norway From Quarks to the Nuclear Many-Body Problem A conference
More informationMolecules in strong magnetic fields
Molecules in strong magnetic fields Trygve Helgaker, Kai Lange, Alessandro Soncini, and Erik Tellgren Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway
More informationDiamagnetism and Paramagnetism in Atoms and Molecules
Diamagnetism and Paramagnetism in Atoms and Molecules Trygve Helgaker Alex Borgoo, Maria Dimitrova, Jürgen Gauss, Florian Hampe, Christof Holzer, Wim Klopper, Trond Saue, Peter Schwerdtfeger, Stella Stopkowicz,
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More information1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More informationBasis sets for electron correlation
Basis sets for electron correlation Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway The 12th Sostrup Summer School Quantum Chemistry
More informationDifferentiable but exact formulation of density-functional theory
Differentiable but exact formulation of density-functional theory Simen Kvaal, 1, a) Ulf Ekström, 1 Andrew M. Teale, 2, 1 and Trygve Helgaker 1 1) Centre for Theoretical and Computational Chemistry, Department
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationElectronic structure theory: Fundamentals to frontiers. 2. Density functional theory
Electronic structure theory: Fundamentals to frontiers. 2. Density functional theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley
More informationSpring College on Computational Nanoscience May Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor
2145-25 Spring College on Computational Nanoscience 17-28 May 2010 Variational Principles, the Hellmann-Feynman Theorem, Density Functional Theor Stefano BARONI SISSA & CNR-IOM DEMOCRITOS Simulation Center
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationElectron Correlation - Methods beyond Hartree-Fock
Electron Correlation - Methods beyond Hartree-Fock how to approach chemical accuracy Alexander A. Auer Max-Planck-Institute for Chemical Energy Conversion, Mülheim September 4, 2014 MMER Summerschool 2014
More informationMolecular Magnetic Properties
Molecular Magnetic Properties Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway European Summer School in Quantum Chemistry
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationAdvanced Quantum Chemistry III: Part 3. Haruyuki Nakano. Kyushu University
Advanced Quantum Chemistry III: Part 3 Haruyuki Nakano Kyushu University 2013 Winter Term 1. Hartree-Fock theory Density Functional Theory 2. Hohenberg-Kohn theorem 3. Kohn-Sham method 4. Exchange-correlation
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
More informationMolecular Magnetism. Molecules in an External Magnetic Field. Trygve Helgaker
Molecular Magnetism Molecules in an External Magnetic Field Trygve Helgaker Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry, University of Oslo, Norway Laboratoire de
More informationDalton Quantum Chemistry Program
1 Quotation from home page: Dalton Quantum Chemistry Program Dalton QCP represents a powerful quantum chemistry program for the calculation of molecular properties with SCF, MP2, MCSCF or CC wave functions.
More informationarxiv: v1 [physics.chem-ph] 29 Sep 2017
Uniform magnetic fields in density-functional theory Erik I. Tellgren, Andre Laestadius, and Simen Kvaal Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O.
More informationMolecular Magnetic Properties. The 11th Sostrup Summer School. Quantum Chemistry and Molecular Properties July 4 16, 2010
1 Molecular Magnetic Properties The 11th Sostrup Summer School Quantum Chemistry and Molecular Properties July 4 16, 2010 Trygve Helgaker Centre for Theoretical and Computational Chemistry, Department
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationOrbital currents in the Colle-Salvetti correlation energy functional and the degeneracy problem. Abstract
Orbital currents in the Colle-Salvetti correlation energy functional and the degeneracy problem S. Pittalis 1, S. Kurth 1, S. Sharma 1,2 and E.K.U. Gross 1 1 Institut für Theoretische Physik, Freie Universität
More informationIndependent electrons in an effective potential
ABC of DFT Adiabatic approximation Independent electrons in an effective potential Hartree Fock Density Functional Theory MBPT - GW Density Functional Theory in a nutshell Every observable quantity of
More informationPreface Introduction to the electron liquid
Table of Preface page xvii 1 Introduction to the electron liquid 1 1.1 A tale of many electrons 1 1.2 Where the electrons roam: physical realizations of the electron liquid 5 1.2.1 Three dimensions 5 1.2.2
More informationIntroduction to density-functional theory. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut
More informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More information3: Density Functional Theory
The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London
More informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
More informationAb initio calculations for potential energy surfaces. D. Talbi GRAAL- Montpellier
Ab initio calculations for potential energy surfaces D. Talbi GRAAL- Montpellier A theoretical study of a reaction is a two step process I-Electronic calculations : techniques of quantum chemistry potential
More informationIntroduction to Computational Quantum Chemistry: Theory
Introduction to Computational Quantum Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC 3108 Course Lectures 2007 Introduction Hartree Fock Theory Configuration Interaction Lectures 1 Introduction
More informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationShort Course on Density Functional Theory and Applications VII. Hybrid, Range-Separated, and One-shot Functionals
Short Course on Density Functional Theory and Applications VII. Hybrid, Range-Separated, and One-shot Functionals Samuel B. Trickey Sept. 2008 Quantum Theory Project Dept. of Physics and Dept. of Chemistry
More informationDensity Func,onal Theory (Chapter 6, Jensen)
Chem 580: DFT Density Func,onal Theory (Chapter 6, Jensen) Hohenberg- Kohn Theorem (Phys. Rev., 136,B864 (1964)): For molecules with a non degenerate ground state, the ground state molecular energy and
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry Chemicum 4th floor vesa.hanninen@helsinki.fi September 10, 2013 Lecture 3. Electron correlation methods September
More informationDensity Functional Theory
Chemistry 380.37 Fall 2015 Dr. Jean M. Standard October 28, 2015 Density Functional Theory What is a Functional? A functional is a general mathematical quantity that represents a rule to convert a function
More informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationWalter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.
Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm
More informationExtending Kohn-Sham density-functional theory
Extending Kohn-Sham density-functional theory Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Email: julien.toulouse@upmc.fr Web page: www.lct.jussieu.fr/pagesperso/toulouse/ January
More informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
More informationDENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY
DENSITY FUNCTIONAL THEORY FOR NON-THEORISTS JOHN P. PERDEW DEPARTMENTS OF PHYSICS AND CHEMISTRY TEMPLE UNIVERSITY A TUTORIAL FOR PHYSICAL SCIENTISTS WHO MAY OR MAY NOT HATE EQUATIONS AND PROOFS REFERENCES
More informationarxiv: v2 [physics.chem-ph] 27 Apr 2018
Generalized Kohn Sham iteration on Banach spaces Andre Laestadius Hylleraas Centre for Quantum Molecular Sciences, Department of Chemistry, University of Oslo, P.O. Box 1033 Blindern, N-0315 Oslo, Norway
More informationTDDFT in Chemistry and Biochemistry III
TDDFT in Chemistry and Biochemistry III Dmitrij Rappoport Department of Chemistry and Chemical Biology Harvard University TDDFT Winter School Benasque, January 2010 Dmitrij Rappoport (Harvard U.) TDDFT
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More information2 Electronic structure theory
Electronic structure theory. Generalities.. Born-Oppenheimer approximation revisited In Sec..3 (lecture 3) the Born-Oppenheimer approximation was introduced (see also, for instance, [Tannor.]). We are
More informationComputational Chemistry I
Computational Chemistry I Text book Cramer: Essentials of Quantum Chemistry, Wiley (2 ed.) Chapter 3. Post Hartree-Fock methods (Cramer: chapter 7) There are many ways to improve the HF method. Most of
More informationThe importance of current contributions to shielding constants in density-functional theory
J our nal Name The importance of current contributions to shielding constants in density-functional theory Sarah Reimann,a Ulf Ekström,a Stella Stopkowicz,a Andrew M. Teale,a,b Alex Borgoo,a and Trygve
More information2.5 Time dependent density functional theory
.5 Time dependent density functional theory The main theorems of Kohn-Sham DFT state that: 1. Every observable is a functional of the density (Hohenger-Kohn theorem).. The density can be obtained within
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationCombining density-functional theory and many-body methods
Combining density-functional theory and many-body methods Julien Toulouse Université Pierre & Marie Curie and CNRS, Paris, France Vrije Universiteit Amsterdam, Netherlands November 2017 Outline 2/23 1
More information3/23/2010 More basics of DFT Kieron Burke and friends UC Irvine Physics and Chemistry References for ground-state DFT ABC of DFT, by KB and Rudy Magyar, http://dft.uci.edu A Primer in Density Functional
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationGeneralized generalized gradient approximation: An improved density-functional theory for accurate orbital eigenvalues
PHYSICAL REVIEW B VOLUME 55, NUMBER 24 15 JUNE 1997-II Generalized generalized gradient approximation: An improved density-functional theory for accurate orbital eigenvalues Xinlei Hua, Xiaojie Chen, and
More informationChemistry 4560/5560 Molecular Modeling Fall 2014
Final Exam Name:. User s guide: 1. Read questions carefully and make sure you understand them before answering (if not, ask). 2. Answer only the question that is asked, not a different question. 3. Unless
More informationAdvanced Electronic Structure Theory Density functional theory. Dr Fred Manby
Advanced Electronic Structure Theory Density functional theory Dr Fred Manby fred.manby@bris.ac.uk http://www.chm.bris.ac.uk/pt/manby/ 6 Strengths of DFT DFT is one of many theories used by (computational)
More informationWave function methods for the electronic Schrödinger equation
Wave function methods for the electronic Schrödinger equation Zürich 2008 DFG Reseach Center Matheon: Mathematics in Key Technologies A7: Numerical Discretization Methods in Quantum Chemistry DFG Priority
More informationDensity-functional theory at noninteger electron numbers
Density-functional theory at noninteger electron numbers Tim Gould 1, Julien Toulouse 2 1 Griffith University, Brisbane, Australia 2 Université Pierre & Marie Curie and CRS, Paris, France July 215 Introduction
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationIntroduction to Computational Chemistry: Theory
Introduction to Computational Chemistry: Theory Dr Andrew Gilbert Rm 118, Craig Building, RSC andrew.gilbert@anu.edu.au 3023 Course Lectures Introduction Hartree Fock Theory Basis Sets Lecture 1 1 Introduction
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
More informationAdiabatic-connection fluctuation-dissipation density-functional theory based on range separation
Adiabatic-connection fluctuation-dissipation density-functional theory based on range separation Julien Toulouse 1 I. Gerber 2, G. Jansen 3, A. Savin 1, W. Zhu 1, J. Ángyán 4 1 Laboratoire de Chimie Théorique,
More informationConvergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems
1 Convergence properties of the coupled-cluster method: the accurate calculation of molecular properties for light systems T. Helgaker Centre for Theoretical and Computational Chemistry, Department of
More informationBuilding a wavefunction within the Complete-Active. Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy
Building a wavefunction within the Complete-Active Active-Space Coupled-Cluster Cluster with Singles and Doubles formalism: straightforward description of quasidegeneracy Dmitry I. Lyakh (Karazin Kharkiv
More information7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig
Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space
More informationQUANTUM CHEMISTRY FOR TRANSITION METALS
QUANTUM CHEMISTRY FOR TRANSITION METALS Outline I Introduction II Correlation Static correlation effects MC methods DFT III Relativity Generalities From 4 to 1 components Effective core potential Outline
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationJack Simons, Henry Eyring Scientist and Professor Chemistry Department University of Utah
1. Born-Oppenheimer approx.- energy surfaces 2. Mean-field (Hartree-Fock) theory- orbitals 3. Pros and cons of HF- RHF, UHF 4. Beyond HF- why? 5. First, one usually does HF-how? 6. Basis sets and notations
More informationElectronic band structure, sx-lda, Hybrid DFT, LDA+U and all that. Keith Refson STFC Rutherford Appleton Laboratory
Electronic band structure, sx-lda, Hybrid DFT, LDA+U and all that Keith Refson STFC Rutherford Appleton Laboratory LDA/GGA DFT is good but... Naive LDA/GGA calculation severely underestimates band-gaps.
More informationKey concepts in Density Functional Theory
From the many body problem to the Kohn-Sham scheme ILM (LPMCN) CNRS, Université Lyon 1 - France European Theoretical Spectroscopy Facility (ETSF) December 12, 2012 Lyon Outline 1 The many-body problem
More information1. Thomas-Fermi method
1. Thomas-Fermi method We consider a system of N electrons in a stationary state, that would obey the stationary Schrödinger equation: h i m + 1 v(r i,r j ) Ψ(r 1,...,r N ) = E i Ψ(r 1,...,r N ). (1.1)
More informationIntroduction to multiconfigurational quantum chemistry. Emmanuel Fromager
Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Notations
More informationPrinciples of Quantum Mechanics
Principles of Quantum Mechanics - indistinguishability of particles: bosons & fermions bosons: total wavefunction is symmetric upon interchange of particle coordinates (space,spin) fermions: total wavefuncftion
More informationElectric properties of molecules
Electric properties of molecules For a molecule in a uniform electric fielde the Hamiltonian has the form: Ĥ(E) = Ĥ + E ˆµ x where we assume that the field is directed along the x axis and ˆµ x is the
More information