Calculating potential energy curves with fixed-node diffusion Monte Carlo: CO and N2

Diffusion quantum Monte Carlo method on diradicals using single- and multi-determinant-Jastrow trial wavefunctions and different orbitals

In this work, the diffusion quantum Monte Carlo (DMC) method is employed to calculate the energies of singlet and triplet states for a series of organic diradicals and diatomic diradicals with π2 configuration. Single-determinant-Jastrow (SDJ) trial wavefunctions for triplet states, two-determinant-Jastrow (2DJ) trial wavefunctions for the singlet states, and multi-determinant-Jastrow (MDJ) trial wavefunctions are employed in DMC calculations using restricted open-shell B3LYP (ROB3LYP) orbitals, complete-active-space self-consistent field (CASSCF) orbitals, state-average CASSCF orbitals, or frozen-CASSCF orbitals. Our results show that DMC energies using either SDJ/2DJ or MDJ with ROB3LYP orbitals are close to or lower than those with the other orbitals for organic diradicals, while they are not very sensitive to the employed orbitals for diatomic diradicals. Furthermore, using MDJ can reduce DMC energies to some extent for most of the investigated organic diradicals and some diatomic diradicals. The importance of MDJ on DMC energies can be estimated based on the percentage of main determinants in the CASCI wavefunction. On the other hand, singlet–triplet gaps can be calculated reasonably with DMC using MDJ with a mean absolute error of less than 2 kcal/mol with all these orbitals. CASCI wavefunctions using density functional theory orbitals are preferred in constructing MDJ trial wavefunctions in practical DMC calculations since it is easier to obtain such wavefunctions than CASSCF methods.

Accurate spectroscopic properties by diffusion quantum Monte Carlo calculations

The capability of Diffusion Quantum Monte Carlo (DMC) to produce high quality potential energy curve (PEC) was evaluated. H₂⁺, HeH⁺ and LiH PECs were built by all-electron fixed-node DMC calculations. Trial wave functions were obtained from Hartree-Fock (HF) (H₂⁺), MCSCF and CI (HeH⁺ and LiH) calculations multiplied by Jastrow factor. The quality of these generated PECs was analyzed throughout equilibrium distance, dissociation energy, vibrational energies and rovibrational spectroscopic constants (ω_e, ω_ex_e, ω_ey_e, α_e, γ_e and B_e). The Discrete Variable Representation (DVR) and the Dunham approaches were used to determine the rovibrational spectroscopic constants. The PECs and the aforementioned properties were also obtained by the following methods: MCSCF/aug-cc-pV5Z (LiH), CCSD(T)/aug-cc-pV5Z (HeH⁺ and LiH) and HF (H₂⁺ and HeH⁺) levels. The results of these DMC computations, specially the DMC-DVR procedure, are the most accurate among others DMC calculations available in the literature for these systems. They suggest that DMC can be used to achieve accurate PECs to produce spectroscopic properties with the same level of accuracy of theoretical benchmarks and experimental data of the literature.

Singlet-triplet gaps in diradicals obtained with diffusion quantum Monte Carlo using a Slater-Jastrow trial wavefunction with a minimum number of determinants

Diradicals are essential species in a wide range of chemical processes, whereas the computational study of their electronic structure often remains a challenge due to near-degeneracy of the frontier molecular orbitals. The fixed-node diffusion quantum Monte Carlo (FN-DMC) method is employed to calculate adiabatic energy gaps of some typical diradicals with the Slater-Jastrow trial wavefunction. The antisymmetrized part of the trial wavefunction is taken to be a linear combination of a minimum number of determinants using RB3LYP orbitals from the closed-shell singlet state or ROB3LYP orbitals from the triplet state. Our results show that using the two-determinant-Jastrow trial wavefunction is necessary to achieve reliable energy differences between closed-shell singlet states. The energy of the triplet state with MS = 1 is calculated to be lower than that with MS = 0 with FN-DMC even using trial wavefunctions with spin-pure states as their antisymmetrized parts and this difference is reduced with better orbitals. This indicates that the fixed-node error is smaller for the triplet state with MS = 1. Adiabatic energy gaps obtained from the present FN-DMC calculations are in reasonable agreement with available experimental values. Compared with results of the high level EOM-SF-CC method, energy gaps of FN-DMC with RB3LYP orbitals are slightly better than those using ROB3LYP orbitals and results of EOM-SF-CCSD. The present FN-DMC calculations using the simplest ansatz for the trial wavefunction can achieve reasonable results for these diradicals and they can readily be applied to large diradicals.

Bias cancellation in one-determinant fixed-node diffusion Monte Carlo: Insights from fermionic occupation numbers

The accuracy of the fixed-node diffusion Monte Carlo (FNDMC) depends on the node location of the supplied trial state Ψ_{T}. The practical FNDMC approaches available for large systems rely on compact yet effective Ψ_{T}, most often containing an explicitly correlated single Slater determinant (SD). However, SD nodes may be better suited to one system than to another, which may possibly lead to inaccurate FNDMC energy differences. It remains a challenge how to estimate nonequivalence or appropriateness of SDs. Here we use the differences of a measure based on the Euclidean distance between the natural orbital occupation number (NOON) vector of the SD and the exact solution in the NOON vector space, which can be viewed as a measure of SD nonequivalence and as a qualitative measure of the expected degree of nondynamic-correlation-related bias in FNDMC energy differences. This is explored on a set of small noncovalent complexes and covalent bond breaking of Si_{2} vs N_{2}. It turns out that NOON-based measures well reflect the magnitude and sign of the bias present in the data available, thus providing insights into the nature of bias cancellation in SD FNDMC energy differences.