Papers

Small particles respond to light through electric and magnetic patterns called multipoles. Their interference determines the direction of scattered light and affects the decay of nearby emitters. Larger nanoantenna arrays, metasurfaces, and other nanophotonic systems demand increasingly costly simulations to meet strict performance and reliability requirements. We show how to represent these responses with quantum bits (qubits) and use quantum circuits to calculate scattering and the Purcell effect, the change in an emitter's decay rate caused by its surroundings. The circuits can encode known multipoles or estimate unknown ones from Maxwell's equations at the particle surface. We test the method against classical theory using an ideal quantum simulator and real IBM and IonQ quantum computers. For a silicon nanoparticle of radius 78 nm illuminated at 679 nm, IonQ gives a forward directivity of 2.984, close to the ideal dipolar value 3, quantifying the concentration of scattered light in the forward direction. For an emitter 10 nm from a gold nanoparticle and emitting at 510 nm, IBM measurements yield a calculated total decay rate 20.48 times its vacuum value versus 20.47 from theory, with correction for errors in reading the qubit states. We also establish how the number of multipoles and the required accuracy set the demands on qubits and measurements. These results motivate quantum-assisted design studies of directional radiation and emitter decay in nanophotonic systems.

Topological insulators lie outside the standard Landau paradigm: Instead of being characterized by local order parameters and spontaneous symmetry breaking, their defining properties are encoded in the global structure of their quantum states and manifested as robust boundary states. Crystalline higher-order topological insulators can extend this principle by hosting localized corner modes in two and higher dimensions. Although these boundary modes can be detected spectroscopically through the local density of states, determining the underlying bulk topology, for example through nested Wilson loops, remains an outstanding experimental challenge. Here we use up to 728 qubits in a superconducting analogue-digital quantum simulator to realize and probe a hierarchy of one-, two-, and three-dimensional topological insulators governed by effective tight-binding Hamiltonians. By reconstructing complex propagators from real-space measurements, we obtain complementary boundary and bulk characterizations: local densities of states reveal localized boundary modes, whereas experimentally reconstructed band projectors yield bulk topological invariants, culminating in a three-dimensional octupole moment extracted through third nested Wilson loops. These results establish a route to realizing, characterizing, and manipulating topologically protected states in programmable three-dimensional quantum matter.

Orbital optimization within the variational quantum eigensolver (VQE) framework is useful for improving active-space quality and accurately capturing electron correlation with shallow quantum circuits. The orbital-optimized unitary coupled cluster (OO-UCC) method [W. Mizukami et al., Phys. Rev. Res. 2, 033421 (2020)], in its original formulation, employs the orbital gradient and Hessian in a Newton–Laphson procedure. In this study, we propose sequential orbital-optimized unitary coupled cluster doubles (SOO-UCCD), extending the gradient-free sequential minimal optimization approach [K. M. Nakanishi et al., Phys. Rev. Res. 2, 043158 (2020)] to orbital rotation parameters. Based on the frequency structure implied by the algebraic properties of the orbital rotation generators, our method determines each orbital rotation angle from a nine-point energy profile. The energy landscape is evaluated by classically transforming the Hamiltonian integrals and contracting them with the reduced density matrices already measured for the VQE energy evaluation. This allows for orbital optimization without evaluating the orbital gradient or Hessian. By building on the core benefits of the OO-UCC framework, our method provides a gradient-free and hyperparameter-free orbital-optimization procedure. State-vector simulations for LiH in the STO-3G basis show that SOO-UCCD reproduces UCCSD-level accuracy along the bond-stretching curve and exhibits macro-iteration convergence comparable to that of L-BFGS-based orbital optimization.

Entanglement witnesses are a key instrument in determining the entanglement of unknown bipartite qubit systems without requiring a full tomography. Riccardi et al. (2020) proposed a set of two-qubit witnesses, denoted $\{W^3\}$, which require a set of local measurements in three bases. In this work we propose an additional set of witnesses, denoted $\{W^5\}$, that build on $\{W^3\}$ and require measurements in two additional bases. We present an adaptive two-step method that uses the measurement results obtained for $\{W^3\}$ to select which $W^5$ to evaluate. We computationally demonstrate that our adaptive two-step protocol can identify 72% of entangled states generated according to the methods of Roik et al. (2021) while requiring measurements in only five of the nine bases needed for full quantum tomography. Additionally, we experimentally verify the behavior of both sets of witnesses using polarization-entangled photon pairs.

