We establish exact scalings for the maximal group velocity of translation-invariant quantum walks in periodic electric fields. Our main result shows that the maximal group velocity decays exponentially with the period of the field in the whole parameter range, thus affirming a conjecture of [3] and at the same time augmenting it to an exact equality. We further demonstrate explicit revival relations and characterize the absolutely continuous spectrum in these models. Our results apply directly also to generalized CMV matrices.

Discrete and continuous variables oftentimes require different treatments in many learning tasks. Identifying the Hamiltonian governing the evolution of a quantum system is a fundamental task in quantum learning theory. While previous works mostly focused on quantum spin systems, where quantum states can be seen as superpositions of discrete bit-strings, relatively little is known about Hamiltonian learning for continuous-variable quantum systems. In this work we focus on learning the Hamiltonian of a bosonic quantum system, a common type of continuous-variable quantum system. This learning task involves an infinite-dimensional Hilbert space and unbounded operators, making mathematically rigorous treatments challenging. We introduce an analytic framework to study the effects of strong dissipation in such systems, enabling a rigorous analysis of cat qubit stabilization via engineered dissipation. This framework also supports the development of Heisenberg-limited algorithms for learning general bosonic Hamiltonians with higher-order terms of the creation and annihilation operators. Notably, our scheme requires a total Hamiltonian evolution time that scales only logarithmically with the number of modes and inversely with the precision of the reconstructed coefficients. On a theoretical level, we derive a new quantitative adiabatic approximation estimate for general Lindbladian evolutions with unbounded generators. Finally, we discuss possible experimental implementations.

We consider the quantum dynamics of a particle on a lattice for large times. Assuming translation invariance, and either discrete or continuous time parameter, the distribution of the ballistically scaled position Q(t)/t converges weakly to a distribution that is compactly supported in velocity space, essentially the distribution of group velocity in the initial state. We show that the total probability of velocities strictly outside the support of the asymptotic measure goes to zero exponentially with t, and we provide a simple method to estimate the exponential rate uniformly in the initial state. Near the boundary of the allowed region the rate function goes to zero like the power 3/2 of the distance to the boundary. The method is illustrated in several examples.

Hamiltonian learning protocols are essential tools to benchmark quantum computers and simulators. Yet rigorous methods for time-dependent Hamiltonians and Lindbladians remain scarce despite their wide use. We close this gap by learning the time-dependent evolution of a locally interacting n-qubit system on a graph of effective dimension D using only preparation of product Pauli eigenstates, evolution under the time-dependent generator for given times, and measurements in product Pauli bases. We assume the time-dependent parameters are well approximated by functions in a known space of dimension m admitting stable interpolation, e.g. by polynomials. Our protocol outputs functions approximating these coefficients to accuracy ε on an interval with success probability 1-δ, requiring only O(ε^-2poly(m)log(nδ^-1)) samples and poly(n,m) pre/postprocessing. Importantly, the scaling in m is polynomial, whereas naive extensions of previous methods scale exponentially. The method estimates time derivatives of observable expectations via interpolation, yielding well-conditioned linear systems for the generator's coefficients. The main difficulty in the time-dependent setting is to evaluate these coefficients at finite times while preserving a controlled link between derivatives and dynamical parameters. Our innovation is to combine Lieb-Robinson bounds, process shadows, and semidefinite programs to recover the coefficients efficiently at constant times. Along the way, we extend state-of-the-art Lieb-Robinson bounds on general graphs to time-dependent, dissipative dynamics, a contribution of independent interest. These results provide a scalable tool to verify state-preparation procedures (e.g. adiabatic protocols) and characterize time-dependent noise in quantum devices.

Quantum hardware has the potential to efficiently solve computationally difficult problems in physics and chemistry to reap enormous practical rewards. Analogue quantum simulation accomplishes this by using the dynamics of a controlled many-body system to mimic those of another system; such a method is feasible on near-term devices. We show that previous theoretical approaches to analogue quantum simulation suffer from fundamental barriers which prohibit scalable experimental implementation. By introducing a new mathematical framework and going beyond the usual toolbox of Hamiltonian complexity theory with an additional resource of engineered dissipation, we show that these barriers can be overcome. This provides a powerful new perspective for the rigorous study of analogue quantum simulators.

