Stephan Eidenbenz

Variational Quantum Simulations of a Two-Dimensional Frustrated Transverse-Field Ising Model on a Trapped-Ion Quantum Computer

Ammar Kirmani [1], Elijah Pelofske [2,3], Andreas Bärtschi, Stephan Eidenbenz [3], Jian-Xin Zhu [1]

Abstract

Quantum computers are an ideal platform to study the ground state properties of strongly correlated systems due to the limitation of classical computing techniques particularly for systems exhibiting quantum phase transitions. While the error rates of Noisy Intermediate-Scale Quantum (NISQ) computers are still high, simulating strongly correlated systems on such devices and extracting information of possible phases may be within reach. The frustrated transverse-field Ising model (TFIM) is such a system with multiple ordered magnetic phases. In this study, we simulate a two-dimensional frustrated TFIM with next-nearest-neighbor spin-exchange interactions at zero temperature. The competition between the nearest-neighbor ferromagnetic and next-nearest-neighbor antiferromagnetic coupling gives rise to frustration in the system. Moreover, the presence of quantum fluctuations makes the ground-state phase profile even richer. We use the Variational Quantum Eigensolver (VQE) to compute the phases on a square lattice with periodic boundary conditions for a system of 16 sites (qubits). The trained VQE circuits are compared to exact diagonalization, allowing us to extract error measures of VQE. We focus on the ground-state phase transitions of this model, where VQE succeeds in finding the dominant magnetic phases. The optimized VQE circuits are then executed on the Quantinuum H1-1 trapped-ion quantum computer without using any error mitigation techniques. Our experiments show near perfect recovery of the magnetic phases of the frustrated model through ground-state energy, the energy derivative, and the spin correlation functions. Thus, we show that the trapped-ion quantum processor is able to achieve reliable simulations of a strongly correlated system within the limitations of the VQE approach.

High-Round QAOA for MAX $k$-SAT on Trapped Ion NISQ Devices

Elijah Pelofske [1], Andreas Bärtschi, John Golden [1], Stephan Eidenbenz [1]

Abstract

The Quantum Alternating Operator Ansatz (QAOA) is a hybrid classical-quantum algorithm that aims to sample the optimal solution(s) of discrete combinatorial optimization problems. We present optimized QAOA circuit constructions for sampling MAX $k$-SAT problems, specifically for $k=3$ and $k=4$. The novel $4$-SAT QAOA circuit construction we present uses measurement based uncomputation, followed by classical feed forward conditional operations. The QAOA circuit parameters for $3$-SAT are optimized via exact classical (noise-free) simulation, using HPC resources to simulate up to $20$ rounds on $10$ qubits. In order to explore the limits of current NISQ devices we execute these optimized QAOA circuits for random $3$-SAT test instances with clause-to-variable ratio $4$ on four trapped ion quantum computers: Quantinuum H1-1 (20 qubits), IonQ Harmony (11 qubits), IonQ Aria 1 (25 qubits), and IonQ Forte (30 qubits). The QAOA circuits that are executed include $n=10$ up to $p=20$, and $n=22$ for $p=1$ and $p=2$. The high round circuits use upwards of 9,000 individual gate instructions, making these some of the largest QAOA circuits executed on NISQ devices. Our main finding is that current NISQ devices perform best at low round counts (i.e., $p = 1,\ldots, 5$) and then -- as expected due to noise -- gradually start returning satisfiability truth assignments that are no better than randomly picked solutions as the number of QAOA rounds are further increased.

Scalable Experimental Bounds for Entangled Quantum State Fidelities

Shamminuj Aktar [1,2], Andreas Bärtschi, Abdel-Hameed A. Badawy [1], Stephan Eidenbenz [2]

Abstract

Estimating the state preparation fidelity of highly entangled states on noisy intermediate-scale quantum (NISQ) devices is important for benchmarking and application considerations. Unfortunately, exact fidelity measurements quickly become prohibitively expensive, as they scale exponentially as $O(3^N)$ for $N$-qubit states, using full state tomography with measurements in all Pauli bases combinations. However, Somma and others [PhysRevA.74.052302] established that the complexity could be drastically reduced when looking at fidelity lower bounds for states that exhibit symmetries, such as Dicke States and GHZ States. These bounds must still be tight enough for larger states to provide reasonable estimations on NISQ devices. For the first time and more than 15 years after the theoretical introduction, we report meaningful lower bounds for the state preparation fidelity of all Dicke States up to $N=10$ and all GHZ states up to $N=20$ on Quantinuum H1 ion-trap systems using efficient implementations of recently proposed scalable circuits for these states. Our achieved lower bounds match or exceed previously reported exact fidelities on superconducting systems for much smaller states. Furthermore, we provide evidence that for large Dicke States $D^N_{N/2}$, we may resort to a GHZ-based approximate state preparation to achieve better fidelity. This work provides a path forward to benchmarking entanglement as NISQ devices improve in size and quality.

Short-Depth Circuits for Dicke State Preparation

Andreas Bärtschi, Stephan Eidenbenz [1]

Abstract

We present short-depth circuits to deterministically prepare any Dicke state |Dn,k>, which is the equal-amplitude superposition of all n-qubit computational basis states with Hamming Weight k. Dicke states are an important class of entangled quantum states with a large variety of applications, and a long history of experimental creation in physical systems. On the other hand, not much is known regarding efficient scalable quantum circuits for Dicke state preparation on realistic quantum computing hardware connectivities. Here we present preparation circuits for Dicke states |Dn,k> with (i) a depth of O(k log(n/k)) for All-to-All connectivity (such as on current ion trap devices); (ii) a depth of O(k sqrt(n/k)) = O(sqrt(nk) for Grid connectivity on grids of size Omega(sqrt(n/s)) x O(sqrt(ns)) with s<=k (such as on current superconducting qubit devices). Both approaches have a total gate count of O(kn), need no ancilla qubits, and generalize to both the preparation and compression of symmetric pure states in which all non-zero amplitudes correspond to states with Hamming weight at most k. Thus our work significantly improves and expands previous state-of-the art circuits which had depth O(n) on a Linear Nearest Neighbor connectivity for arbitrary k (Fundamentals of Computation Theory 2019) and depth O(log n) on All-to-All connectivity for k=1 (Advanced Quantum Technologies 2019).