[QE-users] normalization of wfcs (wfck2r)
Christoph Wolf
wolf.christoph at qns.science
Thu Aug 26 06:24:40 CEST 2021
Dear all,
I am converting a wave-function from k to R space using wfck2r.x. There is
only 1 k-point (particle in a box of 15*15*15 Ang) and the wfc looks like:
==== e( 8) = -4.84501 eV ====
psi = 0.998*[# 2]
|psi|^2 = 0.999
state # 2: atom 1 (Fe ), wfc 2 (l=0 m= 1)
So this is a "pure s state"; i converte the wfc into real space by using:
first_k = 1,
last_k = 1,
first_band = 8
last_band = 8
and I get an ".oct" (plain text) file containing on a grid of 144 x 144
x144:
(* -0.885667802907E+01*, 0.607695317678E-16)
( *-0.627285985227E+01*, -0.418670474801E-15)
( -0.291716024900E+01, -0.436306625693E-15)
....
..
if I understood the manual correctly the order is such that the real part
for example:
re[:,:,0]= [ *-8.85667803, * *-6.02698044, * -2.39305522, ...,
-12.45808728, -11.95557896, -10.79190254]
that is the first 144*144 elements belong to z=0, the next to z=1, z=2,
...z=144.
my volume element would be: (15/0.529)**3/144**3=0.0076 (au^3/grid^3)
if I take the dot product of the complex wfc
<psi*|psi>*0.076=25733.51243066387+0j, but should this not be the same as
|psi|^2 = 0.999?
This is probably a very stupid and obvious mistake but I have been looking
at this for a while now and don't get it...
Best,
Chris
--
Group Leader "Theory of Quantum Systems at Surfaces"
IBS Center for Quantum Nanoscience
Seoul, South Korea
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.quantum-espresso.org/pipermail/users/attachments/20210826/a3a5fb5c/attachment.html>
More information about the users
mailing list