Publications from 2025

We demonstrate an experimentally feasible circuit-QED Bose-Hubbard simulator that reproduces the complex spin dynamics of Heisenberg models. Our method relies on mapping spin-1/2 systems onto bosonic states via the polynomially expanded Holstein-Primakoff (HP) transformation. The HP transformation translates the intricate behavior of spins into a representation that is compatible with bosonic devices like those in a circuit QED setup. For comparison, we also implement the Dyson-Maleev (DM) encoding for spin-1/2 and show that, in this limit, DM and HP are equivalent. We show the equivalence of the DM and the HP transformations for spin-1/2 systems. Rigorous numerical analyses confirm the effectiveness of our HP-based protocol. Specifically, we obtain the concurrence between the spin dynamics and the behavior of microwave photons within our circuit QED-based analog simulator that is designed for the Bose-Hubbard model. By utilizing the microwave photons inherent to circuit QED devices, our framework presents an accessible, scalable avenue for probing quantum spin dynamics in an experimentally viable setting.

Photonic lattices facilitate band structure engineering, supporting both localized and extended modes through their geometric design. However, greater control over these modes can be achieved by taking advantage of the interference effect between external drives with precisely tuned phases and photonic modes within the lattice. In this work, we build on this principle to demonstrate optical switching, directed light propagation and site-specific localization in a one-dimensional photonic lattice of coupled microresonators by resonantly driving the system with a coherent field of controlled phase. Importantly, our experimental results provide direct evidence that increased driving power acts as a tuning parameter enabling nonlinear localization at frequencies previously inaccessible in the linear regime. These findings open new avenues for controlling light propagation and localization in lattices with more elaborate band structures.

Artificial atoms non-locally coupled to waveguides — the so-called giant atoms — offer new opportunities for the control of light and matter. In this work, we show how to use an array of non-locally coupled transmon “molecules” to engineer a passive photonic controlled gate for waveguide photons. In particular, we show that a conditional elastic phase shift between counter-propagating photons arises from the interplay between direction-dependent couplings, engineered through an interplay of non local interactions and molecular binding strength; and the nonlinearity of the transmon array. We analyze the conditions under which a maximal $\pi$-phase shift — and hence a CZ gate — is obtained, and characterize the gate fidelity as a function of key experimental parameters, including finite transmon nonlinearities, emitter spectral inhomogeneities, and limited cooperativity. Our work opens the use of giant atoms as key elements of microwave photonic quantum computing devices.

The field of superconducting qubits is constantly evolving with new circuit designs. However, when it comes to qubit readout, the use of simple transverse linear coupling remains overwhelmingly prevalent. This standard readout scheme has significant drawbacks: in addition to the Purcell effect, it suffers from a limitation on the maximal number of photons in the readout mode, which restricts the signal-to-noise ratio (SNR) and the Quantum Non-Demolition (QND) nature of the readout. Here, we explore the high-power regime by engineering a nonlinear coupling between a transmon qubit and its readout mode. Our approach builds upon previous work by Dassonneville et al. [Physical Review X 10, 011045 (2020)], on qubit readout with a non-perturbative cross-Kerr coupling in a transmon molecule. We demonstrate a readout fidelity of 99.21% with 89 photons utilizing a parametric amplifier. At this elevated photon number, the QND nature remains high at 96.7%. Even with up to 300 photons, the QNDness is only reduced by a few percent. This is qualitatively explained by deriving a critical number of photons associated to the nonlinear coupling, yielding a theoretical value of $\bar{n}_r^\text{crit} = 377$ photons for our sample’s parameters. These results highlight the promising performance of the transmon molecule in the high-power regime for high-fidelity qubit readout.

