Michael J. Biercuk

Large ion crystals in Penning traps provide a platform for quantum simulation and sensing with hundreds of spins, but their continuous rigid-body rotation has so far limited flexible local qubit control. Here we demonstrate programmable site-selective spin control across large rotating ${}^{9}\mathrm{Be}^{+}$ crystals in a Penning trap. A tightly focused off-resonant laser beam drives local $R_z$ phase rotations via differential AC Stark shifts. Beam steering synchronised with crystal rotation enables addressing of arbitrary ions throughout the crystal. Ramsey-based characterisation shows $R_z(π)$ gate fidelity of 94.6% and nearest-neighbour crosstalk of 1.2%. We use this capability to prepare spatially structured spin patterns, generating a biskyrmion spin texture in a single-layer crystal, then extending the method to bilayer crystals we perform layer-selective addressing operations. We further demonstrate dual-quadrature Ramsey sensing by imprinting a relative $π/2$ phase shift between spatial sub-ensembles, enabling simultaneous measurement of orthogonal spin components within a single experimental realisation. These results establish programmable local control in large rotating ion crystals, opening new routes for engineering spatially structured quantum states in multidimensional trapped-ion systems.

Quantum simulation with controllable many-body platforms offers a powerful route to exploring complex phases and dynamics that are difficult to access in natural materials. Among these, topological spin textures such as skyrmions are central to modern condensed-matter physics and play a key role in chiral quantum many-body systems. Their controlled realisation in large, programmable quantum platforms, however, remains an outstanding challenge. Here, we report deterministic generation and site-resolved reconstruction of topological spin textures in a two-dimensional crystal of more than 150 trapped ions. Using globally applied spin-dependent forces, we generate skyrmion configurations and reconstruct the full vector spin field with single-ion resolution, obtaining a winding number of 0.99$\pm$0.02 and a mean local fidelity of 0.87$\pm$0.04. In addition, we implement single-ion-resolved control to deterministically prepare domain-wall states, extending our approach to a broader class of non-uniform spin textures. These results establish trapped-ion crystals as a platform for engineering complex spin textures and open the door to exploring topology-dependent nonequilibrium dynamics in long-range interacting quantum systems.

Bosonic modes are prevalent in all aspects of quantum information processing. However, existing tools for characterizing the quality, stability, and noise properties of bosonic modes are limited, especially in a driven setting. Here, we propose, demonstrate, and analyze a bosonic randomized benchmarking (BRB) protocol that uses randomized displacements of the bosonic modes in phase space to determine their quality. We investigate the impact of common analytic error models, such as heating and dephasing, on the distribution of outcomes over randomized displacement trajectories in phase space. We show that analyzing the distinctive behavior of the mean and variance of this distribution - describable as a gamma distribution - enables identification of error processes, and quantitative extraction of error rates and correlations using a minimal number of measurements. We experimentally validate the analytical models by injecting engineered noise into the motional mode of a trapped ion system and performing the bosonic randomized benchmarking protocol, showing good agreement between experiment and theory. Finally, we investigate the intrinsic error properties in our system, identifying the presence of highly correlated dephasing noise as the dominant process.

Experimental implementations of quantum simulation must balance the controllability of the quantum system under test with decoherence typically introduced through interaction with external control fields. The ratio of coherent interaction strength to decoherence induced by stimulated emission in atomic systems is typically determined by hardware constraints, limiting the flexibility needed to explore different operating regimes. Here, we present an optomechanical system for in-situ tuning of the coherent spin-motion and spin-spin interaction strength in two-dimensional ion crystals confined in a Penning trap. The system introduces active optical positioners into the tightly constrained space of the bore of a superconducting magnet, allowing adjustability of the key hardware parameter which determines the ratio of coherent to incoherent light-matter interaction for fixed optical power. Using precision closed-loop piezo-actuated positioners, the system permits in-situ tuning of the angle-of-incidence of laser beams incident on the ion crystal up to $θ_{\text{ODF}}\approx 28^\circ$. We characterize the system using measurements of the induced mean-field spin precession under the application of an optical dipole force in ion crystals cooled below the Doppler limit through electromagnetically induced transparency cooling. These experiments show approximately a $\times2$ variation in the ratio of the coherent to incoherent interaction strength with changing $θ_{\text{ODF}}$, consistent with theoretical predictions. We characterize system stability over 6000 seconds; rigid mounting of optomechanics to the ion trap structure reduces differential laser movements to approximately $2\times 10^{-3}$ degrees per hour, enabling long-duration experiments. These technical developments will be crucial in future quantum simulations and sensing applications.

