Toward Large-Scale AFQMC Calculations: Large Time Step Auxiliary-Field Quantum Monte Carlo.

We report modifications of the ph-AFQMC algorithm that allow the use of large time steps and reliable time step extrapolation. Our modified algorithm eliminates size-consistency errors present in the standard algorithm when large time steps are employed. We investigate various methods to approximate the exponential of the one-body operator within the AFQMC framework, distinctly demonstrating the superiority of Krylov methods over the conventional Taylor expansion. We assess various propagators within AFQMC and demonstrate that the Split-2 propagator is the optimal method, exhibiting the smallest time-step errors. For the HEAT set molecules, the time-step extrapolated energies deviate on average by only 0.19 kcal/mol from the accurate small time-step energies. For small water clusters, we obtain accurate complete basis-set binding energies using time-step extrapolation with a mean absolute error of 0.07 kcal/mol compared to CCSD(T). Using large time-step ph-AFQMC for the N2 dimer, we show that accurate bond lengths can be obtained while reducing CPU time by an order of magnitude.

Phaseless auxiliary field quantum Monte Carlo with projector-augmented wave method for solids.

We implement the phaseless auxiliary field quantum Monte Carlo method using the plane-wave based projector augmented wave method and explore the accuracy and the feasibility of applying our implementation to solids. We use a singular value decomposition to compress the two-body Hamiltonian and, thus, reduce the computational cost. Consistent correlation energies from the primitive-cell sampling and the corresponding supercell calculations numerically verify our implementation. We calculate the equation of state for diamond and the correlation energies for a range of prototypical solid materials. A down-sampling technique along with natural orbitals accelerates the convergence with respect to the number of orbitals and crystal momentum points. We illustrate the competitiveness of our implementation in accuracy and computational cost for dense crystal momentum point meshes compared to a well-established quantum-chemistry approach, the coupled-cluster ansatz including singles, doubles, and perturbative triple particle-hole excitation operators.

Benchmark Phaseless Auxiliary-Field Quantum Monte Carlo Method for Small Molecules.

We report a scalable Fortran implementation of the phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC) and demonstrate its excellent performance and beneficial scaling with respect to system size. Furthermore, we investigate modifications of the phaseless approximation that can help to reduce the overcorrelation problems common to the ph-AFQMC. We apply the method to the 26 molecules in the HEAT set, the benzene molecule, and water clusters. We observe a mean absolute deviation of the total energy of 1.15 kcal/mol for the molecules in the HEAT set, close to chemical accuracy. For the benzene molecule, the modified algorithm despite using a single-Slater-determinant trial wavefunction yields the same accuracy as the original phaseless scheme with 400 Slater determinants. Despite these improvements, we find systematic errors for the CN, CO2, and O2 molecules that need to be addressed with more accurate trial wavefunctions. For water clusters, we find that the ph-AFQMC yields excellent binding energies that differ from CCSD(T) by typically less than 0.5 kcal/mol.

Approaching the basis-set limit of the dRPA correlation energy with explicitly correlated and projector augmented-wave methods.

The direct random-phase approximation (dRPA) is used to calculate and compare atomization energies for the HEAT set and ten selected molecules of the G2-1 set using both plane waves and Gaussian-type orbitals. We describe detailed procedures to obtain highly accurate and well converged results for the projector augmented-wave method as implemented in the Vienna Ab initio Simulation Package as well as the explicitly correlated dRPA-F12 method as implemented in the TURBOMOLE package. The two approaches agree within chemical accuracy (1 kcal/mol) for the atomization energies of all considered molecules, both for the exact exchange as well as for the RPA. The root mean-square deviation is 0.41 kcal/mol for the exact exchange (evaluated using density functional theory orbitals) and 0.33 kcal/mol for exact exchange plus correlation from the RPA.