We study ansatz-free Lindbladian learning, the problem of reconstructing the generator of an open quantum system without prior knowledge of its Hamiltonian or dissipator structures. This problem exhibits two fundamental precision limits. Hamiltonian components not obscured by dissipation are Heisenberg-limited, while full Lindbladian reconstruction is subject to the quadratically worse standard quantum limit. This creates an apparent algorithmic impasse: achieving the Heisenberg-limited Hamiltonian learning requires suppressing unknown dissipation, which seemingly demands prior reconstruction of the noise and thereby incurs the standard-quantum-limit cost. In this work, we resolve this obstruction by giving an algorithm that learns the Hamiltonian disjoint from dissipator (HDD) at the Heisenberg limit without prior knowledge of either the Hamiltonian or dissipator supports. Our main technical ingredient is a randomized recursive stabilizer-code construction that progressively identifies and suppresses the dominant dissipative terms without reconstructing the full dissipator. Building on this framework, we also introduce an efficient end-to-end algorithm that learns the entire sparse Lindbladian at the standard quantum limit. Finally, we prove that in the ansatz-free setting, the HDD terms constitute the maximal set of Hamiltonian coefficients uniformly learnable at the Heisenberg limit. We show that every term outside HDD is fundamentally standard-quantum limited, extending prior no-go results to the ansatz-free setting. Together, our protocols provide a scalable framework for characterizing open quantum systems, with quantum error correction serving as a key learning primitive.
Comments
Added the HDD lower bound and reorganized the main text. (10 main-text pages, 111 pages in total, 2 figures)
We study ansatz-free Lindbladian learning, the problem of reconstructing the generator of an open quantum system without prior knowledge of its Hamiltonian or dissipator structures. This problem exhibits two fundamental precision limits. Hamiltonian components not obscured by dissipation are Heisenberg-limited, while full Lindbladian reconstruction is subject to the quadratically worse standard quantum limit. This creates an apparent algorithmic impasse: achieving the Heisenberg-limited Hamiltonian learning requires suppressing unknown dissipation, which seemingly demands prior reconstruction of the noise and thereby incurs the standard-quantum-limit cost. In this work, we resolve this obstruction by giving an algorithm that learns the Hamiltonian disjoint from dissipator (HDD) at the Heisenberg limit without prior knowledge of either the Hamiltonian or dissipator supports. Our main technical ingredient is a randomized recursive stabilizer-code construction that progressively identifies and suppresses the dominant dissipative terms without reconstructing the full dissipator. Building on this framework, we also introduce an efficient end-to-end algorithm that learns the entire sparse Lindbladian at the standard quantum limit. Finally, we prove that in the ansatz-free setting, the HDD terms constitute the maximal set of Hamiltonian coefficients uniformly learnable at the Heisenberg limit. We show that every term outside HDD is fundamentally standard-quantum limited, extending prior no-go results to the ansatz-free setting. Together, our protocols provide a scalable framework for characterizing open quantum systems, with quantum error correction serving as a key learning primitive.
Characterizing the dynamics of open quantum systems at the level of microscopic interactions and error mechanisms is essential for calibrating quantum hardware, designing robust simulation protocols, and developing tailored error-correction methods. Under Markovian noise/dissipation, a natural characterization approach is to identify the full Lindbladian generator that gives rise to both coherent (Hamiltonian) and dissipative dynamics. Prior protocols for learning Lindbladians from dynamical data assumed pre-specified interaction structure, which can be restrictive when the relevant noise channels or control imperfections are not known in advance. In this paper, we present a sample-efficient protocol for learning sparse Lindbladians without assuming any a priori structure. Our protocol is ancilla-free, uses only product-state preparations and Pauli-basis measurements, admits provably stable coefficient reconstruction, and achieves near-optimal time resolution, making it compatible with near-term experimental capabilities. Together, this provides a systematic route to scalable characterization of open-system quantum dynamics, especially in settings where the error mechanisms of interest are unknown.
Describing the dynamics of open systems is essential for fault-tolerant quantum computation. Under Markovian assumptions, we can characterize dissipative dynamics via the Lindbladian. Using Bell sampling, we provide an efficient, ansatz-free Lindbladian learning algorithm with polynomial-time classical postprocessing. Motivated by the prevalence of state preparation and measurement (SPAM) noise on near-term devices, we also introduce the first efficient SPAM-robust protocol capable of learning the gauge-independent components of sparse Lindbladians to arbitrary precision in the presence of constant-order SPAM error. In doing so, we provide the first rigorous characterization of the gauge degrees of freedom in noisy Lindbladian learning, precisely identifying which components remain learnable under SPAM noise.
We study the problem of learning an unknown Markovian open-system generator from access to its physical time evolution. This generator, called a Lindbladian, contains Hamiltonian and dissipative coefficients indexed by an exponentially large family of possible Pauli terms. We propose an efficient algorithm that learns arbitrary Lindbladians from time evolution under minimal assumptions. For a Lindbladian of dynamical strength at most Λ, the algorithm estimates every coefficient to error ε using O(Λ2/ε2) experiments and O(Λ/ε2) total evolution time, together with polynomial classical running time. The algorithm consists of two nonadaptive, ancilla-free, and control-free stages: 1. The support-learning stage outputs a candidate support of size poly(Λ/η) that contains every Hamiltonian and dissipative coordinate of magnitude at least η, using O(Λ2/η2) experiments with preparations of product Pauli eigenstates and single-qubit Pauli measurements. 2.The coefficient-learning stage estimates all coefficients in any candidate support of size M to error ε, using O(Λ2logM/ε2) experiments with preparations of random stabilizer states and measurements in random Clifford bases. Composing the two stages identifies and estimates every coefficient of an arbitrary Lindbladian in polynomial time. The experiment-count and total-evolution-time scalings match the lower bounds up to logarithmic factors, so the algorithm is nearly optimal for learning arbitrary Lindbladians.
We design an algorithm for learning the coefficients of an n-qubit constant-local Lindbladian to ε error with O(gd2log(n)/ε2) total evolution time, where g is the single-site energy and d is the (approximate) degree of the interaction graph. Though Lindbladians present new challenges not present in the special case of Hamiltonians, our algorithm achieves the suite of desiderata attained by state-of-the-art Hamiltonian learning algorithms: (1) it uses non-adaptive, ancilla-free randomized Pauli measurement circuits with a time resolution of only Θ(1/g); (2) it works without knowledge of the structure of the unknown Lindbladian; (3) it depends on a smooth form of degree, thereby supporting the learning of quasi-local and power-law Lindbladians. Moreover, we prove a lower bound showing that our algorithm is optimal in each parameter up to logarithmic factors.
Our algorithm is a simple iterative method, where the objective function consists of Fourier coefficients of the Lindbladian restricted to few-site regions. Its analysis identifies the difficulty unique to open systems, which we call "confusing" terms. For settings where the "confusion" is limited, the performance of the algorithm improves. We demonstrate this for the case of structure learning of Hamiltonians from access to real-time evolution, where we obtain a new algorithm that is significantly simpler than previous work. In addition, using the same iterative method, we design the first efficient algorithm for structure learning Hamiltonians from high-temperature Gibbs states.
Learning the unknown interactions that govern a quantum system is crucial for quantum information processing, device benchmarking, and quantum sensing. The problem, known as Hamiltonian learning, is well understood under the assumption that interactions are local, but this assumption may not hold for arbitrary Hamiltonians. Previous methods all require high-order inverse polynomial dependency with precision, unable to surpass the standard quantum limit and reach the gold-standard Heisenberg-limited scaling. Whether Heisenberg-limited Hamiltonian learning is possible without prior assumptions about the interaction structures, a challenge we term , remains an open question. In this work, we present a quantum algorithm to learn arbitrary sparse Hamiltonians without any structure constraints using only black-box queries of the system’s real-time evolution and minimal digital controls to attain Heisenberg-limited scaling in estimation error. Our method is also resilient to state-preparation-and-measurement errors, enhancing its practical feasibility. We numerically demonstrate our ansatz-free protocol for learning physical Hamiltonians and validating analog quantum simulations, benchmarking our performance against the state-of-the-art Heisenberg-limited learning approach. Moreover, we establish a fundamental trade-off between total evolution time and quantum control on learning arbitrary interactions, revealing the intrinsic interplay between controllability and total evolution-time complexity for any learning algorithm. These results pave the way for further exploration into Heisenberg-limited Hamiltonian learning in complex quantum systems under minimal assumptions, potentially enabling new benchmarking and verification protocols.
We study the problem of learning local Lindbladians from black-box access to the physical evolution, where the goal is to estimate all Hamiltonian and dissipative coefficients. For local Lindbladians with low-intersection dissipators and local dynamical strength at most Λ, we design an algorithm that learns the coefficients to target accuracy ε, whose number of channel uses and total evolution time scale as O(Λ2/ε2) and O(Λ/ε2), respectively, with only logarithmic dependence on the number of qubits. The algorithm is built directly from finite-time channel probes. It runs the unknown evolution for short times, estimates the corresponding Pauli transfer matrices from classical shadows, and converts these estimates into Lindbladian coefficients by stable local Fourier inversions.
The algorithm is non-adaptive, uses no ancillas, and requires only random product states as inputs followed by random Pauli measurements. The method does not require knowing the structure of the Lindbladian in advance. We prove matching lower bounds that establish the near-optimality of the algorithm in both resources. By constructing a family of single-qubit dephasing Lindbladians, we show that any algorithm, even an adaptive one with arbitrary ancillas and measurements, requires Ω(Λ2/ε2) channel uses and Ω(Λ/ε2) total evolution time. In particular, the lower bounds imply that the Heisenberg-limited scaling achievable for Hamiltonian learning is information-theoretically impossible once dissipative coefficients must be estimated.