Resolving the locations and discriminating the spin states of individual trapped ions with high fidelity is critical for a large class of applications in quantum computing, simulation, and sensing. We report on a method for high-fidelity state discrimination in large two-dimensional (2D) crystals with over 100 trapped ions in a single trapping region, combining a hardware detector and an artificial neural network. A high-data-rate, spatially resolving, single-photon sensitive timestamping detector performs efficient single-shot detection of 2D crystals in a Penning trap, exhibiting rotation at about $25\,\mathrm{kHz}$. We then train an artificial neural network to process the fluorescence photon data in the rest frame of the rotating crystal in order to identify ion locations with a success rate of $~90\%$, accounting for substantial illumination inhomogeneity across the crystal. Finally, employing a time-binned state detection method, we arrive at an average spin-state detection fidelity of $94(2)\%$. This technique can be used to analyze spatial and temporal correlations in arrays of hundreds of trapped-ion qubits.

Identifying and calibrating quantitative dynamical models for physical quantum systems is important for a variety of applications. Here we present a closed-loop Bayesian learning algorithm for estimating multiple unknown parameters in a dynamical model, using optimised experimental "probe" controls and measurement. The estimation algorithm is based on a Bayesian particle filter, and is designed to autonomously choose informationally-optimised probe experiments with which to compare to model predictions. We demonstrate the performance of the algorithm in both simulated calibration tasks and in an experimental single-qubit ion-trap system. Experimentally, we find that with 60x fewer samples, we exceed the precision of conventional calibration methods, delivering an approximately 93x improvement in efficiency (as quantified by the reduction of measurements required to achieve a target residual uncertainty and multiplied by the increase in accuracy). In simulated and experimental demonstrations, we see that successively longer pulses are selected as the posterior uncertainty iteratively decreases, leading to an exponential improvement in the accuracy of model parameters with the number of experimental queries.

Conical intersections are ubiquitous in chemistry and physics, often governing processes such as light harvesting, vision, photocatalysis, and chemical reactivity. They act as funnels between electronic states of molecules, allowing rapid and efficient relaxation during chemical dynamics. In addition, when a reaction path encircles a conical intersection, the molecular wavefunction experiences a geometric phase, which can affect the outcome of the reaction through quantum-mechanical interference. Past experiments have measured indirect signatures of geometric phases in scattering patterns and spectroscopic observables, but there has been no direct observation of the underlying wavepacket interference. Here, we experimentally observe geometric-phase interference in the dynamics of a wavepacket travelling around an engineered conical intersection in a programmable trapped-ion quantum simulator. To achieve this, we develop a technique to reconstruct the two-dimensional wavepacket densities of a trapped ion. Experiments agree with the theoretical model, demonstrating the ability of analog quantum simulators -- such as those realised using trapped ions -- to accurately describe nuclear quantum effects.

Spectroscopy is one of the most accurate probes of the molecular world. However, predicting molecular spectra accurately is computationally difficult because of the presence of entanglement between electronic and nuclear degrees of freedom. Although quantum computers promise to reduce this computational cost, existing quantum approaches rely on combining signals from individual eigenstates, an approach that is difficult to scale because the number of eigenstates grows exponentially with molecule size. Here, we introduce a method for scalable analog quantum simulation of molecular spectroscopy, by performing simulations in the time domain. Our approach can treat more complicated molecular models than previous ones, requires fewer approximations, and can be extended to open quantum systems with minimal overhead. We present a direct mapping of the underlying problem of time-domain simulation of molecular spectra to the degrees of freedom and control fields available in a trapped-ion quantum simulator. We experimentally demonstrate our algorithm on a trapped-ion device, exploiting both intrinsic electronic and motional degrees of freedom, showing excellent quantitative agreement for a single-mode vibronic photoelectron spectrum of SO$_2$.

