ABINIT Discussion Forums

Dear abinit users,

I have been trying to calculated phonon frequencies, band structures and eventually thermal properties of quartz SiO2 using the parallel version of abinit. During the step where the generation of DDB's is taking place my calculation fails. I do not know whether it has to do with the way the parallel environment is set up or the input file itself. Here I am attaching my input file. The output file is too heavy to attach so I will only copy and paste the end of it. Any input is very appreciated!!
In the output file the MPI setting show the following:
=== MPI ===
Parallel build : yes
Parallel I/O : yes
Time tracing : no
GPU support : no

Should parallel I/O be activated or not?

Also, the following warning :

distrb2: WARNING -
nproc_kpt= 16 >= nkpt= 10* nsppol= 1
The number of processors is larger than nkpt. This is a waste.
invars1: mkmem undefined in the input file. Use default mkmem = nkpt
invars1: With nkpt_me= 1 and mkmem = 10, ground state wf handled in core.
Resetting mkmem to nkpt_me to save memory space.

Here is the end of the output file:

--------------------------------------------------------------------------------
Perturbation wavevector (in red.coord.) 0.000000 0.000000 0.000000
Perturbation : displacement of atom 5 along direction 1
loper3 2 nkpt 16
insy3 : only one element in the set of symmetries for this perturbation :
1 0 0 0 1 0 0 0 1
symatm: atom number 1 is reached starting at atom
1
symatm: atom number 2 is reached starting at atom
2
symatm: atom number 3 is reached starting at atom
3
symatm: atom number 4 is reached starting at atom
4
symatm: atom number 5 is reached starting at atom
5
symatm: atom number 6 is reached starting at atom
6
symatm: atom number 7 is reached starting at atom
7
symatm: atom number 8 is reached starting at atom
8
symatm: atom number 9 is reached starting at atom
9
symkpt : not enough symmetry to change the number of k points.
getmpw: optimal value of mpw= 4402
kpgio: loop on k-points done in parallel
inwffil : before hdr_check
inwffil : examine the header of disk file trf2_1o_DS1_WFK
-P-0000
-P-0000 ================================================================================
-P-0000
-P-0000 - hdr_check: checking restart file header for consistency -
-P-0000
-P-0000
-P-0000 current calculation restart file
-P-0000 ------------------- ------------
-P-0000
-P-0000 calculation expects a wf_planewave | input file contains a wf_planewave
-P-0000 . ABINIT code version 6.12.3 | ABINIT code version 6.12.3
-P-0000 . date 20120719 bantot 416 natom 9 | date 20120719 bantot 260 natom 9
-P-0000 nkpt 16 nsym 2 ngfft 45, 45, 48 | nkpt 10 nsym 2 ngfft 45, 45, 48
-P-0000 ntypat 2 ecut_eff 24.0000000 | ntypat 2 ecut_eff 24.0000000
-P-0000 usepaw 0 | usepaw 0
-P-0000 usewvl 0 | usewvl 0
-P-0000 rprimd: | rprimd:
-P-0000 9.3824902 0.0000000 0.0000000 | 9.3824902 0.0000000 0.0000000
-P-0000 -4.6912451 8.1254749 0.0000000 | -4.6912451 8.1254749 0.0000000
-P-0000 0.0000000 0.0000000 10.2498745 | 0.0000000 0.0000000 10.2498745
-P-0000
-P-0000 hdr_check: WARNING -
-P-0000 input nkpt= 16 not equal disk file nkpt= 10
-P-0000 symafm: | symafm:
-P-0000 1 1 | 1 1
-P-0000 symrel: | symrel:
-P-0000 1 0 0 0 1 0 0 0 1 | 1 0 0 0 1 0 0 0 1
-P-0000 0 1 0 1 0 0 0 0 -1 | 0 1 0 1 0 0 0 0 -1
-P-0000 typat: | typat:
-P-0000 1 1 1 2 2 2 2 2 2 | 1 1 1 2 2 2 2 2 2
-P-0000 so_psp : | so_psp :
-P-0000 1 1 | 1 1
-P-0000 tnons: | tnons:
-P-0000 0.0000000 0.0000000 0.0000000 | 0.0000000 0.0000000 0.0000000
-P-0000 -0.0000000 0.0000000 -0.0000000 | -0.0000000 0.0000000 -0.0000000
-P-0000 znucl: | znucl:
-P-0000 14.00 8.00 | 14.00 8.00
-P-0000 pseudopotential atom type 1: | pseudopotential atom type 1:
-P-0000 pspso 0 pspxc 1 | pspso 0 pspxc 1
-P-0000 pspdat 940714 pspcod 1 zion 4.0 | pspdat 940714 pspcod 1 zion 4.0
-P-0000 pseudopotential atom type 2: | pseudopotential atom type 2:
-P-0000 pspso 0 pspxc 1 | pspso 0 pspxc 1
-P-0000 pspdat 940714 pspcod 1 zion 6.0 | pspdat 940714 pspcod 1 zion 6.0
-P-0000 xred: | xred:
-P-0000 0.5248000 0.0000000 0.3333330 | 0.5248000 0.0000000 0.3333330
-P-0000 -0.0000000 0.5248000 0.6666670 | -0.0000000 0.5248000 0.6666670
-P-0000 0.4752000 0.4752000 0.0000000 | 0.4752000 0.4752000 0.0000000
-P-0000 0.1570000 0.4160000 0.8768000 | 0.1570000 0.4160000 0.8768000
-P-0000 0.5840000 0.7410000 0.2101330 | 0.5840000 0.7410000 0.2101330
-P-0000 0.2590000 0.8430000 0.5434670 | 0.2590000 0.8430000 0.5434670
-P-0000 0.4160000 0.1570000 0.1232000 | 0.4160000 0.1570000 0.1232000
-P-0000 0.7410000 0.5840000 0.7898670 | 0.7410000 0.5840000 0.7898670
-P-0000 0.8430000 0.2590000 0.4565330 | 0.8430000 0.2590000 0.4565330
-P-0000
-P-0000 hdr_check: WARNING -
-P-0000 Restart of self-consistent calculation need translated wavefunctions.
-P-0000 Indeed, critical differences between current calculation and
-P-0000 restart file have been detected in:
-P-0000 * the number, position, or weight of k-points
-P-0000 * the format of wavefunctions (istwfk)
-P-0000 ================================================================================
-inwffil : will read wavefunctions from disk file trf2_1o_DS1_WFK
-P-0000 initwf : disk file gives npw= 4390 nband= 26 for kpt number= 1
-P-0000 initwf : 26 bands have been initialized from disk
-P-0000 leave_test : synchronization done...
wfsinp: loop on k-points and spins done in parallel
pareigocc : MPI_ALLREDUCE
-P-0000 leave_test : synchronization done...
wfsinp: loop on k-points done in parallel
-P-0000 - newkpt: read input wf with ikpt,npw= 1 4390, make ikpt,npw= 1 4390
-P-0000 leave_test : synchronization done...
newkpt: loop on k-points done in parallel
pareigocc : MPI_ALLREDUCE

