Publications 2026

We propose and validate a framework for analog simulation of the Heisenberg spin model using a circuit quantum electrodynamics (circuit-QED) platform. To this end, we develop a continuous family of deformed boson representations of the SU(2) algebra, which includes the Holstein-Primakoff and Dyson-Maleev transformations as special cases. For spin-1/2 systems, we introduce a procedure to circumvent the inherent non-Hermiticity of the representation, showing that this entire family yields the extended Bose-Hubbard (EBH) Hamiltonian. For the experimental realization of this EBH model, we design a scalable circuit-QED architecture based on an engineered Josephson junction array. Numerical simulations confirm that the microwave photon dynamics in this simulator accurately reproduces the original spin dynamics. Our work establishes an experimentally accessible method for investigating complex quantum spin dynamics in a highly controllable bosonic setting.

Recent advances are enabling quantum emitters to interact with photonic crystals, whose electromagnetic modes exhibit complex dispersion relations, spatial mode structure, and polarization textures. However, modeling light-matter behavior in these systems faces a persistent trade-off: electromagnetic approaches based on Maxwell-equation solvers provide realistic vectorial descriptions but are difficult to integrate with quantum many-body and non-perturbative methods, whereas simplified quantum-optical lattice models are tractable but typically rely on scalar and spatially independent light-matter couplings that miss essential features of these structured photonic environments. Here, we introduce a constructive framework to derive quantum-optical lattice descriptions that overcome this trade-off. Combining symmetry-constrained tight-binding constructions with numerically computed photonic band structures and field profiles, our method yields minimal, symmetry-enforced lattice Hamiltonians that reproduce the target photonic dispersion while retaining the mode-resolved (position- and polarization-dependent) structure of the light-matter coupling. We show that these models recover Green’s-function-based emitter dynamics in the perturbative regime, while providing access to non-perturbative quantum dynamical simulations beyond emitter-only descriptions. As a proof of principle, we apply the framework to a two-dimensional photonic crystal and show that it captures polarization-dependent directional emission inaccessible to scalar models, while enabling the analysis of non-Markovian light-matter dynamics and entanglement. Our results provide a practical bridge between classical electromagnetic simulation tools and quantum-optical many-body and non-Markovian modeling in photonic crystal settings.

Combining decoherence protection with directional photon emission in a single waveguide quantum electrodynamics (QED) device remains an open challenge. Here we show that an artificial giant molecule — strongly interacting artificial atoms coupled to a photonic waveguide at multiple spatially separated points — achieves both: a fully operational decoherence-free (DF) qubit and state-dependent chiral single-photon emission, arising from the same photon-interference mechanism. Initialization reduces to a local excitation of a single atom, universal single-qubit gates are implemented by modulating a single atomic frequency, and readout exploits state-dependent chiral emission with directionality reaching 100% and low measurement error of 1.2%. The coexistence of decoherence protection and directional emission in a single device positions giant molecules as protected chiral nodes for modular quantum networks in waveguide QED.

We study the noisy dynamics of periodically driven, discrete-step quantum walks in a one-dimensional photonic lattice. We find that in the bulk, temporal noise that is constant within a Floquet period leads to decoherence-free momentum subspaces, whereas fully random noise destroys coherence in a few time steps. When considering topological edge states, we observe decoherence no matter the type of temporal noise. To explain these results, we derive a non-perturbative master equation to describe the system’s dynamics. We experimentally confirm our findings in a time multiplexed photonic lattice implemented in a double-fiber ring setup subject to laser pulse input states, in which we engineer different types of temporal noise.

