Les piézoélectriques à haute performance sont des composants clés pour les dispositifs agiles. Il a été démontré récemment que les solutions solides (Ba,Ca)(Ti,Zr)O3 (BCTZ) présentent des propriétés électromécaniques prometteuses. Cependant, les mécanismes microscopiques conduisant à de telles caractéristiques restent à éclaircir, et les investigations théoriques de BCTZ demeurent très limitées à ce jour. En conséquence, cette thèse propose d’étudier les propriétés de différentes compositions de solutions solides (Ba,Ca)TiO3-Ba(Ti,Zr)O3 au moyen de calculs de premiers principes, en mettant l’accent sur la dynamiquedu réseau et sur la compétition entre différentes phases ferroélectriques. Nous nous intéressons d’abord aux quatre composés parents BaTiO3, CaTiO3, BaZrO3 et CaZrO3, afin de comparer leurs propriétés et leurs différentes tendances à la ferroélectricité. Ensuite, le coeur de notre étude est une caractérisation systématique des systèmes binaires (Ba,Ca)TiO3 et Ba(Ti,Zr)O3 en utilisant à la fois l’approximation du cristal virtuel (VCA) et des calculs directs sur supercellules.Lorsqu’on passe continument de BaTiO3 à CaTiO3 dans (Ba,Ca)TiO3, la caractéristique principale est une transformation progressive de la ferroélectricité de type B en type A en raison d’effets stériques, lesquels déterminent en grande partie le comportement du système. En particulier, pour les petites concentrations en calcium, nous avons mis en évidence que laferroélectricité guidée par le site B est globalement affaiblie, conduisant à la disparition de la barrière d’énergie entre différents états polaires et à une polarisation quasi-isotrope. Une amélioration considérable de la réponse piézoélectrique résulte de ces caractéristiques. En passant de BaTiO3 à BaZrO3 dans Ba(Ti,Zr)O3, en revanche, le comportement est dominé parles mouvements coopératifs Zr-Ti et l’électrostatique locale. En particulier, la phase R3m est stabilisée significativement pour les faibles concentrations en zirconium. Sous l’effet d’une augmentation de la concentration en zirconium, le système montre une tendance à la réduction de la distorsion polaire; néanmoins, la ferroélectricité peut être préservée localement dans les régionsriches en titane. Grâce à un modèle électrostatique basé sur un super-réseau BaZrO3/mBaTiO3, nous avons également découvert une activation polaire inattendue pour Zr, en fonction d’un ordre atomique spécifique. Un facteur microscopique expliquant la réponse piézoélectrique exaltée dans BCTZ, pour de faibles concentrations en Ca et Zr, peut donc résider dans l’interaction entre la ferroélectricité affaiblie induite par Ti et la ferroélectricité émergente induite par Ca, interaction produisant une anisotropie minimale pour la polarisation. En outre, notre étude comparative révèle que la physique microscopique spécifique de ces solutions solides limite sévèrement l’applicabilité de l’approximation du cristal virtuel (VCA) à ces systèmes. / High-performance piezoelectrics are key components of various smart devicesand, recently, it has been discovered that (Ba,Ca)(Ti,Zr)O3 (BCTZ) solid solutions show appealingelectromechanical properties. Nevertheless, the microscopic mechanisms leading to suchfeatures are still unclear and theoretical investigations of BCTZ remain very limited. Accordingly,this thesis analyzes the properties of various compositions of (Ba,Ca)TiO3-Ba(Ti,Zr)O3solid solutions by means of first-principles calculations, with a focus on the lattice dynamicsand the competition between different ferroelectric phases. We first analyze the four parentcompounds BaTiO3, CaTiO3, BaZrO3 and CaZrO3 in order to compare their properties andtheir different tendency towards ferroelectricity. Then, the core of our study is a systematiccharacterization of the binary systems (Ba,Ca)TiO3 and Ba(Ti,Zr)O3 within both the virtualcrystal approximation (VCA) and direct supercell calculations. When going from BaTiO3 toCaTiO3 in (Ba,Ca)TiO3, the main feature is a gradual transformation from B-type to A-typeferroelectricity due to steric effects that largely determine the behavior of the system. In particular,for low Ca-concentration we found out an overall weakened B-driven ferroelectricitythat produces the vanishing of the energy barrier between different polar states and resultsin a quasi-isotropic polarization. A sizable enhancement of the piezoelectric response resultsfrom these features. When going from BaTiO3 to BaZrO3 in Ba(Ti,Zr)O3, in contrast, thebehavior is dominated by cooperative Zr-Ti motions and the local electrostatics. In particular,low Zr-concentration produces the further stabilization of the R3m-phase. Then, the systemshows the tendency to globally reduce the polar distortion with increasing Zr-concentration.Nevertheless, ferroelectricity can be locally preserved in Ti-rich regions. We also found out anunexpected polar activation of Zr as a function of specific atomic ordering explained via a basicelectrostatic model based on BaZrO3/mBaTiO3 superlattice. A microscopic factor behind theenhanced piezoelectric response in BCTZ, at low concentration of Ca and Zr, can thus be theinterplay between weakened Ti-driven and emerging Ca-driven ferroelectricity, which producesminimal anisotropy for the polarization. In addition, our comparative study reveals that thespecific microscopic physics of these solid solutions sets severe limits to the applicability of thevirtual crystal approximation (VCA) for these systems.

