Papers

We define a nonlinear Fourier transform which maps sequences of contractive $n \times n$ matrices to $SU(2n)$-valued functions on the circle $\mathbb{T}$. We characterize the image of finitely supported sequences and square-summable sequences on the half-line, and construct an inverse for $SU(2n)$-valued functions whose diagonal $n \times n$ blocks are outer matrix functions. As an application, we relate this nonlinear Fourier transform with quantum signal processing over $U(2n)$ and multivariate quantum signal processing.

We report a cavity-enabled solid-state magnetometer based on an NV ensemble coupled with a dielectric cavity, achieving 12 pT/$\sqrt{\rm{Hz}}$ sensitivity and a nearly threefold gain from multispectral features. The features originate from cavity-induced splitting of the NV hyperfine levels and leverages robust quantum coherence in the doubly dressed states of the system to achieve high sensitivity. We project simulated near-term sensitivities approaching 100 fT/$\sqrt{\rm{Hz}}$, close to the Johnson-Nyquist limit. Our results establish frequency multiplexing as a new operational paradigm, offering a robust and scalable quantum resource for metrology under ambient conditions.

Quantum computing (QC) has the potential to revolutionise the future of scientific simulations. To harness the capabilities that QC offers, we can integrate it into hybrid quantum-classical simulations, which can boost the capabilities of supercomputing by leveraging quantum modules that offer speedups over classical counterparts. One example is quantum Monte Carlo integration, which is theorised to achieve a quadratic speedup over classical Monte Carlo, making it suitable for high-energy physics, strong-field QED, and multiple scientific and industrial applications. In this paper, we demonstrate that quantum Monte Carlo can be used to predict the number of pairs created when two photon beams collide head-on, a problem relevant to high-energy physics and intense laser-matter interactions. The results from the quantum simulations demonstrate high accuracy relative to theoretical predictions. The accuracy of the simulations is only constrained by the approximations required to embed polynomials and to initialise the quantum state. We also demonstrate that our algorithm can be used in current quantum hardware, providing up to 90 % accuracy relative to theoretical predictions. Furthermore, we propose pathways towards integrations with classical simulation codes.

Conveyor-mode shuttling in gated Si/SiGe devices enables adiabatic transfer of single electrons, electron patterns and spin qubits confined in quantum dots across several microns with a scalable number of signal lines. To realize their full potential, linear shuttle lanes must connect into a two-dimensional grid with controllable routing. We introduce a T-junction device linking two independently driven shuttle lanes. Electron routing across the junction requires no extra control lines beyond the four channels per conveyor belt. We measure an inter-lane charge transfer fidelity of $F = 100.0000000^{+0}_{-9\times 10^{-7}}\,\%$ at an instantaneous electron velocity of $270\,\mathrm{mm}\,\mathrm{s}^{-1}$. The filling of 54 quantum dots is controlled by simple atomic pulses, allowing us to swap electron patterns, laying the groundwork for a native spin-qubit SWAP gate. This T-junction establishes a path towards scalable, two-dimensional quantum computing architectures with flexible spin qubit routing for quantum error correction.

The potential application of topological superconductivity in quantum transport and quantum information has fueled an intense investigation of hybrid materials with emergent electronic properties, including magnet-superconductor heterostructures. Here, we report evidence of a topological nodal-point superconducting phase in a one-atom-thick in-plane ferromagnet in direct proximity to a conventional $s$-wave superconductor. Low-temperature scanning tunneling spectroscopy data reveal the presence of a double-peak low-energy feature in the local density of states of the hybrid system, which is rationalized via model calculations to be an emergent topological nodal-point superconducting phase with tilted Weyl cones. Our results further establish the combination of in-plane ferromagnetism and conventional superconductivity as a route to design two-dimensional topological quantum phases.