Entanglement entropies of multipartite pure states satisfy linear inequalities; for stabilizer states, the backbone of quantum error correction, the achievable entropy vectors generate a convex cone. We show that, after a linear lift, the closure of this cone is exactly the self-dual slice of a classical object on one more element than there are parties: the closed cone generated by the rank functions of arrangements of subspaces over $\mathrm{GF}(d)$ for a prime local dimension $d$, and of subgroups of abelian groups of exponent dividing $d$ in general. This correspondence brings two features of rank inequalities into play. First, the purifier acts as an extra variable: already for five parties the stabilizer cone is strictly smaller than the quantum linear rank cone, with which earlier evidence suggested it might coincide; an explicit six-variable inequality witnesses the gap. Second, the characteristic of the field enters, so the cone depends on the local dimension. It does not for up to four parties. With enough parties, the cone of $\mathbb{Z}_d$ qudits is contained in that of $\mathbb{Z}_{d'}$ qudits if and only if $d$ divides $d'$, and already for six parties the logical zero state of Steane's code lies outside the stabilizer cone of every odd local dimension. For this last result we show that the seven-variable characteristic-dependent inequalities we checked all follow from Shannon inequalities on the self-dual slice, and instead transfer one of them through a family of subgroups derived from the Fano plane.

Fault-tolerant quantum design relies on simulating noisy quantum error correction (QEC) circuits with hundreds to thousands of physical qubits and large circuit depths. High-throughput stabilizer simulators, such as Stim, are the de facto tools for this regime. However, these simulators typically support only Clifford gates and stochastic Pauli noise. As a result, more realistic device noise models-such as coherent control errors and thermal relaxation-are often approximated using Pauli Twirling. This replacement offers speed but can introduce inaccuracies. To address this, we propose QuasiStim for scalable exact QEC simulation of realistic noise, with QuasiFold as a compiler frontend and DeltaFrame as a sampling backend. QuasiFold lowers sampling overhead by merging equivalent non-Pauli sites in the quantum circuit. DeltaFrame samples the resulting stabilizer trajectories by reusing precomputed predicted references instead of resimulating the full tableau for each shot. QuasiFold can reduce the number of shots by over $100\times$ on bivariate-bicycle (BB) codes. Compared to a per-shot tableau simulator, DeltaFrame is $166\times$–$1331\times$ faster on rotated surface-code memory under coherent $R_Z$ noise as $d$ grows from 5 to 35. QuasiStim also makes the PTA bias measurable: the gap is non-negligible and varies with noise strength, decoder, and circuit schedule. QuasiStim will be open-sourced upon publication at https://github.com/isaachyw/QuasiStim

A multiple-access channel transmits messages from several senders to one receiver. For quantum outputs, the receiver must distinguish states depending on all messages, while tests for different decoding errors may not commute. In this paper, we give a self-contained account of simultaneous decoding, in which one joint measurement recovers the messages at the established achievable rates using independent codebooks. The proof combines Sen's tilted-subspace method with Cheng's analysis of the pretty-good measurement, a standard measurement for distinguishing quantum states. Auxiliary labels close to a common reference state separate the geometric directions associated with different errors. A Gram-matrix bound then allows separately chosen binary tests to certify the performance of one measurement. The labels are removed from the completed decoder, leading to a deterministic packing theorem that applies to finite state ensembles, including unequal priors. We develop the two-sender case first, with worked examples and figures, before treating any fixed number of senders. Sharper finite-error guarantees for the same decoder are developed separately, with sharp bounds for two specified steps of the analysis rather than a claim of globally optimal decoding. Applications recover known achievable rates for unassisted and entanglement-assisted classical communication, distributed classical compression with quantum side information, and coherent-state communication through a beamsplitter or passive multiport device. The common proof makes the geometric mechanism and its quantitative consequences available without repeating the construction for each application.

Quantum decoherence does not permit the operation of (new) miniaturised devices or protocols that can otherwise run unassisted in classical settings, without explicit user agent intervention in the quantum realm. Classical driving signals are inadequate to model a truly quantum device built entirely in the quantum realm. Synchronisation endures miniaturisation without sacrificing quantum coherence, unlike thermalisation. We quantise the Kuramoto model and the Van der Pol oscillator and couple it with qubit(s) to show a robust limit cycle even in typical quantum thermodynamical settings. Novel effects such as Entanglement State Revivals, Zeno locking and thermal current interference are observed. We propose toy model devices that can be implemented harnessing the Synchronisation of a qubit (or an arbitrary collection of qubits) with a nonlinear phase oscillator. The framework theory used here, Spectral Mechanics, is a fresh look at the foundations of conventional quantum mechanics and is extensible to completely new forms of interactions and models in experimental settings.

The protection of quantum information in systems leveraging more than two basis states (qudits) is an ongoing challenge. A tool for transforming the plethora of qubit stabilizer codes into qudit stabilizer codes is the local-dimension-invariant (LDI) extension to the stabilizer formalism. A drawback of the LDI method is that generally in order to ensure the distance of the resultant code is at least preserved the newly selected local-dimension will need to be any prime bigger than a large cutoff value. In this work we consider the repercussions of this large cutoff value, then by focusing on quantum low-density parity-check (qLDPC) codes we are able to bring down this cutoff value. Notably, we are able to reduce it in some cases to the initial local-dimension and for arbitrary CSS qubit-derived codes to singly exponential in $d$ with no dependence on the number of logicals $k$. This reduces the size of the caveat in LDI forms in their distance preservation condition.

Zero-added-loss multiplexing (ZALM)~[Phys. Rev. Appl. 19, 054029 (2023)] is a midpoint-source scheme for entanglement distribution. Its quantum transmitter (QTX) employs two Sagnac sources of polarization-entangled pairs plus partial Bell-state measurements (BSMs) on their idler outputs to herald entanglement in their signal outputs, thus avoiding the switching losses suffered by other multiplexing approaches. When its Sagnac sources' spontaneous parametric downconverters (SPDCs) utilize $\chi^{(2)}$ crystals with phase-matched spectral islands, ZALM avoids the tradeoff between heralding efficiency and output-state fidelity incurred when conventional-crystal SPDCs' outputs are channelized by dense wavelength-division multiplexing filters. This paper reports an end-to-end fidelity and rate analysis for islands-based ZALM that includes: SPDC multipair emissions of all orders and partial-BSM losses in the ZALM QTX; loss in propagation from the ZALM QTX to Alice and Bob's quantum receivers (QRXs); loss and mode mismatch in the mode converters employed at Alice and Bob's QRXs; and loss in the partial-BSM loading of Alice and Bob's emissive quantum memories.

Quantum localization, the phenomenon by which wave diffusion is halted, here understood as the superposition of propagating waves, can have many different origins, such as onsite random potentials in Anderson localization. However, some systems exhibit complex geometries and could potentially reproduce similar results without external potentials. In this work, we study how localization of single particle states is affected by the underlying topology of a simple but highly predictive complex network model: the geometric soft configuration model $\mathbb S^1$. By tuning the model's temperature and heterogeneity parameters, we are able to generate different structural patterns and focus on the transition from a small-world to a non-small-world regime. Our study suggests that small-world networks tend to have extended states, which spread the whole system's size, in contrast to non-small-world networks which tend to localize, showing how complex connectivity can play the role assigned to random energetic disorder.

We numerically study spectral and wave-function properties, as well as time evolution, in finite-size commensurate approximations of the 3D Aubry-André model. We explicitly show, through scaling analysis, the existence of two critical points. One of the critical points is a localization-delocalization transition in real-space, similar to the critical point in the 3D Anderson model. The other is a localization-delocalization transition in a dual momentum space, which does not have any analog in the 3D Anderson model. We estimate their values from spectral measures, such as the gap ratio and the spectral form factor, as well as from wave-function measures, such as the inverse participation ratio and the participation entropy. Additionally, we demonstrate that both critical points exhibit hallmarks of critical points between quantum ergodicity and localization. First, the eigenstates and the eigenspectrum at the critical point fulfill the mutual universal relation between the eigenspectrum compressibility $\chi$ and the eigenstate fractal dimension $d^{(1)}$ [S. Jiricek et al., Proc. Nat. Acad. Sci. 123, e2518027123 (2026)]. Second, the dynamics of the wave packet at the critical point is governed by the critical dynamical exponent $z$, characterized by the spectral fractal dimension and the dimensionality of the lattice [M. Hopjan and L. Vidmar, Phys. Rev. Lett. 135, 060401 (2025)]. Thus, similarly to the 3D Anderson model, the critical dynamical exponent is sub-diffusive with the value of $z=3$.

Ordinal optimization selects designs that are good enough rather than optimal when evaluation is expensive and the budget is small. The alignment probability measures the chance that the selected set contains at least a prescribed number of truly good designs. We formulate quantum ordinal optimization as a class of problems that maximize this probability under a fixed oracle budget. We study exact, coherent pairwise comparisons with one selected design as a first instance, with the budget counted in comparison queries. For every classical comparison budget smaller than the number of designs, we prove that a knockout among uniformly sampled designs attains the optimum. The sample is drawn without replacement and contains one more design than the number of comparisons. For a uniformly random order, the selected design's rank has a polynomial probability law. At the same comparison budget, the maximum quantum degree is at most twice the classical maximum, a factor attained whenever there are more than twice as many designs as comparisons, so a quantum algorithm can realize rank laws that no classical algorithm with the same budget can. We determine the exact quantum optimum with one comparison for every number of designs and every good-enough set. For selecting the best of n designs, the optimal success probability with one quantum comparison is approximately 3.30/n for large n, compared with 2/n for one classical comparison. When the best design is sought, one quantum comparison gives a higher alignment probability than two classical comparisons whenever there are more than nine designs.

The Hanbury Brown-Twiss (HBT) effect marked a significant milestone in quantum optics by illustrating that intensity correlations can reveal inherently quantum features of light. Motivated by the structure of the original 8-mode HBT interferometer, we develop a simplified 4-mode toy-model that captures its essential physical ingredients while significantly reducing analytical complexity. Our toy-model provides an intuitive and well-structured framework for examining the desired correlation functions. We investigate whether our toy-model can exhibit genuinely non-classical behaviour by formulating and applying a Bell-type inequality tailored to its specific measurement configuration. Our analysis shows that the correlations generated by the toy-model violate the inequality, indicating the presence of non-local features that are incompatible with any classical hidden-variable description. The results highlight the power of simplified models in uncovering fundamental aspects of non-classicality in optical systems.

The coherent exchange of pairs of quanta between two bosonic modes is equivalent, up to a numerical factor, to the generator of two-axis twisting, which attains Heisenberg scaling of the squeezing; we determine how much of that equivalence survives in a realistic nonlinear medium. In any such medium the exchange term is accompanied by a local Kerr nonlinearity, and within a sector of fixed particle number the Hamiltonian depends on a single dimensionless ratio of the effective Kerr strength to the exchange strength. We derive two criteria governing this ratio and show that they are independent. The first is the rate at which squeezing is generated: the exponential rate follows in closed form, vanishes exactly when the two strengths are equal in magnitude, and recovers as the square root of the detuning from that point. The second is the depth to which squeezing can be driven: at fixed particle number the attainable depth is a function of the anisotropy of the mean-field saddle alone, reaches the full two-axis-twisting value only for the isotropic saddle, and is independent of the rate, parameter sets sharing the same anisotropy but differing in rate by a sizeable factor giving depths that agree to every printed digit. The closed-form Gaussian dynamics reveals a pitfall: extracting rates from short-time data overestimates them by more than a factor of two unless the saddle is isotropic. An isotropic cubic nonlinearity imposes the value of the ratio where the rate vanishes, so realizing two-axis twisting by pair exchange requires detuning it; a moderate mismatch already recovers almost all of the optimal depth. We also explain why the constraint has not surfaced before: at two quanta the self-Kerr shifts the two coupled states equally, so the two twisting Hamiltonians are rigorously indistinguishable, the degeneracy being lifted only from four quanta.

We consider singular Schrödinger operators with an attractive $\delta'$ interaction of strength $\beta$ supported by a star graph and prove that the ground-state eigenvalue is uniquely maximized by stars of the maximum symmetry. Two different proofs are presented, one for an even number of edges using a comparison to the $\delta$ interaction star, the other general based on Birman-Schwinger technique.

Levitated dielectric sensors are often treated as perturbative probes of a pre-existing optical field. We show that this picture fails for a novel sensor geometry: a dielectric plate whose transverse extent captures the full optical mode. Even an optically thin sensor with modest refractive index can heavily restructure the intracavity fields, effectively partitioning a cavity into two subcavities with substantial cross-coupling and producing new optical resonances and trapping sites with properties which deviate from the standard intuition. Conventional perturbative expressions for the cavity frequency shift and optomechanical coupling depart from the exact solution at thicknesses as small as $\sim\lambda/(10n)$. Furthermore, cavity locking decreases the coupling and increases the trap stiffness by orders of magnitude, and the coupled subcavity structure permits appreciable optical cooling at zero detuning, a distinct behavior unique to this nonperturbative regime. These effects, sharing a common origin, establish a distinct regime of cavity optomechanics in which the sensor is no longer merely a probe of the cavity field, but an element that helps define it.

Strongly interacting fermionic chains obeying a $Z_n$ symmetry can host topological phases with highly degenerate ground-state manifolds. The different ground states are connected by strongly bound fractionalized edge modes called parafermionic zero modes. While exactly preserved under time evolution at zero temperature, interactions lead to a destabilization of these modes at finite energy densities, and, given a generic Hamiltonian, obtaining the lifetime of these modes is often impossible as it requires knowledge of the full spectrum of the Hamiltonian or long time evolutions. Here, we propose a novel quantum-inspired algorithm that can be used to variationally search for these almost strong zero modes in generic parafermionic chain models. We demonstrate analytically and numerically that our ansatz faithfully reproduces the exact results and finds strongly bound parafermionic operators in the non-integrable Kitaev-Hubbard chain and chiral Potts model. Our results further allow us to verify the existence of a true strong zero mode in the symmetrical limit of the chiral Potts model for $Z_3$-parafermions.

Spatial interference can break the inversion symmetry of an effective Josephson potential even in a laterally homogeneous junction with inversion-symmetric local energy. We derive a general criterion based on harmonic-dependent complex spatial form factors: a single local harmonic remains a shifted cosine, whereas multiple harmonics can acquire phases incompatible with any common inversion center, permitting cubic interactions and three-wave mixing. Phase reflection reverses all odd-order expansion coefficients about corresponding minima while preserving the even-order ones. For a capacitively shunted junction at fixed charging energy and arbitrary offset charge, it also preserves the complete quantum spectrum, enabling reversal of the odd nonlinear sector without spectral retuning. A centered cross junction with multiple local harmonics is one concrete example. Within the imposed-field model, a three-component magnetic field can generate an inversion-asymmetric effective potential, and reversing an odd number of field components exactly reflects it without changing the quantum spectrum of the capacitively shunted junction.

In molecules, electronic coherence typically decays within a few hundred femtoseconds, a timescale usually attributed to the environment. We show that it is set to a large extent by which coherences are excited: the molecule can host coherences that the bath barely sees, and optimal control can find and exploit them. Using gradient-based optimal control (GRAPE) with Redfield dynamics, we maximize the total vibronic coherence between the electronic ground and an optically bright electronic excited state in a three-state model with nonadiabatic features, in two regimes of low and high nonadiabatic mixing and for different system-bath coupling strengths. The Liouvillian spectrum reveals, in both regimes, a well-separated coherence mode with substantial ground-bright character that decays on a picosecond timescale. This mode is protected by mixing with the dark state, which brings the geometry of the optically accessible excited state close to that of the ground state and thus weakens its coupling to the bath. Without prior information, GRAPE finds this mode and keeps the photoexcited wave packet near the Franck-Condon region. At equal fluence, it increases the coherence at 800 fs by more than two orders of magnitude over a Gaussian pulse at the weakest coupling, and by many more at stronger coupling. The Liouvillian predicts the bath dependence of the long-lived coherence to within 1%, and a Gaussian pulse tuned to the mode frequency comes within 10% of the best pulse found. Beyond two-level and spin systems, our study puts forward the vibronic structure of molecules as a route to protected coherence subspaces, identifying a dark-state protection mechanism and a design criterion based on molecular parameters. These results lay the groundwork for our goal of using vibronic molecular systems as a resource for long-lived, controllable coherences in molecular quantum information science.

The detection and quantification of quantum entanglement are central problems in quantum information science. Based on the realignment of density matrices and the vectorization of reduced density matrices, we construct a family of enhanced separability criteria for bipartite and multipartite quantum systems. Our criteria recover the previously proposed realignment-based criteria as special cases and detect more entanglement than the earlier ones. Moreover, we derive tighter lower bounds for the concurrence and the convex-roof extended negativity, criteria for detecting the genuine tripartite entanglement, and lower bounds on genuine-multipartite-entanglement concurrence. The advantages of our results are demonstrated with detailed examples.

Stationary scattering admits a dynamical formulation in which the scattering data are extracted from the evolution operator for an effective non-stationary and non-unitary quantum system. This provides a mathematical method of solving scattering problems. We consider the possibility of realizing such a quantum system and relating the scattering information to outcomes of quantum measurements performed on it. To this end, we examine stationary scattering problems in one dimension where the scattering data consist of reflection and transmission amplitudes and the corresponding effective quantum system is a non-unitary two-level system. The states of this system identify points in a Bloch sphere $S^2$ and their evolution yields trajectories in $S^2$. Quantum measurements turn out to be capable of providing information about the reflection amplitudes. We identify configurations displaying reflectionlessness with situations where one or both of the poles of $S^2$ undergo cyclic evolutions. Spectral singularities and coherent perfect absorption correspond to dynamical spin flips which map one of the poles of $S^2$ to the other.

Scaling trapped-ion optical clocks to hundreds of ions remains a major challenge due to inhomogeneous frequency shifts in large ion crystals. Here, we propose a circular surface-trap architecture for a multi-ion optical clock with hundreds of ions. The rotational symmetry of the ion crystal suppresses the main sources of inhomogeneity, quadrupole shift.We investigate the feasibility of a clock based on a ring crystal of 400 $^{40}\mathrm{Ca}^{+}$ ions and propose a suitable experimental setup configuration. We have shown that circular configuration of ions leads to suppression by two orders of magnitude of the quadrupole shift comparing to the linear trap. Combined with stray-field compensation and dynamical decoupling techniques, this architecture provides a pathway toward multi-ion optical clocks with fractional frequency uncertainty at the $10^{-18}$ level and an expected 20-fold improvement in frequency stability compared to single ion clocks.

Can quantum communication determine the size of an anonymous network when no size bound is known? Classical impossibility proofs rely on the covering lemma: an execution can be reproduced on a larger network with the same local structure. Quantum interference invalidates this argument - a measurement record possible on a network can become impossible on a cover. We show that the impossibility of exact counting nevertheless survives. Our main result is a quantum analogue of the covering lemma that preserves local possibilities rather than complete executions. Within a cyclic covering family, every local measurement outcome possible at a node remains possible at a corresponding node in every sufficiently large network. This holds even after information has travelled around the smaller network, beyond the reach of a simple locality argument. We extend the principle to the universal covering tree and finite covers with sufficiently long cycles. The proof combines tensor-network methods with graph coverings: a common representation across network sizes turns normalization into a constraint on which local outcomes can occur. For uniform synchronous protocols with independently prepared nodes, this support-lifting principle rules out zero-error counting and even the production of a guaranteed finite size bound on anonymous rings and the other cyclic covering families considered here. It also rules out exact leader election with termination in finite time with probability one on these families. Measurements, classical control, and local memory growing with time are allowed; no predetermined running-time bound is assumed. Quantum interference thus defeats execution lifting without enabling infallible certification of network size.

Sample-based quantum linear solvers measure quantum Krylov states $e^{-iAT}|b\rangle$ and reconstruct the solution of $Ax=b$ classically on the observed basis states. We develop this approach for positive-definite systems that arise from elliptic partial differential equations. We show that the sampling guarantees hold in any fixed basis with an efficient transform. We then replace the variational reconstruction by a Galerkin solve on the sampled set. For this reconstruction we derive expected and high-probability energy-norm error bounds, a sample complexity in terms of the smallest sampling probability on an essential set, a computable residual certificate, and a support-completeness theorem for the dephased outcome law. We formulate the method for elliptic problems on Riemannian cross sections. In two dimensions the stiffness matrix is metric-independent in conformal coordinates, and curved vessels lead to the same matrix class. Killing symmetries of the cross section give selection rules for the measurement outcomes, and the Lie derivative of the metric gives first-order bounds for geometric perturbations. We separate operator approximations, whose effect we bound, from sampler approximations, which affect coverage but not the target. As examples we treat the implicit step of pulsatile flow in elliptic, lens-shaped, and curved ducts. In simulations of up to 12 qubits, the number of Fourier modes required for 1% error grows from 14 to 60 while the number of unknowns grows from 64 to 4096. A classical residual-greedy selection is nevertheless more efficient than quantum sampling in all cases studied, consistent with the low measured sampling efficiency. On a six-qubit trapped-ion processor, the measured outcome laws follow a replacement-noise model whose strength grows by about 0.6% per two-qubit gate, and a parity selection rule rejects shots at the predicted rate.

The local Hilbert space dimension is a fundamental parameter in quantum many-body systems, but its role in controlling the ergodicity-breaking transition remains less explored. We systematically study this question by introducing toy models of qudits in zero dimensions that exhibit an ergodicity-breaking transition. Using the hybridization condition, we analytically obtain the dependence of the critical value of the interactions on the local Hilbert space dimension. Numerical studies of spectral correlation measures show agreement with the theoretical predictions, hence confirming the existence of the ergodicity-breaking transition in these models. We further introduce a distinct higher spin version of the model and numerically obtain the same critical behavior, demonstrating that the observed critical point is independent of the microscopic structure of interactions considered.

Quantum thermalization of closed systems – i.e. the relaxation of local properties to their thermodynamic value along the Hamiltonian evolution – is generically expected in non-integrable systems at finite energy density. In this work we show that quantum spin systems with at least one spin component coupled through long-range ferromagnetic interactions, and a second spin component coupled to a field, cannot in general thermalize to a paramagnetic state whenever they have a finite, negative energy per spin. The long-range interactions we consider decay as a power law of the distance, $\sim r^{-\alpha}$ with $\alpha < d$, where $d$ is the number of spatial dimensions of the system. The lack of thermalization is generic to an infinite family of Hamiltonians, and it is not associated to any form of integrability. It can instead be understood by the anomalous features of the thermodynamics of the system in the paramagnetic phase, which cannot be reflected in the properties of the eigenstates. This leads to a breakdown of eigenstate typicality and of the eigenstate thermalization hypothesis. Interestingly, the non-thermal states emerging at long times in the dynamics display long-range multipartite entanglement, elevating the breakdown of thermalization to a resource.

We formulate massless tree-level Compton scattering as an effective two-qubit operation on electron and photon helicities. In the Bell basis, the scattering map separates into two parity sectors and admits a simple quantum circuit representation with one nonunitary filter. After Kraus normalization, this filter defines an unsharp incomplete Bell measurement that distinguishes the two Bell-parity sectors but not the relative Bell label within each sector. We further show that the incomplete measurement alone does not surpass the classical teleportation threshold, whereas a subsequent complete Bell-basis readout enables unit conditional teleportation fidelity.

The Gorini-Kossakowski-Lindblad-Sudarshan (GKLS) equation, a cornerstone of open-system theory, describes Markovian semigroup dynamics in which information about the initial preparation cannot revive. Beyond this regime, information revivals are widely interpreted as backflow from a non-Markovian environment. However, a decrease followed by an increase in information retrievable from the system need not imply that information has returned from the environment [PRX Quantum 6, 020316 (2025)]. It may instead arise from noise correlations unrelated to the system's initial preparation, carried by inert degrees of freedom that never interact with the system and therefore cannot mediate backflow. Such revivals are called extrinsic; genuine backflow is reserved for revivals that cannot be explained this way. The framework developed in that work, however, assumes access to arbitrary initial system states, an assumption often incompatible with physical constraints or limited experimental control. This leaves open which properties of the available preparations are necessary for an information revival to certify genuine backflow. Here we extend the framework to allow the description of restricted families of initial states and show that genuine backflow is possible only if the family contains noncommuting states. If all available states commute, as in a classical stochastic process, every revival can be explained by a suitable inert extension and is therefore extrinsic. Noncommutativity of the available system states is thus necessary for information revivals to provide an operational signature of environmental memory.

We show that criticality concentrates extensive sensing information into a local trajectory. In the rapid-readout limit of an all-to-all dissipative feedback network, a single binary trajectory has a Fisher-information rate linear in $N$ for a common weak field, as does the full record. Critical slowing does not remove this gain because the local weight of the slow mode decreases with $N$. An analytic bound guarantees that the local record retains more than half of the full information across the entire channel-parameter range. We also identify how parameter drives and observation dynamics determine local access to collective information beyond permutation symmetry. Furthermore, our sparse-readout theorem establishes a $k^{-1}$ decay of fractional information loss: only $\Theta(\varepsilon^{-1})$ trajectories are needed to retain a fraction $1-\varepsilon$ of the full information rate. Any growing readout set, even of vanishing density, asymptotically recovers the full common-field Fisher-information rate.