Helium computation behaves strangely

Asked by Tommi Höynälänmaa

Helium computation does not get more accurate even though I double the grid point density. Do you have any idea why this is so?

The computations can be found at
http://www.students.tut.fi/~hoynalan/He-20.tar.bz2
and
http://www.students.tut.fi/~hoynalan/He-21.tar.bz2

Yours sincerely,
     Tommi Höynälänmaa

Question information

Language:
English Edit question
Status:
Solved
For:
BigDFT Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Luigi Genovese (luigi-genovese) said :
#1

Dear Tommi,

thanks for your message.
The principal problem in your calculation is the value of crmult: 1.0 is wayy to little for any production result.

Please have a look at the tutorial page on how the basis set should be varied:

http://bigdft.org/Wiki/index.php?title=Basis-set_convergence

Seconldy, for a preiodic calculation as you want to do, the best is always to have a crmult such that the whole box is filled with grid points.
For example, the He20 calculation gives the warning:

- The coarse grid does not fill the entire periodic box

(see at the end of the logfile the list of the warnings.

This means that you are placing the nonzero degrees of freedom only close to the atoms. THis will confine too much the wavefunction and breaks translational invariance, that must be preserved for a 3d periodic calculation.

Do not hesitate in coming back to us should you encounter any problem.

Best

Luigi

Revision history for this message
Luigi Genovese (luigi-genovese) said :
#2

Same reason as Question 254943.