Scaling superconducting quantum computers to the fault-tolerant regime calls for a commensurate scaling of the classical control and readout stack. Today's systems largely rely on room-temperature, rack-based instrumentation connected to dilution-refrigerator cryostats through many coaxial cables. Looking ahead, superconducting fault-tolerant quantum computers (FTQCs) will likely adopt a heterogeneous quantum-classical architecture that places selected electronics at cryogenic stages – for example, cryo-CMOS at 4~K and superconducting digital logic at 4~K and/or mK stages – to curb wiring and thermal-load overheads. This review distills key requirements, surveys representative room-temperature and cryogenic approaches, and provides a transparent first-order accounting framework for cryoelectronics. Using an RSA-2048-scale benchmark as a concrete reference point, we illustrate how scaling targets motivate constraints on multiplexing and stage-wise cryogenic power, and discuss implications for functional partitioning across room-temperature electronics, cryo-CMOS, and superconducting logic.

Both tensor trains (TTs) and quantum states provide compressed representations of grid-structured data with potentially exponential compression power. We present a unified framework for upsampling data encoded in vector amplitudes, with efficient realizations in both classical TT and quantum settings. Starting from an \(n\)-core TT or an \(n\)-qubit state on a coarse grid with \(2^n\) points, the construction produces an \((n+m)\)-core TT or \((n+m)\)-qubit state on a finer grid with \(2^{n+m}\) points. In the TT setting, it supports interpolation, quasi-interpolation, augmentation, and synthesis through efficient low-rank contractions, with the added \(m\) cores retaining constant rank. For function-value encodings, the resulting interpolation satisfies an \(\ell^2\)-error bound independent of the number of added grid points, achieves exponential compression at fixed accuracy, and has a logarithmic complexity in the number of grid points. In the quantum setting, the refined state is prepared by a \(\mathrm{poly}(n,m)\)-size circuit using \(\log(p+1)\) ancillas, where \(p\) controls the smoothness of the quasi-interpolant; the corresponding error scales quadratically with the initial grid spacing. We validate our framework for tensor networks in one-, two-, and three-dimensional examples, including functions, derivatives, airfoil masks, and synthetic random fields such as three-dimensional turbulence. In particular, fractal fields can be generated directly in TT format with logarithmic memory and runtime. These results open a practical route to multiscale solvers, generative models, and geometry-aware algorithms on tensor-network and quantum platforms, with potential applications in scientific simulation, imaging, and real-time graphics.

We present MPM-QIR, a variational-quantum-circuit (VQC) framework for classical image compression and representation whose core objective is to achieve equal or better reconstruction quality at a lower Parameter Compression Ratio (PCR). The method aligns a generative VQC's measurement-probability distribution with normalized pixel intensities and learns positional information implicitly via an ordered mapping to the flattened pixel array, thus eliminating explicit coordinate qubits and tying compression efficiency directly to circuit (ansatz) complexity. A bidirectional convolutional architecture induces long-range entanglement at shallow depth, capturing global image correlations with fewer parameters. Under a unified protocol, the approach attains PSNR $\geq$ 30 dB with lower PCR across benchmarks: MNIST 31.80 dB / SSIM 0.81 at PCR 0.69, Fashion-MNIST 31.30 dB / 0.91 at PCR 0.83, and CIFAR-10 31.56 dB / 0.97 at PCR 0.84. Overall, this compression-first design improves parameter efficiency, validates VQCs as direct and effective generative models for classical image compression, and is amenable to two-stage pipelines with classical codecs and to extensions beyond 2D imagery.

We investigate the $(2+1)$-dimensional $q$-deformed $\mathrm{SU}(N)_k$ Yang-Mills theory in the lattice Hamiltonian formalism, which is characterized by three parameters: the number of colors $N$, the coupling constant $g$, and the level $k$. By treating these as tunable parameters, we explore how key properties of the theory, such as confinement and topological order, emerge in different regimes. Employing a variational mean-field analysis that interpolates between the strong- and weak-coupling regimes, we determine the large-$N$ phase structure in terms of the 't Hooft coupling $\lambda_\mathrm{tH}=g^2N$ and the ratio $k/N$. We find that the topologically ordered phase remains robust at large $N$ under appropriate scalings of these parameters. This result indicates that the continuum limit of large-$N$ gauge theory may be more intricate than naively expected, and motivates studies beyond the mean-field theory, both to achieve a further understanding of confinement in gauge theories and to guide quantum simulations of large-$N$ gauge theories.

