Phonon dispersion curves and atomic mean square displacement for several fcc and bcc materials

The atomic mean square displacement (MSD) and the phonon dispersion curves (PDC's) of a number of face-centred cubic (fcc) and body-centred cubic (bcc) materials have been calclllated from the quasiharmonic (QH) theory, the lowest order (A2 ) perturbation theory (PT) and a recently proposed Green's function (GF) method by Shukla and Hiibschle. The latter method includes certain anharmonic effects to all orders of anharmonicity. In order to determine the effect of the range of the interatomic interaction upon the anharmonic contributions to the MSD we have carried out our calculations for a Lennard-Jones (L-J) solid in the nearest-neighbour (NN) and next-nearest neighbour (NNN) approximations. These results can be presented in dimensionless units but if the NN and NNN results are to be compared with each other they must be converted to that of a real solid. When this is done for Xe, the QH MSD for the NN and NNN approximations are found to differ from each other by about 2%. For the A2 and GF results this difference amounts to 8% and 7% respectively. For the NN case we have also compared our PT results, which have been calculated exactly, with PT results calculated using a frequency-shift approximation. We conclude that this frequency-shift approximation is a poor approximation. We have calculated the MSD of five alkali metals, five bcc transition metals and seven fcc transition metals. The model potentials we have used include the Morse, modified Morse, and Rydberg potentials. In general the results obtained from the Green's function method are in the best agreement with experiment. However, this improvement is mostly qualitative and the values of MSD calculated from the Green's function method are not in much better agreement with the experimental data than those calculated from the QH theory. We have calculated the phonon dispersion curves (PDC's) of Na and Cu, using the 4 parameter modified Morse potential. In the case of Na, our results for the PDC's are in poor agreement with experiment. In the case of eu, the agreement between the tlleory and experiment is much better and in addition the results for the PDC's calclliated from the GF method are in better agreement with experiment that those obtained from the QH theory.

Phonon modes of MgB2 : Super-lattice structures and spectral response

Micrometre-sized MgB2 crystals of varying quality, synthesized at low temperature and autogeneous pressure, are compared using a combination of Raman and Infra-Red (IR) spectroscopy. These data, which include new peak positions in both spectroscopies for high quality MgB2, are interpreted using DFT calculations on phonon behaviour for symmetry-related structures. Raman and IR activity additional to that predicted by point group analyses of the P6/mmm symmetry are detected. These additional peaks, as well as the overall shapes of calculated phonon dispersion (PD) models are explained by assuming a double super-lattice, consistent with a lower symmetry structure for MgB2. A 2x super-lattice in the c-direction allows a simple correlation of the pair breaking energy and the superconducting gap by activation of corresponding acoustic frequencies. A consistent physical interpretation of these spectra is obtained when the position of a phonon anomaly defines a super-lattice modulation in the a-b plane.

Coherent phonon decay and the boron isotope effect for MgB2

Ab-initio DFT calculations for the phonon dispersion (PD) and the Phonon Density Of States (PDOS) of the two isotopic forms (10B and 11B) of MgB2 demonstrate that use of a reduced symmetry super-lattice provides an improved approximation to the dynamical, phonon-distorted P6/mmm crystal structure. Construction of phonon frequency plots using calculated values for these isotopic forms gives linear trends with integer multiples of a base frequency that change in slope in a manner consistent with the isotope effect (IE). Spectral parameters inferred from this method are similar to that determined experimentally for the pure isotopic forms of MgB2. Comparison with AlB2 demonstrates that a coherent phonon decay down to acoustic modes is not possible for this metal. Coherent acoustic phonon decay may be an important contributor to superconductivity for MgB2.

A Revisit to Capture the Entire Behavior of Ultrasonic Wave Dispersion Characteristics of Single Walled Carbon Nanotubes Based on Nonlocal Elasticity Theory and Flugge Shell Model

In this paper, an ultrasonic wave propagation analysis in single-walled carbon nanotube (SWCNT) is re-studied using nonlocal elasticity theory, to capture the whole behaviour. The SWCNT is modeled using Flugge's shell theory, with the wall having axial, circumferential and radial degrees of freedom and also including small scale effects. Nonlocal governing equations for this system are derived and wave propagation analysis is also carried out. The revisited nonlocal elasticity calculation shows that the wavenumber tends to infinite at certain frequencies and the corresponding wave velocity tends to zero at those frequencies indicating localization and stationary behavior. This frequency is termed as escape frequency. This behavior is observed only for axial and radial waves in SWCNT. It has been shown that the circumferential waves will propagate dispersively at higher frequencies in nonlocality. The magnitudes of wave velocities of circumferential waves are smaller in nonlocal elasticity as compared to local elasticity. We also show that the explicit expressions of cut-off frequency depend on the nonlocal scaling parameter and the axial wavenumber. The effect of axial wavenumber on the ultrasonic wave behavior in SWCNTs is also discussed. The present results are compared with the corresponding results (for first mode) obtained from ab initio and 3-D elastodynamic continuum models. The acoustic phonon dispersion relation predicted by the present model is in good agreement with that obtained from literature. The results are new and can provide useful guidance for the study and design of the next generation of nanodevices that make use of the wave propagation properties of single-walled carbon nanotubes.

Phonon scattering limited performance of monolayer MoS2 and WSe2 n-MOSFET

In this paper we show the effect of electron-phonon scattering on the performance of monolayer (1L) MoS2 and WSe2 channel based n-MOSFETs. Electronic properties of the channel materials are evaluated using the local density approximation (LDA) in density functional theory (DFT). For phonon dispersion we employ the small displacement / frozen phonon calculations in DFT. Thereafter using the non-equilibrium Green's function (NEGF) formalism, we study the effect of electron-phonon scattering and the contribution of various phonon modes on the performance of such devices. It is found that the performance of the WSe2 device is less impacted by phonon scattering, showing a ballisticity of 83% for 1L-WSe2 FET for channel length of 10 nm. Though 1L-MoS2 FET of similar dimension shows a lesser ballisticity of 75%. Also in the presence of scattering there exist a a 21-36% increase in the intrinsic delay time (tau) and a 10-18% reduction in peak transconductance (g(m)) for WSe2 and MoS2 devices respectively. (C) 2015 Author(s).

Studies of phonon anharmonicity in solids

Today our understanding of the vibrational thermodynamics of materials at low temperatures is emerging nicely, based on the harmonic model in which phonons are independent. At high temperatures, however, this understanding must accommodate how phonons interact with other phonons or with other excitations. We shall see that the phonon-phonon interactions give rise to interesting coupling problems, and essentially modify the equilibrium and non-equilibrium properties of materials, e.g., thermodynamic stability, heat capacity, optical properties and thermal transport of materials. Despite its great importance, to date the anharmonic lattice dynamics is poorly understood and most studies on lattice dynamics still rely on the harmonic or quasiharmonic models. There have been very few studies on the pure phonon anharmonicity and phonon-phonon interactions. The work presented in this thesis is devoted to the development of experimental and computational methods on this subject.

Modern inelastic scattering techniques with neutrons or photons are ideal for sorting out the anharmonic contribution. Analysis of the experimental data can generate vibrational spectra of the materials, i.e., their phonon densities of states or phonon dispersion relations. We obtained high quality data from laser Raman spectrometer, Fourier transform infrared spectrometer and inelastic neutron spectrometer. With accurate phonon spectra data, we obtained the energy shifts and lifetime broadenings of the interacting phonons, and the vibrational entropies of different materials. The understanding of them then relies on the development of the fundamental theories and the computational methods.

We developed an efficient post-processor for analyzing the anharmonic vibrations from the molecular dynamics (MD) calculations. Currently, most first principles methods are not capable of dealing with strong anharmonicity, because the interactions of phonons are ignored at finite temperatures. Our method adopts the Fourier transformed velocity autocorrelation method to handle the big data of time-dependent atomic velocities from MD calculations, and efficiently reconstructs the phonon DOS and phonon dispersion relations. Our calculations can reproduce the phonon frequency shifts and lifetime broadenings very well at various temperatures.

To understand non-harmonic interactions in a microscopic way, we have developed a numerical fitting method to analyze the decay channels of phonon-phonon interactions. Based on the quantum perturbation theory of many-body interactions, this method is used to calculate the three-phonon and four-phonon kinematics subject to the conservation of energy and momentum, taking into account the weight of phonon couplings. We can assess the strengths of phonon-phonon interactions of different channels and anharmonic orders with the calculated two-phonon DOS. This method, with high computational efficiency, is a promising direction to advance our understandings of non-harmonic lattice dynamics and thermal transport properties.

These experimental techniques and theoretical methods have been successfully performed in the study of anharmonic behaviors of metal oxides, including rutile and cuprite stuctures, and will be discussed in detail in Chapters 4 to 6. For example, for rutile titanium dioxide (TiO₂), we found that the anomalous anharmonic behavior of the B_1g mode can be explained by the volume effects on quasiharmonic force constants, and by the explicit cubic and quartic anharmonicity. For rutile tin dioxide (SnO₂), the broadening of the B_2g mode with temperature showed an unusual concave downwards curvature. This curvature was caused by a change with temperature in the number of down-conversion decay channels, originating with the wide band gap in the phonon dispersions. For silver oxide (Ag₂O), strong anharmonic effects were found for both phonons and for the negative thermal expansion.

Optical-phonon modes and Frohlich potential in quantum dots

By extending the microscopic dipole model on optical-phonon modes as applied in quantum wells and quantum wires, to rectangular quantum dots (QD), optical phonon modes and their accompanying Frohlich potentials in QD are calculated and classified. When the bulk phonon dispersion is ignored, the optical phonon modes in QD can be clearly divided into the confined LO- and TO-bulk-like modes and the extended interface-like modes. Among the interface-like modes, a special attention is given to the corner modes, whose anisotropic behavior is depicted in the long wavelength limit. Based on the numerical results, a set of analytical formula are proposed to approximately describe the bulk-like modes, for which both the optical displacements and Frohlich potentials vanish at the interfaces. (C) 2000 Elsevier Science Ltd. All rights reserved.

Comparison of cooling rates for hot carriers in GaAs AlAs quantum wells based on macroscopic and microscopic phonon models

By using three analytical phonon models in quantum wells-the slab model, the guided-mode model, and the improved version of the Huang-Zhu model [Phys. Rev. B 38, 13 377 (1998)], -and the phonon modes in bulk, the energy-loss rates of hot carriers due to the Frohlich potential scattering in GaAs/AlAs multiple quantum wells (MQW's) are calculated and compared to those obtained based on a microscopic dipole superlattice model. In the study, a special emphasis is put on the effects of the phonon models on the hot-carrier relaxation process when taking the hot-phonon effect into account. Our numerical results show that, the calculated energy-loss rates based on the slab model and on the improved Huang-Zhu model are almost the same when ignoring the hot-phonon effect; however, with the hot phonon effect considered, the calculated cooling rate as well as the hot phonon occupation number do depend upon the phonon models to be adopted. Out of the four analytical phonon models investigated, the improved Huang-Zhu model gives the results most close to the microscopic calculation, while the guided-mode model presents the poorest results. For hot electrons with a sheet density around 10(12)/cm(2), the slab model has been found to overestimate the hot-phonon effect by more than 40% compared to the Huang-Zhu model, and about 75% compared to the microscopic calculation in which the phonon dispersion is fully included. Our calculation also indicates that Nash's improved version [J. Lumin. 44, 315 (1989)] is necessary for evaluating the energy-loss rates in quantum wells of wider well width, because Huang-Zhu's original analytical formulas an only approximately orthogonal for optical phonons associated with small in-plane wave numbers. [S0163-1829(99)08919-5].

OPTICAL-PHONON MODES AND FROHLICH POTENTIAL IN ONE-DIMENSIONAL QUANTUM-WELL WIRES

By extending our microscopic model on optical-phonon modes in quantum wells to one-dimensional (1D) quantum-well wires (QWW), the optical displacements and associated electrostatic potentials of optical-phonon modes in 1D QWW are calculated. The modes can be clearly divided into confined LO bulklike, TO bulklike modes, and extended interfacelike modes provided the bulk phonon dispersion is ignored. The character of each type of mode is illustrated with special attention to the interfacelike modes, which are hybrids of longitudinal- and transverse-optical waves from the corresponding bulk materials. Based on the numerical results, approximate analytical formulas for bulklike modes are presented. As in 2D wells, both the optical displacements and Frohlich potentials for the bulklike modes vanish at the interfaces. The finite dispersion of bulk phonons has a more pronounced effect on the 1D phonon modes because interfacelike modes show mixed characteristics of 2D interface and bulklike modes.

Phonon spectra and thermodynamic properties of rare gas solids based on empirical and semi-empirical (ab initio) two-body potentials : a comparative study /

We study the phonon dispersion, cohesive and thermal properties of raxe gas solids Ne, Ar, Kr, and Xe, using a variety of potentials obtained from different approaches; such as, fitting to crystal properties, purely ab initio calculations for molecules and dimers or ab initio calculations for solid crystalline phase, a combination of ab initio calculations and fitting to either gas phase data or sohd state properties. We explore whether potentials derived with a certain approaxih have any obvious benefit over the others in reproducing the solid state properties. In particular, we study phonon dispersion, isothermal ajid adiabatic bulk moduli, thermal expansion, and elastic (shear) constants as a function of temperatiue. Anharmonic effects on thermal expansion, specific heat, and bulk moduli have been studied using A^ perturbation theory in the high temperature limit using the neaxest-neighbor central force (nncf) model as developed by Shukla and MacDonald [4]. In our study, we find that potentials based on fitting to the crystal properties have some advantage, particularly for Kr and Xe, in terms of reproducing the thermodynamic properties over an extended range of temperatiures, but agreement with the phonon frequencies with the measured values is not guaranteed. For the lighter element Ne, the LJ potential which is based on fitting to the gas phase data produces best results for the thermodynamic properties; however, the Eggenberger potential for Ne, where the potential is based on combining ab initio quantum chemical calculations and molecular dynamics simulations, produces results that have better agreement with the measured dispersion, and elastic (shear) values. For At, the Morse-type potential, which is based on M0ller-Plesset perturbation theory to fourth order (MP4) ab initio calculations, yields the best results for the thermodynamic properties, elastic (shear) constants, and the phonon dispersion curves.

Phonon spectra and thermal properties of some fcc metals using embedded-atom potentials /

By employing the embedded-atom potentials of Mei et ai.[l], we have calculated the dynamical matrices and phonon dispersion curves for six fee metals (Cu,Ag,Au,Ni,Pd and Pt). We have also investigated, within the quasiharmonic approximation, some other thermal properties of these metals which depend on the phonon density of states, such as the temperature dependence of lattice constant, coefficient of linear thermal expansion, isothermal and adiabatic bulk moduli, heat capacities at constant volume and constant pressure, Griineisen parameter and Debye temperature. The computed results are compared with the experimental findings wherever possible. The comparison shows a generally good agreement between the theoretical values and experimental data for all properties except the discrepancies of phonon frequencies and Debye temperature for Pd, Pt and Au. Further, we modify the parameters of this model for Pd and Pt and obtain the phonon dispersion curves which is in good agreement with experimental data.

Phonon softening in Ni-Mn-Ga alloys

The TA2 phonon dispersion curves of Ni-Mn-Ga alloys with different compositions which transform to different martensitic structures have been measured over a broad temperature range covering both paramagnetic and ferromagnetic phases. The branches show an anomaly (dip) at a wave number that depends on the particular martensitic structure, and there is softening of these anomalous phonons with decreasing temperature. This softening is enhanced below the Curie point, as a consequence of spin-phonon coupling. This effect is stronger for systems with higher electronic concentration.

Spin polarized DFT calculations of metal borides

The metal borides, in particular the diborides and hexaborides, contain stoichiometric forms that include insulators, semiconductors and superconductors. In addition, their end-member structures have high symmetry and two atoms although, in general, substitution(s) of multi-valent ions into the metal site occurs consistent with Vegard’s law. These characteristics allow for fundamental comparison of important physical properties such as superconductivity and insulation within a relatively simple structure type. Our early work1,2 has demonstrated this for the hexaborides and this work compares similar attributes across a broader suite of boride structures. In all cases, theoretical calculations are referenced to structures determined via high resolution neutron or X-ray diffraction experiments.