Examples sorted by elements
Hydrogen
Elemental:
H2 molecule with film setup
&input film=t /
11.50 0.00 0.00
0.00 11.50 0.00
0.00 0.00 1.00
1.00
1.00 1.00 1.00
2
1 0.0 0.0 0.75
1 0.0 0.0 -0.75
&film dvac=3.72 dtild=4.86/
&kpt nkpt=1 div1=1 div2=1 /
Bond-length plot:
Elemental:
Hydrogen (APW+lo used for actual calculation)
&input film=f /
&lattice latsys="hP" a0=1.8897269 a=4.00759 c=5.00307 /
4
1 1.0 2.0 0.17498
1 2.0 1.0 -0.17498
1 2.0 1.0 -0.32502
1 1.0 2.0 0.32502
&factor 3.0 3.0 1.0 /
&kpt div1=11 div2=11 div3= 9 tkb=0.0005 /
Surfaces: H / Fe (100)
Compounds:
BH3
Diborane molecule with LDA using film setup
&input film=t /
&lattice latsys="tP" a=15.0 c=5 /
8
5 0.0 0.0 1.65000
5 0.0 0.0 -1.65000
1 1.88 0.0 0.0
1 -1.88 0.0 0.0
1 0.0 1.70 3.0
1 0.0 -1.70 3.0
1 0.0 1.70 -3.0
1 0.0 -1.70 -3.0
&factor 15.0 15.0 1.0/
&kpt nkpt=1 div1=1 div2=1 /
Bond-length plot:
Helium
Elemental :
hcp Helium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.92661 c=4.77913 /
2
2 1.0 1.0 1.0
2 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=27 div2=27 div3=15 tkb=0.0005 /
Lithium
Elemental :
Li
Li
&input film=f /
&lattice latsys='rho' a0=14.437589 alpha=23.22001 /
3
3 0. 0. 0.
3 2. 2. 2.
3 7. 7. 7.
&factor 9. 9. 9. /
&kpt div1=17 div2=17 div3=17 tkb=0.0005 /
Beryllium
Elemental :
hcp Be
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.26276 c=3.57316 /
2
4 1.0 1.0 1.0
4 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=27 div2=27 div3=15 tkb=0.0005 /
Boron
Elemental :
B
B (tricl) (APW+lo used for actual calculation)
&input film=f /
&lattice latsys='aP' a0=1.8897269 a=4.90067 b=4.90067 c=5.05098 alpha=119.02035 beta=60.97965 gamma=120 /
12
5.0 0.98969 0.01031 0.67465
5.0 0.01031 0.98969 0.32535
5.0 0.98969 0.65404 0.67465
5.0 0.01031 0.34596 0.32535
5.0 0.34596 0.01031 0.67465
5.0 0.65404 0.98969 0.32535
5.0 0.77885 0.22115 0.07274
5.0 0.22115 0.77885 0.92726
5.0 0.77885 0.63043 0.07274
5.0 0.22115 0.36957 0.92726
5.0 0.36957 0.22115 0.07274
5.0 0.63043 0.77885 0.92726
&kpt div1=10 div2= 10 div3= 6 tkb=0.0005 /
Compounds:
BH3 See Hydrogen
Carbon
Graphite
grahite C (APW+lo used for actual calculation)
&input film=f /
&lattice latsys='hP' a0=1.8897269 a=2.46857 c=8.84079 /
4
6 0.0 0.0 0.25
6 0.0 0.0 -0.25
6 1.0 2.0 0.25
6 2.0 1.0 -0.25
&factor 3.0 3.0 1.0 /
&kpt div1=16 div2=16 div3= 6 tkb=0.0005 /
Diamond
Diamond
&lattice latsys='cF' a=6.74 /
2
6 0.125 0.125 0.125
6 -0.125 -0.125 -0.125
&kpt div1=9 div2=9 div3=9 /
Bandstructure plot:
Nitrogen
Elemental :
cubic nitrogen (APW+lo used for actual calculation)
&input film=f /
&lattice latsys='cP' a0=1.8897269 a=6.19019 /
8
7 0.05189 0.05189 0.05189
7 0.94811 0.94811 0.94811
7 0.44811 0.94811 0.55189
7 0.55189 0.05189 0.44811
7 0.94811 0.55189 0.44811
7 0.05189 0.44811 0.55189
7 0.55189 0.44811 0.94811
7 0.44811 0.55189 0.05189
&kpt div1= 6 div2= 6 div3= 6 tkb=0.0005 /
Oxygen
Elemental :
monoclinic oxygen (introduce magn. moment with the swsp=T switch!) (APW+lo used for actual calculation)
&input film=f /
&lattice latsys='mP' a0=1.8897269 a=4.57148 b=4.27163 c=4.27682 gamma=66.72128 /
4
8.0 0.14667 0.93782 0 0.14667
8.0 0.85333 0.06218 0 0.85333
8.1 0.14667 0.43782 0.50000 0.14667
8.1 0.85333 0.56218 0.50000 0.85333
&kpt div1= 4 div2= 4 div3= 4 tkb=0.0005 /
Compounds:
H2O: See Hydrogen
TiO2: See Ti
ZnO: See Zn
Fluorine
Elemental :
F (mC) (APW+lo used for actual calculation)
&input film=f cartesian=F /
&lattice latsys='mA' a0=1.8897269 a=6.92543 b=6.16310 c=3.68505 gamma=96.77538 /
4
9 0.09413 -0.0753 0.5853
9 0.90587 0.0753 0.4147
9 0.40587 0.4147 0.0753
9 0.59413 -0.4147 0.9247
&kpt div1=10 div2= 10 div3= 6 tkb=0.0005 /
Neon
Elemental :
fcc Ne
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=4.62262 /
1
10 0.0 0.0 0.0
&kpt div1=15 div2=15 div3=15 tkb=0.0005 /
Sodium
Elemental :
Na
Na
&input film=f /
&lattice latsys='rho' a0=17.797693 alpha=22.92203 /
3
11 0. 0. 0.
11 2. 2. 2.
11 7. 7. 7.
&factor 9. 9. 9. /
&kpt div1=15 div2=15 div3=15 tkb=0.0005 /
Magnesium
Elemental :
hcp Mg
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.19405 c=5.17198 /
2
12 1.0 1.0 1.0
12 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=27 div2=27 div3=15 tkb=0.0005 /
Aluminum
Elemental :
fcc aluminum
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=4.04021 /
1
13 0.0 0.0 0.0
&kpt div1=17 div2=17 div3=17 tkb=0.0005 /
Silicon
Elemental :
alpha Si
&input film=f /
&lattice latsys='cF' a0=2.67247737973 a=3.86709 /
2
14 0.125 0.125 0.125
14 -0.125 -0.125 -0.125
&kpt div1=17 div2=17 div3=17 tkb=0.0005 /
Phosphorus
Black Phosphorus
Elemental :
P (oC)
&input film=f cartesian=T /
&lattice latsys='oC' a0=1.8897269 a=3.30466 b=11.31694 c=4.57327 /
4
15 0 0.09357 0.08719
15 0 0.90643 0.91281
15 0 0.40643 0.58719
15 0 0.59357 0.41281
&kpt div1=20 div2= 6 div3=14 tkb=0.0005 /
scaled cartesian units Note, that all coordinates and symmetry operations are in internal units w.r.t. the (orthogonal) axes of the primitive orthorhombic lattice.
P (black) from ICSD database (http://icsd.fiz-karlsruhe.de/)
&input film=f cartesian=t /
&lattice latsys='oC', a0=1.8897269 a=3.3199 b=10.3678 c=4.3154 /
-1
15 0.0000 0.1044 0.0821
&gen 2 ! as in the int. tables
-1 0 0 -0.5 ! for C m c a (64)
0 -1 0 0.5 ! http://www.cryst.ehu.es
0 0 1 0.5 ! but for primitive unit cell
-1 0 0 0
0 -1 0 0
0 0 -1 0 /
The following plot visualizes the bonding in phosphorus:
Sulfur
Elemental :
S
&input film=f /
&lattice latsys='rho' a0=4.8783866 alpha=91.45541 /
1
16 0. 0. 0.
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Chlorine
Elemental :
Cl (o)
&input film=f cartesian=T /
&lattice latsys='oC' a0=1.8897269 a=7.77875 b=4.34941 c=9.02754 /
4
17 0 0.40882 0.89848
17 0 0.09118 0.39848
17 0 0.90882 0.60152
17 0 0.59118 0.10152
&kpt div1= 6 div2= 10 div3= 4 tkb=0.0005 /
Argon
Elemental :
fcc Ar
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=5.95059 /
1
18 0.0 0.0 0.0
&kpt div1=15 div2=15 div3=15 tkb=0.0005 /
Potasium
Elemental :
bcc potassium
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=5.28589 /
1
19 0.0 0.0 0.0
&kpt div1=19 div2=19 div3=19 tkb=0.0005 /
Calcium
Elemental :
fcc Ca
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=5.52507 /
1
20 0.0 0.0 0.0
&kpt div1=21 div2=21 div3=21 tkb=0.0005 /
Scandium
Elemental :
hcp Scandium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.32145 c=5.16192 /
2
21 1.0 1.0 1.0
21 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=27 div2=27 div3=15 tkb=0.0005 /
========== TiInps ============
Titanium
Elemental :
hcp Titanium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.93664 c=4.65193 /
2
22 1.0 1.0 1.0
22 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=31 div2=31 div3=17 tkb=0.0005 /
Alternative:
Ti
&input film=f cartesian=f /
&lattice latsys='hcp' a0=5.57 a=1.0 c=1.588 /
2 ! number of atoms
22 1.000000 2.000000 1.000000 ! Z,positions
22 2.000000 1.000000 -1.000000 ! Z,positions
&factor 3.0 3.0 4.0 /
Bandstructure-plot:
Compounds:
TiO2:
TiO2 Rutile
&input film=f /
&lattice latsys='tP' a0=1.8897269 a=4.594 c=2.959 /
-2
22 0.0 0.0 0.0
8 0.306 0.306 0.0
&gen 3
-1 0 0 0.0
0 -1 0 0.0
0 0 -1 0.0
-1 0 0 0.5
0 1 0 0.5
0 0 -1 0.5
0 -1 0 0.5
1 0 0 0.5
0 0 1 0.5
The structure looks like this:
Vanadium
Elemental :
bcc V
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=2.99894 /
1
23 0.0 0.0 0.0
&kpt div1=31 div2=31 div3=31 tkb=0.0005 /
Chromium
Elemental :
bcc Cr
&input film=f /
&lattice latsys='cP' a0=1.8897269 a=2.87100 /
2
24.0 0.0 0.0 0.0
24.1 0.5 0.5 0.5
&kpt div1=24 div2=24 div3=24 tkb=0.0005 /
Manganese
Elemental :
tet manganese
&input film=f /
&lattice latsys='tP' a0=1.8897269 a=2.5424166 c=3.59552 /
2
25.1 0.0 0.0 0.0
25.2 0.5 0.5 0.5
&kpt div1=28 div2=28 div3=20 tkb=0.0005 /
Iron
Elemental :
bcc Fe
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=2.83351 /
1
26 0.0 0.0 0.0
&kpt div1=27 div2=27 div3=27 tkb=0.0005 /
Elemental (fcc-structure):
Fe fcc in Cu lattice constant
&input film=f /
&lattice latsys='cF' a0=6.82 /
1
26 0.0 0.0 0.0
Cobalt
Elemental :
hcp Cobalt
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.49680 c=4.03081 /
2
27 1.0 1.0 1.0
27 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Nickel
Elemental :
fcc Ni
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=3.52414 /
1
28 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Copper
Elemental :
fcc copper
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=3.63689 /
1
29 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Alternative
Cu bulk
&input film=f /
&lattice latsys='cF' a=6.82 /
1
29 0.0 0.0 0.0
&kpt div1=13 div2=13 div3=13 /
Fermi-surface:
(111)-Surface
Cu (111) 7 layers
&input film=t /
&lattice latsys='hP' a=4.824 c=10.0 /
7
29 0.0 0.0 3.0
29 1.0 2.0 2.0
29 2.0 1.0 1.0
29 0.0 0.0 0.0
29 1.0 2.0 -1.0
29 2.0 1.0 -2.0
29 0.0 0.0 -3.0
&factor 3.0 3.0 0.2538071 /
&kpt div1=13 div2=13 div3=1 /
Bandstructure-plot:
Zinc
Elemental :
hcp Zinc
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.66169 c=5.00397 /
2
30 1.0 1.0 1.0
30 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=17 div2=17 div3=13 tkb=0.0005 /
Alternative
Zn bulk
&input inistop=t oldfleur=f cartesian=f /
&lattice latsys='hcp' a0=5.0267 a=1.0 c=1.856 /
2 ! number of atoms
30 1.000000 2.000000 0.250000 ! Z,positions
30 2.000000 1.000000 -0.250000 ! Z,positions
&factor 3.0 3.0 1.0
Compounds:
ZnO:
ZnO (Wurtzite)
&input inistop=t oldfleur=f cartesian=f /
&lattice latsys='hcp' a=6.14 c=9.828 /
4 ! number of atoms
30 1.000000 2.000000 0.000000 ! Z,positions
30 2.000000 1.000000 0.500000 ! Z,positions
8 1.000000 2.000000 0.382300 ! Z,positions
8 2.000000 1.000000 0.882300 ! Z,positions
&factor 3.0 3.0 1.0 /
Bandstructure-plot:
Gallium
Elemental :
Ga
Gallium(o)
&input film=f cartesian=T /
&lattice latsys='oC' a0=1.8897269 a=4.56689 b=7.75330 c=4.59978 /
4
31 0 0.15642 0.08120
31 0 0.84358 0.91880
31 0 0.34358 0.58120
31 0 0.65642 0.41880
&kpt div1=12 div2= 8 div3=12 tkb=0.0005 /
Germanium
Elemental :
alpha Ge
&input film=f /
&lattice latsys='cF' a0=2.67247737973 a=4.07419 /
2
32 0.125 0.125 0.125
32 -0.125 -0.125 -0.125
&kpt div1=17 div2=17 div3=17 tkb=0.0005 /
Arsenic
Elemental :
As
&input film=f /
&lattice latsys='rho' a0=7.9729138 alpha=53.843117 /
2
33 0.22657 0.22657 0.22657
33 -0.22657 -0.22657 -0.22657
&kpt div1=13 div2=13 div3=13 tkb=0.0005 /
Selenium
Elemental :
Selenium
&input film=f symor=T /
&lattice latsys='hP' a0=1.8897269 a=4.52464 c=5.04978 /
3
34 0.21903 0 0.33333333333
34 0 0.21903 0.66666666667
34 0.78097 0.78097 0.0
&kpt div1= 9 div2= 9 div3= 6 tkb=0.0005 /
Bromium
Elemental :
Bromium(o)
&input film=f cartesian=T /
&lattice latsys='oC' a0=1.8897269 a=8.22860 b=4.22731 c=9.03323 /
4
35 0 0.37574 0.88284
35 0 0.12426 0.38284
35 0 0.87574 0.61716
35 0 0.62426 0.11716
&kpt div1= 4 div2= 8 div3= 4 tkb=0.0005 /
Krypton
Elemental :
fcc Kr
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=6.4293 /
1
36 0.0 0.0 0.0
&kpt div1=15 div2=15 div3=15 tkb=0.0005 /
Rubidium
Elemental :
bcc rubidium
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=5.67197 /
1
37 0.0 0.0 0.0
&kpt div1=19 div2=19 div3=19 tkb=0.0005 /
Strontium
Elemental :
fcc strontium
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=6.01969 /
1
38 0.0 0.0 0.0
&kpt div1=21 div2=21 div3=21 tkb=0.0005 /
Yttrium
Elemental :
hcp Y
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.66051 c=5.67399 /
2
39 1.0 1.0 1.0
39 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Zirconium
Elemental :
hcp Zirconium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.23603 c=5.18953 /
2
40 1.0 1.0 1.0
40 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Niobium
Elemental :
bcc Niobium
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=3.32239 /
1
41 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Molybdenum
Elemental :
bcc Molybdenum
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=3.16930 /
1
42 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Tecnetium
Elemental :
hcp Technetium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.76187 c=4.41905 /
2
43 1.0 1.0 1.0
43 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Ruthenium
Elemental :
hcp Ruthenium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.72661 c=4.30050 /
2
44 1.0 1.0 1.0
44 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Rhodium
Elemental :
fcc rhodium
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=3.84210 /
1
45 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Palladium
Elemental :
fcc palladium
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=3.95315 /
1
46 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Silver
Elemental :
fcc silver
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=4.16424 /
1
47 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Cadmium
Elemental :
hcp Cadmium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.03716 c=5.77061 /
2
48 1.0 1.0 1.0
48 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=17 div2=17 div3=13 tkb=0.0005 /
Indium
Elemental :
bct In
&input film=f /
&lattice latsys='tI' a0=1.8897269 a=3.29841 c=5.06256 /
1
49 0.0 0.0 0.0
&kpt div1=17 div2=17 div3=13 tkb=0.0005 /
Tin
Elemental :
alpha Sn
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=6.6559537 /
2
50 0.125 0.125 0.125
50 -0.125 -0.125 -0.125
&kpt div1=17 div2=17 div3=17 tkb=0.0005 /
Antimony
Elemental :
Sb
&input film=f /
&lattice latsys='rho' a0=8.6669374 alpha=57.051774 /
2
51 0.23332 0.23332 0.23332
51 -0.23332 -0.23332 -0.23332
&kpt div1=13 div2=13 div3=13 tkb=0.0005 /
Tellurium
Elemental :
hex Te
&input film=f symor=T /
&lattice latsys='hP' a0=1.8897269 a=4.50788 c=5.96346 /
3
52 0.26966 0 0.33333333333
52 0 0.26966 0.66666666667
52 0.73034 0.73034 0.0
&kpt div1= 9 div2= 9 div3= 5 tkb=0.0005 /
Iodine
Elemental :
Iodine (o)
&input film=f /
&lattice latsys='oP' a0=1.8897269 a=8.57701 b=4.55422 c=10.18381 /
8
53 0 0.36410 0.87857
53 0 0.13590 0.37857
53 0 0.86410 0.62143
53 0 0.63590 0.12143
53 0.50000 0.86410 0.87857
53 0.50000 0.63590 0.37857
53 0.50000 0.36410 0.62143
53 0.50000 0.13590 0.12143
&kpt div1= 8 div2=10 div3= 6 tkb=0.0005 /
Xenon
Elemental :
fcc Xe
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=7.05482 /
1
54 0.0 0.0 0.0
&kpt div1=15 div2=15 div3=15 tkb=0.0005 /
Caesium
Elemental :
bcc Cs
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=6.16153 /
1
55 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Barium
Elemental :
bcc barium
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=5.02792 /
1
56 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Hafnium
Elemental :
hcp Hafnium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.20273 c=5.06493 /
2
72 1.0 1.0 1.0
72 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Tantalum
Elemental :
bcc tantalum
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=3.32169 /
1
73 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Compounds:
Ta2D:
Ta2D from ICSD database (http://icsd.fiz-karlsruhe.de/) #61486
&input film=f cartesian=t /
&lattice latsys='oA', a0=1.8897269 a=4.67 b=4.67 c=3.303 /
-2
73 0.2624 0.25 0.25 ! A-centered setting like in #28022
1 0.000 0.00 0.00 ! use H for D
&gen 2 ! as in the int. tables
-1 0 0 0.0 ! for C 2 2 2
0 -1 0 0.0 ! http://www.cryst.ehu.es
0 0 1 0.0
-1 0 0 0.0
0 1 0 0.0
0 0 -1 0.0 / ! we leave out the translation
This is a simple example of an A-centered orthorhombic lattice.
Tungsten
Elemental :
bcc tungsten
&input film=f /
&lattice latsys='cI' a0=1.8897269 a=3.18968 /
1
74 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Rhenium
Elemental :
hcp Re
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.77447 c=4.47993 /
2
75 1.0 1.0 1.0
75 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Osmium
Elemental :
hcp Osmium
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=2.75895 c=4.35694 /
2
76 1.0 1.0 1.0
76 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Iridium
Elemental :
fcc iridium
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=3.87675 /
1
77 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Platinum
Elemental :
fcc platinum
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=3.97675 /
1
78 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Gold
Elemental :
fcc gold
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=4.17410 /
1
79 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Mercury
Elemental :
bct Hg
&input film=f /
&lattice latsys='tI' a0=1.8897269 a=4.10305 c=3.54831 /
1
80 0.0 0.0 0.0
&kpt div1=19 div2=19 div3=19 tkb=0.0005 /
Thallium
Elemental :
hcp Tl
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.59255 c=5.64195 /
2
81 1.0 1.0 1.0
81 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=19 div2=19 div3=13 tkb=0.0005 /
Lead
Elemental :
fcc Pb
&input film=f /
&lattice latsys='cF' a0=1.8897269 a=5.04386 /
1
82 0.0 0.0 0.0
&kpt div1=25 div2=25 div3=25 tkb=0.0005 /
Bismuth
Elemental :
Bi
rho Bi
&input film=f /
&lattice latsys='rho' a0=9.1623962 alpha=56.484592 /
2
83 0.23374 0.23374 0.23374
83 -0.23374 -0.23374 -0.23374
&kpt div1=13 div2=13 div3=13 tkb=0.0005 /
Elemental:
Bi rhl {exp lattice constant, 2 atom unit cell}
&input inistop=t oldfleur=f /
&lattice latsys='rho' a0=8.976, alpha=57.23 /
2 ! number of atoms
83 0.263941 0.263941 0.263941 ! Z,positions
83 -0.263941 -0.263941 -0.263941 ! Z,positions
&soc 0.0 0.0
Bandstructure:
Polonium
Elemental :
cub Po
&input film=f /
&lattice latsys='cP' a0=1.8897269 a=3.34818 /
1
84 0.0 0.0 0.0
&kpt div1=24 div2=24 div3=24 /
Ytterbium
Elemental input:
Yb fcc bulk
&lattice latsys='cF' a0=1.8897269 a=5.49 /
1
70 0.0 0.0 0.0
We treat the 5s and 5p states as local orbitals and the 4f with LDA+U
density matrix file: n_mmp_mat
1.9703395169836 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 -0.0014408308808 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 1.9678035740164 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000 0.0043960476423 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 1.9710835588552 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 -0.0014408308808 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 1.9692234541878
0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000
-0.0014408308808 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 1.9710835588552 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0043960476423 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000 1.9678035740164 0.0000000000000 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 -0.0014408308808 0.0000000000000 0.0000000000000
0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 0.0000000000000 1.9703395169836 0.0000000000000
In this way the lattice constant is approximated quite reasonably:
Density of states:
Lutetium
Elemental:
hcp Lu
&input film=f /
&lattice latsys='hdp' a0=1.8897269 a=3.52388 c=5.48025 /
2
71 1.0 1.0 1.0
71 -1.0 -1.0 -1.0
&factor 3.0 3.0 4.0 /
&kpt div1=25 div2=25 div3=15 tkb=0.0005 /
Proactinium
Elemental:
Pa fcc bulk
&input film=f /
&lattice latsys='cF' a=8.7643 /
1
91 0.0 0.0 0.0
&kpt div1=13 div2=13 div3=13 /