Grover's algorithm is a cornerstone of quantum search algorithm, offering quadratic speedup for unstructured problems. However, limited qubit counts and noise in today's noisy intermediate-scale quantum (NISQ) devices hinder large-scale hardware validation, making efficient classical simulation essential for algorithm development and hardware assessment. We present an iterative Grover simulation framework based on matrix product states (MPS) to efficiently simulate large-scale Grover's algorithm. Within the NVIDIA CUDA-Q environment, we compare iterative and common (non-iterative) Grover's circuits across statevector and MPS backends. On the MPS backend at 29 qubits, the iterative Grover's circuit runs about 15x faster than the common (non-iterative) Grover's circuit, and about 3-4x faster than the statevector backend. In sampling experiments, Grover's circuits demonstrate strong low-shot stability: as the qubit number increases beyond 13, a single-shot measurement still closely mirrors the results from 4,096 shots, indicating reliable estimates with minimal sampling and significant potential to cut measurement costs. Overall, an iterative MPS design delivers speed and scalability for Grover's circuit simulation, enabling practical large-scale implementations.

We compare three proof techniques for composable finite-size security of quantum key distribution under collective attacks, with emphasis on how the resulting secret-key rates behave at practically relevant block lengths. As a benchmark, we consider the BB84 protocol and evaluate finite-size key-rate estimates obtained from entropic uncertainty relations (EUR), from the asymptotic equipartition property (AEP), and from a direct finite-block analysis based on the conditional min-entropy, which we refer to as the finite-size min-entropy (FME) approach. For BB84 we show that the EUR-based bound provides the most favorable performance across the considered parameter range, while the AEP bound is asymptotically tight but can become overly pessimistic at moderate and small block sizes, where it may fail to certify a positive key. The FME approach remains effective in this small-block regime, yielding nonzero rates in situations where the AEP estimate vanishes, although it is not asymptotically optimal for BB84. These results motivate the use of FME-type analyses for continuous-variable protocols in settings where tight EUR-based bounds are unavailable, notably for coherent-state schemes where current finite-size analyses typically rely on AEP-style corrections.

Topological quantum sensing leverages unique topological features to suppress noise and improve the precision of parameter estimation, emerging as a promising tool in both fundamental research and practical application. In this Letter, we propose a sensing protocol that exploits the dynamics of topological quantum walks incorporating localized defects. Unlike conventional schemes that rely on topological protection to suppress disorder and defects, our protocol harnesses the evolution time as a resource to enable precise estimation of the defect parameter. By utilizing topologically nontrivial properties of the quantum walks, the sensing precision can approach the Heisenberg limit. We further demonstrate the performance and robustness of the protocol through Bayesian estimation. Our results show that this approach maintains high precision over a broad range of parameters and exhibits strong robustness against disorder, offering a practical pathway for topologically enhanced quantum metrology.

