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.

I. Current commutators in quantum electrodynamics. II. Electromagnetic mass differences of the pseudoscalar meson and baryon octets

This thesis is in two parts. In the first section, the operator structure of the singular terms in the equal-time commutator of space and time components of the electromagnetic current is investigated in perturbation theory by establishing a connection with Feynman diagrams. It is made very plausible that the singular term is a c number. Some remarks are made about the same problem in the electrodynamics of a spinless particle.

In the second part, an SU(3) symmetric multi-channel calculation of the electromagnetic mass differences in the pseudoscalar meson and baryon octets is carried out with an attempt to include some of the physics of the crossed (pair annihilation) channel along the lines of the recent work by Ball and Zachariasen. The importance of the tensor meson Regge trajectories is emphasized. The agreement with experiment is poor for the isospin one mass differences, but excellent for those with isospin two.

Dynamical models for the Earth's geoid

The Earth's largest geoid anomalies occur at the lowest spherical harmonic degrees, or longest wavelengths, and are primarily the result of mantle convection. Thermal density contrasts due to convection are partially compensated by boundary deformations due to viscous flow whose effects must be included in order to obtain a dynamically consistent model for the geoid. These deformations occur rapidly with respect to the timescale for convection, and we have analytically calculated geoid response kernels for steady-state, viscous, incompressible, self-gravitating, layered Earth models which include the deformation of boundaries due to internal loads. Both the sign and magnitude of geoid anomalies depend strongly upon the viscosity structure of the mantle as well as the possible presence of chemical layering.

Correlations of various global geophysical data sets with the observed geoid can be used to construct theoretical geoid models which constrain the dynamics of mantle convection. Surface features such as topography and plate velocities are not obviously related to the low-degree geoid, with the exception of subduction zones which are characterized by geoid highs (degrees 4-9). Recent models for seismic heterogeneity in the mantle provide additional constraints, and much of the low-degree (2-3) geoid can be attributed to seismically inferred density anomalies in the lower mantle. The Earth's largest geoid highs are underlain by low density material in the lower mantle, thus requiring compensating deformations of the Earth's surface. A dynamical model for whole mantle convection with a low viscosity upper mantle can explain these observations and successfully predicts more than 80% of the observed geoid variance.

Temperature variations associated with density anomalies in the man tie cause lateral viscosity variations whose effects are not included in the analytical models. However, perturbation theory and numerical tests show that broad-scale lateral viscosity variations are much less important than radial variations; in this respect, geoid models, which depend upon steady-state surface deformations, may provide more reliable constraints on mantle structure than inferences from transient phenomena such as postglacial rebound. Stronger, smaller-scale viscosity variations associated with mantle plumes and subducting slabs may be more important. On the basis of numerical modelling of low viscosity plumes, we conclude that the global association of geoid highs (after slab effects are removed) with hotspots and, perhaps, mantle plumes, is the result of hot, upwelling material in the lower mantle; this conclusion does not depend strongly upon plume rheology. The global distribution of hotspots and the dominant, low-degree geoid highs may correspond to a dominant mode of convection stabilized by the ancient Pangean continental assemblage.

Magnetic anisotropy of BaM ferrites

Using the point charge model with weak field scheme, a calculation of the single ion anisotropy of Fe3+ ionsn at all crystal sites in BaFe12O19 was made by perturbation theory. The results show that in addition to the 2b site, all other sites have non-negligible contributions. © 1983.

Abnormal ferroelectric response and enhancement of piezoelectricity of PbTiO3 under uniaxial compression

The ferroelectricity of rhombohedral PbTiO3 under uniaxial compression is investigated from first-principles study. Upon compression, the ferroelectricity decreases until a critical stress of -29 GPa and then increases with a further increase of the magnitude of the uniaxial compressive stress. We also find that uniaxial compression could enhance piezoelectricity and that the maximum piezoelectric coefficient d(33) occurs at sigma(33)=-49 GPa, which supports the experimentally observed piezoelectric behavior in rhombohedral Pb(Mg1/3Nb2/3O3)-0.32PbTiO(3) [Q. Wan, C. Chen, and Y. P. Shen, J. Appl. Phys. 98, 024103 (2005)]. Our calculated results show that the Pb, Ti, and O atoms have different contributions to the total polarization with increasing the magnitude of uniaxial compressive stress, and that when -sigma(33)>55 GPa, the Ti atoms no longer have contributions to the polarization, which leads to the changes of ferroelectricity and piezoelectricity. (C) 2008 American Institute of Physics.

Studies of tetragonal PbTiO3 subjected to uniaxial stress along the c-axis

Tetragonal PbTiO3 under uniaxial stress along the c-axis is investigated from first-principles. The structural parameters, polarization, and squares of the lowest optical phonon frequencies for E(1TO) and A(1)(1TO) modes at Gamma show abrupt changes near a stress sigma(c) of 1.04 GPa, which is related to the dramatic change of elastic constant c(33) resulting from the uniaxial stress applied along the c-axis. We also find that the uniaxial compressive stress could enhance the piezoelectric stress coefficients, whereas the uniaxial tensile stress could enhance the piezoelectric strain coefficients. It is also found that when the magnitude of uniaxial compressive stress sigma(33) is greater than 12 GPa, PbTiO3 is transformed to the paraelectric tetragonal phase.

Elasticity, band-gap bowing, and polarization of AlxGa1-xN alloys

