Search

Lead-based halide perovskite crystals are shown to have strongly anharmonic structural dynamics. This behavior is important because it may be the origin of their exceptional photovoltaic properties. The double perovskite, Cs2AgBiBr6, has been recently studied as a lead-free alternative for optoelectronic applications. However, it does not exhibit the excellent photovoltaic activity of the lead-based halide perovskites. Therefore, to explore the correlation between the anharmonic structural dynamics and optoelectronic properties in lead-based halide perovskites, the structural dynamics of Cs2AgBiBr6 are investigated and are compared to its lead-based analog, CsPbBr3. Using temperature-dependent Raman measurements, it is found that both materials are indeed strongly anharmonic. Nonetheless, the expression of their anharmonic behavior is markedly different. Cs2AgBiBr6 has well-defined normal modes throughout the measured temperature range, while CsPbBr3 exhibits a complete breakdown of the normal-mode picture above 200 K. It is suggested that the breakdown of the normal-mode picture implies that the average crystal structure may not be a proper starting point to understand the electronic properties of the crystal. In addition to our main findings, an unreported phase of Cs2AgBiBr6 is also discovered below approximate to 37 K. Funding Agencies|European Research Concil [850041 - ANHARMONIC]; Alexander von Humboldt FoundationAlexander von Humboldt Foundation; Department of Atomic Energy, Government of IndiaDepartment of Atomic Energy (DAE) [RTI 4007]; SERB IndiaDepartment of Science & Technology (India)Science Engineering Research Board (SERB), India [CRG/2020/003877]; Swedish Research Council (VR)Swedish Research Council [2020-04630, 2021-00486]; Swedish Research CouncilSwedish Research CouncilEuropean Commission [2018-05973]

Understanding heat flow in layered transition metal dichalcogenide (TMD) crystals is crucial for applications exploiting these materials. Despite significant efforts, several basic thermal transport properties of TMDs are currently not well understood, in particular how transport is affected by material thickness and the material's environment. This combined experimental-theoretical study establishes a unifying physical picture of the intrinsic lattice thermal conductivity of the representative TMD MoSe

We put forth an ab initio framework to calculate local moment magnetic interaction parameters, renormalized to treat both the lattice and magnetic systems as a function of temperature $T$. For bcc Fe, magnetic and lattice thermal disorders act in opposition, the former strengthening the Heisenberg-like interactions, while the latter decreasing them. Below ${T}_{C}$, $J$ stays nearly independent of $T$, while around and above ${T}_{C}$, it exhibits a sharp decrease. This remarkable behavior reflects an intricate spin-lattice coupling and its evolution with $T$, in which magnetic interactions and interatomic bonds are each renormalized by the other. This finding is consistent with magnetization data and with the observed softening of magnon and phonon modes at high temperatures. Magnetization as well as magnon and phonon mode softening are discussed.

The proximity to structural phase transitions in IV-VI thermoelectric materials is one of the main reasons for their large phonon anharmonicity and intrinsically low lattice thermal conductivity $\kappa$. However, the $\kappa$ of GeTe increases at the ferroelectric phase transition near $700$ K. Using first-principles calculations with the temperature dependent effective potential method, we show that this rise in $\kappa$ is the consequence of negative thermal expansion in the rhombohedral phase and increase in the phonon lifetimes in the high-symmetry phase. Negative thermal expansion increases phonon group velocities, which counteracts enhanced anharmonicity of phonon modes and boosts $\kappa$ close to the phase transition in the rhombohedral phase. A drastic decrease in the anharmonic force constants in the cubic phase increases the phonon lifetimes and $\kappa$. Strong anharmonicity near the phase transition induces non-Lorentzian shapes of the phonon power spectra. To account for these effects, we implement a novel method of calculating $\kappa$ based on the Green-Kubo approach and find that the Boltzmann transport equation underestimates $\kappa$ near the phase transition. Our findings elucidate the influence of structural phase transitions on $\kappa$ and provide guidance for design of better thermoelectric materials.

Inelastic neutron scattering on a single crystal of silicon was performed at temperatures from 100 to 1500 K. These experimental data were reduced to obtain phonon spectral intensity at all wave vectors $\stackrel{P\vec}{Q}$ and frequencies $\ensuremath{\omega}$ in the first Brillouin zone. Thermal broadenings of the phonon peaks were obtained by fitting and by calculating with an iterative ab initio method that uses thermal atom displacements on an ensemble of superlattices. Agreement between the calculated and experimental broadenings was good, with possible discrepancies at the highest temperatures. Distributions of phonon widths versus phonon energy had similar shapes for computation and experiment. These distributions grew with temperature but maintained similar shapes. Parameters from the ab initio calculations were used to obtain the thermal conductivity from the Boltzmann transport equation, which was in good agreement with experimental data. Despite the high group velocities of longitudinal acoustic phonons, their shorter lifetimes reduced their contribution to the thermal conductivity, which was dominated by transverse acoustic modes.

First-principles calculations were used to study the Fermi surface of body-centered cubic vanadium at elevated temperatures. Supercell calculations accounted for effects of thermal atom displacements on band energies, and band unfolding was used to project the spectral weight of the electron states into the Brillouin zone of a standard bcc unit cell. An electronic topological transition (ETT, or Lifshitz transition) occurred near the $\mathrm{\ensuremath{\Gamma}}$ point with increasing temperature, but the large thermal smearings from the atomic disorder and the Fermi-Dirac distribution reduced the effect of this ETT on the electron-phonon interactions. The phonon dispersions showed thermal stiffening of their Kohn anomalies near the $\mathrm{\ensuremath{\Gamma}}$ point and of the longitudinal N phonon mode. In general the effects of the ETT were overcome by the thermal smearing of the Fermi surface that reduces the spanning vector densities for anomalous phonon modes.

