Dedicated to the centenary of quantum mechanics. In this paper, we introduce a novel and general framework for the variational quantum simulation of Lindblad equations. Building on the close relationship between the unraveled Lindblad dynamics, stochastic Magnus integrators, and variational quantum simulation, we propose a high-order scheme for solving the quantum state diffusion equation using exponential integrators. This formulation facilitates the simulation of wavefunction trajectories within the established framework of variational quantum algorithms for time evolution. Our algorithm significantly enhances robustness in two key aspects: the stability of the simulation with large time steps, and the reduction in the number of quantum trajectories required to accurately simulate the Lindblad dynamics in terms of the ensemble average. We demonstrate the effectiveness of our algorithm through numerical examples in both classical and quantum implementations, including the transverse-field Ising model (TFIM) with damping, the Fenna-Matthews-Olson (FMO) complex, and the radical pair model (RPM). The simulation accuracy can be systematically improved, and the algorithm remains reliable even in highly oscillatory regimes. These methods are expected to be applicable to a broader class of open quantum systems beyond the specific models considered in this study.

Symmetry breaking has been a central theme in classifying quantum phases and phase transitions. Recently, this concept has been extended to the mixed states of open systems, attracting considerable attention due to the emergence of novel physics beyond closed systems. In this work, we reveal a new type of phase transition in mixed states, termed quantum strong-to-weak spontaneous symmetry breaking (SWSSB). Using a combination of field theory calculations and large-scale matrix product state simulations, we map out the global phase diagram of the $XXZ$ critical spin chain under local strong symmetry preserving decoherence, which features an SWSSB phase and a trivial Luttinger liquid phase, separated by a straight critical line that belongs to the boundary Berezinskii-Kosterlitz-Thouless universality class with a varying effective central charge. Importantly, we analyze this transition from two complementary perspectives: on one hand, through the behavior of order parameters that characterize the symmetry breaking; on the other hand, from a quantum information viewpoint by studying entropic quantities and the concept of quantum recoverability. This dual approach allows us to provide a more comprehensive understanding of the phase structure and the nature of the transition. Remarkably, the SWSSB transition in our case is purely quantum in the sense that it can only be driven by tuning the Hamiltonian parameter even under arbitrary decoherence strength, fundamentally distinguishing it from the decoherence-driven SWSSB transitions extensively discussed in previous literature. Importantly, our unified theoretical framework is applicable to a broad class of one-dimensional quantum systems, including spin chains and fermionic systems, whose low-energy physics can be described by Luttinger liquid theory, under arbitrary symmetry-preserving decoherence channels. Finally, we also discuss the experimental relevance of our theory on quantum simulator platforms.

Symmetry inherent in quantum states has been widely used to reduce the effect of noise in quantum error correction and a quantum error-mitigation technique known as symmetry verification. However, these symmetry-based techniques exploit symmetry in quantum states rather than quantum channels, limiting their application to cases where the entire circuit shares the same symmetry. In this work, we propose symmetric channel verification (SCV), a channel purification protocol that leverages the symmetry inherent in quantum channels. By introducing different phases to each symmetric subspace and employing a quantum phase estimation-like circuit, SCV can detect and correct symmetry-breaking noise in quantum channels. We further propose a hardware-efficient implementation of SCV at the virtual level, which requires only a single-qubit ancilla and is robust against the noise in the ancilla qubit. Our protocol is applied to various Hamiltonian simulation circuits and phase estimation circuits, resulting in a significant reduction of errors. Furthermore, in setups where only Clifford unitaries can be used for noise purification, which is relevant in the early fault-tolerant regime, we show that SCV under Pauli symmetry represents the optimal purification method.

Quantum networking technologies use spin qubits and their interface to single photons as core components of a network node. This necessitates the ability to co-design the magnetic- and optical-dipole response of a quantum system—a capability that has been notably absent in solid-state platforms where spin-orbit coupling and the crystalline environment lead to inhomogeneity of electronic $g$-factors and optically active states. Here, we demonstrate the ability to design both the optical and magnetic response of a solid-state quantum emitter a priori. We show that $\mathrm{Ga}\mathrm{As}$ quantum dots (QDs), obtained via local droplet etching epitaxy and already known as exceptionally coherent and efficient quantum light sources, also exhibit spin and optical properties that follow directly from assuming the highest possible system symmetry. Our measurements of electron and hole $g$-tensors—using a new sign-sensitive measurement protocol based on the hyperfine interaction—and of transition dipole moment orientations for charged excitons agree with our predictions from a multiband $\mathit{k}\cdot \mathit{p}$ simulation constrained only by a single atomic-force-microscopy reconstruction of QD morphology. This agreement is verified across multiple wavelength-specific growth runs at different facilities within the range of 730 to 790 nm for the exciton emission. Remarkably, our measurements and simulations track the in-plane electron $g$-factors through a zero-crossing from $-0.1$ to 0.3 and linear optical dipole moment orientations fully determined by an external magnetic field. The robustness and generality of these results establish a fundamentally new paradigm for solid-state spin-photon interfaces: one in which the properties of a spin qubit and its tunable optical interface can be designed—prior to growth—for a target magnetic and photonic environment, with direct applications to scalable and high-fidelity spin-photon entanglement.

