Papers

Coupled electronic and nuclear motions govern chemical reactions, yet resolving how electronic structure evolves during bond dissociation remains a central challenge. Here we investigate the photodissociation of Br2 using correlated photoelectron photoion coincidence measurements. A 400 nm pulse initiates dissociation, while strong field ionization probes the evolving molecular system. Coincident measurement of three dimensional photoion and photoelectron momenta provides simultaneous access to the internuclear separation and the accompanying electronic evolution. We identify multiple distinct stages of electronic evolution during bond dissociation. The reshaping of the ionizing molecular orbital occurs first, followed by redistribution and localization of the electronic charge density, and finally by the gradual decay of residual electronic coherence between the separating atomic centers. Between the molecular and atomic limits, we observe an intermediate bond-breaking state in which localized atomic character coexists with a partially delocalized electronic response. By combining correlated observables with semiclassical modelling, we resolve the temporal ordering of these coupled electronic and nuclear processes and determine their associated dynamical timescales. These results demonstrate how correlated momentum observables can disentangle different aspects of the molecular-to-atomic transition.

Continuous-variable measurement-device-independent quantum key distribution (CV-MDI QKD) can address vulnerabilities on the detection side of a QKD system. The core of this protocol involves continuous-variable Bell measurements performed by an untrusted third party. However, in high-speed systems, spectrum broadening causes Bell measurements to deviate from the ideal single-mode scenario, resulting in mode mismatches, reduced performance, and compromised security. Here, we introduce temporal modes (TMs) to analyze the performance of CV-MDI QKD under continuous-mode scenarios. The mismatch between Bob's transmitting mode and Bell measurement mode has a more significant effect on system performance compared to that on Alice's side. When the Bell receiver is close to Bob and the mismatch is set to just 5%, the transmission distance drastically decreases from 87.96 km to 18.50 km. In comparison, the same mismatch for Alice reduces the distance to 86.83 km. This greater degradation on Bob's side can be attributed to the asymmetry in the data modification step. Furthermore, the mismatch in TM characteristics leads to a significant reduction in the secret key rate by 83% when the transmission distance is set to 15 km, which severely limits the practical usability of the protocol over specific distances. These results indicate that in scenarios involving continuous-mode interference, such as large-scale MDI network setups, careful consideration of each user's TM characteristics is crucial. Rigorous pre-calibration of these modes is essential to ensure the system's reliability and efficiency.

Discrete time crystals (DTCs) are non-equilibrium phases in periodically driven systems that exhibit spontaneous breaking of discrete time-translation symmetry. The stabilization of most DTC phases is achieved via the disorder-induced many-body localization. In this work, we propose an experimental scheme to realize disorder-free DTCs in a periodically driven Rydberg atom array. Our scheme utilizes a linear potential in the atomic detuning to enhance the DTC order, without being tired to (Stark) many-body localization. We numerically demonstrate that the Stark potential enhances the robustness of the DTC against the flip imperfections and extends its lifetime, which are independent of initial states. Thus, our scheme provides a promising way to explore DTCs in Rydberg atom arrays without disorder averaging and special state preparation.

Combining the long-coherence of spin qubits and the capability to transmit information and entanglement through photons, spin-photon interfaces (SPIs) are a promising platform for networked quantum computation and long-distance quantum communication. SPIs that possess local `memory' qubits in addition to the optically coupled `communication' qubit can improve remote entanglement fidelities through brokered entanglement schemes and entanglement purification. In these schemes, it is critical to protect the memory qubit from decoherence during entanglement operations on the communications qubit. Silicon, a platform with mature microelectronic and nanophotonic fabrication, is host to the T centre, an SPI with emission in the telecommunications O-band that directly integrates with silicon nanophotonics. Cavity-coupled T centres are a platform for brokered entanglement distribution in silicon photonic circuits and over long-distance optical fibre links. The T centre's electron and nuclear spin qubits are an intrinsic register of communication and memory qubits respectively, with anisotropic hyperfine coupling. In this work we determine the T centre's hydrogen hyperfine coupling tensor. We also introduce schemes to protect against dephasing or eliminate relaxation of the T centre's hydrogen memory qubit during optical excitation. These results address a key challenge for practical T centre quantum networks.