<p>The elastic constants of Mg_(1-x)Al_xB₂ have been calculated in the regime 0<x<0.25. The calculations show that the ratio, B/G, between the bulk- and the shear-modulus stays well below the empirical ductility limit, 1.75, for all concentrations, indicating that the introduction of Al will not change the brittle behaviour of the material considerably. Furthermore, the tetragonal elastic constant C’ has been calculated for the transition metal alloys Fe-Co, Mo-Tc and W-Re, showing that if a suitable tuning of the alloying is made, these materials have a vanishingly low C'. Thermal expansion calculations of the 4d transition metals have also been performed, showing good agreement with experiment with the exception of Nb and Mo. The calculated phonon dispersions of the 4d metals all give reasonable agreement with experiment. First principles calculations of the thermal expansion of hcp Ti have been performed, showing that this element has a negative thermal expansion along the c-axis which is linked to the closeness of the Fermi level to an electronic topological transition. Calculations of the EOS of fcc Au give support to the suggestion that the ruby pressure scale might underestimate pressures with ~10 GPa at pressures ~150 GPa. The high temperature bcc phase of the group IV metals has been calculated with the novel self-consistent ab-initio dynamical (SCAILD) method. The results show good agreement with experiment, and the free energy resolution of < 1 meV suggests that this method might be suitable for calculating free energy differences between different crystallographic phases as a function of temperature.</p>

The elastic constants of Mg(1-x)AlxB2 have been calculated in the regime 0&lt;x&lt;0.25. The calculations show that the ratio, B/G, between the bulk- and the shear-modulus stays well below the empirical ductility limit, 1.75, for all concentrations, indicating that the introduction of Al will not change the brittle behaviour of the material considerably. Furthermore, the tetragonal elastic constant C’ has been calculated for the transition metal alloys Fe-Co, Mo-Tc and W-Re, showing that if a suitable tuning of the alloying is made, these materials have a vanishingly low C'. Thermal expansion calculations of the 4d transition metals have also been performed, showing good agreement with experiment with the exception of Nb and Mo. The calculated phonon dispersions of the 4d metals all give reasonable agreement with experiment. First principles calculations of the thermal expansion of hcp Ti have been performed, showing that this element has a negative thermal expansion along the c-axis which is linked to the closeness of the Fermi level to an electronic topological transition. Calculations of the EOS of fcc Au give support to the suggestion that the ruby pressure scale might underestimate pressures with ~10 GPa at pressures ~150 GPa. The high temperature bcc phase of the group IV metals has been calculated with the novel self-consistent ab-initio dynamical (SCAILD) method. The results show good agreement with experiment, and the free energy resolution of &lt; 1 meV suggests that this method might be suitable for calculating free energy differences between different crystallographic phases as a function of temperature.

