Search

Quantum Monte Carlo (QMC) methods are powerful approaches for solving electronic structure problems. Although they often provide high-accuracy solutions, the precision of most QMC methods is ultimately limited by a trial wave function that must be used. Recently, an approach has been demonstrated to allow the use of trial wave functions prepared on a quantum computer [Nature 603, 416 (2022)] in the auxiliary-field QMC (AFQMC) method, using classical shadows to estimate the required overlaps. However, this approach has an exponential post-processing step to construct these overlaps, when performing classical shadows obtained using random Clifford circuits. Here, we study an approach to avoid this exponentially scaling step by using a fixed-node Monte Carlo method, based on full configuration interaction quantum Monte Carlo (FCIQMC). This method is applied to the local unitary cluster Jastrow (LUCJ) ansatz. We consider H$_4$, ferrocene and benzene molecules using up to 12 qubits as examples. Circuits are compiled to native gates for typical near-term architectures, and we assess the impact of circuit-level depolarizing noise on the method. We also provide a comparison of AFQMC and fixed-node approaches, demonstrating that AFQMC is more robust to errors, although extrapolations of the fixed-node energy reduce this discrepancy. Although the method can be used to reach chemical accuracy, the sampling cost to achieve this is high even for small active spaces, suggesting caution for the prospect of outperforming conventional QMC approaches., Comment: 16 pages, 8 figures

Quantum error correction (QEC) will be essential to achieve the accuracy needed for quantum computers to realise their full potential. The field has seen promising progress with demonstrations of early QEC and real-time decoded experiments. As quantum computers advance towards demonstrating a universal fault-tolerant logical gate set, implementing scalable and low-latency real-time decoding will be crucial to prevent the backlog problem, avoiding an exponential slowdown and maintaining a fast logical clock rate. Here, we demonstrate low-latency feedback with a scalable FPGA decoder integrated into the control system of a superconducting quantum processor. We perform an 8-qubit stability experiment with up to $25$ decoding rounds and a mean decoding time per round below $1$ ${\mu}s$, showing that we avoid the backlog problem even on superconducting hardware with the strictest speed requirements. We observe logical error suppression as the number of decoding rounds is increased. We also implement and time a fast-feedback experiment demonstrating a decoding response time of $9.6$ ${\mu}s$ for a total of $9$ measurement rounds. The decoder throughput and latency developed in this work, combined with continued device improvements, unlock the next generation of experiments that go beyond purely keeping logical qubits alive and into demonstrating building blocks of fault-tolerant computation, such as lattice surgery and magic state teleportation., Comment: 11 pages, 4 figures, Supplementary Information

In this educational paper, we will discuss calculations on the hydrogen molecule both on classical and quantum computers. In the former case, we will discuss the calculation of molecular integrals that can then be used to calculate potential energy curves at the Hartree--Fock level and to correct them by obtaining the exact results for all states in the minimal basis. Some aspects of spin-symmetry will also be discussed. In the case of quantum computing, we will start out from the second-quantized Hamiltonian and qubit mappings. Using quantum phase estimation, we then provide the circuits for two different algorithms: Trotteization and qubitization. Finally, the significance of quantum error correction will be briefly discussed.

A number of exciting recent results have been seen in the field of quantum error correction. These include initial demonstrations of error correction on current quantum hardware, and resource estimates which improve understanding of the requirements to run large-scale quantum algorithms for real-world applications. In this work, we bridge the gap between these two developments by performing careful estimation of the resources required to fault-tolerantly perform quantum phase estimation (QPE) on a minimal chemical example. Specifically, we describe a detailed compilation of the QPE circuit to lattice surgery operations for the rotated surface code, for a hydrogen molecule in a minimal basis set. We describe a number of optimisations at both the algorithmic and error correction levels. We find that implementing even a simple chemistry circuit requires 1,000 qubits and 2,300 quantum error correction rounds, emphasising the need for improved error correction techniques specifically targeting the early fault-tolerant regime., Comment: 24 pages, 24 figures, 2 tables, source code available at https://github.com/riverlane/h2_compilation/

