The Nuts and Bolts of FirstPrinciples Simulation 23
The Nuts and Bolts of First-Principles Simulation 23: Calculating Phonon Dispersion Relations and other matters! Durham, 6 th-13 th December 2001 Nuts and Bolts 2001 CASTEP Developers’ Group with support from the ESF k Network 23: Calculations using phonons
Outline Introduction phonon dispersion calculations. r Some examples for a range of materials. r Some surfaces. r Density of states. r Free energies. r Phase transitions. r Nuts and Bolts 2001 23: Calculations using phonons� 2
Introduction We will use a linear response method to calculate phonon properties of a crystalline material. r This avoids the need to create a supercell as in the frozen phonon method. r LR theory allows us to calculate phonon properties for any q-vector. Arbitrary q-vector leads to very large supercells (in principle) using frozen phonons. r Nuts and Bolts 2001 23: Calculations using phonons� 3
Examples of LRT Phonons r r r First we will look at Ga. As. The dispersion curve has been calculated by repeated calculation of w(q) for a range of q-vectors. As with band structures, use primitive cell. Nuts and Bolts 2001 23: Calculations using phonons� 4
Dispersion Relation for Ga. As phonon dispersion relation along some lines of high symmetry. r Dots indication neutron-scattering data. r Nuts and Bolts 2001 23: Calculations using phonons� 5
Phonons in Carbon-diamond r r For many nearly all group-V and III-V materials, the maximum phonon frequency occurs at the G point. Anomaly in diamond – G point frequency is a local minimum. Experiment: Inelastic x-ray and neutron scattering data. Nuts and Bolts 2001 23: Calculations using phonons� 6
Anomaly in Diamond Examine q-vectors in all directions around G and we see that there is a local minimum. Shown here in a (100) plane. A similar results is found is other directions. Nuts and Bolts 2001 23: Calculations using phonons� 7
Explanation for Diamond The dispersion of the uppermost branch has a strong over-binding in any direction. r Examination of the phonon eigenvectors indicate direction of atom motion for each q-vector – examine binding for yourself! r This results in a minimum frequency at the zone centre. r This is at variance with other tetrahedral semiconductors. r Nuts and Bolts 2001 23: Calculations using phonons� 8
Wurtzite Structures r r Wide band gap wurtzite semiconductors (eg. Ga. N) are of interest for blue/UV LED’s. Behaviour of carriers (electrons or holes) are influenced by interaction with phonons. Hence lattice-dynamical properties are very important input into various phenomenological models. Experimentally difficult to measure phonon dispersion since decent crystals are hard to grow, hence experimental information is limited. Nuts and Bolts 2001 23: Calculations using phonons� 9
Ga. N – The Wurtzite Structure Semiconducting Ga. N in the hexagonal Wurtzite structure. Nuts and Bolts 2001 The Brillouin Zone constructs and a path of high symmetry. 23: Calculations using phonons� 10
Dispersion Relation in Wurtzite Structure Agreement with experiment is reasonably good. Deviations from experiment exist (eg GM line). What’s better: experiment of theory? Nuts and Bolts 2001 23: Calculations using phonons� 11
Density of States Calculating the density of states for phonons is performed in a similar manner to that of the electronic density of states. r Sample q-vectors over the Brillouin zone (accounting for symmetry) and integrate up. r Interpolation of frequencies between the calculated values at the given q-vectors can be used to increase integration accuracy. r Nuts and Bolts 2001 23: Calculations using phonons� 12
Dispersion and DOS of Si Note: Singularities in DOS is difficult to achieve, but important to get right for free energy calculations! Nuts and Bolts 2001 23: Calculations using phonons� 13
Free Energy r r An important application of ab initio phonon calculations is the ability to evaluate free energies. This is important in the study of phase diagrams. Whichever phase has the lowest free energy will (should!) be the one that’s observed. Free energy of a temperature (T)-volume (V) system is given as follows: Nuts and Bolts 2001 23: Calculations using phonons� 14
Example of Phase Transitions r r r Tin is commonly found in one of two allotropic forms at ambient pressure. The low temperature phase (a-Sn) is a zero gap semiconductor in the diamond structure. Above 13 o. C the crystal transforms into the b phase (white tin which is metallic): body-centred tetragonal. Nuts and Bolts 2001 23: Calculations using phonons� 15
a-b Sn Transition The a-b transition is a simple example of an entropy driven structural phase transition. r It can be examined using the vibrational properties of the two phases of the material. r Examining the free energy equation, each term is easily accessible to the LR-DFPT calculations presented in the previous lecture. r Nuts and Bolts 2001 23: Calculations using phonons� 16
Sn Dispersion Curves Nuts and Bolts 2001 23: Calculations using phonons� 17
Phase Diagram for Sn Theoretical phase transition temperature is slightly above the experimental value. The missing factor is anharmonicity. Nuts and Bolts 2001 23: Calculations using phonons� 18
Surface Phonons r As with the electronic density of states of a surface, we can perform a similar calculation for phonons: This is the surface phonon dispersion curve for In. P (100) surface. Solid lines indicate the phonons that are localised at the surface. Shaded area are bulk phonons projected onto surface BZ. Nuts and Bolts 2001 23: Calculations using phonons� 19
Elastic Constants r r r A full derivation of elastic constants from the Dynamics matrix can be found in A. A. Maradudin, et. al. “Theory of lattice dynamics in the harmonic approximation”, Academic Press, 1971. I will skip a few steps, but basic results are as follows: Allow a, b, g, l, m, n to run over the co-ordinate axes and define: Nuts and Bolts 2001 23: Calculations using phonons� 20
Elastic Constants r If V is the volume of the unit cell, then we can define the following brackets: The elastic constants are then: Nuts and Bolts 2001 23: Calculations using phonons� 21
Definition of Elastic Constants r The elastic constants can then be expressed in their more familiar form by pairing indices by the equivalences: This leads to the standard form, eg. c 11, c 12, c 44, etc These can be used to calculate various modulus, e. g. for a cubic crystal, the bulk modulus is: Nuts and Bolts 2001 23: Calculations using phonons� 22
Structure of Mg. O r r Mg. O (along with Fe. O) is hypothesised to be a major constituent of the lower mantle. Experiments at lower mantle pressures are not possible. Nuts and Bolts 2001 23: Calculations using phonons� 23
Elastic Constants of Mg. O r Elastic constants of Minerals at high pressure is of physical and geological interest, e. g. for interpretation of seismological data Solid Line: Theory. Dashed Lines: Interpolated experimental data using two models. Nuts and Bolts 2001 23: Calculations using phonons� 24
Bulk/Shear Modulus for Mg. O The shear modulus, G, is calculated as follows: Nuts and Bolts 2001 23: Calculations using phonons� 25
Energy Propagation in Mg. O The longitudinal and shearwave velocities, VP and VS of an isotropic aggregate can be calculated using: Nuts and Bolts 2001 23: Calculations using phonons� 26
Propagation Direction in Mg. O Dependence of the longitudinal and two shearwave velocities of Mg. O on propagation direction at pressures of 0 GPa (solid) and 150 GPa (dashed). Nuts and Bolts 2001 23: Calculations using phonons� 27
References P. Pavone, J, . Phys. Condens. Matter 13, 7593 (2001). r B. B. Karki, S. J. Clark, et. al. , American Mineralogist 82, 52 (1997). r J. Frirsch, et. al. , Surface Science 427 -428, 58 (1999). r Nuts and Bolts 2001 23: Calculations using phonons� 28
- Slides: 28