Bosonic quantum error correction enables hardware-efficient protection of quantum information by encoding logical qubits in harmonic oscillators. Bosonic grid states, such as Gottesman-Kitaev-Preskill (GKP) states, are particularly promising due to their potential to correct small displacements and boson loss. However, their generation remains challenging, typically relying on probabilistic protocols or auxiliary qubit systems. Here, we propose deterministic protocols for generating bosonic grid states using programmable nonlinear bosonic circuits composed solely of squeezing, displacement, and Kerr operations. We show that aiming to enforce GKP symmetries in the output of these circuits yields states with competitive performance with respect to current realizations, but whose quality saturates with increasing circuit depth due to imperfect symmetry restoration. Instead, we find that these bosonic circuits naturally give rise to a distinct class of states, that we label as phased-comb states, which are unitarily related to standard grid states but feature an intrinsic phase structure. We demonstrate that these states define a scalable bosonic quantum error-correcting code with near-optimal performance under boson loss comparable to that of approximate GKP states. We further analyze their logical operations and show how to implement a universal gate set for them. Our results establish programmable nonlinear bosonic circuits as a viable route towards the generation of scalable bosonic quantum error-correcting states beyond standard GKP encodings.

We analyze the autonomous stabilization of remote entanglement by driving two distant qubits with the output of a correlated photon source. By treating the qubits as idealized entanglement detectors, we develop a general framework to quantify the maximum amount of entanglement that can be remotely stabilized in this way with a given photon source. We then apply this approach to evaluate the suitability of a single driven two-level system as a minimal resource for autonomous entanglement distribution schemes. While our analysis confirms the presence of distributable entanglement in the Mollow sidebands of a bare two-level system, we show that stabilizing close to maximally entangled states requires additional filter cavities that enhance the relevant correlated emission events compared to other processes. We identify optimized driving and cavity parameters and explain the achievable amount of entanglement in different regimes in terms of an effective two-mode squeezing model. These insights are particularly relevant for quantum networks based on photons or phonons in solid-state systems, where isolated spins, impurity centers, or other two-level defects are readily available, while alternative sources of correlated photons are difficult to realize.

The tunable interaction between stationary quantum bits and propagating modes of light allows for the encoding of quantum information in the state of itinerant photons. This ability fulfills a central requirement for quantum networking, enabling quantum state transfer between distant quantum devices. Conventionally, a symmetric envelope of the photon wavepacket is used for such purposes. Yet, the use of alternative \textit{temporal modes} enables multiple applications in waveguide quantum electrodynamics that remain unexplored experimentally. Here, we use superconducting quantum circuits to generate individual itinerant microwave photons shaped in three mutually orthogonal temporal modes. We transfer the created photons across a 30-m cryogenic link, showing that the orthogonality allows us to decide at the receiver which mode to absorb, reflecting the other two with a selectivity ratio of 40. This experimental capability extends the microwave-frequency quantum communication toolbox, enabling a new photonic degree of freedom.

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.

Quantum technologies in the terahertz (THz) require a coherent interface between addressable qubits and THz quantum channels — a capacity that so far, remains largely underdeveloped. Here, we propose and demonstrate the generation of steady-state entanglement between polar quantum emitters, mediated by THz photons. We exploit strong visible-light driving of the emitters to create Rabi-split dressed eigenstates whose energy separation can be optically tuned into the THz regime. The polar nature of the emitters activates THz transitions within these eigenstates, allowing them to couple to a THz photonic mode that induces collective dissipative dynamics. A coherent driving and control of these effective THz emitters is achieved by using a sideband optical drive with detuning close to the THz transition frequency. The resulting interplay of collective dissipation and driving activates a mechanism to generate steady-state entanglement with high values of the concurrence ($C>0.9$), attainable under experimentally feasible parameters. Crucially, both coherent manipulation and quantum state tomography are implemented entirely through optical means, avoiding direct THz control and detection. This establishes a hybrid visible-THz quantum interface in which a THz channel mediates qubit-qubit entanglement (a key operational requirement for THz quantum technologies) while remaining optically accessible.

Conformal field theory provides universal predictions for Loschmidt amplitudes following quenches from product states to critical Hamiltonians. Building on this observation, we develop a route to extracting conformal data from real-time dynamics without preparing critical low-energy states. After analytic continuation, the Loschmidt amplitude is described by a boundary-CFT partition function on a strip, whose transverse transfer matrix encodes both the boundary operator spectrum and the central charge. Local space-time perturbations of the amplitude are governed by equilibrium correlation functions, and therefore provide access to critical exponents. In parallel, generalized temporal entropies exhibit scaling with time analogous to the equilibrium scaling of spatial entanglement entropy. We show that the low-lying boundary spectrum can be reconstructed from the system-size dependence of finite-chain Loschmidt echoes, whose damped oscillations encode differences of boundary scaling dimensions. Finally, we propose a finite-size scaling protocol that can extract these quantities from simulations or experiments on state-of-the-art quantum platforms.

