[Pw_forum] matdyn.x warning message
Mikiyas Tsegaye
mastermik at gmail.com
Thu Aug 20 21:03:37 CEST 2009
Dear Friends,
I'm trying to plot the phonon dispersion of graphene. I have done the
dispersions for Silicon, and Germanium, and also for Diamond. All of
which worked out well. Graphene, however, is giving me the following
error in the second-to-last matdyn.x step:
--------------------------------------------------------------------
message from routine matdyn
Z* not found in file graphene.fc, TO-LO splitting at q=0 will be absent
A direction for q was not specified: TO-LO splitting will be absent
--------------------------------------------------------------------
The dispersion that it then produces looks very jagged and the numbers
are off. For examples, the acoustic modes don't go to zero at the
gamma point, and the acoustic and optical branches don't touch at the
K-point like they're supposed to. I would really appreciate any help
on this matter. My input file is as follows:
-----------------------------------------------
# SCF step
&control
calculation='scf',
restart_mode='from_scratch',
prefix='graph',
tstress=.true.,
tprnfor=.true.,
pseudo_dir = '$PSEUDO_DIR/',
outdir='$TMP_DIR/'
/
&system
ibrav = 0, nat=2, ntyp= 1,
ecutwfc = 40.0,
ecutrho = 480.0,
occupations = 'smearing',
smearing = 'gaussian',
degauss = 0.02,
/
&electrons
diagonalization = 'cg',
conv_thr = 1.0d-6
/
CELL_PARAMETERS {hexagonal}
4.02 2.321 0.00
4.02 -2.321 0.00
0.00 0.00 69.75
ATOMIC_SPECIES
C 12.01 C.pz-vbc.UPF
ATOMIC_POSITIONS {bohr}
C 0.0 0.0 0.0
C 2.68 0.0 0.0
K_POINTS automatic
16 16 1 0 0 0
$PW_COMMAND < graphene.scf.in > graphene.scf.out
# calculate phonons
phonons of graphene
&INPUTPH
amass(1) = 12.01,
outdir = '$TMP_DIR/',
prefix = 'graph',
tr2_ph = 1.0d-12,
fildyn = 'graph.dyn',
fildrho = 'graph.rho',
ldisp = .true.,
lnscf = .true.,
nq1 = 16,
nq2 = 16,
nq3 = 1
-------------------------------------------------------
Am I missing something in my inputs?
thank you!
- Mikiyas Tsegaye
University of Virginia, ECE
Charlottesville, VA, USA
mst3c at virginia.edu
More information about the Pw_forum
mailing list