Direct evaluation of the force constant matrix in quantum Monte Carlo

Simulating chemical reaction dynamics on quantum computer

The electronic energies of molecules have been successfully evaluated on quantum computers. However, more attention is paid to the dynamics simulation of molecules in practical applications. Based on the variational quantum eigensolver (VQE) algorithm, Fedorov et al. proposed a correlated sampling (CS) method and demonstrated the vibrational dynamics of H2 molecules [J. Chem. Phys. 154, 164103 (2021)]. In this study, we have developed a quantum approach by extending the CS method based on the VQE algorithm (labeled eCS-VQE) for simulating chemical reaction dynamics. First, the CS method is extended to the three-dimensional cases for calculation of first-order energy gradients, and then, it is further generalized to calculate the second-order gradients of energies. By calculating atomic forces and vibrational frequencies for H2, LiH, H+ + H2, and Cl- + CH3Cl systems, we have seen that the approach has achieved the CCSD level of accuracy. Thus, we have simulated dynamics processes for two typical chemical reactions, hydrogen exchange and chlorine substitution, and obtained high-precision reaction dynamics trajectories consistent with the classical methods. Our eCS-VQE approach, as measurement expectations and ground-state wave functions can be reused, is less demanding in quantum computing resources and is, therefore, a feasible means for the dynamics simulation of chemical reactions on the current noisy intermediate-scale quantum-era quantum devices.

Finite-Size Error Cancellation in Diffusion Monte Carlo Calculations of Surface Chemistry

The accurate prediction of reaction mechanisms in heterogeneous (surface) catalysis is one of the central challenges in computational chemistry. Quantum Monte Carlo methods─Diffusion Monte Carlo (DMC) in particular─are being recognized as higher-accuracy, albeit more computationally expensive, alternatives to Density Functional Theory (DFT) for energy predictions of catalytic systems. A major computational bottleneck in the broader adoption of DMC for catalysis is the need to perform finite-size extrapolations by simulating increasingly large periodic cells (supercells) to eliminate many-body finite-size effects and obtain energies in the thermodynamic limit. Here, we show that it is possible to significantly reduce this computational cost by leveraging the cancellation of many-body finite-size errors that accompanies the evaluation of energy differences when calculating quantities like adsorption (binding) energies and mapping potential energy surfaces. We analyze the cancellation and convergence of many-body finite-size errors in two well-known adsorbate/slab systems, H₂O/LiH(001) and CO/Pt(111). Based on this analysis, we identify strategies for obtaining binding energies in the thermodynamic limit that optimally utilize error cancellation to balance accuracy and computational efficiency. Using one such strategy, we then predict the correct order of adsorption site preference on CO/Pt(111), a challenging problem for a wide range of density functionals. Our accurate and inexpensive DMC calculations are found to unambiguously recover the top > bridge > hollow site order, in agreement with experimental observations. We proceed to use this DMC method to map the potential energy surface of CO hopping between Pt(111) adsorption sites. This reveals the existence of an L-shaped top-bridge-hollow diffusion trajectory characterized by energy barriers that provide an additional kinetic justification for experimental observations of CO/Pt(111) adsorption. Overall, this work demonstrates that it is routinely possible to achieve order-of-magnitude speedups and memory savings in DMC calculations by taking advantage of error cancellation in the calculation of energy differences that are ubiquitous in heterogeneous catalysis and surface chemistry more broadly.

Absence of diagonal force constants in cubic Coulomb crystals

The quasi-harmonic model proposes that a crystal can be modelled as atoms connected by springs. We demonstrate how this viewpoint can be misleading: a simple application of Gauss’s law shows that the ion–ion potential for a cubic Coulomb system can have no diagonal harmonic contribution and so cannot necessarily be modelled by springs. We investigate the repercussions of this observation by examining three illustrative regimes: the bare ionic, density tight-binding and density nearly-free electron models. For the bare ionic model, we demonstrate the zero elements in the force constants matrix and explain this phenomenon as a natural consequence of Poisson’s law. In the density tight-binding model, we confirm that the inclusion of localized electrons stabilizes all major crystal structures at harmonic order and we construct a phase diagram of preferred structures with respect to core and valence electron radii. In the density nearly-free electron model, we verify that the inclusion of delocalized electrons, in the form of a background jellium, is enough to counterbalance the diagonal force constants matrix from the ion–ion potential in all cases and we show that a first-order perturbation to the jellium does not have a destabilizing effect. We discuss our results in connection to Wigner crystals in condensed matter, Yukawa crystals in plasma physics, as well as the elemental solids.