Excitation of trapped-ion hyperfine qubits with fast optical Raman pulses enables faster-than-trap-period entangling gates with qubits of long coherence time for practical quantum computation. Achieving high-fidelity fast two-qubit gates requires high-quality spin-dependent kicks (SDKs), which form their fundamental building blocks. Here, we characterize the control parameters (including Raman frequency difference, pulse arrival times, Lamb–Dicke parameter, temperature, pulse width, and SDK time) that maximize the performance of single-ion SDKs for protocols compatible with performed experiments involving a small number of fast pulses. We demonstrate through analytical methods and numerical simulations that, within the model commonly used for infidelity optimization, finite pulse duration is the dominant source of error, exceeding the contribution of secular motion by orders of magnitude for nanosecond-scale SDKs. Low infidelities — below $10^{-3}$ for schemes with $\gtrsim10$ fixed-amplitude, equispaced, picosecond pulses — are achievable in SDK times on the order of nanoseconds. These results provide quantitative design rules for achieving competitive SDK fidelities with current pulsed-laser technology, laying the foundation for sub-microsecond trapped-ion quantum entangling operations.

Characteristic functions of weighted sums of independent random variables exhibit low-rank structure in the quantized tensor train (QTT) representation, also known as matrix product states (MPS), enabling up to exponential compression of their fully non-Gaussian probability distributions. Under variable independence, the global characteristic function factorizes into local terms. Its low-rank QTT structure arises from intrinsic spectral smoothness in continuous models, or from spectral energy concentration as the number of components $D$ grows in discrete models. We demonstrate this on weighted sums of Bernoulli and lognormal random variables. In the former, despite an adversarial, incompressible small-$D$ regime, the characteristic function undergoes a sharp bond-dimension collapse for $D \gtrsim 300$ components, enabling polylogarithmic time and memory scaling. In the latter, the approach reaches high-resolution discretizations of $N = 2^{30}$ frequency modes on standard hardware, far beyond the $N = 2^{24}$ ceiling of dense implementations. These compressed representations enable efficient computation of Value at Risk (VaR) and Expected Shortfall (ES), supporting applications in quantitative finance and beyond.

Interfaces in three-dimensional many-body systems can exhibit rich phenomena beyond the corresponding bulk properties. In particular, they can fluctuate and give rise to massless low energy degrees of freedom even in the presence of a gapped bulk. In this work, we present the first tensor-network study of the paradigmatic interface roughening transition of the 3-D Ising model using highly asymmetric lattices that are infinite in the $(xy)$ direction and finite in $z$. By reducing the problem to an effective 2-D tensor network, we study how truncating the $z$ direction reshapes the physics of the interface. For a truncation based on open boundary conditions, we demonstrate that varying the interface width gives rise to either a $\mathbb{Z}_2$ symmetry breaking transition (for odd $L_z$) or a smooth crossover(for even $L_z$). For antiperiodic boundary conditions, we obtain an effective $\mathbb{Z}_q$ clock model description with $q=2L_z$ that exhibits an intermediate Luttinger liquid phase with an emergent $\U(1)$ symmetry.

The failure of quantum thermalization due to Many-Body Localization (MBL) has evolved from a theoretical concept in spin chains to an experimental reality in synthetic quantum platforms, most notably superconducting circuits based on transmon qubits. Despite its significance, the MBL transition has been studied primarily under purely random disorder and local couplings, leaving more complex and realistic configurations largely unexplored. Here, we investigate the quantum dynamics of transmon networks subject to the unavoidable competing effects of correlated disorder and network-mediated non-local interactions. First, we show how to engineer the physical parameters of the quantum hardware to systematically control the correlated disorder patterns emerging in the system. Then, we demonstrate that the MBL phase transition is robust against such correlations, which is essential for tuning localization properties in realistic transmon devices. This robustness is established through the analysis of the block entanglement entropy variance across disorder realizations, which precisely locates the MBL critical point. Independently, we introduce a local memory parameter, whose dynamics at long evolution times reveals memory retention in the localized phase and yields a critical point consistent with the entropy analysis. These results provide a framework for understanding localization in complex quantum network topologies, with potential implications for multi-qubit processor design.

