[Pw_forum] problem in xspectra calculations
Niharika Joshi
joshiniharika20 at yahoo.in
Wed Sep 28 14:59:02 CEST 2011
Hello,
I am trying the test xspectra calculation for diamond. The scf and fermi energy calculations run fine, but the xanes calculation remains incomplete. In its output file, at the end it says 'Could not find the element C in the table of K edge energies!'
I'm using quantum espresso v4.3.2
Can anyone suggest what possibly is wrong?
Below are the inputs for scf, fermi energy and xanes calculations.
&control
calculation='scf',
pseudo_dir = './Matteo/Gipaw_Pseudo_Generation/C_Hole',
outdir='./tmp/',
prefix='diamondh',
verbosity='high'
/
&system
ibrav = 1,
celldm(1) = 6.740256,
nat=8,
ntyp=2,
nbnd=17,
ecutwfc=40.0,
tot_charge=+1.0,
/
&electrons
mixing_beta = 0.3,
/
ATOMIC_SPECIES
Ch 12.0 C.star1s-pbe-mt_gipaw.UPF
C 12.0 C.pbe-mt_gipaw.UPF
ATOMIC_POSITIONS crystal
Ch 0.0 0.0 0.0
C 0.0 0.5 0.5
C 0.5 0.0 0.5
C 0.5 0.5 0.0
C 0.75 0.75 0.25
C 0.75 0.25 0.75
C 0.25 0.75 0.75
C 0.25 0.25 0.25
K_POINTS automatic
4 4 4 0 0 0
input for calculating fermi energy...
&input_xspectra
calculation='fermi_level',
prefix='diamondh',
outdir='./tmp/',
xread_wf=.true.,
/
&plot
/
&pseudos
filecore='Ch1.wfc',
/
&cut_occ
/
4 4 4 0 0 0
input for xanes....
&input_xspectra
calculation='xanes_dipole',
prefix='diamondh',
outdir='./tmp/',
xniter=1000,
xcheck_conv=50,
xepsilon(1)=1.0,
xepsilon(2)=0.0,
xepsilon(3)=0.0,
xiabs=1,
x_save_file='diamondh.xspectra.sav',
ef_r=0.94036945,
xerror=0.001,
/
&plot
xnepoint=1000,
xgamma=0.8,
xemin=-10.0,
xemax=30.0,
terminator=.true.,
cut_occ_states=.false.,
/
&pseudos
filecore='Ch1.wfc',
r_paw(1)=3.2,
/
&cut_occ
/
4 4 4 1 1 1
Thanking in advance.
-Niharika Joshi
(project student, IISER
Pune)
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20110928/6d334959/attachment.html>
More information about the users
mailing list