Ultrafast chemical reactions are difficult to simulate because they involve entangled, many-body wavefunctions whose computational complexity grows rapidly with molecular size. In photochemistry, the breakdown of the Born-Oppenheimer approximation further complicates the problem by entangling nuclear and electronic degrees of freedom. Here, we show that analog quantum simulators can efficiently simulate molecular dynamics using commonly available bosonic modes to represent molecular vibrations. Our approach can be implemented in any device with a qudit controllably coupled to bosonic oscillators and with quantum hardware resources that scale linearly with molecular size, and offers significant resource savings compared to digital quantum simulation algorithms. Advantages of our approach include a time resolution orders of magnitude better than ultrafast spectroscopy, the ability to simulate large molecules with limited hardware using a Suzuki-Trotter expansion, and the ability to implement realistic system-bath interactions with only one additional interaction per mode. Our approach can be implemented with current technology; e.g., the conical intersection in pyrazine can be simulated using a single trapped ion. Therefore, we expect our method will enable classically intractable chemical dynamics simulations in the near term.

Quantum harmonic oscillators are central to many modern quantum technologies. We introduce a method to determine the frequency noise spectrum of oscillator modes through coupling them to a qubit with continuously driven qubit-state-dependent displacements. We reconstruct the noise spectrum using a series of different drive phase and amplitude modulation patterns in conjunction with a data-fusion routine based on convex optimization. We apply the technique to the identification of intrinsic noise in the motional frequency of a single trapped ion with sensitivity to fluctuations at the sub-Hz level in a spectral range from quasi-DC up to 50 kHz.

We study a class of entangling gates for trapped atomic ions and demonstrate the use of numeric optimization techniques to create a wide range of fast, error-robust gate constructions. Our approach introduces a framework for numeric optimization using individually addressed, amplitude and phase modulated controls targeting maximally and partially entangling operations on ion pairs, complete multi-ion registers, multi-ion subsets of large registers, and parallel operations within a single register. Our calculations and simulations demonstrate that the inclusion of modulation of the difference phase for the bichromatic drive used in the Mølmer-Sørensen gate permits approximately time-optimal control across a range of gate configurations, and when suitably combined with analytic constraints can also provide robustness against key experimental sources of error. We further demonstrate the impact of experimental constraints such as bounds on coupling rates or modulation band-limits on achievable performance. Using a custom optimization engine based on TensorFlow we also demonstrate time-to-solution for optimizations on ion registers up to 20 ions of order tens of minutes using a local-instance laptop, allowing computational access to system-scales relevant to near-term trapped-ion devices.

Manipulating quantum computing hardware in the presence of imperfect devices and control systems is a central challenge in realizing useful quantum computers. Susceptibility to noise limits the performance and capabilities of noisy intermediate-scale quantum (NISQ) devices, as well as any future quantum computing technologies. Fortunately quantum control enables efficient execution of quantum logic operations and algorithms with built-in robustness to errors, without the need for complex logical encoding. In this manuscript we introduce software tools for the application and integration of quantum control in quantum computing research, serving the needs of hardware R&D teams, algorithm developers, and end users. We provide an overview of a set of python-based classical software tools for creating and deploying optimized quantum control solutions at various layers of the quantum computing software stack. We describe a software architecture leveraging both high-performance distributed cloud computation and local custom integration into hardware systems, and explain how key functionality is integrable with other software packages and quantum programming languages. Our presentation includes a detailed mathematical overview of central product features including a flexible optimization toolkit, filter functions for analyzing noise susceptibility in high-dimensional Hilbert spaces, and new approaches to noise and hardware characterization. Pseudocode is presented in order to elucidate common programming workflows for these tasks, and performance benchmarking is reported for numerically intensive tasks, highlighting the benefits of the selected cloud-compute architecture. Finally, we present a series of case studies demonstrating the application of quantum control solutions using these tools in real experimental settings for both trapped-ion and superconducting quantum computer hardware.