Quantum phase estimation (QPE) is a key quantum algorithm, which has been widely studied as a method to perform chemistry and solid-state calculations on future fault-tolerant quantum computers. Recently, several authors have proposed statistical alternatives to QPE that have benefits on early fault-tolerant devices, including shorter circuits and better suitability for error mitigation techniques. However, practical implementations of the algorithm on real quantum processors are lacking. In this paper we practically implement statistical phase estimation on Rigetti's superconducting processors. We specifically use the method of Lin and Tong [PRX Quantum 3, 010318 (2022)] using the improved Fourier approximation of Wan et al. [PRL 129, 030503 (2022)], and applying a variational compilation technique to reduce circuit depth. We then incorporate error mitigation strategies including zero-noise extrapolation and readout error mitigation with bit-flip averaging. We propose a simple method to estimate energies from the statistical phase estimation data, which is found to improve the accuracy in final energy estimates by one to two orders of magnitude with respect to prior theoretical bounds, reducing the cost to perform accurate phase estimation calculations. We apply these methods to chemistry problems for active spaces up to 4 electrons in 4 orbitals, including the application of a quantum embedding method, and use them to correctly estimate energies within chemical precision. Our work demonstrates that statistical phase estimation has a natural resilience to noise, particularly after mitigating coherent errors, and can achieve far higher accuracy than suggested by previous analysis, demonstrating its potential as a valuable quantum algorithm for early fault-tolerant devices., Comment: 24 pages, 13 figures

In this perspective, the various measures of electron correlation used in wavefunction theory, density functional theory and quantum information theory are briefly reviewed. We then focus on a more traditional metric based on dominant weights in the full configuration solution and discuss its behaviour with respect to the choice of the $N$-electron and the one-electron basis. The impact of symmetry is discussed and we emphasize that the distinction between determinants, configuration state functions and configurations as reference functions is useful because the latter incorporate spin-coupling into the reference and should thus reduce the complexity of the wavefunction expansion. The corresponding notions of single determinant, single spin-coupling and single configuration wavefunctions are discussed and the effect of orbital rotations on the multireference character is reviewed by analysing a simple model system. In molecular systems, the extent of correlation effects should be limited by finite system size and in most cases the appropriate choices of one-electron and $N$-electron bases should be able to incorporate these into a low-complexity reference function, often a single configurational one.

Computational chemistry is an essential tool in the pharmaceutical industry. Quantum computing is a fast evolving technology that promises to completely shift the computational capabilities in many areas of chemical research by bringing into reach currently impossible calculations. This perspective illustrates the near-future applicability of quantum computation to pharmaceutical problems. We briefly summarize and compare the scaling properties of state-of-the-art quantum algorithms, and provide novel estimates of the quantum computational cost of simulating progressively larger embedding regions of a pharmaceutically relevant covalent protein-drug complex involving the drug Ibrutinib. Carrying out these calculations requires an error-corrected quantum architecture, that we describe. Our estimates showcase that recent developments on quantum algorithms have dramatically reduced the quantum resources needed to run fully quantum calculations in active spaces of around 50 orbitals and electrons, from estimated over 1000 years using the Trotterisation approach to just a few days with sparse qubitisation, painting a picture of fast and exciting progress in this nascent field.