Altres ajuts: ICN2 is supported by the CERCA Programme/Generalitat de Catalunya. In this study, we use transient thermal gratings-a non-contact, laser-based thermal metrology technique with intrinsically high accuracy-to investigate room-temperature phonon-mediated thermal transport in two nanoporous holey silicon membranes with limiting dimensions of 120 nm and 250 nm, respectively. We compare the experimental results with ab initio calculations of phonon-mediated thermal transport according to the phonon Boltzmann transport equation (BTE) using two different computational techniques. We find that the calculations conducted within the Casimir framework, i.e., based on the BTE with the bulk phonon dispersion and diffuse scattering from surfaces, are in quantitative agreement with the experimental data and thus conclude that this framework is adequate for describing phonon-mediated thermal transport in silicon nanostructures with feature sizes of the order of 100 nm.

Using ab initio calculations, we have analyzed the influence of anharmonic effects on the electronic structure and the phonon-dispersion relations of body-centered-cubic (bcc) niobium (Nb) and investigated the temperature dependence of the Kohn anomaly in this metal. A comparison of the results obtained in the framework of the temperature-dependent effective potential method with those derived within the quasiharmonic approximation demonstrates the importance of the explicit treatment of the finite-temperature effects upon the theoretical description of bcc Nb lattice dynamics. In agreement with experimental results, the inclusion of anharmonic vibrations in our calculations leads to the disappearance of the Kohn anomaly for the acoustic mode in a vicinity of the Gamma point with increasing temperature. Moreover, the calculated phonon self-energy indicates that the origin of the temperature dependence is related to the change of the electronic structure. We have calculated the temperature dependence of the electronic spectral function and analyzed the Fermi surface of Nb. A significant temperature-induced smearing of the electronic states has been identified as the origin of the disappearance of the Kohn anomaly in Nb at elevated temperature. Funding Agencies|Swedish Research Council (VR)Swedish Research Council [2015-04391, 2019-05551]; Swedish Government Strategic Research Areas in Materials Science on Functional Materials at Linkoping University (Faculty Grant SFO-Mat-LiU) [2009-00971]; Knut and Alice Wallenberg Foundation (Wallenberg Scholar Grant) [KAW-2018.0194]; Swedish Foundation for Strategic Research (SSF)Swedish Foundation for Strategic Research [EM16-0004]; Russian Science FoundationRussian Science Foundation (RSF) [18-12-00492]

The lead-free halide double perovskite class of materials offers a promising venue for resolving issues related to toxicity of Pb and long-term stability of the lead-containing halide perovskites. We present a first-principles study of the lattice vibrations in Cs2AgBiBr6, the prototypical compound in this class and show that the lattice dynamics of Cs2AgBiBr6 is highly anharmonic, largely in regards to tilting of AgBr6 and BiBr6 octahedra. Using an energy- and temperature-dependent phonon spectral function, we then show how the experimentally observed cubic-to-tetragonal phase transformation is caused by the collapse of a soft phonon branch. We finally reveal that the softness and anharmonicity of Cs2AgBiBr6 yield an ultralow thermal conductivity, unexpected of high-symmetry cubic structures. Funding Agencies|Swedish Research Council (VR)Swedish Research Council [2019-05551]; Swedish Government Strategic Research Area in Materials Science on Advanced Functional Materials at Linkoping University (Faculty Grant SFO-Mat-LiU) [2009-00971]; Swedish e-science Research Centre (SeRC); Knut and Alice Wallenberg Foundation (Wallenberg Scholar Grant) [KAW-2018.0194]; Ministry of Science and High Education of the Russian Federation of National University of Science and Technology (NUST) "MISIS" [K2-2019-001, 211]; Swedish Research CouncilSwedish Research Council [2016-07213]

Lead chalcogenides have exceptional thermoelectric properties and intriguing anharmonic lattice dynamics underlying their low thermal conductivities. An ideal material for thermoelectric efficiency is the phonon glass–electron crystal, which drives research on strategies to scatter or localize phonons while minimally disrupting electronic-transport. Anharmonicity can potentially do both, even in perfect crystals, and simulations suggest that PbSe is anharmonic enough to support intrinsic localized modes that halt transport. Here, we experimentally observe high-temperature localization in PbSe using neutron scattering but find that localization is not limited to isolated modes – zero group velocity develops for a significant section of the transverse optic phonon on heating above a transition in the anharmonic dynamics. Arrest of the optic phonon propagation coincides with unusual sharpening of the longitudinal acoustic mode due to a loss of phase space for scattering. Our study shows how nonlinear physics beyond conventional anharmonic perturbations can fundamentally alter vibrational transport properties., To optimize the performance of lead chalcogenides for thermoelectric applications, strategies to further reduce the crystal’s thermal conductivity is required. Here, the authors discover anharmonic localized vibrations in PbSe crystals for optimizing the crystal’s vibrational transport properties.

With the examples of the C K-edge in graphite and the B K-edge in hexagonal boron nitride, we demonstrate the impact of vibrational coupling and lattice distortions on the X-ray absorption near-edge structure (XANES) in two-dimensional layered materials. Theoretical XANES spectra are obtained by solving the Bethe–Salpeter equation of many-body perturbation theory, including excitonic effects through the correlated motion of the core hole and excited electron. We show that accounting for zero-point motion is important for the interpretation and understanding of the measured X-ray absorption fine structure in both materials, in particular for describing the σ*-peak structure. Funding agencies: Swedish Government Strategic Research Area in Materials Science on Functional Materials at Linkoping University [2009 00971]; Knut and Alice Wallenbergs Foundation; Swedish Energy Research [43606-1]; Carl Trygger Foundation [CTS16:303, CTS14:310]; JSPS KA