A quantum mechanical setting consisting of a frozen qubit composed with a fully thermalized chaotic system of N states is proposed, with potential relevance to quantum control. Observing the states of the composed system selectively retaining the states leads to the observation of novel localization in the subsystem. At a tuning parameter of 1.0, implying no selection, the system exhibits Wigner-Dyson level spacing statistics, indicative of quantum chaos. As the tuning parameter is reduced and selection occurs at a cutoff, the nearest-neighbor level spacing distribution develops heavier tails, a signature of suppressed spectral mixing and the emergence of non-thermal dynamics. In these regimes, the eigendensity develops a pronounced "cat-ears" structure, reflecting the formation of spatially localized bimodal eigenstates. These topological features persist without transitioning to Poisson statistics, indicating a transition from quantum chaos to a non-thermal, novel many-body localized (MBL) regime-referred to as Wigner Cat Phases. The proposed mixed random matrix ensemble offers a practical probe for sustaining this novel quantum localization setting. Results from our rigorous spectral statistics analysis show how "cat-ears" form in spectral densities based on the degree of selection or disorder and indicate that gap ratio statistics must be used with caution in detecting the full integrable limit due to the possibility of heavy-tailed Wigner-Dyson distributions.

For block-diagonal quantum states with prescribed block weights and dimensions, we determine the minimum von Neumann entropy excess at each trace-norm distance from the blockwise-pure states. The one-block extremal is the classical reverse-Fano profile; the new problem is its constrained assembly when the block weights and dimensions cannot be redistributed. The resulting value function $\Phi_{q,n}(D)$ is attained on the full admissible interval and is the lower envelope of finitely many rank-lattice edges. Every minimizing distance allocation has at most one block between consecutive integer ranks, and the active edge can switch nonmonotonically with $D$. We also determine the lower convex envelope of $\Phi_{q,n}$ explicitly: it is the polygonal curve through synchronized rank vertices and admits a one-parameter dual representation. The exact nonconvex envelope differs from this convex surrogate by at most $q_\diamond\log(5/4)$, with a sharper stage-dependent bound, and the constant is optimal. For binomial charge-sector weights the uniform gap is $O(L^{-1/2})$. The framework yields a symmetry-sector purity certificate; in a three-qubit phase-damping example the fixed-weight certificate remains nontrivial where the universal classical entropy-error envelope is already vacuous.

We explore the correlated dynamics underlying the formation of the quantum Peregrine soliton, a prototypical rogue-wave excitation, utilizing interaction quenches from repulsive to attractive couplings in an ultracold bosonic gas confined in a one-dimensional box trap. The latter emulates the so-called semi-classical initial conditions and the associated gradient catastrophe scenario facilitating the emergence of a high-density, doubly localized waveform. The ensuing multi-orbital variant of the Peregrine soliton features notable deviations from its mean-field sibling, including a reduced peak amplitude, wider core, absence of the side density dips, and earlier formation times. Moreover, Peregrine soliton generation yields coherence losses, while experiencing two-body bunching within each of its sides which show anti-bunching between each other. Controllable seeding of the Peregrine soliton is also demonstrated by tuning the atom number or the box length, while reducing the latter favors the generation of the time-periodic Kuznetsov-Ma breather. Our results highlight that correlations reshape the morphology of rogue-waves in the genuinely quantum, nonintegrable realm, while setting the stage for the emergent field of quantum dispersive hydrodynamics.

Entanglement is one of the key measures of quantum correlations present in nanophotonic systems, with promising applications in quantum optics and beyond. Previous studies have shown that the degree of entanglement between two quantum dot qubits is preserved when a metal nanoparticle is used to mediate the interactions between the qubits. In this work, we investigate long-range plasmonic mediation of qubit–qubit entanglement by studying the impact of the number of mediating metal nanoparticles on stationary concurrence. Collinear and periodically spaced metal nanoparticles that satisfy the weak-coupling approximation are considered. An effective model that enables the derivation of the mediated interactions within the framework of cavity quantum electrodynamics is employed. Under weak driving at the single particle resonance frequency, the model shows that odd-number arrays are more robust to entanglement decay. We attribute this to strong inter-qubit dissipative coupling as a result of a hybridized dipole plasmon resonating with the driving frequency in odd-number arrays. These arrays can sustain non-vanishing stationary entanglement beyond an inter-qubit spacing of one micron, opening up the possibility of independent spatial optical probing of each quantum dot.

