Within the histories formalism the decoherence functional is a formal tool to investigate the emergence of classicality in isolated quantum systems, yet an explicit evaluation of it from first principles has not been reported. We provide such an evaluation for up to five-time histories based on exact numerical diagonalization of the Schrödinger equation. We find a robust emergence of decoherence for slow and coarse observables of a generic random matrix model and extract a finite-size scaling law by varying the Hilbert space dimension over 4 orders of magnitude. Specifically, we conjecture and observe an exponential suppression of coherent effects as a function of the particle number of the system. This suggests a solution to the preferred basis problem of the many-worlds interpretation (or the set selection problem of the histories formalism) within a minimal theoretical framework without relying on environmentally induced decoherence, quantum Darwinism, Markov approximations, low-entropy initial states, or ensemble averages.

Despite the high measurement fidelity that can now be reached, the dispersive qubit readout of circuit quantum electrodynamics is plagued by a loss of its quantum nondemolition character and a decrease in fidelity with increased measurement strength. In this work, we elucidate the nature of this dynamical process, which we refer to as transmon ionization. We develop a comprehensive framework which provides a physical picture of the origin of transmon ionization. This framework consists of three complementary levels of descriptions: a fully quantized transmon-resonator model, a semiclassical model where the resonator is treated as a classical drive on the transmon, and a fully classical model. Crucially, all three approaches preserve the full cosine potential of the transmon and lead to similar predictions. This framework identifies the multiphoton resonances responsible for transmon ionization. It also allows one to efficiently compute numerical estimates of the photon number threshold for ionization, which are in remarkable agreement with recent experimental results. The tools developed within this work are both conceptually and computationally simple, and we expect them to become an integral part of the theoretical underpinning of all circuit QED experiments.

The performance of superconducting quantum circuits is primarily limited by dielectric loss due to interactions with two-level systems (TLSs). State-of-the-art circuits with engineered material interfaces are approaching a limit where dielectric loss from bulk substrates plays an important role. However, a microscopic understanding of dielectric loss in crystalline substrates is still lacking. In this work, we show that boron acceptors in silicon constitute a TLS bath that leads to an energy dissipation channel for superconducting circuits. We discuss how the electronic structure of boron acceptors leads to an effective TLS response in silicon. We sweep the boron concentration in silicon and demonstrate the bulk dielectric loss limit from boron acceptors. We show that boron-induced dielectric loss can be reduced in a magnetic field due to the spin-orbit structure of boron. This work provides the first detailed microscopic description of a TLS bath for superconducting circuits and demonstrates the need for ultrahigh-purity substrates for next-generation superconducting quantum processors.

We introduce a class of quantum non-Markovian processes—dubbed process trees—that exhibit polynomially decaying temporal correlations and memory distributed across timescales. This class of processes is described by a tensor network with treelike geometry whose component tensors are (1) causality-preserving maps (superprocesses) and (2) locality-preserving temporal change-of-scale transformations. We show that the long-range correlations in this class of processes tends to originate almost entirely from memory effects and can accommodate genuinely quantum power-law correlations in time. Importantly, this class allows efficient computation of multitime correlation functions. To showcase the potential utility of this model-agnostic class for numerical simulation of physical models, we show how it can efficiently approximate the strong memory dynamics of the paradigmatic spin-boson model, in terms of arbitrary multitime features. In contrast to an equivalently costly matrix-product-operator representation, the ansatz produces a fiducial characterization of the relevant physics. Finally, leveraging 2D tensor-network renormalization-group methods, we detail an algorithm for deriving a process tree from an underlying Hamiltonian via the Feynmann-Vernon influence functional. Our work lays the foundation for the development of more efficient numerical techniques in the field of strongly interacting open quantum systems, as well as the theoretical development of a temporal renormalization-group scheme.