Long-range quantum systems, in which the interactions decay as $1/r^{\alpha}$, are of increasing interest due to the variety of experimental set-ups in which they naturally appear. Motivated by this, we study fundamental properties of long-range spin systems in thermal equilibrium, focusing on the weak regime of $ \alpha>D$. Our main result is a proof of analiticity of their partition functions at high temperatures, which allows us to construct a classical algorithm with sub-exponential runtime $\exp(\mathcal{O}(\log^2(N/\epsilon)))$ that approximates the log-partition function to small additive error $\epsilon$. As by-products, we establish the equivalence of ensembles and the Gaussianity of the density of states, which we verify numerically in both the weak and strong long-range regimes. This also yields constraints on the appearance of various classes of phase transitions, including thermal, dynamical and excited-state ones. Our main technical contribution is the extension to the quantum long-range regime of the convergence criterion for cluster expansions of Koteck\’y and Preiss.

Quantum Floquet engineering (QFE) seeks to generalize the control of quantum systems with classical external fields, widely known as Semi-Classical Floquet engineering (SCFE), to quantum fields. However, to faithfully capture the physics at arbitrary coupling, a gauge-invariant description of light-matter interaction in cavity-QED materials is required, which makes the Hamiltonian highly non-linear in photonic operators. We provide a non-perturbative truncation scheme of the Hamiltonian, which is valid or arbitrary coupling strength, and use it to investigate the role of light-matter correlations, which are absent in SCFE. We find that even in the high-frequency regime, light-matter correlations can be crucial, in particular for the topological properties of a system. As an example, we show that for a SSH chain coupled to a cavity, light-matter correlations break the original chiral symmetry of the chain, strongly affecting the robustness of its edge states. In addition, we show how light-matter correlations are imprinted in the photonic spectral function and discuss their relation with the topology of the bands.

We present a flip-chip architecture for an array of coupled superconducting qubits, in which circuit components reside inside individual microwave enclosures. In contrast to other flip-chip approaches, the qubit chips in our architecture are electrically floating, which guarantees a simple, fully modular assembly of capacitively coupled circuit components such as qubit, control, and coupling structures, as well as reduced crosstalk between the components. We validate the concept with a chain of three nearest neighbor coupled generalized flux qubits in which the center qubit acts as a frequency-tunable coupler. Using this coupler, we demonstrate a transverse coupling on/off ratio $\approx$ 50, zz-crosstalk $\approx$ 0.7 kHz between resonant qubits and isolation between the qubit enclosures > 60 dB.

Poor man’s Majorana bound states (MBS) arise in minimal Kitaev chains when the parameters are fine-tuned to a sweet spot. We consider an interacting two-site Kitaev chain coupled to a single-mode cavity and show that the sweet spot condition can be controlled with the cavity frequency and the hopping between sites. Furthermore, we demonstrate that photon-mediated effective interactions can be used to screen intrinsic interactions, improving the original quality of the MBS. We describe experimental signatures in the cavity transmission to detect their presence and quality. Our work proposes a new way to tune poor man’s MBS in a quantum dot array coupled to a cavity.

Sum-frequency generation (SFG) allows for coherent upconversion of an electromagnetic signal and has applications in mid-infrared vibrational spectroscopy of molecules. Recent experimental and theoretical studies have shown that plasmonic nanocavities, with their deep sub-wavelength mode volumes, may allow to obtain vibrational SFG signals from a single molecule. In this article, we compute the degree of second order coherence ($g^{(2)}(0)$) of the upconverted mid-infrared field under realistic parameters and accounting for the anharmonic potential that characterizes vibrational modes of individual molecules. On the one hand, we delineate the regime in which the device should operate in order to preserve the second-order coherence of the mid-infrared source, as required in quantum applications. On the other hand, we show that an anharmonic molecular potential can lead to antibunching of the upconverted photons under coherent, Poisson-distributed mid-infrared and visible drives. Our results therefore open a path toward a new kind of bright and tunable source of indistinguishable single photons by leveraging “vibrational blockade” in a resonantly and parametrically driven molecule, without the need for strong light-matter coupling.