Tensor networks provide succinct representations of quantum many-body states and are an important computational tool for strongly correlated quantum systems. Their expressive and computational power is characterized by an underlying entanglement structure, on a lattice or more generally a (hyper)graph, with virtual entangled pairs or multipartite entangled states associated to (hyper)edges. Changing this underlying entanglement structure into another can lead to both theoretical and computational benefits. We study a natural resource theory which generalizes the notion of bond dimension to entanglement structures using multipartite entanglement. It is a direct extension of resource theories of tensors studied in the context of multipartite entanglement and algebraic complexity theory, allowing for the application of the sophisticated methods developed in these fields to tensor networks. The resource theory of tensor networks concerns both the local entanglement structure of a quantum many-body state and the (algebraic) complexity of tensor network contractions using this entanglement structure. We show that there are transformations between entanglement structures which go beyond edge-by-edge conversions, highlighting efficiency gains of our resource theory that mirror those obtained in the search for better matrix multiplication algorithms. We also provide obstructions to the existence of such transformations by extending a variety of methods originally developed in algebraic complexity theory for obtaining complexity lower bounds. The resource theory of tensor networks allows to compare different entanglement structures and should lead to more efficient tensor network representations and contraction algorithms.

Randomized measurements are increasingly appreciated as powerful tools to estimate properties of quantum systems, e.g., in the characterization of hybrid classical-quantum computation. On many platforms they constitute natively accessible measurements, serving as the building block of prominent schemes like shadow estimation. In the real world, however, the implementation of the random gates at the core of these schemes is susceptible to various sources of noise and imperfections, strongly limiting the applicability of protocols. To attenuate the impact of this shortcoming, in this work we introduce an error-mitigated method of randomized measurements, giving rise to a robust shadow estimation procedure. On the practical side, we show that error mitigation and shadow estimation can be carried out using the same session of quantum experiments, hence ensuring that we can address and mitigate the noise affecting the randomization measurements. Mathematically, we develop a picture derived from Fourier-transforms to connect randomized benchmarking and shadow estimation. We prove rigorous performance guarantees and show the functioning using comprehensive numerics. More conceptually, we demonstrate that, if properly used, easily accessible data from randomized benchmarking schemes already provide such valuable diagnostic information to inform about the noise dynamics and to assist in quantum learning procedures.

Characterizing the interactions and dynamics of quantum mechanical systems is an essential task in the development of quantum technologies. We propose an efficient protocol based on the estimation of the time derivatives of few qubit observables using polynomial interpolation for characterizing the underlying Hamiltonian dynamics and Markovian noise of a multi-qubit device. For finite range dynamics, our protocol exponentially relaxes the necessary time resolution of the measurements and quadratically reduces the overall sample complexity compared to previous approaches. Furthermore, we show that our protocol can characterize the dynamics of systems with algebraically decaying interactions. The implementation of the protocol requires only the preparation of product states and single-qubit measurements. Furthermore, we develop a shadow tomography method for quantum channels that is of independent interest. This protocol can be used to parallelize to learn the Hamiltonian, rendering it applicable for the characterization of both current and future quantum devices.

With quantum computing devices increasing in scale and complexity, there is a growing need for tools that obtain precise diagnostic information about quantum operations. However, current quantum devices are only capable of short unstructured gate sequences followed by native measurements. We accept this limitation and turn it into a new paradigm for characterizing quantum gate-sets. A single experiment-random sequence estimation-solves a wealth of estimation problems, with all complexity moved to classical post-processing. We derive robust channel variants of shadow estimation with close-to-optimal performance guarantees and use these as a primitive for partial, compressive and full process tomography as well as the learning of Pauli noise. We discuss applications to the quantum gate engineering cycle, and propose novel methods for the optimization of quantum gates and diagnosing cross-talk.

