Introduction to density functional perturbation theory for lattice dynamics
|
|
- Lewis Lewis
- 5 years ago
- Views:
Transcription
1 Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy)
2 Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion The phonon solution 2 ph.x 3 q2r.x Discrete Fourier transform in 1D matdyn.x
3 Description of a solid Equations of motion The phonon solution Let s consider a periodic solid. We indicate with R I = R µ + d s (1) the equilibrium positions of the atoms. R µ indicate the Bravais lattice vectors and d s the positions of the atoms in one unit cell (s = 1,..., N at ). We take N unit cells with Born-von Karman periodic boundary conditions. Ω is the volume of one cell and V = NΩ the volume of the solid. At time t, each atom is displaced from its equilibrium position. u I (t) is the displacement of the atom I.
4 Description of a solid Equations of motion The phonon solution Within the Born-Oppenheimer adiabatic approximation the nuclei move in a potential energy given by the total energy of the electron system calculated (for instance within DFT) at fixed nuclei. We call E tot (R I + u I ) (2) this energy. The electrons are assumed to be in the ground state for each nuclear configuration. If u I is small, we can expand E tot in a Taylor series with respect to u I. Within the harmonic approximation: E tot (R I +u I ) = E tot (R I )+ Iα E tot u Iα u Iα Iα,Jβ 2 E tot u Iα u Jβ u Iα u Jβ +... where the derivatives are calculated at u I = 0 and α and β indicate the three cartesian coordinates. (3)
5 Description of a solid Equations of motion The phonon solution Equations of motion At equilibrium E tot u Iα = 0, so the Hamiltonian of the ions becomes: H = Iα P 2 Iα 2M I Iα,Jβ 2 E tot u Iα u Jβ u Iα u Jβ (4) where P I are the momenta of the nuclei and M I their masses. The classical motion of the nuclei is given by the N 3 N at functions u Iα (t). These functions are the solutions of the Hamilton equations: u Iα = H P Iα Ṗ Iα = H u Iα (5)
6 Description of a solid Equations of motion The phonon solution Equations of motion-ii With our Hamiltonian: u Iα = P Iα M I or: Ṗ Iα = Jβ M I ü Iα = Jβ 2 E tot u Iα u Jβ u Jβ (6) 2 E tot u Iα u Jβ u Jβ (7)
7 Description of a solid Equations of motion The phonon solution The phonon solution We can search the solution in the form of a phonon. Let s introduce a vector q in the first Brillouin zone. For each q we can write: u µsα (t) = 1 Ms u sα (q)e i(qrµ ωqt) (8) where the time dependence is given by a simple exponential e iωqt and the displacement of the atoms in each cell identified by the Bravais lattice R µ can be obtained from the displacements of the atoms in one unit cell, for instance the one 1 that corresponds to R µ = 0: Ms u sα (q).
8 Description of a solid Equations of motion The phonon solution The phonon solution-ii Inserting this solution in the equations of motion and writing I = (µ, s), J = (ν, s ) we obtain an eigenvalue problem for the 3 N at variables u sα (q): ω 2 qu sα (q) = s β D sαs β(q)u s β(q) (9) where: D sαs β(q) = 1 Ms M s ν is the dynamical matrix of the solid. 2 E tot e iq(rν Rµ) (10) u µsα u νs β
9 ph.x Within DFT the ground state total energy of the solid, calculated at fixed nuclei, is: E tot = ψ i ψ i + V loc (r)ρ(r)d 3 r +E H [ρ]+e xc [ρ]+u II i (11) where ρ(r) is the density of the electron gas: ρ(r) = ψ i (r) 2 (12) i and ψ i are the solution of the Kohn and Sham equations. E H is the Hartree energy, E xc is the exchange and correlation energy and U II is the ion-ion interaction. According to the Hellmann-Feynman theorem, the first order derivative of the ground state energy with respect to an external parameter is: E tot λ = Vloc (r) λ ρ(r)d 3 r + U II λ (13)
10 ph.x Deriving with respect to a second parameter µ: 2 E tot λ = + 2 V loc (r) λ ρ(r)d 3 r + 2 U II λ Vloc (r) ρ(r) λ d 3 r (14) So the new quantity that we need to calculate is the charge density induced, at first order, by the perturbation: ρ(r) = i ψi (r) ψ i(r) + ψi (r) ψ i(r) (15) To fix the ideas we can think that λ = u µsα and µ = u νs β
11 ph.x The wavefunctions obey the following equation: [ 1 ] V KS (r) ψ i (r) = ε i ψ i (r) (16) where V KS = V loc (r) + V H (r) + V xc (r). V KS (r, µ) depends on µ so that also ψ i (r, µ), and ε i (µ) depend on µ. We can expand these quantities in a Taylor series: V KS (r, µ) = V KS (r, µ = 0) + µ V KS(r) +... ψ i (r, µ) = ψ i (r, µ = 0) + µ ψ i(r) +... ε i (µ) = ε i (µ = 0) + µ ε i +... (17)
12 ph.x Inserting these equations and keeping only the first order in µ we obtain: [ 1 ] ψi (r) V KS (r) ε i = V KS ψ i(r) + ε i ψ i(r) (18) where: V KS = V loc + V H + Vxc and V H V xc 1 ρ(r ) = r r d 3 r = dv xc ρ(r) dρ (19) depend self-consistently on the charge density induced by the perturbation.
13 ph.x The induced charge density depends only on P c ψ i where P c = 1 P v is the projector on the conduction bands and P v = i ψ i ψ i is the projector on the valence bands. In fact: ρ(r) = i ψi P (r) c ψ i(r) + ψi (r)p ψ i (r) c + i ψi P (r) v ψ i(r) + ψi (r)p ψ i (r) v (20) ρ(r) = i + ij ψi P (r) c ψ i(r) + ψi (r)p c ψ j (r)ψ i(r) ψ i (r) ( ψ i ψ j + ψ i ψ j ) (21)
14 ph.x DFPT Therefore we can solve the self-consistent linear system: [ 12 ] ψ i (r) 2 + V KS (r) ε i P c = P V KS c ψ i(r) (22) where and ρ(r) = i V KS = V loc + V H + V xc ψi P (r) c ψ i(r) + ψi (r)p ψ i (r) c (23) (24)
15 ph.x The program ph.x solves this self-consistent linear system for 3 N at perturbations at a fixed vector q. With ρ(r) for all the perturbations it calculates the dynamical matrix D sαs β(q) (25) at the given q. Diagonalizing this matrix we obtain 3 N at frequencies ω q. By repeating this procedure for several q we could plot ω q as a function of q and display the phonon dispersions. However, it is more convenient to adopt a different approach that requires the calculation of the dynamical matrix in a small set of points q.
16 q2r.x Discrete Fourier transform in 1D matdyn.x The dynamical matrix of the solid: D sαs β(q) = 1 Ms M s ν 2 E tot e iq(rν Rµ) (26) u µsα u νs β is a periodic function of q with D sαs β(q + G) = D sαs β(q) for any reciprocal lattice vector G. Moreover, since in a solid all Bravais lattice points are equivalent, it does not depend on µ. Eq.26 is a Fourier expansion of a three dimensional periodic function. We have Fourier components only at the discrete values R ν of the Bravais lattice and we can write: 1 Ms M s 2 E tot u µsα u νs β = Ω (2π) 3 d 3 qd sαs β(q)e iq(rν Rµ) (27)
17 q2r.x Discrete Fourier transform in 1D matdyn.x We can use the properties of the discrete Fourier transform and sample the integral in a uniform mesh of points q. This will give the interatomic force constants only for a set of R ν neighbors of R µ. The code q2r.x reads a set of dynamical matrices calculated in a uniform mesh of q points and calculates, using Eq. 27, the interatomic force constants for a few neighbors of the point R µ = 0.
18 q2r.x Discrete Fourier transform in 1D matdyn.x Let us consider a one dimensional periodic function f (x + a) = f (x) with period a. This function can be expanded in a Fourier series and will have a discrete set of Fourier components at the point k n = 2π a n, where n is an integer (positive, negative or zero). f (x) = n c n e iknx (28) where the coefficients of the expansion are: c n = 1 a a 0 f (x)e iknx dx (29) In general, if f (x) is a sufficiently smooth function, c n 0 at large n. Now suppose that we know f (x) only in a uniform set of N points x j = j x where x = a/n and j = 0,..., N 1,
19 q2r.x Discrete Fourier transform in 1D matdyn.x then we can calculate: c n = 1 N N 1 j=0 f (x j )e i 2π N nj (30) c n is a periodic function of n and c n+n = c n. So, if N is sufficiently large that c n = 0 when n N/2 then c n is a good approximation of c n for n < N/2 and the function f (x) = n=n/2 n= N/2 c n e iknx (31) is a good approximation of the function f (x) also on the points x different from x j.
20 q2r.x Discrete Fourier transform in 1D matdyn.x Therefore, if the dynamical matrix is a sufficiently smooth function of q and the interatomic force constants decay sufficiently rapidly in real space, we can use Eq. 26 to calculate the dynamical matrix at arbitrary q, limiting the sum to a few R ν neighbors of R µ = 0. The program matdyn.x reads the interatomic force constants calculated by q2r.x and calculates the dynamical matrices at an arbitrary q using Eq. 26.
21 q2r.x Discrete Fourier transform in 1D matdyn.x This procedure fails in two cases: In metals when there are Kohn anomalies. In this case D sαs β(q) is not a smooth function of q and the interatomic force constants are long range. In polar insulators where the atomic displacements generate long range electrostatic interactions and the dynamical matrix is non analytic for q 0. This case, however, can be dealt with by calculating the Born effective charges and the dielectric constant of the material.
Tutorial on DFPT and TD-DFPT: calculations of phonons and absorption spectra
Tutorial on DFPT and TD-DFPT: calculations of phonons and absorption spectra Iurii Timrov SISSA Scuola Internazionale Superiore di Studi Avanzati, Trieste Italy itimrov@sissa.it Computer modelling of materials
More informationPseudopotential generation and test by the ld1.x atomic code: an introduction
and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)
More informationPhonons In The Elk Code
Phonons In The Elk Code Kay Dewhurst, Sangeeta Sharma, Antonio Sanna, Hardy Gross Max-Planck-Institut für Mikrostrukturphysik, Halle If you are allowed to measure only one property of a material, make
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 informationQuantum Theory of Matter
Imperial College London Department of Physics Professor Ortwin Hess o.hess@imperial.ac.uk Quantum Theory of Matter Spring 014 1 Periodic Structures 1.1 Direct and Reciprocal Lattice (a) Show that the reciprocal
More informationPoisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
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 informationLattice dynamics. Javier Junquera. Philippe Ghosez. Andrei Postnikov
Lattice dynamics Javier Junquera Philippe Ghosez Andrei Postnikov Warm up: a little bit of notation Greek characters ( ) refer to atoms within the unit cell Latin characters ( ) refer to the different
More informationCHEM6085: Density Functional Theory
Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
More informationFundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering
Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)
More informationAb initio phonon calculations in mixed systems
Ab initio phonon calculations in mixed systems Andrei Postnikov apostnik@uos.de Outline: Experiment vs. ab initio theory Ways of theory: linear response and frozen phonon approaches Applications: Be x
More informationCP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
More informationElectronic Structure Methodology 1
Electronic Structure Methodology 1 Chris J. Pickard Lecture Two Working with Density Functional Theory In the last lecture we learnt how to write the total energy as a functional of the density n(r): E
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
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 informationLattice Vibrations. Chris J. Pickard. ω (cm -1 ) 200 W L Γ X W K K W
Lattice Vibrations Chris J. Pickard 500 400 300 ω (cm -1 ) 200 100 L K W X 0 W L Γ X W K The Breakdown of the Static Lattice Model The free electron model was refined by introducing a crystalline external
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More informationALMA: All-scale predictive design of heat management material structures
ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable
More informationTime-dependent density functional theory
Time-dependent density functional theory E.K.U. Gross Max-Planck Institute for Microstructure Physics OUTLINE LECTURE I Phenomena to be described by TDDFT Some generalities on functional theories LECTURE
More informationAdaptively Compressed Polarizability Operator
Adaptively Compressed Polarizability Operator For Accelerating Large Scale ab initio Phonon Calculations Ze Xu 1 Lin Lin 2 Lexing Ying 3 1,2 UC Berkeley 3 ICME, Stanford University BASCD, December 3rd
More information+ V = 0, j = 1,..., 3N (7.1) i =1, 3 m = 1, N, X mi, V X mi. i = 0 in the equilibrium. X mi X nj
7. Lattice dynamics 7.1 Basic consideratons In this section we ll follow a number of well-known textbooks, and only try to keep the notations consistent. Suppose we have a system of N atoms, m=1,...,n,
More informationThe PWcond code: Complex bands, transmission, and ballistic conductance
The PWcond code: Complex bands, transmission, and ballistic conductance SISSA and IOM-DEMOCRITOS Trieste (Italy) Outline 1 Ballistic transport: a few concepts 2 Complex band structure 3 Current of a Bloch
More informationReferences. Documentation Manuals Tutorials Publications
References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =
More informationChemistry 365: Normal Mode Analysis David Ronis McGill University
Chemistry 365: Normal Mode Analysis David Ronis McGill University 1. Quantum Mechanical Treatment Our starting point is the Schrodinger wave equation: Σ h 2 2 2m i N i=1 r 2 i + U( r 1,..., r N ) Ψ( r
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 informationDFT / SIESTA algorithms
DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang
More informationSolid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationSelf Consistent Cycle
Self Consistent Cycle Step 0 : defining your system namelist SYSTEM How to specify the System All periodic systems can be specified by a Bravais Lattice and and atomic basis How to specify the Bravais
More informationWavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts.
Wavefunction and electronic struture in solids: Bloch functions, Fermi level and other concepts. Silvia Casassa Università degli Studi di Torino July 12, 2017 Minnesota Workshop on ab initio MSC Symmetry
More informationThe electronic structure of materials 2 - DFT
Quantum mechanics 2 - Lecture 9 December 19, 2012 1 Density functional theory (DFT) 2 Literature Contents 1 Density functional theory (DFT) 2 Literature Historical background The beginnings: L. de Broglie
More informationElectronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch
Electronic Structure Theory for Periodic Systems: The Concepts Christian Ratsch Institute for Pure and Applied Mathematics and Department of Mathematics, UCLA Motivation There are 10 20 atoms in 1 mm 3
More informationKronig-Penney model. 2m dx. Solutions can be found in region I and region II Match boundary conditions
Kronig-Penney model I II d V( x) m dx E Solutions can be found in region I and region II Match boundary conditions Linear differential equations with periodic coefficients Have exponentially decaying solutions,
More informationMD simulation: output
Properties MD simulation: output Trajectory of atoms positions: e. g. diffusion, mass transport velocities: e. g. v-v autocorrelation spectrum Energies temperature displacement fluctuations Mean square
More informationEtudes ab initio des effets de la température sur le spectre optique des semi-conducteurs
Etudes ab initio des effets de la température sur le spectre optique des semi-conducteurs Woyten Tielesch Laboratoire de Physique de la Matière Condensée et Nanostructures January 27, 2011 Outline Outline
More informationIntegrated Computational Materials Engineering Education
Integrated Computational Materials Engineering Education Lecture on Density Functional Theory An Introduction Mark Asta Dept. of Materials Science and Engineering, University of California, Berkeley &
More informationClassical Theory of Harmonic Crystals
Classical Theory of Harmonic Crystals HARMONIC APPROXIMATION The Hamiltonian of the crystal is expressed in terms of the kinetic energies of atoms and the potential energy. In calculating the potential
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationElectron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals
Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Richard M. Martin UIUC Lecture at Summer School Hands-on introduction to Electronic Structure Materials Computation Center University
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 information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 19: Motion of Electronic Wavepackets Outline Review of Last Time Detailed Look at the Translation Operator Electronic Wavepackets Effective Mass Theorem
More informationMany-Body Perturbation Theory. Lucia Reining, Fabien Bruneval
, Fabien Bruneval Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France European Theoretical Spectroscopy Facility (ETSF) Belfast, 27.6.2007 Outline 1 Reminder 2 Perturbation Theory
More informationCondensed matter physics FKA091
Condensed matter physics FKA091 Ermin Malic Department of Physics Chalmers University of Technology Henrik Johannesson Department of Physics University of Gothenburg Teaching assistants: Roland Jago &
More informationPractical Guide to Density Functional Theory (DFT)
Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between
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 informationPhonon wavefunctions and electron phonon interactions in semiconductors
Phonon wavefunctions and electron phonon interactions in semiconductors Bartomeu Monserrat bm418@cam.ac.uk University of Cambridge Quantum Monte Carlo in the Apuan Alps VII QMC in the Apuan Alps VII Bartomeu
More informationDensity Functional Theory (DFT)
Density Functional Theory (DFT) An Introduction by A.I. Al-Sharif Irbid, Aug, 2 nd, 2009 Density Functional Theory Revolutionized our approach to the electronic structure of atoms, molecules and solid
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 informationFirst-principles Molecular Dynamics Simulations
First-principles Molecular Dynamics Simulations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org MICCoM Computational School, Jul
More informationModel operator approach to calculations of the Lamb shifts in relativistic many-electron atoms
Model operator approach to calculations of the Lamb shifts in relativistic many-electron atoms Vladimir Shabaev a in collaboration with Ilya Tupitsyn a and Vladimir Yerokhin b a St. Petersburg State University
More informationOn the adaptive finite element analysis of the Kohn-Sham equations
On the adaptive finite element analysis of the Kohn-Sham equations Denis Davydov, Toby Young, Paul Steinmann Denis Davydov, LTM, Erlangen, Germany August 2015 Denis Davydov, LTM, Erlangen, Germany College
More informationSolid State Physics II Lattice Dynamics and Heat Capacity
SEOUL NATIONAL UNIVERSITY SCHOOL OF PHYSICS http://phya.snu.ac.kr/ ssphy2/ SPRING SEMESTER 2004 Chapter 3 Solid State Physics II Lattice Dynamics and Heat Capacity Jaejun Yu jyu@snu.ac.kr http://phya.snu.ac.kr/
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 informationBand Structure Calculations; Electronic and Optical Properties
; Electronic and Optical Properties Stewart Clark University of Outline Introduction to band structures Calculating band structures using Castep Calculating optical properties Examples results Some applications
More informationCourse SS18. Organization
0-0 Course SS18 Stefan.Goedecker@unibas.ch Organization Lecture with script: http://comphys.unibas.ch/teaching.htm Simple exercises: Traditional analytic problems and numerical problems (on your own laptop?)
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 8: Lattice Waves in 1D Monatomic Crystals Outline Overview of Lattice Vibrations so far Models for Vibrations in Discrete 1-D Lattice Example of Nearest
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More information2. The Schrödinger equation for one-particle problems. 5. Atoms and the periodic table of chemical elements
1 Historical introduction The Schrödinger equation for one-particle problems 3 Mathematical tools for quantum chemistry 4 The postulates of quantum mechanics 5 Atoms and the periodic table of chemical
More informationQuantum Master Equations for the Electron Transfer Problem
20/01/2010 Quantum Master Equations for the Electron Transfer Problem Seminarvortrag Dekohaerenz und Dissipation in Quantensystemen Antonio A. Gentile The general transport problem in micro/mesoscopic
More informationAnharmonic energy in periodic systems
Anharmonic energy in periodic systems Bartomeu Monserrat University of Cambridge Electronic Structure Discussion Group 13 March 213 Vibrational properties overview Harmonic phonons are a very good approximation.
More informationMolecular Mechanics: The Ab Initio Foundation
Molecular Mechanics: The Ab Initio Foundation Ju Li GEM4 Summer School 2006 Cell and Molecular Mechanics in BioMedicine August 7 18, 2006, MIT, Cambridge, MA, USA 2 Outline Why are electrons quantum? Born-Oppenheimer
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 informationPseudopotentials: design, testing, typical errors
Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 29: Electron-phonon Scattering Outline Bloch Electron Scattering Deformation Potential Scattering LCAO Estimation of Deformation Potential Matrix Element
More informationBorn-Oppenheimer Approximation
Born-Oppenheimer Approximation Adiabatic Assumption: Nuclei move so much more slowly than electron that the electrons that the electrons are assumed to be obtained if the nuclear kinetic energy is ignored,
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 informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
More informationDFT in practice : Part II. Ersen Mete
pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials
More informationNumerical construction of Wannier functions
July 12, 2017 Internship Tutor: A. Levitt School Tutor: É. Cancès 1/27 Introduction Context: Describe electrical properties of crystals (insulator, conductor, semi-conductor). Applications in electronics
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 informationLecture contents. A few concepts from Quantum Mechanics. Tight-binding model Solid state physics review
Lecture contents A few concepts from Quantum Mechanics Particle in a well Two wells: QM perturbation theory Many wells (atoms) BAND formation Tight-binding model Solid state physics review Approximations
More informationElectron levels in a periodic potential: general remarks. Daniele Toffoli January 11, / 46
Electron levels in a periodic potential: general remarks Daniele Toffoli January 11, 2017 1 / 46 Outline 1 Mathematical tools 2 The periodic potential 3 Bloch s theorem and Born-von Karman boundary conditions
More information7 Lattice Vibrations. 7.1 Introduction
7 Lattice Vibrations 7.1 Introduction Up to this point in the lecture, the crystal lattice (lattice vectors & position of basis atoms) was always assumed to be immobile, i.e., atomic displacements away
More informationApplications of Matrix Functions Part III: Quantum Chemistry
Applications of Matrix Functions Part III: Quantum Chemistry Emory University Department of Mathematics and Computer Science Atlanta, GA 30322, USA Prologue The main purpose of this lecture is to present
More informationarxiv: v1 [cond-mat.str-el] 18 Jul 2007
arxiv:0707.2704v1 [cond-mat.str-el] 18 Jul 2007 Determination of the Mott insulating transition by the multi-reference density functional theory 1. Introduction K. Kusakabe Graduate School of Engineering
More informationMicroscopic-macroscopic connection
icroscopic-macroscopic connection Valérie Véniard ETSF France Laboratoire des Solides Irradiés Ecole Polytechniue CEA-DS CNRS 91128 Palaiseau France Outline 1. Introduction: which uantities do we need
More informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
More informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More informationDownloaded 05/09/17 to Redistribution subject to SIAM license or copyright; see
MULTISCALE MODEL. SIMUL. Vol. 15, No. 1, pp. 29 55 c 2017 Society for Industrial and Applied Mathematics ADAPTIVELY COMPRESSED POLARIZABILITY OPERATOR FOR ACCELERATING LARGE SCALE AB INITIO PHONON CALCULATIONS
More informationElectronic Structure of Crystalline Solids
Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem
More informationFirst-principles dynamics of electrons and phonons
First-principles dynamics of electrons and phonons Marco Bernardi 1 Department of Applied Physics and Materials Science, Steele Laboratory, California Institute of Technology, Pasadena, CA 91125, USA.
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 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 informationMODULE 2: QUANTUM MECHANICS. Practice: Quantum ESPRESSO
MODULE 2: QUANTUM MECHANICS Practice: Quantum ESPRESSO I. What is Quantum ESPRESSO? 2 DFT software PW-DFT, PP, US-PP, PAW http://www.quantum-espresso.org FREE PW-DFT, PP, PAW http://www.abinit.org FREE
More informationDFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory
DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,
More informationSecond Quantization Method for Bosons
Second Quantization Method for Bosons Hartree-Fock-based methods cannot describe the effects of the classical image potential (cf. fig. 1) because HF is a mean-field theory. DFF-LDA is not able either
More informationThree Most Important Topics (MIT) Today
Three Most Important Topics (MIT) Today Electrons in periodic potential Energy gap nearly free electron Bloch Theorem Energy gap tight binding Chapter 1 1 Electrons in Periodic Potential We now know the
More informationPhonons II - Thermal Properties (Kittel Ch. 5)
Phonons II - Thermal Properties (Kittel Ch. 5) Heat Capacity C T 3 Approaches classical limit 3 N k B T Physics 460 F 2006 Lect 10 1 Outline What are thermal properties? Fundamental law for probabilities
More informationPhonons and lattice dynamics
Chapter Phonons and lattice dynamics. Vibration modes of a cluster Consider a cluster or a molecule formed of an assembly of atoms bound due to a specific potential. First, the structure must be relaxed
More informationPhysics 541: Condensed Matter Physics
Physics 541: Condensed Matter Physics In-class Midterm Exam Wednesday, October 26, 2011 / 14:00 15:20 / CCIS 4-285 Student s Name: Instructions There are 23 questions. You should attempt all of them. Mark
More informationAb initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More informationThe Hamiltonian. 2.1 The Schrödinger Equation
2 The Hamiltonian Any theoretical description has to start with the definition of the system under consideration and a determination of the fundamental interactions present in the system. This information
More informationConical Intersections. Spiridoula Matsika
Conical Intersections Spiridoula Matsika The Born-Oppenheimer approximation Energy TS Nuclear coordinate R ν The study of chemical systems is based on the separation of nuclear and electronic motion The
More informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
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 informationPhonon Dispersion, Interatomic Force Constants Thermodynamic Quantities
Phonon Dispersion, Interatomic Force Constants Thermodynamic Quantities Umesh V. Waghmare Theoretical Sciences Unit J N C A S R Bangalore ICMR OUTLINE Vibrations and interatomic force constants (IFC) Extended
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 informationLinear response theory and TDDFT
Linear response theory and TDDFT Claudio Attaccalite http://abineel.grenoble.cnrs.fr/ CECAM Yambo School 2013 (Lausanne) Motivations: +- hν Absorption Spectroscopy Many Body Effects!!! Motivations(II):Absorption
More information