Joint Quantum Institute

This tutorial offers a pedagogical guide to hybrid quantum processors that integrate discrete-variable (DV) qubits and continuous-variable (CV) oscillators. Aimed at computer scientists, engineers, and physicists, it provides an overview of the experimental, algorithmic, and architectural aspects of this novel and rapidly developing hardware model. Experimental realizations of this model include superconducting, trapped-ion, and neutral-atom platforms. By combining DV and CV components, hybrid oscillator-qubit processors enable a powerful new paradigm that offers complementary strengths for quantum control, error correction, computation, and simulation. Working toward the goal of a full-stack system connecting applications to CV-DV hardware, we define and formulate abstract machine models and instruction set architectures. These essential abstractions enable codesign of hardware and software, and resource estimation for exploring the potential of current and future hardware for computational and simulation tasks. Using these abstractions, we present both new and existing examples that illustrate the benefits of hybrid CV-DV processors relative to traditional DV-only hardware in computation as well as quantum simulation of physical models. Examples include algorithms for transferring states between DV and CV systems, performing the quantum Fourier transform, and simulation of lattice gauge theories. Relative to qubit-only hardware, the bosonic degrees of freedom natively available in hybrid architectures can substantially reduce the circuit complexity of simulations for physical models containing bosons. A key technique is the extension of quantum signal processing ideas to CV-DV systems. This work is intended to serve as a timely and comprehensive guide to this relatively unexplored yet promising approach to quantum computation and to provide a road map to guide future development.

Ensembles of nitrogen-vacancy (NV) centers in diamond are versatile quantum sensors with broad applications in the physical and life sciences. The concentration of neutral substitutional nitrogen ( [ N s 0 ]) strongly influences NV electronic spin coherence times, sensitivity, and optimal sensing strategies. Diamonds with [ N s 0 ] ∼ 1–10 ppm are a focus of recent material engineering efforts, with higher concentrations being favorable for continuous-wave optically detected magnetic resonance (CW-ODMR) and lower concentrations expected to benefit pulsed magnetometry techniques through extended NV spin coherence times and improved sensing duty cycles. In this work, we synthesize and characterize low- [ N s 0 ] ( ∼ 0.8 ppm), NV-enriched diamond material, engineered through low-strain chemical vapor deposition (CVD) growth on high-quality substrates, 12 C isotopic purification, and controlled electron irradiation and annealing. Our results demonstrate good strain homogeneity in diamonds grown on CVD substrates and spin-bath-limited NV dephasing times. By measuring NV spin and charge properties across a wide range of optical NV excitation intensity, we provide direct comparisons of photon-shot-noise-limited magnetic field sensitivity between the current low- [ N s 0 ] and previously studied higher- [ N s 0 ] ( ∼ 14 ppm) NV-diamond sensors. We show that low- [ N s 0 ] diamond can outperform higher- [ N s 0 ] diamond at moderate and low optical NV excitation intensity. Our results provide practical benchmarks and guidance for selecting NV-diamond sensors tailored to specific experimental constraints and sensing requirements.

Recent experimental progress in controlling open quantum systems enables the pursuit of mixed-state nonequilibrium quantum phases. We investigate whether open quantum systems hosting mixed-state symmetry-protected topological states as steady states retain this property under symmetric perturbations. Focusing on the decohered cluster state – a mixed-state symmetry-protected topological state protected by a combined strong and weak symmetry – we construct a parent Lindbladian that hosts it as a steady state. This Lindbladian can be mapped onto exactly solvable reaction-diffusion dynamics, even in the presence of certain perturbations, allowing us to solve the parent Lindbladian in detail and reveal previously-unknown steady states. Using both analytical and numerical methods, we find that typical symmetric perturbations cause strong-to-weak spontaneous symmetry breaking at arbitrarily small perturbations, destabilize the steady-state mixed-state symmetry-protected topological order. However, when perturbations introduce only weak symmetry defects, the steady-state mixed-state symmetry-protected topological order remains stable. Additionally, we construct a quantum channel which replicates the essential physics of the Lindbladian and can be efficiently simulated using only Clifford gates, Pauli measurements, and feedback.

We characterize the dynamical state of many-body bosonic and fermionic many-body models with intersite Gaussian couplings, on-site non-Gaussian interactions, and local dissipation comprising incoherent particle loss, particle gain, and dephasing. We first establish that, for fermionic systems, if the dephasing noise is larger than the non-Gaussian interactions, irrespective of the Gaussian coupling strength, the system state is a convex combination of Gaussian states at all times. Furthermore, for bosonic systems, we show that if the particle loss and particle gain rates are larger than the Gaussian intersite couplings, the system remains in a separable state at all times. Building on this characterization, we establish that at noise rates above a threshold, there exists a classical algorithm that can efficiently sample from the system state of both the fermionic and bosonic models. Finally, we show that, unlike fermionic systems, bosonic systems can evolve into states that are not convex Gaussian even when the dissipation is much higher than the on-site non-Gaussianity. Similarly, unlike bosonic systems, fermionic systems can generate entanglement even with noise rates much larger than the intersite couplings.

