TUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION

Size: px
Start display at page:

Download "TUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION"

Transcription

1 TUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION 1

2 INVESTIGATED SYSTEM: Silicon, diamond structure Electronic and 0K properties see W. Huhn, Tutorial 2, Wednesday August 2 2

3 THE HARMONIC APPROXIMATION E ({R 0 + R}) E ({R 0 })+ i R 0 R i X 2 j R0 R i R j 3

4 THE HARMONIC APPROXIMATION E ({R 0 + R}) E ({R 0 })+ i R 0 R i X 2 j R0 R i R j Static Equilibrium Energy from DFT Hessian Φij 4

5 THE HARMONIC APPROXIMATION E ({R 0 + R}) E ({R 0 })+ i R 0 R i X 2 j R0 R i R j Static Equilibrium Energy from DFT Hessian Φij Determine Hessian aka the Harmonic Force Constants Φij: from Density-Functional Perturbation Theory S. Baroni, P. Giannozzi, and A. Testa, Phys. Rev. Lett. 58, 1861 (1987) & S. Baroni, et al., Rev. Mod. Phys. 73, 515 (2001). from Finite Differences K. Kunc, and R. M. Martin, Phys. Rev. Lett. 48, 406 (1982) & K. Parlinski, Z. Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 78, 4063 (1997). 5

6 THE HARMONIC APPROXIMATION E ({R 0 + R}) E ({R 0 })+ i R 0 R i X 2 j R0 R i R j Static Equilibrium Energy from DFT Hessian Φij Determine Hessian aka the Harmonic Force Constants Φij: from Density-Functional Perturbation Theory phonopy-fhi-aims S. Baroni, P. Giannozzi, and A. Testa, Phys. Rev. Lett. 58, 1861 (1987) & S. A. Baroni, Togo, F. et Oba, al., Rev. and Mod. I. Tanaka, Phys. Phys. 73, 515 Rev. (2001). B 78, (2008). from Finite Differences K. Kunc, and R. M. Martin, Phys. Rev. Lett. 48, 406 (1982) & K. Parlinski, Z. Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 78, 4063 (1997). 6

