Fonones y Elasticidad bajo presión ab initio .

1. Fonones y Elasticidad bajo presión ab initio . Alfonso Muñoz Dpto. de Física Fundamental II Universidad de La Laguna. Tenerife. MALTA Consolider Team Canary Islands, SPAIN

2. Plan de la charla: • Introducción : Ab initio methods • Fonones. Propiedades dinámicas. Ejemplos • Elasticidad estabilidad mecánica bajo presión • Conclusiones.

3. Ab initio methods • State of the art Ab Initio Total Energy Pseudopotential calculations are useful to study many properties of materials. • No experimental input required (even the structure). Only Z is required. They can provide and “predict” many properties of the material if the approximations are correct! • DFT is the standar theory applied, it is “exact” but one need to use approximations, XC functional (LDA, GGA etc…), BZ integration with k-special points, etc. (some problems in high correlated systems, f-electrons etc..). DFPT also available, allows to study phonons, elastic constants etc… • More elaborated approximations are also available, like LDA + U, MD, etc.. • Many computer programs available, some times free (Abinit, quantum espresso, VASP, CASTEP, etc…) • Ab initio methods provide and alternative and complimentary technique to the experiments under extreme conditions.

4. "Those who are enamoured of Practice without Theory are like a pilot who goes into a ship without rudder or compass and never has any certainty of where he is going. Practice should always be based upon a sound knowledge of Theory.“ Leonardo da Vinci, (1452-1519)

6. Well tested:

7. . “Prediction is very difficult, especially about the future”. Niels Bohr (1885-1962)

8. Fonones y espectroscopía, ¿paraque? Thermal Expansion Superconductivity Thermal Conductivity Elasticity - deformation

9. Hooke’s law! IFC Lattice Dynamics • Lattice Potential: • Harmonic approximation: ISVS: A hands-on introduction to ABINIT

10. Phonons: linear chain of atoms Harmonic approximation: Ansatz: => ISVS: A hands-on introduction to ABINIT

11. Phonons: linear chain of atoms Linear chain of atoms

12. Linear chain with two different "spring constants" two atoms per unit cell Ansatz: ISVS: A hands-on introduction to ABINIT

13. Phonons Linear chain with two different "spring constants" Two solutions: acoustic (-) and optic (+) branches ISVS: A hands-on introduction to ABINIT

14. THz cm-1 Si 3C-SiC LO TO G. Nilsson and G. Nelin, PRB 6, 3777 (1972) W. Weber, PRB 15, 4789 (1977) J. Serrano et al., APL 80, 23 (2002) 3D Phonon Dispersion Relations 3THz ~ 100 cm-1 ; 1meV ~ 8 cm-1 j = 3N branches

15. Polar crystals: LST relation • Ionic crystals: Macroscopic electric field Lyddane-Sachs-Teller relation Born effective charges X. Gonze and C. Lee, PRB 55, 10355 (1997)

16. Anisotropy: crystal field GaN T. Ruf et al., PRL 86, 906 (2001)

17. Anisotropy: Selection Rules • Not all modes are visible with the same technique! B1: SILENT modes • Not all allowed modes are visible at the same time! J.M. Zhang et al., PRB 56, 14399 (1997)

18. Christoffel equations Elasticity h-BN A. Bosak et al., PRB 73 041402(R) (2006)

19. Vibrational spectroscopies • Probes: • Light: photons • Particles • Brillouin spectroscopy • Raman spectroscopy • Infrared absorption spect. • Inelastic X-ray Scattering • electrons:High Resolution e- Energy Loss • He:He atom scattering • neutrons • Time-of-flight spectroscopy • Inelastic Neutron Scattering ISVS: A hands-on introduction to ABINIT

20. Vibrational spectroscopies

21. X-ray scattering • Excitations ~ meV, ~ Å-1 • whole BZ available • No kinematics restrictions • Dispersion + r() • Energy resolution ~1meV • Raman spec. • 1meV-eV Excitations • Optic phonons at the center of • the Brillouin zone • High resolution • Different selection rules • Neutron scattering • Excitations ~ meV, ~ Å-1 • whole BZ available • Dispersion + r(w) Kinematical limit: vs < 3000 m/s Vibrational spectroscopies Optic b. Acoustic branches • Brillouin spec. • Excitations of 2meV-0.6meV • Acoustic phonon branches at low q • (sound waves) • Information: E Vs (sound speed) • linewidth Atenuation

22. Infrared Spectroscopy • Absorption spectroscopy: dipolar selection rules • Target: polar molecular vibrations, determination of functional groups in organic compounds, polar modes in crystals ISVS: A hands-on introduction to ABINIT

23. BILBAO CRYSTAL… SERVER SAM

24. ISOTROPY PACKAGE (STOKES ET AL.) SMODES, FINDSYM, ETC……..

25. Interatomic force constants Born-Oppenheimer approximation In the harmonic approximation, the total energy of a crystal with small atomic position deviations is: where the matrix of IFC’s is defined as:

26. Physical Interpretation of the Interatomic Force Constants The force conjugate to the position of a nucleus, can always be written: We can thus rewrite the IFC’s in more physically descriptive fashion: The IFC’s are the rate of change of the atomic forces when we displace another atom in the crystal.

27. Dynamical properties under pressure.. The construction of the dynamical matrix at gamma point is very simple: Phonon dispersion, DOS, PDOS requires supercell calculations. Also DFPT allows to include T effects, Thermod. properties, etc…

