Bulk Modulus calculation using Quantum ESPRESSO


Quantum ESPRESSO has a package called ev.x that can be used to calculate the bulk modulus of a material. ev.x basically offers calculation of 4 equation of states, namely: These equation of states give the equilibrium volume (lattice constant in case of cubic) as well as the bulk modulus and it’s first derivative.

The procedure is pretty simple. You first need to calculate the energy of the system at different volumes of unit cell. So you would basically need to run several SCF calculations at different values of lattice parameters (thereby different volumes) and note down the corresponding energies in a text file. For cubic systems you can have the lattice parameter in the first column and the Energy in the second column, while for non-cubic systems, you’d note down the volumes in the first column.

A typical file would look like this:

Let’s say this text file is called: volvsE.txt

Now, we can use this to get the Murnaghan [1] equation of state which is very popular among the scientific community.
Murnaghan EOS is given by:

\boxed{E(V)=E(V_0)+\frac{B_0V}{B'_0}\left[\frac{\left(V_0/V\right)^{B_0'}}{B'_0-1}+1\right]-\frac{B_0V_0}{B_0'-1}  }

Now you just need to run ev.x from the command-line (terminal) and supply inputs to the following prompts:

[email protected]:~$ /home/manas/qe-6.1/bin/ev.x
Lattice parameter or Volume are in (au, Ang) > Angstrom
Assuming Angstrom
Enter type of bravais lattice (fcc, bcc, sc, noncubic) > sc
Enter type of equation of state :
1=birch1, 2=birch2, 3=keane, 4=murnaghan > 4
Input file > volvsE.txt
Output file > eos

Once, the calculation is done you should have file (eos.txt) that looks like this:

Example output file generated from ev.x

This contains information of the equilibrium volume, bulk modulus and it’s derivative and the minimum energy.

You can use the above generated data to create a neat plot of the equation of state (as you’d see in many publications), using gnuplot. The following shell script, plots the calculated data-points from file volvsE.txt as well as the fitted equation of state curve. [NOTE: This shell script is for cubic systems only. You’d need to modify the code to implement it for non-cubic systems]

Create a file called script.sh and add the following code to it:

V=$(echo "$a*$a*$a"|bc)
echo "
set terminal png size 1000,500 
set output 'EOSplot.png'
set xlabel 'Volume (Ang^3)'
set ylabel 'Energy (Ry)'
set title 'Equation of state'
#set key box linestyle 1
f(x)= "$emin"+(x*10**(-30)*"$k"*10**9/"$dk"*(("$V"/x)**"$dk"/("$dk"-1)+1)-"$k"*10**9*"$V"*10**(-30)/("$dk"-1))/(1.6*10**(-19)*13.6)
plot f(x) w l title 'Murnaghan EOS fit', 'volvsE.txt' u 1:2 w points pointtype 7 pointsize 2 title 'Calculated data-points'
set terminal postscript enhanced color solid 22
set output 'EOSplot.eps'
set xlabel 'Volume (Ang^3)'
set ylabel 'Energy (Ry)'
set title 'Equaiton of state'
set autoscale
plot f(x) w l title 'Murnaghan EOS fit', 'volvsE.txt' u 1:2 w points pointtype 7 pointsize 2 title 'Calculated data-points'
set term x11" >EOSplotScript.p

Now save the file and open the terminal and run the command:
chmod u+x script.sh

Now simply run the script by providing the four arguments, in the following order: lattice parameter, bulk modulus, derivative of BM, and minimum energy.

Example: ./script.sh 5.5 81 3.9 -588.656

The above command will create a gnuplot script file called: EOSplotScript.p

This can then be run using gnuplot using the following command:
gnuplot ./EOSplotScript.p

When you run this command, it will generate a PNG image as well as an EPS file.

The output looks like this:

Example of the PNG and EPS file generated by the above script


[1] F.D. Murnaghan, Proceedings of the National Academy of Sciences 30, 244 (1944).

I'm a physicist specializing in theoretical, computational and experimental condensed matter physics. I like to develop Physics related apps and softwares from time to time. Can code in most of the popular languages. Like to share my knowledge in Physics and applications using this Blog and a YouTube channel.

Leave a Reply

Your email address will not be published. Required fields are marked *