7 THE FINITE DIFFERENCE APPROACH K. Parlinski, Z. Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 78, 4063 (1997). A. Togo, F. Oba, and I. Tanaka, Phys. Rev. B 78, (2008). Finite differences using normalized displacements d: ij j R 0 F j F j(r 0 i + " d i R 0 " 7

8 THE FINITE DIFFERENCE APPROACH K. Parlinski, Z. Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 78, 4063 (1997). A. Togo, F. Oba, and I. Tanaka, Phys. Rev. B 78, (2008). Finite differences using normalized displacements d: ij j R 0 F j F j(r 0 i + " d i R 0 " Example: Diamond Si (2 atoms in the basis): 0 xx yx zx xx yx zx xy yy zy xy yy zy xz yz zz xz yz zz xx yx zx xx yx zx xy yy zy xy yy zy xz yz zz xz yz zz 1 C A Hessian has 36 entries: 6 displacements d required 8

9 THE FINITE DIFFERENCE APPROACH K. Parlinski, Z. Q. Li, and Y. Kawazoe, Phys. Rev. Lett. 78, 4063 (1997). A. Togo, F. Oba, and I. Tanaka, Phys. Rev. B 78, (2008). Finite differences using normalized displacements d: ij j R 0 F j F j(r 0 i + " d i R 0 " 0 xx yx zx xx yx zx xy yy zy xy yy zy xz yz zz xz yz zz Example: Diamond Si (2 atoms in the basis): xx yx zx xx yx zx xy yy zy xy yy zy Space Group Analysis xz yz zz xz yz zz 1 C A 0 xx xy xx 0 xy 0 xx xy xx 0 xy 0 xz yz xx xz yz xx xx yz xz xx yz xz xy 0 xx 0 xy xx xy 0 xx 0 xy xx 1 C A Hessian has 5 unique, non-zero entries: Only 1 displacement d required 9

10 VIBRATIONS IN A CRYSTAL 101 e.g. N. W Ashcroft and N. D. Mermin, Solid State Physics (1976) also see Oliver Hofmann, Tuesday August 1 Dynamical matrix: D i 0 j 0 (q) = X j e i (q (R 0 j R 0 j 0 )) p Mi 0M j 0 i 0 j Equation of Motion becomes an Eigenvalue Problem: D(q) [ (q)] =! 2 (q) [ (q)] Analytical Solution in Real Space: Superposition of Harmonic Oscillations R j (t) = R 0 j + Re X s A s p Mi e i (q (R 0 j R 0 j 0 )! s (q)t) [ s (q)] j 0!

11 (a) Edit your file control.in so that it contains the following lines phonon displacement 0.01 (b) Run phonopy-fhi-aims by typing phonopy-fhi-aims (c) Change into the directory phonopy-fhi-aims-displacement-01 and run FHI-aims: cd phonopy-fhi-aims-displacement-01 mpirun -np 4 aims.x > phonopy-fhi-aims-displacement-01.out (d) Change into parent directory and run phonopy-fhi-aims again cd.. phonopy-fhi-aims

12 CONVERGING THE SUPERCELL Fourier Transform can be truncated since Φij = 0 D i 0 j 0 (q) = X j e i (q (R 0 j R 0 j 0 )) p Mi 0M j 0 i 0 j for large Rjj 0 = Rj 0 - Rj 0 phonon displacement 0.01 phonon supercell k_grid phonon displacement 0.01 phonon supercell k_grid 2 2 2

13 CONVERGING THE SUPERCELL Fourier Transform can be truncated since Φij = 0 D i 0 j 0 (q) = X j e i (q (R 0 j R 0 j 0 )) p Mi 0M j 0 i 0 j for large Rjj 0 = Rj 0 - Rj 0 phonon displacement 0.01 phonon supercell k_grid phonon displacement 0.01 phonon supercell k_grid To achieve convergence, it is essential to have a consistent description of the electronic structure for all supercell sizes: #atoms #k-points constant 13

14 CONVERGING THE SUPERCELL Fourier Transform can be truncated since Φij = 0 D i 0 j 0 (q) = X j e i (q (R 0 j R 0 j 0 )) p Mi 0M j 0 i 0 j for large Rjj 0 = Rj 0 - Rj 0 Oblique Cell: Not all Cartesian directions are treated consistently! Cubic ( spherical ) Cell: Consistent assessment of all cartesian directions! 14

15 VIBRATIONAL BAND STRUCTURE # control.in : Plot vibrational band structure phonon band Gamma Delta phonon band Delta X phonon band X W phonon band W K phonon band K Gamma phonon band Gamma Lambda phonon band Lambda L ω (cm -1 ) Γ X W K Γ Λ L

16 VIBRATIONAL DENSITY OF STATES g(!) = X s Z dq (2 ) 3 (!!(q)) = X s Z!(q)=! ds 1 (2 ) 3 r!(q) # control.in : Plot vibrational density of states phonon dos g(ω) (a.u.) ω (cm -1 )

17 VIBRATIONAL DENSITY OF STATES g(!) = X s Z dq (2 ) 3 (!!(q)) = X s 600 Z!(q)=! ds 1 (2 ) 3 r!(q) ω (cm -1 ) g(ω) (a.u.) Γ X W K Γ Λ L ω (cm -1 )

18 THE HARMONIC FREE ENERGY F ha (T ) = E({R 0 }) Z Static Equilibrium Energy + + Z d d g( ) ~ 2 Zero-point vibration g( ) k B T ln 1 e ~ k B T Thermally induced vibrations 18

19 FREE ENERGY AND HEAT CAPACITY F ha (ev) c V (k B ) F (T ) C V = T = V 2 V T (K) 19

20 THE QUASI-HARMONIC APPROXIMATION 20

21 THE HARMONIC APPROXIMATION H = X i T i X i,j ij R i R =0 Lattice expansion vanishes in the harmonic approximation. THE QUASI-HARMONIC APPROACH H = X i T i X i,j ij(v ) R i R 6=0 Assess lattice expansion by explicitly accounting for the volume dependence of the Hessian.

22 THE QUASI-HARMONIC APPROACH (free) energy Lattice constant a0 can be determined from Birch-Murnaghan fit of E(a0) cf. William Huhn, Practical Session 2 EDFT lattice constant a0

23 THE QUASI-HARMONIC APPROACH (free) energy +Fha(0K) Lattice constant a0 can be determined from Birch-Murnaghan fit of E(a0) cf. William Huhn, Practical Session 2 Add vibrational free energy for each individual value of a0 EDFT lattice constant a0 23

24 THE QUASI-HARMONIC APPROACH (free) energy a0(t) EDFT lattice constant a0 +Fha(T2) +Fha(T1) +Fha(0K) Lattice constant a0 can be determined from Birch-Murnaghan fit of E(a0) cf. William Huhn, Practical Session 2 Add vibrational free energy for each individual value of a0 Repeat for each temperature 0K < T1 < T2 Birch-Murnaghan fits for each individual temperature allow to determine temperature dependence of lattice constant a0(t). 24

25 EXERCISE 3 LATTICE 5e-06 (T )= 1 p 4e-06 3e-06 α (1/K) 2e-06 1e e Temperature (K) 25

26 EXERCISE 3 LATTICE 5e-06 (T )= 1 p 4e-06 3e-06 α (1/K) 2e-06 1e-06 < 0? 0-1e Temperature (K) 26

27 ELECTRON-PHONON COUPLING 27

28 BAND GAP RENORMALIZATION Electronic band gaps often exhibit a distinct temperature dependence Linear extrapolation yields the bare gap at 0K, i.e., the gap for immobile nuclei (classical limit) Linear Extrapolation Actual band gap at 0K differs from the bare gap: Band gap renormalization 28 M. Cardona, Solid State Comm. 133, 3 (2005).

29 ELECTRON-PHONON COUPLING conduction band "(k) bare band gap valence band 29

30 ELECTRON-PHONON COUPLING conduction band conduction band "(k) bare band gap valence band "(k) band gap valence band 30

31 BAND GAP TEMPERATURE DEPENDENCE What is the physical mechanism? Exercise 4: Lattice Expansion? Exercise 5: Atomic Motion? Use results from exercise 3 Molecular Dynamics 31

32 HARMONIC MOLECULAR DYNAMICS M I RI (t) = r Ri E DFT M I RI (t) = X j ij R j Harmonic Approximation # Molecular Dynamics MD_MB_init MD_time_step MD_schedule MD_segment 5.0 NVT_parrinello harmonic_potential_only fc_constants.dat MD_segment 20.0 NVT_parrinello harmonic_potential_only fc_constants.dat # Equilibration # Sample phase space 32

33 WARNING: In the following exercises, the computational settings, in particular the reciprocal space grid (tag k_grid), the basis set, and supercells, have been chosen to allow a rapid computation of the exercises in the limited time and within the CPU resources available during the tutorial session. In a real production calculation, the reciprocal space grid, the basis set, and the supercells would all have to be converged with much more care, but the qualitative trends hold already with the present settings Happy Computing! Christian Carbogno Mohsen Yarmohammadi 33 Maja Lenz Florian Knoop

34 LATTICE EXPANSION S. Biernacki and M. Scheffler, Phys. Rev. Lett. 63, 290 (1989). Z Free energy definition: F ha (T! 0) / d! g(!) ~! equilibrium V Frequencies lowered at frequenncy (cm -1 ) expanded V larger volumes!

35 LATTICE EXPANSION S. Biernacki and M. Scheffler, Phys. Rev. Lett. 63, 290 (1989). Z Free energy definition: F ha (T! 0) / d! g(!) ~! equilibrium V Frequencies lowered at frequenncy (cm -1 ) expanded V larger volumes! Acoustic frequencies increased at 0 larger volumes! 35

TUTORIAL 6: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION

TUTORIAL 6: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION Hands-On Tutorial Workshop, July 29 th 2014 TUTORIAL 6: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION Christian Carbogno & Manuel Schöttler Fritz-Haber-Institut der Max-Planck-Gesellschaft,

More information

Density functional theory and beyond: Computational materials science for real materials Los Angeles, July 21 August 1, 2014

Density functional theory and beyond: Computational materials science for real materials Los Angeles, July 21 August 1, 2014 Density functional theory and beyond: Computational materials science for real materials Los Angeles, July 21 August 1, 2014 R L e y E N e x Tutorial VI: Phonons, Lattice Expansion, and Band-gap Renormalization

More information

Hands-on Workshop: First-principles simulations of molecules and materials: Berlin, July 13 - July 23, 2015

Hands-on Workshop: First-principles simulations of molecules and materials: Berlin, July 13 - July 23, 2015 Hands-on Workshop: First-principles simulations of molecules and materials: Berlin, July 13 - July 23, 2015 R L e y E N e x Tutorial IV: Phonons, Lattice Expansion, and Band-gap Renormalization Manuscript

More information

Crystal Relaxation, Elasticity, and Lattice Dynamics

Crystal Relaxation, Elasticity, and Lattice Dynamics http://exciting-code.org Crystal Relaxation, Elasticity, and Lattice Dynamics Pasquale Pavone Humboldt-Universität zu Berlin http://exciting-code.org PART I: Structure Optimization Pasquale Pavone Humboldt-Universität

More information

ab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 30, 2009 This afternoon s plan

ab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 30, 2009 This afternoon s plan ab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 3, 29 This afternoon s plan introductory talk Phonons: harmonic vibrations for solids Phonons:

More information

Before we start: Important setup of your Computer

Before we start: Important setup of your Computer Before we start: Important setup of your Computer change directory: cd /afs/ictp/public/shared/smr2475./setup-config.sh logout login again 1 st Tutorial: The Basics of DFT Lydia Nemec and Oliver T. Hofmann

More information

Supporting information. The Unusual and the Expected in the Si/C Phase Diagram. Guoying Gao, N. W. Ashcroft and Roald Hoffmann.

Supporting information. The Unusual and the Expected in the Si/C Phase Diagram. Guoying Gao, N. W. Ashcroft and Roald Hoffmann. Supporting information The Unusual and the Expected in the Si/C Phase Diagram Guoying Gao, N. W. Ashcroft and Roald Hoffmann Table of Contents Computational Methods...S1 Hypothetical Structures for Si

More information

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any

A DARK GREY P O N T, with a Switch Tail, and a small Star on the Forehead. Any Y Y Y X X «/ YY Y Y ««Y x ) & \ & & } # Y \#$& / Y Y X» \\ / X X X x & Y Y X «q «z \x» = q Y # % \ & [ & Z \ & { + % ) / / «q zy» / & / / / & x x X / % % ) Y x X Y $ Z % Y Y x x } / % «] «] # z» & Y X»

More information

Neatest and Promptest Manner. E d i t u r ami rul)lihher. FOIt THE CIIILDIIES'. Trifles.

