To obtain estimates of electronic energies, the Variational Quantum Eigensolver (VQE) technique performs separate measurements for multiple parts of the system Hamiltonian. Current quantum hardware is restricted to projective single-qubit measurements, and, thus, only parts of the Hamiltonian that form mutually qubit-wise commuting groups can be measured simultaneously. The number of such groups in the electronic structure Hamiltonians grows as N 4, where N is the number of qubits, thereby putting serious restrictions on the size of the systems that can be studied. Using a partitioning of the system Hamiltonian as a linear combination of unitary operators, we found a circuit formulation of the VQE algorithm that allows one to measure a group of fully anticommuting terms of the Hamiltonian in a single series of single-qubit measurements. Numerical comparison of the unitary partitioning to previously used grouping of Hamiltonian terms based on their qubit-wise commutativity is consistent with an N -fold reduction in the number of measurable groups.
Variational quantum algorithms are promising applications of noisy intermediate-scale quantum (NISQ) computers. These algorithms consist of a number of separate prepare-and-measure experiments that estimate terms in a Hamiltonian. The number of terms can become overwhelmingly large for problems at the scale of NISQ hardware that may soon be available. We use unitary partitioning (developed independently by Izmaylov et al. [J. Chem. Theory Comput. 16, 190 (2020)]) to define variational quantum eigensolver procedures in which additional unitary operations are appended to the ansatz preparation to reduce the number of terms. This approach may be scaled to use all coherent resources available after ansatz preparation. We also study the use of asymmetric qubitization to implement the additional coherent operations with lower circuit depth. Using this technique, we find a constant factor speedup for lattice and random Pauli Hamiltonians. For electronic structure Hamiltonians, we prove that linear term reduction with respect to the number of orbitals, which has been previously observed in numerical studies, is always achievable. For systems represented on 10--30 qubits, we find that there is a reduction in the number of terms by approximately an order of magnitude. Applied to the plane-wave dual-basis representation of fermionic Hamiltonians, however, unitary partitioning offers only a constant factor reduction. Finally, we show that noncontextual Hamiltonians may be reduced to effective commuting Hamiltonians using unitary partitioning.
Current implementations of the Variational Quantum Eigensolver (VQE) technique for solving the electronic structure problem involve splitting the system qubit Hamiltonian into parts whose elements commute within their single qubit subspaces. The number of such parts rapidly grows with the size of the molecule, this increases the uncertainty in the measurement of the energy expectation value because elements from different parts need to be measured independently. To address this problem we introduce a more efficient partitioning of the qubit Hamiltonian using fewer parts that need to be measured separately. The new partitioning scheme is based on two ideas: 1) grouping terms into parts whose eigenstates have a single-qubit product structure, and 2) devising multi-qubit unitary transformations for the Hamiltonian or its parts to produce less entangled operators. The first condition allows the new parts to be measured in the number of involved qubit consequential one-particle measurements. Advantages of the new partitioning scheme resulting in several-fold reduction of separately measured terms are illustrated on the H2 and LiH problems.
The contextual subspace variational quantum eigensolver (CS-VQE) is a hybrid quantum-classical algorithm that approximates the ground-state energy of a given qubit Hamiltonian. It achieves this by separating the Hamiltonian into contextual and noncontextual parts. The ground-state energy is approximated by classically solving the noncontextual problem, followed by solving the contextual problem using VQE, constrained by the noncontextual solution. In general, computation of the contextual correction needs fewer qubits and measurements compared with solving the full Hamiltonian via traditional VQE. We simulate CS-VQE on different tapered molecular Hamiltonians and apply the unitary partitioning measurement reduction strategy to further reduce the number of measurements required to obtain the contextual correction. Our results indicate that CS-VQE combined with measurement reduction is a promising approach to allow feasible eigenvalue computations on noisy intermediate-scale quantum devices. We also provide a modification to the CS-VQE algorithm; the CS-VQE algorithm previously could cause an exponential increase in Hamiltonian terms but with this modification now at worst will scale quadratically.
Near-term quantum computers will be limited in the number of qubits on which they can process information as well as the depth of the circuits that they can coherently carry out. To date, experimental demonstrations of algorithms such as the Variational Quantum Eigensolver (VQE) have been limited to small molecules using minimal basis sets for this reason. In this work we propose incorporating an orbital optimization scheme into quantum eigensolvers wherein a parametrized partial unitary transformation is applied to the basis functions set in order to reduce the number of qubits required for a given problem. The optimal transformation is found by minimizing the ground state energy with respect to this partial unitary matrix. Through numerical simulations of small molecules up to 16 spin orbitals, we demonstrate that this method has the ability to greatly extend the capabilities of near-term quantum computers with regard to the electronic structure problem. We find that VQE paired with orbital optimization consistently achieves lower ground state energies than traditional VQE when using the same number of qubits and even frequently achieves lower ground state energies than VQE methods using more qubits.
Estimating the expectation value of an observable
O for a quantum state evolved under a general, possibly non-unitary, operator
A is a fundamental task in quantum simulation, chemistry, and variational algorithms. Existing linear-combination-of-unitaries (LCU) approaches implement
A on hardware, requiring ancilla qubits, controlled unitaries and, in randomized variants, fresh circuit executions for every Monte Carlo sample. These costs are prohibitive on NISQ and early fault-tolerant devices. We propose HERMES (Hardware-Efficient Reusable Measurements for Estimating LCU-based expectation valueS), which moves the processing of
A entirely to the classical domain. Expanding
A and
O in the Pauli basis and sampling pairs of Pauli terms classically, HERMES rewrites
μ=Tr(OAρA†)/Tr(AρA†) as a ratio of weighted sums of Pauli expectation values on the unmodified state
ρ, with no ancillas, controlled unitaries, or quantum implementation of
A. A two-level measurement reduction scheme (Pauli deduplication and commuting-group partitioning) bounds the number of distinct measurement settings independently of the sample count, unlike randomized LCU methods whose circuit executions grow linearly with it. On ground-state energy estimation of
H2 and
LiH via imaginary-time evolution, HERMES matches or beats Pauli single-ancilla LCU baseline using shallow circuits with no two-qubit gates beyond state preparation under qubit-wise commuting grouping. It is also more robust to noise and completes a noisy
LiH simulation infeasible for the baseline, at the cost of classical overhead growing with the number of Pauli terms in
A.