Encoding quantum information in quantum states with disjoint wave-function support and noise-insensitive energies is the key behind the idea of qubit protection. While fully protected qubits are expected to offer exponential protection against both energy relaxation and pure dephasing, simpler circuits may grant partial protection with currently achievable parameters. Here, we study a fluxonium circuit in which the wave functions are engineered to minimize their overlap while benefiting from a first-order-insensitive flux sweet spot. Taking advantage of a large superinductance ($L\sim 1\mathrm{\mu }\mathrm{H}$), our circuit incorporates a resonant tunneling mechanism at zero external flux that couples states with the same fluxon parity, thus enabling bifluxon tunneling. The states $|0⟩$ and $|1⟩$ are encoded in wave functions with parities 0 and 1, respectively, ensuring a minimal form of protection against relaxation. Two-tone spectroscopy reveals the energy-level structure of the circuit and the presence of $4\pi$ quantum-phase slips between different potential wells corresponding to $m=\pm 1$ fluxons, which can be precisely described by a simple fluxonium Hamiltonian or by an effective bifluxon Hamiltonian. Despite suboptimal fabrication, the measured relaxation ($T_1=177\pm 3\mathrm{\mu }\mathrm{s}$) and dephasing ($T_2^{\mathrm{E}}=75\pm 5\mathrm{\mu }\mathrm{s}$) times not only demonstrate the relevance of our approach but also open an alternative direction toward quantum computing using partially protected fluxonium qubits.

During a quantum measurement, superpositions of states with different observable properties probabilistically collapse into one with a sharp value of the measured observable. In macroscopic quantum systems, this collapse arises via a continuous measurement-induced phase transition (MIPT) at a critical value of the strength of interaction with the measurement apparatus. MIPTs lie outside established paradigms for equilibrium or nonequilibrium critical phenomena and delineate distinct, stable dynamical and computational phases of matter. Quantum computers enable programmable simulation of the interaction of a measurement apparatus with a dynamical quantum system, to explore MIPT phenomena over a range of system sizes while retaining quantum coherence. Yet, existing experimental protocols rely on fundamentally nonscalable postselection techniques or direct classical simulation of quantum circuits. Here, we report the scalable observation of finite-size scaling evidence for an observable-sharpening MIPT in monitored quantum circuits in a chain of ${}^{171}{\mathrm{Yb}}^{+}$ ions in Quantinuum’s H1-1 trapped-ion quantum processor. By leveraging an equivalent description as a statistical physics problem, we implement scalable classical algorithms to infer the value of the measured observable from a single experimental shot. This technique enables a truly scalable protocol to observe observable-sharpening MIPTs in generic classes of circuits that cannot be directly classically simulated and also provides enhanced means to detect and suppress errors in the quantum simulation.

The negatively charged tin-vacancy center in diamond (${\mathrm{SnV}}^{-}$) is an emerging platform for building the next generation of long-distance quantum networks. This is due to the ${\mathrm{SnV}}^{-}$’s favorable optical and spin properties including bright emission, insensitivity to electronic noise, and long spin coherence times at temperatures above 1 K. Here, we demonstrate measurement of a single ${\mathrm{SnV}}^{-}$ electronic spin with a single-shot readout fidelity of 87.4%, which can be further improved to 98.5% by conditioning on multiple readouts. In the process, we develop understanding of the relationship between strain, magnetic field, spin readout, and microwave spin control. We show that high-fidelity readout is compatible with rapid microwave spin control, demonstrating a favorable parameter regime for use of the ${\mathrm{SnV}}^{-}$ center as a high-quality spin-photon interface. Finally, we use weak quantum measurement to study measurement-induced dephasing; this illuminates the fundamental interplay between measurement and decoherence in quantum mechanics, and provides a universal method to characterize the efficiency of color-center spin readout. Taken together, these results overcome an important hurdle in the development of the ${\mathrm{SnV}}^{-}$-based quantum technologies and, in the process, develop techniques and understanding broadly applicable to the study of solid-state quantum emitters.

To control and measure the state of a quantum system, it must necessarily be coupled to external degrees of freedom. This inevitably leads to spontaneous emission via the Purcell effect, photon-induced dephasing from measurement backaction, and errors caused by unwanted interactions with nearby quantum systems. To tackle this fundamental challenge, we make use of the design flexibility of superconducting quantum circuits to form a multimode element—an artificial molecule—with symmetry-protected modes. The proposed circuit consists of three superconducting islands coupled to a central island via Josephson junctions. It exhibits two essential nonlinear modes, one of which is flux insensitive and used as the protected qubit mode. The second mode is flux tunable and serves via a cross-Kerr-type coupling as a mediator to control the dispersive coupling of the qubit mode to the readout resonator. We demonstrate the Purcell protection of the qubit mode by measuring relaxation times that are independent of the mediated dispersive coupling. We show that the coherence of the qubit is not limited by photon-induced dephasing when detuning the mediator mode from the readout resonator and thereby reducing the dispersive coupling. The resulting highly protected qubit, which we refer to as P-mon, with tunable interactions may serve as a basic building block of a scalable quantum processor architecture, in which qubit decoherence is strongly suppressed.