Motivated by the applications of secure multi-party computation as a privacy-protecting data analysis tool, and identifying oblivious transfer as one of its main practical enablers, we propose a practical realization of randomized quantum oblivious transfer. By using only symmetric cryptography primitives to implement commitments, we construct computationally secure randomized oblivious transfer without the need for public-key cryptography or assumptions imposing limitations on the adversarial devices. We show that the protocol is secure under an indistinguishability-based notion of security and demonstrate an experimental implementation to test its real-world performance. Its security and performance are then compared to both quantum and classical alternatives, showing potential advantages over existing solutions based on the noisy storage model and public-key cryptography.

We investigate the finite-temperature phase diagram of polar molecules shielded by dual microwave fields using the path integral Monte Carlo method combined with the worm algorithm. We determine the critical temperature $T_c$ for Bose-Einstein condensations (BECs) and identify two distinct phases below $T_c$: the expanding gas (EG) phase and the self-bound gas (SBG) phase. We further analyze the temperature and interaction-strength dependence of the condensate and superfluid fractions. Notably, in contrast to dilute atomic BECs, the SBG phase displays a low condensate fraction and a high superfluid fraction, resembling the behavior of strongly correlated ${}^4\mathrm{He}$ superfluids. These significant many-body correlations arise from the interplay between long-range dipole-dipole interactions and the short-range shielding potential. Furthermore, we demonstrate that the aspect ratio of the gas provides a characteristic geometric signature to accurately determine the EG-to-SBG transition, robust against external trapping potentials. Our findings provide unbiased and numerically exact results to guide upcoming experiments with polar molecules.

When experimentally learning the action of a continuous-variable quantum process by probing it with inputs, there will often be some restriction on the input states used. One experimentally simple way to probe a quantum channel is to use low-energy coherent states. Learning a quantum channel in this way presents difficulties, due to the fact that two channels may act similarly on low-energy inputs but very differently for high-energy inputs. They may also act similarly on coherent-state inputs but differently on nonclassical inputs. Extrapolating the behavior of a channel for more general input states from its action on the far more limited set of low-energy coherent states is a case of out-of-distribution generalization. To be sure that such generalization gives meaningful results, one needs to relate error bounds for the training set to bounds that are valid for all inputs. We show that for any pair of channels that act sufficiently similarly on low-energy coherent-state inputs, one can bound how different the input-output relations are for any (high-energy or highly nonclassical) input. This proves that out-of-distribution generalization is always possible for learning quantum channels using low-energy coherent states, as long as enough samples are used.

The surface code, one of the leading candidates for quantum error correction, is known to protect encoded quantum information against stochastic, i.e., incoherent errors. The protection against coherent errors, such as from unwanted gate rotations, is however understood only for special cases, such as rotations about the $X$ or $Z$ axes. Here we consider generic single-qubit coherent errors in the surface code, i.e., rotations by angle $\alpha$ about an axis that can be chosen arbitrarily. We develop a statistical mechanical mapping for such errors and perform entanglement analysis in transfer matrix space to numerically establish the existence of an error-correcting phase, which we chart in a subspace of rotation axes to estimate the corresponding maximum-likelihood thresholds ${\alpha }_{\mathrm{th}}$. The classical statistical mechanics model we derive is a random-bond Ising model with complex couplings and four-spin interactions (i.e., a complex-coupled Ashkin-Teller model). The error-correcting phase, $\alpha <{\alpha }_{\mathrm{th}}$, where the logical error rate decreases exponentially with code distance, is shown to correspond in transfer matrix space to a gapped one-dimensional quantum Hamiltonian exhibiting spontaneous breaking of a ${\mathbb{Z}}_2$ symmetry. Our numerical results rest on two key ingredients: (i) we show that the state evolution under the transfer matrix, a nonunitary $(1+1)$-dimensional quantum circuit, can be efficiently numerically simulated using matrix product states; and (ii) based on this approach, we also develop an algorithm to (approximately) sample syndromes based on their Born probability. The ${\alpha }_{\mathrm{th}}$ values we find show that the maximum-likelihood thresholds for coherent errors are larger than those for the corresponding incoherent errors (from the Pauli twirl), and significantly exceed the values found using minimum weight perfect matching.

