Comparison of 15.09 FLO-SIC code against Perfect code. Energies of first and second iterations are compared.
Discrepancy at the first iteration tells you SIC energy calculation is wrong, and discrepancy at the second iteration tells you SIC matrix is calculated wrong.
DFT Energy Comparison | |||||||
DFT | 1509 code | perfect | |||||
System | IT | ETOT | EDFT | Difference | |||
O2 | 1 | -149.270313466 | -149.270313466 | 0.00E+00 | |||
2 | -149.320322245 | -149.320322245 | 0.00E+00 | ||||
C | 1 | -37.466448632 | -37.466448632 | 0.00E+00 | |||
2 | -37.467806082 | -37.467806082 | 0.00E+00 | ||||
Ne | 1 | -128.229726561 | -128.229726561 | 0.00E+00 | |||
2 | -128.229726888 | -128.229726888 | 0.00E+00 | ||||
Si | 1 | -288.217076863 | -288.217076863 | 0.00E+00 | |||
2 | -288.217404963 | -288.217404963 | 0.00E+00 | ||||
Ar | 1 | -525.939481957 | -525.939481957 | 0.00E+00 | |||
2 | -525.93948397 | -525.93948397 | 0.00E+00 | ||||
Cu | 1 | -1637.778538616 | -1637.778538616 | 0.00E+00 | |||
2 | -1637.758717915 | -1637.758717917 | 2.00E-09 | ||||
Zn | 1 | -1776.561063172 | -1776.561063173 | 1.00E-09 | |||
2 | -1776.56106501 | -1776.56106501 | 0.00E+00 | ||||
As | 1 | -2232.571800009 | -2232.571800009 | 0.00E+00 | |||
2 | -2232.572666378 | -2232.572666378 | 0.00E+00 | ||||
Kr | 1 | -2750.132934998 | -2750.132934998 | 0.00E+00 | |||
2 | -2750.132962212 | -2750.132962212 | 0.00E+00 | ||||
N2 | 1 | -108.666809409 | -108.666809409 | 0.00E+00 | |||
2 | -108.691396857 | -108.691396857 | 0.00E+00 | ||||
SIC Energy Comparison | |||||||
SIC | 1509 code | perfect+fermocc+newwaveour | |||||
IT | ETOT | EDFT+SIC | EDFT | EDFT+SIC | difference | ||
O2 | 1 | -149.270313466 | -150.43551959 | -149.270313466 | -150.435519589 | 0.00E+00 | -1.00E-09 |
2 | -149.313961185 | -150.528246116 | -149.313961186 | -150.528246114 | 1.00E-09 | -2.00E-09 | |
C | 1 | -37.465962808 | -37.938618958 | -37.465962808 | -37.939900385 | 0.00E+00 | 1.28E-03 |
2 | -37.462300658 | -37.947549537 | -37.46229779 | -37.949067952 | -2.87E-06 | 1.52E-03 | |
Ne | 1 | -128.229726561 | -129.267624633 | -128.229726561 | -129.267624631 | 0.00E+00 | -2.00E-09 |
2 | -128.212353051 | -129.278524947 | -128.212353052 | -129.278524945 | 1.00E-09 | -2.00E-09 | |
Si | 1 | -288.215517308 | -289.996877871 | -288.215517308 | -289.997759511 | 0.00E+00 | 8.82E-04 |
2 | -288.206194342 | -290.008093002 | -288.206175697 | -290.008691422 | -1.86E-05 | 5.98E-04 | |
Ar | 1 | -525.939481957 | -528.522213441 | -525.939481957 | -528.522213436 | 0.00E+00 | -5.00E-09 |
2 | -525.925261427 | -528.535665809 | -525.925261427 | -528.535665803 | 0.00E+00 | -6.00E-09 | |
Cu | 1 | -1637.778538616 | -1642.93362471 | -1637.778538616 | -1642.933624715 | 0.00E+00 | 5.00E-09 |
2 | -1637.678478171 | -1642.941458513 | -1637.678478168 | -1642.941458518 | -3.00E-09 | 5.00E-09 | |
Zn | 1 | -1776.561063172 | -1782.054634963 | -1776.561063173 | -1782.054634979 | 1.00E-09 | 1.60E-08 |
2 | -1776.466620381 | -1782.064793862 | -1776.466620378 | -1782.064793876 | -3.00E-09 | 1.40E-08 | |
As | 1 | -2232.571800009 | -2239.071003081 | -2232.571800009 | -2239.071003081 | 0.00E+00 | 0.00E+00 |
2 | -2232.527365829 | -2239.090088444 | -2232.52736583 | -2239.090088443 | 1.00E-09 | -1.00E-09 | |
Kr | 1 | -2750.132934998 | -2757.579879852 | -2750.132934998 | -2757.579879855 | 0.00E+00 | 3.00E-09 |
2 | -2750.096332313 | -2757.602179869 | -2750.096332313 | -2757.602179872 | 0.00E+00 | 3.00E-09 | |
N2 | 1 | -108.666809409 | -109.819197311 | -108.666809409 | -109.819197311 | 0.00E+00 | 0.00E+00 |
2 | -108.67679902 | -109.853899069 | -108.67679902 | -109.853899069 | 0.00E+00 | 0.00E+00 | |
MPI vs Serial | HPC Ifort | ||||||
1509 code | MPI np 12 | Serial | |||||
IT | ETOT | EDFT+SIC | ETOT | EDFT+SIC | difference | ||
O2 | 1 | -149.270313466 | -150.43551959 | -149.270313466 | -150.43551959 | 0.00E+00 | 0.00E+00 |
2 | -149.313961186 | -150.528246116 | -149.313961185 | -150.528246116 | -1.00E-09 | 0.00E+00 | |
C | 1 | -37.465962808 | -37.939879844 | -37.465962808 | -37.939879844 | 0.00E+00 | 0.00E+00 |
2 | -37.46229871 | -37.949044701 | -37.46229871 | -37.949044701 | 0.00E+00 | 0.00E+00 | |
Ne | 1 | -128.229726561 | -129.267624633 | -128.229726561 | -129.267624633 | 0.00E+00 | 0.00E+00 |
2 | -128.212353052 | -129.278524947 | -128.212353051 | -129.278524947 | -1.00E-09 | 0.00E+00 | |
Si | 1 | -288.215517308 | -289.996464126 | -288.215517308 | -289.996464126 | 0.00E+00 | 0.00E+00 |
2 | -288.206181268 | -290.00777049 | -288.206181268 | -290.007770489 | 0.00E+00 | -1.00E-09 | |
Ar | 1 | -525.939481957 | -528.522213441 | -525.939481957 | -528.522213441 | 0.00E+00 | 0.00E+00 |
2 | -525.925261428 | -528.535665809 | -525.925261427 | -528.535665809 | -1.00E-09 | 0.00E+00 | |
Cu | 1 | -1637.778538616 | -1642.933624709 | -1637.778538616 | -1642.933624709 | 0.00E+00 | 0.00E+00 |
2 | -1637.678478203 | -1642.941458525 | -1637.67847817 | -1642.941458513 | -3.30E-08 | -1.20E-08 | |
Zn | 1 | -1776.561063172 | -1782.054634963 | -1776.561063172 | -1782.054634963 | 0.00E+00 | 0.00E+00 |
2 | -1776.466620406 | -1782.064793871 | -1776.466620381 | -1782.064793862 | -2.50E-08 | -9.00E-09 | |
As | 1 | -2232.571800009 | -2239.071003081 | -2232.571800009 | -2239.071003081 | 0.00E+00 | 0.00E+00 |
2 | -2232.527365849 | -2239.090088449 | -2232.527365829 | -2239.090088444 | -2.00E-08 | -5.00E-09 | |
Kr | 1 | -2750.132934998 | -2757.579879857 | -2750.132934998 | -2757.579879857 | 0.00E+00 | 0.00E+00 |
2 | -2750.096332345 | -2757.602179876 | -2750.096332313 | -2757.60217987 | -3.20E-08 | -6.00E-09 | |
N2 | 1 | -108.666809409 | -109.819197311 | -108.666809409 | -109.819197311 | 0.00E+00 | 0.00E+00 |
2 | -108.67679902 | -109.853899069 | -108.67679902 | -109.853899069 | 0.00E+00 | 0.00E+00 | |
MPI vs Serial | HPC gfortran | ||||||
1509 code | MPI np 8 | ||||||
IT | ETOT | EDFT+SIC | ETOT | EDFT+SIC | difference | ||
O2 | 1 | -149.270313466 | -150.43551959 | -149.270313466 | -150.435519589 | 0.00E+00 | -1.00E-09 |
2 | -149.313961186 | -150.528246116 | -149.313961186 | -150.528246114 | 0.00E+00 | -2.00E-09 | |
C | 1 | -37.465962808 | -37.939791134 | -37.465962808 | -37.938618958 | 0.00E+00 | -1.17E-03 |
2 | -37.462297359 | -37.948949253 | -37.462300658 | -37.947549537 | 3.30E-06 | -1.40E-03 | |
Ne | 1 | -128.229726561 | -129.267624633 | -128.229726561 | -129.267624633 | 0.00E+00 | 0.00E+00 |
2 | -128.212353052 | -129.278524947 | -128.212353051 | -129.278524947 | -1.00E-09 | 0.00E+00 | |
Si | 1 | -288.215517308 | -289.996605824 | -288.215517308 | -289.996877871 | 0.00E+00 | 2.72E-04 |
2 | -288.206185843 | -290.007884028 | -288.206194342 | -290.008093002 | 8.50E-06 | 2.09E-04 | |
Ar | 1 | -525.939481957 | -528.522213441 | -525.939481957 | -528.522213441 | 0.00E+00 | 0.00E+00 |
2 | -525.925261428 | -528.535665809 | -525.925261427 | -528.535665809 | -1.00E-09 | 0.00E+00 | |
Cu | 1 | -1637.778538613 | -1642.933624709 | -1637.778538616 | -1642.93362471 | 3.00E-09 | 1.00E-09 |
2 | -1637.678478201 | -1642.941458527 | -1637.678478171 | -1642.941458513 | -3.00E-08 | -1.40E-08 | |
Zn | 1 | -1776.561063172 | -1782.054634963 | -1776.561063172 | -1782.054634963 | 0.00E+00 | 0.00E+00 |
2 | -1776.466620406 | -1782.064793871 | -1776.466620381 | -1782.064793862 | -2.50E-08 | -9.00E-09 | |
As | 1 | -2232.571800003 | -2239.071003076 | -2232.571800009 | -2239.071003081 | 6.00E-09 | 5.00E-09 |
2 | -2232.527365843 | -2239.090088443 | -2232.527365829 | -2239.090088444 | -1.40E-08 | 1.00E-09 | |
Kr | 1 | -2750.132934998 | -2757.579879852 | -2750.132934998 | -2757.579879852 | 0.00E+00 | 0.00E+00 |
2 | -2750.096332345 | -2757.602179876 | -2750.096332313 | -2757.602179869 | -3.20E-08 | -7.00E-09 | |
N2 | 1 | -108.666809409 | -109.819197311 | -108.666809409 | -109.819197311 | 0.00E+00 | 0.00E+00 |
2 | -108.67679902 | -109.853899069 | -108.67679902 | -109.853899069 | 0.00E+00 | 0.00E+00 |