getmpw: optimal value of mpw= 4402
kpgio: loop on k-points done in parallel
inwffil : before hdr_check
inwffil : examine the header of disk file trf2_1o_DS1_WFK
-P-0000
-P-0000 ================================================================================
-P-0000
-P-0000 - hdr_check: checking restart file header for consistency -
-P-0000
-P-0000
-P-0000 current calculation restart file
-P-0000 ------------------- ------------
-P-0000
-P-0000 calculation expects a wf_planewave | input file contains a wf_planewave
-P-0000 . ABINIT code version 6.12.3 | ABINIT code version 6.12.3
-P-0000 . date 20120719 bantot 416 natom 9 | date 20120719 bantot 260 natom 9
-P-0000 nkpt 16 nsym 2 ngfft 45, 45, 48 | nkpt 10 nsym 2 ngfft 45, 45, 48
-P-0000 ntypat 2 ecut_eff 24.0000000 | ntypat 2 ecut_eff 24.0000000
-P-0000 usepaw 0 | usepaw 0
-P-0000 usewvl 0 | usewvl 0
-P-0000 rprimd: | rprimd:
-P-0000 9.3824902 0.0000000 0.0000000 | 9.3824902 0.0000000 0.0000000
-P-0000 -4.6912451 8.1254749 0.0000000 | -4.6912451 8.1254749 0.0000000
-P-0000 0.0000000 0.0000000 10.2498745 | 0.0000000 0.0000000 10.2498745
-P-0000
-P-0000 hdr_check: WARNING -
-P-0000 input nkpt= 16 not equal disk file nkpt= 10
-P-0000 symafm: | symafm:
-P-0000 1 1 | 1 1
-P-0000 symrel: | symrel:
-P-0000 1 0 0 0 1 0 0 0 1 | 1 0 0 0 1 0 0 0 1
-P-0000 0 1 0 1 0 0 0 0 -1 | 0 1 0 1 0 0 0 0 -1
-P-0000 typat: | typat:
-P-0000 1 1 1 2 2 2 2 2 2 | 1 1 1 2 2 2 2 2 2
-P-0000 so_psp : | so_psp :
-P-0000 1 1 | 1 1
-P-0000 tnons: | tnons:
-P-0000 0.0000000 0.0000000 0.0000000 | 0.0000000 0.0000000 0.0000000
-P-0000 -0.0000000 0.0000000 -0.0000000 | -0.0000000 0.0000000 -0.0000000
-P-0000 znucl: | znucl:
-P-0000 14.00 8.00 | 14.00 8.00
-P-0000 pseudopotential atom type 1: | pseudopotential atom type 1:
-P-0000 pspso 0 pspxc 1 | pspso 0 pspxc 1
TACC: Cleaning up after job: 2698534
TACC: Done.

Could that be also a reason why this calculation fails?
Should I include such variables as npkpt npband in my input file?? Please help!

well, I think your shiftk is not compatible with a primitive hexagonal cell as you have here. Do an initial run with kptrlen 60 (say), and prtkpt 1, so abinit can tell you what are good values for kptrlatt and shiftk for your system.Also, by the way, your kpt grid looks VERY coarse right now. Before doing a phonon calc I would make sure that the ground state stress is accurately converged with both ecut and kpt mesh.