TUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION
|
|
- Ami Bridges
- 5 years ago
- Views:
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
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 informationDensity 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 informationHands-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 informationCrystal 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 informationab 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 informationBefore 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 informationSupporting 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 informationA 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 informationNeatest 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 informationTHERMOSTATS 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 informationUnit 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 informationProgramming 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 informationPhonon 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 informationQ 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 informationPredicting 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 informationLOWELL 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 informationA. 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 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 informationLOWELL 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 informationPhysics 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 informationMultiple 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 informationStructural 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 informationTheoretical 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 informationMECH 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 informationHomework 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 informationTheoretical 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 informationAn 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 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 informationarxiv: 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 informationLOWELL 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 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 informationDFT 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 informationAdvanced 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 informationTwo 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 informationAnswers 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 informationTutorial 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 informationIntroduction 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 informationMultiple 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 informationQuasi-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 informationStress 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 informationGeometry 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 "
3 5? #< q 8 2 / / ) 9 ) 2 ) > < _ / ] > ) 2 ) ) 5 > x > [ < > < ) > _ ] ]? <
More informationModule #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 informationChapter 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 informationBasic 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 informationPEAT 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 informationProperties 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 informationMANY 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 informationFrom 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 informationFinite-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 informationLecture 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 informationRigid 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 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 informationexciting 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 informationOWELL 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 informationStrain-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 informationThermodynamics 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 information1 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 informationNDT&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 informationAlgebraic 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
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 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]
(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 informationClosed-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 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 informationAbinit 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 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 informationUniversity 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 informationWave 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 informationMP464: 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 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 informationIf 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 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 informationStress, 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 informationHomogenization 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 informationVibrational 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 informationChapter 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 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 informationTopological 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 informationSIMULATION 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 informationChapter 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 informationProject 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 informationApplications 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 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 informationChemistry 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 informationWinmostar 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 informationequation 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 informationRotations 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 informationRotational 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 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 informationMicroscopic-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 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 informationarxiv: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 informationRepresenting 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 informationElectromagnetism 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 informationVibrational 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 informationPotentials, 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 informationKeble 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 informationMATH 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