Entanglement is an extraordinary feature of quantum mechanics. Sources of entangled optical photons were essential to test the foundations of quantum physics through violations of Bell’s inequalities. More recently, entangled many-body states have been realized via strong nonlinear interactions in microwave circuits with superconducting qubits. Here, we demonstrate a chip-scale source of entangled optical and microwave photonic qubits. Our device platform integrates a piezo-optomechanical transducer with a superconducting resonator which is robust under optical illumination. We drive a photon-pair generation process and employ a dual-rail encoding intrinsic to our system to prepare entangled states of microwave and optical photons. We place a lower bound on the fidelity of the entangled state by measuring microwave and optical photons in two orthogonal bases. This entanglement source can directly interface telecom wavelength time-bin qubits and gigahertz frequency superconducting qubits, two well-established platforms for quantum communication and computation, respectively.

We study the statistics of matrix elements of local operators in the basis of energy eigenstates in a paradigmatic, integrable, many-particle quantum theory, the Lieb-Liniger model of bosons with repulsive delta-function interactions. Using methods of quantum integrability, we determine the scaling of matrix elements with system size. As a consequence of the extensive number of conservation laws, the structure of matrix elements is fundamentally different from, and much more intricate than, the predictions of the eigenstate thermalization hypothesis for generic models. We uncover an interesting connection between this structure for local operators in interacting integrable models and the one for local operators that are not local with respect to the elementary excitations in free theories. We find that typical off-diagonal matrix elements $⟨\mathit{\mu }|\mathcal{O}|\mathit{\lambda }⟩$ in the same macrostate scale as $\exp \mathbf{(}-c^{\mathcal{O}}L\ln (L)-L{M}_{\mathit{\mu },\mathit{\lambda }}^{\mathcal{O}}\mathbf{)}$, where the probability distribution function for $M_{\mathit{\mu },\mathit{\lambda }}^{\mathcal{O}}$ is well described by Fréchet distributions and $c^{\mathcal{O}}$ depends only on macrostate information. In contrast, typical off-diagonal matrix elements between two different macrostates scale as $\exp (-d^{\mathcal{O}}{L}^2)$, where $d^{\mathcal{O}}$ depends only on macrostate information. Diagonal matrix elements depend only on macrostate information up to finite-size corrections.

In this paper, we develop a novel method to solve problems involving quantum optical systems coupled to coherent quantum feedback loops featuring time delays. Our method is based on exact mappings of such non-Markovian problems to equivalent Markovian driven dissipative quantum many-body problems. In this work, we show that the resulting Markovian quantum many-body problems can be solved (numerically) exactly and efficiently using tensor network methods for a series of paradigmatic examples, consisting of driven quantum systems coupled to waveguides at several distant points. In particular, we show that our method allows solving problems in so far inaccessible regimes, including problems with arbitrary long time delays and arbitrary numbers of excitations in the delay lines. We obtain solutions for the full real-time dynamics as well as the steady state in all these regimes. Finally, motivated by our results, we develop a novel mean-field approach, which allows us to find the solution semianalytically, and we identify parameter regimes where this approximation is in excellent agreement with our tensor network results.

By applying a microwave drive to a specially designed Josephson circuit, we have realized a double-well model system: a Kerr oscillator submitted to a squeezing force. We have observed, for the first time, the spectroscopic fingerprint of a quantum double-well Hamiltonian when its barrier height is increased: a pairwise level kissing (coalescence) corresponding to the exponential reduction of tunnel splitting in the excited states as they sink under the barrier. The discrete levels in the wells also manifest themselves in the activation time across the barrier which, instead of increasing smoothly as a function of the barrier height, presents steps each time a pair of excited states is captured by the wells. This experiment illustrates the quantum regime of Arrhenius’s law, whose observation is made possible here by the unprecedented combination of low dissipation, time-resolved state control, 98.5% quantum nondemolition single shot measurement fidelity, and complete microwave control over all Hamiltonian parameters in the quantum regime. Direct applications to quantum computation and simulation are discussed.