We present a first-principles theoretical approach to calculate temperature dependent phonon dispersions in bcc Fe, which captures finite temperature spin-lattice coupling by treating thermal disorder in both the spin and lattice systems simultaneously. With increasing temperature, thermal atomic displacements are found to induce increasingly large fluctuations in local magnetic moment magnitudes. The calculated phonon dispersions of bcc Fe show excellent agreement with measured data over a wide range of temperatures both above and below the magnetic and structural transition temperatures, suggesting the applicability of the developed approach to other magnetic materials. Funding Agencies|National Science Foundation EFRI 2-DARE program [1433467]; Boston College dissertation fellowship award; National Science FoundationNational Science Foundation (NSF) [ACI-1548562]; NSFNational Science Foundation (NSF) [ACI-1445606]

The prototypical antiferroelectric PbZrO3 has several unsettled questions, such as the nature of the antiferroelectric transition, a possible intermediate phase, and the microscopic origin of the Pbam ground state. Using first-principles calculations, we show that no phonon becomes truly soft at the cubic-to-Pbam transition temperature, and the order-disorder character of this transition is clearly demonstrated based on molecular dynamics simulations and potential energy surfaces. The out-of-phase octahedral tilting is an important degree of freedom, which can collaborate with other phonon distortions and form a complex energy landscape with multiple minima Candidates of the possible intermediate phase are suggested based on the calculated kinetic barriers between energy minima, and the development of a first-principles-based effective Hamiltonian. The use of this latter scheme further reveals that specific bilinear interactions between local dipoles and octahedral tiltings play a major role in the formation of the Pbam ground state, which contrasts with most of the previous explanations. Funding Agencies|ONR [N00014-17-1-2818]; DARPA under the MATRIX program [HR0011-15-2-0038]; Priority Academic Program Development (PAPD) of Jiangsu Higher Education Institutions; Arkansas High Performance Computer Center at the University of Arkansas; Cloud Pilot project

Structural phase transitions and soft phonon modes pose a long-standing challenge to computing electron-phonon ($e$-ph) interactions in strongly anharmonic crystals. Here we develop a first-principles approach to compute $e$-ph scattering and charge transport in materials with anharmonic lattice dynamics. Our approach employs renormalized phonons to compute the temperature-dependent $e$-ph coupling for all phonon modes, including the soft modes associated with ferroelectricity and phase transitions. We show that the electron mobility in cubic ${\mathrm{SrTiO}}_{3}$ is controlled by scattering with longitudinal optical phonons at room temperature and with ferroelectric soft phonons below 200 K. Our calculations can accurately predict the temperature dependence of the electron mobility in ${\mathrm{SrTiO}}_{3}$ between 150--300 K, and reveal the microscopic origin of its roughly ${T}^{\ensuremath{-}3}$ trend. Our approach enables first-principles calculations of $e$-ph interactions and charge transport in broad classes of crystals with phase transitions and strongly anharmonic phonons.

This work presents an integrated approach to study the crystal chemistry and phonon heat capacity of complex layered oxides. Two quaternary delafossites are synthesized from ternary parent compounds and copper monohalides via a topochemical exchange reaction that preserves the honeycomb ordering of the parent structures. For each compound, Rietveld refinement of the powder X-ray diffraction patterns is examined in both monoclinic C2/ c and rhombohedral R3̅ m space groups. Honeycomb ordering occurs only in the monoclinic space group. Bragg peaks associated with honeycomb ordering acquire an asymmetric broadening known as the Warren line shape that is commonly observed in layered structures with stacking disorder. Detailed TEM analysis confirms honeycomb ordering within each layer in both title compounds and establishes a twinning between the adjacent layers instead of the more conventional shifting or skipping stacking faults. The structural model is then used to calculate phonon dispersions and heat capacity from first principles. In both compounds, the calculated heat capacity accurately describes the experimental data. The integrated approach presented here offers a platform to carefully analyze the phonon heat capacity in complex oxides where the crystal structure can produce magnetic frustration. Isolating phonon contribution from total heat capacity is a necessary and challenging step toward a quantitative study of spin liquid materials with exotic magnetic excitations such as spinons and Majorana fermions. A quantitative understanding of phonon density of states based on crystal chemistry as presented here also paves the way toward higher efficiency thermoelectric materials.

Despite the widespread use of silicon in modern technology, its peculiar thermal expansion is not well understood. Adapting harmonic phonons to the specific volume at temperature, the quasiharmonic approximation, has become accepted for simulating the thermal expansion, but has given ambiguous interpretations for microscopic mechanisms. To test atomistic mechanisms, we performed inelastic neutron scattering experiments from 100 K to 1,500 K on a single crystal of silicon to measure the changes in phonon frequencies. Our state-of-the-art ab initio calculations, which fully account for phonon anharmonicity and nuclear quantum effects, reproduced the measured shifts of individual phonons with temperature, whereas quasiharmonic shifts were mostly of the wrong sign. Surprisingly, the accepted quasiharmonic model was found to predict the thermal expansion owing to a large cancellation of contributions from individual phonons.