Loophole-free quantum nonlocality often demands experiments with high complexity (defined by all parties' settings and outcomes) and multiple efficient detectors. Here, we identify the fundamental efficiency and complexity thresholds for quantum steering using two-qubit entangled states. Remarkably, it requires only one photon detector on the untrusted side, with efficiency $\epsilon > 1/X$, where $X \geq 2$ is the number of settings on that side. This threshold applies to all pure entangled states, in contrast to analogous Bell-nonlocality tests, which require almost unentangled states to be loss-tolerant. We confirm these predictions in a minimal-complexity ($X = 2$ for the untrusted party and a single three-outcome measurement for the trusted party), detection-loophole-free photonic experiment with $\epsilon = (51.6 \pm 0.4)\% $.

In this paper, we present a reproducible benchmarking framework that systematically compares QML models with architecture-matched classical counterparts across three financial tasks: (i) directional return prediction on U.S. and Turkish equities, (ii) live-trading simulation with Quantum LSTMs versus classical LSTMs on the S&P 500, and (iii) realized volatility forecasting using Quantum Support Vector Regression. By standardizing data splits, features, and evaluation metrics, our study provides a fair assessment of when current-generation QML models can match or exceed classical methods. Our results reveal that quantum approaches show performance gains when data structure and circuit design are well aligned. In directional classification, hybrid quantum neural networks surpass the parameter-matched ANN by +3.8 AUC and +3.4 accuracy points on AAPL stock and by +4.9 AUC and +3.6 accuracy points on Turkish stock KCHOL. In live trading, the QLSTM achieves higher risk-adjusted returns in two of four S&P~500 regimes. For volatility forecasting, an angle-encoded QSVR attains the lowest QLIKE on KCHOL and remains within $\sim$0.02-0.04 QLIKE of the best classical kernels on S&P~500 and AAPL. Our benchmarking framework clearly identifies the scenarios where current QML architectures offer tangible improvements and where established classical methods continue to dominate.

Signal amplification by reversible exchange (SABRE) is a nuclear spin hyperpolarization technique in which the transient interaction of parahydrogen (pH2) and a target substrate with an iridium complex leads to polarization transfer to the substrate. Here, we use a parahydrogen-enhanced, spin-selective NMR method to investigate pyruvate binding, which is combined with exchange-model fitting and DFT calculations. Our study reveals several key findings that reshape the current understanding of SABRE: (a) intramolecular hydrogen exchange of the hydrides, occurring faster than pyruvate or H2 loss; (b) the discovery of a novel stable [Ir(H)2(IMes)(ąppa}1-pyr)(DMSO)2] complex; and (c) the potential role of counterions (here Na+) in Ir-pyruvate binding. Previously unknown insights into complex kinetics and distributions as a function of temperature, [DMSO], [pyruvate], and hydrogen pressure are presented. The methods demonstrated here, exemplified by SABRE, provide a framework that is expected to guide future research in the field.

We investigate the computational hardness of estimating the quantum $\alpha$-Rényi entropy ${\rm S}^{\tt R}_{\alpha}(\rho) = \frac{\ln {\rm Tr}(\rho^\alpha)}{1-\alpha}$ and the quantum $q$-Tsallis entropy ${\rm S}^{\tt T}_q(\rho) = \frac{1-{\rm Tr}(\rho^q)}{q-1}$, both of which converge to the von Neumann entropy as the order approaches $1$. The promise problems Quantum $\alpha$-Rényi Entropy Approximation (RényiQEA$_\alpha$) and Quantum $q$-Tsallis Entropy Approximation (TsallisQEA$_q$) ask whether $ {\rm S}^ {\tt R}_{\alpha}(\rho)$ or ${\rm S}^{\tt T}_q(\rho)$, is at least $\tau_1$ or at most $\tau_2$, where $\tau_1 - \tau_2$ is typically a positive constant. Previous hardness results cover only the von Neumann entropy (order $1$) and some cases of the quantum $q$-Tsallis entropy, while existing approaches do not readily extend to other orders. We establish that for all positive real $\alpha$ and $q$, and also for $\alpha=\infty$, the rank-$2$ variants Rank2RényiQEA$_\alpha$ and Rank2TsallisQEA$_q$ are BQP-hard. Combined with prior (rank-dependent) quantum query algorithms in Wang, Guan, Liu, Zhang, and Ying (TIT 2024), Wang, Zhang, and Li (TIT 2024), and Liu and Wang (SODA 2025), as well as the one derived from O'Donnell and Wright (STOC 2016), our results imply: - For all real orders $\alpha > 0$ or $\alpha=\infty$, and for all real orders $0 < q \leq 1$, LowRankRényiQEA$_\alpha$ and LowRankTsallisQEA$_q$ are BQP-complete, where both are restricted versions of RényiQEA$_\alpha$ and TsallisQEA$_q$ with $\rho$ of polynomial rank. - For all real order $q>1$, TsallisQEA$_q$ is BQP-complete. Our hardness results stem from reductions based on new inequalities relating the $\alpha$-Rényi or $q$-Tsallis binary entropies of different orders. These reductions differ substantially from previous approaches, and the inequalities are of independent interest.