Neatest and Promptest Manner. E d i t u r ami rul)lihher. FOIt THE CIIILDIIES'. Trifles. » ~ $ ) 7 x X ) / ( 8 2 X 39 ««x» ««! «! / x? \» «({? «» q «(? (?? x! «? 8? ( z x x q? ) «q q q ) x z x 69 7( X X ( 3»«! ( ~«x ««x ) (» «8 4 X «4 «4 «8 X «x «(» X) ()»» «X «97 X X X 4 ( 86) x) ( ) z z

More information

THERMOSTATS AND THERMAL TRANSPORT IN SOLIDS

THERMOSTATS AND THERMAL TRANSPORT IN SOLIDS Hands-On Tutorial Workshop, July 19 th 2011 THERMOSTATS AND THERMAL TRANSPORT IN SOLIDS Christian Carbogno FRITZ-HABER-INSTITUT DER MAX-PLANCK-GESELLSCHAFT, BERLIN - GERMANY Hands-On Tutorial Workshop,

More information

Unit IV State of stress in Three Dimensions

Unit IV State of stress in Three Dimensions Unit IV State of stress in Three Dimensions State of stress in Three Dimensions References Punmia B.C.,"Theory of Structures" (SMTS) Vol II, Laxmi Publishing Pvt Ltd, New Delhi 2004. Rattan.S.S., "Strength