The fragility of quantum systems makes them ideally suited for sensing applications at the nanoscale. However, interpreting the output signal of a qubit-based sensor is generally complicated by background clutter due to out-of-band spectral leakage, as well as ambiguity in signal origin when the sensor is operated with noisy hardware. We present a sensing protocol based on optimally band-limited "Slepian functions" that can overcome these challenges, by providing narrowband sensing of ambient dephasing noise, coupling additively to the sensor along the ${z}$-axis, while permitting isolation of the target noise spectrum from other contributions coupling along a different axis. This is achieved by introducing a finite-difference control modulation, which linearizes the sensor's response and affords tunable band-limited "windowing" in frequency. Building on these techniques, we experimentally demonstrate two spectral estimation capabilities using a trapped-ion qubit sensor. We first perform efficient experimental reconstruction of a "mixed" dephasing spectrum, composed of a broadband $1/f$-type spectrum with discrete spurs. We then demonstrate the simultaneous reconstruction of overlapping dephasing and control noise spectra from a single set of measurements, in a setting where the two noise sources contribute equally to the sensor's response. Our approach provides a direct means to augment quantum-sensor performance in the presence of both complex broadband noise environments and imperfect control signals, by optimally complying with realistic time-bandwidth constraints.

New quantum computing architectures consider integrating qubits as sensors to provide actionable information useful for decoherence mitigation on neighboring data qubits, but little work has addressed how such schemes may be efficiently implemented in order to maximize information utilization. Techniques from classical estimation and dynamic control, suitably adapted to the strictures of quantum measurement, provide an opportunity to extract augmented hardware performance through automation of low-level characterization and control. In this work, we present an autonomous learning framework, Noise Mapping for Quantum Architectures (NMQA), for adaptive scheduling of sensor-qubit measurements and efficient spatial noise mapping (prior to actuation) across device architectures. Via a two-layer particle filter, NMQA receives binary measurements and determines regions within the architecture that share common noise processes; an adaptive controller then schedules future measurements to reduce map uncertainty. Numerical analysis and experiments on an array of trapped ytterbium ions demonstrate that NMQA outperforms brute-force mapping by up-to $18$x ($3$x) in simulations (experiments), calculated as a reduction in the number of measurements required to map a spatially inhomogeneous magnetic field with a target error metric. As an early adaptation of robotic control to quantum devices, this work opens up exciting new avenues in quantum computer science.

Entangling operations are among the most important primitive gates employed in quantum computing and it is crucial to ensure high-fidelity implementations as systems are scaled up. We experimentally realize and characterize a simple scheme to minimize errors in entangling operations related to the residual excitation of mediating bosonic oscillator modes that both improves gate robustness and provides scaling benefits in larger systems. The technique employs discrete phase shifts in the control field driving the gate operation, determined either analytically or numerically, to ensure all modes are de-excited at arbitrary user-defined times. We demonstrate an average gate fidelity of 99.4(2)% across a wide range of parameters in a system of $^{171}\text{Yb}^{+}$ trapped ion qubits, and observe a reduction of gate error in the presence of common experimental error sources. Our approach provides a unified framework to achieve robustness against both static and time-varying laser amplitude and frequency detuning errors. We verify these capabilities through system-identification experiments revealing improvements in error-susceptibility achieved in phase-modulated gates.