Elastic constants, the bulk modulus, Young's modulus, band-gap bowing coefficients, spontaneous and piezoelectric polarizations, and piezoelectric coefficients of hexagonal AlxGa1-xN ternary alloys are calculated using first-principles methods. The fully relaxed structures and the structures subjected to homogeneous biaxial and uniaxial tension are investigated. We show that the biaxial tension in the plane perpendicular to the c axis and the uniaxial tension along the c axis all reduce the bulk modulus, whereas they reduce and enhance Young's modulus, respectively. We find that the biaxial and uniaxial tension can enhance the bowing coefficients. We also find that the biaxial tension can enhance the total polarization, while the uniaxial tension will suppress the total polarization. (C) 2008 American Institute of Physics.

Spin Hall effect of excitons with spin-orbit coupling

The center-of-mass motion of a quasi-two-dimensional exciton with spin-orbit coupling (SOC) in the presence of a perpendicular electric field is calculated by perturbation theory. The results indicate that a quasi-two-dimensional exciton with SOC can exhibit the spin Hall effect (SHE), which is similar to two-dimensional electrons and holes. A likely way to establish exciton SHE in experiments and a possible phase transition from dark to bright state driven by SOC are suggested. (c) 2007 American Institute of Physics.

The influence of internal electric fields on the transition energy of InGaN/gaN quantum well

We present a new way to meet the amount of strain relaxation in an InGaN quantum well layer grown on relaxed GaN by calculating and measuring its internal field. With perturbation theory, we also calculate the transition energy of InGaN/GaN SQWs as affected by internal fields. The newly reported experimental data by Graham et al. fit our calculations well on the assumption that the InGaN well layer suffered a 20% strain relaxation, we discuss the differences between our calculated results and the experimental data. Our calculation suggests that with the increase of indium mole fraction in the InGaN/GaN quantum well, the effect of polarization fields on the luminescence of the quantum well will increase. Moreover, our calculation also suggests that an increase in the quantum well width by only one monolayer can result in a large reduction in the transition energy. (c) 2006 Elsevier B.V. All rights reserved.

Hydrate dissociation conditions for gas mixtures containing carbon dioxide, hydrogen, hydrogen sulfide, nitrogen, and hydrocarbons using SAFT

A new method, a molecular thermodynamic model based on statistical mechanics, is employed to predict the hydrate dissociation conditions for binary gas mixtures with carbon dioxide, hydrogen, hydrogen sulfide, nitrogen, and hydrocarbons in the presence of aqueous solutions. The statistical associating fluid theory (SAFT) equation of state is employed to characterize the vapor and liquid phases and the statistical model of van der Waals and Platteeuw for the hydrate phase. The predictions of the proposed model were found to be in satisfactory to excellent agreement with the experimental data.

Observation of quantum beat in Rb by parametric four-wave mixing

Two coupled parametric four-wave-mixing processed in Rb atoms are studied using perturbation theory, which reveals clear evidence of the appearance of quantum beat at 608 cm(-1), corresponding to the energy difference of the 7s - 5d states of Rb atoms, in the parametric four-wave-mixing signals. A pump-probe technique is utilized to observe the quantum beat. Time-varying characteristics of the quantum beat are investigated using time-dependent Fourier transform. The results show that the time-varying characteristics of the quantum beat potential tool for monitoring the dissociation of molecules.

Photoluminescence of GaInP under high pressure

Photoluminescence of GaInP under hydrostatic pressure is investigated. The Gamma valley of disordered GaInP shifts sublinearly upwards with respect to the top of the valence band with increasing pressure and this sublinearity is caused by the nonlinear relationship between lattice constant and hydrostatic pressure. The Gamma valleys of ordered GaInP rise more slowly than that of the disordered one and the relationship between the band gap and the pressure can not be explained in the same way. Taking into account the interactions between the Gamma valley and the folded L valleys, as well as, the X valleys, the experimental pressure dependences of the band gap of ordered GaInP epilayers are calculated and fitted quite well using first order perturbation theory. The results indicate that simultaneous ordering along [111] and [100] directions can occur in ordered GaInP. (C) 1996 American Institute of Physics.

A DOUBLE-QUANTUM-WELL IN A DRIVING LASER FIELD

The dynamic effect of electrons in a double quantum well under the influence of a monochromatic driving laser field is investigated. Closed-form solutions for the quasienergy and Floquet states are obtained with the help of SU(2) symmetry. For the case of weak interlevel coupling, explicit expressions of the quasienergy are presented by the use of perturbation theory, from which it is found that as long as the photon energy is not close to the tunnel splitting, the electron will be confined in an initially occupied eigenstate of the undriven system during the whole evolution process. Otherwise, it will transit between the lowest two levels in an oscillatory behavior.

Strain effects on optical polarisation properties in (11(2)over-bar2) plane GaN films

We present the theoretical results of the electronic band structure of wurtzite GaN films under biaxial strains in the (11 (2) over bar2)-plane The calculations are performed by the kappa p perturbation theory approach through using the effective-mass Hamiltonian for an arbitrary direction The results show that the transition energies decrease with the biaxial strains changing from -0 5% to 0 5% For films of (11 (2) over bar2)-plane, the strains are expected to be anisotropic in the growth plane Such anisotropic strains give rise to valence band mixing which results in dramatic change in optical polarisation property The strain can also result in optical polarisation switching phenomena Finally, we discuss the applications of these properties to the (11 (2) over bar2) plane GaN based light emitting diode and lase diode

Delta-resonance effective mass in medium

Based on the relativistic chiral effective field theory, we study the effective mass of the Delta-resonance in medium by investigating the self-energy of the Delta-resonance related to the pi N decay channel in symmetric nuclear matter. We find that the effective mass of Delta-resonance decreases evidently with increasing nuclear density rho. In our calculation, we also consider the influence of the shifts of the nucleon mass, pion mass and its decay constant due to the restoration of chiral symmetry in medium. The results are roughly consistent with the data given by Lawrence Berkeley National Laboratory.