28. Relation between the IFC’s and the dynamical matrix The Fourier transform of the IFC’s is directly related to the dynamical matrix, The phonon frequencies are then obtained by diagonalization of the dynamical matrix or equivalently by the solution of this eigenvalue problem: masses phonon displacementpattern square ofphonon frequencies

29. Phonon band structure of α-Quartz SiO2 9 atoms per unit cell Nb. of phonon bands: Nb. of acoustic bands: 3 Nb. of optical bands: Polar crystal : LO non-analyticity Directionality ! [X.Gonze, J.-C.Charlier, D.C.Allan, M.P.Teter, PRB50, 13055 (1994)]

30. High - temperature : Fluorite structure ( , one formula unit per cell )  Supercell calculation + interpolation ! Long-range dipole-dipoleinteraction not taken into account LO-TO splitting Wrong behaviour Calculated phonon dispersions of ZrO2 in the cubic structure at the equilibrium lattice constant a0 = 5.13 Å. [From Parlinski K., Li Z.Q., and Kawazoe Y.,Phys. Rev. Lett. 78, 4063 (1997)] • DFPT (Linear-response)with = 5.75 = -2.86 and = 5.75 • LO - TO splitting 11.99 THz • Non-polar mode is OK

31. Thermodynamic properties In the harmonic approximation, the phonons can be treated as an independent boson gas. They obey the Bose-Einstein distribution: The total energy of the gas can be calculated directly using the standard formula: Phonon DOS Energy of the harmonic oscillator All thermodynamic properties can be calculated in this manner. Note: ISVS 2008: Phonon Bands and Thermodynamic Properties

32. Evenwith f electrons (PRB 85, 024317 (2012) TbPO4 , DyPO4

33. ZnS, Phys. Rev B 81 075207 (2010) Cardona, …Muñoz . et al.

35. Spin-orbit, phonon dispersión, temperature effects, etc…. DFPT Inverted s-o interaction. Contribution of the negative splitting of 5d G15 states of Hg wich overcompesate the positive splitting of the S 3p.

36. CuGaS2 electronic and phononic properties Eficient photovoltaic materials. (Phys.Rev. B 83,195208 (2011) ) Chalcopyrite tetragonal SG I-42d Few it is know about this compounds. We did a structural , electronic and phononic study of the thermodynamical properties • Two main groups of chalcopyrites: • I-III-VI2 derived from II-VI zb compounds (CuGaS2, AgGaS2..) • II-IV-V2 derived from III-V zb comp.ounds (ZnGaAs2,….) • Two formula units per primitive cell • We will focus on the study of some • thermodynamics properties, like the • specific heat with emphasis in the low-T • region where appear strong desviation • of the Debye T3 law, phonons, etc…

37. silicon zincblende, ZnSS chalcopyrite, CuGaS2

39. Elastic Constants Cij (no experimental data available)

40. 5 Phonons Starting from the electronic structure we calculate the phonon dispersion relations with density functional perturbation theory. We compare them with Raman and IR measurements at the center of the BZ (see Figure). comparison with Raman and IR measurements () shows good agreement inelastic neutron scattering data are not available as yet CuGaS2

41. Phonon Density of States Through BZ integration of the phonon dispersion relations the phonon density of states (total or projected on the individual atoms Cu, Ga, S)) are obtained (see Figure). 7 PDOS (states / formula unit) sulphur-like below 120 cm-1: essentially Cu- and Ga-like phonons above 280 cm-1: essentially S-like phonons midgap feature at ~180 cm-1: Ga-, Cu-like Cu-Ga Ga-Cu The partial density of states are useful for calculating the effect of isotope disorder on the phonon linewidths

42. Two-phonon No second-order Raman spectra available. The calculated sum an difference densities will help to interpret future measured spectra. It is posible to establish a correspondence between the calculated two-phonon Raman spectra of CuGaS2 and other two-phonon measured spectra of binary compounds.

43. The effect of phonons the on Vo(T) for a (cubic) crystal can be expressed in terms of mode Grüneisen parameters γqj : Due to the large number of phonons bands, a first approximation is to use only the values at the Zone center for the evaluation of the termal expansion coefficient. The temperature dependence of Vo for q= 0 is: Or from thermodyn… using S(P,T)

46. Heat capacity The phonon DOS allows to calculate the Free Energy F(T), and the specific heat at constant volume And the constant pressure Cp can be obtained:

48. 8 qDebye(0) = 355 K Comparison of calculated and measured specific heat CuGaS2 versus AgGaS2 peak at ~ 20 K in CP/T 3 representation from Cu/Ga like phonons (ratio 1:6 to low-frequency peak in phonon DOS) ABINIT LDA reproduces peak position, but absolute value at peak ~20% lower VASP GGA reproduces peak position and magnitude

49. 9 Comparison of calculated and measured specific heat Extension to AgGaS2 and AgGaTe2 lattice softening by Cu  Ag replacement lattice softening by S  Te replacement

50. Even more properties? Inclusion of Temperature effects is computationally very expensive, e-ph interaction,… Many experimental results of T dependence of the gap in binary and ternary compounds. The degree of cation-anion hybridization on the electronic an vibrational properties, leads to anomalous dependence of the band gaps with temperature. The presence of d-electrons in upper VB lead to anomalies, like negative s-o splitting. For example in Cu or Ag chalcopyrite, the other constituents correspond to decrease the gap, but Cu or Ag tends to increase. The sum of boths effects generates a non monotonic dependence of gaps with T. It can be fitted using two Einstein oscillator according to: E0 is the zero-point un-renormalized gap energy, A1 is the contribution to the zero-point renormal., nB is the Bose-Eisntein function .