We investigate the combined effects of spatial curvature and topology on the properties of the vacuum state for a charged scalar field localized on the (2+1)-dimensional Beltrami pseudosphere, assuming that the field obeys quasiperiodicity condition with constant phase. As important local characteristics of the vacuum state the vacuum expectation values (VEVs) of the field squared and energy-momentum tensor are evaluated. The contributions in the VEVs coming from geometry with an uncompactified azimuthal coordinate are divergent, whereas the compact counterparts are finite and are analysed both numerically and asymptotically. For small values of proper radius of the compactified dimension, the leading terms of topological contributions are independent of the field mass and curvature coupling parameter, increasing by a power-law. In the opposite limit, the VEVs decay following a power-law in the general case. In the special case of a conformally coupled massless field the behavior is different. Unlike the VEV of field squared and vacuum energy density, the radial and azimuthal stresses are increasing by absolute value. As a consequence, the effects of nontrivial topology are strong for the stresses in this case at small values of radial coordinate.

The coherent mobility of doublons, arising from second-order virtual dissociation-recombination processes, fundamentally limits their use as information carriers in the strongly interacting Bose-Hubbard model. We propose a disorder-free suppression mechanism by introducing an optimized nearest-neighbor pair-hopping term that destructively interferes with the dominant virtual hopping channel. Using the third-order Schrieffer-Wolff transformation, we derive an analytical optimal condition that accounts for lattice geometry corrections. Exact numerical simulations demonstrate that this optimized scheme achieves near-complete dynamical arrest and entanglement preservation in one-dimensional chains, while in two-dimensional square lattices, it significantly suppresses ballistic spreading yet permits a slow residual expansion. Furthermore, in the many-body regime, finite-size scaling analysis identifies the observed long-lived density-wave order as a prethermal plateau emerging from the dramatic separation of microscopic and thermalization timescales.

We investigate mixed-state entanglement measures, namely reflected entropy, mutual information and logarithmic negativity, for two disjoint intervals in one-dimensional systems excited by a finite number of quasiparticles. While whole system is in a pure state, the two disjoint intervals are in a generically mixed state. To address the problem that natural subsystem bases are generically non-orthonormal in such excited states, we use a general and efficient algorithm that computes these measures directly from the density matrix expressed in an arbitrary non-orthonormal basis. Applying this method to classical, bosonic, and fermionic quasiparticle excitations on a circle, we discover a universal additivity property: in the limit of large momentum differences, the mixed-state entanglement of a multi-quasiparticle state decomposes exactly into the sum of independent contributions. This additivity unifies the entanglement behavior across classical and quantum statistics, with the classical result emerging naturally as a special case. Our findings establish a robust computational framework for mixed-state entanglement in excited many-body systems and reveal a generic decoupling mechanism that governs entanglement distribution beyond the ground state.

Non-Hermitian skin effects are conventionally manifested as boundary localization of eigenstates in linear systems. In nonlinear settings, however, where eigenstates are no longer well defined, it becomes unclear how skin effects should be faithfully characterized. Here, we propose a Koopman-based characterization of nonlinear skin effects, in which localization is defined in terms of Koopman eigenfunctions in a lifted observable space, rather than physical states. Using a minimal nonlinear extension of the Hatano-Nelson model, we show that dominant Koopman eigenfunctions localize sharply on higher-order observables, in stark contrast to linear skin effects confined to linear observables. This lifted-space localization governs the sensitivity to boundary amplitude perturbations, providing a distinct dynamical signature of the nonlinear skin effect. Our results establish the Koopman framework as a natural setting in which skin effects unique to nonlinear non-Hermitian systems can be identified.