We present the design, construction and characterization of an experimental system capable of supporting a broad class of quantum simulation experiments with hundreds of spin qubits using Be-9 ions in a Penning trap. This article provides a detailed overview of the core optical and trapping subsystems, and their integration. We begin with a description of a dual-trap design separating loading and experimental zones and associated vacuum infrastructure design. The experimental-zone trap electrodes are designed for wide-angle optical access (e.g. for lasers used to engineer spin-motional coupling across large ion crystals) while simultaneously providing a harmonic trapping potential. We describe a near-zero-loss liquid-cryogen-based superconducting magnet, employed in both trapping and establishing a quantization field for ion spin-states, and equipped with a dual-stage remote-motor LN2LHe recondenser. Experimental measurements using a nuclear magnetic resonance (NMR) probe demonstrate part-per-million homogeneity over 7 mm-diameter cylindrical volume, with no discernible effect on the measured NMR linewidth from pulse-tube operation. Next we describe a custom-engineered inbore optomechanical system which delivers ultraviolet (UV) laser light to the trap, and supports multiple aligned optical objectives for top- and sideview imaging in the experimental trap region. We describe design choices including the use of non-magnetic goniometers and translation stages for precision alignment. Further, the optomechanical system integrates UV-compatible fiber optics which decouple the system's alignment from remote light sources. Using this system we present site-resolved images of ion crystals and demonstrate the ability to realize both planar and three-dimensional ion arrays via control of rotating wall electrodes and radial laser beams. Looking to future work, we include interferometric..

We demonstrate stabilisation of an ultraviolet diode laser via Doppler free spectroscopy of Ytterbium ions in a discharge. Our technique employs polarization spectroscopy, which produces a natural dispersive lineshape whose zero-crossing is largely immune to environmental drifts, making this signal an ideal absolute frequency reference for Yb$^+$ ion trapping experiments. We stabilise an external-cavity diode laser near 369 nm for cooling Yb$^+$ ions, using amplitude-modulated polarisation spectroscopy and a commercial PID feedback system. We achieve stable, low-drift locking with a standard deviation of measured laser frequency ~400 kHz over 10 minutes, limited by the instantaneous linewidth of the diode laser. These results and the simplicity of our optical setup makes our approach attractive for stabilization of laser sources in atomic-physics applications.

The presence of long-range quantum spin correlations underlies a variety of physical phenomena in condensed matter systems, potentially including high-temperature superconductivity. However, many properties of exotic strongly correlated spin systems (e.g., spin liquids) have proved difficult to study, in part because calculations involving N-body entanglement become intractable for as few as N~30 particles. Feynman divined that a quantum simulator - a special-purpose "analog" processor built using quantum particles (qubits) - would be inherently adept at such problems. In the context of quantum magnetism, a number of experiments have demonstrated the feasibility of this approach. However, simulations of quantum magnetism allowing controlled, tunable interactions between spins localized on 2D and 3D lattices of more than a few 10's of qubits have yet to be demonstrated, owing in part to the technical challenge of realizing large-scale qubit arrays. Here we demonstrate a variable-range Ising-type spin-spin interaction J_ij on a naturally occurring 2D triangular crystal lattice of hundreds of spin-1/2 particles (9Be+ ions stored in a Penning trap), a computationally relevant scale more than an order of magnitude larger than existing experiments. We show that a spin-dependent optical dipole force can produce an antiferromagnetic interaction J_ij ~ 1/d_ij^a, where a is tunable over 0<a<3; d_ij is the distance between spin pairs. These power-laws correspond physically to infinite-range (a=0), Coulomb-like (a=1), monopole-dipole (a=2) and dipole-dipole (a=3) couplings. Experimentally, we demonstrate excellent agreement with theory for 0.05<a<1.4. This demonstration coupled with the high spin-count, excellent quantum control and low technical complexity of the Penning trap brings within reach simulation of interesting and otherwise computationally intractable problems in quantum magnetism.

We demonstrate spectroscopy and thermometry of individual motional modes in a mesoscopic 2D ion array using entanglement-induced decoherence as a method of transduction. Our system is a $\sim$400 $μ$m-diameter planar crystal of several hundred $^9$Be$^+$ ions exhibiting complex drumhead modes in the confining potential of a Penning trap. Exploiting precise control over the $^9$Be$^+$ valence electron spins, we apply a homogeneous spin-dependent optical dipole force to excite arbitrary transverse modes with an effective wavelength approaching the interparticle spacing ($\sim$20 \nolinebreak$μ$m). Center-of-mass displacements below 1 nm are detected via entanglement of spin and motional degrees of freedom.

