The CASINO forum

Also, you can have a fixed node error. If you are not working with the full wave function that gave the reference result in the first place (and I think that the value you cite was obtained by CCSDT), DMC may well also give an energy that is higher than the value from the "standard" quant...

Vladimir, I'm not sure what you are asking here. You're ground state energy is: -109.542(1) a.u., whilst the exact value (from Ref. 38 of the quoted PRE article) is -109.5421 a.u. These two are consistent with each other (at the level of statistical significance you are working), and your VMC energ...

Hello CASINO developers I'm using the QChem program to get trial WFN Qchem has “valence” active space methods VOD which allows you to get WFN and T2-amplitudes to construct CASINO input https://www.q-chem.com/qchem-website/manual/qchem43_manual/sect-ccactive.html When I create VOD single determinant...

Hello Cyrus I just read the article and want to note that in the case of single determinant WFN for Be atom and Be2 molecule which are known to give the wrong topology of nodal surface DMC time step bias no longer linear I applied polynomial interpolation be_be2.png E(Be) = -0.022459*x^1.5 - 0.01044...

If you look at Fig. 13 of Umrigar, Nightingale, Runge, J. Chem. Phys. 1993, you will see that it is nearly quadratic near tau=0. I do not know of a theoretical argument for the behaviour of the time-step error. What is the x-scale on your plots? Cyrus X scale is time step divided by τ = 1/(3Z^2), w...

Hello everyone who does DMC calculations. I calculated the time step bias dependence for simple full electron systems like H-Cl atoms and some diatomic. I've found that small time step asymptotic is rather quadratic than linear whereas large time step asymptotic is linear in some cases. For example ...

Hello CASINO developers. I wonder why the CASINO uses a spherical representation for the gaussian basic functions, while the сartesian representation requires a smaller number of computations to find both the values, derivatives and Laplacian etc. There is no difference only for S and P the orbitals...

Vladimir_Konjkov wrote: Yes, I also wanted to improve significantly your inelegant code, and compiled profiling version that gave me 'empty' gmon.out.

FFLAGS_prof = -pg -no-pie fix generation of the gmon.out
Also mpif90 is the deprecated OpenMPI wrapper in favor of mpifort.

Hello Mike. In the gm2_bf SUBROUTINE evaluates the elements of the prim array in CASE of g-orbitals, it seems to me that in this function such calculations are not required. I apologize, I just automatically made the code review. Some guy tried to do that once to the Gaussian evaluators (he was talk...

Hello all It is well known that for effective DMC calculation one need, among other things, to choose the optimal time-step. Fort all-electron atomic systems, the linear regime to occur for time steps less than of the order τ= 1/(3Z^2), where Z is the largest atomic number occurring in the system. h...