The effects of temperature and pressure on the phonons of GaN were calculated for both the wurtzite and zinc-blende structures. The quasiharmonic approximation (QHA) gave reasonable results for the temperature dependence of the phonon DOS at zero pressure but unreliably predicted the combined effects of temperature and pressure. Pressure was found to change the explicit anharmonicity, altering the thermal shifts of phonons and more notably qualitatively changing the evolution of phonon lifetimes with increasing temperature. These effects were largest for the optical modes, and phonon frequencies below approximately 5 THz were adequately predicted with the QHA. The elastic anisotropies of GaN in both wurtzite and zinc-blende structures were calculated from the elastic constants as a function of pressure at 0 K. The elastic anisotropy increased with pressure until reaching elastic instabilities at 40 GPa (zinc blende) and 65 GPa (wurtzite). The calculated instabilities are consistent with proposed transformation pathways to rocksalt GaN and place upper bounds on the pressures at which wurtzite and zinc-blende GaN can be metastable. Funding Agencies|Office of Science of the U.S. Department of Energy [DE-AC02-05CH11231]; Department of Energy through the Carnegie-DOE Alliance Centers Stewardship Sciences Academic Alliance Program; Department of Energy through Office of Science Graduate Fellowship Program (DOE SCGF) [DE-AC05-06OR23100]

International audience; This work presents an integrated approach to study the crystal chemistry and phonon heat capacity of complex layered oxides. Two quaternary delafossites are synthesized from ternary parent compounds and copper monohalides via a topochemical exchange reaction that preserves the honeycomb ordering of the parent structures. For each compound, Rietveld refinement of the powder X-ray diffraction patterns is examined in both monoclinic C2/c and rhombohedral R (3) over barm space groups. Honeycomb ordering occurs only in the monoclinic space group. Bragg peaks associated with honeycomb ordering acquire an asymmetric broadening known as the Warren line shape that is commonly observed in layered structures with stacking disorder. Detailed TEM analysis confirms honeycomb ordering within each layer in both title compounds and establishes a twinning between the adjacent layers instead of the more conventional shifting or skipping stacking faults. The structural model is then used to calculate phonon dispersions and heat capacity from first principles. In both compounds, the calculated heat capacity accurately describes the experimental data. The integrated approach presented here offers a platform to carefully analyze the phonon heat capacity in complex oxides where the crystal structure can produce magnetic frustration. Isolating phonon contribution from total heat capacity is a necessary and challenging step toward a quantitative study of spin liquid materials with exotic magnetic excitations such as spinons and Majorana fermions. A quantitative understanding of phonon density of states based on crystal chemistry as presented here also paves the way toward higher efficiency thermoelectric materials.

We discuss theoretical description of pressure-induced phase transitions by means of first-principles calculations in the framework of density functional theory. We illustrate applications of theoretical tools that allow one to take into account configurational and vibrational disorders, considering Ti-V alloys as a model system. The universality of the first-principles theory allows us to apply it in studies of different phenomena that occur in the Ti-V system upon compression. Besides the transitions between different crystal structures, we discuss isostructural transitions in bcc Ti-V alloys. Moreover, we present arguments for possible electronic transitions in this system, which may explain peculiar behaviour of elastic properties of V upon compression. Funding Agencies|Ministry of Education and Science of the Russian Federation [14.Y26.31.0005]; SSF [SRL 10-0026]; Russian Foundation for Basic Researches [13-02-00606a]; Tomsk State University Academic D.I. Mendeleev Fund Program; Program of Fundamental Research of State Academies of Sciences

Molecular crystals such as polyethylene are of intense interest as flexible thermal conductors, yet their intrinsic upper limits of thermal conductivity remain unknown. Here, we report a study of the vibrational properties and lattice thermal conductivity of a polyethylene molecular crystal using an ab initio approach that rigorously incorporates nuclear quantum motion and finite temperature effects. We obtain a thermal conductivity along the chain direction of around 160 W m^{-1} K^{-1} at room temperature, providing a firm upper bound for the thermal conductivity of this molecular crystal. Furthermore, we show that the inclusion of quantum nuclear effects significantly impacts the thermal conductivity by altering the phase space for three-phonon scattering. Our computational approach paves the way for ab initio studies and computational material discovery of molecular solids free of any adjustable parameters.

Phonon partial densities of states (pDOS) of ^(57)Fe_3 C were measured from cryogenic temperatures through the Curie transition at 460 K using nuclear resonant inelastic x-ray scattering. The cementite pDOS reveal that low-energy acoustic phonons shift to higher energies (stiffen) with temperature before the magnetic transition. This unexpected stiffening suggests strongly nonharmonic vibrational behavior that impacts the thermodynamics and elastic properties of cementite. Density functional theory calculations reproduced the anomalous stiffening observed experimentally in cementite by accounting for phonon-phonon interactions at finite temperatures. The calculations show that the low-energy acoustic phonon branches with polarizations along the [010] direction are largely responsible for the anomalous thermal stiffening. The effect was further localized to the motions of the Fe_(II) site within the orthorhombic structure, which participates disproportionately in the anomalous phonon stiffening.

The knowledge of lattice thermal conductivity of materials under realistic conditions is vitally important since many modern technologies require either high or low thermal conductivity. Here, we propose a theoretical model for determining lattice thermal conductivity, which takes into account the effect of microstructure. It is based on ab initio description that includes the temperature dependence of the interatomic force constants and treats anharmonic lattice vibrations. We choose ScN as a model system, comparing the computational predictions to the experimental data by time-domain thermoreflectance. Our experimental results show a trend of reduction in lattice thermal conductivity with decreasing domain size predicted by the theoretical model. These results suggest a possibility to control thermal conductivity by microstructural tailoring and provide a predictive tool for the effect of the microstructure on the lattice thermal conductivity of materials based on ab initio calculations. Funding Agencies|European Research Council under the European Communitys Seventh Framework Programme [FP/2007-2013]; ERC [335383]; Swedish Government Strategic Research Area in Materials Science on Functional Materials at Linkoping University [2009 00971]; Swedish Research Council [2012-4430, 2016-03365, 330-2014-6336, 2014-4750, 637-2013-7296]; Linnaeus Environment LiLi-NFM; Swedish Foundation for Strategic Research (SSF) through the Future Research Leaders 5 Program; NanoCaTe project (FP7) [604647]; National University of Singapore Startup Grant

