PseudoDojo notebook for Ar

In [1]:
from __future__ import print_function, division, unicode_literals
%matplotlib inline
In [2]:
from IPython.display import HTML
HTML('''<script>
code_show=true;
function code_toggle() {
 if (code_show){
 $('div.input').hide();
 } else {
 $('div.input').show();
 }
 code_show = !code_show
}
$( document ).ready(code_toggle);
</script>
The raw code for this IPython notebook is by default hidden for easier reading.
To toggle on/off the raw code, click <a href="javascript:code_toggle()">here</a>.''')
Out[2]:
The raw code for this IPython notebook is by default hidden for easier reading. To toggle on/off the raw code, click here.
In [3]:
from pseudo_dojo.core.pseudos import dojopseudo_from_file
pseudo = dojopseudo_from_file('/Users/setten/software/conda_install/my_stack/pseudo_dojo/pseudo_dojo/pseudos/website/pw_h_sr/Ar.psp8')
report = pseudo.dojo_report

ONCVPSP Input File

In [4]:
input_file = pseudo.filepath.replace(".psp8", ".in")
%cat $input_file

# atsym z nc nv iexc psfile
Ar 18 3 2 -001012 psp8
# n l f
1 0 2.0
2 0 2.0
2 1 6.0
3 0 2.0
3 1 6.0
# lmax
2
# l rc ep ncon nbas qcut
0 1.55 0.0 4 8 8.2
1 1.55 0.0 4 8 8.2
2 1.55 0.1 4 8 8.2
# lloc lpopt rc5 dvloc0
4 5 1.5 0.0
# l nproj debl
0 2 1.5
1 2 1.5
2 2 1.0
# icmod fcfact rcfact
3 7.0 1.25
# epsh1 epsh2 depsh
-12.0 12.0 0.02
# rlmax drl
4.0 0.01
 0
In [5]:
# Get data from the oncvpsp output file
from pseudo_dojo.ppcodes.oncvpsp import OncvOutputParser, PseudoGenDataPlotter
onc_parser = OncvOutputParser(pseudo.filepath.replace(".psp8", ".out"))

# Parse the file and build the plotter
onc_parser.scan()
plotter = onc_parser.make_plotter()

AE and PS radial wavefunctions $\phi(r)$

In [6]:
fig = plotter.plot_radial_wfs(show=False)

Arctan of the logarithmic derivatives

From oncvpsp documentation: The plots show $\phi(E) = \arctan(R * d \psi(r)/dr |_R)$ for some $R$ greater than the core radius, where $\psi$ is the solution of the non-local radial equation regular at the origin (i.e., the outward-integrated solution). For a well-designed pseudopotential, $\phi(E)$ will closely track that of the all-electron potential over a wide range of energies from well-below to well-above the valence semi-core states of interest. The steps of $\pi$ indicate localized pseudo wave functions. Spurious steps of $\pi$ indicate "ghost" states, which are localized states than on investigation turn out to have more nodes than appropriate for their energies.

For $GW$ pseudos, no significant deviation should be present up to 8 Hartree.

In [7]:
fig = plotter.plot_atan_logders(show=False)

Convergence in $G$-space estimated by ONCVPSP

These results are obtained in the atomic configuration and should give a reasonable estimate of the convergence behaviour wrt to ecut in crystalline systems.

In [8]:
fig = plotter.plot_ene_vs_ecut(show=False)

Projectors

In general the second projector in any channel should have one node more that the first one. Pushing the energy of the second projector too high may cause an additional node. This will most likely introduce ghosts.

In [9]:
fig = plotter.plot_projectors(show=False)

Core-Valence-Model charge densities

Much better convergence properties can been achieved with icmod 3. In this case, fcfact mainly determines the height of the model core charge while rcfact mainly determines the width of the model core charge.

In [10]:
fig = plotter.plot_densities(show=False)

Local potential and $l$-dependent potentials

In [11]:
fig = plotter.plot_potentials(show=False)

Ghosts Test

Self-consistent band structure calculation performed on a regular mesh. The algorithm to detect ghosts is just an indication usually on the side of false positives. Zoom in on the band plot to see if an actual ghost is there.

In [12]:
fig = report.plot_ebands(with_soc=False, show=False); fig