Quantum computers are special purpose machines that are expected to be particularly useful in simulating strongly correlated chemical systems. The quantum computer excels at treating a moderate number of orbitals within an active space in a fully quantum mechanical manner. We present a quantum phase estimation calculation on F$_2$ in a (2,2) active space on Rigetti's Aspen-11 QPU. While this is a promising start, it also underlines the need for carefully selecting the orbital spaces treated by the quantum computer. In this work, a scheme for selecting such an active space automatically is described and simulated results obtained using both the quantum phase estimation (QPE) and variational quantum eigensolver (VQE) algorithms are presented and combined with a subtractive method to enable accurate description of the environment. The active occupied space is selected from orbitals localized on the chemically relevant fragment of the molecule, while the corresponding virtual space is chosen based on the magnitude of interactions with the occupied space calculated from perturbation theory. This protocol is then applied to two chemical systems of pharmaceutical relevance: the enzyme [Fe] hydrogenase and the photosenzitizer temoporfin. While the sizes of the active spaces currently amenable to a quantum computational treatment are not enough to demonstrate quantum advantage, the procedure outlined here is applicable to any active space size, including those that are outside the reach of classical computation., Comment: Submitted to Journal of Computational Chemistry

We present a study of fixed and partial-node approximations in Slater determinant basis sets, using full configuration interaction quantum Monte Carlo (FCIQMC) to perform sampling. Walker annihilation in the FCIQMC method allows partial-node simulations to be performed, relaxing the nodal constraint to converge to the FCI solution. This is applied to ab initio molecular systems, using symmetry-projected Jastrow mean-field wave functions for complete active space (CAS) problems. Convergence and the sign problem within the partial-node approximation are studied, which is shown to eventually be limited in its use due to the large walker populations required. However the fixed-node approximation results in an accurate and practical method. We apply these approaches to various molecular systems and active spaces, including ferrocene and acenes. This also provides a test of symmetry-projected Jastrow mean-field wave functions in variational Monte Carlo (VMC) for a new set of problems. For trans-polyacetylene molecules and acenes we find that the time to perform a constant number of fixed-node FCIQMC iterations scales as O(N^{1.44}) and O(N^{1.75}) respectively, resulting in an efficient method for CAS-based problems that can be applied accurately to large active spaces., Comment: Supporting information is included

We present a stochastic approach to perform strongly contracted n-electron valence state perturbation theory (SC-NEVPT), which only requires one- and two-body reduced density matrices, without introducing approximations. We use this method to perform SC-NEVPT2 for CASSCF wave functions obtained from selected configuration interaction, although the approach is applicable to a larger class of wave functions, including those from variational Monte Carlo (VMC). The accuracy of this approach is demonstrated for small test systems, and the scaling is investigated with the number of virtual orbitals and the molecule size. We also find the SC-NEVPT2 energy to be relatively insensitive to the quality of the reference wave function. Finally, the method is applied to the Fe(II)-Porphyrin system with a (32e, 29o) active space, and to the isomerization of Cu2O2 in a (28e, 32o) active space., Comment: 13 pages, 3 figures

We present NECI, a state-of-the-art implementation of the Full Configuration Interaction Quantum Monte Carlo algorithm, a method based on a stochastic application of the Hamiltonian matrix on a sparse sampling of the wave function. The program utilizes a very powerful parallelization and scales efficiently to more than 24000 CPU cores. In this paper, we describe the core functionalities of NECI and recent developments. This includes the capabilities to calculate ground and excited state energies, properties via the one- and two-body reduced density matrices, as well as spectral and Green's functions for ab initio and model systems. A number of enhancements of the bare FCIQMC algorithm are available within NECI, allowing to use a partially deterministic formulation of the algorithm, working in a spin-adapted basis or supporting transcorrelated Hamiltonians. NECI supports the FCIDUMP file format for integrals, supplying a convenient interface to numerous quantum chemistry programs and it is licensed under GPL-3.0., Comment: 68 pages, 8 figures. To be published in the Journal of Chemical Physics, full supplementary files are to be published together with the article

PYSCF is a Python-based general-purpose electronic structure platform that both supports first-principles simulations of molecules and solids, as well as accelerates the development of new methodology and complex computational workflows. The present paper explains the design and philosophy behind PYSCF that enables it to meet these twin objectives. With several case studies, we show how users can easily implement their own methods using PYSCF as a development environment. We then summarize the capabilities of PYSCF for molecular and solid-state simulations. Finally, we describe the growing ecosystem of projects that use PYSCF across the domains of quantum chemistry, materials science, machine learning and quantum information science.

