I am a new user in QMC and I use this program for a little time, for this reason I don't have much experience.
This last week I tried to obtain a dissociation energy with the casino program. First, I optimized the Jastrow parameters and then I calculated the DMC energy. To obtain the energy I used 100000 dmc steps and 2000 configurations for the statistics and I repeated this for different time step(0.001, 0.0025 and 0.005). The problem is when the time step is small(0.001) I represent the reblock graph for the dmc and I don't find the "plateau" at small values of dtdmc in the reblock.
Block len Std error Err in error
1 1.190778E04 2.663150E07
2 1.579759E04 4.996579E07
4 2.119098E04 9.478781E07
8 2.847128E04 1.801075E06
16 3.808577E04 3.407375E06
32 5.042606E04 6.380601E06
64 6.531861E04 1.168938E05
128 8.417852E04 2.131036E05
256 1.056928E03 3.785678E05
512 1.269994E03 6.439546E05
1024 1.481055E03 1.064076E04
2048 1.718183E03 1.754235E04
4096 1.815489E03 2.638888E04 *** BEST ***
8192 1.489433E03 3.122125E04
16384 1.728491E03 5.359214E04
32768 1.592540E03 7.764297E04
Then I try to increment the number of steps (600000) but I found a value very similar.
When I increase the value of dtdmc(0.0025 and 0.005) a "plateau" appear at high values(Blocklen=4096) even if I increase the number of steps to 800000. For this reason I want to know which value of dtdmc I have to use, because I don't know how to choose the right value. Should I increase more the number of steps to obtain a better reblock or I use the lowest value of dtdmc in which appear the "plateau"?.
Question about dtdmc
Re: Question about dtdmc
Ivan,
Could you give more information about your specific system? Most importantly, are you using pseudo potentials or are you performing all electron calculations?
Thanks, Ryan.
Could you give more information about your specific system? Most importantly, are you using pseudo potentials or are you performing all electron calculations?
Thanks, Ryan.

 Posts: 3
 Joined: Thu Nov 17, 2016 5:55 pm
Re: Question about dtdmc
Dear Ryan,
Now I work with two benzenes and I performing a all electron calculation. Do you need more information?
Ivan
Now I work with two benzenes and I performing a all electron calculation. Do you need more information?
Ivan
Re: Question about dtdmc
Ivan,
I would say that your time steps seem reasonable then. Barring any malfunctions the only advice I can give is to accumulate more statistics. 800k steps for an AE dmc calculation is not massively big.
How many equilibration steps have you performed?
Are you running with an uptodate version of CASINO? (I. E. Is it an old version that was compiled with gfortran?)
Thanks, Ryan.
I would say that your time steps seem reasonable then. Barring any malfunctions the only advice I can give is to accumulate more statistics. 800k steps for an AE dmc calculation is not massively big.
How many equilibration steps have you performed?
Are you running with an uptodate version of CASINO? (I. E. Is it an old version that was compiled with gfortran?)
Thanks, Ryan.

 Posts: 3
 Joined: Thu Nov 17, 2016 5:55 pm
Re: Question about dtdmc
Dear Ryan,
I use the version 2.12 of the CASINO program.
I used 5000 equilibration steps, because I read that is enough with a small number of steps. Anyway I will try to run the calculation with more statistics steps and I will check the results, but the problem is that I don't know how to choose the value for dtdmc.
Ivan.
I use the version 2.12 of the CASINO program.
I used 5000 equilibration steps, because I read that is enough with a small number of steps. Anyway I will try to run the calculation with more statistics steps and I will check the results, but the problem is that I don't know how to choose the value for dtdmc.
Ivan.
Re: Question about dtdmc
Dear Ivan,
If you have compiled v2.12 of CASINO with gfortran, then be careful with your AE DMC calculations  we recently noticed a problem with the default gfortran compiler flags that leads to additional noise (and sometimes disastrous breaking of the variational principle) in AE DMC. This is fixed as of v2.13.599, but your version should still be bugged (check if the "Ofast" flag appears in <CASINO BASE DIRECTORY>/arch/data/f90/gnu.arch  if it does, I would replace it with "O3 ffastmath fstackarrays" and recompile (and rerun your DMC...).
If you use a different compiler, then this wont be an issue .
With regards to choosing a sensible dtdmc and a reasonable equilibration period, I would recommend the following advice.
1) dtdmc should be chosen such that the root meansquared distance diffused by an electron in one move (which is = sqrt(3 x dtdmc x p) with p = the step acceptance probability  usually modulated such that it is close to 1) is of order the smallest physical length scale in the problem (for AE calculations, this will be small, and hence dtdmc will typically be smaller than it would be if one used PPs for the same system).
2) the equilibration period should be chosen such that, on average, each electron in the system has chance to "properly explore the system"  if L is the longest length scale in your system, then you would want to have L ~ sqrt(3N_eq x dtdmc x p). Once you have decided on dtdmc via part (1), you can then find N_eq  hopefully a reasonable number of equilibration steps (something of order thousands is typical).
Hope this helps you!
Thanks, Ryan.
If you have compiled v2.12 of CASINO with gfortran, then be careful with your AE DMC calculations  we recently noticed a problem with the default gfortran compiler flags that leads to additional noise (and sometimes disastrous breaking of the variational principle) in AE DMC. This is fixed as of v2.13.599, but your version should still be bugged (check if the "Ofast" flag appears in <CASINO BASE DIRECTORY>/arch/data/f90/gnu.arch  if it does, I would replace it with "O3 ffastmath fstackarrays" and recompile (and rerun your DMC...).
If you use a different compiler, then this wont be an issue .
With regards to choosing a sensible dtdmc and a reasonable equilibration period, I would recommend the following advice.
1) dtdmc should be chosen such that the root meansquared distance diffused by an electron in one move (which is = sqrt(3 x dtdmc x p) with p = the step acceptance probability  usually modulated such that it is close to 1) is of order the smallest physical length scale in the problem (for AE calculations, this will be small, and hence dtdmc will typically be smaller than it would be if one used PPs for the same system).
2) the equilibration period should be chosen such that, on average, each electron in the system has chance to "properly explore the system"  if L is the longest length scale in your system, then you would want to have L ~ sqrt(3N_eq x dtdmc x p). Once you have decided on dtdmc via part (1), you can then find N_eq  hopefully a reasonable number of equilibration steps (something of order thousands is typical).
Hope this helps you!
Thanks, Ryan.