Files prepared and density convergence / ELF calculation put into the queue on 32 nodes with 8 ppn and a walltime of 36 hours.

Density successfully converged after 16 SCF steps and ELF printed.
Final etot = -505.54887994153 Hartrees.
Final diff in etot = -2.309E-09.
nkpt = 258.
natom = 80.
nband = 144.
ngfft = 120 45 100.
ecut = 25 Hartrees.
Run on 32 nodes, 8 ppn.
It took 7804.1 seconds per processor = 2 hours 10 minutes
Size of abinit ELF output file: 4.6 mb
Size of ELF vesta file: 12.5 mb

Next step is to change the input file to read in the density, converge the wavefunction, and print the PDOS. I'm going to let iatsph and ratsph be set to default. Default iatsph goes from 1 to natsph, exactly as we want. Default ratsph sets all r to 2.0 Bohr for norm conserving pseudopotentials, which I believe we have. If the calculation is successful, we can tweak the ratsph and see what happens. The calculation is now in the queue.

When run on 32 nodes, 8 ppn = 256 processors, calculation runs out of memory. Will rerun it with 64 nodes, 4 ppn = 256 processors while adding #PBS -l pvmem=5GB to the submission file.

The memory problems have been avoided and the PDOS has successfully calculated, however before the calculation can finish, the following error occurs:

At line 121 of file hdr_skip.F90 (unit = 12, file = 'mz84o_DOS_AT0004')
Fortran runtime error: Missing format for FORMATTED data transfer

Furthermore, the PDOS files for atoms 3 (the first of a group of oxygens) and 9 (the first of a group of carbons) are nowhere to be found, though in the log file it says they were written. The other nine printed PDOS files do appear to be good though, therefore I'm trying to rerun the same calculation for mz133 while only printing the PDOS file for atom 3 to see if we can "fill in the gaps", so to speak.

I've successfully calculated the PDOS for atom 9 of mz133, therefore I'm going to try and calculate PDOS for atoms 3 and 9 for this compound to see if I can do them together. Calculation is in the queue now.

PDOS for atoms 3 and 9 calculated successfully. Here are the specs:
nkpt = 258.
natom = 80.
nband = 180.
ngfft = 120 45 100.
ecut = 25 Hartrees.
Run on 64 nodes, 4 ppn.
It took 10493.2 seconds per processor = 2 hours 55 minutes

Therefore the PDOS is complete. Hurrah!

Alex took a look at the abinit generated ELFs and the whole induced inert pair stuff is not at all present. Because we trust the ELFs from abinit a lot more than those from LMTO this means that this part of the story is simply gone. Kinda sucks, but it's better to be right.

I haven't really dug into the PDOS yet, however I did use it to calculate the band gap = 3.4022 eV compared to mz133 = 3.3599 eV.

Te coordination examined with Shape. See Shape Calculations for more information.

Unless otherwise stated, the content of this page is licensed under Creative Commons Attribution-ShareAlike 3.0 License