Autonomous quantum memories are a way to passively protect quantum information using engineered dissipation that creates an “always-on'' decoder. We analyze Markovian autonomous decoders that can be implemented with a wide range of qubit and bosonic error-correcting codes, and derive several upper bounds and a lower bound on the logical error rate in terms of correction and noise rates. These bounds suggest that, in general, there is always a correction rate, possibly size-dependent, above which autonomous memories exhibit arbitrarily long coherence times. For any given autonomous memory, size dependence of this correction rate is difficult to rule out: we point to common scenarios where autonomous decoders that stochastically implement active error correction must operate at rates that grow with code size. For codes with a threshold, we show that it is possible to achieve faster-than-polynomial decay of the logical error rate with code size by using superlogarithmic scaling of the correction rate. We illustrate our results with several examples. One example is an exactly solvable global dissipative toric code model that can achieve an effective logical error rate that decreases exponentially with the linear lattice size, provided that the recovery rate grows proportionally with the linear lattice size.

In analogy to the development of classical computing, the true potential of quantum devices at scale lies in their interconnection within a network. The “quantum internet” [1] is still a distant aspiration with the main limitation being the remote entanglement generation rate dropping off sharply with distance due to the fiber losses of the visible and UV photons involved. Still, there have been a number of successful demonstrations of remote entanglement with different platforms [2]–[5]. Ongoing efforts have focused on converting entangled photons emitted from atomic systems into fiber-compatible wavelengths, but this conversion process incurs efficiency loss, noise, and system complexity. In this conference paper, we present initial results from direct transmission of infrared ion-photons emitted from Strontium ions at 1,092 nm through a municipal dark fiber link. This wavelength falls into a low-loss window in fiber transmission that will allow for local area quantum networks without the frequency-conversion overhead. This comes at the cost of a smaller branching ratio, limiting the transmission rate. However, there is a distance range of a few kilometers where this is an advantageous choice, making this scheme highly promising for city-scale quantum networks. We stimulate trapped Strontium ions to emit single photons at 1,092 nm. After emission, the ion is left in one of two states, which form a metastable qubit (see diagram in figure 1).

We study the intravalley spin-polarized superconductivity in rhombohedral tetralayer graphene, which has been discovered experimentally in Han et al. [arXiv:2408.15233]. We construct a minimal model for the intravalley spin-polarized superconductivity, assuming a simplified anisotropic interaction that depends only on the angle between the incoming and outgoing momenta. Despite the absence of Fermi surface nesting, we show that superconductivity can emerge near the Van Hove singularity with the maximal ${T}_{c}$ near a bifurcation point of the peaks in the density of states. We identify the $p+ip$, $h+ih$, and the nodal $f$-wave pairings as the possible states, which are all pair density wave orders due to the intravalley nature. Furthermore, these pair density wave orders require a finite attractive threshold for superconductivity, resulting in a narrow stripe shape of superconducting region, consistent with experimental findings. We point out that the Kohn-Luttinger mechanism is a plausible explanation with a dominant $p+ip$ pairing. The possibility of realizing intravalley spin-polarized superconductivity in other rhombohedral graphene systems is also discussed.

Density functional theory (DFT) offers a desirable balance between quantitative accuracy and computational efficiency in practical many-electron calculations. Its central component, the exchange-correlation energy functional, has been approximated with increasing levels of complexity ranging from strictly local approximations to nonlocal and orbital dependent expressions with many tuned parameters. In this paper, we formulate a general way of rewriting complex density functionals using deep neural networks in a way that allows for simplified computation of Kohn-Sham potentials as well as higher functional derivatives through automatic differentiation, enabling access to highly nonlinear response functions and forces. These goals are achieved by using a recently developed class of robust neural network models capable of modeling functionals, as opposed to functions, with explicitly enforced spatial symmetries. Functionals treated in this way are then called and can be seamlessly integrated with existing DFT workflows. Tests are performed for a dataset featuring a large variety of molecular structures and popular meta–generalized gradient approximation density functionals, where we successfully eliminate orbital dependencies coming from the kinetic energy density, and discover a high degree of transferability to a variety of physical systems. The presented framework is general and could be extended to more complex orbital and energy dependent functionals as well as refined with specialized datasets. Published by the American Physical Society 2025

Abstract We study spin-dynamics and quantum magnetism with ultracold highly-magnetic atoms. In particular, we focus on the interactions among rare-earth atoms localized in a site of an optical-lattice potential, modeled as a cylindrically symmetric harmonic oscillator in the presence of a weak external magnetic field. The interactions between the atoms are modeled using multi-channel Hamiltonian containing multiple spin-spin and anisotropic spin-orbit interactions with strengths that depend on the separation between the atoms. We studied the eigenenergies of the atom pair in a site for different lattice-site geometries and magnetic field strengths. In parallel, we compared these energies to those found from a simplified approach, where the complex-collisional physics is replaced by a two-length-scale pseudo-potential containing a contact and the magnetic dipole-dipole interaction. The eigenenergies within this model can be computed analytically within the Born approximation as well as non-perturbatively for stronger contact interactions. We have shown that the pseudo-potential model can accurately represent the multi-channel Hamiltonian in certain parameter regimes of the shape or geometry of the site of an optical lattice. The pseudo-potential forms the starting point for many-body, condensed matter simulations involving many atom pairs in different sites of an optical lattice.