The thermoelectric properties of 54 different group 4 half-Heusler (HH) alloys have been studied from first principles. Electronic transport was studied with density functional theory using hybrid functionals facilitated by the $\mathbf{k} \cdot \mathbf{p}$ method, while the temperature dependent effective potential method was used for the phonon contributions to the figure of merit $ZT$. The phonon thermal conductivity was calculated including anharmonic phonon-phonon, isotope, alloy and grain-boundary scattering. HH alloys have an ${\it XYZ}$ composition and those studied here are in the group 4-9-15 (Ti,Zr,Hf)(Co,Rh,Ir)(As,Sb,Bi) and group 4-10-14 (Ti,Zr,Hf)(Ni,Pd,Pt)(Ge,Sn,Pb). The electronic part of the thermal conductivity was found to significantly impact $ZT$ and thus the optimal doping level. Furthermore, the choice of functional was found to significantly affect thermoelectric properties, particularly for structures exhibiting band alignment features. The intrinsic thermal conductivity was significantly reduced when alloy and grain boundary scattering were accounted for, which also reduced the spread in thermal conductivity. It was found that sub-lattice disorder on the ${\it Z}$-site, i.e. the site occupied by group 14 or 15 elements, was more effective than ${\it X}$-site substitution, occupied by group 4 elements. The calculations confirmed that ZrNiSn, ZrCoSb and ZrCoBi based alloys display promising thermoelectric properties. A few other n-type and p-type compounds were also predicted to be potentially excellent thermoelectric materials, given that sufficiently high charge carrier concentrations can be achieved. This study provides insight into the thermoelectric potential of HH alloys and casts light on strategies to optimize thermoelectric performance of multicomponent alloys.

Progress in materials science through thermodynamic modelling may rest crucially on access to a database, such as that developed by Scientific Group Thermodata Europe (SGTE) around 1990. It gives the Gibbs energy G(T)of the elements in the form of series as a function of temperature, i.e. essentially a curve fitting to experimental data. In the light of progress in theoretical understanding and first-principles calculation methods, the possibility for an improved database description of the thermodynamics of the elements has become evident. It is the purpose of this paper to provide a framework for such work. Lattice vibrations, which usually give the major contribution to G(T), are treated in some detail with a discussion of neutron scattering studies of anharmonicity in aluminium, first-principles calculations including ab initio molecular dynamics (AIMD), and the strength and weakness of analytic model representations of data. Similarly, electronic contributions to G(T) are treated on the basis of the density of states N(E) for metals, with emphasis on effects at high T. Further, we consider G(T) below 300 K, which is not covered by SGTE. Other parts in the paper discuss metastable and dynamically unstable lattices, G(T) in the region of superheated solids and the requirement on a database in the calculation of phase diagrams.

We review basic ideas behind state-of-the-art techniques for first-principles theoretical simulations of the phase stabilities and properties of alloys. We concentrate on methods that allow for an efficient treatment of compositional and thermal disorder effects. In particular, we present novel approach to evaluate free energy for strongly anharmonic systems. Theoretical tools are then employed in studies of two materials systems relevant for nuclear energy applications: Fe-Cr and Zr-based alloys. In particular, we investigate the effect of hydrostatic pressure and multicomponent alloying on the mixing enthalpy of Fe-Cr alloys, and show that in the ferromagnetic state both of them reduce the alloy stability at low Cr concentration. For Zr-Nb alloys, we demonstrate how microscopic parameters calculated from first-principles can be used in higher-level models.

Wurtzite aluminium nitride (AlN) is a technologically important wide-band-gap semiconductor with an unusually high thermal conductivity, used in optical applications and as a heatsink substrate. Explaining many of its properties depends on an accurate description of its lattice dynamics, which have thus far only been captured in the quasiharmonic approximation. In this work, we show that anharmonic effects have a considerable impact on its phase stability and transport properties, since they are much stronger in the rocksalt phase. We construct a theoretical pressure-temperature phase diagram of AlN, demonstrating that the rocksalt phase is stabilized by increasing temperature, with respect to the wurtzite phase. We recover the thermal conductivity of the wurtzite phase (320Wm^(-1)K^(-1) under ambient conditions) and compute the hitherto unknown thermal conductivity of the rocksalt phase (81Wm^(-1)K^(-1)). We also show that the electronic band gap decreases with temperature. These findings provide further evidence that anharmonic effects cannot be ignored in simulations of materials intended for high-temperature applications.

Ab initio molecular dynamics, supported by inelastic neutron scattering and nuclear resonant inelastic x-ray scattering, showed an anomalous thermal softening of the M_{5}^{-} phonon mode in B2-ordered FeTi that could not be explained by phonon-phonon interactions or electron-phonon interactions calculated at low temperatures. A computational investigation showed that the Fermi surface undergoes a novel thermally driven electronic topological transition, in which new features of the Fermi surface arise at elevated temperatures. The thermally induced electronic topological transition causes an increased electronic screening for the atom displacements in the M_{5}^{-} phonon mode and an adiabatic electron-phonon interaction with an unusual temperature dependence.