We develop a general theoretical framework for computing the time-resolved magneto-optical Kerr effect in ultrafast pump-probe setups, formulated within the Dynamical Projective Operatorial Approach (DPOA) and its application to the generalized linear-response theory for pumped systems. Furthermore, we exploit this formalism to express the post-pump optical conductivity and consequently the Kerr rotation in terms of the time-evolved single-particle density matrix (SPDM), providing a transparent and computationally efficient description of photo-excited multi-band systems. This extension, in addition to its lower computational cost, has the advantage of allowing the inclusion of phenomenological damping. We illustrate the formalism using both (i) a two-band tight-binding model, which captures the essential physics of ultrafast spin-charge dynamics and the Kerr rotation, and (ii) weakly spin-polarized germanium, as a realistic playground with a complex band structure. The results demonstrate that, by exploiting DPOA and/or its SPDM extension, one can reliably reproduce both the short-time features under the pump-pulse envelope and the long-time dynamics after excitation, offering a versatile framework for analyzing time-resolved magneto-optical Kerr effect experiments in complex materials. Moreover, this analysis clearly shows that the Kerr rotation can be used to deduce experimentally the relevant n-photon resonances for a given specific material.

Quantum Key Distribution (QKD) schemes are key exchange protocols based on the physical properties of quantum channels. They avoid the computational-hardness assumptions that underlie the security of classical key exchange. Continuous-Variable QKD (CVQKD), in contrast to qubit-based discrete-variable (DV) schemes, makes use of quadrature measurements of the electromagnetic field. CVQKD has the advantage of being compatible with standard telecom equipment, but at long distances has to deal with very low signal to noise ratios, which necessitates labour-intensive error correction. It is challenging to implement the error correction decoding in realtime. In this paper we introduce a random-codebook error correction method that is suitable for long range Gaussian-modulated CVQKD. We use likelihood ratio scoring with block rejection based on thresholding. For proof-technical reasons, the accept/reject decisions are communicated in encrypted form; in this way we avoid having to deal with non-Gaussian states in the analysis of the leakage. The error correction method is highly parallelisable, which is advantageous for realtime implementation. Under conservative assumptions on the computational resources, we predict a realtime key ratio of at least 8% of the Devetak-Winter value, which outperforms existing reconciliation schemes.

Analog coprocessors for neural networks are an intensively developing field. They provide approximate results of computations for relatively low energy cost and at high speed. We show that a set of ring resonators with Brillouin interaction between photons and phonons, being coupled to a waveguide, can be used to implement matrix-vector multiplication. The input vector is formed by occupancies of the anti-Stokes optical modes pumped via spontaneous Brillouin scattering, i.e, scattering on thermal phonons. Brillouin scattering rates and coupling constants between ring resonators and the waveguide form the matrix. The system allows for parallel computations in frequency band.

We propose and implement an aberration correction method for the creation of extended arrays of spots well beyond the isoplanatic region of any optical system. The method relies on an extensive calibration of aberrations in terms of Zernike polynomials over the full accessible field of an optical system. We introduce a modified Gerchberg-Saxton algorithm for generating holographic phase masks creating fully corrected arbitrary arrays of spots. By applying the method to an aspherical lens, and using a liquid-crystal spatial light modulator (SLM), we increase the aberration-free field of view from 50 to 500 $\mu$m, only limited by the largest diffraction angles accessible to the SLM. This opens new perspectives for the generation of large arrays of optical tweezers, especially for neutral atom based quantum processors and simulators.

