quantum-espresso/KCW/examples/example06/reference/SPIN_DOWN/cri3.kcw-wann2kcw.out

434 lines
22 KiB
Plaintext

=--------------------------------------------------------------------------------=
::: ::: :::::::: ::: :::
:+: :+: :+: :+: :+: :+:
+:+ +:+ +:+ +:+ +:+
+#++:++ +#+ +#+ +:+ +#+
+#+ +#+ +#+ +#+ +#+#+ +#+
#+# #+# #+# #+# #+#+# #+#+#
### ### ######## ### ###
Koopmans functional implementation based on DFPT; please cite this program as
N.Colonna, R. De Gennaro, E. Linscott, and N. Marzari, JCTC 18, 5435 (2022)
If you use the non-collinear mode (with/without spin-orbit coupling) please cite
A. Marrazzo and N. Colonna, arXiv:2402.14575 (2024)
=--------------------------------------------------------------------------------=
Program KCW v.7.3.1 starts on 5Jun2024 at 10:53:56
This program is part of the open-source Quantum ESPRESSO suite
for quantum simulation of materials; please cite
"P. Giannozzi et al., J. Phys.:Condens. Matter 21 395502 (2009);
"P. Giannozzi et al., J. Phys.:Condens. Matter 29 465901 (2017);
"P. Giannozzi et al., J. Chem. Phys. 152 154105 (2020);
URL http://www.quantum-espresso.org",
in publications or presentations arising from this work. More details at
http://www.quantum-espresso.org/quote
Parallel version (MPI), running on 8 processors
MPI processes distributed on 1 nodes
K-points division: npool = 4
R & G space division: proc/nbgrp/npool/nimage = 2
240404 MiB available memory on the printing compute node when the environment starts
KCW INPUT SUMMARY
============================================
CONTROL
# title = W2K Silicon
# out_dir =/home/colonna_n/CO
# prefix = cri3
# calculation = wann2kcw
# kcw_iverbosity = 2
# kcw_at_ks = F
# MP grid = 2 2 2
# spin_component = 2
# homo_only = F
# read_unitary_matrix = T
# check_ks = T
# l_vcut = F
# assume_isolated = none
# io_sp = F
# io_real_space = F
WANNIER
# seedname = wann
# num_wann_occ = 32
# num_wann_emp = 10
# have_empty = T
# has_disentangle = F
# l_unique_manifold = F
============================================
INFO: Reading pwscf data
Reading xml data from directory:
/home/colonna_n/CODES/q-e-kcw-soc/tempdir/cri3.save/
IMPORTANT: XC functional enforced from input :
Exchange-correlation= PW
( 1 4 0 0 0 0 0)
Any further DFT definition will be discarded
Please, verify this is what you really want
Parallelization info
--------------------
sticks: dense smooth PW G-vecs: dense smooth PW
Min 1079 1079 295 31069 31069 4530
Max 1080 1080 296 31070 31070 4533
Sum 2159 2159 591 62139 62139 9063
Using Slab Decomposition
Check: negative core charge= -0.000182
Reading collected, re-writing distributed wavefunctions
REPORT # of electrons
nelec= 70.00000000
nelup= 38.00000000
neldw= 32.00000000
nkstot= 16
nspin= 2
INFO: Buffer for KS wfcs, OPENED
INFO: Optimal Matrix READ
INFO: total number of Wannier functions 42
INFO: Unitary matrix, READ from file
INFO: Buffer for WFs, OPENED
INFO: Buffer for WFs ALL-k, OPENED
INFO: Minimizing orbitals from Unitary Matrix Rotation
INFO: Performing a check on the eigenvalues of the rotated KS Hamilotnian ...
INFO: Performing a check on the eigenvalues of the rotated KS Hamiltonian ... DONE
INFO: Minimizing orbitals DEFINED
INFO: Coulomb q+G=0 treatment:
INFO: Divergence none
INFO: q-grid dimension 2 2 2
INFO: cell volume 1819.528131144983
INFO: Gamma Extrapolation F
INFO: extrapolation q->0 term not estimated
INFO: Bare Coulomb q+G=0 0.00000E+00
INFO: PREPARING THE KCW CALCULATION ...
INFO: Compute Wannier-orbital Densities ...
==============================================================================
iq = 1
The Wannier density at q = 0.0000000 0.0000000 0.0000000 [Cart ]
The Wannier density at q = 0.0000000 0.0000000 0.0000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 1 1 ) 1 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 1 2 ) 2 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 1 3 ) 3 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 1 4 ) 4 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 1 5 ) 5 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 1 6 ) 6 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 1 7 ) 7 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 1 8 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 1 Structure Factor S(q) [Re, Im] = 8.00000000 0.00000000
==============================================================================
iq = 2
The Wannier density at q = 0.0000000 -0.6474914 0.1944075 [Cart ]
The Wannier density at q = -0.0000000 -0.0000000 0.5000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 2 1 ) 2 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 2 2 ) 1 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 2 3 ) 4 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 2 4 ) 3 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 2 5 ) 6 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 2 6 ) 5 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 2 7 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 2 8 ) 7 + 0.0000 0.0000 1.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 2 Structure Factor S(q) [Re, Im] = 0.00000000 0.00000000
==============================================================================
iq = 3
The Wannier density at q = -0.5607440 0.3237457 0.1944075 [Cart ]
The Wannier density at q = 0.0000000 0.5000000 -0.0000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 3 1 ) 3 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 3 2 ) 4 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 3 3 ) 1 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 3 4 ) 2 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 3 5 ) 7 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 3 6 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 3 7 ) 5 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 3 8 ) 6 + 0.0000 1.0000 0.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 3 Structure Factor S(q) [Re, Im] = 0.00000000 0.00000000
==============================================================================
iq = 4
The Wannier density at q = -0.5607440 -0.3237457 0.3888149 [Cart ]
The Wannier density at q = -0.0000000 0.5000000 0.5000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 4 1 ) 4 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 4 2 ) 3 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 4 3 ) 2 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 4 4 ) 1 + 0.0000 1.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 4 5 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 4 6 ) 7 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 4 7 ) 6 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 4 8 ) 5 + 0.0000 1.0000 1.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 4 Structure Factor S(q) [Re, Im] = 0.00000000 0.00000000
==============================================================================
iq = 5
The Wannier density at q = 0.5607440 0.3237457 0.1944075 [Cart ]
The Wannier density at q = 0.5000000 0.0000000 -0.0000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 5 1 ) 5 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 5 2 ) 6 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 5 3 ) 7 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 5 4 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 5 5 ) 1 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 5 6 ) 2 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 5 7 ) 3 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 5 8 ) 4 + 1.0000 0.0000 0.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 5 Structure Factor S(q) [Re, Im] = 0.00000000 0.00000000
==============================================================================
iq = 6
The Wannier density at q = 0.5607440 -0.3237457 0.3888149 [Cart ]
The Wannier density at q = 0.5000000 -0.0000000 0.5000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 6 1 ) 6 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 6 2 ) 5 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 6 3 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 6 4 ) 7 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 6 5 ) 2 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 6 6 ) 1 + 1.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 6 7 ) 4 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 6 8 ) 3 + 1.0000 0.0000 1.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 6 Structure Factor S(q) [Re, Im] = 0.00000000 0.00000000
==============================================================================
iq = 7
The Wannier density at q = 0.0000000 0.6474914 0.3888149 [Cart ]
The Wannier density at q = 0.5000000 0.5000000 -0.0000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 7 1 ) 7 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 7 2 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 7 3 ) 5 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 7 4 ) 6 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 7 5 ) 3 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 7 6 ) 4 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 7 7 ) 1 + 1.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 7 8 ) 2 + 1.0000 1.0000 0.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 7 Structure Factor S(q) [Re, Im] = 0.00000000 -0.00000000
==============================================================================
iq = 8
The Wannier density at q = 0.0000000 0.0000000 0.5832224 [Cart ]
The Wannier density at q = 0.5000000 0.5000000 0.5000000 [Cryst]
==============================================================================
The map (iq,ik) --> ip + G ( 8 1 ) 8 + 0.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 8 2 ) 7 + 0.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 8 3 ) 6 + 0.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 8 4 ) 5 + 0.0000 1.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 8 5 ) 4 + 1.0000 0.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 8 6 ) 3 + 1.0000 0.0000 1.0000 [Cryst]
The map (iq,ik) --> ip + G ( 8 7 ) 2 + 1.0000 1.0000 0.0000 [Cryst]
The map (iq,ik) --> ip + G ( 8 8 ) 1 + 1.0000 1.0000 1.0000 [Cryst]
INFO: Map k+q -> p in 1BZ DONE
INFO: rho_q(r) DONE
INFO: iq = 8 Structure Factor S(q) [Re, Im] = 0.00000000 0.00000000
INFO: Wannier density number
iwann= 1 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 2 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 3 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 4 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 5 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 6 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 7 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 8 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 9 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 10 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 11 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 12 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 13 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 -0.00000000
iwann= 14 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 15 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 16 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 17 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 18 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 19 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 20 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 21 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 22 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 23 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 24 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 25 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 26 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 27 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 -0.00000000
iwann= 28 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 -0.00000000
iwann= 29 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 -0.00000000
iwann= 30 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 31 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 32 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 33 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 34 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 35 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 36 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 37 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 38 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 39 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 40 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 41 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
iwann= 42 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.00000000 0.00000000
INFO: Orbital Self-Hartree (SH)
orb 1 SH 0.899339
orb 2 SH 0.899339
orb 3 SH 0.905985
orb 4 SH 0.906151
orb 5 SH 0.906161
orb 6 SH 0.905985
orb 7 SH 0.906151
orb 8 SH 0.906161
orb 9 SH 0.329559
orb 10 SH 0.329559
orb 11 SH 0.329559
orb 12 SH 0.329559
orb 13 SH 0.329559
orb 14 SH 0.329559
orb 15 SH 0.252512
orb 16 SH 0.257142
orb 17 SH 0.256175
orb 18 SH 0.257095
orb 19 SH 0.256262
orb 20 SH 0.252488
orb 21 SH 0.252514
orb 22 SH 0.256167
orb 23 SH 0.257146
orb 24 SH 0.252512
orb 25 SH 0.257142
orb 26 SH 0.256175
orb 27 SH 0.257095
orb 28 SH 0.256262
orb 29 SH 0.252488
orb 30 SH 0.252514
orb 31 SH 0.256167
orb 32 SH 0.257146
orb 33 SH 0.367563
orb 34 SH 0.367565
orb 35 SH 0.367563
orb 36 SH 0.367559
orb 37 SH 0.367562
orb 38 SH 0.367560
orb 39 SH 0.290361
orb 40 SH 0.290362
orb 41 SH 0.290362
orb 42 SH 0.290361
INFO: PREPARING THE KCW CALCULATION ... DONE
KCW : 9.08s CPU 10.46s WALL
INITIALIZATION:
map : 0.00s CPU 0.01s WALL ( 8 calls)
rho_of_q : 1.82s CPU 2.73s WALL ( 8 calls)
init_vloc : 0.01s CPU 0.01s WALL ( 1 calls)
init_us_1 : 0.00s CPU 0.00s WALL ( 1 calls)
General routines
fft : 0.00s CPU 0.01s WALL ( 5 calls)
ffts : 2.97s CPU 2.98s WALL ( 2352 calls)
davcio : 0.00s CPU 0.02s WALL ( 12 calls)
KCW : 9.08s CPU 10.47s WALL
This run was terminated on: 10:54: 7 5Jun2024
=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=