Dispersive readout of superconducting qubits is often limited by readout-drive-induced transitions between qubit levels. While there is a growing understanding of such effects in transmon qubits, the case of highly nonlinear fluxonium qubits is more complex. We theoretically analyze measurement-induced state transitions (MISTs) during the dispersive readout of a fluxonium qubit. We focus on a new mechanism: a simultaneous transition/excitation involving the qubit and an internal mode of the Josephson-junction array in the fluxonium circuit. Using an adiabatic Floquet approach, we show that these new kinds of MIST processes can be relevant when using realistic circuit parameters and relatively low readout drive powers. They also contribute to excess qubit dephasing even after a measurement is complete. In addition to outlining basic mechanisms, we also investigate the dependence of such transitions on the circuit parameters. We find that with a judicious choice of frequency allocations or coupling strengths, these parasitic processes can most likely be avoided.

Learning faithful representations of quantum states is crucial to fully characterizing the variety of many-body states created on quantum processors. While various tomographic methods, such as classical shadow and matrix product state (MPS) tomography have shown promise in characterizing a wide class of quantum states, they face unique limitations in detecting topologically ordered two-dimensional states. To address this problem, we implement and study a heuristic tomographic method that combines variational optimization on tensor networks with randomized measurement techniques. Using this approach, we demonstrate its ability to learn the ground state of the surface-code Hamiltonian as well as an experimentally realizable quantum spin liquid state. In particular, we perform numerical experiments using MPS ansätze and systematically investigate the sample complexity required to achieve high fidelities for systems with sizes of up to 48 qubits. In addition, we provide theoretical insights into the scaling of our learning algorithm by analyzing the statistical properties of maximum-likelihood estimation. Notably, our method is sample-efficient and experimentally friendly, only requiring snapshots of the quantum state measured randomly in the $X$ or $Z$ bases. Using this subset of measurements, our approach can effectively learn any real pure states represented by tensor networks, and we rigorously prove that random-$XZ$ measurements are tomographically complete for such states.

One of the critical challenges solid-state quantum processors face is the presence of fabrication imperfections and two-level systems, which render certain qubits and gates either inoperable or much noisier than tolerable by quantum error-correction protocols. To address this challenge, we develop a suite of novel and highly performant methods for adapting surface code patches in the presence of defective qubits and gates, which we call Snakes and Ladders. We explain how our algorithm generates and compares several strategies in order to find the optimal one for any given configuration of defective components, as well as introduce heuristics to improve runtime and minimize computing resources required by our algorithm. In addition to memory storage we also show how to apply our methods to lattice surgery protocols. Compared to prior works, our methods significantly improve the code distance of the adapted surface code patches for realistic defect rates, resulting in a logical performance similar to that of the defect-free patches.

The interplay between coherence and system-environment interactions is at the basis of a wide range of phenomena, from quantum information processing to charge and energy transfer in molecular systems, biomolecules, and photochemical materials. In this work, we use a Frenkel exciton model with long-range interacting qubits coupled to a damped collective bosonic mode to investigate vibrationally assisted transfer processes in donor-acceptor systems featuring internal substructures analogous to light-harvesting complexes. We find that certain delocalized excitonic states maximize the transfer rate and that the entanglement is preserved during the dissipative transfer over a wide range of parameters. We investigate the reduction in transfer caused by static disorder, white noise, and finite temperature and study how transfer efficiency scales as a function of the number of dimerized monomers and the component number of each monomer, finding which excitonic states lead to optimal transfer. Finally, we provide a realistic experimental setting to realize this model in analog trapped-ion quantum simulators. Analog quantum simulation of systems comprising many and increasingly complex monomers could offer valuable insights into the design of light-harvesting materials, particularly in the nonperturbative intermediate parameter regime examined in this study, where classical simulation methods are resource intensive.

