### Comparing pwscf and Abinit MPI SCF calculations

Posted:

**Tue Mar 30, 2021 5:59 pm**Hi,

I want to compare the energies and forces found by using a non-local pseudopotential and a local pseudopotential. The non-local pseudopotential calculation was performed in pwscf (Quantum Espresso) using a UPF norm-conserving pseudopotential; however, due to the limitations on pseudopotential format in pwscf, I need to use Abinit for the local pseudopotential calculation.

Bear in mind I am new to Abinit, so there may be something obvious I have missed, but the Abinit convergence seems far far slower. I have attached both input and output files (note that the Abinit calculation ran out of time so is incomplete and is in fact the "log" file), to hopefully shed some light on what is going on.

The system is a 108 atom Al cell which has been equilibrated at 100 kK using MD. Therefore, I would also like to perform the "ground-state" (perhaps "thermal equilibrium state" is a more accurate term) minimisation at an electron temperature of 100 kK (= 0.6334 Ry = 0.3167 Ha) with Fermi-Dirac smearing. This is the reason for the large (4500) number of bands.

pwscf seems reasonably simple to run in parallel; it automatically detected the number of available cores and parallelised accordingly. Abinit is a bit more complicated, but I have used the autoparal to hopefully find the best allocation of resources. This is where I believe the biggest bottleneck is probably found.

Each calculation was given 12 hrs with 16GB per core and 40 cores. pwscf predicted it would need 367GB of RAM, whereas Abinit predicts it only needs 1GB of memory, which implies either I don't understand how to read the Abinit log file, or it is not making use of the available memory fully. Abinit also used all 12 hrs and didn't complete, compared to pwscf completing the calculation in ~5 hrs.

I've tried to keep all other parameters (e.g. ecut, etol, nkpt, etc.) as similar between the two runs as possible to minimise any other differences, and I am not convinced that the local pseudopotential could cause such a slow down in comparison to the non-local pseudopotential. I am, however, going to do another identical run but with a non-local pseudopotential in Abinit to check this.

Any ideas on what I am doing incorrectly or could change would be much appreciated.

Thanks.

I want to compare the energies and forces found by using a non-local pseudopotential and a local pseudopotential. The non-local pseudopotential calculation was performed in pwscf (Quantum Espresso) using a UPF norm-conserving pseudopotential; however, due to the limitations on pseudopotential format in pwscf, I need to use Abinit for the local pseudopotential calculation.

Bear in mind I am new to Abinit, so there may be something obvious I have missed, but the Abinit convergence seems far far slower. I have attached both input and output files (note that the Abinit calculation ran out of time so is incomplete and is in fact the "log" file), to hopefully shed some light on what is going on.

The system is a 108 atom Al cell which has been equilibrated at 100 kK using MD. Therefore, I would also like to perform the "ground-state" (perhaps "thermal equilibrium state" is a more accurate term) minimisation at an electron temperature of 100 kK (= 0.6334 Ry = 0.3167 Ha) with Fermi-Dirac smearing. This is the reason for the large (4500) number of bands.

pwscf seems reasonably simple to run in parallel; it automatically detected the number of available cores and parallelised accordingly. Abinit is a bit more complicated, but I have used the autoparal to hopefully find the best allocation of resources. This is where I believe the biggest bottleneck is probably found.

Each calculation was given 12 hrs with 16GB per core and 40 cores. pwscf predicted it would need 367GB of RAM, whereas Abinit predicts it only needs 1GB of memory, which implies either I don't understand how to read the Abinit log file, or it is not making use of the available memory fully. Abinit also used all 12 hrs and didn't complete, compared to pwscf completing the calculation in ~5 hrs.

I've tried to keep all other parameters (e.g. ecut, etol, nkpt, etc.) as similar between the two runs as possible to minimise any other differences, and I am not convinced that the local pseudopotential could cause such a slow down in comparison to the non-local pseudopotential. I am, however, going to do another identical run but with a non-local pseudopotential in Abinit to check this.

Any ideas on what I am doing incorrectly or could change would be much appreciated.

Thanks.