[Pw_forum] error: npert > max_irr_dim
Aritz Leonardo
swblelia at sw.ehu.es
Wed Feb 22 10:14:22 CET 2006
Hello:
This is the job that crashes. From a list of Q vectors
q=0.30,0.211695,0.0 gives the mentioned error
I am running the job in a parallel cluster of OPTERONs with portlan 6.0
compiler.
I repeat that the same calculation with E=32 ended ok!
Thanks a lot.
E=35.0
ca=6.5
a=4.65
# %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%5
for PARAM in `cat qELPH1`;
do
QX="`echo $PARAM | awk -F: '{print $1}'`";
QY="`echo $PARAM | awk -F: '{print $2}'`";
QZ="`echo $PARAM | awk -F: '{print $3}'`";
echo "$QX , $QY , $QZ ,$a";
#
# self-consistent calculation
#
cat > $a.scf.in << EOF
&control
calculation='scf',
title='C MONOLAYER surface phonons'
restart_mode='from_scratch',
prefix='C',
pseudo_dir = './',
outdir='$TMP_DIR/',
tprnfor = .true.
tstress = .true.
/
&system
ibrav=4,celldm(1)=$a celldm(3)=$ca,nat=2,ntyp= 1,
ecutwfc=$E,occupations='smearing',
smearing='m-p',degauss=0.0075
/
&electrons
conv_thr = 1.0d-8
mixing_beta= 0.5
/
ATOMIC_SPECIES
C 12.01 C.pz-vbc.UPF
ATOMIC_POSITIONS
C 0.0000 0.00000000 0.000
C 0.0000 0.577350269 0.000
K_POINTS {automatic}
20 20 1 1 1 0
EOF
#
# non self-consistent calculation
#
cat > $a.nscf.in << EOF
&control
calculation='phonon'
restart_mode='from_scratch',
pseudo_dir = './',
outdir='$TMP_DIR/',
prefix='C'
/
&system
ibrav=4,celldm(1)=$a, celldm(3)=$ca,nat=2,ntyp=1,
ecutwfc=$E,occupations='smearing',
smearing='m-p',degauss=0.0075
/
&electrons
conv_thr = 1.0d-8
mixing_beta = 0.5
/
&phonon
xqq(1)=$QX, xqq(2)=$QY, xqq(3)=$QZ
/
ATOMIC_SPECIES
C 12.01 C.pz-vbc.UPF
ATOMIC_POSITIONS
C 0.0000 0.000000000 0.000
C 0.0000 0.577350269 0.000
K_POINTS {automatic}
20 20 1 1 1 0
EOF
#
# phonon calculation
#
cat > $a.ph.in << EOF
ESTA LINEA ES NECESARIA
&inputph
alpha_mix(1)=0.2
tr2_ph=1.0d-13,
prefix='C',
amass(1)=12.01,
outdir='$TMP_DIR/'
fildyn='C.dyn'
fildvscf='Cdv'
/
$QX $QY $QZ
EOF
More information about the Pw_forum
mailing list