DFT is a valuable tool for calculating adsorption energies toward designing materials for hydrogen storage. However, dispersion forces being absent from the theory, it remains unclear how the consideration of van der Waals (vdW) interactions affects
such calculations. For the first time, we applied diffusion Monte Carlo (DMC) to evaluate the adsorption characteristics of a hydrogen molecule on a (5,5) armchair silicon-carbide nanotube (H$_2$-SiCNT). Within the framework of density functional theory (DFT), we also benchmarked various exchange-correlation functionals, including those recently developed for treating dispersion or vdW interactions. We found that the vdW-corrected DFT methods agree well with DMC, whereas the local (semilocal) functional significantly over (under)-binds. Furthermore, we fully optimized the H$_2$-SiCNT geometry within the DFT framework and investigated the correlation between structure and charge density. The vdW contribution to adsorption was found to be non-negligible at approximately 1 kcal/mol per hydrogen molecule, which amounts to 9-29 % of the ideal adsorption energy required for hydrogen storage applications.
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 tr
ial 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.
While Diffusion Monte Carlo (DMC) is in principle an exact stochastic method for textit{ab initio} electronic structure calculations, in practice the fermionic sign problem necessitates the use of the fixed-node approximation and trial wavefunctions
with approximate nodes (or zeros) must be used. This approximation introduces a variational error in the energy that potentially can be tested and systematically improved. Here, we present a computational method that produces trial wavefunctions with systematically improvable nodes for DMC calculations of periodic solids. These trial wavefunctions are efficiently generated with the configuration interaction using a perturbative selection made iteratively (CIPSI) method. A simple protocol in which both exact and approximate results for finite supercells are used to extrapolate to the thermodynamic limit is introduced.
Quantum Monte Carlo (QMC) methods are some of the most accurate methods for simulating correlated electronic systems. We investigate the compatibility, strengths and weaknesses of two such methods, namely, diffusion Monte Carlo (DMC) and auxiliary-fi
eld quantum Monte Carlo (AFQMC). The multi-determinant trial wave functions employed in both approaches are generated using the configuration interaction using a perturbative selection made iteratively (CIPSI) technique. Complete basis set full configuration interaction (CBS-FCI) energies estimated with CIPSI are used as a reference in this comparative study between DMC and AFQMC. By focusing on a set of canonical finite size solid state systems, we show that both QMC methods can be made to systematically converge towards the same energy once basis set effects and systematic biases have been removed. AFQMC shows a much smaller dependence on the trial wavefunction than DMC while simultaneously exhibiting a much larger basis set dependence. We outline some of the remaining challenges and opportunities for improving these approaches.
We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been
greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing accuracy. Advances in real space methods include techniques for accurate computation of band gaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods including GW and density functional based techniques. To provide an improved foundation for these calculations we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.
Diffusion Monte Carlo (DMC) calculations were performed for an accurate description of the nature of the O$_2$ adsorption on a single layer graphene. We investigated the stable orientation of O$_2$ at a specific adsorption site as well as its equilib
rium adsorption energy. At equilibrium adsorption distances, an O$_2$ molecule was found to prefer a horizontal orientation, where the O-O bond is parallel to the graphene surface, to the vertical orientation. However, the vertical orientation is favored at the O$_2$-graphene distances shorter than the equilibrium distance, which could be understood by the steric repulsion between O and C atoms. Contrary to previous DFT calculations, our DMC calculations show that the midpoint of a C-C bond (a bridge site) is energetically preferred for the O$_2$ adsorption to a center of a hexagonal ring (a hollow site). The lowest DMC adsorption energy was found at an intermediate point between a hollow and a bridge site, where the O$_2$ adsorption energy was estimated to be -0.142(4) eV that was in very good agreement with the recently-reported experimental value. Finally, we have found that O$_2$ is very diffusive on the surface of graphene with the diffusion barrier along a bridge-hollow-bridge path being as small as ~ 11 meV.
Cubic hafnia (HfO$_2$) is of great interest for a number of applications in electronics because of its high dielectric constant. However, common defects in such applications degrade the properties of hafina. We have investigated the electronic proper
ties of oxygen vacancies and nitrogen substitution in cubic HfO$_2$ using first principles calculations based on density functional theory (DFT) and many-body diffusion Monte Carlo (DMC) methods. We investigate five different charge defect states of oxygen vacancies, as well as substitutional N defects which can lead to local magnetic moments. Both DMC and DFT calculations shows that an oxygen vacancy induces strong lattice relaxations around the defect. Finally, we compare defect formation energies, charge and spin densities obtained from DMC with results obtained using DFT. While the obtained formation energies from DMC are 0.6~eV -- 1.5~eV larger than those from GGA+U, the agreement for the most important defects, neutral and positively charged oxygen vacancies, and nitrogen substitutional defect, under oxygen-poor conditions are in reasonably good agreement. Our work confirms that nitrogen can act to passivate cubic hafnia for applications in electronics.
Zirconia (zirconium dioxide) and hafnia (hafnium dioxide) are binary oxides used in a range of applications. Because zirconium and hafnium are chemically equivalent, they have three similar polymorphs, and it is important to understand the properties
and energetics of these polymorphs. However, while density functional theory calculations can get the correct energetic ordering, the energy differences between polymorphs depend very much on the specific density functional theory approach, as do other quantities such as lattice constants and bulk modulus. We have used highly accurate quantum Monte Carlo simulations to model the three zirconia and hafnia polymorphs. We compare our results for structural parameters, bulk modulus, and cohesive energy with results obtained from density functional theory calculations. We also discuss comparisons of our results with existing experimental data, in particular for structural parameters where extrapolation to zero temperature can be attempted. We hope our results of structural parameters as well as for cohesive energy and bulk modulus can serve as benchmarks for density-functional theory based calculations and as a guidance for future experiments.