Bias-tailored codes such as the XZZX surface code and the domain wall color code achieve high dephasing-biased thresholds because, in the infinite-bias limit, their $Z$ syndromes decouple into one-dimensional repetition-like chains; the $X^3Z^3$ Floquet code shows an analogous strip-wise structure for detector events in spacetime. We capture this common mechanism by defining strip-symmetric biased codes, a class of static stabilizer and dynamical (Floquet) codes for which, under pure dephasing and perfect measurements, each elementary $Z$ fault is confined to a strip and the Z-detector–fault incidence matrix is block diagonal. For such codes the Z-detector hypergraph decomposes into independent strip components and maximum-likelihood $Z$ decoding factorizes across strips, yielding complexity savings for matching-based decoders. We characterize strip symmetry via per-strip stabilizer products, viewed as a $\mathbb{Z}_2$ 1-form symmetry, place XZZX, the domain wall color code, and $X^3Z^3$ in this framework, and introduce synthetic strip-symmetric detector models and domain-wise Clifford constructions that serve as design tools for new bias-tailored Floquet codes.

We present a quantum algorithm for simulating dissipative diffusion dynamics generated by positive semidefinite operators of the form $A=L^\dagger L$, a structure that arises naturally in standard discretizations of elliptic operators. Our main tool is the Kannai transform, which represents the diffusion semigroup $e^{-AT}$, where $T$ is the final simulation time, as a Gaussian-weighted superposition of unitary wave propagators. For target accuracy $\varepsilon$, this representation leads to a linear-combination-of-unitaries implementation with a Gaussian tail and yields query complexity $\widetilde{O}(\sqrt{\|A\|\,T\,\log(1/\varepsilon)})$, up to the standard dependence on state-preparation and output-norm factors, improving the scaling in $\|A\|$, $T$, and $\varepsilon$ compared with generic Hamiltonian-simulation-based methods. We instantiate the method for the heat equation and biharmonic diffusion under non-periodic physical boundary conditions, and further use it as a subroutine for constant-coefficient linear parabolic surrogates arising in entropy-penalization schemes for the viscous Hamilton–Jacobi equations. In the long-time regime, under a spectral-gap assumption, the same framework gives a structured quantum linear solver by exploiting convergence to the steady state. For normalized positive definite systems $\sigma(A)\subset[1,\kappa]$, the solver outputs an $\varepsilon$-approximation to the state proportional to $\mathbf{x}=A^{-1}\mathbf{b}$ with query complexity $\widetilde{O}\left(\frac{\|\mathbf{b}\|}{\|\mathbf{x}\|}\sqrt{\kappa}\log^2\frac{\|\mathbf{b}\|}{\varepsilon\|\mathbf{x}\|}\right)$.

The Rashba, Dresselhaus, and Weyl Hamiltonians form a foundational framework for modeling spin-orbit interactions across condensed matter systems. Although they describe distinct material classes and produce seemingly different spin textures, they are conventionally treated as separate, unrelated theoretical frameworks. Here, this work demonstrates that the linear 2D Rashba and Weyl models are connected by a specific unitary transformation that maps one Hamiltonian exactly onto the other. The same unitary can be applied to map the linear Dresselhaus-1 model onto the Dresselhaus-2 models and vice versa. Such hidden correspondence establishes a unified theoretical foundation for spin-orbit interactions, deepening our conceptual understanding of spin-orbit coupling and opening new avenues for exploring complex spin textures. To illustrate the application, this work introduces a unique, improved, and more realistic model Hamiltonian H_MKM combining all known foundational spintronic models, where the stringent condition of equal spin-orbit coupling strength of Rashba and Dresselhaus may not be required to observe persistent spin texture under MKM transformation.