Reliable quantum technology requires knowledge of the dynamics governing the underlying system. This problem of characterizing and benchmarking quantum devices or experiments in continuous time is referred to as the Hamiltonian learning problem. In contrast to multi-qubit systems, learning guarantees for the dynamics of bosonic systems have hitherto remained mostly unexplored. For $m$-mode Hamiltonians given as polynomials in annihilation and creation operators with modes arranged on a lattice, we establish a simple moment criterion in terms of the particle number operator which ensures that learning strategies from the finite-dimensional setting extend to the bosonic setting, requiring only coherent states and heterodyne detection on the experimental side. We then propose an enhanced procedure based on added dissipation that even works if the Hamiltonian time evolution violates this moment criterion: With high success probability it learns all coefficients of the Hamiltonian to accuracy $\varepsilon$ using a total evolution time of $\mathcal{O}(\varepsilon^{-2}\log(m))$. Our protocol involves the experimentally reachable resources of projected coherent state preparation, dissipative regularization akin to recent quantum error correction schemes involving cat qubits stabilized by a nonlinear multi-photon driven dissipation process, and heterodyne measurements. As a crucial step in our analysis, we establish our moment criterion and a new Lieb-Robinson type bound for the evolution generated by an arbitrary bosonic Hamiltonian of bounded degree in the annihilation and creation operators combined with photon-driven dissipation. Our work demonstrates that a broad class of bosonic Hamiltonians can be efficiently learned from simple quantum experiments, and our bosonic Lieb-Robinson bound may independently serve as a versatile tool for studying evolutions on continuous variable systems.

We study quantum dichotomies and the resource theory of asymmetric distinguishability using a generalization of Strassen's theorem on preordered semirings. We find that an asymptotic variant of relative submajorization, defined on unnormalized dichotomies, is characterized by real-valued monotones that are multiplicative under the tensor product and additive under the direct sum. These strong constraints allow us to classify and explicitly describe all such monotones, leading to a rate formula expressed as an optimization involving sandwiched Renyi divergences. As an application we give a new derivation of the strong converse error exponent in quantum hypothesis testing.

Discrete lattice models are a cornerstone of quantum many-body physics. They arise as effective descriptions of condensed matter systems and lattice-regularized quantum field theories. Lieb-Robinson bounds imply that if the degrees of freedom at each lattice site only interact locally with each other, correlations can only propagate with a finite group velocity through the lattice, similarly to a light cone in relativistic systems. Here we show that Lieb-Robinson bounds are equivalent to the locality of the interactions: a system with k-body interactions fulfills Lieb-Robinson bounds in exponential form if and only if the underlying interactions decay exponentially in space. In particular, our result already follows from the behavior of two-point correlation functions for single-site observables and generalizes to different decay behaviours as well as fermionic lattice models. As a side-result, we thus find that Lieb-Robinson bounds for single-site observables imply Lieb-Robinson bounds for bounded observables with arbitrary support.

Randomized benchmarking refers to a collection of protocols that in the past decade have become central methods for characterizing quantum gates. These protocols aim at efficiently estimating the quality of a set of quantum gates in a way that is resistant to state preparation and measurement errors. Over the years many versions have been developed, however a comprehensive theoretical treatment of randomized benchmarking has been missing. In this work, we develop a rigorous framework of randomized benchmarking general enough to encompass virtually all known protocols as well as novel, more flexible extensions. Overcoming previous limitations on error models and gate sets, this framework allows us, for the first time, to formulate realistic conditions under which we can rigorously guarantee that the output of any randomized benchmarking experiment is well described by a linear combination of matrix exponential decays. We complement this with a detailed analysis of the fitting problem associated with randomized benchmarking data. We introduce modern signal processing techniques to randomized benchmarking, prove analytical sample complexity bounds, and numerically evaluate performance and limitations. In order to reduce the resource demands of this fitting problem, we introduce novel, scalable postprocessing techniques to isolate exponential decays, significantly improving the practical feasibility of a large set of randomized benchmarking protocols. These postprocessing techniques overcome shortcomings in efficiency of several previously proposed methods such as character benchmarking and linear-cross entropy benchmarking. Finally, we discuss, in full generality, how and when randomized benchmarking decay rates can be used to infer quality measures like the average fidelity. On the technical side, our work substantially extends the recently developed Fourier-theoretic perspective on randomized benchmarking by making use of the perturbation theory of invariant subspaces, as well as ideas from signal processing.