We propose a novel experimental protocol to measure generalized temporal entropies in many-body quantum systems. Our approach involves using local operators as probes to characterize the out-of-equilibrium dynamics induced by a geometric double quench on a replicated system. Such protocol mimics the path-integral on the corresponding Riemann surface encoding generalized temporal entanglement. We present the results of tensor network simulations of one-dimensional systems which validate the protocol and demonstrate the experimental feasibility of measuring generalized temporal entropies, and we outline the experimental requirements for implementing these quenches using state-of-the-art quantum simulators. Therefore, our results provide a physical interpretation of the meaning of generalized temporal entropies. Furthermore, they reveal that the dynamics induced on two replicas of the Ising model in a transverse field differ qualitatively from the ones of its non-integrable extension, suggesting that generalized temporal entropies can be used as a tool for identifying different dynamical classes in quantum systems.

We study temporal correlations in interacting quantum systems through the influence functional of a half-infinite quantum Ising chain. Using R\’enyi entropies and temporal mutual information, we confirm that integrable dynamics is captured by the quasiparticle picture. In contrast, generic ergodic Hamiltonian dynamics exhibits pronounced deviations from random-circuit universality, and its generalization including a symmetry accounting for energy conservation. Instead, we find a long mesoscopic regime suggestive of a slow spectral reorganization of the influence functional. Our results reveal a rich temporal structure in generic Hamiltonian dynamics and point to limitations of existing random-circuit paradigms at experimentally and numerically relevant timescales.

We describe a novel scheme for the generation of stationary entanglement between two separated qubits that are driven by a purely thermal photon source. While in this scenario the qubits remain in a separable state at all times when the source is broadband, i.e. Markovian, the qubits relax into an entangled steady state once the bandwidth of the thermal source is sufficiently reduced. We explain this phenomenon by the appearance of a quasiadiabatic dark state and identify the most relevant nonadiabatic corrections that eventually lead to a breakdown of the entangled state, once the temperature is too high. This effect demonstrates how the non-Markovianity of an otherwise incoherent reservoir can be harnessed for quantum communication applications in optical, microwave, and phononic networks. As two specific examples, we discuss the use of filtered room-temperature noise as a passive resource for entangling distant superconducting qubits in a cryogenic quantum link or solid-state spin qubits in a phononic quantum channel.

Cooperative emission is well understood for idealized symmetric systems, but its limits in spatially extended, free-space ensembles remain an open question. Here, we derive a universal law for the scaling of the maximum photon emission rate with system size that unifies both ordered arrays and disordered atomic clouds in arbitrary dimensions at fixed density. We demonstrate that, for a fixed atomic density, the maximum emission rate scales universally as the product of the atom number and the system’s optical depth, with the latter encoding the dimensional scaling across all regimes from independent emission to the Dicke limit. Furthermore, we establish a scaling law for directional detection, revealing that the observed rate depends on the detector’s numerical aperture: small apertures yield Dicke-like quadratic scaling, whereas large apertures recover our integrated universal bound. Our results establish optical depth as the parameter governing many-body cooperative emission in both ordered and disordered ensembles, and reveal that directional and total-emission scalings must be carefully distinguished in experimental settings.