We show that Weyl's abandoned idea of local scale invariance has a natural realization at the quantum level in pilot-wave (de Broglie-Bohm) theory. We obtain the Weyl covariant derivative by complexifying the electromagnetic gauge coupling parameter. The resultant non-hermiticity has a natural interpretation in terms of local scale invariance in pilot-wave theory. The conserved current density is modified from $|\psi|^2$ to the local scale invariant, trajectory-dependent ratio $|\psi|^2/ \mathbf 1^2[\mathcal C]$, where $\mathbf 1[\mathcal C]$ is a scale factor that depends on the pilot-wave trajectory $\mathcal C$ in configuration space. All physical predictions are local scale invariant, even in the presence of mass terms. Our approach is general, and we implement it for the Schrödinger and Pauli equations, and for the Dirac equation in curved spacetime, each coupled to an external electromagnetic field. We also implement it in quantum field theory for the case of a quantized axion field interacting with a quantized electromagnetic field. We discuss the equilibrium probability density and show that the corresponding trajectories are unique. Our results provide a pivotal understanding of local scale invariance in quantum theory.

We propose a protocol to tailor dynamical quantum phase transitions (DQPTs) by double-mode squeezing onto the initial state in the XY chain. The effect of squeezing depends critically on the system's symmetry and parameters. When the squeezing operator breaks particle-hole symmetry (PHS), DQPTs become highly tunable, allowing one to either induce transitions within a single phase or suppress them. Remarkably, when PHS is preserved and the squeezing strength reaches $r=\pi/4$, a universal class of DQPTs emerges, independent of the quench path. This universality is characterized by two key features: (i) the collapse of all Fisher zeros onto the real-time axis, and (ii) the saturation of intermode entanglement to its maximum in each $(k,-k)$ modes. Moreover, the critical momenta governing the DQPTs coincide exactly with the modes attaining the maximal entanglement. At this universal point, the dynamical phase vanishes, leading to a purely geometric evolution marked by $\pi$-jumps in the Pancharatnam geometric phase. Our work establishes initial-state squeezing as a versatile tool for tailoring far-from-equilibrium criticality and reveals a direct link between entanglement saturation and universal nonanalytic dynamics.

We propose the Many-Body Quantum Score (MBQS), a practical and scalable application-level benchmark protocol designed to evaluate how accurately quantum processing units (QPUs)—both gate-based and analog—can simulate many-body quantum dynamics. MBQS quantifies performance by identifying the maximum number of qubits with which a QPU can reliably reproduce correlation functions of the transverse-field Ising model following a specific quantum quench. This quench reveals a global correlation surge, a collective peak of connected correlations at all distances, which to our knowledge had not been noticed so far in this context. This paper presents the MBQS protocol and highlights its design principles, supported by analytical insights, classical simulations, and experimental data. It also reports an implementation on Ruby, an analog QPU based on Rydberg atoms developed by the Pasqal company. These findings demonstrate MBQS's potential as a robust and informative tool for benchmarking near-term quantum devices for many-body physics. While the present experimental validation has been performed on an analog device, the protocol has the potential to be implemented on digital platforms as well.

We present an optical study of near-surface Er$^{3+}$ ensembles in waveguide-integrated CaWO$_4$ and YVO$_4$, investigating how nanophotonic coupling modifies rare-earth spectroscopy. In particular, we compare bulk excitation with evanescently coupled TE and TM waveguide modes. In Er$^{3+}$:CaWO$_4$, we observe a pronounced polarization-dependent surface effect. TE-coupled spectra closely reproduce bulk behavior. In contrast, TM coupling induces strong inhomogeneous broadening and an asymmetric low-energy shoulder of the site S1 Y1Z1 transition, with linewidths exceeding those of the bulk by more than a factor of four. Temperature-dependent measurements and surface termination studies indicate that surface charges are the dominant mechanism. Er$^{3+}$:YVO$_4$ remains largely unaffected by mode polarization, and surface termination leads only to minor spectral shifts. These observations suggest that non-charge-neutral rare-earth systems are more susceptible to surface-induced decoherence sources than charge-neutral hosts.

