-
Notifications
You must be signed in to change notification settings - Fork 0
Expand file tree
/
Copy pathupf-spec.html
More file actions
685 lines (646 loc) · 22.2 KB
/
upf-spec.html
File metadata and controls
685 lines (646 loc) · 22.2 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
182
183
184
185
186
187
188
189
190
191
192
193
194
195
196
197
198
199
200
201
202
203
204
205
206
207
208
209
210
211
212
213
214
215
216
217
218
219
220
221
222
223
224
225
226
227
228
229
230
231
232
233
234
235
236
237
238
239
240
241
242
243
244
245
246
247
248
249
250
251
252
253
254
255
256
257
258
259
260
261
262
263
264
265
266
267
268
269
270
271
272
273
274
275
276
277
278
279
280
281
282
283
284
285
286
287
288
289
290
291
292
293
294
295
296
297
298
299
300
301
302
303
304
305
306
307
308
309
310
311
312
313
314
315
316
317
318
319
320
321
322
323
324
325
326
327
328
329
330
331
332
333
334
335
336
337
338
339
340
341
342
343
344
345
346
347
348
349
350
351
352
353
354
355
356
357
358
359
360
361
362
363
364
365
366
367
368
369
370
371
372
373
374
375
376
377
378
379
380
381
382
383
384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
473
474
475
476
477
478
479
480
481
482
483
484
485
486
487
488
489
490
491
492
493
494
495
496
497
498
499
500
501
502
503
504
505
506
507
508
509
510
511
512
513
514
515
516
517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
538
539
540
541
542
543
544
545
546
547
548
549
550
551
552
553
554
555
556
557
558
559
560
561
562
563
564
565
566
567
568
569
570
571
572
573
574
575
576
577
578
579
580
581
582
583
584
585
586
587
588
589
590
591
592
593
594
595
596
597
598
599
600
601
602
603
604
605
606
607
608
609
610
611
612
613
614
615
616
617
618
619
620
621
622
623
624
625
626
627
628
629
630
631
632
633
634
635
636
637
638
639
640
641
642
643
644
645
646
647
648
649
650
651
652
653
654
655
656
657
658
659
660
661
662
663
664
665
666
667
668
669
670
671
672
673
674
675
676
677
678
679
680
681
682
683
684
685
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html lang="en">
<head>
<meta charset="utf-8" />
<title>Unified Pseudopotential Format - QUANTUMESPRESSO</title>
<meta name="generator" content="GravCMS" />
<meta name="description" content="Grav is an easy to use, yet powerful, open source flat-file CMS" />
<meta http-equiv="X-UA-Compatible" content="IE=edge" />
<meta name="robots" content="follow, all" />
<meta name="google-site-verification" content="g4CG61y_kLfkqAsfUPkF2MefchgN00GozvDUZMeK4z8" /> <link rel="canonical" href="https://pseudopotentials.quantum-espresso.org/home/unified-pseudopotential-format" />
<link href="/user/themes/quantum/css/style.css" type="text/css" rel="stylesheet">
<link href="/user/themes/quantum/css/style-red.css" type="text/css" rel="stylesheet">
<link href="/user/themes/quantum/css/superfish.css" type="text/css" rel="stylesheet">
<link href="/user/themes/quantum/css/bootstrap.css" type="text/css" rel="stylesheet">
<link href="/user/themes/quantum/css/prettyPhoto.css" type="text/css" rel="stylesheet">
<link href="https://fonts.googleapis.com/css?family=Lato:300" type="text/css" rel="stylesheet">
<link href="/user/plugins/markdown-notices/assets/notices.css" type="text/css" rel="stylesheet">
<link href="/user/plugins/form/assets/form-styles.css" type="text/css" rel="stylesheet">
<link href="/user/plugins/mathjax/assets/css/mathjax.css" type="text/css" rel="stylesheet">
<link href="/user/plugins/login/css/login.css" type="text/css" rel="stylesheet">
<script src="/user/themes/quantum/js/jquery.js"></script>
<script src="/user/plugins/mathjax/assets/js/mathjax.js"></script>
<script src="https://cdnjs.cloudflare.com/ajax/libs/mathjax/2.7.0/MathJax.js?config=TeX-AMS-MML_HTMLorMML"></script>
<script src="/user/themes/quantum/js/jquery.cycle.all.js"></script>
<script src="/user/themes/quantum/js/custom.js"></script>
<script src="/user/themes/quantum/js/superfish.js"></script>
<script src="/user/themes/quantum/js/jquery.prettyPhoto.js"></script>
<script src="/user/themes/quantum/js/functions.js"></script>
<script>
</script>
</head>
<body class="">
<div id="outer-container">
<div id="inner-container">
<!-- BEGIN HEADER -->
<div id="header">
<div id="logo">
<a href="/">
<img src="/user/themes/quantum/images/logo_header.jpg" alt="QUANTUMESPRESSO" />
</a>
</div>
<div id="nav-menu">
<ul id="menu-menu" class="sf-menu">
<li class="menu-item" >
<a class="sf-with-ul" href="http://www.quantum-espresso.org">
HOME
</a>
</li>
<li class="current_page_item menu-item ">
<a href="/" >
Pseudopotentials
</a>
</li>
</ul>
</div>
</div>
<!-- END OF HEADER -->
<!-- BEGIN CONTENT -->
<div id="content-inner">
<div id="left-bar">
<div class="left-bar-content">
<h2>Resources</h2>
<ul class="leftbar-list">
<li class="page_item page-item-36">
<a href="https://www.materialscloud.org/discover/sssp/table/efficiency" title="SSSP">SSSP</a>
</li>
<li class="page_item page-item-36">
<a href="https://dalcorso.github.io/pslibrary/" title="PSlibrary">PSlibrary</a>
</li>
<li class="page_item page-item-36">
<a href="/legacy_tables" title="clickable tables">clickable tables</a>
</li>
<li class="page_item page-item-36">
<a href="unified-pseudopotential-format" title="Unified Pseudopotential Format">Unified Pseudopotential Format</a>
</li>
</ul>
</div>
</div>
<div id="content-right" >
<div id="page-title">
<h1>Unified Pseudopotential Format</h1>
</div>
<div class="maincontent">
<p><em>IMPORTANT NOTICE</em>: the UPF format is undergoing a major overhaul and will
become a true XML format with a schema. The format described here may become
obsolete.</p>
<p>The Unified Pseudopotential Format (UPF), currently at v.2.0.1, is designed to store different kinds of pseudopotentials:</p>
<ul>
<li>Norm-conserving (NC) pseudopotentials (PP) in nonlocal form</li>
<li>As above, in both semilocal (SL) and nonlocal (NL) form</li>
<li>Ultrasoft (US) PP (aka Vanderbilt)</li>
<li>Projector Augmented Waves (PAW) datasets</li>
<li>Additional data for the GIPAW reconstruction of all-electron (AE) charge</li>
</ul>
<p>into a flexible file format, having a XML-like structure. The PP is stored in numerical form on a radial grid.</p>
<h2>General structure</h2>
<ul>
<li>
<p>The file is formatted, starts with a line</p>
<pre><UPF version="2.0.1"></pre>
<p>ends with a line</p>
<pre></UPF></pre>
</li>
<li>
<p>The file contains “fields”. A field whose name is “FOO” is delimited by a starting line containing <FOO> and an ending line containing </FOO>, as in the example below:</p>
<pre><FOO>
(content of field FOO)
</FOO></pre>
<p><FOO> and </FOO> are “delimiters” of field FOO.</p>
</li>
<li>
<p>A field may have “attributes”, as in the following example:</p>
<pre><FOO attr="bar">
(content of field FOO)
</FOO></pre>
</li>
<li>
<p>A field name is case-sensitive: use all caps in case of doubt!</p>
</li>
<li>
<p>A field name can contain only letters and digits</p>
</li>
<li>
<p>Spaces are not allowed between the <> and the field name</p>
</li>
<li>
<p>A delimiter need not to start at the beginning of a line</p>
</li>
<li>
<p>Trailing characters in the line after the > of a delimiter are ignored</p>
</li>
<li>
<p>Fields may contain numeric data, character strings, or other fields (“subfields”)</p>
</li>
<li>
<p>Blank lines in a field are ignored.</p>
</li>
<li>
<p>Comments are introduced by <– and terminated by –> (on the same line)</p>
</li>
<li>
<p>Maximum line length is 80 characters</p>
</li>
<li>
<p>Data in fields must be readable using fortran free format</p>
</li>
</ul>
<h2>Defined fields</h2>
<p>First-level fields defined in v.2.0.1 of the UPF format are:</p>
<ul>
<li>PP_INFO</li>
<li>PP_HEADER</li>
<li>PP_MESH</li>
<li>PP_NLCC (optional)</li>
<li>PP_LOCAL</li>
<li>PP_NONLOCAL</li>
<li>PP_SEMILOCAL (optional, only for norm-conserving)</li>
<li>PP_PSWFC (optional)</li>
<li>PP_FULL_WFC (only for PAW)</li>
<li>PP_RHOATOM</li>
<li>PP_PAW (only for PAW)</li>
</ul>
<p>PP_INFO should be the first field, but only for human readers: it is meant to contain info that would allow to generate again the pseudopotential, not data to be read. PP_HEADER must precede PP_MESH that must precede all the remaining fields.</p>
<p>Fields that are not defined are ignored.</p>
<h2>Field specifications</h2>
<p>The meaning of variables is explained at the end of each field. Do loops, indicated by lines containing only dots (…), can be written with as many numbers per line as desired, within the limit of 80 characters per line.</p>
<p>All quantities are in atomic Rydberg units: e^2=2, m=1/2, hbar=1. Lengths are in Bohr (0.529177 A), energies in Ry (13.6058 eV). Potentials are multiplied by e so they have the units of energy. Beware: some quantities computed on the radial mesh contain a multiplying <em>r</em> factor.</p>
<h3>PP_INFO</h3>
<p>The PP_INFO field may contain any piece of information that is deemed useful to re-generate the pseudopotential. The recommanded structure is the following:</p>
<pre><PP_INFO>
Generated using XXX code v.N
Author: Jon Doe
Generation date: 32Oct1976
Pseudopotential type: SL|NC|1/r|US|PAW
Element: Tc
Functional: SLA PW PBX PBC
Suggested minimum cutoff for wavefunctions: N Ry
Suggested minimum cutoff for charge density: M Ry
Non-/scalar-/fully-relativistic pseudopotential
Local potential generation info (L, rcloc, pseudization)
Pseudopotential is spin-orbit/contains GIPAW data
Valence configuration:
nl, pn, l, occ, Rcut, Rcut US, E pseu
els(1), nns(1), lchi(1), oc(1), rcut(1), rcutus(1), epseu(1)
...
els(n), nns(n), lchi(n), oc(n), rcut(n), rcutus(n), epseu(n)
Generation configuration:
as above, including all states used in generation
Pseudization used: Martins-Troullier/RRKJ
<PP_INPUTFILE>
Copy of the input file used in generation
</PP_INPUTFILE>
</PP_INFO>
</pre>
<hr />
<p><strong>PP_HEADER</strong></p>
<p>Structure:</p>
<pre>
<PP\_HEADER attr1="value1" ... attrN="valueN"> ... </PP\_HEADER>
attr value
generated "Generation code"
author "Author"
date "Generation date"
comment "Brief description"
element "Chemical Symbol"
pseudo_type "NC | SL | 1/r | US | PAW"
relativistic "scalar | full | nonrelativistic"
is_ultrasoft .F. | .T.
is_paw .F. | .T.
is_coulomb .F. | .T.
has_so .F. | .T.
has_wfc .F. | .T.
has_gipaw .F. | .T.
paw\_as\_gipaw .F. | .T.
core_correction .F. | .T.
functional "dft"
z_valence Zval
total_psenergy etotps
wfc_cutoff ecutwfc
rho_cutoff ecutrho
l_max lmax
l\_max\_rho lmax_rho
l_local lloc
mesh_size mesh
number\_of\_wfc nwfc
number\_of\_proj nbeta
</PP_HEADER>
</pre>
<hr />
<p>For cases in which different values of attributes are listed, the first one is the default one. Meaning of variables:</p>
<ul>
<li>
<p>“NC” = Norm-Conserving PP, fully nonlocal form only</p>
</li>
<li>
<p>“SL” = Norm-Conserving PP, nonlocal and semilocal forms available</p>
</li>
<li>
<p>“US” = Ultrasoft (Vanderbilt) PP. Implies: is_uspp=.T.</p>
</li>
<li>
<p>“1/r”= Coulomb potential:. Implies: is_coulomb=.T.</p>
</li>
<li>
<p>“PAW”=Projector-Augmented Wave set. Implies: is_paw=.T</p>
</li>
<li>
<p>has_so: fully relativistic PP with spin-orbit terms</p>
</li>
<li>
<p>has_wfc: contains all-electron orbitaks in field PP_FULL_WFC</p>
</li>
<li>
<p>has_gipaw: contains data for GIPAW reconstructions in PP_GIPAW</p>
</li>
<li>
<p>paw_as_gipaw= ?</p>
</li>
<li>
<p>nlcc: non-linear core correction is included</p>
</li>
<li>
<p>“dft”: a label identifying the exchange-correlation functional</p>
</li>
<li>
<p>Zval: valence charge</p>
</li>
<li>
<p>etotps: total pseudo-valence energy of PP</p>
</li>
<li>
<p>ecutwfc: suggested cutoff for orbital expansion into plane waves</p>
</li>
<li>
<p>ecutrho: suggested cutoff for charge density expansion</p>
</li>
<li>
<p>lmax: max angular: momentum component in PP</p>
</li>
<li>
<p>lmax_rho: as above, in atomic charge density (PAW only)</p>
</li>
<li>
<p>mesh: number of points in radial grid</p>
</li>
<li>
<p>wfc: is the number of atomic (pseudo-)orbitals in section PP_PSWFC. May not coincide with the number of atomic states used in PP generation</p>
</li>
<li>
<p>nbeta: number of Kleinman-Bylander projectors (“beta functions”) included in the PP (field PP_NONLOCAL)</p>
</li>
<li>
<ul>
<li>*</li>
</ul>
</li>
</ul>
<p><strong>PP_MESH</strong></p>
<p>Structure:</p>
<pre>
<PP_MESH dx="dx" mesh="mesh" xmin="xmin" rmax="rmax" zmesh="Zmesh">
<PP_R>
r(1) r(2) ... r(mesh)
</PP_R>
<PP_RAB>
rab(1) rab(2) ... rab(mesh)
</PP_RAB>
</PP_MESH>
</pre>
<hr />
<p>Meaning of the variables:</p>
<ul>
<li>dx, mesh, xmi, rmax, Zmesh: radial grid arameters</li>
<li>r (1:mesh): radial grid points (a.u.). Can be one of the following:<em><br />
r(i) = exmin+i*dx/Zmesh</em> or <em>r(i) = (exmin+i*dx-1)/Zmesh</em>, with <em>r(mesh)=rmax</em></li>
<li>rab(mesh): factor required for discrete integration: <em>∫ f(r) dr =∑ifi rabi</em>.</li>
</ul>
<p><strong>PP_NLCC</strong></p>
<pre>
<PP_NLCC>
rho\_atc(1) rho\_atc(2) ... rho_atc(mesh)
</PP_NLCC>
</pre>
<hr />
<p>Optional, needed only for PP with core corrections. Meaning of variables:</p>
<ul>
<li>rho_atc(mesh) : core charge for nonlinear core correction (true charge, not multiplied by 4π<em>r2</em>)</li>
</ul>
<p><strong>PP_SEMILOCAL</strong></p>
<pre>
<PP_SEMILOCAL>
<PP_VNL1 L="l1" J="j1">
V\_1(1) V\_1(2) ... V_1(mesh)
</PP_VNL1 L="l" J="j">
<PP_VNL2 L="l2" J="j2">
V\_2(1) V\_2(2) ... V_2(mesh)
</PP_VNL2 L="l" J="j">
...
</PP_SEMILOCAL>
</pre>
<hr />
<p>Optional, NC PP with semilocal form only. There are nbeta terms, each one containing potential V(1:mesh) for the specified value of l (also j in fully relativistic case only)</p>
<p><strong>PP_LOCAL</strong></p>
<pre>
<PP_LOCAL>
vloc(1) vloc(2) ... vloc(mesh)
</PP_LOCAL>
</pre>
<hr />
<p>vloc(mesh) : local potential (Ry a.u.) sampled on the radial grid</p>
<p><strong>PP_NONLOCAL</strong></p>
<p>Structure:</p>
<pre>
<PP_NONLOCAL>
<PP_BETA>
1 lll(1) "Beta L"
kkbeta(1)
beta(1,1) beta(2,1) ... beta(kkbeta(1),1)
</PP_BETA>
...
<PP_BETA>
nbeta lll(nbeta) "Beta L"
kkbeta(nbeta)
beta(1,nbeta) beta(2,nbeta) ... beta(kkbeta(nbeta),nbeta)
</PP_BETA>
<PP_DIJ>
nd, "Number of nonzero Dij"
do nb=1,nbeta
do mb=nb,nbeta
if (abs (dion (nb, mb) ) > 0) then
nb mb dion(nb,mb) "Q_int"
end if
end do
end do
</PP_DIJ>
<PP_QIJ>
nqf "nqf"
<PP_RINNER>
rinner(0) rinner(1) ... rinner(2*lmax)
</PP_RINNER>
do nb=1,nbeta
do mb=nb,nbeta
nb mb lll(mb) "i j (l)"
qqq(nb,mb) "Q_int"
qfunc(1, nb,mb) qfunc(2, nb,mb) ... qfunc(mesh, nb,mb)
<PP_QFCOEF>
do l=0,2*lmax
do i=1,nqf
qfcoef(i,l,nb,mb)
end do
end do
</PP_QFCOEF>
end do
end do
</PP_QIJ>
</PP_NONLOCAL>
</pre>
<hr />
<ul>
<li>lll(i): angular momentum of projector i</li>
<li>kkbeta(i): number of mesh points for projector i (must be ≤ mesh )</li>
<li>beta(i): projector <em>riβ(ri)</em> (note the factor <em>r</em>!)</li>
<li>dion(i,j): the <em>Dij</em> factors of the nonlocal PP: V<em>NL = </em>∑i,j Di,j |βi><βj_</li>
<li>nqf: number of expansion coefficients for <em>qij</em> (may be zero)</li>
<li>rinner(i): for <em>r < rinner(i)</em> Q functions are pseudized (not read if nqf=0)</li>
<li>qqq(i,j) : <em>Qij = ∫ qij(r) d3r</em></li>
<li>qfunc: <em>r2 qij(r)</em> for <em>r > rinner(i)</em></li>
<li>qfcoef: expansion coefficients of <em>r2 qij(r)</em> for <em>r < rinner(i)</em> (not read if nqf=0)</li>
<li>Note on units: the “beta” and “dion” are defined as in Vanderbilt’s USPP code and should have Bohr^{-1/2} and Ry units, respectively. Since they enter the calculation only as (beta*D*beta), some converters may actually produce “dion” in Ry^{-1} units and “beta” in Ry*Bohr^{-1/2} units instead, as suggested by the Kleinman-Bylader transformation.</li>
</ul>
<p><strong>PP_PSWFC</strong></p>
<pre>
<PP_PSWFC>
els(1) lchi(1) oc(1) "Wavefunction"
chi(1,1) chi(2,1) ... chi(mesh,1)
..........
els(natwfc) lchi(natwfc) oc(natwfc) "Wavefunction"
chi(1,natwfc) chi(2,natwfc) ... chi(mesh,natwfc)
</PP_PSWFC>
</pre>
<hr />
<ul>
<li>chi(mesh,i): <em>χi(r)</em>, <em>i</em>-th radial atomic (pseudo-)orbital (radial part of the KS equation, multiplied by <em>r</em>)</li>
<li>els(natwf), lchi(natwf), oc(natwf): as in PP_HEADER</li>
</ul>
<p><strong>PP_RHOATOM</strong></p>
<pre>
<PP_RHOATOM>
rho\_at(1) rho\_at(2) ... rho_at(mesh)
</PP_RHOATOM>
</pre>
<hr />
<ul>
<li>rho_at(mesh): radial atomic (pseudo-)charge .This is <em>4π r2</em> times the true charge.</li>
</ul>
<h2>Additional Fields</h2>
<h3>PAW</h3>
<p>If a PAW dataset is contained in the UPF file then the additional structure <PAW> is present; it contains the fields listed in the following sections.</p>
<p><strong>PP_PAW_FORMAT_VERSION</strong></p>
<pre>
<PP\_PAW\_FORMAT_VERSION>
version number
</PP\_PAW\_FORMAT_VERSION>
</pre>
<hr />
<p>Contains version number, current version is 0.1.</p>
<p><strong>PP_AUGMENTATION</strong></p>
<pre>
<PP_AUGMENTATION>
Shape of augmentation charge:
BESSEL | GAUSS | PSQ | ...
r\_match\_augfun, irc "augmentation max radius"
lmax_aug "augmentation max angular momentum"
"Augmentation multipoles:"
nb = 1,nbeta
nb1 = 1,nbeta
l = 0,lmax_aug
augmom(nb,nb1,l)
enddo
enddo
enddo
"Augmentation functions:"
do l = 0,lmax_aug
do nb = 1,nbeta
do nb1 = 1,nbeta
if (abs(augmom(nb,nb1,l)) > 0)
"label of the augmentation function"
augfun(k), k = 1, mesh
endif
enddo
enddo
enddo
</PP_AUGMENTATION>
</pre>
<hr />
<p>Data required to build augmentation functions:</p>
<ul>
<li>BESSEL|GAUSS|PSQ|…: function used to pseudize the augmentation</li>
<li>r_match_augfun: range beyond which all the augmentation functions are zero (a.u.)</li>
<li>irc: index of the radial grid closer to, and greater than, r_match_augfun</li>
<li>augmom: multipole of augmentation channel (nb,nb1); it is computed as:<br />
<em>mlnb,nb1 =∫ dr r2 rl (χAEnb(r) χAEnb1(r) – χPSnb(r) χPSnb1(r))</em></li>
<li>augfun: the augmentation function, it is stored only if the corresponding augmentation multipole is different from zero.</li>
</ul>
<p><strong>PP_AE_RHO_ATC</strong></p>
<pre>
<PP\_AE\_RHO_ATC>
aeccharg(k), k = 1,mesh
</PP\_AE\_RHO_ATC>
</pre>
<hr />
<p>All-electron atomic density on the radial grid.</p>
<p><strong>PP_AEWFC</strong></p>
<pre>
<PP_AEWFC>
do nb = 1,nbeta
aewfc(k, nb), k = 1,mesh
end do
</PP_AEWFC>
</pre>
<hr />
<p>All-electron wavefunctions used for the generation of the dataset; there is one wavefunction for each beta projector.</p>
<p><strong>PP_PSWFC_FULL</strong></p>
<pre>
<PP\_PSWFC\_FULL>
do nb = 1,nbeta
pswfc(k, nb), k = 1,mesh
end do
</PP\_PSWFC\_FULL>
</pre>
<hr />
<p>Pseudo wavefunction used for the generation of the dataset; note that in the PP_PSWFC field only the occupied wavefunctions are stored while for PAW you need a wavefunction for each projector.</p>
<p><strong>PP_AEVLOC</strong></p>
<pre>
<PP_AEVLOC>
do nb = 1,nbeta
aewfc(k, nb), k = 1,mesh
end do
</PP_AEVLOC>
</pre>
<hr />
<p>All-electron local potential</p>
<p><strong>PP_KDIFF</strong></p>
<pre>
<PP_KDIFF>
nb = 1,nbeta
nb1 = 1,nbeta
kdiff(nb, nb1)
enddo
enddo
<PP_KDIFF>
</pre>
<hr />
<p>Kinetic energy difference between all-electron and pseudo component of each augmentation channel.</p>
<p><strong>PP_OCCUP</strong></p>
<pre>
<PP_OCCUP>
do nb = 1,nbeta
occ(nb)
end do
</PP_OCCUP>
</pre>
<hr />
<p>Occupations of atomic orbitals.</p>
<p><strong>PP_GRID_RECON</strong></p>
<pre>
<PP\_GRID\_RECON>
"Minimal info to reconstruct the radial grid:"
grid%dx, " dx"
grid%xmin, " xmin"
grid%rmax, " rmax"
grid%zmesh," zmesh"
<PP\_SQRT\_R>
grid%sqr(k), k=1,mesh
</PP\_SQRT\_R>
</PP\_GRID\_RECON>
</pre>
<hr />
<p>Addition data necessary to accurately reconstruct the radial grid used for the dataset generation.</p>
<p><strong>GIPAW</strong></p>
<p>GIPAW additional data is necessary to reconstruct all-electron charge density using the gipaw.x program included in QE distribution.</p>
<p><strong>PP_GIPAW_FORMAT_VERSION</strong></p>
<pre>
<PP\_GIPAW\_FORMAT_VERSION>
version number
</PP\_GIPAW\_FORMAT_VERSION>
</pre>
<hr />
<p>Contains version number, current version is 0.1.</p>
<p><strong>GIPAW_CORE_ORBITALS</strong></p>
<pre>
<GIPAW\_CORE\_ORBITALS>
n\_core\_orbitals "number of core orbitals"
<GIPAW\_CORE\_ORBITAL>
n, l "orbital n and l quantum numbers"
core_orbital(k), k = 1,mesh
</GIPAW\_CORE\_ORBITAL>
</pre>
<p>Repeated for each core orbital</p>
<pre>
</GIPAW\_CORE\_ORBITALS>
</pre>
<hr />
<p>Core orbitals.</p>
<p><strong>GIPAW_LOCAL_DATA</strong></p>
<pre>
<GIPAW\_LOCAL\_DATA>
<GIPAW\_VLOCAL\_AE>
vlocal_ae(k), k = 1,mesh
</GIPAW\_VLOCAL\_AE>
<GIPAW\_VLOCAL\_PS>
vlocal_ps(k), k = 1,mesh
</GIPAW\_VLOCAL\_PS>
</GIPAW\_LOCAL\_DATA>
</pre>
<hr />
<p>All electron and pseudo local potentials, sampled on the radial grid.</p>
<p><strong>GIPAW_ORBITALS</strong></p>
<pre>
<GIPAW_ORBITALS>
<GIPAW\_AE\_ORBITAL>
el(nb), ll(nb)
wfs_ae(k,nb), k = 1, mesh
</GIPAW\_AE\_ORBITAL>
<GIPAW\_PS\_ORBITAL>
rcut(nb), rcutus(nb)
wfs_ae(k,nb), k = 1, mesh
</GIPAW\_PS\_ORBITAL>
Repeated for each valence orbital.
</GIPAW_ORBITALS>
</pre>
<hr />
<ul>
<li>el: principal quantum number (0,1,2..)</li>
<li>ll: angular momentum quantum number</li>
<li>rcut: inner cutof radius (a.u.)</li>
<li>rcutus: outer cutoff radius (a.u.)</li>
<li>wfc_ae: all-electron wavefunction sample on radial grid</li>
<li>wfc_ps: pseudo wavefunction sample on radial grid</li>
</ul>
<p>Last change: July 10, 2018</p>
</div>
</div>
<!-- BEGIN FOOTER -->
<div id="footer">
<div class="footer-logo">
<a href="https://foundation.quantum-espresso.org">
<img src="/user/themes/quantum/images/logo_qef_cropped.jpg" alt="QUANTUM Espresso Foundation" />
</a>
<a href="http://www.max-centre.eu/" align="right">
<img src="/user/themes/quantum/images/MAX_logo.png" alt="Quantum Espresso is a flaghip code of the MaX CoE">
</a>
</div>
<div class="footer-cite">
What I cannot <span class="red">compute</span>, I do not understand (adapted from Richard P. Feynman)<br />
</div>
<div class="footer-text">
</div>
</div>
<div class="footer-logo">
</div>
<!-- END OF FOOTER -->
</div>
</div>
</body>
</html>