skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.

Explore Research Products in the PAR
It may take a few hours for recently added research products to appear in PAR search results.


  • NSF Public Access
  • Search Results
  • Enhancing efficiency and scope of first-principles quasiharmonic approximation methods through the calculation of third-order elastic constants
Title: Enhancing efficiency and scope of first-principles quasiharmonic approximation methods through the calculation of third-order elastic constants
A novel computational strategy is presented to calculate from first principles the coefficient of thermal expansion and the elastic constants of a material over meaningful intervals of temperature and pressure. This strategy combines a novel implementation of the quasiharmonic approximation to calculate the isothermal-isochoric linear and nonlinear elastic constants of a material, with elementary equations of nonlinear continuum mechanics. Our implementation of the quasiharmonic approximation relies on finite deformations, the use of nonprimitive supercells to describe a material, a recently proposed technique to calculate generalized mode Grüneisen parameters, and the numerical differentiation of the stress tensor to calculate both second- and third-order elastic constants. The combination of this method with nonlinear continuum mechanics is shown to yield accurate predictions of lattice parameters and linear elastic constants of a material over finite intervals of temperature and pressure, at the cost of calculating isothermal second- and third-order elastic constants for a single reference state. Here, the validity and limits of our novel methods are assessed by carrying out calculations of MgO based on classical interatomic potentials. To demonstrate potential, our methods are then used in conjunction with a density functional theory approach to calculate thermal expansion and elastic properties of silicon, lithium hydrate, and graphite.  more » « less
Award ID(s):
2036176
PAR ID:
10337842
Author(s) / Creator(s):
;
Date Published:
Journal Name:
Physical review materials
Volume:
6
Issue:
4
ISSN:
2476-0455
Page Range / eLocation ID:
043803
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. ; ; ; ; (, High Pressure Research)
    null (Ed.)
    Thermoelastic behavior of transition metal boride Os2B3 was studied under quasi-hydrostatic and isothermal conditions in a Paris-Edinburgh cell to 5.4 GPa and 1273 K. In-situ Energy Dispersive X-ray diffraction was used to determine interplanar spacings of the hexagonal crystal structure and thus the volume and axial compression. P-V-T data were fitted to a 3rd Order Birch-Murnaghan equation of state with a temperature modification to determine thermal elastic constants. The bulk modulus was shown to be K0 = 402 ± 21 GPa when the first pressure derivative was held to K0’ = 4.0 from the room temperature P-V curve. Under a quadratic fit α=α_0+α_1 T-α_2 T^(-2), the thermal expansion coefficients were determined to be α_0=1.862×10^(-5) K-1, α_1=0.841×10^(-9) K-2, and α_2=-0.525 K. Density functional theory (DFT) with the quasi-harmonic approximation (QHA) were further employed to study Os2B3, including its P-V-T curves, phonon spectra, bulk modulus, specific heat, thermal expansion, and the Grüneisen parameter. A good agreement between the first-principle theory and experimental observations was achieved, highlighting the success of the Armiento-Mattsson 2005 generalized gradient approximation functional employed in this study and QHA for describing thermodynamic properties of Os2B3. 
    more » « less
  2. ; ; ; ; ; ; ; ; (, 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. ; ; ; (, Physical Review Materials)
    While anomalous diffusion coefficients with non-Arrhenius-like temperature dependence are observed in a number of metals, a conclusive comprehensive framework of explanation has not been brought forward to date. Here, we use first-principles calculations based on density functional theory to calculate self-diffusion coefficients in the bcc metals Mo and β-Ti by coupling quasiharmonic transition state theory and large-displacement phonon calculations and show that anharmonicity from thermal expansion is the major reason for the anomalous temperature dependence. We use a modified Debye approach to quantify the thermal expansion over the entire temperature range and introduce a method to relax the vacancy structure in a mechanically unstable crystal such as β-Ti. The effect of thermal expansion is found to be crucial for the nonlinear, non-Arrhenius “anomalous” self-diffusion in both bcc systems, with β-Ti showing a 60% larger relative nonlinearity parameter than Mo. Our results point to temperature dependence in the diffusion prefactor from thermal expansion as the major origin of anomalous self-diffusion. The methodology proposed for β-Ti is general and simple enough to be applicable to other mechanically unstable crystals. 
    more » « less
  4. ; ; (, Review of Scientific Instruments)
    Time-domain thermoreflectance and frequency-domain thermoreflectance (FDTR) have been widely used for non-contact measurement of anisotropic thermal conductivity of materials with high spatial resolution. However, the requirement of a high thermoreflectance coefficient restricts the choice of metal coating and laser wavelength. The accuracy of the measurement is often limited by the high sensitivity to the radii of the laser beams. We describe an alternative frequency-domain pump-probe technique based on probe beam deflection. The beam deflection is primarily caused by thermoelastic deformation of the sample surface, with a magnitude determined by the thermal expansion coefficient of the bulk material to measure. We derive an analytical solution to the coupled elasticity and heat diffusion equations for periodic heating of a multilayer sample with anisotropic elastic constants, thermal conductivity, and thermal expansion coefficients. In most cases, a simplified model can reliably describe the frequency dependence of the beam deflection signal without knowledge of the elastic constants and thermal expansion coefficients of the material. The magnitude of the probe beam deflection signal is larger than the maximum magnitude achievable by thermoreflectance detection of surface temperatures if the thermal expansion coefficient is greater than 5 × 10 −6  K −1 . The uncertainty propagated from laser beam radii is smaller than that in FDTR when using a large beam offset. We find a nearly perfect matching of the measured signal and model prediction, and measure thermal conductivities within 6% of accepted values for materials spanning the range of polymers to gold, 0.1–300 W/(m K). 
    more » « less
  5. ; ; ; (, New Journal of Physics)
    Abstract Parametric amplification of an elastic wave and a framework for using elastic waves that could enable a new generation of high performance, low noise acoustic amplifiers, mixers and circulators are presented. Using a novel approach with nonlinear materials produces highly desirable non-reciprocal characteristics. Parametric amplification of a weak elastic signal wave is achieved by an elastic pump wave of higher intensity. By careful selection of material orientation together with precise excitation of signal and pump waves, ‘up frequency conversion’ is suppressed and selective amplification of the elastic signal wave occurs at its original frequency. In addition, a general mathematical framework is developed and used for analytical studies of coupled wave equations in nonlinear anisotropic materials. The results obtained from the analytical studies are verified using a finite element implementation. 
    more » « less
  • Citation Formats
  • MLA
  • APA
  • Chicago
  • BibTeX
  • Save / Share this Record
  • Send to Email