Found 133 dispersionless states
[0] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1, eig=-19.063, occ=2.000
[1] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=2, eig=-19.063, occ=2.000
[2] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=3, eig=-19.063, occ=2.000
[3] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=311, eig=99.976, occ=0.000
[4] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=312, eig=99.976, occ=0.000
[5] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=355, eig=110.506, occ=0.000
[6] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=356, eig=110.506, occ=0.000
[7] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=357, eig=110.506, occ=0.000
[8] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=418, eig=121.380, occ=0.000
[9] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=419, eig=121.380, occ=0.000
[10] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=420, eig=121.380, occ=0.000
[11] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=552, eig=152.865, occ=0.000
[12] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=553, eig=152.865, occ=0.000
[13] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=554, eig=152.865, occ=0.000
[14] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=555, eig=152.869, occ=0.000
[15] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=556, eig=152.869, occ=0.000
[16] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=557, eig=152.869, occ=0.000
[17] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=571, eig=155.560, occ=0.000
[18] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=572, eig=155.560, occ=0.000
[19] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=573, eig=155.560, occ=0.000
[20] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=574, eig=155.569, occ=0.000
[21] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=581, eig=157.389, occ=0.000
[22] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=582, eig=157.389, occ=0.000
[23] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=606, eig=161.076, occ=0.000
[24] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=607, eig=161.076, occ=0.000
[25] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=608, eig=161.076, occ=0.000
[26] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=609, eig=161.084, occ=0.000
[27] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=610, eig=161.084, occ=0.000
[28] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=636, eig=167.017, occ=0.000
[29] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=637, eig=167.017, occ=0.000
[30] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=638, eig=167.017, occ=0.000
[31] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=646, eig=168.545, occ=0.000
[32] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=647, eig=168.545, occ=0.000
[33] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=648, eig=168.545, occ=0.000
[34] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=687, eig=174.543, occ=0.000
[35] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=688, eig=174.543, occ=0.000
[36] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=689, eig=174.543, occ=0.000
[37] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=699, eig=175.812, occ=0.000
[38] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=700, eig=175.812, occ=0.000
[39] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=701, eig=175.812, occ=0.000
[40] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=705, eig=176.889, occ=0.000
[41] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=706, eig=176.889, occ=0.000
[42] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=707, eig=176.889, occ=0.000
[43] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=714, eig=178.576, occ=0.000
[44] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=715, eig=178.576, occ=0.000
[45] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=716, eig=178.576, occ=0.000
[46] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=730, eig=182.518, occ=0.000
[47] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=731, eig=182.518, occ=0.000
[48] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=732, eig=182.518, occ=0.000
[49] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=744, eig=184.010, occ=0.000
[50] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=745, eig=184.010, occ=0.000
[51] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=746, eig=184.010, occ=0.000
[52] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=747, eig=185.224, occ=0.000
[53] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=759, eig=187.086, occ=0.000
[54] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=760, eig=187.086, occ=0.000
[55] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=761, eig=187.086, occ=0.000
[56] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=765, eig=188.263, occ=0.000
[57] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=766, eig=188.263, occ=0.000
[58] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=767, eig=188.263, occ=0.000
[59] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=774, eig=192.592, occ=0.000
[60] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=775, eig=192.592, occ=0.000
[61] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=795, eig=197.727, occ=0.000
[62] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=796, eig=197.727, occ=0.000
[63] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=797, eig=197.727, occ=0.000
[64] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=798, eig=198.007, occ=0.000
[65] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=799, eig=198.007, occ=0.000
[66] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=800, eig=198.007, occ=0.000
[67] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=849, eig=205.750, occ=0.000
[68] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=850, eig=205.750, occ=0.000
[69] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=851, eig=205.750, occ=0.000
[70] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=853, eig=205.988, occ=0.000
[71] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=854, eig=205.988, occ=0.000
[72] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=855, eig=205.988, occ=0.000
[73] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=876, eig=208.819, occ=0.000
[74] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=877, eig=208.819, occ=0.000
[75] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=878, eig=208.819, occ=0.000
[76] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=879, eig=208.857, occ=0.000
[77] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=882, eig=209.685, occ=0.000
[78] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=883, eig=209.685, occ=0.000
[79] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=884, eig=209.685, occ=0.000
[80] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=889, eig=210.448, occ=0.000
[81] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=890, eig=210.448, occ=0.000
[82] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=897, eig=211.844, occ=0.000
[83] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=898, eig=211.868, occ=0.000
[84] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=899, eig=211.868, occ=0.000
[85] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=915, eig=215.185, occ=0.000
[86] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=916, eig=215.185, occ=0.000
[87] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=917, eig=215.185, occ=0.000
[88] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=921, eig=216.373, occ=0.000
[89] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=922, eig=216.373, occ=0.000
[90] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=923, eig=216.373, occ=0.000
[91] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=939, eig=219.019, occ=0.000
[92] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=940, eig=219.019, occ=0.000
[93] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=941, eig=219.019, occ=0.000
[94] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=963, eig=222.572, occ=0.000
[95] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=964, eig=222.572, occ=0.000
[96] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=965, eig=222.572, occ=0.000
[97] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1002, eig=228.838, occ=0.000
[98] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1003, eig=228.838, occ=0.000
[99] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1004, eig=228.838, occ=0.000
[100] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1017, eig=231.589, occ=0.000
[101] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1018, eig=231.589, occ=0.000
[102] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1053, eig=237.633, occ=0.000
[103] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1059, eig=238.547, occ=0.000
[104] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1060, eig=238.547, occ=0.000
[105] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1061, eig=238.547, occ=0.000
[106] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1062, eig=238.585, occ=0.000
[107] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1063, eig=238.718, occ=0.000
[108] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1064, eig=238.718, occ=0.000
[109] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1065, eig=238.718, occ=0.000
[110] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1074, eig=241.707, occ=0.000
[111] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1075, eig=241.707, occ=0.000
[112] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1076, eig=241.707, occ=0.000
[113] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1077, eig=242.404, occ=0.000
[114] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1078, eig=242.404, occ=0.000
[115] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1079, eig=242.404, occ=0.000
[116] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1080, eig=242.422, occ=0.000
[117] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1081, eig=242.422, occ=0.000
[118] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1082, eig=242.422, occ=0.000
[119] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1089, eig=243.699, occ=0.000
[120] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1090, eig=243.699, occ=0.000
[121] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1091, eig=243.699, occ=0.000
[122] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1117, eig=246.857, occ=0.000
[123] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1146, eig=250.905, occ=0.000
[124] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1147, eig=250.907, occ=0.000
[125] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1148, eig=250.909, occ=0.000
[126] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1149, eig=250.963, occ=0.000
[127] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1150, eig=250.965, occ=0.000
[128] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1151, eig=250.967, occ=0.000
[129] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1152, eig=251.168, occ=0.000
[130] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1157, eig=252.003, occ=0.000
[131] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1158, eig=252.035, occ=0.000
[132] spin=0, kpt=[0.063, 0.063, 0.063], weight=0.016, band=1159, eig=252.078, occ=0.000
Out[12]:

Convergence of the total energy wrt ecut

The energies are obtained from the deltafactor calculations performed at the Wien2K equilibrium volume

In [13]:
fig = report.plot_etotal_vs_ecut(show=False)
In [14]:
fig = report.plot_etotal_vs_ecut(inv_ecut=True, show=False)

Convergence of the Delta-Gauge results

The Delta-gauge compares the Equation Of State (EOS) of the elemental solid of the element calculated using the pseudo potential to reference curves calculated using an all electron method. The Delta-gauge was introduced by K. Lejaeghere, V. Van Speybroeck, G. Van Oost, and&S. Cottenier in Critical Reviews in Solid State and Materials Sciences 39, 1

A comparison using the Delta-gauge between many codes and many pseudo tables can be found at the center for molecular modeling and in Science 351, 1394-1395

In [15]:
fig = report.plot_deltafactor_convergence(xc=pseudo.xc, what=("dfact_meV", "dfactprime_meV"), show=False)

Reference data: DeltaFactorEntry(symbol='Ar', v0=30.12052, b0=0.04405881291960992, b1=5.731, xc=PW)

Convergence of $\Delta v_0$, $\Delta b_0$, and $\Delta b_1$ (deltafactor tests)

In [16]:
# Absolute difference with respect to Wien2k results.
fig = report.plot_deltafactor_convergence(xc=pseudo.xc, what=("-dfact_meV", "-dfactprime_meV"), show=False)

Reference data: DeltaFactorEntry(symbol='Ar', v0=30.12052, b0=0.04405881291960992, b1=5.731, xc=PW)

Delta-gauge EOS for the different cutoff energies

In [17]:
fig = report.plot_deltafactor_eos(show=False)

Convergence of the GBRV lattice parameters

The GBRV tests compare the lattice parameter of a FCC and BCC lattice of the element to all electron reference data. The test was introduced by Kevin F. Garrity, Joseph W. Bennett, Karin M. Rabe, and David Vanderbilt in developing th GBRV pseudo potential table. More information can be found in Computational Materials Science 81, 446-452.

In [18]:
fig = report.plot_gbrv_convergence(show=False)

dojo report does not contain trial: gbrv_fcc

Convergence of the phonon frequencies at $\Gamma$

The calculation is performed with the Wien2k relaxed parameters obtained from the deltafactor CIF files.

In [19]:
fig = report.plot_phonon_convergence(show=False)

dojo report does not contain trial: phgamma

GBRV EOS for the FCC structure

In [20]:
fig = report.plot_gbrv_eos(struct_type="fcc", show=False)

dojo report does not contain trial: gbrv_fcc

GBRV EOS for the BCC structure

In [21]:
fig = report.plot_gbrv_eos(struct_type="bcc", show=False)

dojo report does not contain trial: gbrv_bcc