deepmodeling / abacus-develop

An electronic structure package based on either plane wave basis or numerical atomic orbitals.
http://abacus.ustc.edu.cn
GNU Lesser General Public License v3.0
169 stars 129 forks source link

I have the data-HR and data-SR, I need an rR file without scf calculation #4045

Open iduygnay opened 5 months ago

iduygnay commented 5 months ago

Background

rR is the output when I turn on "out_mat_r" tag. b6f443f71a7d3f630f2a2357042d787

Describe the solution you'd like

Is it possible to get the rR matrix without scf calculation?

Task list only for developers

Notice Possible Changes of Behavior (Reminder only for developers)

No response

Notice any changes of core modules (Reminder only for developers)

No response

Notice Possible Changes of Core Modules (Reminder only for developers)

No response

Additional Context

No response

Task list for Issue attackers (only for developers)

iduygnay commented 5 months ago

In detail, now I have the data-HR and data-SR, I need an rR file without scf calculation, and the number of R should be the same as the HR and SR I provided. data-SR-sparse_SPIN0.zip

iduygnay commented 5 months ago

data-HR-sparse_SPIN0.zip

dyzheng commented 5 months ago

You can wait for the merging of PR#4084 and try to set "calculation nscf" with "out_mat_r".

Or you can set "calculation scf" with "scf_nmax 1" and "out_mat_r" in current version.

iduygnay commented 5 months ago

rR.zip This is the rR I got with the calculation you mentioned above. But pyatb reported that the shape is different from the SR or HR. I offered the SR and HR files above. Can you check the shape of them? I think the main issue is the number of R. image If you check the shapes for them are the same, I will contact with pyatb developers.

iduygnay commented 5 months ago

This could be solved by pyatb developers. Thanks!

iduygnay commented 4 months ago

I found that for the large system, it always goes to oom for scf calculation or I have to wait a long time for wfc initialization. But the only thing I want to get is the rR matrix. Is there any suggestion to get this matrix without a large memory? Or can it be quick?

I also tried nscf calculation, but I don't have the charge or wfc files, tho it can not give me the rR matrix.

mohanchen commented 11 hours ago

We will try to have a better solution to this feature.