# Arbitrary Ground State Observables from Quantum Computed Moments
**Authors**: Harish J. Vallury1 and Lloyd C. L. Hollenberg
> 0009-0007-3355-54580000-0001-7672-6965School of Physics,
The University of MelbourneParkville, VIC 3010, AustraliaEmail: 1hvallury@student.unimelb.edu.au,
2lloydch@unimelb.edu.au
Abstract
The determination of ground state properties of quantum systems is a fundamental problem in physics and chemistry, and is considered a key application of quantum computers. A common approach is to prepare a trial ground state on the quantum computer and measure observables such as energy, but this is often limited by hardware constraints that prevent an accurate description of the target ground state. The quantum computed moments (QCM) method has proven to be remarkably useful in estimating the ground state energy of a system by computing Hamiltonian moments with respect to a suboptimal or noisy trial state. In this paper, we extend the QCM method to estimate arbitrary ground state observables of quantum systems. We present preliminary results of using QCM to determine the ground state magnetisation and spin-spin correlations of the Heisenberg model in its various forms. Our findings validate the well-established advantage of QCM over existing methods in handling suboptimal trial states and noise, extend its applicability to the estimation of more general ground state properties, and demonstrate its practical potential for solving a wide range of problems on near-term quantum hardware.
I Introduction
Many difficult problems with practical application in physics and chemistry involve determining the ground state properties of a quantum system. This has been a popular road to follow towards the development of near-term applications of quantum computers, with much of the recent focus primarily on variational quantum algorithms [1, 2, 3]. These hybrid approaches involve the preparation of complicated states designed to represent the ground state of the desired quantum system, via the measurement and minimisation of energy [4, 5, 6]. Most interesting problems at scale require high-depth circuits to prepare such states, and are thus subject to noise which is often severely detrimental to the quality of results [7, 8].
Some time ago, an analytic cluster expansion of the Lanczos recursion was uncovered [9]. Following from this, it was found that for a given quantum system Hamiltonian $H$ , one could compute Hamiltonian moments $\langle H^{k}\rangle$ with respect to a trial state with some overlap with the ground state to produce a powerful approximation to the ground state energy [10, 11], shown below to fourth order in the moments:
$$
E_{\rm 0}^{\rm{L}(4)}=c_{1}-\frac{c_{2}^{2}}{c_{3}^{2}-c_{2}c_{4}}\left[\sqrt{%
3c_{3}^{2}-2c_{2}c_{4}}-c_{3}\right], \tag{4}
$$
where the cumulants $c_{n}$ are defined as:
$$
c_{n}=\langle H^{n}\rangle-\sum_{p=0}^{n-2}\binom{n-1}{p}c_{p+1}\langle H^{n-1%
-p}\rangle. \tag{2}
$$
The rapid development of quantum computer technology has provided a fitting application for this result, allowing for advanced state preparation and measurement of these Hamiltonian moments in order to solve interesting problems. We refer to this as the quantum computed moments (QCM) method, and recent work applied to problems in quantum magnetism [12, 13] and chemistry [14] has demonstrated the effectiveness of QCM in estimating the ground state energy of a system, even when moments are measured with respect to a suboptimal or noisy trial state. This technique can correct for a low complexity ansatz and a high degree of noise, proving to be remarkably useful on present-day quantum hardware for a variety of different problem instances.
In recent times, these QCM studies have been accompanied elsewhere in the literature by a general renewed interest in Hamiltonian moments in a quantum computing context [15, 16, 17, 18], with the primary focus being on the ground state energy problem. However, a problem which is much more practically useful and indeed more computationally complex is that of determining arbitrary ground state observables [19]. It turns out that Hamiltonian moments can also be used to estimate general ground state properties via a Hellman-Feynman approach [20, 21]. Here we develop this technique in the QCM context. When applied to models of interest, our results show that the crucial robustness property of QCM for energy problems indeed carries over to this more general ground state observable procedure.
II Method
Suppose we wish to find the expectation value of some general observable (given by Hermitian operator $A$ ) with respect to the ground state $\ket{\psi_{0}}$ of a Hamiltonian $H$ . Consider the following parameterised Hamiltonian:
$$
\displaystyle H_{\lambda}=H+\lambda A \tag{3}
$$
The parameterised Hamiltonian $H_{\lambda}$ has ground state $\ket{\psi_{\lambda,0}}$ with energy $E_{\lambda,0}$ which we can approximate with:
$$
E_{\lambda,0}\approx E^{{\rm L}(4)}_{\lambda,0}|_{\psi_{t}}. \tag{4}
$$
The approximation in Equation (4) is expressed entirely in terms of Hamiltonian moments $\langle\psi_{t}|(H_{\lambda})^{k}|\psi_{t}\rangle$ , up to $k=4$ , with respect to some trial state $\ket{\psi_{t}}$ which has nonzero overlap with $\ket{\psi_{\lambda,0}}$ .
Starting with the Hellmann-Feynman theorem,
$$
\langle\psi_{\lambda,0}|\,\diff{H_{\lambda}}{\lambda}\,|\psi_{\lambda,0}%
\rangle=\diff{E_{\lambda,0}}{\lambda}, \tag{5}
$$
one can substitute the derivative of the Hamiltonian in Equation (3) and use the approximation in Equation (4) to obtain:
$$
\langle\psi_{\lambda,0}|\,A\,|\psi_{\lambda,0}\rangle\approx\frac{\rm d}{%
\mathrm{d}\lambda}\left(E^{{\rm L}(4)}_{\lambda,0}|_{\psi_{t}}\right). \tag{4}
$$
Thus if we want to compute the observable $A$ with respect to the ground state $\ket{\psi_{0}}$ of $H$ , we can choose $\ket{\psi_{t}}$ to have some overlap with $\ket{\psi_{0}}$ and evaluate the derivative in Equation (6) around $\lambda=0$ :
$$
\langle\psi_{0}|\,A\,|\psi_{0}\rangle\approx\frac{\rm d}{\mathrm{d}\lambda}%
\left(E^{{\rm L}(4)}_{\lambda,0}|_{\psi_{t}}\right)_{\lambda=0}. \tag{4}
$$
This leads to a practically useful final expression for the ground state observable:
$$
\langle A\rangle_{0}^{\rm L(4)}\equiv\frac{1}{2\varepsilon}\left(E_{+%
\varepsilon,0}^{\rm L(4)}-{E_{-\varepsilon,0}^{\rm L(4)}}\right)\approx\langle%
\psi_{0}|\,A\,|\psi_{0}\rangle, \tag{4}
$$
where $\varepsilon>0$ is sufficiently small as to find the $\lambda=0$ derivative, and the $E_{±\varepsilon,0}^{\rm L(4)}$ are calculated via Equation (1) from the moments $\langle\,(H±\varepsilon A)^{k}\,\rangle$ with respect to some trial state $\ket{\psi_{t}}$ up to $k=4$ .
We expand upon the analytic studies of this method, first introduced in [20, 21], by now considering the approach in a quantum computing context. Note that, in practice, the $E_{+\varepsilon,0}^{\rm L(4)}$ and ${E_{-\varepsilon,0}^{\rm L(4)}}$ in Equation 8 are not computed separately, but rather come from the same set of measurements on the quantum computer. The derivative is taken with respect to a parameter only introduced in the classical post-processing stage, so remains unaffected by the stochastic nature of the quantum measurement. We can thus expect to see similar robustness behaviour as has been observed in the QCM energy estimate, and in the following section, we test the effectiveness of this arbitrary observables technique by performing noisy simulations on a variety of trial states and Hamiltonians.
III Results
<details>
<summary>x1.png Details</summary>