Many symmetry protected or symmetry enriched phases of quantum matter have the property that every ground state in a given such phase endows measurement based quantum computation with the same computational power. Such phases are called computational phases of quantum matter. Here, we experimentally verify four theoretical predictions for them on an IBM superconducting quantum device. We comprehensively investigate how symmetric imperfections of the resource states translate into logical decoherence, and how this decoherence is mitigated. In particular, the central experiment probes the scaling law from which the uniformity of computational power follows. We also analyze the correlated regime, where local measurements give rise to logical operations collectively. We test the prediction that densest packing of a measurement-based algorithms remains the most efficient, in spite of the correlations. Our experiments corroborate the operational stability of measurement based quantum computation in quantum phases of matter with symmetry.

We introduce an implementation of a non-degenerate parametric amplifier in which the signal and idler modes, respectively, a microwave mode and an ensemble of spins (e.g., nitrogen-vacancy centers in diamond), are operated in their linear regime. This paramp, which amplifies signals in both parts at room and cryogenic temperatures, can be used to generate both the two-mode and single-mode squeezing of either system. It requires merely modulating the frequency of the spin ensemble at the sum of the cavity and spin frequencies (providing the classical pump) with the two systems sufficiently detuned. This effect is remarkable given that modulating a spin ensemble by itself produces neither amplification nor squeezing, unlike modulating an oscillator, and that an off-resonant perturbative analysis would suggest that modulating the spin ensemble merely parametrically drives the cavity mode. With typical cavity parameters including a cavity quality factor~$Q=10^4$, and a 1 GHz modulation amplitude, the microwave signal can be amplified by approximately $18~\mbox{dB}$ in $1.7~\mbox{$\mu$s}$, with a resonant bandwidth of about $0.5~\mbox{MHz}$. At $10~\mbox{mK}$ with the same modulation amplitude and a cavity and spin $Q=5\times 10^4$ it generates approximately $5~\mbox{dB}$ of squeezing. We also examine the experimental requirements for implementation.

We examine multiphoton interference with a symmetric $SU(N)$ beam splitter $S_N$, an extension of features of the $SU(2)$ 50/50 beam splitter extended Hong-Ou-Mandel (eHOM) effect, whereby one obtains a zero amplitude (probability) for the output coincidence state (defined by equal number of photons $n/N$ in each output port), when a total number $n$ of photons impinges on the $N$-port device. These are transitions of the form $|n_1,n_2,\ldots,n_N\rangle\overset{S_N}{\to}|n/N\rangle^{\otimes N}$, where $n=\sum_{i=1}^N n_i$, which generalize the Hong-Ou-Mandel (HOM) effect $|1,1\rangle \overset{S_2}{\to}|1,1\rangle $, the eHOM effect $|n_1,n_2\rangle \overset{S_2}{\to}|\tfrac{n_1+n_2}{2},\tfrac{n_1+n_2}{2}\rangle $, and the generalized HOM effect (gHOM) $|1\rangle^{\otimes N}\overset{S_N}{\to}|1\rangle^{\otimes N}$, which have previously been studied in the literature. The emphasis of this work is on illuminating how the overall destructive interference occurs in separate groups of destructive interferences of sub-amplitudes of the total zero amplitude. We develop symmetry properties for the generalized eHOM effect (geHOM) $|n_1,n_2,\ldots,n_N\rangle\overset{S_N}{\to}|n/N\rangle^{\otimes N}$ involving a zero amplitude governed by Perm($\Lambda$)=0, for an appropriately constructed matrix $\Lambda(S_N)$ built from the matrix elements of $S_N$. We develop an analytical constraint equation for Perm$(\Lambda)$ for arbitrary $N$ that allows us to determine when it is zero. We generalize the SU(2) beam splitter feature of central nodal line (CNL), which has a zero diagonal along the output probability distribution when one of the input states is of odd parity (containing only odd number of photons), to the general case of $N = 2 * N'$ where $N'\in odd$.