The interest in improving the thermoelectric response of bulk materials has received a boost after it has been recognized that layered materials, in particular SnSe, show a very large thermoelectric figure of merit. This result has received great attention while it is now possible to conceive other similar materials or experimental methods to improve this value. Before we can now think of engineering this material it is important we understand the basic mechanism that explains this unusual behavior, where very low thermal conductivity and a high thermopower result from a delicate balance between the crystal and electronic structure. In this Letter, we present a complete temperature evolution of the Seebeck coefficient as the material undergoes a soft crystal transformation and its consequences on other properties within SnSe by means of first-principles calculations. Our results are able to explain the full range of considered experimental temperatures. Funding Agencies|NSF [1434897]; American Chemical Society Petroleum Research Fund [54075-ND10]; ULg; Actions De Recherche Concertees (ACR) grants from Communaute Francaise de Belgique [10/15-03, 15/19-09]; COST network EUSpec [MP1306]; COST network XLIC [CM1204]; EU [RI-312763]; Knut & Alice Wallenberg Foundation (KAW) project "Isotopic Control for Ultimate Material Properties"; Swedish Foundation for Strategic Research (SSF) program [SRL10-002]

A first-principles nonequilibrium molecular dynamics (NEMD) study employing the color-diffusion algorithm has been conducted to obtain the bulk ionic conductivity and the diffusion constant of gadolinium-doped cerium oxide (GDC) in the 850-1150 K temperature range. Being a slow process, ionic diffusion in solids usually requires simulation times that are prohibitively long for ab initio equilibrium molecular dynamics. The use of the color-diffusion algorithm allowed us to substantially speed up the oxygen-ion diffusion. The key parameters of the method, such as field direction and strength as well as color-charge distribution, have been investigated and their optimized values for the considered system have been determined. The calculated ionic conductivity and diffusion constants are in good agreement with available experimental data. Funding Agencies|Swedish Energy Agency (STEM) [355151]; Carl Tryggers Foundation [CTS 14:433]; Swedish Research Council (VR) [2014-5993, 2014-4750, 637-2013-7296]; LiLi-NFM; Swedish Government Strategic Research Area Grant in Materials Science

We develop a method to accurately and efficiently determine the vibrational free energy as a function of temperature and volume for substitutional alloys from first principles. Taking Ti1-xAlxN alloy as a model system, we calculate the isostructural phase diagram by finding the global minimum of the free energy corresponding to the true equilibrium state of the system. We demonstrate that the vibrational contribution including anharmonicity and temperature dependence of the mixing enthalpy have a decisive impact on the calculated phase diagram of a Ti1-xAlxN alloy, lowering the maximum temperature for the miscibility gap from 6560 to 2860 K. Our local chemical composition measurements on thermally aged Ti0.5Al0.5N alloys agree with the calculated phase diagram. Funding Agencies|Swedish Foundation for Strategic Research Programs (Stiftelsen for Strategisk Forskning) [SRL10-0026, RMA08-0069]; FUNCASE; Swedish Research Council (Vetenskapsradet) [2012-4401, 621-2011-4426, 2015-04391, 637-2013-7296, 621-2011-4417, 330-2014-6336]; Swedish Government Strategic Research Area in Materials Science on Functional Materials at Linkoping University (Faculty Grant SFO-Mat-LiU) [2009 00971]; Swedish Governmental Agency for Innovation Systems (Vinnova) through the M-ERA.net project MC2; SECO Tools AB; Marie Sklodowska Curie Actions, Cofund [INCA 600398]; Erasmus Mundus Joint European Doctoral Programme DocMASE; Ministry of Education and Science of the Russian Federation [K2-2016-013]; Deutsche Forschungsgemeinschaft; Federal State Government of Saarland [INST 256/298-1 FUGG]

Lead chalcogenides such as PbS, PbSe, and PbTe are of interest for their exceptional thermoelectric properties and strongly anharmonic lattice dynamics. Although PbTe has received the most attention, PbSe has a lower thermal conductivity despite being stiffer, a trend that prior first-principles calculations have not reproduced. Here, we use ab-initio calculations that explicitly account for strong anharmonicity to identify the origin of this low thermal conductivity as an anomalously large anharmonic interaction, exceeding in strength that in PbTe, between the transverse optic and longitudinal acoustic branches. The strong anharmonicity is reflected in the striking observation of an intrinsic localized mode that forms in the acoustic frequencies. Our work shows the deep insights into thermal phonons that can be obtained from ab-initio calculations that are not confined to the weak limit of anharmonicity., Comment: 18 pages, 5 figures

We revisit the color-diffusion algorithm [Aeberhard et al., Phys. Rev. Lett. 108, 095901 (2012)] in non equilibrium ab initio molecular dynamics (NE-AIMD) and propose a simple efficient approach for the estimation of monovacancy jump rates in crystalline solids at temperatures well below melting. Color-diffusion applied to monovacancy migration entails that one lattice atom (colored atom) is accelerated toward the neighboring defect site by an external constant force F. Considering bcc molybdenum between 1000 and 2800 K as a model system, NE-AIMD results show that the colored-atom jump rate k(NE) increases exponentially with the force intensity F, up to F values far beyond the linear-fitting regime employed previously. Using a simple model, we derive an analytical expression which reproduces the observed k(NE)(F) dependence on F. Equilibrium rates extrapolated by NE-AIMD results are in excellent agreement with those of unconstrained dynamics. The gain in computational efficiency achieved with our approach increases rapidly with decreasing temperatures and reaches a factor of 4 orders of magnitude at the lowest temperature considered in the present study. Funding Agencies|Knut and Alice Wallenberg Foundation [2011.0094]; Swedish Research Council (VR) [621-2011-4417, 2015-04391, 637-2013-7296, 330-2014-336]; Linkoping Linnaeus Initiative LiLi-NFM [2008-6572]; Swedish Government Strategic Research Area Grant in Materials Science on Advanced Functional Materials [MatLiU 2009-00971]; Ministry of Education and Science of the Russian Federation [14.Y26.31.0005]; Tomsk State University Academic D. I. Mendeleev Fund Program

