VMC energy doesn't equal SCF one.

General discussion of the Cambridge quantum Monte Carlo code CASINO; how to install and setup; how to use it; what it does; applications.
Post Reply
Tomohiro Ichibha
Posts: 10
Joined: Sun May 24, 2015 8:14 am

VMC energy doesn't equal SCF one.

Post by Tomohiro Ichibha »

Dear all,

I performed VMC calculations with the CI trial wave function
,which is made without all pseudo potential and
of which basis set is gaussian basis basis set.
But I was not able to get correct vmc energy:
First, I performed the calculation with 'cusp correction', and
vmc energy was about 10000 Hartree higher than SCF energy.
Second, I performed the calculation without 'cusp correction', and
vmc energy was about 7 Hartree higher than SCF energy.

Could you tell me how to get correct vmc energy?

Ichibha
Vladimir_Konjkov
Posts: 165
Joined: Wed Apr 15, 2015 3:14 pm

Re: VMC energy doesn't equal SCF one.

Post by Vladimir_Konjkov »

Dear Tomohiro Ichibha.

It would be very useful to describe the steps to reproduce your results, including the version of software used, input files, and possibly additional instructions.
Otherwise it is impossible to understand what happened in your case.

Vladimir.
In Soviet Russia Casino plays you.
Neil Drummond
Posts: 113
Joined: Fri May 31, 2013 10:42 am
Location: Lancaster
Contact:

Re: VMC energy doesn't equal SCF one.

Post by Neil Drummond »

One obvious point: a difference of 10,000 Hartree (=0.27 MeV) indicates that something has gone completely wrong in the calculations!

Best wishes,

Neil.
Vladimir_Konjkov
Posts: 165
Joined: Wed Apr 15, 2015 3:14 pm

Re: VMC energy doesn't equal SCF one.

Post by Vladimir_Konjkov »

Neil Drummond wrote:One obvious point: a difference of 10,000 Hartree (=0.27 MeV) indicates that something has gone completely wrong in the calculations!

Best wishes,

Neil.
I managed several times to get such a big difference in energy by setting jastrow U TERM cutoff to >100 au in varmin_linjas optimisation, but sample variance was even greater in this case.
I do not know the details of the Ichibha calculation, so do not want anything to suppose.

Vladimir.
In Soviet Russia Casino plays you.
Tomohiro Ichibha
Posts: 10
Joined: Sun May 24, 2015 8:14 am

Re: VMC energy doesn't equal SCF one.

Post by Tomohiro Ichibha »

Dear all,

Sorry for my late response.
By substituting gpcc for normal cusp_correction, The calculation runs correctly.

Thank you very much, Vladimir and Neil.

Ichibha
Mike Towler
Posts: 239
Joined: Thu May 30, 2013 11:03 pm
Location: Florence
Contact:

Re: VMC energy doesn't equal SCF one.

Post by Mike Towler »

It would be interesting to see why that is. Can you send me the input files?

Cheers
Mike
Mdeible
Posts: 13
Joined: Tue Jun 11, 2013 2:46 pm

Re: VMC energy doesn't equal SCF one.

Post by Mdeible »

Hi,

Can you also state what molecule you are using? I have seen this problem before as well, where the gaussian_cusp_correction gives clearly erroneous values for the VMC calculation but the general purpose scheme gives reasonable values. I don't have specific data on hand, but I seem to recall it happening on molecules that contain several hydrogen atoms, such as water clusters or a benzene or anthracene molecule.

Thanks,
Mike
Mike Towler
Posts: 239
Joined: Thu May 30, 2013 11:03 pm
Location: Florence
Contact:

Re: VMC energy doesn't equal SCF one.

Post by Mike Towler »

Hi guys,

I'm now been given access to these files - I'll try to find time to investigate the alleged problem this week.

Cheers,
Mike
Post Reply