View source for GSFE Curve Generation Script

[[Al P6-DFT|< Back to GSFE Curve]]

=Purpose=
This BASH script will generate run [[Code: VASP]] while the shearing the FCC cell generated by a python script.

=Implementation=
* The shearing displacement is governed by this line <pre>for a in `seq -w 0.0 0.05 1.0`</pre>Edit the sequence to edit the shearing.
* Alter this line <pre>./<name>.py 4.05 1 4 1 $a p</pre> with an appropriate POSCAR genertion file name, lattice constant, periodicity, and output format.
* The formatting for this line is:<pre><name>.py <a> <nx> <ny> <nz> <t> <p/d></pre> where the arguments are:

<pre>

        Usage: fcc_111_gsfe_curve.py a nx ny nz t p/d        
               a  - equilibrium lattice parameter
               nx,ny,nz - periodicity in x,y and z
               t - displacement of top block of atoms
               p/d - output filetype p(POSCAR)/d(LAMMPS datafile)
        Example: fcc_111_gsfe_curve.py 4.05 1 4 1 0.1 p

</pre>

=Source Code=
<pre>
#!/bin/sh

rm EvsA 
nproc=$1
totT=0

for a in `seq -w 0.0 0.05 1.0`
do
echo "a= $a"

./<name>.py 4.05 1 4 1 $a p
mpiexec -np $nproc ./job.sh

E=`tail -n1 OSZICAR | awk '{ print $5}'`

echo $a $E >> EvsA

cleanvaspfiles;rm C* W*

done
</pre>



[[Category: Script]]