The short quantum link regime, where the photon travel time $\tau$ is comparable to the emitter lifetime $1/\gamma$, is experimentally relevant but theoretically underexplored: existing few-mode descriptions lose validity as retardation and multimode effects become significant. Using a Delay Differential Equation (DDE) framework that admits exact analytical solutions from the single-mode cavity limit to the multimode waveguide continuum, we show that emitters coupled to a short link spontaneously lock into self-synchronized Rabi oscillations driven by coherent photon echoes, breaking the link’s discrete time-displacement symmetry. The resulting spectral structure — persistent quasi-dark states and vacuum Rabi splitting, including in the superstrong coupling regime — enables efficient quantum state transfer (QST): benchmarking three protocols across the full $\gamma\tau$ parameter space, we find that STIRAP exploits the quasi-dark-state structure to achieve a quadratic infidelity floor $\mathcal{O}((\gamma\tau)^2)$, outperforming both SWAP (linear error $\mathcal{O}(\gamma\tau)$) and wavepacket engineering for $\gamma\tau \lesssim 1.44$, even in regimes where retardation cannot be neglected. These results establish photon-echo synchronization as an engineering resource for quantum state transfer, with DDE modeling providing the exact analytical predictions needed to design and optimize short-link experiments on current circuit-QED hardware.

Subwavelength arrays of atoms trapped in optical lattices or tweezers are inherently susceptible to deformations: Optomechanical forces displace atoms within their trapping potential and produce lattice distortions, which in turn modify the optical response of the array. We show that this optomechanical coupling hybridizes collective atomic excitations (polaritons) with phonons, forming polaron-polaritons — the fundamental quasiparticles governing light-matter interactions in arrays of trapped atoms. Using analytical polaron theory and numerical simulations, we find that: (1) phonons can strongly enhance the decay of subradiant states, but also enable their efficient excitation; (2) transport of dark excitations remains remarkably robust even at low trap frequencies, except when a polariton can resonantly scatter phonons; and (3) motion reduces the reflectivity of a two-dimensional atomic mirror; by identifying design principles that mitigate this degradation, we recover reflectivity above 99% under realistic conditions. Our findings lay the foundation for analyzing motional effects in key applications and suggest new ways to harness them in state-of-the-art experiments.

Quantum dot based platforms offer a promising route towards realizing the Kitaev chain Hamiltonian hosting Majorana bound states (MBSs). Poor man’s MBSs arise in a two-site Kitaev chain when the parameters of the system are fine-tuned to the sweet spot. Based on our previous work [Phys. Rev. B 111, 155410 (2025)], we consider a microscopic model for the Kitaev chain based on quantum dots with proximity effect embedded in a photonic cavity. We find that the photon coupling in the microscopic model yields an effective Hamiltonian where the cavity affects the pairing term. However, we demonstrate that even in this case, it is possible to screen particle interactions and reach the sweet spot condition for the emergence of the poor man’s MBSs. In particular, we find that attractive particle interactions can be canceled for the cavity prepared in the zero-photon state, while repulsive ones can be screened with a cavity prepared in the one-photon state. Furthermore, in case of a large number of photons in the cavity, we find that the hopping amplitudes are suppressed resulting in a degenerate spectrum. This motivates the use of quantum light for engineering poor man’s MBSs with cavity embedding.

Long-range interactions are the source of many equilibrium and out-of-equilibrium quantum many-body phenomena. Analog simulators based on ionic, atomic, superconducting, and molecular systems provide a natural platform to obtain these interactions using vibration- and photon-mediated processes. Recent experimental advances, such as their integration in multi-mode cavities and waveguides, or the use of Raman-assisted transitions, enable dynamical control over both the strength and the spatial range of these interactions, thereby rendering them programmable. Here, we develop a hybrid classical-quantum toolbox that exploits this tunability to enhance many-body state preparation in analog simulators beyond fixed-connectivity architectures. Our approach is based on classical pre-compilation in homogeneous small systems, whose optimized parameters are extrapolated iteratively to larger system sizes, and then refined on the quantum hardware using noise-aware hybrid re-optimization and error-mitigation techniques. We benchmark this strategy across several fermionic, spin-1/2, and spin-1 models, demonstrating orders-of-magnitude improvements in fidelity and energy estimation for system sizes ranging from 100 to 1000 particles. Finally, we show that the combination of such high-fidelity programmable state preparation techniques with tunable-range out-of-equilibrium dynamics enables controlled studies of many-body thermalization in regimes accessible to current experimental platforms. Our results establish programmable long-range interactions as a powerful resource for next-generation analog quantum simulators.