The computationally expensive evaluation and storage of high-rank reduced density matrices (RDMs) has been the bottleneck in the calculation of dynamic correlation for multireference wave functions in large active spaces. We present a stochastic formulation of multireference configuration interaction (MRCI) and perturbation theory (MRPT) that avoids the need for these expensive RDMs. The algorithm presented here is flexible enough to incorporate a wide variety of active space reference wave functions, including selected configuration interaction, matrix product states, and symmetry-projected Jastrow mean field wave functions. It enjoys the usual attractive features of Monte Carlo methods, such as embarrassing parallelizability and low memory costs. We find that the stochastic algorithm is already competitive with the deterministic algorithm for small active spaces, containing as few as 14 orbitals. We illustrate the utility of our stochastic formulation using benchmark applications., Comment: 10 pages, 4 figures

We present an approach to combining selected configuration interaction (SCI) and initiator full configuration interaction quantum Monte Carlo (i-FCIQMC). In the current i-FCIQMC scheme, the space of initiators is chosen dynamically by a population threshold. Here, we instead choose initiators as the selected space (V) from a prior SCI calculation, allowing substantially larger initiator spaces for a given walker population. While SCI+PT2 adds a perturbative correction in the first-order interacting space (FOIS) beyond V, the approach presented here allows a variational calculation in the same space, and a perturbative correction in the second-order interacting space. The use of a fixed initiator space reintroduces population plateaus into FCIQMC, but it is shown that the plateau height is typically only a small multiple of the size of V. Thus, for a comparable fundamental memory cost to SCI+PT2, a substantially larger space can be sampled. The resulting method can be seen as a complementary approach to SCI+PT2, which is more accurate but slower for a common selected/initiator space. More generally, our results show that approaches exist to significantly improve initiator energies in i-FCIQMC, while still ameliorating the fermion sign problem relative to the original FCIQMC method., Comment: 13 pages, 6 figures

We propose the use of preconditioning in FCIQMC which, in combination with perturbative estimators, greatly increases the efficiency of the algorithm. The use of preconditioning allows a time step close to unity to be used (without time-step errors), provided that multiple spawning attempts are made per walker. We show that this approach substantially reduces statistical noise on perturbative corrections to initiator error, which improve the accuracy of FCIQMC but which can suffer from significant noise in the original scheme. Therefore, the use of preconditioning and perturbatively-corrected estimators in combination leads to a significantly more efficient algorithm. In addition, a simpler approach to sampling variational and perturbative estimators in FCIQMC is presented, which also allows the variance of the energy to be calculated. These developments are investigated and applied to benzene (30e,108o), an example where accurate treatment is not possible with the original method., Comment: 15 pages, 7 figures

Building on the success of Quantum Monte Carlo techniques such as diffusion Monte Carlo, alternative stochastic approaches to solve electronic structure problems have emerged over the last decade. The full configuration interaction quantum Monte Carlo (FCIQMC) method allows one to systematically approach the exact solution of such problems, for cases where very high accuracy is desired. The introduction of FCIQMC has subsequently led to the development of coupled cluster Monte Carlo (CCMC) and density matrix quantum Monte Carlo (DMQMC), allowing stochastic sampling of the coupled cluster wave function and the exact thermal density matrix, respectively. In this article we describe the HANDE-QMC code, an open-source implementation of FCIQMC, CCMC and DMQMC, including initiator and semi-stochastic adaptations. We describe our code and demonstrate its use on three example systems; a molecule (nitric oxide), a model solid (the uniform electron gas), and a real solid (diamond). An illustrative tutorial is also included.