More information

Programming Project 2: Harmonic Vibrational Frequencies

Programming Project 2: Harmonic Vibrational Frequencies Programming Project 2: Harmonic Vibrational Frequencies Center for Computational Chemistry University of Georgia Athens, Georgia 30602 Summer 2012 1 Introduction This is the second programming project

More information

Phonon calculations with SCAN

Phonon calculations with SCAN Workshop on the SCAN density functional: Fundamentals, practices, and extensions Temple university, Philadelphia May 18th, 2017 Hands-on tutorial 3 Phonon calculations with SCAN Yubo Zhang and Jianwei

More information

Q SON,' (ESTABLISHED 1879L

Q SON,' (ESTABLISHED 1879L ( < 5(? Q 5 9 7 00 9 0 < 6 z 97 ( # ) $ x 6 < ( ) ( ( 6( ( ) ( $ z 0 z z 0 ) { ( % 69% ( ) x 7 97 z ) 7 ) ( ) 6 0 0 97 )( 0 x 7 97 5 6 ( ) 0 6 ) 5 ) 0 ) 9%5 z» 0 97 «6 6» 96? 0 96 5 0 ( ) ( ) 0 x 6 0

More information

Predicting the Structure of Solids by DFT

Predicting the Structure of Solids by DFT Questions? Hudson Hall 235 or Hudson Hall 1111 Predicting the Structure of Solids by DFT Hands-On Instructions Contents 1 Cohesive Energy for Bulk Phases of Si 11 Setting up the Structures 12 Structure-Dependent

More information

LOWELL WEEKLY JOURNAL.

LOWELL WEEKLY JOURNAL. Y $ Y Y 7 27 Y 2» x 7»» 2» q» ~ [ } q q $ $ 6 2 2 2 2 2 2 7 q > Y» Y >» / Y» ) Y» < Y»» _»» < Y > Y Y < )»» >» > ) >» >> >Y x x )»» > Y Y >>»» }> ) Y < >» /» Y x» > / x /»»»»» >» >» >»» > > >» < Y /~ >

More information

A. H. Hall, 33, 35 &37, Lendoi

A. H. Hall, 33, 35 &37, Lendoi 7 X x > - z Z - ----»»x - % x x» [> Q - ) < % - - 7»- -Q 9 Q # 5 - z -> Q x > z»- ~» - x " < z Q q»» > X»? Q ~ - - % % < - < - - 7 - x -X - -- 6 97 9

More information

Lattice dynamics. Javier Junquera. Philippe Ghosez. Andrei Postnikov

Lattice 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 information

LOWELL WEEKLY JOURNAL

LOWELL WEEKLY JOURNAL Y -» $ 5 Y 7 Y Y -Y- Q x Q» 75»»/ q } # ]»\ - - $ { Q» / X x»»- 3 q $ 9 ) Y q - 5 5 3 3 3 7 Q q - - Q _»»/Q Y - 9 - - - )- [ X 7» -» - )»? / /? Q Y»» # X Q» - -?» Q ) Q \ Q - - - 3? 7» -? #»»» 7 - / Q

More information

Physics with Neutrons I, WS 2015/2016. Lecture 11, MLZ is a cooperation between:

Physics with Neutrons I, WS 2015/2016. Lecture 11, MLZ is a cooperation between: Physics with Neutrons I, WS 2015/2016 Lecture 11, 11.1.2016 MLZ is a cooperation between: Organization Exam (after winter term) Registration: via TUM-Online between 16.11.2015 15.1.2015 Email: sebastian.muehlbauer@frm2.tum.de

More information

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015

Multiple Integrals and Vector Calculus (Oxford Physics) Synopsis and Problem Sets; Hilary 2015 Multiple Integrals and Vector Calculus (Oxford Physics) Ramin Golestanian Synopsis and Problem Sets; Hilary 215 The outline of the material, which will be covered in 14 lectures, is as follows: 1. Introduction

More information

Structural Calculations phase stability, surfaces, interfaces etc

Structural Calculations phase stability, surfaces, interfaces etc Structural Calculations phase stability, surfaces, interfaces etc Keith Refson STFC Rutherford Appleton Laboratory September 19, 2007 Phase Equilibrium 2 Energy-Volume curves..................................................................

More information

Theoretical Concepts of Spin-Orbit Splitting

