Quantum Simulation Lab

Monolithic three-dimensional (3D) Paul traps combine the high-precision microfabrication of two-dimensional (2D) chip traps with the deep trapping potentials and low heating rates characteristic of macroscopic 3D Paul traps, which are typically machined by traditional means and mechanically assembled. However, achieving low motional heating rates and optical access with a high numerical aperture (NA) while maintaining the high radio-frequency (RF) voltages required for trapping heavy ionic species, such as Yb$^{+}$ and Ba$^{+}$, remains a significant technical challenge. In this work, we present a fused-silica, monolithic segmented 3D Paul trap with an ion-electrode distance of 250 $μ$m, and stable operation at high RF voltages. We benchmark the performance of the trap using Yb$^{+}$ ions, demonstrating axially homogeneous trapping potentials spanning over 200 $μ$m about the axial center of the trap, high multi-directional optical access (up to 0.7 NA), and radial motional heating as low as $\dot{\bar n}=1.1 \pm 0.1 $ quanta/s at radial trap frequencies about 3 MHz near room temperature. Furthermore, we observe a motional Ramsey coherence time, ${T}_{2}$, of about 95 ms for the radial center-of-mass mode. We demonstrate the generation of a two-qubit Bell state with a parity contrast of ${99.3}^{+0.7} _{-1.5}$% with state preparation and measurement correction. These results establish fused-silica monolithic 3D Paul traps as a scalable, modular platform for quantum simulation, computation, metrology, and networking with heavy ionic species.

We propose and demonstrate an experimental scheme to engineer thermal baths with independently tunable temperatures and dissipation rates for the motional modes of a trapped-ion system. This approach enables robust thermal-state preparation and quantum simulations of open-system dynamics in bosonic and spin-boson models at well-controlled finite temperatures. We benchmark our protocol by experimentally realizing out-of-equilibrium dynamics of a charge-transfer model at different temperatures. We observe that, when the process occurs at a higher temperature, the transfer rate spectrum broadens, with reduced rates at small donor-acceptor energy gaps and enhanced rates at large gaps. We then employ our scheme to study local-temperature effects in a two-mode vibrationally assisted exciton transfer system, where we observe thermally activated interference pathways for excitation transfer.

We study a cost Hamiltonian compilation problem for the quantum approximate optimization algorithm (QAOA) applied to the Max-Cut problem, focusing on trapped-ion quantum computers. Instead of standard compilation with CNOT and Rz gates, we employ global coupling operations and single-qubit bit flips. Prior work by Rajakumar et al. established that such a compilation is always possible. Minimizing operational error requires short operation sequences. The problem reduces to a low-rank semi-discrete decomposition of the graph's adjacency matrix, where the minimum achievable rank, the graph coupling number gc(G), represents the number of global control layers. Rajakumar et al. introduced the Union of Stars construction, proving gc(G) <= 3n - 2 for unweighted graphs with n vertices, and gave an O(m)-rank construction for weighted graphs. We concentrate on unweighted graphs. We derive structural properties of the compilation problem and show the Union of Stars method is order-optimal by proving a lower bound of gc(G) >= n - 1 for a family of graphs. We also improve the general upper bound to 2.5n + 2. For particular graph families -- cliques, perfect matchings, paths, and cycles -- we provide sharper bounds. Further, we reveal a link between the problem and Hadamard matrix theory. Finally, we introduce a compact mixed-integer programming (MIP) formulation that outperforms the previously studied exponential-size MIP.

Quantum simulation offers a route to study open-system molecular dynamics in non-perturbative regimes by programming the interactions among electronic, vibrational, and environmental degrees of freedom on similar energy scales. Trapped-ion systems possess this capability, with their native spins, phonons, and tunable dissipation integrated within a single platform. Here, we demonstrate an open-system quantum simulation of charge and exciton transfer in a multi-mode linear vibronic coupling model. Employing tailored spin-phonon interactions alongside reservoir engineering techniques, we emulate a system with two dissipative vibrational modes coupled to donor and acceptor electronic sites and follow its non-equilibrium dynamics. We continuously tune the system from the charge transfer (CT) regime to the vibrationally assisted exciton transfer (VAET) regime by controlling the vibronic coupling strengths. We find that degenerate modes enhance CT and VAET rates at large energy gaps, while non-degenerate modes activate slow-mode pathways that reduce the energy-gap dependence, thus enlarging the window for efficient transfer. These results show that the presence of one additional vibration introduces interfering vibrationally assisted pathways and reshapes non-perturbative quantum excitation transfer. Our work establishes a scalable and hardware-efficient route to simulating chemically relevant, many-mode vibronic processes with engineered environments, guiding the design of next-generation organic photovoltaics and molecular electronics.

Electron transfer is at the heart of many fundamental physical, chemical, and biochemical processes essential for life. The exact simulation of these reactions is often hindered by the large number of degrees of freedom and by the essential role of quantum effects. Here, we experimentally simulate a paradigmatic model of molecular electron transfer using a multispecies trapped-ion crystal, where the donor-acceptor gap, the electronic and vibronic couplings, and the bath relaxation dynamics can all be controlled independently. By manipulating both the ground-state and optical qubits, we observe the real-time dynamics of the spin excitation, measuring the transfer rate in several regimes of adiabaticity and relaxation dynamics. Our results provide a testing ground for increasingly rich models of molecular excitation transfer processes that are relevant for molecular electronics and light-harvesting systems.

Combinatorial optimization problems on graphs have broad applications in science and engineering. The Quantum Approximate Optimization Algorithm (QAOA) is a method to solve these problems on a quantum computer by applying multiple rounds of variational circuits. However, there exist several challenges limiting the real-world applications of QAOA. In this paper, we demonstrate on a trapped-ion quantum computer that QAOA results improve with the number of rounds for multiple problems on several arbitrary graphs. We also demonstrate an advanced mixing Hamiltonian that allows sampling of all optimal solutions with predetermined weights. Our results are a step towards applying quantum algorithms to real-world problems.

The conventional framework for defining and understanding phases of matter requires thermodynamic equilibrium. Extensions to non-equilibrium systems have led to surprising insights into the nature of many-body thermalization and the discovery of novel phases of matter, often catalyzed by driving the system periodically. The inherent heating from such Floquet drives can be tempered by including strong disorder in the system, but this can also mask the generality of non-equilibrium phases. In this work, we utilize a trapped-ion quantum simulator to observe signatures of a non-equilibrium driven phase without disorder: the prethermal discrete time crystal (PDTC). Here, many-body heating is suppressed not by disorder-induced many-body localization, but instead via high-frequency driving, leading to an expansive time window where non-equilibrium phases can emerge. We observe a number of key features that distinguish the PDTC from its many-body-localized disordered counterpart, such as the drive-frequency control of its lifetime and the dependence of time-crystalline order on the energy density of the initial state. Floquet prethermalization is thus presented as a general strategy for creating, stabilizing and studying intrinsically out-of-equilibrium phases of matter.