### Visual Description
## Combined Chart: Energy and Correlation vs. Position
### Overview
The image presents two line charts, (a) and (b), displaying the relationship between position (x) and energy (H, E0^(L(4)), E0) in chart (a), and position (x) and correlation (C, C0^(L(4)), C0) in chart (b). Both charts share the same x-axis and are divided into three regions, corresponding to different states denoted as |ψ0^(x=-1)>, |ψ0^(x=0)>, and |ψ0^(x=1)>. The charts also include diagrams of a 3x3 grid, representing a lattice structure, with specific elements highlighted in chart (b).
### Components/Axes
**Chart (a): Energy vs. Position**
* **Title:** (a)
* **Y-axis:** Energy, with a scale from -0.55 to -0.35 in increments of 0.05.
* **X-axis:** x, with a scale from -1.0 to 1.0 in increments of 0.5.
* **Vertical Dividers:** Vertical black lines at x = -0.5 and x = 0.5, separating the regions.
* **Regions:**
* Left: |ψ0^(x=-1)>
* Middle: |ψ0^(x=0)>
* Right: |ψ0^(x=1)>
* **Lattice Diagram:** A 3x3 grid of nodes connected by lines, centered below the x-axis.
* **Formula:** H = (1/4q) * Σ(i,j) (ZiZj + x(XiXj + YiYj)) located at the top of the chart.
* **Legend (Top-Left):**
* Blue: <H>
* Orange: E0^(L(4))
* Green: E0
**Chart (b): Correlation vs. Position**
* **Title:** (b)
* **Y-axis:** Correlation, with a scale from 0.40 to 1.00 in increments of 0.20.
* **X-axis:** x, with a scale from -1.0 to 1.0 in increments of 0.5.
* **Vertical Dividers:** Vertical black lines at x = -0.5 and x = 0.5, separating the regions.
* **Regions:**
* Left: |ψ0^(x=-1)>
* Middle: |ψ0^(x=0)>
* Right: |ψ0^(x=1)>
* **Lattice Diagram:** A 3x3 grid of nodes connected by lines, with two nodes highlighted in red.
* **Formula:** Cij = ZiZj located near the lattice diagram.
* **Legend (Top-Left):**
* Blue: <C>
* Orange: <C>0^(L(4))
* Green: <C>0
### Detailed Analysis
**Chart (a): Energy vs. Position**
* **<H> (Blue):** The line is constant at approximately -0.36 in the region -0.5 < x < 0.5. It decreases linearly from approximately -0.36 at x = -0.5 to approximately -0.56 at x = -1.0, and from approximately -0.36 at x = 0.5 to approximately -0.56 at x = 1.0.
* **E0^(L(4)) (Orange):** The line forms a curve that is approximately -0.36 at x = -0.5 and x = 0.5, and reaches a maximum of approximately -0.35 at x = 0.0. It decreases linearly from approximately -0.36 at x = -0.5 to approximately -0.56 at x = -1.0, and from approximately -0.36 at x = 0.5 to approximately -0.56 at x = 1.0.
* **E0 (Green):** The line forms a curve that is approximately -0.39 at x = -0.5 and x = 0.5, and reaches a maximum of approximately -0.36 at x = 0.0. It decreases linearly from approximately -0.39 at x = -0.5 to approximately -0.56 at x = -1.0, and from approximately -0.39 at x = 0.5 to approximately -0.56 at x = 1.0.
**Chart (b): Correlation vs. Position**
* **<C> (Blue):** The line is constant at approximately 1.0 in the region -0.5 < x < 0.5. It decreases linearly from approximately 1.0 at x = -0.5 to approximately 0.3 at x = -1.0, and from approximately 1.0 at x = 0.5 to approximately 0.3 at x = 1.0.
* **<C>0^(L(4)) (Orange):** The line forms a curve that is approximately 1.0 at x = -0.5 and x = 0.5, and reaches a maximum of approximately 0.98 at x = 0.0. It decreases linearly from approximately 1.0 at x = -0.5 to approximately 0.3 at x = -1.0, and from approximately 1.0 at x = 0.5 to approximately 0.3 at x = 1.0.
* **<C>0 (Green):** The line forms a curve that is approximately 0.85 at x = -0.5 and x = 0.5, and reaches a maximum of approximately 0.95 at x = 0.0. It decreases linearly from approximately 0.85 at x = -0.5 to approximately 0.3 at x = -1.0, and from approximately 0.85 at x = 0.5 to approximately 0.3 at x = 1.0.
### Key Observations
* Both charts show a distinct change in behavior at x = -0.5 and x = 0.5, corresponding to the boundaries of the |ψ0^(x=0)> region.
* Within the |ψ0^(x=0)> region, the energy values are relatively constant for <H>, while the correlation values are constant for <C>.
* Outside the |ψ0^(x=0)> region, all data series exhibit a linear decrease as x moves away from 0.
* The lattice diagram in chart (b) highlights two nodes, suggesting a specific correlation (Cij) between them.
### Interpretation
The charts illustrate how energy and correlation change with position (x) under different conditions, represented by the states |ψ0^(x=-1)>, |ψ0^(x=0)>, and |ψ0^(x=1)>. The constant values within the |ψ0^(x=0)> region suggest a stable or uniform state, while the linear decrease outside this region indicates a change or decay in energy and correlation. The lattice diagrams provide a visual representation of the system being studied, and the highlighted nodes in chart (b) emphasize the correlation between specific elements within the lattice. The formulas provided give the mathematical relationships being explored. The differences between the blue, orange, and green lines in each chart likely represent different levels of approximation or different theoretical models used to calculate the energy and correlation.
</details>
Figure 1: (a) Ground state energy estimates for the Heisenberg XXZ model as a function of $x$ in the Hamiltonian. Green shows the exact ground state energy, $E_{0}$ . The energy estimates, $\langle H\rangle$ (blue) and $E_{0}^{\rm L(4)}$ (orange), are evaluated with respect to three trial states: the $x=-1$ (ferromagnetic) ground state for $-1<x<-\frac{1}{2}$ , the $x=0$ (Néel) ground state for $-\frac{1}{2}<x<\frac{1}{2}$ , and the $x=1$ (antiferromagnetic) ground state for $\frac{1}{2}<x<1$ . (b) Estimates of ground state $ZZ$ correlation between two next-nearest-neighbour spins (marked in red) versus $x$ , using the same three trial states. Like before, the exact correlation is shown in green, the direct expectation value in blue, and the moments-based estimate in orange.
<details>
<summary>x2.png Details</summary>

