More Like this

1. Extreme sensitivity of higher-order interatomic force constants and thermal conductivity to the energy surface roughness of exchange-correlation functionals

   https://doi.org/10.1063/5.0173762  

   Zhou, Hao; Zhou, Shuxiang; Hua, Zilong; Bawane, Kaustubh; Feng, Tianli (November 2023, Applied Physics Letters)

   In this Letter, we report that the fourth-order interatomic force constants (4th-IFCs) are significantly sensitive to the energy surface roughness of exchange-correlation (XC) functionals in density functional theory calculations. This sensitivity, which is insignificant for the second- (2nd-) and third-order (3rd-) IFCs, varies for different functionals in different materials and can cause misprediction of thermal conductivity by several times of magnitude. As a result, when calculating the 4th-IFCs using the finite difference method, the atomic displacement needs to be taken large enough to overcome the energy surface roughness, in order to accurately predict phonon lifetime and thermal conductivity. We demonstrate this phenomenon on a benchmark material (Si), a high-thermal conductivity material (BAs), and a low thermal conductivity material (NaCl). For Si, we find that the LDA, PBE, and PBEsol XC functionals are all smooth to the 2nd- and 3rd-IFCs but all rough to the 4th-IFCs. This roughness can lead to a prediction of nearly one order of magnitude lower thermal conductivity. For BAs, all three functionals are smooth to the 2nd- and 3rd-IFCs, and only the PBEsol XC functional is rough for the 4th-IFCs, which leads to a 40% underestimation of thermal conductivity. For NaCl, all functionals are smooth to the 2nd- and 3rd-IFCs but rough to the 4th-IFCs, leading to a 70% underprediction of thermal conductivity at room temperature. With these observations, we provide general guidance on the calculation of 4th-IFCs for an accurate thermal conductivity prediction. 

   more » « less
2. ALATDYN: A set of Anharmonic LATtice DYNamics codes to compute thermodynamic and thermal transport properties of crystalline solids

   https://doi.org/10.1016/j.cpc.2025.109575  

   Esfarjani, Keivan; Stokes, Harold; Nayeb_Sadeghi, Safoura; Liang, Yuan; Timalsina, Bikash; Meng, Han; Shiomi, Junichiro; Liao, Bolin; Sun, Ruoshi (July 2025, Computer Physics Communications)

   We introduce a lattice dynamics package which calculates elastic, thermodynamic and thermal transport properties of crystalline materials from data on their force and potential energy as a function of atomic positions. The data can come from density functional theory (DFT) calculations or classical molecular dynamics runs performed in a supercell. First, the model potential parameters, which are anharmonic force constants are extracted from the latter runs. Then, once the anharmonic model is defined, thermal conductivity and equilibrium properties at finite temperatures can be computed using lattice dynamics, Boltzmann transport theories, and a variational principle respectively. In addition, the software calculates the mechanical properties such as elastic tensor, Gruneisen parameters and the thermal expansion coefficient within the quasi-harmonic approximation (QHA). Phonons, elastic constants and thermodynamic properties results applied to the germanium crystal will be illustrated. Using the force constants as a force field, one may also perform molecular dynamics (MD) simulations in order to investigate the combined effects of anharmonicity and defect scattering beyond perturbation theory. 

   more » « less
3. Peculiar anharmonicity of Ruddlesden Popper metal halides: temperature-dependent phonon dephasing

   https://doi.org/10.1039/D1MH01010G  

   Rojas-Gatjens, Esteban; Silva-Acuña, Carlos; Kandada, Ajay Ram (January 2022, Materials Horizons)

   The anharmonicity of the Ruddlesden Popper metal-halide lattice, and its consequences for their electronic and optical properties, are paramount in their basic semiconductor physics. It is thus critical to identify specific anharmonic optical phonons that govern their photophysics. Here, we address the nature of phonon–phonon scattering probabilities of the resonantly excited optical phonons that dress the electronic transitions in these materials. Based on the temperature dependence of the coherent phonon lifetimes, we isolate the dominant anharmonic phonon and quantify its phonon–phonon interaction strength. Intriguingly, we also observe that the anharmonicity is distinct for different phonons, with a few select modes exhibiting temperature-independent coherence lifetimes, indicating their predominantly harmonic nature. However, the population and dephasing dynamics of excitons are dominated by the anharmonic phonon. 

   more » « less
4. Density functional theory study of bulk properties of transition metal nitrides

   https://doi.org/10.1039/D2CP06082E  

   Lynn, Michael O.; Ologunagba, Damilola; Dangi, Beni B.; Kattel, Shyam (February 2023, Physical Chemistry Chemical Physics)

   Density functional theory (DFT) calculations are performed to compute the lattice constants, formation energies and vacancy formation energies of transition metal nitrides (TMNs) for transition metals (TM) ranging from 3d–5d series. The results obtained using six different DFT exchange and correlation potentials (LDA, AM05, BLYP, PBE, rPBE, and PBEsol) show that the experimental lattice constants are best predicted by rPBE, while the values obtained using AM05, PBE, rPBE and PBEsol lie between the LDA and BLYP calculated values. A linear relationship is observed between the lattice constants and formation energies with the mean radii of TM and the difference in the electronegativity of TM and N in TMNs, respectively. Our calculated vacancy formation energies, in general, show that N-vacancies are more favorable than TM-vacancies in most TMNs. We observe that N-vacancy formation energies are linearly correlated with the calculated bulk formation energies indicating that TMNs with large negative formation energies are less susceptible to the formation of N-vacancies. Thus, our results from this extensive DFT study not only provide a systematic comparison of various DFT functionals in calculating the properties of TMNs but also serve as reference data for the computation-driven experimental design of materials. 

   more » « less
5. The Lieb-Oxford Lower Bounds on the Coulomb Energy, Their Importance to Electron Density Functional Theory, and a Conjectured Tight Bound on Exchange

   John P. Perdew; Jianwei Sun (July 2022, The Elliott Lieb Anniversary Volume)

   M. Lewin, Rupert L. (Ed.)

   Abstract: Lieb and Oxford (1981) derived rigorous lower bounds, in the form of local functionals of the electron density, on the indirect part of the Coulomb repulsion energy. The greatest lower bound for a given electron number N depends monotonically upon N, and the N→∞ limit is a bound for all N. These bounds have been shown to apply to the exact density functionals for the exchange- and exchange-correlation energies that must be approximated for an accurate and computationally efficient description of atoms, molecules, and solids. A tight bound on the exact exchange energy has been derived therefrom for two-electron ground states, and is conjectured to apply to all spin-unpolarized electronic ground states. Some of these and other exact constraints have been used to construct two generations of non-empirical density functionals beyond the local density approximation: the Perdew-Burke-Ernzerhof (PBE) generalized gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA. 

   more » « less