Building robust instruments capable of making interferometric measurements with precision beyond the standard quantum limit remains an important goal in many metrology laboratories. We describe here the basic concepts underlying spin squeezing experiments that allow one to surpass this limit. In priniciple it is possible to reach the so-called Heisenberg limit, which constitutes an improvement in precision by a factor $\sqrt{N}$, where $N$ is the number of particles on which the measurement is carried out. In particular, we focus on recent progress toward implementing spin squeezing with a cloud of beryllium ions in a Penning ion trap, via the geometric phase gate used more commonly for performing two-qubit entangling operations in quantum computing experiments.

We discuss the use of two-dimensional $^{9}$Be$^{+}$ ion crystals for experimental tests of quantum control techniques. Our primary qubit is the 124 GHz ground-state electron spin flip transition, which we drive using microwaves. An ion crystal represents a spatial ensemble of qubits, but the effects of inhomogeneities across a typical crystal are small, and as such we treat the ensemble as a single effective spin. We are able to initialize the qubits in a simple state and perform a projective measurement on the system. We demonstrate full control of the qubit Bloch vector, performing arbitrary high-fidelity rotations ($τ_π\sim$200 $μ$s). Randomized Benchmarking demonstrates an error per gate (a Pauli-randomized $π/2$ and $π$ pulse pair) of $8\pm1\times10^{-4}$. Ramsey interferometry and spin-locking measurements are used to elucidate the limits of qubit coherence in the system, yielding a typical free-induction decay coherence time of $T_{2}\sim$2 ms, and a limiting $T_{1ρ}\sim$688 ms. These experimental specifications make ion crystals in a Penning trap ideal candidates for novel experiments in quantum control. As such, we briefly describe recent efforts aimed at studying the error-suppressing capabilities of dynamical decoupling pulse sequences, demonstrating an ability to extend qubit coherence and suppress phase errors. We conclude with a discussion of future avenues for experimental exploration, including the use of additional nuclear-spin-flip transitions for effective multiqubit protocols, and the potential for Coulomb crystals to form a useful testbed for studies of large-scale entanglement.

One approach to maintaining phase coherence of qubits through dynamical decoupling consists of applying a sequence of Hahn spin-echo pulses. Recent studies have shown that, in certain noise environments, judicious choice of the delay times between these pulses can greatly improve the suppression of phase errors compared to traditional approaches. By enforcing a simple analytical condition, we obtain sets of dynamical decoupling sequences that are designed for optimized noise filtration and are spectrum-independent up to a single scaling factor set by the coherence time of the system. We demonstrate the efficacy of these sequences in suppressing phase errors through measurements on a model qubit system, $^{9}$Be$^{+}$ ions in a Penning trap. Our combined theoretical and experimental studies show that in high-frequency-dominated noise environments this approach may suppress phase errors orders of magnitude more efficiently than comparable techniques can.

We present a detailed experimental study of the Uhrig Dynamical Decoupling (UDD) sequence in a variety of noise environments. Our qubit system consists of a crystalline array of $^{9}$Be$^{+}$ ions confined in a Penning trap. We use an electron-spin-flip transition as our qubit manifold and drive qubit rotations using a 124 GHz microwave system. We study the effect of the UDD sequence in mitigating phase errors and compare against the well known CPMG-style multipulse spin echo as a function of pulse number, rotation axis, noise spectrum, and noise strength. Our results agree well with theoretical predictions for qubit decoherence in the presence of classical phase noise, accounting for the effect of finite-duration $π$ pulses. Finally, we demonstrate that the Uhrig sequence is more robust against systematic over/underrotation and detuning errors than is multipulse spin echo, despite the precise prescription for pulse-timing in UDD.