In this thesis I have studied some problems of condensed matter at high pressures and temperatures by means of numerical simulations based on Density Functional Theory (DFT). The stability of MgCO3 and CaCO3 carbonates at the Earth's mantle conditions may play an important role in the global carbon cycle through the subduction of the oceanic crust. By performing ab initio electronic structure calculations, we observed a new high-pressure phase transition within the Pmcn structure of CaCO3. This transformation is characterized by the change of the sp-hybridization state of carbon atom and indicates a change to a new crystal-chemical regime. By performing ab initio Molecular Dynamics simulations we show the new phase to be stable at 250 GPa and 1000K. Thus, the formation of sp3 hybridized bonds in carbonates can explain the stability of MaCO3 and CaCO3 at pressures corresponding to the Earth's lower mantle conditions. We have also calculated phase transition sequence in CaCO3, SrCO3 and BaCO3, and have found that, despite the fact that these carbonates are isostructural and undergo the same type of aragonite to post-aragonite transition, their phase transformation sequences are different at high pressures. The continuous improvement of the high-pressure technique led to the discovery of new composite structures at high pressures and complex phases of many elements in the periodic table have been determined as composite host-guest incommensurate structures. We propose a procedure to accurately describe the structural parameters of an incommensurate phase using ab initio methods by approximating it with a set of analogous commensurate supercells and exploiting the fact that the total energy of the system is a function of structural parameters. By applying this method to the Sc-II phase, we have determined the incommensurate ratio, lattice parameters and Wyckoff positions of Sc-II in excellent agreement with the available experimental data. Moreover, we predict the occurrence of an incommensurate high-pressure phase in Ca from first-principle calculations within this approach. The implementation of DFT in modern electronic structure calculation methods proved to be very successful in predicting the physical properties of a solid at low temperature. One can rigorously describe the thermodynamics of a crystal via the collective excitation of the ionic lattice, and the ab initio calculations give an accurate phonon spectra in the quasi-harmonic approximation. Recently an elegant method to calculate phonon spectra at finite temperature in a self-consistent way by using first principles methods has been developed. Within the framework of self-consistent ab initio lattice dynamics approach (SCAILD) it is possible to reproduce the observed stable phonon spectra of high-temperature bcc phase of Ti, Zr and Hf with a good accuracy. We show that this method gives also a good description of the thermodynamics of hcp and bcc phases of Ti, Zr and Hf at high temperatures, and we provide a procedure for the correct estimation of the hcp to bcc phase transition temperature.

This thesis presents a theoretical study of materials under high pressure using ab initio lattice dynamics based on density functional theory and density functional perturbation theory using both super-cell and linear response approach. Ab initio lattice dynamics using super-cell approach is applied to compare our theoretical predictions with experimental findings. Phonon dispersion curves of fcc α-γ cerium are calculated and compared with inelastic X-ray scattering data. Pressure dependency of phonon density of states in two cubic phases TiO2 allows us to assign the observed cubic phase in experiments to be of fluorite rather than pyrite structure. Dynamical stability of cotunnite TiO2 phase at low pressure can explain the observed quenching phenomena in experiments. Our calculated O2 vibron mode in both ε-ζ phases of solid oxygen supports the hypothesis that both phases are iso-structural. Hydrogen-dense materials attract great attention not only because they open a path to study phenomena related to metallization (superconductivity) of solid hydrogen but also because they are closely related to important industrial applications (hydrogen storage). Using linear response method, we find that metallic fcc-AlH3 is dynamically stabilized in the range of 72-106 GPa and can persist at ambient pressure if finite temperature effects are considered. For SiH4, we test dynamical stability, Raman spectra, zero point energy, and utilize GW calculations for self energy correction. We find that a metallic tetragonal phase of SiH4 can be assigned to the experimentally observed one. Our ab initio lattice dynamics calculations based on density functional perturbation theory predict that fcc-YH3 is a pressure-induced superconductor with a high transition temperature of 40 K at 17.7 GPa. With increasing pressure this material undergoes a superconductor-metal-superconductor transition and the underlying mechanism of this transition can simultaneously explains also the observed metal-insulator transition at 25 GPa in YH3-δ.