Although the rutile structure of TiO$_2$ is stable at high temperatures, the conventional quasiharmonic approximation predicts that several acoustic phonons decrease anomalously to zero frequency with thermal expansion, incorrectly predicting a structural collapse at temperatures well below 1000\,K. Inelastic neutron scattering was used to measure the temperature dependence of the phonon density of states (DOS) of rutile TiO$_2$ from 300 to 1373\,K. Surprisingly, these anomalous acoustic phonons were found to increase in frequency with temperature. First-principles calculations showed that with lattice expansion, the potentials for the anomalous acoustic phonons transform from quadratic to quartic, stabilizing the rutile phase at high temperatures. In these modes, the vibrational displacements of adjacent Ti and O atoms cause variations in hybridization of $3d$ electrons of Ti and $2p$ electrons of O atoms. With thermal expansion, the energy variation in this "phonon-tracked hybridization" flattens the bottom of the interatomic potential well between Ti and O atoms, and induces a quarticity in the phonon potential., Comment: 7 pages, 6 figures, supplemental material (3 figures)

We investigate the harmonic and anharmonic contributions to the phonon spectrum of lead telluride and perform a complete characterization of how thermal properties of PbTe evolve as temperature increases. We analyze the thermal resistivitys variationwith temperature and clarify misconceptions about existing experimental literature. The resistivity initially increases sublinearly because of phase space effects and ultra strong anharmonic renormalizations of specific bands. This effect is the strongest factor in the favorable thermoelectric properties of PbTe, and it explains its limitations at higher T. This quantitative prediction opens the prospect of phonon phase space engineering to tailor the lifetimes of crucial heat carrying phonons by considering different structure or nanostructure geometries. We analyze the available scattering volume between TO and LA phonons as a function of temperature and correlate its changes to features in the thermal conductivity. Funding Agencies|Marie Curie Actions from the European Union [PIIFR-GA-2011-911070]; American Chemical Society Petroleum Research Fund [54075-ND10]; National Science Foundation [1434897]; Action de Recherches Concertees (ARC) from the Communaute Francaise de Belgique [10/15-03]; EU FP7 [RI-283493, RI-312763]; Swedish Research Council (VR) program [637-2013-7296]

Structural phase transitions in epitaxial stoichiometric VN/MgO(011) thin films are investigated using temperature-dependent synchrotron x-ray diffraction (XRD), selected-area electron diffraction (SAED), resistivity measurements, high-resolution cross-sectional transmission electron microscopy, and ab initio molecular dynamics (AIMD). At room temperature, VN has the B1 NaCl structure. However, below T-c = 250 K, XRD and SAED results reveal forbidden (00l) reflections of mixed parity associated with a noncentrosymmetric tetragonal structure. The intensities of the forbidden reflections increase with decreasing temperature following the scaling behavior I proportional to (T-c - T)(1/2). Resistivity measurements between 300 and 4 K consist of two linear regimes resulting from different electron/phonon coupling strengths in the cubic and tetragonal-VN phases. The VN transport Eliashberg spectral function alpha F-2(tr)(h omega), the product of the phonon density of states F(h omega) and the transport electron/phonon coupling strength alpha(2)(tr)(h omega), is determined and used in combination with AIMD renormalized phonon dispersion relations to show that anharmonic vibrations stabilize the NaCl structure at T greater than T-c. Free-energy contributions due to vibrational entropy, often neglected in theoretical modeling, are essential for understanding the room-temperature stability of NaCl-structure VN, and of strongly anharmonic systems in general. Funding Agencies|Swedish Research Council (VR) program [637-2013-7296, 2014-5790, 2009-00971, 2013-4018]; Swedish Government Strategic Research Area Grant in Materials Science [SFO Mat-LiU 2009-00971]; DOE Office of Science by Argonne National Laboratory [DE-AC02-06CH11357]

A recently developed temperature-dependent effective potential method is employed to study the martensitic phase transformation in AuZn. This method is based on ab initio molecular dynamics and allows to obtain finite-temperature lattice vibrational properties. We show that the transversal acoustic TA(2)[110] mode associated with the phase transition is stabilized at 300 K. Temperature evolution of single-phonon dynamic structure factor at the wave vector q= 1/3[1,1,0], associated with phonon softening and Fermi surface nesting, was also studied. (C) 2015 The Authors. Published by Elsevier Ltd.

Synchrotron reflection x-ray thermal diffuse scattering (TDS) measurements, rather than previously reported transmission TDS, are carried out at room temperature and analyzed using a formalism based upon second-order interatomic force constants and long-range Coulomb interactions to obtain quantitative determinations of MgO phonon dispersion relations (h) over bar omega(j) (q), phonon densities of states g((h) over bar omega), and isochoric temperature-dependent vibrational heat capacities c(v) (T). We use MgO as a model system for investigating reflection TDS due to its harmonic behavior as well as its mechanical and dynamic stability. Resulting phonon dispersion relations and densities of states are found to be in good agreement with independent reports from inelastic neutron and x-ray scattering experiments. Temperature-dependent isochoric heat capacities cv (T), computed within the harmonic approximation from (h) over bar omega(j) (q) values, increase with temperature from 0.4 x 10(-4) eV/atom K at 100 K to 1.4 x 10(-4) eV/atom K at 200 K and 1.9 x 10(-4) eV/atom K at 300 K, in excellent agreement with isobaric heat capacity values c(p) (T) between 4 and 300 K. We anticipate that the experimental approach developed here will be valuable for determining vibrational properties of heteroepitaxial thin films since the use of grazing-incidence (theta less than or similar to theta(c), where theta(c) is the density-dependent critical angle) allows selective tuning of x-ray penetration depths to less than or similar to 10 nm. Funding Agencies|Swedish Research Council (VR) [2014-5790, 2013-4018]; Swedish Government Strategic Research Area (SFO) in Materials Science on Advanced Functional Material (MatLiU AFM) [2009-00971]; US Department of Energy (DOE), Office of Science, Office of Basic Energy Sciences [DE-FG02-07ER46383]; DOE Office of Science, Argonne National Laboratory [DE-AC02-06CH11357]

