Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userTotal forces in the diffusion Monte Carlo method with nonlocal pseudopotentials
113
0
0.0
(
0
)
Ask ChatGPT about the research
No Arabic abstract
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.
rate research
Read More
We present a version of the T-moves approach for treating nonlocal pseudopotentials in diffusion Monte Carlo which has much smaller time-step errors than the existing T-moves approaches, while at the same time preserving desirable features such as the upper-bound property for the energy. In addition, we modify the reweighting factor of the projector used in diffusion Monte Carlo to reduce the time-step error. The latter is applicable not only to pseudopotential calculations but to all-electron calculations as well.
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.
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.
We introduce a semistochastic implementation of the power method to compute, for very large matrices, the dominant eigenvalue and expectation values involving the corresponding eigenvector. The method is semistochastic in that the matrix multiplication is partially implemented numerically exactly and partially with respect to expectation values only. Compared to a fully stochastic method, the semistochastic approach significantly reduces the computational time required to obtain the eigenvalue to a specified statistical uncertainty. This is demonstrated by the application of the semistochastic quantum Monte Carlo method to systems with a sign problem: the fermion Hubbard model and the carbon dimer.
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.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