Fermionic atoms in optical lattices provide a native implementation of Fermi-Hubbard (FH) models that can be used as analog quantum simulators of many-body fermionic systems. Recent experimental advances include the time-dependent local control of chemical potentials and tunnelings, and thus enable to operate this platform digitally as a programmable quantum simulator. Here, we explore these opportunities and develop ground-state preparation algorithms for different fermionic models, based on the ability to implement both single-particle and many-body, high-fidelity fermionic gates, as provided by the native FH Hamiltonian. In particular, we first design variational, pre-compiled quantum circuits to prepare the ground state of the natively implemented FH model, with significant speedups relative to competing adiabatic protocols. Besides, the versatility of this variational approach enables to target extended FH models, i.e., including terms that are not natively realized on the platform. As an illustration, we include next-nearest-neighbor tunnelings at finite dopings, relevant in the context of $d$-wave superconductivity. Furthermore, we discuss how to approximate the imaginary-time evolution using variational fermionic circuits, both as an alternative state-preparation strategy, and as a subroutine for the Quantum Lanczos algorithm to further improve the energy estimation. We benchmark our protocols for ladder geometries, though they can be readily applied to 2D experimental setups to address regimes beyond the capabilities of current classical methods. These results pave the way for more efficient and comprehensive explorations of relevant many-body phases with existing programmable fermionic quantum simulators.

We present a quantum-inspired solver for the one-dimensional Gross-Pitaevskii equation in the Quantics Tensor-Train (QTT) representation. By evolving the system entirely within a low-rank tensor manifold, the method sidesteps the memory and runtime barriers that limit conventional finite-difference and spectral schemes. Two complementary algorithms are developed: an imaginary-time projector that drives the condensate toward its variational ground state and a rank-adapted fourth-order Runge-Kutta integrator for real-time dynamics. The framework captures a broad range of physical scenarios – including barrier-confined condensates, quasi-random potentials, long-range dipolar interactions, and multicomponent spinor dynamics – without leaving the compressed representation. Relative to standard discretizations, the QTT approach achieves an exponential reduction in computational resources while retaining quantitative accuracy, thereby extending the practicable regime of Gross-Pitaevskii simulations on classical hardware. These results position tensor networks as a practical bridge between high-performance classical computing and prospective quantum hardware for the numerical treatment of nonlinear Schrodinger-type partial differential equations.

We report a systematic investigation of universal quantum chaotic signatures in the transverse field Ising model on an Erd\H{o}s-R\’enyi network. This is achieved by studying local spectral measures such as the level spacing and the level velocity statistics. A spectral form factor analysis is also performed as a global measure, probing energy level correlations at arbitrary spectral distances. Our findings show that these measures capture the breakdown of chaotic behavior upon varying the connectivity and strength of the transverse field in various regimes. We demonstrate that the level spacing statistics and the spectral form factor signal this breakdown for sparsely and densely connected networks. The velocity statistics capture the surviving chaotic signatures in the sparse limit. However, these integrable-like regimes extend over a vanishingly small segment in the full range of connectivity.

The last two decades has seen quantum thermodynamics become a well established field of research in its own right. In that time, it has demonstrated a remarkably broad applicability, ranging from providing foundational advances in the understanding of how thermodynamic principles apply at the nano-scale and in the presence of quantum coherence, to providing a guiding framework for the development of efficient quantum devices. Exquisite levels of control have allowed state-of-the-art experimental platforms to explore energetics and thermodynamics at the smallest scales which has in turn helped to drive theoretical advances. This Roadmap provides an overview of the recent developments across many of the field’s sub-disciplines, assessing the key challenges and future prospects, providing a guide for its near term progress.

The study of many-body quantum systems out of equilibrium remains a significant challenge with complexity barriers arising in both state and operator-based representations. In this work, we review recent approaches based on finding better contraction strategies for the full spatio-temporal tensor networks that encode the path integral of the dynamics, as well as the conceptual integration of influence functionals, process tensors, and transfer matrices within the tensor network formalism. We discuss recent algorithmic developments, highlight the complexity of influence functionals in various dynamical regimes and present consistent results of different communities, showing how ergodic dynamics render these functionals exponentially difficult to compress. Finally, we provide an outlook on strategies to encode complementary influence functional overlaps, paving the way for accurate descriptions of open and closed quantum systems with tensor networks.

