<div dir="ltr">Hi everyone!<div>I'm new in Quantum-Espresso and i having problems to minimize que lattice parameter of gold</div><div>The minimum energy is with 7.0A and in literature we find 4.07 A and in other programs i find 4.2</div>
<div>I'm using this code to calculate the curve</div><div><br></div><div><div>#!/bin/sh</div><div>####################################################################</div><div>#</div><div># define the following variables according to your needs</div>
<div>#</div><div>#espresso_dir=top_directory_of_espresso_package</div><div>outdir=/home'</div><div>pseudo_dir=/home/pseudo'</div><div>####################################################################</div><div>
<br></div><div>rm -f au.scf.out au.etot_vs_alat</div><div>touch au.etot_vs_alat</div><div><br></div><div>for alat in  4.05 4.06 4.07 4.08 4.09 4.1 4.2 5.0 6.0 7.0 8.0 9.0 ; do</div><div><br></div><div># self-consistent calculation</div>
<div>cat > <a href="http://au.scf.in">au.scf.in</a> << EOF</div><div>&control</div><div>    calculation = 'scf'</div><div>    prefix='gold',</div><div>    pseudo_dir = pseudo_dir,</div><div>    outdir= outdir,</div>
<div> /</div><div> &system</div><div>    ibrav=  2, a=$alat , nat=  4, ntyp= 1,</div><div>    ecutwfc = 40.0, ecutrho = 480.0,</div><div>    occupations='smearing', smearing = 'mv', degauss = 0.002</div>
<div> /</div><div> &electrons</div><div>    mixing_mode = 'local-TF'</div><div>    mixing_beta = 0.7</div><div>    conv_thr = 1.0d-6</div><div> /</div><div>ATOMIC_SPECIES</div><div> Au  196.96654  Au.bp-van_ak.UPF</div>
<div>ATOMIC_POSITIONS (crystal)</div><div> Au 0.00 0.00 0.00</div><div> Au 0.5  0.5  0.0</div><div> Au 0.5  0.0  0.5</div><div> Au 0.0  0.5  0.5</div><div>K_POINTS automatic</div><div>   8 8 8 1 1 1</div><div><br></div><div>
EOF</div><div><br></div><div>pw.x < <a href="http://au.scf.in">au.scf.in</a> > au.scf.out</div><div><br></div><div>grep -e 'lattice parameter' -e ! au.scf.out | \</div><div>      awk '/lattice/{a=$(NF-1)}/!/{print a, $(NF-1)}' >> au.etot_vs_alat</div>
<div><br></div><div>done</div><div><br></div><div>I find this in a tutorial of Professor Narasimhan. 7.0 A i think is large for gold, so i try to find for Silver and the alat was almost the same. I tested all the PP for gold that i find in the QE page, and has the same result with different energies values. If some one could help i will be grateful.</div>
<div><br></div><div><br></div>-- <br>Aline A. Bidoul<div>Engenheira Física</div><div>Mestranda em Engenharia Biomédica</div>
</div></div>