This paper uncovers and exploits a link between a central object in harmonic analysis, the so-called Schur functions, and the very hot topic of symmetry protected topological phases of quantum matter. This connection is found in the setting of quantum walks, i.e. quantum analogs of classical random walks. We prove that topological indices classifying symmetry protected topological phases of quantum walks are encoded by matrix Schur functions built out of the walk. This main result of the paper reduces the calculation of these topological indices to a linear algebra problem: calculating symmetry indices of finite-dimensional unitaries obtained by evaluating such matrix Schur functions at the symmetry protected points ± 1 . The Schur representation fully covers the complete set of symmetry indices for 1D quantum walks with a group of symmetries realizing any of the symmetry types of the tenfold way. The main advantage of the Schur approach is its validity in the absence of translation invariance, which allows us to go beyond standard Fourier methods, leading to the complete classification of non-translation invariant phases for typical examples.

The presence of chiral modes on the edges of quantum Hall samples is essential to our understanding of the quantum Hall effect. In particular, these edge modes should support ballistic transport and therefore, in a single particle picture, be supported in the absolutely continuous spectrum of the singleparticle Hamiltonian. We show in this note that if a free fermion system on the two-dimensional lattice is gapped in the bulk, and has a nonvanishing Hall conductance, then the same system put on a half-space geometry supports edge modes whose spectrum fills the entire bulk gap and is absolutely continuous.

With quantum computing devices increasing in scale and complexity, there is a growing need for tools that obtain precise diagnostic information about quantum operations. However, current quantum devices are only capable of short unstructured gate sequences followed by native measurements. We accept this limitation and turn it into a new paradigm for characterizing quantum gate-sets. A single experiment - random sequence estimation - solves a wealth of estimation problems, with all complexity moved to classical post-processing. We derive robust channel variants of shadow estimation with close-to-optimal performance guarantees and use these as a primitive for partial, compressive and full process tomography as well as the learning of Pauli noise. We discuss applications to the quantum gate engineering cycle, and propose novel methods for the optimization of quantum gates and diagnosing cross-talk.

Tensor network states form a variational ansatz class widely used, both analytically and numerically, in the study of quantum many-body systems. It is known that if the underlying graph contains a cycle, e.g., as in projected entangled pair states, then the set of tensor network states of given bond dimension is not closed. Its closure is the tensor network variety. Recent work has shown that states on the boundary of this variety can yield more efficient representations for states of physical interest, but it remained unclear how to systematically find and optimize over such representations. We address this issue by defining an ansatz class of states that includes states at the boundary of the tensor network variety of given bond dimension. We show how to optimize over this class in order to find ground states of local Hamiltonians by only slightly modifying standard algorithms and code for tensor networks. We apply this method to different models and observe favorable energies and runtimes when compared with standard tensor network methods.

We classify chiral symmetric periodically driven quantum systems on a one-dimensional lattice. The driving process is local, can be continuous, or discrete in time, and we assume a gap condition for the corresponding Floquet operator. The analysis is in terms of the unitary operator at a half-period, the half-step operator. We give a complete classification of the connected classes of half-step operators in terms of five integer indices. On the basis of these indices, it can be decided whether the half-step operator can be obtained from a continuous Hamiltonian driving, or not. The half-step operator determines two Floquet operators, obtained by starting the driving at zero or at half-period, respectively. These are called timeframes and are chiral symmetric quantum walks. Conversely, we show under which conditions two chiral symmetric walks determine a common half-step operator. Moreover, we clarify the connection between the classification of half-step operators and the corresponding quantum walks. Within this theory, we prove bulk-edge correspondence and show that a second timeframe allows to distinguish between symmetry protected edge states at +1 and -1 which is not possible for a single timeframe.

We consider the spectral and dynamical properties of one-dimensional quantum walks placed into homogenous electric fields according to a discrete version of the minimal coupling principle. We show that for all irrational fields the absolutely continuous spectrum of these systems is empty, and prove Anderson localization for almost all (irrational) fields. This result closes a gap which was left open in the original study of electric quantum walks: a spectral and dynamical characterization of these systems for typical fields. Additionally, we derive an analytic and explicit expression for the Lyapunov exponent of this model. Making use of a connection between quantum walks and CMV matrices our result implies Anderson localization for CMV matrices with a particular choice of skew-shift Verblunsky coefficients as well as for quasi-periodic unitary band matrices.