We study the effect of a time-varying solenoidal vector potential for a quantum particle confined to a ring. The setup appears to be a time-varying version of the Aharonov-Bohm effect, but since the particle moves in the presence of fields, it is not strictly an Aharonov-Bohm effect. The results are similar to the ac Stark effect, but with a time-varying electric field coming from the vector potential, rather than the scalar potential. We compare and contrast the present effect with the standard ac Stark effect. The signature of this setup is the generation of quasi-energy sidebands which are observable via spectroscopy.

The probabilistic interference fringes observed in the double slit experiment vividly demonstrate the quantum superposition principle, yet they also highlight a fundamental conceptual challenge: the relationship between a system before and after the measurement. According to Copenhagen interpretation, an unobserved quantum system evolves continuously based on the Schrodinger equation, whereas observation induces an instantaneous collapse of the wave function to an eigenstate. This contrast between continuous evolution and sudden collapse renders the single particle behavior particularly enigmatic, especially given that quantum mechanics itself is constructed upon the statistical behavior of ensembles rather than individual entities. In this study, we introduce a Double Slit Diffraction Surrogate Model DSM based on deep learning, designed to capture the mapping between wave functions and probability distributions. The DSM explores multiple potential propagation paths and adaptively selects optimal transmission channels using gradient descent, forming a backbone for the information through the network. By comparing the interpretability of paths and interference, we propose an intuitive physical analogy: the particle behaves like a stretchable dough, extending across both slits, reconnecting after transmission, allowing detachment before the barrier. Monte Carlo simulations confirm that this framework can naturally reproduce the characteristic interference and diffraction probability patterns. Our approach offers a novel, physically interpretable perspective on quantum superposition and measurement induced collapse. The dough analogy is expected to extend to other quantum phenomena. Finally, we provide a dough based picture, attempting to unify interference, entanglement, and tunneling as manifestations of the same underlying phenomenon.

Fluctuation theorems provide universal constraints on nonequilibrium energy and entropy fluctuations, making them a natural framework to assess how and to what extent quantum resources become thermodynamically relevant. We develop a unified framework for incorporating a generic quantum resource, including athermality, quantum coherence, and entanglement, into fluctuation theorems. We work within the end point measurement scheme, which avoids an initial energy measurement and allows quantum resources in the initial state to affect nonequilibrium energy statistics. We derive a family of quantum fluctuation theorems, including generalized Jarzynski equalities and Crooks type fluctuation relations, in which corrections decompose into resource resolved contributions. For single systems, we introduce the concept of weight of athermality, and combine it with the weight of coherence to isolate distinct thermodynamic effects of these quantum resources. For bipartite systems, we furthermore obtain two families of entanglement-resolved fluctuation theorems using an appended correlation operator and the best separable approximation, respectively. Finally, we introduce the concepts of coherence and entanglement fluctuation distances, as Kullback Leibler divergences, which quantify the thermodynamic relevance of quantum resources in a process-dependent and operational manner.

The dynamics of plasmas are governed by a set of non-linear differential equations which remain challenging to solve directly for large 2D and 3D problems. Here we investigate how tensor networks could be applied to plasmas described by the Vlasov-Maxwell system of equations and investigate parameter regimes which show promise for efficient simulations. We show for low-dimensional problems that the simplest form of tensor networks known as a Matrix Product State performs sufficiently well, however in regimes with a strong permanent magnetic field or high-dimensional problems one may need to consider alternative tensor network geometries. We conclude the study of the Vlasov-Maxwell system with the application of tensor networks to an industrially relevant test case and validate our results against state of the art plasma solvers based on Particle-In-Cell codes. We also extend the application of tensor networks to the alternative plasma description of Magnetohydrodynamics and outline how this can be encoded using Matrix Product States.

