Optimization problem of the AE CUTOFFS

General discussion of the Cambridge quantum Monte Carlo code CASINO; how to install and setup; how to use it; what it does; applications.

Optimization problem of the AE CUTOFFS

Postby Braulio_Brito » Tue Sep 12, 2017 2:09 pm

I am facing an optimization problem of the AE CUTOFFS (They are not being optmized) using the CASINO version 2.13.584. I have been previously working with the version 2.12.1 which works properly. Can anyone help me with this?
Braulio_Brito
 
Posts: 7
Joined: Tue Mar 04, 2014 11:07 am

Re: Optimization problem of the AE CUTOFFS

Postby Pablo_Lopez_Rios » Tue Sep 12, 2017 2:43 pm

Hi Braulio,

If I understand correctly you are trying to optimize the backflow AE cutoffs but their values do not change. Could you post your input / correlation.data / out files to see if there are any obvious problems?

Best,
Pablo
Hey there! I am using CASINO.
Pablo_Lopez_Rios
 
Posts: 39
Joined: Thu Jan 30, 2014 1:25 am

Re: Optimization problem of the AE CUTOFFS

Postby Braulio_Brito » Tue Sep 12, 2017 4:52 pm

Hi Pablo,

The files asked are attached.

Best regards,

Braulio Gabriel.
Attachments
ProblemCasino.tar.gz
(246.22 KiB) Downloaded 105 times
Braulio_Brito
 
Posts: 7
Joined: Tue Mar 04, 2014 11:07 am

Re: Optimization problem of the AE CUTOFFS

Postby Pablo_Lopez_Rios » Tue Sep 12, 2017 5:18 pm

Hi Braulio,

For some reason the AE cutoffs are being initialized to 1.1E-8 a.u., which is essentially impossible to optimize. Add the following to your initial correlation.data just before the "END BACKFLOW" label:

Code: Select all
 START AE CUTOFFS
 Nucleus ; Set ; Cutoff length     ;  Optimizable (0=NO; 1=YES)
  1         1         1.d0                         1
  2         1         1.d0                         1
  3         1         1.d0                         1
 END AE CUTOFFS


which will not only set the cutoffs to reasonable initial values but will also use the same L_AE for the three nuclei in your system.

I will look into why the code came up with such poor initial guesses for the cutoffs.

Hope this helps,
Pablo
Hey there! I am using CASINO.
Pablo_Lopez_Rios
 
Posts: 39
Joined: Thu Jan 30, 2014 1:25 am

Re: Optimization problem of the AE CUTOFFS

Postby Braulio_Brito » Wed Sep 13, 2017 4:47 pm

Thank you Pablo. I'll do some testing.
Braulio_Brito
 
Posts: 7
Joined: Tue Mar 04, 2014 11:07 am

Re: Optimization problem of the AE CUTOFFS

Postby Braulio_Brito » Wed Sep 13, 2017 4:57 pm

I found something different in the optimized coefficients of the AE CUTOFFS. They have the same values.

START AE CUTOFFS
Nucleus ; Set ; Cutoff length ; Optimizable (0=NO; 1=YES)
1 1 0.878846061910215 1
2 1 0.878846061910215 1
3 1 0.878846061910215 1
END AE CUTOFFS
Braulio_Brito
 
Posts: 7
Joined: Tue Mar 04, 2014 11:07 am

Re: Optimization problem of the AE CUTOFFS

Postby Pablo_Lopez_Rios » Wed Sep 13, 2017 5:08 pm

Hi Braulio,

Yes, I mentioned this in my previous post - the label "1" in the specification of the AE CUTOFFS forces them to be the same, so there are two fewer parameters to worry about. In your correlation.data file you force the parameters in the MU and PHI backflow terms to be the same for the three atoms in your system (by the specification of "Atoms in this set: 1 2 3"), so it makes sense for the AE CUTOFFS to be forced to be equal too.

In any case, if this is not what you want you can just set a different label for each cutoff:

Code: Select all
 START AE CUTOFFS
 Nucleus ; Set ; Cutoff length     ;  Optimizable (0=NO; 1=YES)
  1         1         1.d0                         1
  2         2         1.d0                         1
  3         3         1.d0                         1
 END AE CUTOFFS


Best,
Pablo
Hey there! I am using CASINO.
Pablo_Lopez_Rios
 
Posts: 39
Joined: Thu Jan 30, 2014 1:25 am

Re: Optimization problem of the AE CUTOFFS

Postby Braulio_Brito » Wed Sep 13, 2017 5:12 pm

Thanks Pablo for the explanation.
Braulio_Brito
 
Posts: 7
Joined: Tue Mar 04, 2014 11:07 am

Re: Optimization problem of the AE CUTOFFS

Postby Vladimir_Konjkov » Thu Oct 11, 2018 3:48 am

Hello all.

I noticed that when one optimize backflow parameters using VMC emin, the AE CUTOFFS parameters is bad or not optimized at all i.e. accepts random values or not optimal in terms of subsequent DMC calculation.

As stated in APPENDIX B: ZEROING THE BACKFLOW DISPLACEMENT AT AE ATOMS in article Inhomogeneous backflow transformations in quantum Monte Carlo P. Lopez Rıos, A. Ma, N. D. Drummond, M. D. Towler, and R. J. Needs
AE cutoff limit the radii of non-intersecting spheres around the atoms inside which an artificial multiplicative cutoff function is applied to all contributions to the backflow displacement of particle that do not depend on the distance to the central AE atom. This includes the homogeneous backflow displacement and the inhomogeneous contributions from neighboring atoms.


Manually selecting the values of the AE CUTOFFS parameters, I found that the optimal values fall into a range of 0.2-0.5 au.

With large values of the parameter the energy becomes worse in the subsequent DMC calculation, and with smaller values, the variance in VMC emin optimization increases and it is difficult to obtain reproducible values for optimal backflow parameters.

Although if one increase the number VMC configurations up to 10^6-10^7 with decreasing AE CUTOFFS to 0.01 and optimize subsequently one can get a better DMC energy value. This is true for both atoms and molecules and DMC energy get better roughly 1mHa per atom when AE cutoff decrease from 0.5 to 0.01 (for hydrogen atoms did not check).

I have a question how effective the formula (B1 from APPENDIX B) is for small AE CUTOFFS (order of 0.01 au) and whether it can be improved without complicating the calculations?

Best Vladimir.

PS.
I found http://www.lancaster.ac.uk/staff/drummo ... apers.html is very usefull for citing

PPS.
I use ETA term of order = 9
In Soviet Russia Casino play you.
Vladimir_Konjkov
 
Posts: 119
Joined: Wed Apr 15, 2015 3:14 pm


Return to The CASINO program

Who is online

Users browsing this forum: No registered users and 5 guests

cron