Theoretical Concepts of Spin-Orbit Splitting Chapter 9 Theoretical Concepts of Spin-Orbit Splitting 9.1 Free-electron model In order to understand the basic origin of spin-orbit coupling at the surface of a crystal, it is a natural starting point

More information

MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso

MECH 5312 Solid Mechanics II. Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso MECH 5312 Solid Mechanics II Dr. Calvin M. Stewart Department of Mechanical Engineering The University of Texas at El Paso Table of Contents Preliminary Math Concept of Stress Stress Components Equilibrium

More information

Homework 1/Solutions. Graded Exercises

Homework 1/Solutions. Graded Exercises MTH 310-3 Abstract Algebra I and Number Theory S18 Homework 1/Solutions Graded Exercises Exercise 1. Below are parts of the addition table and parts of the multiplication table of a ring. Complete both

More information

Theoretical Material Science: Electronic structure theory at the computer Exercise 14: Brillouin zone integration

Theoretical Material Science: Electronic structure theory at the computer Exercise 14: Brillouin zone integration Theoretical Material Science: Electronic structure theory at the computer Exercise 14: Brillouin zone integration Prepared by Lydia Nemec and Volker Blum Berlin, May 2012 Some rules on expected documentation

More information

An introduction to Solid State NMR and its Interactions

An introduction to Solid State NMR and its Interactions An introduction to Solid State NMR and its Interactions From tensor to NMR spectra CECAM Tutorial September 9 Calculation of Solid-State NMR Parameters Using the GIPAW Method Thibault Charpentier - CEA

More information

Phonon wavefunctions and electron phonon interactions in semiconductors

Phonon 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 information

arxiv: v1 [cond-mat.mtrl-sci] 1 Mar 2011

arxiv: v1 [cond-mat.mtrl-sci] 1 Mar 2011 Phonon-phonon interactions in transition metals Laurent Chaput, 1 Atsushi Togo,,3 Isao Tanaka, 3,4 and Gilles Hug 1 Institut Jean Lamour, UMR CNRS 7198, Nancy Université, Bd. des Aiguillettes, BP 3, 5456

More information

LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort

LOWELL JOURNAL. MUST APOLOGIZE. such communication with the shore as Is m i Boimhle, noewwary and proper for the comfort - 7 7 Z 8 q ) V x - X > q - < Y Y X V - z - - - - V - V - q \ - q q < -- V - - - x - - V q > x - x q - x q - x - - - 7 -» - - - - 6 q x - > - - x - - - x- - - q q - V - x - - ( Y q Y7 - >»> - x Y - ] [

More information

DFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory

DFT 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 information

DFT EXERCISES. FELIPE CERVANTES SODI January 2006

DFT EXERCISES. FELIPE CERVANTES SODI January 2006 DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a

More information

Advanced Spectroscopy. Dr. P. Hunt Rm 167 (Chemistry) web-site:

Advanced Spectroscopy. Dr. P. Hunt Rm 167 (Chemistry) web-site: Advanced Spectroscopy Dr. P. Hunt p.hunt@imperial.ac.uk Rm 167 (Chemistry) web-site: http://www.ch.ic.ac.uk/hunt Maths! Coordinate transformations rotations! example 18.1 p501 whole chapter on Matrices

More information

Two Posts to Fill On School Board

Two Posts to Fill On School Board Y Y 9 86 4 4 qz 86 x : ( ) z 7 854 Y x 4 z z x x 4 87 88 Y 5 x q x 8 Y 8 x x : 6 ; : 5 x ; 4 ( z ; ( ) ) x ; z 94 ; x 3 3 3 5 94 ; ; ; ; 3 x : 5 89 q ; ; x ; x ; ; x : ; ; ; ; ; ; 87 47% : () : / : 83

More information

Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014

Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Answers Quantum Chemistry NWI-MOL406 G. C. Groenenboom and G. A. de Wijs, HG00.307, 8:30-11:30, 21 jan 2014 Question 1: Basis sets Consider the split valence SV3-21G one electron basis set for formaldehyde

More information

Tutorial on DFPT and TD-DFPT: calculations of phonons and absorption spectra

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 information

Introduction to Condensed Matter Physics

Introduction to Condensed Matter Physics Introduction to Condensed Matter Physics Elasticity M.P. Vaughan Overview Overview of elasticity Classical description of elasticity Speed of sound Strain Stress Young s modulus Shear modulus Poisson ratio

More information

Multiple Integrals and Vector Calculus: Synopsis

Multiple Integrals and Vector Calculus: Synopsis Multiple Integrals and Vector Calculus: Synopsis Hilary Term 28: 14 lectures. Steve Rawlings. 1. Vectors - recap of basic principles. Things which are (and are not) vectors. Differentiation and integration

More information

Quasi-Harmonic Theory of Thermal Expansion

Quasi-Harmonic Theory of Thermal Expansion Chapter 5 Quasi-Harmonic Theory of Thermal Expansion 5.1 Introduction The quasi-harmonic approximation is a computationally efficient method for evaluating thermal properties of materials. Planes and Manosa

More information

Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation

Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Stress equilibrium in southern California from Maxwell stress function models fit to both earthquake data and a quasi-static dynamic simulation Peter Bird Dept. of Earth, Planetary, and Space Sciences

More information

Geometry Optimisation

Geometry Optimisation Geometry Optimisation Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Motivation Overview of Talk Background

More information

' Liberty and Umou Ono and Inseparablo "

' Liberty and Umou Ono and Inseparablo 3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <

More information

Module #3. Transformation of stresses in 3-D READING LIST. DIETER: Ch. 2, pp Ch. 3 in Roesler Ch. 2 in McClintock and Argon Ch.

Module #3. Transformation of stresses in 3-D READING LIST. DIETER: Ch. 2, pp Ch. 3 in Roesler Ch. 2 in McClintock and Argon Ch. HOMEWORK From Dieter -3, -4, 3-7 Module #3 Transformation of stresses in 3-D READING LIST DIETER: Ch., pp. 7-36 Ch. 3 in Roesler Ch. in McClintock and Argon Ch. 7 in Edelglass The Stress Tensor z z x O

More information

Chapter 5 Phonons II Thermal Properties

Chapter 5 Phonons II Thermal Properties Chapter 5 Phonons II Thermal Properties Phonon Heat Capacity < n k,p > is the thermal equilibrium occupancy of phonon wavevector K and polarization p, Total energy at k B T, U = Σ Σ < n k,p > ħ k, p Plank

More information

Basic Equations of Elasticity

Basic Equations of Elasticity A Basic Equations of Elasticity A.1 STRESS The state of stress at any point in a loaded bo is defined completely in terms of the nine components of stress: σ xx,σ yy,σ zz,σ xy,σ yx,σ yz,σ zy,σ zx,andσ

More information

PEAT SEISMOLOGY Lecture 2: Continuum mechanics

PEAT SEISMOLOGY Lecture 2: Continuum mechanics PEAT8002 - SEISMOLOGY Lecture 2: Continuum mechanics Nick Rawlinson Research School of Earth Sciences Australian National University Strain Strain is the formal description of the change in shape of a

More information

Properties of the stress tensor

Properties of the stress tensor Appendix C Properties of the stress tensor Some of the basic properties of the stress tensor and traction vector are reviewed in the following. C.1 The traction vector Let us assume that the state of stress

More information

MANY BILLS OF CONCERN TO PUBLIC

MANY BILLS OF CONCERN TO PUBLIC - 6 8 9-6 8 9 6 9 XXX 4 > -? - 8 9 x 4 z ) - -! x - x - - X - - - - - x 00 - - - - - x z - - - x x - x - - - - - ) x - - - - - - 0 > - 000-90 - - 4 0 x 00 - -? z 8 & x - - 8? > 9 - - - - 64 49 9 x - -

