Do you want to publish a course? Click here

Non-orthogonal determinants in multi-Slater-Jastrow trial wave functions for fixed-node diffusion Monte Carlo

75   0   0.0 ( 0 )
 Added by Lucas Wagner
 Publication date 2018
  fields Physics
and research's language is English




Ask ChatGPT about the research

The accuracy and efficiency of ab-initio quantum Monte Carlo (QMC) algorithms benefits greatly from compact variational trial wave functions that accurately reproduce ground state properties of a system. We investigate the possibility of using multi-Slater-Jastrow trial wave functions with non-orthogonal determinants by optimizing identical single particle orbitals independently in separate determinants. As a test case, we compute variational and fixed-node diffusion Monte Carlo (FN-DMC) energies of a C$_2$ molecule. For a given multi-determinant expansion, we find that this non-orthogonal orbital optimization results in a consistent improvement in the variational energy and the FN-DMC energy on the order of a few tenths of an eV. Our calculations indicate that trial wave functions with non-orthogonal determinants can improve computed energies in a QMC calculation when compared to their orthogonal counterparts.



rate research

Read More

125 - Sam Azadi , , Thomas D. Kuhne 2016
We use the diffusion quantum Monte Carlo to revisit the enthalpy-pressure phase diagram of the various products from the different proposed decompositions of H$_2$S at pressures above 150~GPa. Our results entails a revision of the ground-state enthalpy-pressure phase diagram. Specifically, we find that the C2/c HS$_2$ structure is persistent up to 440~GPa before undergoing a phase transition into the C2/m phase. Contrary to density functional theory, our calculations suggest that the C2/m phase of HS is more stable than the I4$_1$/amd HS structure over the whole pressure range from 150 to 400 GPa. Moreover, we predict that the Im-3m phase is the most likely candidate for H$_3$S, which is consistent with recent experimental x-ray diffraction measurements.
124 - A. Badinski , R. J. Needs 2008
We report exact expressions for atomic forces in the diffusion Monte Carlo (DMC) method when using nonlocal pseudopotentials. We present approximate schemes for estimating these expressions in both mixed and pure DMC calculations, including the pseudopotential Pulay term which has not previously been calculated and the Pulay nodal term which has not been calculated for real systems in pure DMC simulations. Harmonic vibrational frequencies and equilibrium bond lengths are derived from the DMC forces and compared with those obtained from DMC potential energy curves. Results for four small molecules show that the equilibrium bond lengths obtained from our best force and energy calculations differ by less than 0.002 Angstrom.
By combining density-functional theory (DFT) and wave function theory (WFT) via the range separation (RS) of the interelectronic Coulomb operator, we obtain accurate fixed-node diffusion Monte Carlo (FN-DMC) energies with compact multi-determinant trial wave functions. In particular, we combine here short-range exchange-correlation functionals with a flavor of selected configuration interaction (SCI) known as emph{configuration interaction using a perturbative selection made iteratively} (CIPSI), a scheme that we label RS-DFT-CIPSI. One of the take-home messages of the present study is that RS-DFT-CIPSI trial wave functions yield lower fixed-node energies with more compact multi-determinant expansions than CIPSI, especially for small basis sets. Indeed, as the CIPSI component of RS-DFT-CIPSI is relieved from describing the short-range part of the correlation hole around the electron-electron coalescence points, the number of determinants in the trial wave function required to reach a given accuracy is significantly reduced as compared to a conventional CIPSI calculation. Importantly, by performing various numerical experiments, we evidence that the RS-DFT scheme essentially plays the role of a simple Jastrow factor by mimicking short-range correlation effects, hence avoiding the burden of performing a stochastic optimization. Considering the 55 atomization energies of the Gaussian-1 benchmark set of molecules, we show that using a fixed value of $mu=0.5$~bohr$^{-1}$ provides effective error cancellations as well as compact trial wave functions, making the present method a good candidate for the accurate description of large chemical systems.
Spin crossover molecules have recently emerged as a family of compounds potentially useful for implementing molecular spintronics devices. The calculations of the electronic properties of such molecules is a formidable theoretical challenge as one has to describe the spin ground state of a transition metal as the legand field changes. The problem is dominated by the interplay between strong electron correlation at the transition metal site and charge delocalization over the ligands, and thus it fits into a class of problems where density functional theory may be inadequate. Furthermore, the crossover activity is extremely sensitive to environmental conditions, which are difficult to fully characterize. Here we discuss the phase transition of a prototypical spin crossover molecule as obtained with diffusion Monte Carlo simulations. We demonstrate that the ground state changes depending on whether the molecule is in the gas or in the solid phase. As our calculation provides a solid benchmark for the theory we then assess the performances of density functional theory. We find that the low spin state is always over-stabilized, not only by the (semi-)local functionals, but even by the most commonly used hybrids (such as B3LYP and PBE0). We then propose that reliable results can be obtained by using hybrid functionals containing about 50% of exact-exchange.
Quantum Monte Carlo simulations of interacting electrons in solids often use Slater-Jastrow trial wave functions with Jastrow factors containing one- and two-body terms. In uniform systems the long-range behavior of the two-body term may be deduced from the random-phase approximation (RPA) of Bohm and Pines. Here we generalize the RPA to nonuniform systems. This gives the long-range behavior of the inhomogeneous two-body correlation term and provides an accurate analytic expression for the one-body term. It also explains why Slater-Jastrow trial wave functions incorporating determinants of Hartree-Fock or density-functional orbitals are close to optimal even in the presence of an RPA Jastrow factor. After adjusting the inhomogeneous RPA Jastrow factor to incorporate the known short-range behavior, we test it using variational Monte Carlo calculations. We find that the most important aspect of the two-body term is the short-range behavior due to electron-electron scattering, although the long-range behavior described by the RPA should become more important at high densities.
comments
Fetching comments Fetching comments
mircosoft-partner

هل ترغب بارسال اشعارات عن اخر التحديثات في شمرا-اكاديميا