[Pw_forum] Electron-phonon interaction ...
Amit Kumar
amit76.india at gmail.com
Sat Dec 16 13:31:53 CET 2006
Respected Sir,
I think the error is not so simple.
Because I did all what you think a student is supposed to
do before posting a problem.
Here is my input ------>
cat > hexagonal_2D_sample.scf.fit.in << EOF
&control
calculation='scf'
restart_mode='from_scratch',
prefix='hexa',
pseudo_dir = '$PSEUDO_DIR/',
outdir='$TMP_DIR/'
tstress = .true.
tprnfor = .true.
/
&system
ibrav = 4,
celldm(1) = 4.16,
celldm(3) = 16.0,
nat = 2,
ntyp = 1,
ecutwfc =40.0,
ecutrho = 240.0,
occupations='smearing',
smearing='methfessel-paxton',
degauss=0.01,
la2F = .true.,
/
&electrons
conv_thr = 1.0d-8
mixing_beta = 0.7
ATOMIC_POSITIONS (angstrom)
A 0.000000000 0.001189495 0.000000000
A 0.000000000 1.418810505 0.000000000
K_POINTS {automatic}
32 32 4
****************************************************************************
Same input file for SCF calculation only with out la2F flag and
with
K_POINTS {automatic}
8 8 1
*****************************************************************************
cat > hexagonal_2D_sample.elph.in << EOF
Electron-phonon coefficients for 2D hexagoanl sample
&inputph
tr2_ph=1.0d-10,
prefix='hexa',
fildvscf='hexadv',
amass(1)=58.9332,
outdir='$TMP_DIR/',
fildyn='hexa.dyn',
elph=.true.,
trans=.true.,
ldisp=.true.
nq1=4, nq2=4, nq3=1
/
EOF
$ECHO " running the el-ph calculation...\c"
*******************************************************************************************************
I have checked for several combinations of K_POINTS and nq1, nq2 and nq3
also.
But I' still getting the same error.
Looking forward to your kind help.
Regards,
Amit
As for the answer to the question below and similar ones
(e.g. when the error is followed by a message from the code)
I can advice the first debugging/input_check step:
-- goto the code and find sobroutine which gives the error
or grep the error message in the source by making for instance
grep -i "cannot remap grid on k-point list" ./*/*.f90 (or PH/*.f
90)
-- if you find that place, you may see "a hint"
-- You think what are "k1,k2,k3" in the code. You find that you sh
ould
specify these numbers as zero in the input (otherwise there
are problems with the symmetry... because somebody -- read Gosi
a
in this case -- was lazy to solve it).
-- You change your input, and if it still doesn't work.... you sen
d us
your input..
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20061216/bd3c4c98/attachment.html>
More information about the users
mailing list