More information

From Electrons to Materials Properties

From Electrons to Materials Properties From Electrons to Materials Properties: DFT for Engineers and Materials Scientists Funk Fachgebiet Werkstoffe des Bauwesens und Bauchemie Mar-17, 2016 From Electrons to Materials Properties Density Functional

More information

Finite-temperature equation of state. T ln 2sinh h" '-

Finite-temperature equation of state. T ln 2sinh h '- Finite-temperature equation of state * $ F(V,T) = U 0 + k B T ln 2sinh h" '- #, & )/ + % 2k B T (. 0 # Compute vibrational modes, frequencies Evaluate at a given volume V Compute F at various temperatures

More information

Lecture 8. Stress Strain in Multi-dimension

Lecture 8. Stress Strain in Multi-dimension Lecture 8. Stress Strain in Multi-dimension Module. General Field Equations General Field Equations [] Equilibrium Equations in Elastic bodies xx x y z yx zx f x 0, etc [2] Kinematics xx u x x,etc. [3]

More information

Rigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient

Rigid body simulation. Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body dynamics Rigid body simulation Once we consider an object with spatial extent, particle system simulation is no longer sufficient Rigid body simulation Unconstrained system no contact Constrained

More information

Electronic Structure Theory for Periodic Systems: The Concepts. Christian Ratsch

Electronic 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 information

exciting in a nutshell

exciting in a nutshell http://exciting-code.org exciting in a nutshell Pasquale Pavone Humboldt-Universität zu Berlin http://exciting-code.org exciting in a (coco)nutshell Pasquale Pavone Humboldt-Universität zu Berlin Outline

More information

OWELL WEEKLY JOURNAL

OWELL WEEKLY JOURNAL Y \»< - } Y Y Y & #»»» q ] q»»»>) & - - - } ) x ( - { Y» & ( x - (» & )< - Y X - & Q Q» 3 - x Q Y 6 \Y > Y Y X 3 3-9 33 x - - / - -»- --

More information

Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity

Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity Torino, Italy, September 4-9, 2016 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it

More information

Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations

Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations, USA, July 9-14, 2017 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it 2017 Outline -

More information

1 Hooke s law, stiffness, and compliance

1 Hooke s law, stiffness, and compliance Non-quilibrium Continuum Physics TA session #5 TA: Yohai Bar Sinai 3.04.206 Linear elasticity I This TA session is the first of three at least, maybe more) in which we ll dive deep deep into linear elasticity

More information

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16.