Quantum battery, as a novel energy storage device, offers the potential for unprecedented efficiency and performance beyond the capabilities of classical systems, with broad implications for future quantum technologies. Here, we experimentally \RefC{demonstrate quantum charging advantage (QCA)} in a scalable solid-state quantum battery. More specifically, we show how double-excitation Hamiltonians for two-level systems promote scalable QCA \RefB{with standard methods.} We effectively implement the collective evolution of quantum systems with 2 up to 12 battery cells in a superconducting quantum processor, and study the performance of quantum charging compared to its uncorrelated classical counterpart. The model considered is a linear chain of superconducting transmon qubits with only \textit{nearest-neighbor} and \textit{pairwise} interactions, which constitute the simplest model of a multi-cell quantum battery. Our results empirically realize substantial QCA without the necessity of adopting long-range and many-body interactions \RefB{ and showcase the quantum features of the QB charging processes with measurements of non-zero coherent ergotropy, incoherent ergotropy and entanglement,} revealing a promising prospect for further developments of efficient and experimentally feasible protocols for QCA.

Continuously monitored quantum systems are emerging as promising platforms for quantum metrology, where a central challenge is to identify measurement strategies that optimally extract information about unknown parameters encoded in the complex quantum state of emitted radiation. Different measurement strategies effectively access distinct temporal modes of the emitted field, and the resulting choice of mode can strongly impact the information available for parameter estimation. While a ubiquitous approach in quantum optics is to select frequency modes through spectral filtering, the metrological potential of this technique has not yet been systematically quantified. We develop a theoretical framework to assess this potential by modeling spectral detection as a cascaded quantum system, allowing us to reconstruct the full density matrix of frequency-filtered photonic modes and to compute their associated Fisher information. This framework provides a minimal yet general method to benchmark the performance of spectral measurements in quantum optics, allowing to identify optimal filtering strategies in terms of frequency selection, detector linewidth, and metrological gain accessible through higher-order frequency-resolved correlations and mean-field engineering. These results lay the groundwork for identifying and designing optimal sensing strategies in practical quantum-optical platforms.

We, as researchers in quantum science and technology, are publishing this manifesto to express our deep concerns about the current geopolitical situation and the global race to rearm. We firmly oppose all forms of militarization in our societies and, in particular, within the academic world. We categorically reject the use of our research for military applications, population control, or surveillance. We stand against the practice of military funding for research. This manifesto is a call to action: to confront the elephant in the room of quantum research, and to unite all researchers who share our views. Our main goals are: i) To express, as a unified collective, our rejection of the use of our research for military purposes; ii) To open a debate in our community about the ethical implications of quantum research for military purposes; iii) To create a forum where concerned scientists can share their opinions and join forces in support of demilitarized research; iv) To advocate for the establishment of a public database listing all research projects at public universities funded by military or defense agencies. In what follows, we lay out our concerns and the rationale behind our opposition to the militarization of quantum research.

Digitized adiabatic quantum factorization is a hybrid algorithm that exploits the advantage of digitized quantum computers to implement efficient adiabatic algorithms for factorization through gate decompositions of analog evolutions. In this paper, we harness the flexibility of digitized computers to derive a digitized adiabatic algorithm able to reduce the gate-demanding costs of implementing factorization. To this end, we propose a new approach for adiabatic factorization by encoding the solution of the problem in the kernel subspace of the problem Hamiltonian, instead of using ground-state encoding considered in the standard adiabatic factorization proposed by Peng $et$ $al$. [Phys. Rev. Lett. 101, 220405 (2008)]. Our encoding enables the design of adiabatic factorization algorithms belonging to the class of Quadratic Unconstrained Binary Optimization (QUBO) methods, instead the Polinomial Unconstrained Binary Optimization (PUBO) used by standard adiabatic factorization. We illustrate the performance of our QUBO algorithm by implementing the factorization of integers $N$ up to 8 bits. The results demonstrate a substantial improvement over the PUBO formulation, both in terms of reduced circuit complexity and increased fidelity in identifying the correct solution.

