Jastrow function

Typically QMC calculations are carried out in the FNA using a Slater-Jastrow wave function which is written as a product of an antisymmetric function, and an exponentiated Jastrow function, F Wsj =. The antisymmetric function... 

Although the product form of a SJ wave function has typically been used with a single global Jastrow factor, recently [111] the benefit of separate Jastrow terms for each molecular orbital of the antisymmetric function was demonstrated. The latter approach improves treatment of local electron correlation by facilitating adjustment to the local molecular environment. Also, the nodal structure of the trial wave function can better reflect the parameters of multiple Jastrow functions than the single global Slater-Jastrow wave function. This modification notably improves the nodal structure of trial wave functions and FN-DMC energies. 

The Fade function has a cusp at r = 0 that can be adjusted to match the Coulomb cusp conditions by adjusting the a parameter. The Sun form also has a cusp, but approaches its asymptotic value far more quickly than the Fade function, which is useful for the linear scaling methods. An exponential form proposed by Manten and Luchow is similar to the Sun form, but shifted by a constant. By itself, the shift affects only the normalization of the Slater-Jastrow function, but has other consequences when the function is used to construct more elaborate correlation functions. The polynomial Fade function does not have a cusp, but its value goes to zero at a finite distance. 

This ansatz was first proposed by Boys and Handy (BH) [142] and Schmidt and Moskowitz (SM) later arrived at the same form by considering averaged backflow effects [143], The SMBH form includes both the e-n function (via the m, n, 0 = m, 0, 0 terms) and the homogeneous e-e Jastrow function (via the m, n, 0 = 0, 0, 0 terms). The remaining e-e-n terms modulate the e-e correlation function according to the e-n distances. 

The unreweighted variance does not minimize the variance of (A), but it is nevertheless a valid minimization criterion because the variance of the exact wave function will be zero regardless of the distribution of walker coordinates. The renewed interest of Drummond and Needs [154] stems from their observation that the unreweighted variance is sometimes especially easy to minimize. When the only parameters being optimized are linear parameters of the Jastrow function the local energy is... 

The Slater-Jastrow function is the standard form of the guide function (j) in QMC, a product of one or several Slater determinants and a Jastrow correlation function... 

Diedrich et al. could demonstrate with calculations on the dimers of methane, ammonia, and water, as well as the benzene dimer, that DMC performs very well on the whole range of interactions from pure dispersive to mainly electrostatic. " They used pseudopotentials and HF orbitals. With a similar approach, Korth et al. calculated the full S22 test set of dimers and the pairs of nucleic adds both in the Watson-Crick and the stacked conformation. The benchmark calculations revealed a mean absolute deviation for the binding energy of only 0.68kcal/mol. Very accurate results for the parallel displaced benzene dimer were obtained by Sorella et al. who obtained a binding energy of 2.2kcal/mol. These authors used their AGP approach with a Jastrow function and carefully optimized wave function parameters. 

The typical trial wavefunction for QMC calculations on molecular systems consists of the product of a Slater determinant multiplied by a second function, which accounts to some extent for electron correlation with use of interelectron distances. The trial wavefunctions are most often taken from relatively simple analytic variational calculations, in most cases from calculations at the SCF level. Thus, for the 10-electron system methane," the trial function may be the product of the SCF function, which is a 10 x 10 determinant made up of two 5x5 determinants, and a Jastrow function for each pair of electrons. 

The values of the Jastrow constants b and c may be specified as Vi for pairs of electrons with opposite spins and as % for pairs with identical spins. This avoids infinities in the local energy for two electrons at the same position. The Jastrow functions incorporate the main effects of electron-electron interactions and give a significant improvement over simple SCF trial functions. 

Of course the cusp can be represented by including the interparticle distances in a trial wave functions, most simply by means of Jastrow factors, exponentials of the inter-particle distances. But the problems of integral evaluation with such fac-... 

The function T is not provided by the QMC method. It need not be highly accurate in all regions in space, but it must have accurate nodal surfaces in 3n-dimensional space, and it must be accurate near singularities of V. In practice, simple SCF wave functions provide sufficiently accurate nodal surfaces, and symmetrical Jastrow factors [ll] that serve to satisfy... 