We study minimum-error identification of an unknown single-bit Boolean function given black-box (oracle) access with one allowed query. Rather than stopping at an abstract optimal measurement, we give a fully constructive solution: an explicit state preparation and an explicit measurement unitary whose computational-basis readout achieves the Helstrom-optimal success probability 3/4 for distinguishing the four possible functions. The resulting circuit is low depth, uses a fixed gate set, and (in this smallest setting) requires no entanglement in the input state. Beyond the specific example, the main message is operational. It highlights a regime in which optimal oracle discrimination is not only well-defined but implementably explicit: the optimal POVM collapses to a compact gate-level primitive that can be compiled, verified, and composed inside larger routines. Motivated by this, we discuss a "what if" question that is open in spirit: for fixed (n,m,k), could optimal k-query identification (possibly for large hypothesis classes) admit deterministic, closed-form descriptions of the inter-query unitaries and the final measurement unitary acting on the natural n+m-qubit input–output registers (and, if needed, small work registers)? Even when such descriptions are not compact and do not evade known circuit-complexity barriers for generic Boolean functions, making the optimum constructive at the circuit level would be valuable for theory-to-hardware translation and for clarifying which forms of "oracle access" are physically meaningful.

Entanglement asymmetry is a relative entropy that faithfully diagnoses symmetry breaking in quantum states, possibly within a spatial subregion. In this work, we extend such framework to higher-form symmetries and compute entanglement asymmetry in theories with spontaneously-broken continuous zero- and higher-form symmetries. One of our central results is an entropic Coleman-Mermin-Wagner theorem, for 0- and $p$-form symmetries, valid also on subregions, which forbids spontaneous breaking of continuous $p$-form symmetries in spacetime dimensions $d\leq p+2$. Our theorem not only qualifies symmetry breaking, it also quantifies it: spontaneous breaking triggers a nonvanishing entanglement asymmetry that grows monotonically towards the infrared, and counts the number of Goldstone fields. Along the way, we derive standalone results concerning the entanglement entropy and asymmetry of Goldstone bosons and gauge fields. In particular, we find a closed-form expression for the Rényi asymmetries of a compact scalar field on spherical subregions in three and four spacetime dimensions, and for higher-form gauge fields in higher dimensions.

Exploring universal entanglement structure in many-body systems is both fundamental and challenging, particularly when the system undergoes non-unitary operations. In this work, we uncover a universal mechanism for realizing maximal entanglement swapping in fermionic Gaussian states subjected to projective Bell measurements. We consider two initially decoupled, half-filled copies of a free-fermion system in arbitrary dimensions and perform post-selective Bell measurements on half of the corresponding sites across the two copies. Remarkably, the post-measurement state factorizes into a product of Bell pairs, establishing maximal interlayer entanglement entirely independent of the initial Gaussian state. We derive this post-measurement state exactly for general particle-number-conserving fermionic Gaussian states, establishing both the validity and universality of the mechanism, with numerical simulations serving as consistency checks. This phenomenon arises from a robust interplay between fermionic statistics and Gaussianity, revealing a distinct fermionic route to measurement-induced maximal entanglement.

Photodynamic therapy (PDT) is a targeted cancer treatment that uses light-activated photosensitizers to generate reactive oxygen species that selectively destroy tumor cells, generally causing less collateral damage than conventional treatments. However, its clinical success hinges on the availability of photosensitizers with strong optical sensitivity and high efficiency in generating reactive oxygen species. While classical computational methods have provided useful insights into photosensitizer design, they struggle to scale and often lack the accuracy needed for these simulations. In this work, we show how fault-tolerant quantum algorithms can be used to identify promising photosensitizer candidates for PDT. To predict photosensitizer performance, we assess two computational properties. First, we quantify light sensitivity by calculating the cumulative absorption in the therapeutic window with a threshold projection algorithm. Second, we determine the efficiency of reactive oxygen generation by estimating intersystem crossing (ISC) rates using the evolution-proxy approach, complemented by a vibronic dynamic treatment where appropriate. We apply these algorithms to a clinically relevant and actively pursued class of photosensitizers, BODIPY derivatives, including heavy-atom and transition-metal-substituted systems that are challenging for classical methods. Our resource estimates, obtained with PennyLane, suggest that systems with active spaces ranging from 11 to 45 spatial orbitals can be simulated using $180$-$350$ logical qubits and Toffoli gate depths between $10^7$ and $10^9$, placing our algorithms within reach of realistic fault-tolerant quantum devices. This paves the way to an efficient quantum-based workflow for designing photosensitizers that can accelerate the discovery of new PDT agents.