Ti1−xAlxN is a technologically important alloy that undergoes a process of high temperature age-hardening that is strongly influenced by its elastic properties. We have performed first principles calculations of the elastic constants and anisotropy using the newly developed symmetry imposed force constant temperature dependent effective potential method, that include lattice vibrations and therefore the effects of temperature, including thermal expansion and intrinsic anharmonicity. These are compared with in situ high temperature x-ray diffraction measurements of the lattice parameter. We show that anharmonic effects are crucial to the recovery of finite temperature elasticity. The effects of thermal expansion and intrinsic anharmonicity on the elastic constants are of the same order, and cannot be considered separately. Furthermore, the effect of thermal expansion on elastic constants is such that the volume change induced by zero point motion has a significant effect. For TiAlN, the elastic constants soften non-uniformly with temperature: C11 decreases substantially when the temperature increases for all compositions, resulting in an increased anisotropy. These findings suggest that an increased Al content and annealing at higher temperatures will result in a harder alloy. Funding agencies: Swedish Research Council (VR) [621-2011-4426, 621-2012-4401, 637-2013-7296]; Swedish Foundation for Strategic Research (SSF) [RMA08-0069, SRL10-0026]; VINNOVA [2013-02355(MC2)]; Erasmus Mundus Joint European Doctoral Program DocMASE; Ministry of EducationVid tidpunkten för disputation förelåg publikation som manuskript

We present first-principles calculations of the thermal and thermal transport properties of ${\mathrm{Bi}}_{2}{\mathrm{Te}}_{3}$ that combine an ab initio molecular dynamics (AIMD) approach to calculate interatomic force constants (IFCs) along with a full iterative solution of the Peierls-Boltzmann transport equation for phonons. The newly developed AIMD approach allows determination of harmonic and anharmonic interatomic forces at each temperature, which is particularly appropriate for highly anharmonic materials such as ${\mathrm{Bi}}_{2}{\mathrm{Te}}_{3}$. The calculated phonon dispersions, heat capacity, and thermal expansion coefficient are found to be in good agreement with measured data. The lattice thermal conductivity, ${\ensuremath{\kappa}}_{l}$, calculated using the AIMD approach nicely matches measured values, showing better agreement than the ${\ensuremath{\kappa}}_{l}$ obtained using temperature-independent IFCs. A significant contribution to ${\ensuremath{\kappa}}_{l}$ from optic phonon modes is found. Already at room temperature, the phonon line shapes show a notable broadening and onset of satellite peaks reflecting the underlying strong anharmonicity.

The anharmonic lattice dynamics of rock-salt thermoelectric compounds SnTe and PbTe are investigated with inelastic neutron scattering (INS) and first-principles calculations. The experiments show that, surprisingly, although SnTe is closer to the ferroelectric instability, phonon spectra in PbTe exhibit a more anharmonic character. This behavior is reproduced in first-principles calculations of the temperature-dependent phonon self-energy. Our simulations reveal how the nesting of phonon dispersions induces prominent features in the self-energy, which account for the measured INS spectra and their temperature dependence. We establish that the phase-space for three-phonon scattering processes, rather than just the proximity to the lattice instability, is the mechanism determining the complex spectrum of the transverse-optical ferroelectric mode.

We present a theoretical first-principles method to calculate the free energy of a magnetic system in its high-temperature paramagnetic phase, including vibrational, electronic, and magnetic contributions. The method for calculating free energies is based on ab-initio molecular dynamics and combines a treatment of disordered magnetism using disordered local moments molecular dynamics (DLM-MD) with the temperature dependent effective potential (TDEP) method to obtain the vibrational contribution to the free energy. We illustrate the applicability of the method by obtaining the anharmonic free energy for the paramagnetic cubic and the antiferromagnetic orthorhombic phases of chromium nitride. The influence of lattice dynamics on the transition between the two phases is demonstrated by constructing the temperature-pressure phase diagram., Comment: 10 pages, 9 figures

We have developed a thorough and accurate method of determining anharmonic free energies, the temperature dependent effective potential technique (TDEP). It is based on ab initio molecular dynamics followed by a mapping onto a model Hamiltonian that describes the lattice dynamics. The formalism and the numerical aspects of the technique are described in detail. A number of practical examples are given, and results are presented, which confirm the usefulness of TDEP within ab initio and classical molecular dynamics frameworks. In particular, we examine from first principles the behavior of force constants upon the dynamical stabilization of the body centered phase of Zr, and show that they become more localized. We also calculate the phase diagram for 4He modeled with the Aziz et al. potential and obtain results which are in favorable agreement both with respect to experiment and established techniques. Funding Agencies|Knut & Alice Wallenberg Foundation (KAW) project "Isotopic Control for Ultimate Material Properties"||Swedish Research Council (VR)|621-2011-4426|LiLi-NFM||Swedish Foundation for Strategic Research (SSF) program|SRL10-0026|Ministry of Education and Science of the Russian Federation within the framework of Program Scientific and Scientific-Pedagogical Personnel for Innovative Russia|14.B37.21.089014.A18.21.0893|Swedish Government Strategic Research Area Grant in Materials Science