NDT&E Methods: UT. VJ Technologies CAVITY INSPECTION. Nondestructive Testing & Evaluation TPU Lecture Course 2015/16. CAVITY INSPECTION NDT&E Methods: UT VJ Technologies NDT&E Methods: UT 6. NDT&E: Introduction to Methods 6.1. Ultrasonic Testing: Basics of Elasto-Dynamics 6.2. Principles of Measurement 6.3. The Pulse-Echo

More information

Algebraic Expressions

Algebraic Expressions Algebraic Expressions 1. Expressions are formed from variables and constants. 2. Terms are added to form expressions. Terms themselves are formed as product of factors. 3. Expressions that contain exactly

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

+ 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 information

The quasi-harmonic approximation (QHA)

The quasi-harmonic approximation (QHA) The quasi-harmonic approximation (QHA) M. Palumbo 19/01/2017 Trieste, Italy Limitations of the harmonic approximation E tot(r I, u I )=E tot(r I )+ I,α E tot u u Iα + 1 2 E tot Iα 2 I,α u u Iα u Iα u Jβ

More information

(2) A two-dimensional solid has an electron energy band of the form, . [1]

(2) A two-dimensional solid has an electron energy band of the form, . [1] (1) The figure shows a two-dimensional periodic lattice, containing A atoms (white) and B atoms (black). The section of lattice shown makes a 3a 4a rectangle, as shown (measured from A atom to A atom).

More information

Closed-Form Solution Of Absolute Orientation Using Unit Quaternions

Closed-Form Solution Of Absolute Orientation Using Unit Quaternions Closed-Form Solution Of Absolute Orientation Using Unit Berthold K. P. Horn Department of Computer and Information Sciences November 11, 2004 Outline 1 Introduction 2 3 The Problem Given: two sets of corresponding

More information

Phonon Dispersion, Interatomic Force Constants Thermodynamic Quantities

Phonon 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 information

Abinit electron phonon interaction calculations for geniuses

Abinit electron phonon interaction calculations for geniuses Abinit electron phonon interaction calculations for geniuses Matthieu Verstraete, Bin Xu, and Momar Diakhate Unité Nanomat, Dept of Physics, Université de Liège, Belgium 1 I. INTRODUCTION This document

More information

Phonons In The Elk Code

Phonons 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 information

University of Chinese Academy of Sciences, Beijing , People s Republic of China,

University of Chinese Academy of Sciences, Beijing , People s Republic of China, SiC 2 Siligraphene and Nanotubes: Novel Donor Materials in Excitonic Solar Cell Liu-Jiang Zhou,, Yong-Fan Zhang, Li-Ming Wu *, State Key Laboratory of Structural Chemistry, Fujian Institute of Research

More information

Wave and Elasticity Equations

Wave and Elasticity Equations 1 Wave and lasticity quations Now let us consider the vibrating string problem which is modeled by the one-dimensional wave equation. Suppose that a taut string is suspended by its extremes at the points

More information

MP464: Solid State Physics Problem Sheet

MP464: Solid State Physics Problem Sheet MP464: Solid State Physics Problem Sheet 1 Write down primitive lattice vectors for the -dimensional rectangular lattice, with sides a and b in the x and y-directions respectively, and a face-centred rectangular

More information

6.730 Physics for Solid State Applications

6.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 information

If the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc.

