Phonon-Phonon Interactions in the Polarization Dependence of Raman Scattering

We have found that the polarization dependence of Raman scattering in organic crystals at finite temperatures can only be described by a fourth-rank tensor formalism. This generalization of the second-rank Raman tensor stems from the effect of off-diagonal components in the crystal self-energy on the light scattering mechanism. We thus establish a novel manifestation of phonon-phonon interaction in inelastic light scattering, markedly separate from the better-known phonon lifetime.

Finite Difference Interpolation for Reduction of Grid-Related Errors in Real-Space Pseudopotential Density Functional Theory

The real-space pseudopotential approach is a well-known method for large-scale density functional theory (DFT) calculations. One of its main limitations, however, is the introduction of errors associated with the positioning of the underlying real-space grid, a phenomenon usually known as the "egg-box" effect. The effect can be controlled by using a finer grid, but this raises the cost of the calculations or even undermines their feasibility altogether. Therefore, there is ongoing interest in the reduction of the effect per a given real-space grid. Here, we present a finite difference interpolation of electron orbitals as a means of exploiting the high resolution of the pseudopotential to reduce egg-box effects systematically. We implement the method in PARSEC, a finite difference real-space pseudopotential DFT code, and demonstrate error mitigation and improved convergence at a low additional computational cost.

Diverging Expressions of Anharmonicity in Halide Perovskites

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, Cs₂ AgBiBr₆ , 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 Cs₂ AgBiBr₆ are investigated and are compared to its lead-based analog, CsPbBr₃ . 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. Cs₂ AgBiBr₆ has well-defined normal modes throughout the measured temperature range, while CsPbBr₃ 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 Cs₂ AgBiBr₆ is also discovered below ≈37 K.

Unraveling Heat Transport and Dissipation in Suspended MoSe2 from Bulk to Monolayer

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₂ . Thermal conductivity measurements using Raman thermometry on a large set of clean, crystalline, suspended crystals with systematically varied thickness are combined with ab initio simulations with phonons at finite temperature. The results show that phonon dispersions and lifetimes change strongly with thickness, yet the thinnest TMD films exhibit an in-plane thermal conductivity that is only marginally smaller than that of bulk crystals. This is the result of compensating phonon contributions, in particular heat-carrying modes around ≈0.1 THz in (sub)nanometer thin films, with a surprisingly long mean free path of several micrometers. This behavior arises directly from the layered nature of the material. Furthermore, out-of-plane heat dissipation to air molecules is remarkably efficient, in particular for the thinnest crystals, increasing the apparent thermal conductivity of monolayer MoSe₂ by an order of magnitude. These results are crucial for the design of (flexible) TMD-based (opto-)electronic applications.

Anharmonicity and Ultralow Thermal Conductivity in Lead-Free Halide Double Perovskites

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 Cs_{2}AgBiBr_{6}, the prototypical compound in this class and show that the lattice dynamics of Cs_{2}AgBiBr_{6} is highly anharmonic, largely in regards to tilting of AgBr_{6} and BiBr_{6} 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 Cs_{2}AgBiBr_{6} yield an ultralow thermal conductivity, unexpected of high-symmetry cubic structures.

Intrinsic anharmonic localization in thermoelectric PbSe

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.

Electron-Phonon Scattering in the Presence of Soft Modes and Electron Mobility in SrTiO_{3} Perovskite from First Principles

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 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 SrTiO_{3} between 150-300 K, and reveal the microscopic origin of its roughly T^{-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.

Metastable silica high pressure polymorphs as structural proxies of deep Earth silicate melts

Modelling of processes involving deep Earth liquids requires information on their structures and compression mechanisms. However, knowledge of the local structures of silicates and silica (SiO2) melts at deep mantle conditions and of their densification mechanisms is still limited. Here we report the synthesis and characterization of metastable high-pressure silica phases, coesite-IV and coesite-V, using in situ single-crystal X-ray diffraction and ab initio simulations. Their crystal structures are drastically different from any previously considered models, but explain well features of pair-distribution functions of highly densified silica glass and molten basalt at high pressure. Built of four, five-, and six-coordinated silicon, coesite-IV and coesite-V contain SiO6 octahedra, which, at odds with 3rd Pauling's rule, are connected through common faces. Our results suggest that possible silicate liquids in Earth's lower mantle may have complex structures making them more compressible than previously supposed.

Nuclear quantum effect with pure anharmonicity and the anomalous thermal expansion of silicon

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.

Lattice Thermal Conductivity of Polyethylene Molecular Crystals from First-Principles Including Nuclear Quantum Effects

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.

Two-Step Phase Transition in SnSe and the Origins of its High Power Factor from First Principles

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.

Lattice Vibrations Change the Solid Solubility of an Alloy at High Temperatures

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 Ti_{1-x}Al_{x}N 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 Ti_{1-x}Al_{x}N alloy, lowering the maximum temperature for the miscibility gap from 6560 to 2860 K. Our local chemical composition measurements on thermally aged Ti_{0.5}Al_{0.5}N alloys agree with the calculated phase diagram.

Highly Efficient Free Energy Calculations of the Fe Equation of State Using Temperature-Dependent Effective Potential Method

Free energy calculations at finite temperature based on ab initio molecular dynamics (AIMD) simulations have become possible, but they are still highly computationally demanding. Besides, achieving simultaneously high accuracy of the calculated results and efficiency of the computational algorithm is still a challenge. In this work we describe an efficient algorithm to determine accurate free energies of solids in simulations using the recently proposed temperature-dependent effective potential method (TDEP). We provide a detailed analysis of numerical approximations employed in the TDEP algorithm. We show that for a model system considered in this work, hcp Fe, the obtained thermal equation of state at 2000 K is in excellent agreement with the results of standard calculations within the quasiharmonic approximation.

Thermally Driven Electronic Topological Transition in FeTi

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.

Phonon self-energy and origin of anomalous neutron scattering spectra in SnTe and PbTe thermoelectrics

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, combined with the proximity to the lattice instability, is the mechanism determining the complex spectrum of the transverse-optic ferroelectric mode.

Charge redistribution mechanisms of ceria reduction

Charge redistribution at low oxygen vacancy concentrations in ceria have been studied in the framework of the density functional theory. We propose a model to approach the dilute limit using the results of supercell calculations. It allows one to reproduce the characteristic experimentally observed behavior of composition versus oxygen pressure dependency. We show that in the dilute limit the charge redistribution is likely to be driven by a mechanism different from the one involving electron localization on cerium atoms. We demonstrate that it can involve charge localization on light element impurities.

Impact of thermal vasomotor control on the heart rate variability of newborn infants

The reactivity of the periodic heart rate variability (HRV) to the thermal vasomotor control was studied during quiet sleep in a total of 20 neonates which were classified in 5 groups according to their maturity. Thermal stimulation on the peripheral skin was applied by an air blower at the rates of 5, 2.5 and 1.5 cycles/min. The oscillation of successive R-R intervals of the ECG was determined by a hybrid computer technique during thermal stimulations and compared to spontaneous HRV activity. The indices of the HRV and the total power of the periodic HRV remained unchanged during thermal stimulations. The response of the periodic HRV was selective at the spectral band corresponding to the frequency of stimulation. It was maximal at 5 cycles/min. The relative response was highly significant (P less than 0.001) and it increased with increasing maturity. The state of SGA (small for gestational age) seemed to decrease the vasomotor response of the HRV. The results suggested that the thermal vasomotor control was functional already in small preterm infants and it tended to increase with increasing maturity as does the autonomic cardiac control in general.