Classical shadow estimation (CSE) is a powerful tool for learning the properties of quantum states and quantum processes. Here we consider the CSE task for quantum unitary channels. By querying an unknown unitary channel $\mathcal{U}$ multiple times in quantum experiments, the goal is to learn a classical description from which one can accurately predict many different linear properties of the channel, i.e., the expectation values of arbitrary observables measured on the output of $\mathcal{U}$ upon arbitrary input states. Based on collective measurements on multiple systems, we propose a query-efficient protocol for this task, whose query complexity has a quadratic advantage over the previous best approach for this problem, and almost saturates the information-theoretic lower bound. To further enhance practicality, we also present a variant protocol using only single-copy measurements, which still offers much better query performance than previous protocols that do not use quantum memory, and can serve as a key subroutine for learning an arbitrary unknown Hamiltonian from dynamics. In addition to linear properties of unitary channels, our protocol can also be applied to simultaneously predict many nonlinear properties, such as out-of-time-ordered correlators.

The Majorana tetron is a prototypical topological qubit stemming from the ground-state degeneracy of a superconducting island hosting four Majorana modes. This degeneracy manifests as an effective nonlocal spin degree of freedom, the most paradigmatic signature of which is the topological Kondo effect. Degeneracies of states with different fermionic parities also characterize minimal Kitaev chains, which have lately emerged as a platform to realize and study unprotected versions of Majorana modes, dubbed poor man’s Majorana modes. Here, we introduce the “poor man’s Majorana tetron,” comprising four quantum dots coupled via a floating superconducting island. Its charging energy yields nontrivial correlations among the dots, although, unlike a standard tetron, it is not directly determined by the fermionic parity of the Majorana modes. The poor man’s tetron displays parameter regions with a twofold-degenerate ground state with odd fermionic parity, which gives rise to an effective Anderson impurity model when coupled to external leads. We show that this system can approach a regime featuring the topological Kondo effect under a suitable tuning of experimental parameters. Therefore, the poor man’s tetron is a promising device to observe the nonlocality of Majorana modes and their related fractional conductance.

Ordered atomic arrays feature an enhanced collective optical response compared to random atomic ensembles due to constructive interference in resonant dipole-dipole interactions. One consequence of this is the existence of a large shift of the transition with respect to the bare atomic frequency. In the linear optics regime (low light intensity), one observes a spectroscopic shift of the Lorentzian atomic line often called the collective Lamb shift. For stronger driving, many excitations are present in the system rendering the calculation of this shift theoretically challenging, but its understanding is important for instance when performing Ramsey spectroscopy in optical clocks. Here we report on the study of the collective optical response of a one-dimensional array of 30 dysprosium atoms. We drive the atoms on the narrow intercombination transition isolating a two-level system, and measure the atomic state with single-shot state readout using a broad transition. In the linear optics regime, we measure the shift of the resonance in steady state due to dipole interactions, and measure how this shift depends on the interatomic distance. We further resolve at the single atom level how the excitation is distributed over the array. Then, on the same transition we perform Ramsey spectroscopy i.e., away from the linear regime. We observe a time-dependent shift, that allows us to draw the connection between the collective Lamb shift observed in the linear optics regime and in the large-excitation case.

An important open question for the current generation of highly controllable quantum devices is understanding which phases can be realized as stable steady states under local quantum dynamics. In this work, we show how robust steady-state phases with both Abelian and non-Abelian mixed-state topological order can be stabilized, in two spatial dimensions, against generic “heralded” noise using active dynamics that incorporate measurement and feedback, modeled as a fully local Lindblad master equation. These topologically ordered steady states are two-way connected to pure topologically ordered ground states using local quantum channels, and preserve quantum information for a time that is exponentially large in the system size. Specifically, we present explicit constructions of families of local Lindbladians for both Abelian (${\mathbb{Z}}_2$) and non-Abelian ($D_4$) topological order whose steady states host mixed-state topological order when the noise is below a threshold strength. As the noise strength is increased, these models exhibit first-order transitions to intermediate mixed-state phases where they encode robust classical memories, followed by (first-order) transitions to a trivial steady state at high noise rates. When the noise is imperfectly heralded, steady-state order disappears but our active dynamics significantly enhances the lifetime of the encoded logical information. To carry out the numerical simulations for the non-Abelian $D_4$ case, we introduce a generalized stabilizer tableau formalism that permits efficient simulation of the non-Abelian Lindbladian dynamics.

