Noncovalent Interactions by QMC: Speedup by One-Particle Basis-Set Size Reduction

While it is empirically accepted that the fixed-node diffusion Monte-Carlo (FN-DMC) depends only weakly on the size of the one-particle basis sets used to expand its guiding functions, limits of this observation are not settled yet. Our recent work indicates that under the FN error cancellation conditions, augmented triple zeta basis sets are sufficient to achieve a benchmark level of 0.1 kcal/mol in a number of small noncovalent complexes. Here we report on a possibility of truncation of the one-particle basis sets used in FN-DMC guiding functions that has no visible effect on the accuracy of the production FN-DMC energy differences. The proposed scheme leads to no significant increase in the local energy variance, indicating that the total CPU cost of large-scale benchmark noncovalent interaction energy FN-DMC calculations may be reduced.Comment: ACS book chapter, accepte

Diffusion Monte Carlo Study of Para -Diiodobenzene Polymorphism Revisited

We revisit our investigation of the diffusion Monte Carlo (DMC) simulation of p-DIB molecular crystal polymorphism. [J. Phys. Chem. Lett. 2010, 1, 1789-1794] We perform, for the first time, a rigorous study of finite-size effects and choice of nodal surface on the prediction of polymorph stability in molecular crystals using fixed-node DMC. Our calculations are the largest which are currently feasible using the resources of the K computer and provide insights into the formidable challenge of predicting such properties from first principles. In particular, we show that finite-size effects can influence the trial nodal surface of a small (1×1×1) simulation cell considerably. We therefore repeated our DMC simulations with a 1×3×3 simulation cell, which is the largest such calculation to date. We used a DFT nodal surface generated with the PBE functional and we accumulated statistical samples with ∼6.4×105 core-hours for each polymorph. Our final results predict a polymorph stability consistent with experiment, but indicate that results in our previous paper were somewhat fortuitous. We analyze the finite-size errors using model periodic Coulomb (MPC) interactions and kinetic energy corrections, according to the CCMH scheme of Chiesa, Ceperley, Martin, and Holzmann. We investigate the dependence of the finite-size errors on different aspect ratios of the simulation cell (k-mesh convergence) in order to understand how to choose an appropriate ratio for the DMC calculations. Even in the most expensive simulations currently possible, we show that the finite size errors in the DMC total energies are far larger than the energy difference between the two polymorphs, although error cancellation means that the polymorph prediction is accurate. Finally, we found that the T-move scheme is essential for these massive DMC simulations in order to circumvent population explosions and large time-step biases.Chemistry and Chemical Biolog

Disentanglement of triplet and singlet states of azobenzene: direct EELS detection and QMC modeling

Singlet and triplet excited states of trans-azobenzene have been measured in the gas phase by electron energy loss spectroscopy (EELS). In order to interpret the strongly overlapping singlet and triplet bands in the spectra a set of large-scale correlated quantum Monte-Carlo (QMC) simulations was performed. The EELS/QMC combination of methods yields an excellent agreement between theory and experiment and for the two low-lying excited singlet and two low-lying triplet states permitted their unambiguous assignment. In addition, EELS revealed two overlapping electronic states in the band commonly assigned as S2, the lower one with a pronounced vibrational structure, the upper one structureless. Finally, the agreement between theory and experiment was shown to further increase by taking computationally into account the finite temperature effects

Practical Diffusion Monte Carlo Simulations for Large Noncovalent Systems

Fixed-node diffusion Monte Carlo (FNDMC) simulations are one of the most promising methods for describing the noncovalent systems to high accuracy within reasonable computational times. The advent of massively parallel computers enables one to apply FNDMC to various noncovalent systems such as supramolecules and molecular crystals. It is, however, to be noted that a reliable description of subtle noncovalent interactions requires a much higher accuracy than that of typical chemical bindings, e.g., the subchemical accuracy of 0.1 kcal/mol for small noncovalent complexes. This is a severe requirement for FNDMC based on stochastic approaches and raises the computational issues of reliable estimates of not only error bar, but also energy itself. Firstly, our recent works on several noncovalent systems are demonstrated. Then we address the issues and propose a new strategy for statistical estimates to meet the subchemical accuracy