Understanding the dynamical consequences of quantum phase transitions on thermodynamical quantities, such as work statistics and entropy production, is one of the most intriguing aspect of quantum many-body systems, pinpointing the emergence of irreversibility to critical features. In this work, we investigate the critical fingerprints appearing in these key thermodynamical quantities for a mean-field critical system undergoing a finite-time cycle, starting from a thermal state at a generic inverse temperature. In contrast to spatially extended systems, the presence of a mean-field critical point in a finite-time cycle leads to constant irreversible work even in the limit of infinitely slow driving. This links with the fact that a slow finite-time cycle results in a constant amount of squeezing, which enables us to derive analytical expressions for the work statistics and irreversible entropy, depending solely on the mean-field critical exponents and the functional form of the control parameter near the critical point. We find that the probability of observing negative work values, corresponding to negative irreversible entropy, is inversely proportional to the time the system remains near to the critical point, and this trend becomes less pronounced the lower the temperature of the initial thermal state. Finally, we determine the irreversibility traits under squeezing generation at zero-temperature using the relative entropy of coherence.

This work introduces a theoretical framework to model the collective dynamics of quantum emitters in highly non-Markovian environments, interacting through the exchange of photons with significant retardations. The formalism consists on a set of coupled delay differential equations for the emitter’s polarizations $\sigma^\pm_i$, supplemented by input-output relations that describe the field mediating the interactions. These equations capture the dynamics of both linear (bosonic) and nonlinear (two-level) emitter arrays. It is exact in some limits$-$e.g., bosonic emitters or generic systems with up to one collective excitation$-$and can be integrated to provide accurate results for larger numbers of photons. These equations support a study of collective spontaneous emission of emitter arrays in open waveguide-QED environments. This study uncovers an effect we term cascaded super- and sub-radiance, characterized by light-cone-limited propagation and increasingly correlated photon emission across distant emitters. The collective nature of this dynamics for two-level systems is evident both in the enhancement of collective emission rates, as well as in a superradiant burst with a faster than linear growth. While these effects should be observable in existing circuit QED devices or slight generalizations thereof, the formalism put forward in this work can be extended to model other systems such as network of quantum emitters or the generation of correlated photon states.

The Quantum Approximate Optimization Algorithm (QAOA) is a variational ansatz that resembles the Trotterized dynamics of a Quantum Annealing (QA) protocol. This work formalizes this connection formally and empirically, showing the angles of a multilayer QAOA circuit converge to universal QA trajectories. Furthermore, the errors in both QAOA circuits and QA paths act as thermal excitations in pseudo-Boltzmann probability distributions whose temperature decreases with the invested resource — i.e. integrated angles or total time — and which in QAOA also contain a higher temperature arising from the Trotterization. This also means QAOA and QA are cooling protocols and simulators of partition functions whose target temperature can be tuned by rescaling the universal trajectory. The average cooling power of both methods exhibits favorable algebraic scalings with respect to the target temperature and problem size, whereby in QAOA the coldest temperature is inversely proportional to the number of layers, $T\sim 1/p$, and to the integrated angles — or integrated interactions in QA.

We consider a linear array of trapped ions subjected to local parametric modulation of the trapping potential and continuous laser cooling. In our model, the phase of the parametric modulation varies linearly along the array, breaking time-reversal symmetry and inducing non-trivial topological effects. The linear response to an external force is investigated with the Green’s function formalism. We predict the appearance of topological amplification regimes in which the trapped ion array behaves as a directional amplifier of vibrational excitations. The emergence of topological phases is determined by a winding number related to non-Hermitian point-gap topology. Beyond its fundamental interests as a topological driven-dissipative system, our setup can be used for quantum sensing of ultra-weak forces and electric fields. We consider a scheme in which a trapped ion at one edge of the array acts as a sensor of an ultra-weak force, and the vibrational signal gets amplified towards the last trapped ion, which acts as a detector. We consider arrays of 2-30 $^{25}$Mg$^+$ ions, assuming that the detector ion’s displacement is measured via fluorescence with a spatial resolution of 200-500 nm, and predict sensitivities as small as 1 yN $\cdot$ Hz$^{-1/2}$. Our system has the advantage that the detected force frequency can be tuned by adjusting the frequency of the periodic drive.