### Visual Description
## Chart: Quantum Measurement Expectation Values vs. Parameter g
### Overview
The image presents a series of line plots arranged in a 3x3 grid. Each plot displays the expectation values of quantum mechanical operators as a function of a parameter 'g'. The plots are organized by two parameters: 'p' (columns) and 'F' (rows). The x-axis represents 'g', and the y-axis represents the expectation value. Three different expectation values are plotted in each subplot, distinguished by color.
### Components/Axes
* **X-axis:** 'g', ranging from 0.00 to 2.00, with tick marks at intervals of 0.25.
* **Y-axis:** Expectation value, ranging from -0.5 to 0.0, with tick marks at intervals of 0.1.
* **Titles:**
* Top row: p = 0.01, p = 0.10, p = 0.50
* Left column: F ≈ 0.99, F ≈ 0.80, F ≈ 0.40
* **Legend (present in each subplot):**
* Blue: ⟨ψg t|M|ψg t⟩
* Orange: ⟨M⟩(4)0|ψg t⟩
* Green: ⟨ψg 0|M|ψg 0⟩
### Detailed Analysis
The data is presented in a 3x3 grid, where each cell contains a plot of three different expectation values as a function of 'g'. The columns represent different values of 'p' (0.01, 0.10, 0.50), and the rows represent different values of 'F' (approximately 0.99, 0.80, and 0.40).
**Row 1: F ≈ 0.99**
* **p = 0.01:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.45.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
* **p = 0.10:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.4.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
* **p = 0.50:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately -0.1 and decreases to around -0.2.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
**Row 2: F ≈ 0.80**
* **p = 0.01:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately 0.0 and fluctuates, ending around -0.4.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
* **p = 0.10:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately 0.0 and fluctuates, ending around -0.1.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
* **p = 0.50:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately -0.1 and fluctuates, ending around -0.2.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
**Row 3: F ≈ 0.40**
* **p = 0.01:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately 0.0 and fluctuates significantly, ending around -0.1.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
* **p = 0.10:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately 0.0 and fluctuates significantly, ending around -0.1.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
* **p = 0.50:**
* Blue (⟨ψg t|M|ψg t⟩): Starts at approximately -0.1 and fluctuates significantly, ending around -0.2.
* Orange (⟨M⟩(4)0|ψg t⟩): Starts at approximately 0.0 and decreases to around -0.5, with significant fluctuations.
* Green (⟨ψg 0|M|ψg 0⟩): Starts at approximately 0.0 and decreases to around -0.5.
### Key Observations
* The blue line (⟨ψg t|M|ψg t⟩) shows the most variation with changes in 'p' and 'F'. It generally fluctuates more as 'F' decreases.
* The orange (⟨M⟩(4)0|ψg t⟩) and green (⟨ψg 0|M|ψg 0⟩) lines tend to converge and remain relatively stable, decreasing from 0.0 to approximately -0.5 as 'g' increases.
* Higher values of 'p' seem to result in less change in the blue line as 'g' increases, especially when 'F' is low.
### Interpretation
The plots illustrate how the expectation values of certain quantum mechanical operators change with the parameter 'g' under different conditions defined by 'p' and 'F'. The parameter 'F' seems to influence the stability of the system, with lower values leading to more fluctuations in the expectation value represented by the blue line. The parameter 'p' appears to affect the overall magnitude of the expectation value, particularly for the blue line. The convergence of the orange and green lines suggests that the corresponding expectation values are less sensitive to changes in 'p' and 'F' compared to the blue line. The data suggests that the system's behavior is complex and depends on the interplay between the parameters 'g', 'p', and 'F'.
</details>
Figure 2: Staggered magnetisation for a 6-site 1D Heisenberg antiferromagnet with staggered external field, as a function of field parameter $g$ in the Hamiltonian (Eq. 11). Comparison of methods for obtaining ground state staggered magnetisation – direct $\langle M\rangle$ (blue) and QCM $\langle M\rangle_{0}^{\rm L(4)}$ (orange) estimates vs. exact value (green). To compute the estimates, trial states $\ket{\psi_{t}^{g}}$ were generated at each value of $g$ by starting at the ground state $\ket{\psi_{0}^{g}}$ and applying a random rotation as $\ket{\psi_{t}^{g}}=e^{-i\theta M}\ket{\psi_{0}^{g}}$ , where $M$ is a random Hermitian matrix sampled from the Gaussian unitary ensemble. By tuning the parameter $\theta$ , trial states of varying degrees of overlap with the ground state $F=\left|\langle\psi_{t}^{g}|\psi_{0}^{g}\rangle\right|^{2}$ were generated. These trial states ${\rho}=\ket{\psi_{t}^{g}}\!\bra{\psi_{t}^{g}}$ were then passed through a depolarising error channel ${\rho}\mapsto(1-\frac{3p}{4}){\rho}+\frac{p}{4}(X{\rho}X+Y{\rho}Y+Z{\rho}Z)$ with noise parameter $p$ , which was varied from 1% up to 50%.
Of particular interest in the study of quantum magnetism is the Heisenberg XXZ model, given by the Hamiltonian:
$$
H=\frac{1}{4q}\sum_{\langle i\,j\rangle}\left(Z_{i}Z_{j}+x(X_{i}X_{j}+Y_{i}Y_{%
j})\right). \tag{9}
$$
Here, $q$ is the number of qubits, the $S_{i}$ are the Pauli operators $S∈\{X,Y,Z\}$ acting on the $i^{\rm th}$ qubit, the $\langle i\,j\rangle$ are the nearest neighbour couplings in the lattice. The parameter $x$ controls the strength of the transverse spin components relative to the Ising-like interaction. As a first test case of the moments-based ground state observable estimate (Eq. 8), we take the XXZ model over a 2D square lattice, and consider the $ZZ$ correlation operator,
$$
C_{ij}=Z_{i}Z_{j}, \tag{10}
$$
acting on a pair of next-nearest-neighbour qubits $i$ and $j$ . Focusing on these longer range interactions provides additional insights into the system behaviour beyond what is captured explicitly by the Hamiltonian.
On a quantum computer, to simulate and study the ground state behaviour (e.g. the energy and spin-spin correlations) of this system as the parameter $x$ varies, one might need to prepare the ground state using a different circuit at each value of $x$ . However, one can perform far fewer computations by leveraging the fact that the QCM approach has reduced sensitivity to the choice of trial state. A well-chosen state that accurately describes the solution at a given value of $x$ will also produce an accurate QCM estimate for nearby values of $x$ . Figure 1 (a) shows that the ground state energy of the Heisenberg XXZ model can be fully characterised for all $-1<x<1$ using the QCM estimate $E_{0}^{\rm L(4)}$ with respect to just three different trial states (ground states at $x∈\{-1,0,1\}$ ), in stark contrast to $\langle H\rangle$ from the standard variational approach.
In Figure 1 (b), the difference between the two methods is made even clearer when these states are used to find estimates of other ground state properties. QCM correctly estimates the correlation in the ground state between the corner site and its next-nearest-neighbour in the middle of the $4× 3$ square lattice, remaining remarkably close to the exact value as $x$ varies, diverging only slightly at $x=±\frac{1}{2}$ where the ground state overlap is at a minimum.
In practice on a real device, measuring the moments up to $H^{4}$ requires one to multiply out the terms (or Pauli strings) in the Hamiltonian and determine their expectation values. These Pauli strings can be grouped into simultaneously measurable tensor product basis (TPB) sets, the total number of which scales logarithmically with the number of Pauli strings [12]. Each TPB corresponds to one measurement on the quantum computer – thus when determining a general ground state observable using QCM, if the operator describes interactions present explicitly in the Hamiltonian, no additional measurements need to be made over the energy calculation. In our Heisenberg XXZ case, to go from the energy calculation in Figure 1 (a) to the next-nearest-neighbour correlation in Figure 1 (b), the total number of Pauli strings increases from 66 343 to 68 960, and the total number of TPB measurements goes from 1 906 to 1 973. Note that these measurements need to only be made once for each of the three trial states, after which the entire curve for $-1<x<1$ can be filled out via simple classical post-processing.
So far, we have shown that the QCM estimate of a ground state observable other than energy shows the same resilience to an imperfect choice of trial state, however the presence of the noise robustness property still remains to be verified. For this, we turn our attention to the antiferromagnetic Heisenberg model with an external (staggered) field:
$$
H=\frac{1}{4}\sum_{\langle i\,j\rangle}(X_{i}X_{j}+Y_{i}Y_{j}+Z_{i}Z_{j})+%
\frac{g}{2}\sum_{k}(-1)^{k}Z_{k}, \tag{11}
$$
where $g$ is the external field parameter. An operator of interest in this case is the staggered magnetisation $M$ , given by:
$$
M=\frac{1}{2}\sum_{k}(-1)^{k}Z_{k}. \tag{12}
$$
We look at this model over 6 qubits in 1D. To simulate the results one might typically see in a variational quantum algorithm, given that preparing ansaẗze that can fully represent the target state is often challenging, we approximate an intermediate step of such an algorithm by starting with the exact ground state vector $\ket{\psi_{0}^{g}}$ for each $g$ and applying random rotations to obtain a “trial state” $\ket{\psi_{t}^{g}}$ . Both the direct and moments-based estimates of the staggered magnetisation are then calculated with respect to this trial state after applying various levels of depolarising noise. Figure 2 shows the interplay between these two phenomena in our simulations for varying average ground state fidelity, $F=\left|\langle\psi_{t}^{g}|\psi_{0}^{g}\rangle\right|^{2}$ , and depolarising noise, ${\rho}\mapsto(1-\frac{3p}{4}){\rho}+\frac{p}{4}(X{\rho}X+Y{\rho}Y+Z{\rho}Z)$ with parameter $p$ . We observe that the robustness of the QCM technique to low fidelity trial states and noise, previously only seen in the context of energy problems [12, 13, 14], persists when considering arbitrary ground state observables.
We see even for trial states with $\sim$ 40% ground state fidelity and 50% depolarising error, at each value of $g$ , the moments based estimate $\langle M\rangle_{0}^{\rm L(4)}$ provides a powerful correction to the direct measurement of $\langle M\rangle$ , remaining remarkably close to the exact value for the ground state staggered magnetisation. This is more or less the level of noise one would expect to see in a typical variational circuit on real hardware with 1% CNOT error, so the result bodes well for the prospect of applying this technique on present-day devices.
IV Conclusion
We have presented an extension of the QCM method to estimate arbitrary ground state observables of quantum systems, with a focus on the Heisenberg model. We observe the robustness of the QCM technique to low fidelity trial states and noise, previously only seen in the context of ground-state energy problems. In particular, the results presented in this paper demonstrate the ability of QCM to estimate ground state properties such as magnetisation and spin-spin correlations of the Heisenberg model to a high precision, even with suboptimal trial states and levels of noise present in near-term quantum hardware.
Further advancement of this technique will lead to improved understanding of the conditions under which various QCM observable estimates hold greater or lesser validity. We reiterate that the results presented here are from simplified noisy simulations, and the key next steps will involve performing the technique on real quantum hardware and under improved noise models – previous analysis of the QCM energy estimate in [13] in addition to the preliminary results presented here for the generalisation to arbitrary observables suggests that we should also expect this technique to work quite well on a real noisy device. Furthermore, the examples we have studied here have only scratched the surface of possible Hamiltonians and observables that could be suitable for this approach. Many intriguing problems invite future application of this technique, such as calculating ground state observables for the Fermi-Hubbard model [22] and detecting phase transitions in more unconventional systems by measuring nonlocal order parameters [23].
QCM has so far proven to be an exciting and potentially useful technique for obtaining the ground-state of many body systems on quantum computers, even in the presence of noise and with suboptimal trial states. Our findings here generalise the method to other observables of interest and further highlight its practical potential for solving a wide range of problems on near-term quantum hardware.
Acknowledgements
This research was supported by the University of Melbourne through the establishment of the IBM Quantum Network Hub at the University. HJV is supported by an Australian Government Research Training Program Scholarship.
References
- [1] Yudong Cao et al. “Quantum chemistry in the age of quantum computing” In Chemical reviews 119.19 ACS Publications, 2019, pp. 10856–10915 DOI: 10.1021/acs.chemrev.8b00803
- [2] Sam McArdle et al. “Quantum computational chemistry” In Reviews of Modern Physics 92.1 APS, 2020, pp. 015003 DOI: 10.1103/RevModPhys.92.015003
- [3] Marco Cerezo et al. “Variational quantum algorithms” In Nature Reviews Physics 3.9 Nature Publishing Group UK London, 2021, pp. 625–644 DOI: 10.1038/s42254-021-00348-9
- [4] Alberto Peruzzo et al. “A variational eigenvalue solver on a photonic quantum processor” In Nature communications 5.1 Nature Publishing Group, 2014, pp. 1–7 DOI: 10.1038/ncomms5213
- [5] Abhinav Kandala et al. “Hardware-efficient variational quantum eigensolver for small molecules and quantum magnets” In Nature 549.7671 Nature Publishing Group, 2017, pp. 242–246 DOI: 10.1038/nature23879
- [6] Yunseong Nam et al. “Ground-state energy estimation of the water molecule on a trapped-ion quantum computer” In npj Quantum Information 6.1 Nature Publishing Group UK London, 2020, pp. 33 DOI: 10.1038/s41534-020-0259-3
- [7] Samson Wang et al. “Noise-induced barren plateaus in variational quantum algorithms” In Nature communications 12.1 Nature Publishing Group, 2021, pp. 1–11 DOI: 10.1038/s41467-021-27045-6
- [8] Daniel Stilck França and Raul Garcia-Patron “Limitations of optimization algorithms on noisy quantum devices” In Nature Physics 17.11 Nature Publishing Group UK London, 2021, pp. 1221–1227 DOI: 10.1038/s41567-021-01356-3
- [9] Lloyd C L Hollenberg “Plaquette expansion in lattice Hamiltonian models” In Physical Review D 47.4 APS, 1993, pp. 1640 DOI: 10.1103/PhysRevD.47.1640
- [10] Lloyd C L Hollenberg and N S Witte “General nonperturbative estimate of the energy density of lattice Hamiltonians” In Physical Review D 50.5 APS, 1994, pp. 3382 DOI: 10.1103/PhysRevD.50.3382
- [11] Lloyd C L Hollenberg and N S Witte “Analytic solution for the ground-state energy of the extensive many-body problem” In Physical Review B 54.23 APS, 1996, pp. 16309 DOI: 10.1103/PhysRevB.54.16309
- [12] Harish J Vallury, Michael A Jones, Charles D Hill and Lloyd C L Hollenberg “Quantum computed moments correction to variational estimates” In Quantum 4 Verein zur Förderung des Open Access Publizierens in den Quantenwissenschaften, 2020, pp. 373 DOI: 10.22331/q-2020-12-15-373
- [13] Harish J. Vallury et al. “Noise-robust ground state energy estimates from deep quantum circuits” In Quantum 7 Verein zur Förderung des Open Access Publizierens in den Quantenwissenschaften, 2023, pp. 1109 DOI: 10.22331/q-2023-09-11-1109
- [14] Michael A Jones, Harish J Vallury, Charles D Hill and Lloyd C L Hollenberg “Chemistry beyond the Hartree–Fock energy via quantum computed moments” In Scientific Reports 12.1 Nature Publishing Group, 2022, pp. 1–9 DOI: 10.1038/s41598-022-12324-z
- [15] Karol Kowalski and Bo Peng “Quantum simulations employing connected moments expansions” In The Journal of Chemical Physics 153.20 AIP Publishing LLC, 2020, pp. 201102 DOI: 10.1063/5.0030688
- [16] Kazuhiro Seki and Seiji Yunoki “Quantum power method by a superposition of time-evolved states” In PRX Quantum 2.1 APS, 2021, pp. 010333 DOI: 10.1103/PRXQuantum.2.010333
- [17] Philippe Suchsland et al. “Algorithmic error mitigation scheme for current quantum processors” In Quantum 5 Verein zur Förderung des Open Access Publizierens in den Quantenwissenschaften, 2021, pp. 492 DOI: 10.22331/q-2021-07-01-492
- [18] Joseph C Aulicino, Trevor Keen and Bo Peng “State preparation and evolution in quantum computing: A perspective from Hamiltonian moments” In International Journal of Quantum Chemistry 122.5 Wiley Online Library, 2022, pp. e26853 DOI: 10.1002/qua.26853
- [19] Ruizhe Zhang, Guoming Wang and Peter Johnson “Computing ground state properties with early fault-tolerant quantum computers” In Quantum 6 Verein zur Förderung des Open Access Publizierens in den Quantenwissenschaften, 2022, pp. 761 DOI: 10.22331/q-2022-07-11-761
- [20] Lloyd C L Hollenberg and Michael J Tomlinson “Staggered magnetisation in the Heisenberg antiferromagnet” In Australian journal of physics 47.2 CSIRO Publishing, 1994, pp. 137–144 DOI: 10.1071/PH940137
- [21] N S Witte, Lloyd C L Hollenberg and Zheng Weihong “Two-dimensional XXZ model ground-state properties using an analytic Lanczos expansion” In Physical Review B 55.16 APS, 1997, pp. 10412 DOI: 10.1103/PhysRevB.55.10412
- [22] Stasja Stanisic et al. “Observing ground-state properties of the Fermi-Hubbard model using a scalable algorithm on a quantum computer” In Nature Communications 13.1 Nature Publishing Group UK London, 2022, pp. 5743 DOI: 10.1038/s41467-022-33335-4
- [23] Ken N Okada, Keita Osaki, Kosuke Mitarai and Keisuke Fujii “Identification of topological phases using classically-optimized variational quantum eigensolver” In arXiv preprint arXiv:2202.02909, 2022 DOI: 10.48550/arXiv.2202.02909