scnc
/
quantum-espresso
mirror of https://gitlab.com/QEF/q-e.git
1
0
Fork 0
Code Issues Packages Projects Releases Wiki Activity
quantum-espresso/test-suite/kcw_wann-nc/benchmark.out.git.inp=Si_nc...

555 lines
30 KiB
Plaintext
Raw Permalink Blame History

=--------------------------------------------------------------------------------=
::: ::: :::::::: ::: :::
:+: :+: :+: :+: :+: :+:
+:+ +:+ +:+ +:+ +:+
+#++:++ +#+ +#+ +:+ +#+
+#+ +#+ +#+ +#+ +#+#+ +#+
#+# #+# #+# #+# #+#+# #+#+#
### ### ######## ### ###
Koopmans functional implementation based on DFPT; please cite this program as
N.Colonna, R. De Gannaro, E. Linscott, and N. Marzari, JCTC 18, 5435 (2022)
=--------------------------------------------------------------------------------=
Program KCW v.7.1 starts on 16Feb2024 at 15:28:26
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 1 processors
MPI processes distributed on 1 nodes
5015 MiB available memory on the printing compute node when the environment starts
KCW INPUT SUMMARY
============================================
CONTROL
# title = KI Silicon
# out_dir = ./out/
# prefix = Si_nc
# calculation = wann2kcw
# kcw_iverbosity = 2
# kcw_at_ks = F
# MP grid = 2 2 2
# spin_component = 1
# homo_only = F
# read_unitary_matrix = T
# check_ks = T
# l_vcut = F
# assume_isolated = none
WANNIER
# seedname = Si_nc
# num_wann_occ = 8
# have_empty = F
# has_disentangle = F
# l_unique_manifold = F
============================================
INFO: Reading pwscf data
Reading xml data from directory:
./out/Si_nc.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
G-vector sticks info
--------------------
sticks: dense smooth PW G-vecs: dense smooth PW
Sum 295 295 91 3287 3287 609
Using Slab Decomposition
Reading collected, re-writing distributed wavefunctions
Message from routine kcw_readin:
Non-collinear KCW calculation.
INFO: Buffer for KS wfcs, OPENED
INFO: total number of Wannier functions 8
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 ...
Rotated Occupation Matrix (ROM) ik= 1 xk = 0.0000 0.0000 0.0000
1.0000 -0.0000 0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000
-0.0000 1.0000 -0.0000 -0.0000 0.0000 0.0000 -0.0000 0.0000
0.0000 -0.0000 1.0000 -0.0000 0.0000 0.0000 -0.0000 -0.0000
0.0000 -0.0000 -0.0000 1.0000 0.0000 0.0000 -0.0000 0.0000
-0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 0.0000 0.0000 1.0000 -0.0000 0.0000
-0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 1.0000 0.0000
-0.0000 0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 1.0000
Trace 8.000000000176044
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= 0.0000 0.0000 0.0000 ... DONE
WANN -5.7047 -5.7047 6.3115 6.3115 6.3611 6.3611 6.3611 6.3611
PWSCF -5.7047 -5.7047 6.3115 6.3115 6.3611 6.3611 6.3611 6.3611
Rotated Occupation Matrix (ROM) ik= 2 xk = -0.5000 0.5000 -0.5000
1.0000 0.0000 0.0000 -0.0000 0.0000 0.0000 -0.0000 -0.0000
0.0000 1.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000
0.0000 -0.0000 1.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 1.0000 -0.0000 0.0000 0.0000 0.0000
0.0000 0.0000 0.0000 -0.0000 1.0000 -0.0000 0.0000 -0.0000
0.0000 0.0000 -0.0000 0.0000 -0.0000 1.0000 0.0000 -0.0000
-0.0000 0.0000 -0.0000 0.0000 0.0000 0.0000 1.0000 -0.0000
-0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 -0.0000 1.0000
Trace 7.999999999998614
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= -0.5000 0.5000 -0.5000 ... DONE
WANN -3.3363 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
PWSCF -3.3363 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
Rotated Occupation Matrix (ROM) ik= 3 xk = 0.5000 0.5000 0.5000
1.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 0.0000
-0.0000 1.0000 0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000
-0.0000 0.0000 1.0000 0.0000 0.0000 0.0000 -0.0000 -0.0000
-0.0000 -0.0000 0.0000 1.0000 -0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 0.0000 -0.0000 1.0000 0.0000 0.0000 0.0000
-0.0000 -0.0000 0.0000 0.0000 0.0000 1.0000 0.0000 -0.0000
0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000 1.0000 -0.0000
0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 1.0000
Trace 7.999999999961760
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= 0.5000 0.5000 0.5000 ... DONE
WANN -3.3363 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
PWSCF -3.3363 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
Rotated Occupation Matrix (ROM) ik= 4 xk = 0.0000 1.0000 0.0000
1.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 0.0000
-0.0000 1.0000 0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000
-0.0000 0.0000 1.0000 -0.0000 0.0000 0.0000 -0.0000 0.0000
-0.0000 0.0000 -0.0000 1.0000 -0.0000 -0.0000 0.0000 -0.0000
-0.0000 -0.0000 0.0000 -0.0000 1.0000 -0.0000 0.0000 -0.0000
-0.0000 -0.0000 0.0000 -0.0000 -0.0000 1.0000 0.0000 -0.0000
0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 0.0000 0.0000 -0.0000 -0.0000 -0.0000 0.0000 1.0000
Trace 8.000000000158813
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= 0.0000 1.0000 0.0000 ... DONE
WANN -1.5308 -1.5308 -1.5308 -1.5308 3.4058 3.4058 3.4058 3.4058
PWSCF -1.5309 -1.5309 -1.5308 -1.5308 3.4057 3.4057 3.4058 3.4058
Rotated Occupation Matrix (ROM) ik= 5 xk = -0.5000 -0.5000 0.5000
1.0000 0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000
0.0000 1.0000 0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000
0.0000 0.0000 1.0000 0.0000 -0.0000 -0.0000 -0.0000 -0.0000
-0.0000 0.0000 0.0000 1.0000 -0.0000 -0.0000 -0.0000 0.0000
-0.0000 -0.0000 -0.0000 -0.0000 1.0000 0.0000 0.0000 0.0000
0.0000 0.0000 -0.0000 -0.0000 0.0000 1.0000 -0.0000 0.0000
-0.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000 1.0000 0.0000
-0.0000 -0.0000 -0.0000 0.0000 0.0000 0.0000 0.0000 1.0000
Trace 8.000000000063805
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= -0.5000 -0.5000 0.5000 ... DONE
WANN -3.3363 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
PWSCF -3.3363 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
Rotated Occupation Matrix (ROM) ik= 6 xk = -1.0000 0.0000 0.0000
1.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 0.0000
-0.0000 1.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 -0.0000
-0.0000 -0.0000 1.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
-0.0000 -0.0000 -0.0000 1.0000 0.0000 0.0000 0.0000 0.0000
-0.0000 0.0000 -0.0000 0.0000 1.0000 0.0000 0.0000 -0.0000
-0.0000 -0.0000 -0.0000 0.0000 0.0000 1.0000 0.0000 -0.0000
0.0000 0.0000 -0.0000 0.0000 0.0000 0.0000 1.0000 0.0000
0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 0.0000 1.0000
Trace 8.000000000104521
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= -1.0000 0.0000 0.0000 ... DONE
WANN -1.5308 -1.5308 -1.5308 -1.5308 3.4058 3.4058 3.4058 3.4058
PWSCF -1.5309 -1.5309 -1.5308 -1.5308 3.4057 3.4057 3.4058 3.4058
Rotated Occupation Matrix (ROM) ik= 7 xk = 0.0000 0.0000 1.0000
1.0000 -0.0000 0.0000 0.0000 0.0000 -0.0000 0.0000 0.0000
-0.0000 1.0000 0.0000 0.0000 -0.0000 -0.0000 0.0000 0.0000
0.0000 0.0000 1.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000
0.0000 0.0000 -0.0000 1.0000 -0.0000 -0.0000 0.0000 -0.0000
0.0000 -0.0000 -0.0000 -0.0000 1.0000 -0.0000 -0.0000 0.0000
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 1.0000 -0.0000 -0.0000
0.0000 0.0000 -0.0000 0.0000 -0.0000 -0.0000 1.0000 -0.0000
0.0000 0.0000 -0.0000 -0.0000 0.0000 -0.0000 -0.0000 1.0000
Trace 8.000000000233962
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= 0.0000 0.0000 1.0000 ... DONE
WANN -1.5308 -1.5308 -1.5308 -1.5308 3.4058 3.4058 3.4058 3.4058
PWSCF -1.5309 -1.5309 -1.5308 -1.5308 3.4057 3.4057 3.4058 3.4058
Rotated Occupation Matrix (ROM) ik= 8 xk = -0.5000 0.5000 0.5000
1.0000 -0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000 -0.0000
-0.0000 1.0000 0.0000 0.0000 0.0000 0.0000 -0.0000 -0.0000
-0.0000 0.0000 1.0000 -0.0000 -0.0000 -0.0000 0.0000 -0.0000
0.0000 0.0000 -0.0000 1.0000 -0.0000 -0.0000 -0.0000 -0.0000
-0.0000 0.0000 -0.0000 -0.0000 1.0000 -0.0000 0.0000 -0.0000
0.0000 0.0000 -0.0000 -0.0000 -0.0000 1.0000 0.0000 -0.0000
0.0000 -0.0000 0.0000 -0.0000 0.0000 0.0000 1.0000 0.0000
-0.0000 -0.0000 -0.0000 -0.0000 -0.0000 -0.0000 0.0000 1.0000
Trace 7.999999999968267
ROM eig 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000
KS Hamiltonian calculation at k= -0.5000 0.5000 0.5000 ... DONE
WANN -3.3363 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
PWSCF -3.3364 -3.3363 -0.7601 -0.7601 5.0942 5.0942 5.1269 5.1269
INFO: Performing a check on the eigenvalues of the rotated KS Hamiltonian ... DONE
INFO: Minimizing orbitals DEFINED
INFO: PREPARING THE KCW CALCULATION ...
INFO: Compute Wannier-orbital Densities ...
INFO: Coulomb q+G=0 treatment:
INFO: Divergence none
INFO: q-grid dimension 2 2 2
INFO: cell volume 270.256360799503
INFO: Gamma Extrapolation F
INFO: extrapolation q->0 term not estimated
INFO: Bare Coulomb q+G=0 0.00000E+00
==============================================================================
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.5000000 0.5000000 -0.5000000 [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.5000000 0.5000000 0.5000000 [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.0000000 1.0000000 0.0000000 [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.5000000 -0.5000000 0.5000000 [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 = -1.0000000 0.0000000 0.0000000 [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.0000000 1.0000000 [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.5000000 0.5000000 0.5000000 [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: Orbital Self-Hartree (SH)
orb 1 SH 0.272873
orb 2 SH 0.272874
orb 3 SH 0.272874
orb 4 SH 0.272873
orb 5 SH 0.272873
orb 6 SH 0.272874
orb 7 SH 0.272873
orb 8 SH 0.272874
iwann= 1 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 -0.000000
iwann= 1 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 1 ipol= 2 int rho_wann[ipol](r) [Re, Im] = -0.000001 0.000000
iwann= 1 ipol= 2 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 1 ipol= 3 int rho_wann[ipol](r) [Re, Im] = -0.000001 0.000000
iwann= 1 ipol= 3 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 1 ipol= 4 int rho_wann[ipol](r) [Re, Im] = 0.999997 -0.000000
iwann= 1 ipol= 4 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 2 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 -0.000000
iwann= 2 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 2 ipol= 2 int rho_wann[ipol](r) [Re, Im] = 0.000001 -0.000000
iwann= 2 ipol= 2 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 2 ipol= 3 int rho_wann[ipol](r) [Re, Im] = 0.000001 0.000000
iwann= 2 ipol= 3 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 2 ipol= 4 int rho_wann[ipol](r) [Re, Im] = -0.999997 0.000000
iwann= 2 ipol= 4 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 3 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 0.000000
iwann= 3 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 3 ipol= 2 int rho_wann[ipol](r) [Re, Im] = 0.000000 -0.000000
iwann= 3 ipol= 2 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 3 ipol= 3 int rho_wann[ipol](r) [Re, Im] = -0.000001 0.000000
iwann= 3 ipol= 3 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 3 ipol= 4 int rho_wann[ipol](r) [Re, Im] = 0.999997 0.000000
iwann= 3 ipol= 4 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 4 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 0.000000
iwann= 4 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 4 ipol= 2 int rho_wann[ipol](r) [Re, Im] = -0.000000 0.000000
iwann= 4 ipol= 2 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 4 ipol= 3 int rho_wann[ipol](r) [Re, Im] = 0.000001 0.000000
iwann= 4 ipol= 3 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 4 ipol= 4 int rho_wann[ipol](r) [Re, Im] = -0.999997 -0.000000
iwann= 4 ipol= 4 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 5 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 -0.000000
iwann= 5 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 5 ipol= 2 int rho_wann[ipol](r) [Re, Im] = -0.000000 -0.000000
iwann= 5 ipol= 2 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 5 ipol= 3 int rho_wann[ipol](r) [Re, Im] = 0.000001 0.000000
iwann= 5 ipol= 3 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 5 ipol= 4 int rho_wann[ipol](r) [Re, Im] = 0.999997 -0.000000
iwann= 5 ipol= 4 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 6 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 -0.000000
iwann= 6 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 6 ipol= 2 int rho_wann[ipol](r) [Re, Im] = 0.000000 0.000000
iwann= 6 ipol= 2 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 6 ipol= 3 int rho_wann[ipol](r) [Re, Im] = -0.000001 -0.000000
iwann= 6 ipol= 3 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 6 ipol= 4 int rho_wann[ipol](r) [Re, Im] = -0.999997 0.000000
iwann= 6 ipol= 4 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 7 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 -0.000000
iwann= 7 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 7 ipol= 2 int rho_wann[ipol](r) [Re, Im] = 0.000001 0.000000
iwann= 7 ipol= 2 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 7 ipol= 3 int rho_wann[ipol](r) [Re, Im] = 0.000000 -0.000000
iwann= 7 ipol= 3 int Im[rho_wann[ipol](r)] = -0.000000
iwann= 7 ipol= 4 int rho_wann[ipol](r) [Re, Im] = 0.999997 -0.000000
iwann= 7 ipol= 4 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 8 ipol= 1 int rho_wann[ipol](r) [Re, Im] = 1.000000 -0.000000
iwann= 8 ipol= 1 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 8 ipol= 2 int rho_wann[ipol](r) [Re, Im] = -0.000001 0.000000
iwann= 8 ipol= 2 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 8 ipol= 3 int rho_wann[ipol](r) [Re, Im] = -0.000000 0.000000
iwann= 8 ipol= 3 int Im[rho_wann[ipol](r)] = 0.000000
iwann= 8 ipol= 4 int rho_wann[ipol](r) [Re, Im] = -0.999997 0.000000
iwann= 8 ipol= 4 int Im[rho_wann[ipol](r)] = -0.000000
INFO: PREPARING THE KCW CALCULATION ... DONE
KCW : 0.86s CPU 1.01s WALL
INITIALIZATION:
map : 0.00s CPU 0.00s WALL ( 8 calls)
rho_of_q : 0.49s CPU 0.49s WALL ( 8 calls)
init_vloc : 0.00s CPU 0.00s WALL ( 1 calls)
init_us_1 : 0.01s CPU 0.02s WALL ( 1 calls)
h_psi : 0.04s CPU 0.04s WALL ( 8 calls)
h_psi : 0.04s CPU 0.04s WALL ( 8 calls)
add_vuspsi : 0.00s CPU 0.00s WALL ( 8 calls)
General routines
calbec : 0.00s CPU 0.00s WALL ( 8 calls)
fft : 0.01s CPU 0.04s WALL ( 7 calls)
ffts : 0.10s CPU 0.10s WALL ( 896 calls)
fftw : 0.24s CPU 0.24s WALL ( 2304 calls)
davcio : 0.01s CPU 0.01s WALL ( 176 calls)
KCW : 0.86s CPU 1.02s WALL
This run was terminated on: 15:28:27 16Feb2024
=------------------------------------------------------------------------------=
JOB DONE.
=------------------------------------------------------------------------------=
Reference in New Issue View Git Blame Copy Permalink