We introduce SeeMPS, a Python library dedicated to implementing tensor network algorithms based on the well-known Matrix Product States (MPS) and Quantized Tensor Train (QTT) formalisms. SeeMPS is implemented as a complete finite precision linear algebra package where exponentially large vector spaces are compressed using the MPS/TT formalism. It enables both low-level operations, such as vector addition, linear transformations, and Hadamard products, as well as high-level algorithms, including the approximation of linear equations, eigenvalue computations, and exponentially efficient Fourier transforms. This library can be used for traditional quantum many-body physics applications and also for quantum-inspired numerical analysis problems, such as solving PDEs, interpolating and integrating multidimensional functions, sampling multivariate probability distributions, etc.

We revisit the U(1) quantum link model in a ladder geometry, finding, by finite-size scaling, that the critical exponent $\nu=1$ and the central charge $c=1/2$ are consistent with the Ising universality class for all phase transitions observed. A blind application of the Harris criterion would suggest that this criticality is lost in the presence of the disorder. It turns out not to be the case. For the disorder affecting ladder’s rung hoppings only, we have found that the transitions survive disappearing only for quite strong disorder. The disorder in the ladder’s legs destroys the nonzero mass phase criticality, while the symmetry-protected topological phase for zero mass survives a small disorder. The observed robustness against disorder is explained qualitatively using field-theoretic arguments.

Driven-dissipative quantum systems can exhibit robust transport and amplification in topological regimes, yet the connection between topology and the extent of correlations remains largely unexplored. In this work, we develop a general framework that links topological phases in driven-dissipative systems to bosonic correlations via the singular value decomposition (SVD). In essence, we claim that non-Hermitian topology in quadratic Liouvillians is directly encoded in steady-state correlations, providing an intrinsic characterization of topology without external probes. We show that topological amplification induces disorder-robust long-range order (LRO) in steady-state correlations at fixed frequency, establishing frequency-resolved correlations as direct signatures of non-Hermitian topological phases. We introduce a vector-valued topological invariant that captures the total number of singular-value gap closings across the frequency axis, extending the concept of adiabatic deformation from topological insulators to the case of topological phases of quadratic Liouvillians. Within this framework, we further demonstrate that the spatial structure of equal-time correlations encodes global topological information, manifested as a Gaussian spatial decay with distance in the topological phase, in contrast to the exponential decay characteristic of trivial phases. These findings open new avenues for quantum sensing and correlation engineering in non-Hermitian systems, with feasible implementations in platforms such as trapped ions and superconducting circuits.

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.

Variational ans\”atze are a cornerstone of quantum many-body physics, providing compact approximations to complex ground states using finite resources. Recent quantum-technology advances have introduced a new class based on layered parameterized evolutions. Assessing whether they can represent critical ground states is challenging: correlations span all length scales, while finite circuit depth limits how far they extend. Building on finite-resource scaling from tensor networks, we assign each ansatz an emergent correlation length $\xi_D$, the longest range over which it faithfully captures critical correlations. Its growth with refinement parameter $D$, $\xi_D \propto D^\kappa$, defines an exponent $\kappa$ measuring how efficiently an architecture converts resources into long-distance correlations. Applying this framework to the critical transverse-field Ising model, with $D$ the circuit depth of parameterized evolutions, we compare ans\”atze with nearest-neighbor and long-range generators, and layers where generators act separately or combined. All ans\”atze are compatible with algebraic growth, but fitted exponents range from $\kappa\simeq1$ to $\kappa\simeq3$. Exponential interactions give the largest exponents, while power-law interactions stay close to nearest-neighbor behavior, showing long-range support alone gives no scaling advantage. Combined-generator layers generally outperform separable ones, so layer organization matters alongside interaction range. Finally, a quasiparticle analysis shows finite depth leaves an unresolved window of width $\xi_D^{-1}$ around the low-energy modes responsible for long-distance correlations. The emergent correlation length thus acts as an infrared resolution scale, providing a benchmark for critical-state preparation and a guide for designing resource-efficient variational architectures.