TransWikia.com

Temperature effect on elastic constant using VASP

Matter Modeling Asked by Niraja moharana on February 15, 2021

I want to incorporate the temperature effect on elastic constant (EC) of a system using VASP.
From literature, there are two ways which can be used to get the temperature effect on EC i) ab initio molecular dynamics ii) standard DFT with quasi harmonic approximation (as far as I know). Based on the second approach, standard DFT calculation with IBRION = 5 or 6 with ISIF = 3 tags can generate EC at 0 Kelvin for a particular volume. Again, volume vs. temperature data is one of the outputs of the phonon codes (QHA). But somehow I am not able to correlate this two quantities.

Can somebody please explain a detailed way for getting the temperature effect by considering both theories and code (VASP)?

One Answer

It is always better to write your own code to get elastic constants. For example , In case of cubic system, we need three types of distorsion to unit cell. Now elastic constant is simply slope of second order fit of energies at different value of distorsion (see different publications) normalized with Volume minima (volume corresponds to minimum energy) . To calculate elastic constants at higher temperature, we need two things

  1. Free energy of crystal at different temperature( use phonopy)
  2. Volume minima at higher temp( we can get from v-t)

Second order fit of free energy (at higher temp) normalized with volume minima of that temp will give elastic constants at higher temperature.

Answered by pranav kumar on February 15, 2021

Add your own answers!

Ask a Question

Get help from others!

© 2024 TransWikia.com. All rights reserved. Sites we Love: PCI Database, UKBizDB, Menu Kuliner, Sharing RPP