In this work, we introduce a technique for reducing the length of a quantum stabilizer code, and we call this deflation of the code. Deflation can be seen as a generalization of the well-known puncturing and shortening techniques in cases where more than a single qudit is removed. We show that the parameters of the deflated quantum code can be controlled, and argue that a similar approach is not as beneficial when applied to classical linear codes. Furthermore, it is shown that deflation introduces additional freedom compared to applying just puncturing and shortening consecutively. We exemplify that it is possible to obtain better parameters by deflating a code rather than consecutively using puncturing and shortening.

Existing protocols for quantum communication networks usually assume an initial allocation of quantum entanglement resources, which are then manipulated through local operations and classical communication (LOCC) to establish high-fidelity entanglement between distant parties. It is generally held that the resulting fidelity would increase monotonically with the entanglement budget. Here, we show that for noncooperative LOCC protocols, the resulting fidelity may decrease as more entanglement is added to a network with non-pure states. This effect results from a quantum analog of selfish routing and constitutes a potential obstacle to the optimal use of resources in large quantum networks.

We investigate the statistical properties of orthogonal, or real, Clifford circuits doped with magic and imaginary resources. By developing the Weingarten calculus for the real Clifford group, we derive the exact overlap distribution of real stabilizer states, identifying a new universality class: the orthogonal Clifford Porter-Thomas distribution. We prove that local real architectures recover this global statistic in logarithmic depth. Furthermore, we uncover a sharp hierarchy in resource requirements: while retrieving Haar statistics necessitates a polylogarithmic amount of magic states, recovering the full unitary Clifford statistics requires only a single phase gate.

We show that a two-proton emitter with a diproton-correlated initial state can act as a source of spin-correlated proton pairs. Using a time-dependent three-body model, we investigate the two-proton emission of $^{16}$Ne ($^{14}$O$+2p$) and analyze the spin correlation of the emitted protons. We find that, when the emission proceeds as a democratic three-body process from an initial state containing a spin-singlet diproton correlation, the emitted protons exhibit a pronounced spin-correlation pattern exceeding the local-hidden-variable bound. This spin correlation closely resembles that of a pure spin-singlet pair. In contrast, this pattern is lost when the process is dominated by the sequential emission or when the initial diproton correlation is absent. These results demonstrate that a certain class of two-proton emitters can deliver spin-entangled proton pairs, and their spin correlation reflects the diproton correlation embedded in the initial state.

We show that a single moving quantum sensor provides complete access to spatially correlated scalar fields. We demonstrate that with either trajectory or internal state control, one can selectively measure any linear functional, e.g. a gradient or a spatial Fourier series coefficient, while successfully eliminating all noise signals with orthogonal spatial correlation. This even exceeds the capabilities of a sensor network consisting of multiple entangled, yet spatially fixed, quantum sensors, where the number of suppressed noise signals is limited by the number of sensor positions. We show that one can achieve an improved scaling of the quantum Fisher information for moving sensors beyond the static fundamental limit of $T^2$.

We develop and employ general Tree Tensor Networks (TTNs) to compute the vibrational spectra for two model systems: a set of 64-dimensional coupled oscillators and acetonitrile. We explore various tree architectures, ranging from the simple linear structure of Matrix Product States (MPS), to trees where only the leaf nodes carry a physical leg – as commonly seen in the underlying ansatz of the Multilayer Multiconfiguration Time-Dependent Hartree (ML-MCTDH) method – and further to more general trees in which all nodes are allowed to possess a physical leg. In addition, we implement Locally Optimal Block Preconditioned Conjugate Gradient (LOBPCG) methods and Inverse Iteration methods as eigensolvers. Benchmarking runtime and accuracy shows that all tested topologies can reach high accuracy. For acetonitrile, inverse-iteration refinement brings all 84 computed states below 1~cm$^{-1}$ error, while the fork-4 tree, a comb-like tree with four backbone nodes, provides the best overall balance between accuracy and cost. MPS remains computationally attractive, whereas more connected trees generally improve accuracy at fixed bond dimension. All numerical simulations were performed using PyTreeNet, a Python package designed for flexible tensor network computations.