We perform excited-state variational Monte Carlo and diffusion Monte Carlo calculations using a simple and efficient wave function ansatz. This ansatz follows the recent variation-after-response formalism, accurately approximating a configuration interaction singles wave function as a sum of only two non-orthogonal Slater determinants, and further including important orbital relaxation. The ansatz is used to perform diffusion Monte Carlo calculations with large augmented basis sets, comparing to benchmarks from near-exact quantum chemical methods. The significance of orbital optimization in excited-state diffusion Monte Carlo is demonstrated, and the excited-state optimization procedure is discussed in detail. Diffuse excited states in water and formaldehyde are studied, in addition to the formaldimine and benzonitrile molecules., Comment: 12 pages, 3 figures

In this article we study examples of systematic biases that can occur in quantum Monte Carlo methods due to the accumulation of non-linear expectation values, and approaches by which these errors can be corrected. We begin with a study of the Krylov-projected FCIQMC (KP-FCIQMC) approach, which was recently introduced to allow efficient, stochastic calculation of dynamical properties. This requires the solution of a sampled effective Hamiltonian, resulting in a non-linear operation on these stochastic variables. We investigate the probability distribution of this eigenvalue problem to study both stochastic errors and systematic biases in the approach, and demonstrate that such errors can be significantly corrected by moving to a more appropriate basis. This is lastly expanded to include consideration of the correlation function QMC approach of Ceperley and Bernu, showing how such an approach can be taken in the FCIQMC framework., Comment: 12 pages, 7 figures

We formulate a general, arbitrary-order stochastic response formalism within the Full Configuration Interaction Quantum Monte Carlo framework. This modified stochastic dynamic allows for the exact response properties of correlated multireference electronic systems to be systematically converged upon for systems far out of reach of traditional exact treatments. This requires a simultaneous coupled evolution of a response state alongside the zeroth-order state, which is shown to be stable, non-transient and unbiased. We demonstrate this with application to the static dipole polarizability of molecular systems, and in doing so, resolve a discrepancy between restricted and unrestricted high-level coupled-cluster linear response results which were the high-accuracy benchmark in the literature.

We present a perturbative correction within initiator full configuration interaction quantum Monte Carlo (i-FCIQMC). In the existing i-FCIQMC algorithm, a significant number of spawned walkers are discarded due to the initiator criteria. Here we show that these discarded walkers have a form that allows calculation of a second-order Epstein-Nesbet correction, that may be accumulated in a trivial and inexpensive manner, yet substantially improves i-FCIQMC results. The correction is applied to the Hubbard model, the uniform electron gas and molecular systems., Comment: 5 pages, 1 figure, with Supplementary material also included

QMCPACK is an open source quantum Monte Carlo package for ab-initio electronic structure calculations. It supports calculations of metallic and insulating solids, molecules, atoms, and some model Hamiltonians. Implemented real space quantum Monte Carlo algorithms include variational, diffusion, and reptation Monte Carlo. QMCPACK uses Slater-Jastrow type trial wave functions in conjunction with a sophisticated optimizer capable of optimizing tens of thousands of parameters. The orbital space auxiliary field quantum Monte Carlo method is also implemented, enabling cross validation between different highly accurate methods. The code is specifically optimized for calculations with large numbers of electrons on the latest high performance computing architectures, including multicore central processing unit (CPU) and graphical processing unit (GPU) systems. We detail the program's capabilities, outline its structure, and give examples of its use in current research calculations. The package is available at http://www.qmcpack.org .

PySCF is a general-purpose electronic structure platform designed from the ground up to emphasize code simplicity, both to aid new method development, as well as for flexibility in computational workflow. The package provides a wide range of tools to support simulations of finite size systems, extended systems with periodic boundary conditions, low dimensional periodic systems, and custom Hamiltonians, using mean-field and post-mean-field methods with standard Gaussian basis functions. To ensure easy of extensibility, PySCF uses the Python language to implement almost all its features, while computationally critical paths are implemented with heavily optimized C routines. Using this combined Python/C implementation, the package is as efficient as the best existing C or Fortran based quantum chemistry programs. In this paper we document the capabilities and design philosophy of the current version of the PySCF package.