Now let us look at the paper. Eqn. (1) gives the form of the transcorrelated wave function C0, where C = li >jfiri,rj) is a Jastrow factor, and is a determinant. This compact wave function includes the effects of electron correlation through the introduction of r in C. The form C0 is taken as the trial wave function in quantum Monte Carlo (QMC) molecular computations today. Indeed the explicit form for/(r r,j is most often used by the QMC community. The transcorrelated wave function was obtained by solving (C //C - W) = 0, which Boys called the transcorrelation wave equation. Because C //C is a non-Hermitian operator, it was important to devise independent assessments of the accuracy of the wavefunction C0. 

Over the years there have been important progress in finding trial functions substantially more accurate then the pair product form for homogeneous systems [12,13]. Within the generalized Feynman-Kac formalism, it is possible to systematically improve a given trial function [13,14]. The first corrections to the pair product action with plane wave orbitals are a three-body correlation term which modifies the correlation part of the trial function (Jastrow) and a backfiow transformation which changes the orbitals and therefore the nodal structure (or the phase) of the trial function [14]. The new trial function has the form... 

For metallic hydrogen we have described a parameter-free trial function which does not need optimization. However, if we use the pair proton action both for molecular or LDA orbitals, we are left with free parameters in the Jastrow factor and with the width of the gaussians for molecular orbitals. Optimization of the parameters in a trial function is crucial for the success of VMC. Bad upper bounds do not give much physical information. Good trial functions will be needed in the Projector Monte Carlo method. First, we must decide on what to optimize and then how to perform the optimization. There are several possibilities for the quantity to optimize and depending on the physical system, one or other of the criteria may be best. 

Fig. 2. Ne = Np = 16, Vs = 1.31. Dependence of total energy, variance and energy difference for a pair of proton configurations S, S ) on the RQMC projection time. The study is performed for Te = 0.02Dotted lines represent the variational estimates with their error bars. In panel b) and c) the lines are exponential fits to data and in panel d) the continnons line is a linear fit in the region < 0.005. Black circles (3BF-A) are resnlts obtained with the analitical three-body and backflow trial wave functions discnssed earlier, the red triangle is a variational resnlt with a Slater-Jastrow trial function with simple plane wave orbitals and the blue squares are results from a trial function with LDA orbitals and an optmized two-body Jastrow...

If the r-space wave function is a linear combination of Slater determinants constructed from a set of spin-orbitals 0/, then its p-space counterpart is the same linear combination of Slater determinants constructed from the spin-momentals j obtained as Fourier transforms, Eq. (9), of the spin-orbitals. The overwhelming majority of contemporary r-space wave functions can be expressed as a linear combination of Slater determinants, and in these cases only three-dimensional Eourier transforms, Eq. (9), of the spin-orbitals are necessary to obtain the corresponding A -electron wave function in p space. Use of the full Eq. (5) becomes necessary only for wave functions, such as Hylleraas- or Jastrow-type wave functions, that are not built from a one-electron basis set. Examples of transformations to momentum space of such wave functions for He and H, can be found elsewhere [20-22]. 

Currently, the ubiquitous choice for the trial function is of the Slater-Jastrow or pair-product form. It is a linear combination of spin-up and spin-down determinants of one-body orbitals multiplied by a correlation factor represented by an exponential of one-body, two-body, and so on, terms [16, 30] ... 

A suitable trial wave function can be constructed from a Hartree-Fock wave function multiplied with a suitable correlation function, often taken as a Jastrow factor 7(r). 

One of the most recent theoretical advances has been to correlate almost all n-p and p-p data below 300 Mev in terms of boundary conditions on S, P and D states with only one energy dependent parameter (see Breit and Feshbach ). However, there is still no answer to whether it is possible to write an analytical expression as a function of energy for the law of force between two nucleons. It may be that the need of a Lorentz invariant expression for the interaction is the fundamental reason why potentials (even of the repulsive core type (Jastrow )] specifically have not been able to fit all of the data. 

The antisymmetric wave functions in the previous section account for electron correlation indirectly through correlation among the coefficients of the geminal or Cl expansions. More compact descriptions of electron correlation are achieved by Jastrow correlation functions that depend explicitly on interelectronic distance. A Jastrow correlation function F can be parameterized in an infinite number of ways. F can be partitioned into a hierarchy of terms. .. TV in which F describes correlations among n electrons. 

The QMC method is thought to be promising for the treatment of dynamical and static eiectron correiation effects with the compact functional form of wave functions. One standard form of the wave function is the Jastrow-Slater wave function. The Jastrow-Slater wave function is defined by... 

---