In recent years dual-unitary circuits and their multi-unitary generalizations have emerged as exactly solvable yet chaotic models of quantum many-body dynamics. However, a systematic picture for the solvability of multi-unitary dynamics remains missing. We present a framework encompassing a large class of such non-integrable models with exactly solvable dynamics, which we term completely reducible circuits. In these circuits, the entanglement membrane determining operator growth and entanglement dynamics can be characterized analytically. Completely reducible circuits extend the notion of space-time symmetry to more general lattice geometries, breaking dual-unitarity globally but not locally, and allow for a rich phenomenology going beyond dual-unitarity. As example, we introduce circuits that support four and five directions of information flow. We derive a general expression for the entanglement line tension in terms of the pattern of information flow in spacetime. The solvability is shown to be related to the absence of knots of this information flow, connecting entanglement dynamics to the Kauffman bracket as knot invariant. Building on these results, we propose that in general non-integrable dynamics the curvature of the entanglement line tension can be interpreted as a density of information transport. Our results provide a new and unified framework for exactly solvable models of many-body quantum chaos, encompassing and extending known constructions.

Floquet systems display rich phenomena, such as time crystals, with many-body localisation (MBL) protecting the phases from heating. While several types of Floquet phases have been classified, a unified picture of Floquet MBL is still emerging. Static phases have been fruitfully studied via "symmetry topological field theory" (SymTFT), wherein the universal features of $G$-symmetric systems are elucidated by placing them on the boundary of a topological order of one dimension higher. In this work, we provide a SymTFT approach to classifying $G$-symmetric Floquet MBL phases in 1D, for $G$ a finite Abelian group with on-site unitary action. In the SymTFT, these 1D systems correspond to the boundaries of the quantum double associated to $G$, and the classification naturally arises from considering the Lagrangian subgroups and boundary excitations of the quantum double. The classification covers all known Floquet phases while uncovering others previously unexplored, along with bulk features of phases thought to have only boundary signatures. We refer to the latter phases as "dual" time crystals. For static phases, we show how anyons of the quantum double and (string) order parameters provide a natural and simple interpretation of known classification schemes. By extending our framework to the boundaries of twisted quantum doubles, we uncover a new time-crystalline phase with non-onsite symmetry, which cannot be obtained through local, symmetric Hamiltonian drives. We numerically demonstrate evidence for the absolute stability of this phase, and observe that for open boundary conditions it has greater stability to symmetric perturbations. We finally discuss perspectives on using programmable quantum devices to realise and probe the phases we discuss. Our results show that SymTFT provides a powerful approach to unifying phases and features of Floquet systems.

We establish a direct correspondence between the Lanczos approach and the orthogonal polynomials approach in random matrix theory. In the large-$N$ and continuum limits, the average Lanczos coefficients and the recursion coefficients become equivalent, with the precise mapping $b(1-x)=\sqrt{R(x)}$ and $a(1-x)=S(x)$. As a result, the two formalisms yield identical expressions for the leading density of states. We further analyze the Krylov dynamics associated with the recursion coefficients and show that the orthogonal polynomials admit a natural interpretation as Krylov polynomials. This picture is realized explicitly in the Gaussian Unitary Ensemble, where all quantities can be computed analytically.

Efficient simulation of interacting fermionic systems is a key application of near-term quantum computers, but is hindered by the overhead required to encode fermionic operators on qubit hardware. Here, we consider models with $N$ fermionic modes in which each participates in at most a constant number $d$ of interactions and study the circuit depth required to implement the Trotterized time evolution on qubit hardware with all-to-all connectivity. We introduce an encoding that augments each physical fermionic mode with a small number of auxiliary fermions, enabling the removal of Jordan–Wigner strings. Although the preparation of the auxiliary fermion state incurs an initial overhead, this state remains invariant under time evolution. As a result, long-time evolution can be implemented with asymptotically optimal circuit depth, reducing a previously multiplicative $O(\log N)$ overhead to an additive overhead. Our results thus establish that the simulation of sparse fermionic models on qubit hardware matches the performance achievable on ideal fermionic hardware up to constant factors and $O(dN)$ ancillary qubits.