ab initio Lattice Vibrations: Calculating the Thermal Expansion Coeffcient Felix Hanke & Martin Fuchs June 30, 2009 This afternoon s plan
|
|
- Harvey Cross
- 6 years ago
- Views:
Transcription
1 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: how Thermodynamics with phonons practical exercise Obtain & assess phonons for Si Calculate & understand thermal expansion
2 Recap: Molecular vibrations (non-periodic) Newton s equations for small displacements of atom!, coordinate i about PES minimum: m α ẍ α,i = F α,i ({x β,j }) m α ẍ α,i + 2 E (δx α,i ) (δx β,j ) δx β,j β,j ω 2 [ m α δx α,i ]= β,j In symmetric matrix form 1 mα m β 2 E (δx α,i ) (δx β,j ) [ m j δx β,j ] The harmonic lattice Translational symmetry (PBC s) leads to k-dependence One eigenvalue problem for each k, giving 3Natoms modes ω 2 n(k)ξ n (k) =D(k)ξ n (k) D is the dynamic matrix, contains information from all supercells R, all pairs of atoms! and ", and all pairs of coordinates i and j: D αβ,ij = R e ik R mi m j 2 E r (R=),αi r (R),βj
3 Phonon dispersion relation Diamond fcc conventional cell 2 atoms per primitive cell = 6 phonon branches a = 3.5Å Recall diatomic molecules: 2 atoms = 6 degrees of freedom! =! =! =!vib 3 translations 2 rotations 1 vibration How about solids? 1 5 Frequency " (cm -1 )15 Phonon dispersion relation Diamond fcc conventional cell 2 atoms per primitive cell = 6 phonon branches optical branches! X W K! L k-vector acoustic branches a = 3.5Å DOS
4 Phonon density of states DOS = number of phonon modes per unit frequency per unit cell g(ω) = δ(ω ω n (k)) exp [ (ω ω n(k)) 2 ] 2σ 2 n,k n,k smearing # makes it feasible to Converge smearing vs plot (same as electronic DOS!) k-point sampling! numerical computation: frequency interval [$+%$/2,$-%$/2] with discrete k-point sampling g(ω) = w ω+ ω/2 k dω exp [ (ω ω n(k)) 2 ] 2π 2σ 2 n,k ω ω/2 The Direct Method approximation: finite interaction distance, use finite supercells D αβ,ij = R e ik R mi m j 2 E r (R=),αi r (R),βj
5 Harmonic approximation: Thermodynamics Approximate thermal effects with independent harmonic oscillators: one for each mode F = k,n = k,n ln Z(k,N) [ ωk,n 2 ] + k B T ln[1 exp( ω k,n /k B T )] Practically, one tends to use the density of states: = dω g(ω) [ ] ω 2 + k BT ln[1 exp( ω/k B T )] Quasiarmonic approximation: Changing lattice constants Alternative to explicit treatment of each phonon at all a s Anharmonic properties implicitly through dependence on lattice constant Treat many lattice constants harmonically & minimize free energy over all a s
6 Lattice Constants: Zero Point Vibrations how does the DOS change as fct of lattice constant? does this have any effect on the lattice? Density of states a = 3.3 Å a = 3.5 Å a = 3.7 Å diamond a Phonon frequency! (cm -1 ) Energy (ev) 1.5 Lattice Constants: Zero Point Vibrations Decreasing phonon free energy at larger a: Lattice constant including ZPE is larger than T= result ZPE = Å Å dω g(ω) ω 2 diamond Lattice constant (Å) no ZPE ZPE only sum
7 Lattice constants: Temperature effects Murnaghan fits of Ftot(a,T) at different temperatures lattice constant & Bulk modulus a (Å) diamond quasiharmonic no phonons B (MBar) quasiharmonic no phonons Temperature (K) " (K -1 ) Thermal expansion coefficient need to minimize F(a,T) at a given temperature to find lattice constant F (a, T )=E el (a)+f ph (a, T ) Differentiate to find α(t )= 1 a Watch the scale!!! 4!1-6 3!1-6 2!1-6 1!1-6 da dt diamond Temperature (K)
8 Heat capacity: cv Computed from free energy c v (T ) = T ds dt = T 2 F V T 2 V = dω g(ω) ( ω)2 exp( ω/k B T ) k B T 2 (exp( ω/k B T ) 1) 2 Grüneisen parameter: coupling to! α = γc v 3B c v (k B /unit cell) diamond Temperature (K) Breakdown of the harmonic approximation High temperatures: anharmonic excitations soft modes - dynamically stabilized structures V ( ξ) imaginary frequency Quantum effects, eg in hydrogen bonded crystals ξ
9 Beyond the quasiharmonic approximation MD to describe complete PES beyond quasiharmonic approximation F = E el + F qh + F anh + F vac thermodynamic integration: slowly switching on anh terms from qh solution Determine from self-consistently optimizing vacancy volume Grabowski, Ismer, Hickel & Neugebauer, Phys. Rev. B (29) Today s tutorial: Thermal properties of Silicon this talk afternoon exercise
10 Schedule for this tutorial (A) simple phonon dispersion relation (B) converging the supercell 45 min 45 min (C) zero point vibrations 8 min (D) thermal expansion 4 min Exercise A: calculations of phonons with FHI-aims use Si (diamond structure) lattice constant = 5.44 Å, light species defaults Supercell size 2x2x2 Compute: phonon dispersion relation for the directions "-X-W-K-"-L density of states, specific heat cv converge DOS see also: FHI-aims manual, section 4.5
11 Exercise A: calculations of phonons with FHI-aims Phonon calculation works similar to vibrations: run aims.phonons.workshop.mpi.pl in the directory containing your control.in and geometry.in files. Output phonon_band_structure.dat - same format as e-bands phonon_dos.dat - frequency, density of states phonon_free_energy.dat - T, F(T), U(T), cv(t), Svib(T) output stream - status reports phonon_workdir/ - working files & restart info Attention Please specify a SINGLE PRIMITIVE CELL ONLY! The script does all the necessary copying & displacing. Exercise A: calculations of phonons with FHI-aims Phonon dispersion calculation completely driven by phonon keyword in control.in (see FHI-aims manual, section 3.5) supercell size - phonon supercell (for now) DOS specification - phonon dos Free energy - phonon free_energy k-grid: hand set to match supercell size, i.e. k_grid Band structure: phonon band <start> <end> <Npoints> <sname> <ename> working directory: /usr/local/aimsfiles/tutorial6/exercise_a
12 Exercise A: Solutions Frequency " (cm -1 ) ! X W K! L c v (k B /unit cell) k-vector Temperature (K) DOS 4 2 Exercise B: Supercell convergence Compute phonon dispersion and DOS for the supercell sizes 4x4x4, 6x6x6 Use the converged DOS settings from the last exercise. All necessary files are in directory /usr/local/aimsfiles/tutorial6/exercise_b NOTE: We have provided partial control.in and geometry.in files as well as ALL of the required DFT output for the remainder of this tutorial
13 Exercise B: solutions Phonon frequency! (cm -1 ) x2x2 4x4x4 6x6x6! X W K! L k-vector Exercise B: solutions DOS (arb units) 2x2x2 4x4x4 6x6x6 2 4 Phonon frequency (cm -1 )
14 Exercise C: zero point vibrations Using data for the lattice constants provided in the directory /usr/local/aimsfiles/tutorial6/exercise_c+d For each lattice constant provided, calculate the total energy of a single unit cell using a 12x12x12 k-point grid. Store the result in the format 4x4x4_a*.***/output.single_point_energy Again, for each lattice constant, calculate the phonon free energy & specific heat for a temperature range from K to 8K in 81 steps (***) Extract the ZPE from the T= free energies from file phonon_free_energy.dat in each directory produce a Murnaghan fit with and without ZPE and plot the two fits & the resulting lattice constants (***) These exact settings are important for exercise D Exercise C: solutions Energy (mev) ZPE E(a), no ZPE E(a), incl ZPE Lattice constant a (Å)
15 Exercise D: compute a(t), cv &! Again, using the data provided directory /usr/local/aimsfiles/tutorial6/exercise_c+d use the script eval_alpha.sh to calculate the lattice constant at different temperatures & the thermal expansion coefficient eval_alpha.sh has to be started in the directory exercise_c+d and requires the EXACT specifications of phonon free_energy from the last exercise. find cv(t) at the optimal lattice constant a=5.44å Compare with previously calculated values Exercise D: solutions c v (k B /unit cell) Specific heat at constant volume Temperature (K)
16 Exercise D: solutions lattice constant a(t) & thermal expansion coefficient! " (K -1 ) a (Å) ! Temperature (K) Density of states Density oof states (arb units) a = 5.4 Å a = 5.45 Å a = 5.5 Å frequency reordering 25 5 Phonon frequency! (cm -1 ) a = 5.4 Å a = 5.45 Å a = 5.5 Å
17 Phonon frequency " (cm -1 ) Phonon band structure! X W K! L k-vector DOS Negative thermal expansion coefficient: Why optical phonon bands become softer acoustic bands become harder in some part of the BZ
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 informationTUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION
TUTORIAL 8: PHONONS, LATTICE EXPANSION, AND BAND-GAP RENORMALIZATION 1 INVESTIGATED SYSTEM: Silicon, diamond structure Electronic and 0K properties see W. Huhn, Tutorial 2, Wednesday August 2 2 THE HARMONIC
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 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 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 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 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 informationIntroduction to solid state physics
PHYS 342/555 Introduction to solid state physics Instructor: Dr. Pengcheng Dai Professor of Physics The University of Tennessee (Room 407A, Nielsen, 974-1509) Chapter 5: Thermal properties Lecture in pdf
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 informationPHONON HEAT CAPACITY
Solid State Physics PHONON HEAT CAPACITY Lecture 11 A.H. Harker Physics and Astronomy UCL 4.5 Experimental Specific Heats Element Z A C p Element Z A C p J K 1 mol 1 J K 1 mol 1 Lithium 3 6.94 24.77 Rhenium
More informationAn Introduction to Lattice Vibrations
An Introduction to Lattice Vibrations Andreas Wacker 1 Mathematical Physics, Lund University November 3, 2015 1 Introduction Ideally, the atoms in a crystal are positioned in a regular manner following
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 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 informationSMARTMET project: Towards breaking the inverse ductility-strength relation
SMARTMET project: Towards breaking the inverse ductility-strength relation B. Grabowski, C. Tasan SMARTMET ERC advanced grant 3.8 Mio Euro for 5 years (Raabe/Neugebauer) Adaptive Structural Materials group
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 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 informationLecture 12: Phonon heat capacity
Lecture 12: Phonon heat capacity Review o Phonon dispersion relations o Quantum nature of waves in solids Phonon heat capacity o Normal mode enumeration o Density of states o Debye model Review By considering
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 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 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 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 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 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 informationPhonons II: Thermal properties
Phonons II: Thermal properties specific heat of a crystal density of state Einstein mode Debye model anharmonic effect thermal conduction A technician holding a silica fibre thermal insulation tile at
More informationPhysics 211B : Problem Set #0
Physics 211B : Problem Set #0 These problems provide a cross section of the sort of exercises I would have assigned had I taught 211A. Please take a look at all the problems, and turn in problems 1, 4,
More informationPhonons (Classical theory)
Phonons (Classical theory) (Read Kittel ch. 4) Classical theory. Consider propagation of elastic waves in cubic crystal, along [00], [0], or [] directions. Entire plane vibrates in phase in these directions
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 informationBasics of Electronic-Structure Theory
A guideline through the tutorial 23.06.2009 For each exercise in this tutorial, a separate directory should be created in your home directory in which the input files are generated and the calculation
More informationPhonons I - Crystal Vibrations (Kittel Ch. 4)
Phonons I - Crystal Vibrations (Kittel Ch. 4) Displacements of Atoms Positions of atoms in their perfect lattice positions are given by: R 0 (n 1, n 2, n 3 ) = n 10 x + n 20 y + n 30 z For simplicity here
More information4. Thermal properties of solids. Time to study: 4 hours. Lecture Oscillations of the crystal lattice
4. Thermal properties of solids Time to study: 4 hours Objective After studying this chapter you will get acquainted with a description of oscillations of atoms learn how to express heat capacity for different
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 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 informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2009 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 Spring 2009 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 Spring 2009 POP QUIZ
More informationSoutheast University, Nanjing, China 2 Department of Applied Physics, Aalto University,
Supplementary Information to Solubility of Boron, Carbon and Nitrogen in Transition Metals: Getting Insight into Trends from First-Principles Calculations Xiaohui Hu, 1,2 Torbjörn Björkman 2,3, Harri Lipsanen
More information3.091 Introduction to Solid State Chemistry. Lecture Notes No. 5a ELASTIC BEHAVIOR OF SOLIDS
3.091 Introduction to Solid State Chemistry Lecture Notes No. 5a ELASTIC BEHAVIOR OF SOLIDS 1. INTRODUCTION Crystals are held together by interatomic or intermolecular bonds. The bonds can be covalent,
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechanics Rajdeep Sensarma sensarma@theory.tifr.res.in Quantum Dynamics Lecture #2 Recap of Last Class Schrodinger and Heisenberg Picture Time Evolution operator/ Propagator : Retarded
More information5.62 Physical Chemistry II Spring 2008
MI OpenCourseWare http://ocw.mit.edu 5.6 Physical Chemistry II Spring 008 For information about citing these materials or our erms of Use, visit: http://ocw.mit.edu/terms. 5.6 Spring 008 Lecture Summary
More informationMatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing
MatSci 331 Homework 4 Molecular Dynamics and Monte Carlo: Stress, heat capacity, quantum nuclear effects, and simulated annealing Due Thursday Feb. 21 at 5pm in Durand 110. Evan Reed In this homework,
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 informationTransition Theory Abbreviated Derivation [ A - B - C] # E o. Reaction Coordinate. [ ] # æ Æ
Transition Theory Abbreviated Derivation A + BC æ Æ AB + C [ A - B - C] # E A BC D E o AB, C Reaction Coordinate A + BC æ æ Æ æ A - B - C [ ] # æ Æ æ A - B + C The rate of reaction is the frequency of
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 informationNon-Continuum Energy Transfer: Phonons
Non-Continuum Energy Transfer: Phonons D. B. Go Slide 1 The Crystal Lattice The crystal lattice is the organization of atoms and/or molecules in a solid simple cubic body-centered cubic hexagonal a NaCl
More informationHands-On Tutorial on Ab Initio Molecular Simulations Berlin, July 12 21, Tutorial II: Periodic systems Manuscript for Exercise Problems
Hands-On Tutorial on Ab Initio Molecular Simulations Berlin, July 12 21, 2011 Tutorial II: Periodic systems Manuscript for Exercise Problems Prepared by Jürgen Wieferink, Lydia Nemec, and Volker Blum.
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 22 Dec 2004
High-accuracy first-principles determination of the structural, vibrational and thermodynamical properties of diamond, graphite, and derivatives arxiv:cond-mat/412643v1 [cond-mat.mtrl-sci] 22 Dec 24 Nicolas
More informationintroduction of thermal transport
Subgroup meeting 2010.12.07 introduction of thermal transport members: 王虹之. 盧孟珮 introduction of thermal transport Phonon effect Electron effect Lattice vibration phonon Debye model of lattice vibration
More informationThermodynamics Molecular Model of a Gas Molar Heat Capacities
Thermodynamics Molecular Model of a Gas Molar Heat Capacities Lana Sheridan De Anza College May 3, 2018 Last time modeling an ideal gas at the microscopic level rms speed of molecules equipartition of
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 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 informationThermodynamics of hexagonal-close-packed iron under Earth s core conditions
PHYSICAL REVIEW B, VOLUME 64, 045123 Thermodynamics of hexagonal-close-packed iron under Earth s core conditions D. Alfè, 1,2 G. D. Price, 1 and M. J. Gillan 2 1 Research School of Geological and Geophysical
More informationThe electronic structure of materials 1
Quantum mechanics 2 - Lecture 9 December 18, 2013 1 An overview 2 Literature Contents 1 An overview 2 Literature Electronic ground state Ground state cohesive energy equilibrium crystal structure phase
More informationPHYSICS 219 Homework 2 Due in class, Wednesday May 3. Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235.
PHYSICS 219 Homework 2 Due in class, Wednesday May 3 Note: Makeup lectures on Friday May 12 and 19, usual time. Location will be ISB 231 or 235. No lecture: May 8 (I m away at a meeting) and May 29 (holiday).
More informationPH575 Spring Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5
PH575 Spring 2014 Lecture #26 & 27 Phonons: Kittel Ch. 4 & 5 PH575 POP QUIZ Phonons are: A. Fermions B. Bosons C. Lattice vibrations D. Light/matter interactions PH575 POP QUIZ Phonon dispersion relation:
More informationSupplementary Online Materials: Formation of Stoichiometric CsF n Compounds
1 2 3 4 5 6 7 8 9 1 11 12 13 Supplementary Online Materials: Formation of Stoichiometric CsF n Compounds Qiang Zhu, 1, a) Artem R. Oganov, 1, 2, 3 and Qingfeng Zeng 4 1) Department of Geosciences, Stony
More informationLecture 9 Anharmonic effects in crystals.
Lecture 9 Anharmonic effects in crystals. 1 Introduction Most of the effects related to lattice dynamics that you have so far encountered in this course and in previous courses lattice specific heat, Debye-Waller
More informationRelative stability of α and β boron
Relative stability of α and β boron Michael Widom 1 and Marek Mihalkovič 2 1 Department of Physics, Carnegie Mellon University, Pittsburgh PA 15213, USA 2 Department of Physics, Slovakian Academy of Sciences,
More informationSummary: Thermodynamic Potentials and Conditions of Equilibrium
Summary: Thermodynamic Potentials and Conditions of Equilibrium Isolated system: E, V, {N} controlled Entropy, S(E,V{N}) = maximum Thermal contact: T, V, {N} controlled Helmholtz free energy, F(T,V,{N})
More informationPhonons II. Thermal Properties
Chapter 5. Phonons II. Thermal Properties Thermal properties of phonons As mentioned before, we are now going to look at how what we know about phonons will lead us to a description of the heat capacity
More informationFYS Vår 2015 (Kondenserte fasers fysikk)
FYS410 - Vår 015 (Kondenserte fasers fysikk) http://www.uio.no/studier/emner/matnat/fys/fys410/v15/index.html Pensum: Introduction to Solid State Physics by Charles Kittel (Chapters 1-9 and 17, 18, 0)
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 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 informationPhonons Thermal energy Heat capacity Einstein model Density of states Debye model Anharmonic effects Thermal expansion Thermal conduction by phonons
3b. Lattice Dynamics Phonons Thermal energy Heat capacity Einstein model Density of states Debye model Anharmonic effects Thermal expansion Thermal conduction by phonons Neutron scattering G. Bracco-Material
More informationSolutions for Homework 4
Solutions for Homework 4 October 6, 2006 1 Kittel 3.8 - Young s modulus and Poison ratio As shown in the figure stretching a cubic crystal in the x direction with a stress Xx causes a strain e xx = δl/l
More informationLab 1: Empirical Energy Methods Due: 2/14/18
Lab 1: Empirical Energy Methods Due: 2/14/18 General remarks on scientific scripting Scientific scripting for managing the input and output data is an important component of modern materials computations,
More informationThe Toric-Boson model and quantum memory at finite temperature
The Toric-Boson model and quantum memory at finite temperature A.H., C. Castelnovo, C. Chamon Phys. Rev. B 79, 245122 (2009) Overview Classical information can be stored for arbitrarily long times because
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 informationQuantum Condensed Matter Physics Lecture 5
Quantum Condensed Matter Physics Lecture 5 detector sample X-ray source monochromator David Ritchie http://www.sp.phy.cam.ac.uk/drp2/home QCMP Lent/Easter 2019 5.1 Quantum Condensed Matter Physics 1. Classical
More informationChapter 5: Thermal Properties of Crystal Lattices
Chapter 5: Thermal Properties of Crystal Lattices Debye January 30, 07 Contents Formalism. The Virial Theorem............................. The Phonon Density of States...................... 5 Models of
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 informationLecture 3: Optical Properties of Insulators, Semiconductors, and Metals. 5 nm
Metals Lecture 3: Optical Properties of Insulators, Semiconductors, and Metals 5 nm Course Info Next Week (Sept. 5 and 7) no classes First H/W is due Sept. 1 The Previous Lecture Origin frequency dependence
More informationLecture 10 Planck Distribution
Lecture 0 Planck Distribution We will now consider some nice applications using our canonical picture. Specifically, we will derive the so-called Planck Distribution and demonstrate that it describes two
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 informationLecture 6 Quantum Mechanical Systems and Measurements
Lecture 6 Quantum Mechanical Systems and Measurements Today s Program: 1. Simple Harmonic Oscillator (SHO). Principle of spectral decomposition. 3. Predicting the results of measurements, fourth postulate
More informationDFT calculation of pressure induced phase transition in silicon
DFT calculation of pressure induced phase transition in silicon Michael Scherbela scherbela@student.tugraz.at 2015-12-07 Contents 1 Introduction 2 2 Structure of the input file 2 3 Calculation of phase
More informationGeneration of Thermal Scattering Laws for YH 2 using Ab Initio Methods
Generation of Thermal Scattering Laws for YH 2 using Ab Initio Methods Michael L. Zerkle Bettis Atomic Power Laboratory WPEC SG42 Meeting May 18, 2015 May 18-19, 2015 WPEC SG42 Slide 1 Outline Motivation
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 informationAb initio Berechungen für Datenbanken
J Ab initio Berechungen für Datenbanken Jörg Neugebauer University of Paderborn Lehrstuhl Computational Materials Science Computational Materials Science Group CMS Group Scaling Problem in Modeling length
More informationThe Dulong-Petit (1819) rule for molar heat capacities of crystalline matter c v, predicts the constant value
I believe that nobody who has a reasonably reliable sense for the experimental test of a theory will be able to contemplate these results without becoming convinced of the mighty logical power of the quantum
More informationLecture 11 - Phonons II - Thermal Prop. Continued
Phonons II - hermal Properties - Continued (Kittel Ch. 5) Low High Outline Anharmonicity Crucial for hermal expansion other changes with pressure temperature Gruneisen Constant hermal Heat ransport Phonon
More information6.730 Physics for Solid State Applications
6.730 Physics for Solid State Applications Lecture 5: Specific Heat of Lattice Waves Outline Review Lecture 4 3-D Elastic Continuum 3-D Lattice Waves Lattice Density of Modes Specific Heat of Lattice Specific
More informationTopic 5-1: Introduction to Phonons Kittel pages: 91, 92
Topic 5-1: Introduction to Phonons Kittel pages: 91, 92 Summary: In this video we introduce the concept that atoms are not rigid, fixed points within the lattice. Instead we treat them as quantum harmonic
More informationC. Show your answer in part B agrees with your answer in part A in the limit that the constant c 0.
Problem #1 A. A projectile of mass m is shot vertically in the gravitational field. Its initial velocity is v o. Assuming there is no air resistance, how high does m go? B. Now assume the projectile is
More informationSupplementary Information
Supplementary Information Ballistic Thermal Transport in Carbyne and Cumulene with Micron-Scale Spectral Acoustic Phonon Mean Free Path Mingchao Wang and Shangchao Lin * Department of Mechanical Engineering,
More informationVibrational Spectroscopy
Vibrational Spectroscopy Keith Refson STFC Rutherford Appleton Laboratory August 28, 2009 Density Functional Methods for Experimental Spectroscopy 2009: Oxford 1 / 22 Two similar structures Zincblende
More informationPhysics 607 Final Exam
Physics 607 Final Exam Please be well-organized, and show all significant steps clearly in all problems. You are graded on your work, so please do not just write down answers with no explanation! Do all
More information9.3. Total number of phonon modes, total energy and heat capacity
Phys50.nb 6 E = n = n = exp - (9.9) 9... History of the Planck distribution or the Bose-Einstein distribution. his distribution was firstly discovered by Planck in the study of black-body radiation. here,
More informationRoad map (Where are we headed?)
Road map (Where are we headed?) oal: Fairly high level understanding of carrier transport and optical transitions in semiconductors Necessary Ingredients Crystal Structure Lattice Vibrations Free Electron
More informationBasics of DFT applications to solids and surfaces
Basics of DFT applications to solids and surfaces Peter Kratzer Physics Department, University Duisburg-Essen, Duisburg, Germany E-mail: Peter.Kratzer@uni-duisburg-essen.de Periodicity in real space and
More informationOutline. Introduction: graphene. Adsorption on graphene: - Chemisorption - Physisorption. Summary
Outline Introduction: graphene Adsorption on graphene: - Chemisorption - Physisorption Summary 1 Electronic band structure: Electronic properties K Γ M v F = 10 6 ms -1 = c/300 massless Dirac particles!
More informationSupplementary Information
Supplementary Information Supplementary Figure 1: After structural optimization of the CH 3 NH 3 PbI 3 unit cell, the resulting relaxed volumes for three distinct orientation of the MA molecules are shown.
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department Statistical Physics I Spring Term 2013
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Physics Department 8.044 Statistical Physics I Spring Term 2013 Problem 1: The Big Bang Problem Set #9 Due in hand-in box by 4;00 PM, Friday, April 19 Early in the
More informationPHYS3113, 3d year Statistical Mechanics Tutorial problems. Tutorial 1, Microcanonical, Canonical and Grand Canonical Distributions
1 PHYS3113, 3d year Statistical Mechanics Tutorial problems Tutorial 1, Microcanonical, Canonical and Grand Canonical Distributions Problem 1 The macrostate probability in an ensemble of N spins 1/2 is
More informationMolecular energy levels and spectroscopy
Molecular energy levels and spectroscopy 1. Translational energy levels The translational energy levels of a molecule are usually taken to be those of a particle in a three-dimensional box: n x E(n x,n
More informationQuasiharmonic models for the calculation of thermodynamic properties of crystalline silicon under strain
JOURNAL OF APPLIED PHYSICS 99, 64314 26 Quasiharmonic models for the calculation of thermodynamic properties of crystalline silicon under strain H. Zhao and Z. Tang Department of Mechanical and Industrial
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
More informationCrystals. Peter Košovan. Dept. of Physical and Macromolecular Chemistry
Crystals Peter Košovan peter.kosovan@natur.cuni.cz Dept. of Physical and Macromolecular Chemistry Lecture 1, Statistical Thermodynamics, MC26P15, 5.1.216 If you find a mistake, kindly report it to the
More information1+e θvib/t +e 2θvib/T +
7Mar218 Chemistry 21b Spectroscopy & Statistical Thermodynamics Lecture # 26 Vibrational Partition Functions of Diatomic Polyatomic Molecules Our starting point is again the approximation that we can treat
More information10.40 Lectures 23 and 24 Computation of the properties of ideal gases
1040 Lectures 3 and 4 Computation of the properties of ideal gases Bernhardt L rout October 16 003 (In preparation for Lectures 3 and 4 also read &M 1015-1017) Degrees of freedom Outline Computation of
More informationIntroduction to Nonlinear Optics
Introduction to Nonlinear Optics Prof. Cleber R. Mendonca http://www.fotonica.ifsc.usp.br Outline Linear optics Introduction to nonlinear optics Second order nonlinearities Third order nonlinearities Two-photon
More information