I have a question regarding the partial density of state information obtained from NWChem. When I do a calculation with a material involving partially filled f shells (so 5f U), the density of states output includes s, p, and d densities for alpha and beta electrons (great), but the f state occupation is reported as being uniformly zero. (This is using the newer 14-electron U psp, not the older 6-electron one which did not have 5f electrons.) Physically, I expect that a great deal of the electronic DOS should be due to the 5f electrons, and this is borne out by other researchers. Am I doing something incorrect, or is this a latent shortcoming of the pseudopotential or routine? I have tried this on two different machines with 6.0 and 6.1, so I don't think it's a compilation problem (although I can't a priori rule that out).
My input is:
title "Gamma-uranium 2-atom bcc cell - electronic structure"
echo
permanent_dir .
scratch_dir .
start gamma
memory 1500 mb
geometry center noautosym noautoz print
system crystal
lat_a 3.506
lat_b 3.506
lat_c 3.506
alpha 90.0d0
beta 90.0d0
gamma 90.0d0
end
U1 0 0 0
U2 0.5 0.5 0.5
end
nwpw
simulation_cell
ngrid 64 64 64
end
# Spin multiplicity: 2 * 5 (U) = 10
mult 10
monkhorst-pack 5 5 5
ewald_rcut 3.0
ewald_ncut 8
xc hse
cutoff 120.0
#mapping 2
end
set nwpw:minimizer 2 # Grassman LMBFGS minimizer
set nwpw:cif_filename gamma.opt # create a CIF file containing optimization history
set nwpw:lcao_skip .true.
set nwpw:kbpp_ray .true.
set nwpw:kbpp_filter .true.
set includestress .true.
set includelattice .true.
# Get the electronic characteristics of this geometry.
nwpw
dos
efield
mulliken
end
task pspw energy
task band energy
nwpw
virtual 16
brillouin_zone
zone_name bccpath
path bcc gamma h n gamma p h n
end
zone_structure_name bccpath
end
task band structure
nwpw
virtual 16
dos-grid 5 5 5
end
task band dos
Typical output is of the following types, which mostly seem to make sense.
gamma.alpha_band:
0.000000E+00 -0.771206E+00 -0.728674E+00 0.115314E+00 0.115562E+00 0.115594E+00 0.237159E+
00 0.237443E+00 0.237483E+00 0.652761E+00 0.673993E+00 0.674313E+00 0.817007E+00 0.817133
E+00 0.817151E+00 0.819150E+00 0.864800E+00 0.864854E+00 0.865273E+00 0.905204E+00 0.9052
25E+00 0.905514E+00 0.906242E+00 0.906313E+00 0.906887E+00 0.926702E+00 0.926756E+00 0.92
6771E+00 0.926784E+00 0.104918E+01 0.104931E+01 0.126065E+01 0.126068E+01 0.126119E+01 0.
136659E+01
0.483115E-01 -0.770907E+00 -0.728905E+00 0.115004E+00 0.115198E+00 0.117224E+00 0.233673E+
00 0.236920E+00 0.237128E+00 0.656933E+00 0.675980E+00 0.676297E+00 0.815910E+00 0.816007
E+00 0.817525E+00 0.819236E+00 0.863958E+00 0.865117E+00 0.865453E+00 0.894464E+00 0.9007
58E+00 0.901096E+00 0.910887E+00 0.911125E+00 0.912313E+00 0.926518E+00 0.926593E+00 0.92
6656E+00 0.935386E+00 0.104965E+01 0.104979E+01 0.123823E+01 0.125614E+01 0.125646E+01 0.
137268E+01
...
gamma.DOS_alpha_all:
-0.84830E+00 0.91855E-01
-0.84461E+00 0.97942E-01
-0.84091E+00 0.10471E+00
-0.83722E+00 0.11226E+00
-0.83352E+00 0.12074E+00
...
gamma.DOS_alpha_d:
-0.84830E+00 0.13291E-02
-0.84461E+00 0.13352E-02
-0.84091E+00 0.13414E-02
-0.83722E+00 0.13476E-02
...
gamma.DOS_alpha_f:
-0.84830E+00 0.00000E+00
-0.84461E+00 0.00000E+00
-0.84091E+00 0.00000E+00
-0.83722E+00 0.00000E+00
...
0.99184E+00 0.00000E+00
0.99554E+00 0.00000E+00
(Beta results are similar but opposite in sign, if appropriate.)
|