Papers

The resource-theoretic approach to quantum thermodynamics typically assumes perfect knowledge of the thermal equilibrium state, an idealization incompatible with finite experimental precision. We develop a framework for equilibrium uncertainty by representing the equilibrium reference as a set of candidate states. Under a generic geometric condition, we prove a no-go theorem that sharply limits athermality “purification”: converting an uncertain athermal state into a definite target is either trivial or impossible. We then derive exact one-shot entropic characterizations of work extraction and formation for two work-storage models, a clean battery with known equilibrium and a dirty battery with uncertain equilibrium. Both models exhibit strong asymptotic irreversibility even under arbitrarily small uncertainty. An explicit example reveals two distinct extremes: clean batteries display a bound-entanglement-like phenomenon, with positive formation cost but zero extractable work, whereas dirty batteries allow positive work extraction but require infinite formation cost. These phenomena show that equilibrium uncertainty is not a minor perturbation of the standard theory, but a structural ingredient that fundamentally reshapes the limits of quantum thermodynamics.

Quantum simulation is a cornerstone application of quantum computing, yet how fundamental quantum resources–entanglement and non-stabilizerness (“magic")–shape simulation fidelity remains an open question. In this work, we establish a rigorous connection between these resources and the statistical behavior of algorithmic errors arising in Hamiltonian simulation based on the Trotter-Suzuki formula. By analyzing ensembles of states with fixed entanglement entropy or magic, we make two key discoveries: First, the variance of the Trotter error decreases with increasing entanglement entropy, indicating a stronger concentration of error for entangled states. Moreover, we find that the kurtosis of the error exhibits a negative linear dependence on magic, implying that states with high magic possess lighter-tailed error distributions and thus a reduced probability of large deviations. These findings reveal a subtle phenomenon: quantum resources that obstruct classical emulation may, paradoxically, enhance the intrinsic robustness of quantum simulation, highlighting a constructive interplay between complexity and stability in quantum computation.

Bright squeezed light has entered strong-field atomic physics, creating an immediate need to characterize the quantum field delivered to the target on subcycle timescales. Here we show that attosecond streaking maps the coherent displacement and phase-sensitive covariance of a displaced squeezed field onto distinct harmonics: the mean photoelectron momentum follows an $\omega$-periodic shift of the spectral center, whereas the momentum variance exhibits a $2\omega$-periodic breathing. A Feynman–Vernon formulation represents Gaussian quantum light by stochastic vector-potential trajectories, enabling Coulomb-resolved TDSE simulations with a finite XUV gate and without explicit photon-state propagation. A coherent-state reference calibrates the streaking phase and variance background, enabling retrieval of the coherent amplitude and phase together with the squeezed-noise amplitude and squeezing phase. This establishes attosecond streaking as an in situ, gas-phase diagnostic of bright squeezed light.

We theoretically investigate the statistical properties of photon emission of a driven two-level single-molecule system undergoing spectral diffusion induced by nonequilibrium environmental fluctuations. Within the framework of the generating function method and the stochastic Liouville equation, we analyze the influence of the nonequilibrium characteristics of environmental fluctuations respectively governed by nonstationary Ornstein-Uhlenbeck noise and random telegraph noise on the photon counting statistics of the driven single-molecule system. In the slow modulation limit of spectral diffusion, the intensity and statistical fluctuations of photon emission depend on the environmental nonequilibrium characteristics at short time scales, whereas they become independent of the nonequilibrium characteristics of environmental fluctuations in the steady state. In the fast modulation limit of spectral diffusion, neither the line shape nor the Mandel's parameter depends on the environmental nonequilibrium characteristics owing to the rapid relaxation of environmental fluctuations. These findings not only shed light on the role of nonequilibrium environmental fluctuations in shaping the photon emission properties of single-molecule systems but also provide a basis for distinguishing between equilibrium and nonequilibrium characteristics of environmental fluctuations in experimental measurements.

Variational quantum algorithms on bosonic quantum processors are an emerging paradigm for quantum chemistry calculations, exploiting the natural alignment between molecular structure and harmonic oscillator-based hardware. We introduce the qumode-based variational quantum deflation framework (QumVQD) for finding both electronic and vibrational excited state energies on qumode-based architectures. We validate the approach through electronic structure calculations on H$_{\text{2}}$ and linear H$_{\text{4}}$, where we introduce Hamming-weight filtering of the Fock basis to enforce particle number conservation and eliminate spurious eigenstates by reducing the required Hilbert space, which reduces the required number of qumodes in turn. We achieve agreement with full configuration interaction (FCI) using the STO-3G basis set within the chemical accuracy threshold at most points along the potential energy surfaces. Extending to the vibrational structure, we combine QumVQD with an existing Hamiltonian fragmentation approach based on Cartan subalgebra, allowing us to compute the vibrational eigenenergies of CO$_{\text{2}}$ and H$_{\text{2}}$S to spectroscopic accuracy with per-fragment circuits that scale as $O(N)$ in single-qumode gates and $O(N^2)$ in beam-splitter gates for $N$ qumodes. For the case of CO$_{\text{2}}$, we get total gate counts more than an order of magnitude smaller than those reported for qubit-based vibrational algorithms at this system size. These results demonstrate that bosonic quantum devices are a viable platform for excited-state quantum chemistry, particularly for vibrational problems where qubit-based methods incur substantial boson-to-qubit mapping overhead.

In recent years, the growing scale of programmable neutral-atom arrays has led to an increasing demand for higher-power Rydberg excitation light. Although pulsed amplifiers deliver higher peak power than continuous-wave lasers, their use for efficient coherent Rydberg excitation of single atoms in arrays has been limited by challenges such as pulse distortion, synchronization with excitation sequences, and spectral linewidth broadening. Here, we address these issues using a fiber-based master-oscillator power-amplifier system. We demonstrate efficient coherent Rydberg excitation of single atoms in a rubidium atom array, achieving performance comparable to continuous-wave methods. This study provides a potentially new technical pathway toward future large-scale quantum simulation and computation with Rydberg atom arrays.

In Si(111) crystals, a strong biaxial tensile strain applied within the (111) plane is considered to shift the lowest energy point of the conduction band from the $\Delta$ valley to the L valley. Electrons confined in this L valley experience a splitting of their quadruply degenerate energy levels into an undegenerate single-level ground state (L1) and a triply degenerate excited state (L3). The energy of the single-level ground state is sufficiently low relative to the energies of the L3 valley and the $\Delta$ valley, making it optimal as a two-level system for a qubit. Using deformation potential theory and incorporating quantum effects from electron confinement in the SiGe/Si(111)/SiGe structure, we determine the value of the biaxial tensile strain causing the shift of the conduction band energy minimum from the $\Delta$ valley to the L valley, along with the corresponding Ge concentration. We also calculate the critical thickness for the plastic relaxation of the Si quantum well under this large biaxial tensile strain and examine the feasibility of realizing it as a SiGe/Si(111)/SiGe heterostructure.

Nuclear lattice effective field theory has become an important framework for quantum many-body calculations in nuclear physics, yet its classical implementation remains increasingly challenging for more general interactions and larger systems. In this work, we develop a quantum-computing framework for a three-dimensional nuclear lattice model. We construct a variational quantum eigensolver framework and systematically compare the Jordan-Wigner and Gray code encodings. Our analysis shows that for the few-body systems considered here, Gray code combined with symmetry reduction yields a substantially more compact qubit representation. Based on this framework, we perform numerical studies for $^{2}\mathrm{H}$, $^{3}\mathrm{H}$, and $^{4}\mathrm{He}$ on finite lattices. The calculated ground-state energies exhibit a clear approach toward the corresponding experimental binding energies as the lattice size increases. These results provide a proof-of-principle foundation for future quantum simulations of nuclear many-body problems.

We investigate a many-body interacting system of quantum kicked rotors, where each rotor resides in its respective quantum resonance. Rich many-body dynamics are found to emerge from the interplay between the principal and secondary resonances. In particular, for both the wavepacket and bipartite entanglement entropy, we analytically demonstrate three distinct dynamical regimes – quadratic spreading (growth), period-2 oscillation, and their hybrid – governed by the respective symmetries of the relevant potentials. Based on these symmetries, the connection between the wavepacket and the entanglement dynamics is illustrated. Other related issues are also discussed, including higher-order resonance effects, the robustness of the predicted dynamical behaviors, extension to many-body kicked tops, and relevance to experimental studies.

The phenomena of subthreshold photoemission and absorption under coherent and entangled-photon-pair illumination are reviewed, and the generation and properties of entangled-photon pairs are surveyed. Three prominent forms of subthreshold photoemission are examined: one-photon Fermi-tail photoemission (FTP), two-photon photoemission (TPP), and entangled-two-photon photoemission (ETPP). Experimental methods for measuring subthreshold photocurrents and photoelectron count rates are discussed, along with strategies for enhancing selected contributions. Experimental observations of FTP from a CsK$_2$Sb photocathode in a photomultiplier tube (PMT), under both coherent and entangled-photon-pair illumination, are reviewed, and the role of FTP as a noise source in two-photon measurements is elucidated. TPP from Na and CsK$_2$Sb photocathodes in a PMT under classical-light illumination is considered, as are TPP and ETPP from a CsK$_2$Sb photocathode in a channel photomultiplier (CPM) under coherent and entangled-photon-pair illumination. The observation of ETPP is facilitated by the use of a CPM, which suppresses FTP, and by low-intensity illumination, which minimizes TPP. Quantum models of TPP and ETPP accord well with experiment. Entangled-two-photon absorption (ETPA) is analyzed, as are its applications in entangled-two-photon fluorescence microscopy (ETPFM) and entangled-two-photon spectroscopy (ETPS). The three principal forms of subthreshold absorption parallel those of subthreshold photoemission: singleton-induced Boltzmann-tail absorption; cousin-induced/singleton-pair-induced two-photon absorption; and twin-induced ETPA. Heuristic particle and fully quantum models of these processes are compared, and experimental studies of ETPA and ETPFM, together with methods for enhancing their observability, are summarized.

Experiments on attractive Bose-Einstein condensates (BECs) have unlocked many intriguing out-of-equilibrium dynamics through the interplay between matter-wave dispersion and nonlinear attractive interaction. Competition between these effects leads to fascinating phenomena such as wave collapse, modulational instability, and formation of multidimensional bright solitons. This chapter reviews experimental studies on attractive condensates, with a primary focus on alkali atoms featuring two-body contact interactions. We review recent experimental advances in optical trapping and interaction control techniques, which have enabled new studies on attractive condensates in three and also in lower dimensions. Specifically, we discuss pioneering and recent experimental observations on the dynamics and stability of attractive BECs, including the formation of bright solitons, their collisions, and excitations in quasi-one-dimensional traps. Recent observations of the elusive two-dimensional Townes solitons and vortex solitons are also discussed in this Chapter. We then highlight an experimental technique revealing the nonclassical signatures of modulational instability in an attractive condensate.

Superconducting quantum processors have largely converged on transmon-based architectures, while alternative qubit modalities with intrinsic error protection have lacked a demonstrated path to scalable system integration. In particular, although tunable-coupler-mediated interactions have been validated for small fluxonium systems, it remains unclear whether such designs can be scaled to a multi-qubit lattice. Here, we establish a scalable fluxonium processor architecture based on a modular qubit-coupler unit cell engineered to suppress residual interactions and spectator errors in a many-qubit lattice. The system enables parallel single-qubit gate fidelities approaching 99.99% and two-qubit CZ gate fidelities around 99%. With an optimized gate duration of 32 ns, the best CZ gate fidelity reaches 99.9%. We further validate this architecture in a 22-qubit processor based on the same configuration, where parallel operations enable the deterministic generation of Greenberger-Horne-Zeilinger states involving up to 10 qubits. Together, these results demonstrate that the fluxonium-tunable-coupler unit cell composes without emergent interaction pathologies and establish fluxonium as a scalable superconducting qubit platform.

Consciousness and quantum mechanics are among the most puzzling phenomena studied in the sciences. Some scholars suggest they are related, though others think this claim commits a "minimization of mystery" fallacy. The aim of this programmatic paper is to draw attention to a less widely discussed parallel between consciousness and quantum mechanics: both challenge the classical objectivist worldview of science. Under certain assumptions, they are each in tension with a package of metaphysical theses – "non-relationalism", "non-fragmentation", and "one world" – that jointly make up that worldview. This points to three distinct non-objectivist responses: the "relationalist", "fragmentalist", and "many-subjective-worlds" ones. We will map out their pros and cons.

Synthetic Aperture Radar (SAR) data is inherently complex-valued, while quantum machine learning (QML) models operate in complex Hilbert spaces. This similarity suggests that using both the magnitude and phase of SAR data in quantum encoding should help automatic target recognition in SAR images. In this study, we test this assumption by comparing five encoding strategies for quantum models: magnitude-only encoding, joint magnitude-phase encoding, in-phase and quadrature encoding, preprocessed phase encoding, and a purely quantum architecture. All approaches are evaluated under a unified experimental setup on the MSTAR benchmark dataset. Surprisingly, we find that magnitude-only encoding performs better than phase-inclusive encodings in hybrid quantum-classical models. It achieves 99.57 percent accuracy on the 3-class task and 71.19 percent accuracy on the 8-class task, outperforming complex-valued alternatives under the same framework. Adding phase information provides little or no improvement and can sometimes degrade performance. However, in purely quantum models with only 184 to 224 trainable parameters and no classical neural-network layers, phase information becomes much more important, improving accuracy by up to 21.65 percentage points. These findings show that the usefulness of phase information depends not only on the data, but also on the architecture used to process it. Hybrid models can compensate for missing phase information through their classical components, while pure quantum models rely more strongly on phase information for class discrimination. The results provide practical guidance for encoding complex-valued data in quantum machine learning and highlight the importance of jointly designing encoding strategies and model architectures for current quantum systems.

A physical relationship for the combination of scrambling via spread complexity and entanglement is characterized by the fidelity of quantum unitary dynamics. Time-integrated quantities for the degree of quantum ergodicity capture state spreading from initial to late times under physically plausible operator perturbations. Selecting maximally entangled state as a initial condition evaluates the scrambling of entanglement in this setting. For this reason, we utilize the Rosenzweig-Porter ensembles across different ergodic regimes. Computed integrated spread complexity and associated integrated quantum state fidelity display a monotonic inverse relationship for maximally entangled states from integrability to chaos. Information bounds in non-local high-energy systems and holographic scrambling can be evaluated with this approach.

If an operator $H$ anticommutes with a chirality operator $\Gamma_*$ such that $\Gamma_*^2=1$, the null space of $H$ can be decomposed in a direct sum of two spaces having positive and negative chiralities, respectively. When both spaces are finite dimensional, one can define an index, $\mathrm{Ind}(\Gamma_*,H)$, as the difference of dimensions of these two spaces. The key issue is whether $\mathrm{Ind}(\Gamma_*,H)$ is topologically protected, i.e., whether it remains constant under smooth variations of the parameters and background fields entering $H$. For Hermitian Dirac operators, topological protection of the index is guaranteed by the Atiyah–Singer theorem. In this paper, by using the heat kernel methods, we show that $\mathrm{Ind}(\Gamma_*,H)$ is topologically protected also for non-hermitian operators $H$ as long as they are diagonalizable and satisfy some ellipticity conditions. The index is given by a heat kernel coefficient which is an integral of a local expression depending on background fields.

We evaluate the effective action for the Dynamical Casimir Effect (DCE) for a real scalar field in $d+1$ dimensions within the worldline formulation of quantum field theory. The scalar field is coupled to a spacetime-dependent mass term, which here plays the role of the moving medium and imposes imperfect boundary conditions on time-dependent surfaces. Expanding in powers of the departure of the geometry from a planar configuration, the worldline path integral factorizes into simpler, lower-dimensional ones. In the limit of a strong coupling to the surface, we recover the Dirichlet result and derive the systematic corrections in inverse powers of the coupling, calculating the imaginary part of the effective action up to arbitrary order of said powers. Finally, we also apply the method to a two-surface configuration.

A fully fault-tolerant implementation of the quantum error-detecting Iceberg $[[2m, 2m-2, 2]]$ code applied to a Toffoli circuit achieved beyond-break-even error detection on a leading trapped-ion quantum computer, where the effect of encoding a circuit with a quantum error-detection code enables increased fidelity compared to an unencoded circuit. This code was also applied to Bell state preparation circuits, where a lean non-fault-tolerant implementation of the Iceberg code enables a fidelity gain as well. This highlights the important point that, at least for small-scale circuits with a substantial portion of error-free runs, it can be effective simply to use error detection to filter out the runs with errors. Furthermore, experiments performed in this work highlight the necessity for judicious compilation of circuits not only for a given hardware but also within a quantum error detection code.

We discuss the semiclassical approximation to transport problems in quantum chaotic systems. The figures of merit are moments of the transmission matrix and of the time delay matrix. After reviewing a few results obtained by treating these matrices are random matrices, we show how expressions for their elements in terms of sums over trajectories lead to diagrammatic formulations that correspond to perturbative calculations. This semiclassical approach agrees with random matrix theory when it should, and allows further elements to be incorporated, like tunnel barriers, superconductors, absorption effects. We also discuss how this approach can be encoded in matrix integrals, resulting in a powerful and versatile theory that is amenable to algebraic solutions.

Quantum phase measurements offer a complementary route to axion searches. We show that axion-photon interactions can imprint both Aharonov-Bohm (AB) and Berry phases in experimentally motivated quantum setups. For a coherently oscillating axion dark matter background, the induced effective current generates a time dependent magnetic flux in an rf-SQUID, leading to a measurable voltage signal through the Josephson phase. For representative benchmarks, this AB phase search reaches the minimum axion-photon coupling $g_{a\gamma\gamma}^{\mathrm{min}}\sim 7.8\times10^{-14}~\mathrm{GeV}^{-1}$ at axion mass $m_a\sim 10^{-10}~\mathrm{eV}$, with projected sensitivity that can improve on existing limits in that parameter space by roughly one to two orders of magnitude. We also identify a geometric phase observable in a Mach-Zehnder interferometer with an adiabatically rotating magnetic field, providing a proof-of-principle phase-based probe of meV-scale axions even when they do not constitute the dark matter, although sensitivity on the coupling remains weaker than current bounds with conservative tabletop benchmarks. Extending the analysis to a three level photon-axion quasiparticle (AQP)-axion system, with the AQP realized in a topological magnetic insulator, we find a potentially measurable THz Berry phase dominated by the AQP sector, furnishing a nontrivial validation of the formalism in a richer coupled system. These setups establish quantum phase observables as a useful new framework for axion searches, with immediate phenomenological promise in superconducting circuits and longer term potential in quantum enhanced interferometry.

Quantum computational sensing (QCS) combines quantum sensing with quantum computing to extract task-relevant information from the physical world. QCS can in principle achieve an accuracy advantage for specific tasks versus the alternative of raw-signal estimation using conventional quantum sensing followed by task-specific classical postprocessing. Here we report the experimental demonstration of quantum computational displacement sensing (QCDS) with a superconducting circuit comprising a qubit coupled to an oscillator. We consider binary classification sensing tasks, where the goal is to predict the class label of a single complex-valued displacement sensed once by the oscillator. Rather than estimating the displacement, our computational-sensing protocol – using parameterized quantum circuits before and after sensing – attempts to determine the binary class label using quantum processing and map it onto the ground or excited state of the qubit. A single measurement of the qubit directly outputs the prediction. We implemented circuits with up to 24 entangling gates and 38 free parameters, which were trained in silico. We show that increasing the circuit depth systematically improves expressivity and classification accuracy. We experimentally obtained an accuracy advantage over a suite of protocols that first use conventional quantum sensing to estimate the displacement before using classical postprocessing to perform prediction. For certain tasks, our protocol achieves a 15-percentage-points higher classification accuracy than most conventional approaches considered. Our results establish the feasibility of quantum computational sensing with noisy superconducting hardware and illustrate how integrating quantum computation with quantum sensing can enhance performance when the goal is to estimate a property or function of a signal rather than to estimate the signal.

Quasiparticle poisoning following particle impacts poses a significant challenge to the development of fault-tolerant superconducting quantum computers, as a sudden excess of quasiparticles can simultaneously degrade the coherence of multiple qubits across large device arrays. In this work, we present a statistical analysis that models the time evolution of radiation-induced qubit energy relaxation through quasiparticle density dynamics. This study provides insight into quasiparticle loss processes by distinguishing between recombination and trapping decay channels and assessing their respective impact on qubit performance. We precisely measure quasiparticle recombination in multiple transmon qubits and uncover an unexpected dependence of qubit relaxation dynamics on deposited energy. By linking correlated relaxation events across qubits to ballistic phonon propagation, we introduce a statistical localization approach to extract the energy deposited in the substrate, which is in good agreement with Monte Carlo simulation. This work establishes the quantitative framework for using an arbitrary subset of superconducting transmon qubits in a QPU as energy-resolving witness particle detectors.

We establish a rigorous operator-space condition for the thermalization of typical initial states. We introduce a concept of simple slow operators (SSOs), defined as operators that evolve slowly under Heisenberg-picture dynamics and have significant small-sized components. We show that the non-thermalization of typical initial states drawn from a weakly entangled state ensemble implies the existence of SSOs; equivalently, the absence of SSOs implies that typical initial states thermalize. We establish these results by introducing the concept of an ensemble variance norm of an operator, defined as the typical magnitude of the expectation value of that operator with respect to states in the ensemble. For low-entanglement ensembles, the norm is related to operator sizes, allowing us to establish a direct link between operator dynamics and thermalization.

In quantum systems, correlations are subject to monogamy relations between the whole and its parts. In particular, entanglement between two sets of particles constrains the entanglement within each set and viceversa. This raises the following question: In a multipartite quantum system, can global correlations disentangle all subsystems of certain size? By exploiting unitary symmetries, we formulate this question as a spectral marginal problem, namely the existence of $n$-qubit pure states whose marginals below a threshold cannot be entangled by any unitary evolution. We call such states threshold entanglement sharing (TES) states. We provide constructive existence proofs and bound the marginal sizes for up to 200 qubits. Further, we show that TES states contain both entanglement and magic, the main resources for quantum computing. Our proof strategy bridges techniques that were treated separately up to date, and independently resolves two open problems in quantum information theory: maximal multipartite entanglement and maximum purity of absolutely separable states.

Randomized measurements access nonlinear functionals without full tomography, yet turning third-order local single-copy data into a strong entanglement test remains difficult. We convert the reduction criterion into an experimentally measurable separability criterion by testing it on squared affine combinations of the identity, the local marginals, and the state itself. This yields a $4\times4$ matrix $\bar{\mathfrak{M}}(\rho)$ built from experimentally accessible second- and third-order local invariants. Entanglement is certified when its minimum eigenvalue $\mathcal{E}_4(\rho)$ becomes negative. We prove that all separable states satisfy $\bar{\mathfrak{M}}(\rho)\succeq0$, and that the sign of $\mathcal{E}_4(\rho)$ can be inferred from single-copy randomized measurements with dimension-independent sample complexity. For isotropic states on $d\times d$, the second-order purity criterion detects entanglement only for $p\sim d^{-1/2}$, whereas our third-order witness reaches $p\sim 2/d$, close to the separability threshold $p\sim 1/d$. A complementary nonisotropic benchmark shows that the affine marginal directions become essential once the local states are not maximally mixed.

Unstable periodic orbits act as organizing structures for classical chaotic systems and underpin quantum scarring. Long known in single-particle systems, genuine quantum scars based on unstable periodic orbits have been recently extended to isolated many-body systems for time-independent Hamiltonians. Their fate under periodic driving, however, remains largely uncharted, challenged by the expectation that these systems should in general heat to infinite temperature. Here, we investigate how genuine scarring competes with the drive in a Floquet many-body system. Using chaotic spin chains, we demonstrate that Floquet states remain scarred in the high-frequency limit. Beyond this static correspondence, we uncover additional, driving-induced Floquet scars with no static analog. We construct a rich dynamical stability diagram with intermediate-frequency regimes of enhanced and quenched scarring, which we understand with a classical analysis of the Lyapunov exponent. Our results position Floquet systems as a natural platform for tuning the scarring behavior of quantum many-body systems.

We introduce a framework for realizing universal fermionic quantum processing with globally controlled itinerant fermionic particles. Our approach is tailored to the example of neutral atoms in optical lattices, but transposes to other setups with similar capabilities. We give constructive protocols to realize arbitrary fermionic processes, with time-dependent control over global parameters of the experimental setup, such as tunneling and interaction in a Fermi-Hubbard type model. We first prove the universality of our framework and then discuss implementation variants, such as hybrid analog-digital simulation of extended Fermi-Hubbard models, e.g., with long-range couplings.

Measurement time represents a critical bottleneck limiting the operational speed of neutral atom quantum computers, as it cannot be accelerated through parallelization like other quantum operations. We present a protocol for fast measurement of neutral atoms based on a new, fast multi-atom Rydberg gate that significantly reduces the measurement integration time and improves the measurement fidelity. Our approach employs a multi-qubit register of $N$ ancilla atoms within a single Rydberg blockade region to measure a single data qubit. This enables an $N$-fold enhancement in photon emission collections, while reducing the measurement's sensitivity to loss. The scheme requires spectral separation between the data qubit and the ancillae, achievable through either a dual-species architecture or a targeted light shift. Beyond this, the scheme is straightforward to implement: it relies only on global pulses, global photon collection, and avoids both atom shuttling and numerically optimized pulses. Simulations of a Cs–Rb platform demonstrate that with only five ancillae ($N=5$), measurement infidelity below $10^{-3}$ within $6\ \mu\text{s}$ is achievable.

We introduce the flavoured lattice Schwinger model, a $(1{+}1)$-dimensional $U(1)$ lattice gauge theory in which the fermion doubling problem is resolved by staggering a $\mathbb{Z}_{2}$ flavour degree of freedom rather than staggering chirality. Unlike the standard approaches, this construction preserves an exact axial $U(1)$ symmetry at finite lattice spacing. We derive the continuum limit, showing that the model reduces to the two-flavour massless Schwinger model, with flavours $\alpha\in\{0,1\}$ sharing one dynamical $U(1)$ gauge field. The central result is a well-defined, regularised, gauge-invariant lattice axial charge $Q_{G}^{A}$ whose continuum non-conservation $\langle dQ_{G}^{A}/dt\rangle = -(2g/\pi)\!\int\! dx\,\langle E(x)\rangle$ arises as a direct dynamical consequence of minimal gauge coupling. A particle-hole transformation on the $\chi$ flavour exposes a hidden $U_{L}(2)\times U_{R}(2)$ chiral symmetry; non-Abelian bosonisation then identifies the model with a massive abelian Schwinger sector tensored with the level-$1$ $SU(2)$ Wess–Zumino–Witten model. Finally, we show that embedding the flavoured fermions in a ribbon-shaped $(2{+}1)$D Bernevig–Hughes–Zhang topological insulator and gauging the bulk in a constant background field factorises the boundary theory into two decoupled single-flavour Schwinger models, one on each edge, identifying the lattice factor of $2$ as one quantum of Schwinger anomaly per edge.

Tensor-network simulations of quantum many-body dynamics are fundamentally limited by entanglement build-up, which leads to exponentially growing computational costs. Furthermore, these classical simulation algorithms are inherently sequential as typically a tensor network representation of the quantum state is updated incrementally at each time step. We build on recently introduced randomized quantum algorithms for time evolution (TE-PAI), and adapt them to the classical simulation context with the purpose of enabling massive parallelisation. Our MPS TE-PAI approach achieves exact time evolution on average (unbiased estimator) and proceeds by representing an ensemble of randomized shallow Trotter-variant circuits as tensor networks. As each circuit instance yields a deterministic quantum state (or observable expecation value), the only source of randomness is the sampling of circuit variants; the absence of shot noise therefore yields a reduced estimator variance relative to quantum hardware implementations of TE-PAI. We simulate representative disordered one-dimensional spin-ring Hamiltonians, and numerically observe reductions in the per-sample gate-count by a factor of up to $10^3$ relative to Trotterized MPS evolution, yielding orders of magnitude reduction in the time-to-solution under realistic levels of parallelisation. Finally, we numerically observe that MPS TE-PAI is substantially more robust against severe bond-dimension truncation than product formulas, potentially making it useful for the simulation of strongly correlated systems where truncation is necessary in practice. We also demonstrate that the approach can be used naturally in combination with existing time evolution algorithms, effectively extending their time depth via parallelisation.