Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Is the ccECP Br.upf correct? #74

Open
zenandrea opened this issue May 11, 2022 · 1 comment
Open

Is the ccECP Br.upf correct? #74

zenandrea opened this issue May 11, 2022 · 1 comment

Comments

@zenandrea
Copy link

I noticed that the Br.upf file for the ccECP pseudo:
https://pseudopotentiallibrary.org/recipes/Br/ccECP/Br.upf
has one local component and two non local,
while the parametrised pseudo (e.g., https://pseudopotentiallibrary.org/recipes/Br/ccECP/Br.ccECP.gamess)
do have three non local components other than the local component.
Is it a mistake or is it intentional?

@aannabe
Copy link
Contributor

aannabe commented May 11, 2022

The reason for having different non-local components in the quantum chemistry formats vs *.upf files is that the former is a semi-local (SL) potential while the latter is a transformed Kleinman-Bylander (KB) potential. The transformation is carried out for efficiency purposes, it reduces the number of integrals from O(NM2) to O(NM) (M=number of basis functions, N=number of atoms). Therefore, these two potentials are actually slightly different except for the atomic configuration which was used to construct the KB potential in which case they produce numerically the same results.

In the KB transformation, one has the freedom to choose the local channel to be any l <= l_max + 1 number. So choosing l <= l_max will effectively reduce the number of non-local channels by 1. Of course, one can choose l_max + 1 to be local, then the number of non-local channels remains the same. The reason for not doing this is to achieve the best transferability (going from SL to KB) and to avoid the ghost states (i.e, in KB form, for a given l, eigenfunctions with more radial nodes can have lower energies than the nodeless ground state eigenfunctions, which is not physical). In addition, in this case, choosing local = l_max + 1 would require f-occupations which increases difficulties with DFT convergences for the reference configuration.

So actually this reduction is there for all semi-local ccECPs. So far, we didn't observe problems due this SL->KB transformation, especially given that in QMC we use the original semi-local potential (*.xml) with the correct number of non-local channels and without transformation biases.

This is my understanding, others can correct me if something looks off.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants