QMCPACK / QMCPACK/qmcpack

pw2qmcpack: wrong orbital norm due to use of QE bandgroups

Open
#3,095 6 comments 0 reactions 0 assignees View on GitHub
bug converter
Dominant language
C++
Stars
403
Forks
154
Avg merge
1d 13h
Merged PRs (30d)
90

Description

I am running into an issue with pw2qmcpack. I attempted to write orbitals from a hybrid calculation.
pw2qmcpack prints:

```
Parallelization info
--------------------
sticks: dense smooth PW G-vecs: dense smooth PW
Min 2529 2529 668 185725 185725 25283
Max 2530 2530 669 185726 185726 25285
Sum 40477 40477 10699 2971607 2971607 404541

Generating pointlists ...
new r_m : 0.2884 (alat units) 2.3621 (a.u.) for type 1
new r_m : 0.2824 (alat units) 2.3127 (a.u.) for type 2
new r_m : 0.2824 (alat units) 2.3127 (a.u.) for type 3

negative rho (up, down): 0.000E+00 5.934E-06
The wrong norm of k-point 15 band 1 , after collection before writing, is 10.000000000000030
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 5 band 1 , after collection before writing, is 10.000000000000050
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 7 band 1 , after collection before writing, is 10.000000000000089
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 11 band 1 , after collection before writing, is 10.000000000000091
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 21 band 1 , after collection before writing, is 9.999999999999972
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 9 band 1 , after collection before writing, is 10.000000000000002
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 19 band 1 , after collection before writing, is 10.000000000000012
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 13 band 1 , after collection before writing, is 9.999999999999947
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 3 band 1 , after collection before writing, is 9.999999999999899
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
The wrong norm of k-point 17 band 1 , after collection before writing, is 9.999999999999957
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
esh5 create pwscf_output/pwscf.pwscf.h5_part9
esh5 create pwscf_output/pwscf.pwscf.h5_part1
esh5 create pwscf_output/pwscf.pwscf.h5_part5
esh5 create pwscf_output/pwscf.pwscf.h5_part7
esh5 create pwscf_output/pwscf.pwscf.h5_part2
esh5 create pwscf_output/pwscf.pwscf.h5_part6
esh5 create pwscf_output/pwscf.pwscf.h5_part10
esh5 create pwscf_output/pwscf.pwscf.h5_part8
esh5 create pwscf_output/pwscf.pwscf.h5_part4
esh5 create pwscf_output/pwscf.pwscf.h5_part3
The wrong norm of k-point 1 band 1 , after collection before writing, is 10.000000000000075
The orbitals went wrong before being written to h5 file. Please first add debug=.true. in the pw2qmcpack input file to check if the orbitals can be read from QE files correctly.
esh5 create pwscf_output/pwscf.pwscf.h5
```

I found it interesting that the norm is integer valued, but not 1.
This is related to a previous issue, #2538 , however in my case the pseudos should be norm conserving -- and indeed, DFT+U orbitals can be generated using these same pseudos (though DFT+U conversions were run on a different machine but with same version of QE).

I have submitted a run with `debug=.true.` to the queue, but thought I would submit the issue in the interim.

Relevant files are here:
[p2q_run.zip](https://github.com/QMCPACK/qmcpack/files/6304017/p2q_run.zip)

**To Reproduce**
Steps to reproduce the behavior:
I am running QE 6.4.1 on NERSC Cori. Patched and built by hand.

**System:**
- NERSC Cori

Contributor guide

Open the contributing guide

Assessment

This issue has not been assessed yet.

Get new issues in your inbox

A short digest of beginner-friendly GitHub issues.