If the symmetry axes of a uniform symmetric body coincide with the coordinate axes, the products of inertia (Ixy etc. Prof. O. B. Wright, Autumn 007 Mechanics Lecture 9 More on rigid bodies, coupled vibrations Principal axes of the inertia tensor If the symmetry axes of a uniform symmetric body coincide with the coordinate

More information

Electronic Structure of Crystalline Solids

Electronic 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 information

Stress, Strain, Mohr s Circle

Stress, Strain, Mohr s Circle Stress, Strain, Mohr s Circle The fundamental quantities in solid mechanics are stresses and strains. In accordance with the continuum mechanics assumption, the molecular structure of materials is neglected

More information

Homogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections. Andrew Hall October 16, 2015

Homogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections. Andrew Hall October 16, 2015 Homogenization Methods for Full Core Solution of the Pn Transport Equations with 3-D Cross Sections Andrew Hall October 16, 2015 Outline Resource-Renewable Boiling Water Reactor (RBWR) Current Neutron

More information

Vibrational states of molecules. Diatomic molecules Polyatomic molecules

Vibrational states of molecules. Diatomic molecules Polyatomic molecules Vibrational states of molecules Diatomic molecules Polyatomic molecules Diatomic molecules V v 1 v 0 Re Q Birge-Sponer plot The solution of the Schrödinger equation can be solved analytically for the

More information

Chapter 2 Governing Equations

Chapter 2 Governing Equations Chapter Governing Equations Abstract In this chapter fundamental governing equations for propagation of a harmonic disturbance on the surface of an elastic half-space is presented. The elastic media is

More information

Ab initio phonon calculations in mixed systems

Ab 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 information

Topological insulator part I: Phenomena

Topological insulator part I: Phenomena Phys60.nb 5 Topological insulator part I: Phenomena (Part II and Part III discusses how to understand a topological insluator based band-structure theory and gauge theory) (Part IV discusses more complicated

More information

SIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS

SIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS SIMULATION AND INTERPRETATION OF BOREHOLE GEOPHYSICAL MEASUREMENTS USING hp FINTE ELEMENTS hp-fem team: D. Pardo, M. J. Nam, L. Demkowicz, C. Torres-Verdín, V. M. Calo, M. Paszynski, and P. J. Matuszyk

More information

Chapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set

Chapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless

More information

Project Report: Band Structure of GaAs using k.p-theory

Project Report: Band Structure of GaAs using k.p-theory Proect Report: Band Structure of GaAs using k.p-theory Austin Irish Mikael Thorström December 12th 2017 1 Introduction The obective of the proect was to calculate the band structure of both strained and

More information

Applications of Eigenvalues & Eigenvectors

Applications of Eigenvalues & Eigenvectors Applications of Eigenvalues & Eigenvectors Louie L. Yaw Walla Walla University Engineering Department For Linear Algebra Class November 17, 214 Outline 1 The eigenvalue/eigenvector problem 2 Principal

More information

7 Lattice Vibrations. 7.1 Introduction

7 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 information

Chemistry 431. NC State University. Lecture 17. Vibrational Spectroscopy

Chemistry 431. NC State University. Lecture 17. Vibrational Spectroscopy Chemistry 43 Lecture 7 Vibrational Spectroscopy NC State University The Dipole Moment Expansion The permanent dipole moment of a molecule oscillates about an equilibrium value as the molecule vibrates.

More information

Winmostar tutorial LAMMPS Melting point V X-Ability Co,. Ltd. 2017/8/17

Winmostar tutorial LAMMPS Melting point V X-Ability Co,. Ltd. 2017/8/17 Winmostar tutorial LAMMPS Melting point V7.025 X-Ability Co,. Ltd. question@winmostar.com Contents Configure I. Build solid phase II. Equilibration of solid phase III. Equilibration of liquid phase IV.

More information

equation of state and bulk modulus of AlP, AlAs and AlSb semiconductor compounds

equation of state and bulk modulus of AlP, AlAs and AlSb semiconductor compounds PRAMANA c Indian Academy of Sciences Vol. 64, No. 1 journal of January 2005 physics pp. 153 158 Total energy, equation of state and bulk modulus of AlP, AlAs and AlSb semiconductors A R JIVANI, H J TRIVEDI,

More information

Rotations and vibrations of polyatomic molecules

Rotations and vibrations of polyatomic molecules Rotations and vibrations of polyatomic molecules When the potential energy surface V( R 1, R 2,..., R N ) is known we can compute the energy levels of the molecule. These levels can be an effect of: Rotation

More information

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep.

Rotational Motion. Chapter 4. P. J. Grandinetti. Sep. 1, Chem P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. Rotational Motion Chapter 4 P. J. Grandinetti Chem. 4300 Sep. 1, 2017 P. J. Grandinetti (Chem. 4300) Rotational Motion Sep. 1, 2017 1 / 76 Angular Momentum The angular momentum of a particle with respect

More information

ALMA: All-scale predictive design of heat management material structures

ALMA: 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 information

Microscopic-Macroscopic connection. Silvana Botti

Microscopic-Macroscopic connection. Silvana Botti relating experiment and theory European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre for Computational

More information

Lecture contents. A few concepts from Quantum Mechanics. Tight-binding model Solid state physics review

Lecture 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 information

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 14 May 2003

arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 14 May 2003 LA-UR-3-239 arxiv:cond-mat/35331v1 [cond-mat.mtrl-sci] 14 May 23 Thermal Stabilization of the HCP Phase in Titanium Sven P. Rudin 1, M. D. Jones 2, and R. C. Albers 1 1 Los Alamos National Laboratory,

More information

Representing High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks

Representing High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks Energy Landscapes, Chemnitz 200 Representing High-Dimensional Potential-Energy Surfaces by Artificial Neural Networks Jörg Behler Lehrstuhl für Theoretische Chemie Ruhr-Universität Bochum D-44780 Bochum,

More information

Electromagnetism II Lecture 7

Electromagnetism II Lecture 7 Electromagnetism II Lecture 7 Instructor: Andrei Sirenko sirenko@njit.edu Spring 13 Thursdays 1 pm 4 pm Spring 13, NJIT 1 Previous Lecture: Conservation Laws Previous Lecture: EM waves Normal incidence

More information

Vibrational Analysis in Gaussian

Vibrational Analysis in Gaussian Vibrational Analysis in Gaussian Joseph W. Ochterski, Ph.D. help@gaussian.com October 9, 1999 Minor corrections 14 June 018 Abstract One of the most commonly asked questions about Gaussian is What is the

More information

Potentials, periodicity

Potentials, periodicity Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium

More information

Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering

Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering Tomi Johnson Keble College - Hilary 2012 Section VI: Condensed matter physics Tutorial 2 - Lattices and scattering Please leave your work in the Clarendon laboratory s J pigeon hole by 5pm on Monday of

More information

MATH 19520/51 Class 5

MATH 19520/51 Class 5 MATH 19520/51 Class 5 Minh-Tam Trinh University of Chicago 2017-10-04 1 Definition of partial derivatives. 2 Geometry of partial derivatives. 3 Higher derivatives. 4 Definition of a partial differential

More information