Group-IV color centers in diamond are promising candidates for quantum networks due to their dominant zero-phonon line and symmetry-protected optical transitions that connect to coherent spin levels. The negatively charged tin-vacancy (SnV) center possesses long electron spin lifetimes due to its large spin-orbit splitting. However, the magnetic dipole transitions required for microwave spin control are suppressed, and strain is necessary to enable these transitions. Recent work has shown spin control of strained emitters using microwave lines that suffer from Ohmic losses, restricting coherence through heating. We utilize a superconducting coplanar waveguide to measure SnV centers subjected to strain, observing substantial improvement. A detailed analysis of the SnV center electron spin Hamiltonian based on the angle-dependent splitting of the ground and excited states is performed. We demonstrate coherent spin manipulation and obtain a Hahn echo coherence time of up to $T_2=430\mathrm{\mu }\mathrm{s}$. With dynamical decoupling, we can prolong coherence to $T_2=10\mathrm{ms}$, about a sixfold improvement compared to earlier works. We also observe a nearby coupling ${}^{13}\mathrm{C}$ spin, which may serve as a quantum memory, thus substantiating the potential of SnV centers in diamond and demonstrates the benefit of superconducting microwave structures.

We consider the statistical properties of eigenstates of the time-evolution operator in chaotic many-body quantum systems. Our focus is on correlations between eigenstates that are specific to spatially extended systems and that characterize entanglement dynamics and operator spreading. In order to isolate these aspects of dynamics from those arising as a result of local conservation laws, we consider Floquet systems in which there are no conserved densities. The correlations associated with scrambling of quantum information lie outside the standard framework established by the eigenstate thermalization hypothesis (ETH). In particular, ETH provides a statistical description of matrix elements of local operators between pairs of eigenstates, whereas the aspects of dynamics we are concerned with arise from correlations among sets of four or more eigenstates. We establish the simplest correlation function that captures these correlations and discuss features of its behavior that are expected to be universal at long distances and low energies. We also propose a maximum-entropy ansatz for the joint distribution of a small number $n$ of eigenstates. In the case $n=2$, this ansatz reproduces ETH. For $n=4$ it captures both the growth with time of entanglement between subsystems, as characterized by the purity of the time-evolution operator, and also operator spreading, as characterized by the behavior of the out-of-time-order correlator. We test these ideas by comparing results from Monte Carlo sampling of our ansatz with exact diagonalization studies of Floquet quantum circuits.

Dynamical decoupling and Hamiltonian engineering are well-established techniques that have been used to control qubit systems. However, designing the corresponding methods for qudit systems has been challenging due to the lack of a Bloch sphere representation, more complex interactions, and additional control constraints. By identifying several general structures associated with such problems, we develop a formalism for the robust dynamical decoupling and Hamiltonian engineering of strongly interacting qudit systems. Our formalism significantly simplifies qudit pulse-sequence design while naturally incorporating robustness conditions necessary for experimental practicality. We experimentally demonstrate these techniques in a strongly interacting, disordered ensemble of spin-1 nitrogen-vacancy centers, achieving more than an order-of-magnitude improvement in coherence time over existing pulse sequences. We further describe how our techniques enable the engineering of exotic many-body phenomena such as quantum many-body scars, and open up new opportunities for quantum metrology with enhanced sensitivities. These results enable wide-reaching new applications for dynamical decoupling and Hamiltonian engineering in many-body physics and quantum metrology.

Fault-tolerant quantum computation with bosonic qubits often necessitates the use of noisy discrete-variable ancillae. In this work, we establish a comprehensive and practical fault-tolerance framework for such a hybrid system and synthesize it with fault-tolerant protocols by combining bosonic quantum error correction (QEC) and advanced quantum control techniques. We introduce essential building blocks of error-corrected gadgets by leveraging ancilla-assisted bosonic operations using a generalized variant of path-independent quantum control. Using these building blocks, we construct a universal set of error-corrected gadgets that tolerate a single-photon loss and an arbitrary ancilla fault for four-legged cat qubits. Notably, our construction requires only dispersive coupling between bosonic modes and ancillae, as well as beam-splitter coupling between bosonic modes, both of which have been experimentally demonstrated with strong strengths and high accuracy. Moreover, each error-corrected bosonic qubit is comprised of only a single bosonic mode and a three-level ancilla, featuring the hardware efficiency of bosonic QEC in the full fault-tolerant setting. We numerically demonstrate the feasibility of our schemes using current experimental parameters in the circuit-QED platform. Finally, we present a hardware-efficient architecture for fault-tolerant quantum computing by concatenating the four-legged cat qubits with an outer qubit code utilizing only beam-splitter couplings. Our estimates suggest that the overall noise threshold can be reached using existing hardware. These developed fault-tolerant schemes extend beyond their applicability to four-legged cat qubits and can be adapted for other rotation-symmetrical codes, offering a promising avenue toward scalable and robust quantum computation with bosonic qubits.

