[Pw_forum] SCF convergence test - iteration
Ker Park
kerpark at hotmail.com
Mon Oct 22 04:46:43 CEST 2012
Hello all,
I am new to quantum espresso. I am trying to do the scf convergence test (ecutwfc) with silicon.I was able to run an scf calculation with a single 'ecutwfc', but I want the calculation with multiple 'ecutwfc's using iteration method (for-loop) to save me the trouble of running each calculation one by one.
Below is the input file that I used, but I've got an error saying "ecutwfc value is out of range"I don't think my cutoff energy was out of range but a syntax error.I executed the simulation with a command "~/pw.x < si.scf.in > si.scf.out"
--------------------------------------------------------------------------------------------! This code does the convergence test
for CUTOFF in 15 20docat >> conv$CUTOFF.in << EOF
&control calculation='scf' restart_mode='from_scratch', !pseudo_dir='directory where pseudopotentials are stored/', !outdir='directory where large files are written/' pseudo_dir='./', outdir='./tmp' prefix='PH', / &system ibrav=2, celldm(1)=10.20, nat=2, ntyp=1, ecutwfc =$CUTOFF / &electrons conv_thr = 1.0d-8 mixing_beta = 0.7 /ATOMIC_SPECIES Si 28.0855 Si.pz-vbc.UPFATOMIC_POSITIONS Si 0.00 0.00 0.00 Si 0.25 0.25 0.25 K_POINTS 10 0.1250000 0.1250000 0.1250000 1.00 0.1250000 0.1250000 0.3750000 3.00 0.1250000 0.1250000 0.6250000 3.00 0.1250000 0.1250000 0.8750000 3.00 0.1250000 0.3750000 0.3750000 3.00 0.1250000 0.3750000 0.6250000 6.00 0.1250000 0.3750000 0.8750000 6.00 0.1250000 0.6250000 0.6250000 3.00 0.3750000 0.3750000 0.3750000 1.00 0.3750000 0.3750000 0.6250000 3.00
EOF
pw.x < conv$CUTOFF.in > conv$CUTOFF.out----------------------------------------------------------------------------------------------
Any correction or example would be appreciated. Also, based on what language has this code been written?
Thank you.
Best,Ker
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20121021/7d3e6407/attachment.html>
More information about the users
mailing list