The mechanisms of thermal transport in defect-free silicon nanostructures are examined using a combination of lattice dynamics (LD) calculations and the Boltzmann transport equation (BTE). To begin, the thermal conductivity reduction in thin films is examined using a hierarchical method that first predicts phonon transport properties using LD calculations, and then solves the phonon BTE using the lattice Boltzmann method. This approach, which considers all of the phonons in the first Brillouin-zone, is used to assess the suitability of common assumptions used to reduce the computational effort. Specifically, we assess the validity of: (i) neglecting the contributions of optical modes, (ii) the isotropic approximation, (iii) assuming an averaged bulk mean-free path (i.e., the Gray approximation), and (iv) using the Matthiessen rule to combine the effect of different scattering mechanisms. Because the frequency-dependent contributions to thermal conductivity change as the film thickness is reduced, assumptions that are valid for bulk are not necessarily valid for thin films. Using knowledge gained from this study, an analytical model for the length-dependence of thin film thermal conductivity is presented and compared to the predictions of the LD-based calculations. The model contains no fitting parameters and only requires the bulk lattice constant, bulk thermal conductivity, and an acoustic phonon speed as inputs. By including the mode-dependence of the phonon lifetimes resulting from phonon-phonon and phonon-boundary scattering, the model predictions capture the approach to the bulk thermal conductivity better than predictions made using Gray models based on a single lifetime. Both the model and the LD-based method are used to assess a procedure commonly used to extract bulk thermal conductivities from length-dependent molecular dynamics simulation data. Because the mode-dependence of thermal conductivity is not included in the derivation of this extrapolation procedure, using it can result in significant error. Finally, phonon transport across a silicon/vacuum-gap/silicon structure is modelled using lattice dynamics and Landauer theory. The phonons transmit thermal energy across the vacuum gap via atomic interactions between the leads. Because the incident phonons do not encounter a classically impenetrable potential barrier, this mechanism is not a tunneling phenomenon. The heat flux due to phonon transport can be 4 orders of magnitude larger than that due to photon transport predicted from near-field radiation theory.

Thesis (Ph.D.)--Mississippi State University. Department of Electrical and Computer Engineering. / Title from title screen. Includes bibliographical references.

The mechanisms of thermal transport in defect-free silicon nanostructures are examined using a combination of lattice dynamics (LD) calculations and the Boltzmann transport equation (BTE). To begin, the thermal conductivity reduction in thin films is examined using a hierarchical method that first predicts phonon transport properties using LD calculations, and then solves the phonon BTE using the lattice Boltzmann method. This approach, which considers all of the phonons in the first Brillouin-zone, is used to assess the suitability of common assumptions used to reduce the computational effort. Specifically, we assess the validity of: (i) neglecting the contributions of optical modes, (ii) the isotropic approximation, (iii) assuming an averaged bulk mean-free path (i.e., the Gray approximation), and (iv) using the Matthiessen rule to combine the effect of different scattering mechanisms. Because the frequency-dependent contributions to thermal conductivity change as the film thickness is reduced, assumptions that are valid for bulk are not necessarily valid for thin films. Using knowledge gained from this study, an analytical model for the length-dependence of thin film thermal conductivity is presented and compared to the predictions of the LD-based calculations. The model contains no fitting parameters and only requires the bulk lattice constant, bulk thermal conductivity, and an acoustic phonon speed as inputs. By including the mode-dependence of the phonon lifetimes resulting from phonon-phonon and phonon-boundary scattering, the model predictions capture the approach to the bulk thermal conductivity better than predictions made using Gray models based on a single lifetime. Both the model and the LD-based method are used to assess a procedure commonly used to extract bulk thermal conductivities from length-dependent molecular dynamics simulation data. Because the mode-dependence of thermal conductivity is not included in the derivation of this extrapolation procedure, using it can result in significant error. Finally, phonon transport across a silicon/vacuum-gap/silicon structure is modelled using lattice dynamics and Landauer theory. The phonons transmit thermal energy across the vacuum gap via atomic interactions between the leads. Because the incident phonons do not encounter a classically impenetrable potential barrier, this mechanism is not a tunneling phenomenon. The heat flux due to phonon transport can be 4 orders of magnitude larger than that due to photon transport predicted from near-field radiation theory.

The calculation of the thermal conductivity of condensed matter has posed a significant challenge to engineers and scientists for almost a century. Thermal conductivity models have been successfully applied to many materials however many challenges still remain. One serious challenge is the inability of current thermal conductivity models to calculate the thermal conductivity of highly complex materials. Another challenge is managing error introduced by using an effective interatomic potential, for many materials this problem is exacerbated because their effective potentials have not been extensively used or characterized. Recent interest in nanostructures has initiated a new set of challenges and unanswered questions. This work addresses different aspects of the aforementioned challenges by using zeolite MFI and gallium nitride as case studies.

Thesis (M.S.)--Mississippi State University. Department of Physics and Astronomy. / Title from title screen. Includes bibliographical references.