A relevant problem in the theory of open quantum systems is the lack of complete positivity of dynamical maps obtained after weak-coupling approximations, a famous example being the Redfield master equation. A number of approaches exist to recover well-defined evolutions under additional Markovian assumptions, but much less is known beyond this regime. Here, we propose a numerical method to cure the complete-positivity violation issue while preserving the non-Markovian features of an arbitrary original dynamical map. The idea is to replace its unphysical Choi operator with its closest physical one, mimicking recent work on quantum process tomography. We also show that the regularized dynamics is more accurate in terms of reproducing the exact dynamics, which allows us to heuristically push the utilization of these master equations in moderate coupling regimes, where the loss of positivity can have a relevant impact.

A ubiquitous problem in quantum physics is to understand the ground-state properties of many-body systems. Confronted with the fact that exact diagonalization quickly becomes impossible when increasing the system size, variational approaches are typically employed as a scalable alternative: Energy is minimized over a subset of all possible states and then different physical quantities are computed over the solution state. Despite remarkable success, rigorously speaking, all that variational methods offer are upper bounds on the ground-state energy. On the other hand, so-called relaxations of the ground-state problem based on semidefinite programming represent a complementary approach, providing lower bounds to the ground-state energy. However, in their current implementation, neither variational nor relaxation methods offer provable bound on other observables in the ground state beyond the energy. In this work, we show that the combination of the two classes of approaches can be used to derive certifiable bounds on the value of any observable in the ground state, such as correlation functions of arbitrary order, structure factors, or order parameters. We illustrate the power of this approach in paradigmatic examples of 1D and 2D spin-$1/2$ Heisenberg models. To improve the scalability of the method, we exploit the symmetries and sparsity of the considered systems to reach sizes of hundreds of particles at much higher precision than previous works. Our analysis therefore shows how to obtain certifiable bounds on many-body ground-state properties beyond energy in a scalable way.

Interfacing cold atoms with integrated nanophotonic devices could offer new paradigms for engineering atom-light interactions and provide a potentially scalable route for quantum sensing, metrology, and quantum information processing. However, it remains a challenging task to efficiently trap a large ensemble of cold atoms on an integrated nanophotonic circuit. Here, we demonstrate direct loading of an ensemble of up to 70 atoms into an optical microtrap on a nanophotonic microring circuit. Efficient trap loading is achieved by employing degenerate Raman-sideband cooling in the microtrap, where a built-in spin-motion coupling arises directly from the vector light shift of the evanescent-field potential on a microring. Atoms are cooled into the trap via optical pumping with a single free space beam. We have achieved a trap lifetime approaching 700 ms under continuous cooling. We show that the trapped atoms display large cooperative coupling and superradiant decay into a whispering-gallery mode of the microring resonator, holding promise for explorations of new collective effects. Our technique can be extended to trapping a large ensemble of cold atoms on nanophotonic circuits for various quantum applications.

Ultracold polar molecules are promising for quantum information processing and searches for physics beyond the standard model. Laser cooling to ultracold temperatures is an established technique for trapped diatomic and triatomic molecules. Further cooling of the molecules to near the motional ground state is crucial for reducing various dephasings in quantum and precision applications. In this work, we demonstrate Raman sideband cooling (RSC) of CaF molecules in optical tweezers to near their motional ground state, with average motional occupation quantum numbers of ${\overline{n}}_x=0.16(12)$, ${\overline{n}}_y=0.17(17)$ (radial directions), and ${\overline{n}}_z=0.22(16)$ (axial direction), and a 3-D motional-ground-state probability of $54\pm 18\%$ of the molecules that survive the RSC. This process paves the way to increase molecular coherence times in optical tweezers for robust quantum computation and simulation applications.