problem with phonon calculations
Posted: Mon Aug 06, 2012 2:20 pm
Someone please help me out. These are my input files,
# Crystalline CaB6 : computation of the phonon spectrum
ndtset 12
#Set 1 : ground state self-consistency
getwfk1 0 # Cancel default
kptopt1 1 # Automatic generation of k points, taking
# into account the symmetry
nqpt1 0 # Cancel default
tolvrs1 1.0d-18 # SCF stopping criterion (modify default)
rfphon1 0 # Cancel default
#Q vectors for all datasets
#Complete set of symmetry-inequivalent qpt chosen to be commensurate
# with kpt mesh so that only one set of GS wave functions is needed.
#Generated automatically by running GS calculation with kptopt=1,
# nshift=0, shiftk=0 0 0 (to include gamma) and taking output kpt set
# file as qpt set. Set nstep=1 so only one iteration runs.
nqpt 1 # One qpt for each dataset (only 0 or 1 allowed)
# This is the default for all datasets and must
# be explicitly turned off for dataset 1.
qpt2 0.00000000E+00 0.00000000E+00 0.00000000E+00
qpt3 0.00000000E+00 0.00000000E+00 0.00000000E+00
qpt4 2.50000000E-01 0.00000000E+00 0.00000000E+00
qpt5 5.00000000E-01 0.00000000E+00 0.00000000E+00
qpt6 2.50000000E-01 2.50000000E-01 0.00000000E+00
qpt7 5.00000000E-01 2.50000000E-01 0.00000000E+00
qpt8 5.00000000E-01 5.00000000E-01 0.00000000E+00
qpt9 2.50000000E-01 2.50000000E-01 2.50000000E-01
qpt10 5.00000000E-01 2.50000000E-01 2.50000000E-01
qpt11 5.00000000E-01 5.00000000E-01 2.50000000E-01
qpt12 5.00000000E-01 5.00000000E-01 5.00000000E-01
#Set 2 : Response function calculation of d/dk wave function
iscf2 -3 # Need this non-self-consistent option for d/dk
kptopt2 2 # Modify default to use time-reversal symmetry
rfphon2 0 # Cancel default
rfelfd2 2 # Calculate d/dk wave function only
tolwfr2 1.0d-22 # Use wave function residual criterion instead
#Set 3 : Response function calculation of Q=0 phonons and electric field pert.
getddk3 2 # d/dk wave functions from last dataset
kptopt3 2 # Modify default to use time-reversal symmetry
rfelfd3 3 # Electric-field perturbation response only
#Sets 4-12 : Finite-wave-vector phonon calculations (defaults for all datasets)
getwfk 1 # Use GS wave functions from dataset1
kptopt 3 # Need full k-point set for finite-Q response
rfphon 1 # Do phonon response
rfatpol 1 7 # Treat displacements of all atoms
rfdir 1 1 1 # Do all directions (symmetry will be used)
tolvrs 1.0d-8 # This default is active for sets 3-10
#########################################################################################################
#Data common to all datasets
nshiftk 1
shiftk 0.5 0.5 0.5
ngkpt 4 4 4 ! This is much too low : should be at least 24x24x24
acell 3*4.19
diemac 10.0
ecut 45.00 ! This is also too low
iscf 5
nstep 25
ixc 11
natom 7
ntypat 2
rprim 1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0
# Atomic positions
xred 0.00 0.00 0.00 #Ca1
0.50 0.50 0.20 #B2
0.50 0.50 -0.20 #B2
0.50 0.20 0.50 #B2
0.50 -0.20 0.50 #B2
0.20 0.50 0.50 #B2
0.80 0.50 0.50 #B2
nband 16
typat 1 2 2 2 2 2 2
znucl 20 5
# add to conserve old < 6.7.2 behavior for calculating forces at each SCF step
optforces 1
And
!Input file for the anaddb code. Analysis of the CaB6 DDB
!Flags
ifcflag 1 ! Interatomic force constant flag
!Wavevector grid number 1 (coarse grid, from DDB)
brav 1 ! Bravais Lattice : 1-S.C., 2-F.C., 3-B.C., 4-Hex.)
ngqpt 4 4 4 ! Monkhorst-Pack indices
nqshft 1 ! number of q-points in repeated basic q-cell
q1shft 3*0.0
!Effective charges
asr 1 ! Acoustic Sum Rule. 1 => imposed asymetrically
chneut 1 ! Charge neutrality requirement for effective charges.
!Interatomic force constant info
dipdip 1 ! Dipole-dipole interaction treatment
!Phonon band structure output for band2eps - See note near end for
! dealing with gamma LO-TO splitting issue.
eivec 4
!Wavevector list number 1 (Reduced coordinates and normalization factor)
nph1l 51 ! number of phonons in list 1
qph1l 0.000 0.000 0.000 1.0 !(gamma point)
0.000 0.050 0.000 1.0
0.000 0.100 0.000 1.0
0.000 0.150 0.000 1.0
0.000 0.200 0.000 1.0
0.000 0.250 0.000 1.0
0.000 0.300 0.000 1.0
0.000 0.350 0.000 1.0
0.000 0.400 0.000 1.0
0.000 0.450 0.000 1.0
0.000 0.500 0.000 1.0 !(X point)
0.050 0.050 0.000 1.0
0.100 0.100 0.000 1.0
0.150 0.150 0.000 1.0
0.200 0.200 0.000 1.0
0.250 0.250 0.000 1.0
0.300 0.300 0.000 1.0
0.350 0.350 0.000 1.0
0.400 0.400 0.000 1.0
0.450 0.450 0.000 1.0
0.500 0.500 0.000 1.0 !(M point)
0.450 0.450 0.000 1.0
0.400 0.400 0.000 1.0
0.350 0.350 0.000 1.0
0.300 0.300 0.000 1.0
0.250 0.250 0.000 1.0
0.200 0.200 0.000 1.0
0.150 0.150 0.000 1.0
0.100 0.100 0.000 1.0
0.050 0.050 0.000 1.0
0.000 0.000 0.000 1.0 !(GAMMA point)
0.050 0.050 0.050 1.0
0.100 0.100 0.100 1.0
0.150 0.150 0.150 1.0
0.200 0.200 0.200 1.0
0.250 0.250 0.250 1.0
0.300 0.300 0.300 1.0
0.350 0.350 0.350 1.0
0.400 0.400 0.400 1.0
0.450 0.450 0.450 1.0
0.500 0.500 0.500 1.0 !(R point)
0.450 0.500 0.450 1.0
0.400 0.500 0.400 1.0
0.350 0.500 0.350 1.0
0.300 0.500 0.300 1.0
0.250 0.500 0.250 1.0
0.200 0.500 0.200 1.0
0.150 0.500 0.150 1.0
0.100 0.500 0.100 1.0
0.050 0.500 0.050 1.0
0.000 0.500 0.000 1.0 !(X point)
!Wavevector list number 2 (Cartesian directions for non-analytic gamma phonons)
!The output for this calculation must be cut-and-pasted into the
! t59_out.freq file to be used as band2eps input to get proper LO-TO
! splitting at gamma. Note that gamma occurrs twice.
nph2l 1 ! number of directions in list 2
qph2l 1.0 0.0 0.0 0.0
# This line added when defaults were changed (v5.3) to keep the previous, old behaviour
symdynmat 0
Now my frequency file has 21 columns, but the in the abinit example has 6 columns. Someone should please correct my input files. Thanks.
# Crystalline CaB6 : computation of the phonon spectrum
ndtset 12
#Set 1 : ground state self-consistency
getwfk1 0 # Cancel default
kptopt1 1 # Automatic generation of k points, taking
# into account the symmetry
nqpt1 0 # Cancel default
tolvrs1 1.0d-18 # SCF stopping criterion (modify default)
rfphon1 0 # Cancel default
#Q vectors for all datasets
#Complete set of symmetry-inequivalent qpt chosen to be commensurate
# with kpt mesh so that only one set of GS wave functions is needed.
#Generated automatically by running GS calculation with kptopt=1,
# nshift=0, shiftk=0 0 0 (to include gamma) and taking output kpt set
# file as qpt set. Set nstep=1 so only one iteration runs.
nqpt 1 # One qpt for each dataset (only 0 or 1 allowed)
# This is the default for all datasets and must
# be explicitly turned off for dataset 1.
qpt2 0.00000000E+00 0.00000000E+00 0.00000000E+00
qpt3 0.00000000E+00 0.00000000E+00 0.00000000E+00
qpt4 2.50000000E-01 0.00000000E+00 0.00000000E+00
qpt5 5.00000000E-01 0.00000000E+00 0.00000000E+00
qpt6 2.50000000E-01 2.50000000E-01 0.00000000E+00
qpt7 5.00000000E-01 2.50000000E-01 0.00000000E+00
qpt8 5.00000000E-01 5.00000000E-01 0.00000000E+00
qpt9 2.50000000E-01 2.50000000E-01 2.50000000E-01
qpt10 5.00000000E-01 2.50000000E-01 2.50000000E-01
qpt11 5.00000000E-01 5.00000000E-01 2.50000000E-01
qpt12 5.00000000E-01 5.00000000E-01 5.00000000E-01
#Set 2 : Response function calculation of d/dk wave function
iscf2 -3 # Need this non-self-consistent option for d/dk
kptopt2 2 # Modify default to use time-reversal symmetry
rfphon2 0 # Cancel default
rfelfd2 2 # Calculate d/dk wave function only
tolwfr2 1.0d-22 # Use wave function residual criterion instead
#Set 3 : Response function calculation of Q=0 phonons and electric field pert.
getddk3 2 # d/dk wave functions from last dataset
kptopt3 2 # Modify default to use time-reversal symmetry
rfelfd3 3 # Electric-field perturbation response only
#Sets 4-12 : Finite-wave-vector phonon calculations (defaults for all datasets)
getwfk 1 # Use GS wave functions from dataset1
kptopt 3 # Need full k-point set for finite-Q response
rfphon 1 # Do phonon response
rfatpol 1 7 # Treat displacements of all atoms
rfdir 1 1 1 # Do all directions (symmetry will be used)
tolvrs 1.0d-8 # This default is active for sets 3-10
#########################################################################################################
#Data common to all datasets
nshiftk 1
shiftk 0.5 0.5 0.5
ngkpt 4 4 4 ! This is much too low : should be at least 24x24x24
acell 3*4.19
diemac 10.0
ecut 45.00 ! This is also too low
iscf 5
nstep 25
ixc 11
natom 7
ntypat 2
rprim 1.0 0.0 0.0
0.0 1.0 0.0
0.0 0.0 1.0
# Atomic positions
xred 0.00 0.00 0.00 #Ca1
0.50 0.50 0.20 #B2
0.50 0.50 -0.20 #B2
0.50 0.20 0.50 #B2
0.50 -0.20 0.50 #B2
0.20 0.50 0.50 #B2
0.80 0.50 0.50 #B2
nband 16
typat 1 2 2 2 2 2 2
znucl 20 5
# add to conserve old < 6.7.2 behavior for calculating forces at each SCF step
optforces 1
And
!Input file for the anaddb code. Analysis of the CaB6 DDB
!Flags
ifcflag 1 ! Interatomic force constant flag
!Wavevector grid number 1 (coarse grid, from DDB)
brav 1 ! Bravais Lattice : 1-S.C., 2-F.C., 3-B.C., 4-Hex.)
ngqpt 4 4 4 ! Monkhorst-Pack indices
nqshft 1 ! number of q-points in repeated basic q-cell
q1shft 3*0.0
!Effective charges
asr 1 ! Acoustic Sum Rule. 1 => imposed asymetrically
chneut 1 ! Charge neutrality requirement for effective charges.
!Interatomic force constant info
dipdip 1 ! Dipole-dipole interaction treatment
!Phonon band structure output for band2eps - See note near end for
! dealing with gamma LO-TO splitting issue.
eivec 4
!Wavevector list number 1 (Reduced coordinates and normalization factor)
nph1l 51 ! number of phonons in list 1
qph1l 0.000 0.000 0.000 1.0 !(gamma point)
0.000 0.050 0.000 1.0
0.000 0.100 0.000 1.0
0.000 0.150 0.000 1.0
0.000 0.200 0.000 1.0
0.000 0.250 0.000 1.0
0.000 0.300 0.000 1.0
0.000 0.350 0.000 1.0
0.000 0.400 0.000 1.0
0.000 0.450 0.000 1.0
0.000 0.500 0.000 1.0 !(X point)
0.050 0.050 0.000 1.0
0.100 0.100 0.000 1.0
0.150 0.150 0.000 1.0
0.200 0.200 0.000 1.0
0.250 0.250 0.000 1.0
0.300 0.300 0.000 1.0
0.350 0.350 0.000 1.0
0.400 0.400 0.000 1.0
0.450 0.450 0.000 1.0
0.500 0.500 0.000 1.0 !(M point)
0.450 0.450 0.000 1.0
0.400 0.400 0.000 1.0
0.350 0.350 0.000 1.0
0.300 0.300 0.000 1.0
0.250 0.250 0.000 1.0
0.200 0.200 0.000 1.0
0.150 0.150 0.000 1.0
0.100 0.100 0.000 1.0
0.050 0.050 0.000 1.0
0.000 0.000 0.000 1.0 !(GAMMA point)
0.050 0.050 0.050 1.0
0.100 0.100 0.100 1.0
0.150 0.150 0.150 1.0
0.200 0.200 0.200 1.0
0.250 0.250 0.250 1.0
0.300 0.300 0.300 1.0
0.350 0.350 0.350 1.0
0.400 0.400 0.400 1.0
0.450 0.450 0.450 1.0
0.500 0.500 0.500 1.0 !(R point)
0.450 0.500 0.450 1.0
0.400 0.500 0.400 1.0
0.350 0.500 0.350 1.0
0.300 0.500 0.300 1.0
0.250 0.500 0.250 1.0
0.200 0.500 0.200 1.0
0.150 0.500 0.150 1.0
0.100 0.500 0.100 1.0
0.050 0.500 0.050 1.0
0.000 0.500 0.000 1.0 !(X point)
!Wavevector list number 2 (Cartesian directions for non-analytic gamma phonons)
!The output for this calculation must be cut-and-pasted into the
! t59_out.freq file to be used as band2eps input to get proper LO-TO
! splitting at gamma. Note that gamma occurrs twice.
nph2l 1 ! number of directions in list 2
qph2l 1.0 0.0 0.0 0.0
# This line added when defaults were changed (v5.3) to keep the previous, old behaviour
symdynmat 0
Now my frequency file has 21 columns, but the in the abinit example has 6 columns. Someone should please correct my input files. Thanks.