Boson sampling stands out as a promising approach toward experimental demonstration of quantum computational advantage. However, the presence of physical noise in near-term experiments hinders the realization of quantum computational advantage with boson sampling. Since physical noise in near-term boson-sampling devices is inevitable, precise characterization of the boundary of noise rates where the classical intractability of boson sampling is maintained is crucial for quantum computational advantage using near-term devices. In this work, we identify the level of partial-distinguishability noise that upholds the classical intractability of boson sampling. We find that boson sampling with on average $O(\log N)$ distinguishable photons out of $N$ input photons maintains the equivalent complexity to the ideal boson-sampling case. By providing strong complexity-theoretical evidence for the classical intractability of noisy boson sampling, we expect that our findings will ultimately facilitate the demonstration of quantum computational advantage with noisy boson-sampling experiments in the near future.

A complex but important challenge in understanding quantum mechanical phenomena is the simulation of quantum many-body dynamics. Although quantum computers offer significant potential to accelerate these simulations, their practical application is currently limited by noise and restricted scalability. In this work, we address these problems by introducing a hybrid ansatz combining the strengths of quantum and classical computational methods. Using Trotterization, we propose to evolve an initial state on the quantum computer according to a simplified Hamiltonian, focusing on terms that are difficult to simulate classically. A classical model then corrects the simulation by including the terms omitted in the quantum circuit. While the classical ansatz is optimized during the time evolution, the quantum circuit has no variational parameters. Derivatives can thus be calculated purely classically, avoiding challenges arising in the optimization of parametrized quantum circuits. We demonstrate three applications of this hybrid method. First, our approach allows us to avoid swap gates in the quantum circuit by restricting the quantum part of the ansatz to hardware-efficient terms of the Hamiltonian. Second, we can mitigate errors arising from the Trotterization of the time-evolution unitary. Finally, we can extend the system size while keeping the number of qubits in the quantum circuit constant by including additional degrees of freedom in the classical ansatz.

Semiconductor-based spin qubits embedded into a superconducting microwave cavity constitute a fast-progressing and promising platform for realizing fast and fault-tolerant qubit control with long-range two-qubit coupling. The flopping-mode spin qubit consists of a single electron in a double quantum dot; it combines a charge qubit with a spin qubit. With its strong and tunable cavity coupling, the flopping-mode qubit is proven to be well suited for low-power qubit control and cavity-mediated long-range quantum gates. The singlet-triplet (ST) and exchange-only (EO) qubits are multielectron realizations that go without broadband control and are protected from some types of noise, but are challenging to couple to each other and to microwave cavities. We combine the flopping-mode concept with the ST and EO qubits and propose two new flopping-mode qubits that consist of three (four) quantum dots, occupied by two (three) electrons near the $(1,0,1)\leftrightarrow (0,1,1)[(1,0,1,1)\leftrightarrow (0,1,1,1)]$ charge transition. The two-electron system augments the ${\mathrm{ST}}_0$ spin qubit with a charge qubit that interacts transversally and longitudinally with a cavity. Both couplings are highly tunable, and the longitudinal coupling distinguishes the flopping-mode ST qubit from the regular flopping-mode qubit. The longitudinal coupling allows for nondissipative universal control similar to superconducting transmon qubits. The EO flopping-mode qubit comprises four dots occupied by three electrons and opens a new possibility to perform two-qubit gates for EO qubits that are challenging to perform directly with the exchange coupling. We use input-output theory to provide means of extracting the coupling strengths from cavity transmission data.

In classical thermodynamics, heat must spontaneously flow from hot to cold systems. In quantum thermodynamics, the same law applies when considering multipartite product thermal states evolving unitarily. If initial correlations are present, anomalous heat flow can happen, temporarily making cold thermal states colder and hot thermal states hotter. Such an effect can happen due to entanglement but also because of classical randomness, hence lacking a direct connection with nonclassicality. In this work, we introduce scenarios in which anomalous heat flow does have a direct link to nonclassicality, defined as the failure of noncontextual models to explain experimental data. We start by extending known noncontextuality inequalities to a setup in which sequential transformations are considered. We then show a class of quantum prepare-transform-measure protocols, characterized by a time interval $(0,{\tau }_c)$ for a given critical time ${\tau }_c$, where anomalous heat flow happens only if a noncontextuality inequality is violated. We also analyze a recent experiment from Micadei et al. [Nat. Commun. 10, 2456 (2019)] and find the critical time ${\tau }_c$ based on the authors’ experimental parameters. We conclude by investigating heat flow in the evolution of two-qudit systems, showing that our findings are not an artifact of using two-qubit systems.