Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:1912.07636

Scalable Bayesian Hamiltonian learning

T. J. Evans, Robin Harper, Steven T. Flammia

Published Dec 16, 2019• 33 citations
arXivHTMLPDF

Abstract

As the size of quantum devices continues to grow, the development of scalable methods to characterise and diagnose noise is becoming an increasingly important problem. Recent methods have shown how to efficiently estimate Hamiltonians in principle, but they are poorly conditioned and can only characterize the system up to a scalar factor, making them difficult to use in practice. In this work we present a Bayesian methodology, called Bayesian Hamiltonian Learning (BHL), that addresses both of these issues by making use of any or all, of the following: well-characterised experimental control of Hamiltonian couplings, the preparation of multiple states, and the availability of any prior information for the Hamiltonian. Importantly, BHL can be used online as an adaptive measurement protocol, updating estimates and their corresponding uncertainties as experimental data become available. In addition, we show that multiple input states and control fields enable BHL to reconstruct Hamiltonians that are neither generic nor spatially local. We demonstrate the scalability and accuracy of our method with numerical simulations on up to 100 qubits. These practical results are complemented by several theoretical contributions. We prove that a $k$-body Hamiltonian $H$ whose correlation matrix has a spectral gap $Δ$ can be estimated to precision $\varepsilon$ with only $\tilde{O}\bigl(n^{3k}/(\varepsilon Δ)^{3/2}\bigr)$ measurements. We use two subroutines that may be of independent interest: First, an algorithm to approximate a steady state of $H$ starting from an arbitrary input that converges factorially in the number of samples; and second, an algorithm to estimate the expectation values of $m$ Pauli operators with weight $\le k$ to precision $ε$ using only $O(ε^{-2} 3^k \log m)$ measurements, which quadratically improves a recent result by Cotler and Wilczek.

DOI
10.48550/arxiv.1912.07636

Related papers

Robust and Efficient Hamiltonian Learning

Wenjun Yu, Jinzhao Sun, Zeyao Han +1 more

With the fast development of quantum technology, the sizes of both digital and analog quantum systems increase drastically. In order to have better control and understanding of the quantum hardware, an important task is to characterize the interaction, i.e., to learn the Hamiltonian, which determines both static and dynamic properties of the system. Conventional Hamiltonian learning methods either require costly process tomography or adopt impractical assumptions, such as prior information on the Hamiltonian structure and the ground or thermal states of the system. In this work, we present a robust and efficient Hamiltonian learning method that circumvents these limitations based only on mild assumptions. The proposed method can efficiently learn any Hamiltonian that is sparse on the Pauli basis using only short-time dynamics and local operations without any information on the Hamiltonian or preparing any eigenstates or thermal states. The method has a scalable complexity and a vanishing failure probability regarding the qubit number. Meanwhile, it performs robustly given the presence of state preparation and measurement errors and resiliently against a certain amount of circuit and shot noise. We numerically test the scaling and the estimation accuracy of the method for transverse field Ising Hamiltonian with random interaction strengths and molecular Hamiltonians, both with varying sizes and manually added noise. All these results verify the robustness and efficacy of the method, paving the way for a systematic understanding of the dynamics of large quantum systems.

arXiv:2201.00190•Quantum•Jun 29, 2023•40 citations

Optimised Bayesian system identification in quantum devices

Thomas M. Stace, Jiayin Chen, Li Li +6 more

Identifying and calibrating quantitative dynamical models for physical quantum systems is important for a variety of applications. Here we present a closed-loop Bayesian learning algorithm for estimating multiple unknown parameters in a dynamical model, using optimised experimental "probe" controls and measurement. The estimation algorithm is based on a Bayesian particle filter, and is designed to autonomously choose informationally-optimised probe experiments with which to compare to model predictions. We demonstrate the performance of the algorithm in both simulated calibration tasks and in an experimental single-qubit ion-trap system. Experimentally, we find that with 60x fewer samples, we exceed the precision of conventional calibration methods, delivering an approximately 93x improvement in efficiency (as quantified by the reduction of measurements required to achieve a target residual uncertainty and multiplied by the increase in accuracy). In simulated and experimental demonstrations, we see that successively longer pulses are selected as the posterior uncertainty iteratively decreases, leading to an exponential improvement in the accuracy of model parameters with the number of experimental queries.

arXiv:2211.09090•Nov 16, 2022•1 citation

Adaptive Bayesian phase estimation for quantum error correcting codes

Fernando Martínez-García, Davide Vodola, Markus Müller

Abstract Realisation of experiments even on small and medium-scale quantum computers requires an optimisation of several parameters to achieve high-fidelity operations. As the size of the quantum register increases, the characterisation of quantum states becomes more difficult since the number of parameters to be measured grows as well and finding efficient observables in order to estimate the parameters of the model becomes a crucial task. Here we propose a method relying on application of Bayesian inference that can be used to determine systematic, unknown phase shifts of multi-qubit states. This method offers important advantages as compared to Ramsey-type protocols. First, application of Bayesian inference allows the selection of an adaptive basis for the measurements which yields the optimal amount of information about the phase shifts of the state. Secondly, this method can process the outcomes of different observables at the same time. This leads to a substantial decrease in the resources needed for the estimation of phases, speeding up the state characterisation and optimisation in experimental implementations. The proposed Bayesian inference method can be applied in various physical platforms that are currently used as quantum processors.

arXiv:1904.06166•New Journal of Physics•Nov 27, 2019•14 citations

Optimal Ansatz-free Hamiltonian Learning In Situ

Taiqi Zhou, Weiyuan Gong

Characterizing the features of a Hamiltonian that governs a quantum system serves as a fundamental subroutine of quantum device calibration, signal sensing, and error correction. Recent works have proposed protocols achieving the optimal Heisenberg-limited scaling learning ansatz-free Hamiltonians from their real-time evolutions without fully specifying interaction structures. However, these protocols rely on both deep circuits with interleaving probes and control, and extremely short time resolution, making them difficult to implement on near- and intermediate-term in situ quantum experiments. In this work, we propose a computationally efficient, control-free, and ancilla-free algorithm that uses only Pauli product state preparation and measurement, and learns an ansatz-free Hamiltonian $H$ with $||H||\leqΛ$ in total evolution time of $Θ(\fracΛ{ε^2}\log(\fracΛε))$. The evolution time cost of our algorithm is optimal for any control-free protocols as we further prove a lower bound of $Ω(\fracΛ{ε^2}\log(\fracΛε))$. Technically, our method introduces a randomized-sampling framework that combines band-limited kernel-based time sampling with a displacement sieve for Hamiltonian structure learning. The characteristic probe time resolution depends only on $Λ$ instead of $\varepsilon$, which makes our protocol especially appealing in the high-precision regime for sensing and calibration applications. We also show that the algorithm maintains the same asymptotic total evolution time in the presence of state-preparation-and-measurement (SPAM) noise when the Hamiltonian is local after calibration. Our results demonstrate the fundamental cost of experimentally friendly Hamiltonian learning and provide a practical route to rigorous in situ characterization of near-term quantum platforms.

quant-pharXiv:2606.19486•Jun 19, 2026•0 citations

Multidimensional Quantum Estimation and Model Learning Framework Based on Variational Bayesian Inference

Federico Belliardo, Erik M. Gauger, M. H. Abobeih +3 more

The advancement and scaling of quantum technology has made the learning and identification of quantum systems and devices in highly multidimensional parameter spaces a pressing task for a variety of applications. In many cases, the integration of real-time feedback control and adaptive choice of measurement settings places strict demands on the speed of this task. Here we present a joint model selection and parameter estimation algorithm that is fast and operable on a large number of model parameters. The algorithm is based on variational Bayesian inference, which approximates the target posterior distribution by optimizing a tractable family of distributions, making it more scalable than exact inference methods relying on sampling and that generally suffer from high variance and computational cost in high-dimensional spaces. We show how a regularizing prior can be used to select between competing models, each comprising a different number of parameters, identifying the simplest model that explains the experimental data. The regularization can further separate the degrees of freedom, e.g., quantum systems in the environment or processes, which contribute to major features in the observed dynamics, with respect to others featuring small coupling, which only contribute to a background. As an application of the introduced framework, we consider the problem of the identification of multiple individual nuclear spins with a single electron spin quantum sensor, relevant for nanoscale nuclear magnetic resonance and for the implementation of multiqubit quantum networking nodes. With the number of environmental spins unknown a priori, our Bayesian approach is able to correctly identify the model, i.e., the number of spins and their couplings. We benchmark the algorithm on both simulated and experimental data, using standard figures of merit, and demonstrate that we can estimate dozens of parameters within minutes. Our methodology is compatible with the implementation of real-time adaptive choice of experimental settings in multiparameter estimation, which has the potential to greatly decrease the estimation time in a variety of metrology, tomography, and characterization tasks by prioritizing the measurements with the highest content of information.

arXiv:2507.23130•PRX Quantum•Apr 16, 2026•0 citations

Optimized Bayesian system identification in quantum devices

Thomas M. Stace, Jiayin Chen, Li Li +6 more

Identifying and calibrating quantitative dynamical models for physical quantum systems is important for a variety of applications. Here we present a closed-loop Bayesian learning algorithm for estimating multiple unknown parameters in a dynamical model, using optimized experimental ``probe'' controls and measurement. The estimation algorithm is based on a Bayesian particle filter, and is designed to autonomously choose informationally optimized probe experiments with which to compare to model predictions. We demonstrate the performance of the algorithm in both simulated calibration tasks and in an experimental single-qubit ion-trap system. Experimentally, we find that, with $60$ times fewer samples, we exceed the precision of conventional calibration methods, delivering an approximately $93$ times improvement in efficiency (as quantified by the reduction of measurements and resets required to achieve a target residual uncertainty and multiplied by the increase in accuracy). In simulated and experimental demonstrations, we see that successively longer pulses are selected as the posterior uncertainty iteratively decreases, leading to an exponential improvement in the accuracy of model parameters with the number of experimental queries, and a commensurate increase in the per-query time.

Physical Review Applied•Jan 9, 2024•6 citations
As the size of quantum devices continues to grow, the development of scalable methods to characterise and diagnose noise is becoming an increasingly important problem. Recent methods have shown how to efficiently estimate Hamiltonians in principle, but they are poorly conditioned and can only characterize the system up to a scalar factor, making them difficult to use in practice. In this work we present a Bayesian methodology, called Bayesian Hamiltonian Learning (BHL), that addresses both of these issues by making use of any or all, of the following: well-characterised experimental control of Hamiltonian couplings, the preparation of multiple states, and the availability of any prior information for the Hamiltonian. Importantly, BHL can be used online as an adaptive measurement protocol, updating estimates and their corresponding uncertainties as experimental data become available. In addition, we show that multiple input states and control fields enable BHL to reconstruct Hamiltonians that are neither generic nor spatially local. We demonstrate the scalability and accuracy of our method with numerical simulations on up to 100 qubits. These practical results are complemented by several theoretical contributions. We prove that a kkk-body Hamiltonian HHH whose correlation matrix has a spectral gap ΔΔΔ can be estimated to precision ε\varepsilonε with only O~(n3k/(εΔ)3/2)\tilde{O}\bigl(n^{3k}/(\varepsilon Δ)^{3/2}\bigr)O~(n3k/(εΔ)3/2) measurements. We use two subroutines that may be of independent interest: First, an algorithm to approximate a steady state of HHH starting from an arbitrary input that converges factorially in the number of samples; and second, an algorithm to estimate the expectation values of mmm Pauli operators with weight ≤k\le k≤k to precision εεε using only O(ε−23klog⁡m)O(ε^{-2} 3^k \log m)O(ε−23klogm) measurements, which quadratically improves a recent result by Cotler and Wilczek.
With the fast development of quantum technology, the sizes of both digital and analog quantum systems increase drastically. In order to have better control and understanding of the quantum hardware, an important task is to characterize the interaction, i.e., to learn the Hamiltonian, which determines both static and dynamic properties of the system. Conventional Hamiltonian learning methods either require costly process tomography or adopt impractical assumptions, such as prior information on the Hamiltonian structure and the ground or thermal states of the system. In this work, we present a robust and efficient Hamiltonian learning method that circumvents these limitations based only on mild assumptions. The proposed method can efficiently learn any Hamiltonian that is sparse on the Pauli basis using only short-time dynamics and local operations without any information on the Hamiltonian or preparing any eigenstates or thermal states. The method has a scalable complexity and a vanishing failure probability regarding the qubit number. Meanwhile, it performs robustly given the presence of state preparation and measurement errors and resiliently against a certain amount of circuit and shot noise. We numerically test the scaling and the estimation accuracy of the method for transverse field Ising Hamiltonian with random interaction strengths and molecular Hamiltonians, both with varying sizes and manually added noise. All these results verify the robustness and efficacy of the method, paving the way for a systematic understanding of the dynamics of large quantum systems.
Identifying and calibrating quantitative dynamical models for physical quantum systems is important for a variety of applications. Here we present a closed-loop Bayesian learning algorithm for estimating multiple unknown parameters in a dynamical model, using optimised experimental "probe" controls and measurement. The estimation algorithm is based on a Bayesian particle filter, and is designed to autonomously choose informationally-optimised probe experiments with which to compare to model predictions. We demonstrate the performance of the algorithm in both simulated calibration tasks and in an experimental single-qubit ion-trap system. Experimentally, we find that with 60x fewer samples, we exceed the precision of conventional calibration methods, delivering an approximately 93x improvement in efficiency (as quantified by the reduction of measurements required to achieve a target residual uncertainty and multiplied by the increase in accuracy). In simulated and experimental demonstrations, we see that successively longer pulses are selected as the posterior uncertainty iteratively decreases, leading to an exponential improvement in the accuracy of model parameters with the number of experimental queries.
Abstract Realisation of experiments even on small and medium-scale quantum computers requires an optimisation of several parameters to achieve high-fidelity operations. As the size of the quantum register increases, the characterisation of quantum states becomes more difficult since the number of parameters to be measured grows as well and finding efficient observables in order to estimate the parameters of the model becomes a crucial task. Here we propose a method relying on application of Bayesian inference that can be used to determine systematic, unknown phase shifts of multi-qubit states. This method offers important advantages as compared to Ramsey-type protocols. First, application of Bayesian inference allows the selection of an adaptive basis for the measurements which yields the optimal amount of information about the phase shifts of the state. Secondly, this method can process the outcomes of different observables at the same time. This leads to a substantial decrease in the resources needed for the estimation of phases, speeding up the state characterisation and optimisation in experimental implementations. The proposed Bayesian inference method can be applied in various physical platforms that are currently used as quantum processors.
Characterizing the features of a Hamiltonian that governs a quantum system serves as a fundamental subroutine of quantum device calibration, signal sensing, and error correction. Recent works have proposed protocols achieving the optimal Heisenberg-limited scaling learning ansatz-free Hamiltonians from their real-time evolutions without fully specifying interaction structures. However, these protocols rely on both deep circuits with interleaving probes and control, and extremely short time resolution, making them difficult to implement on near- and intermediate-term in situ quantum experiments. In this work, we propose a computationally efficient, control-free, and ancilla-free algorithm that uses only Pauli product state preparation and measurement, and learns an ansatz-free Hamiltonian HHH with ∣∣H∣∣≤Λ||H||\leqΛ∣∣H∣∣≤Λ in total evolution time of Θ(Λε2log⁡(Λε))Θ(\fracΛ{ε^2}\log(\fracΛε))Θ(ε2Λ​log(εΛ​)). The evolution time cost of our algorithm is optimal for any control-free protocols as we further prove a lower bound of Ω(Λε2log⁡(Λε))Ω(\fracΛ{ε^2}\log(\fracΛε))Ω(ε2Λ​log(εΛ​)). Technically, our method introduces a randomized-sampling framework that combines band-limited kernel-based time sampling with a displacement sieve for Hamiltonian structure learning. The characteristic probe time resolution depends only on ΛΛΛ instead of ε\varepsilonε, which makes our protocol especially appealing in the high-precision regime for sensing and calibration applications. We also show that the algorithm maintains the same asymptotic total evolution time in the presence of state-preparation-and-measurement (SPAM) noise when the Hamiltonian is local after calibration. Our results demonstrate the fundamental cost of experimentally friendly Hamiltonian learning and provide a practical route to rigorous in situ characterization of near-term quantum platforms.
The advancement and scaling of quantum technology has made the learning and identification of quantum systems and devices in highly multidimensional parameter spaces a pressing task for a variety of applications. In many cases, the integration of real-time feedback control and adaptive choice of measurement settings places strict demands on the speed of this task. Here we present a joint model selection and parameter estimation algorithm that is fast and operable on a large number of model parameters. The algorithm is based on variational Bayesian inference, which approximates the target posterior distribution by optimizing a tractable family of distributions, making it more scalable than exact inference methods relying on sampling and that generally suffer from high variance and computational cost in high-dimensional spaces. We show how a regularizing prior can be used to select between competing models, each comprising a different number of parameters, identifying the simplest model that explains the experimental data. The regularization can further separate the degrees of freedom, e.g., quantum systems in the environment or processes, which contribute to major features in the observed dynamics, with respect to others featuring small coupling, which only contribute to a background. As an application of the introduced framework, we consider the problem of the identification of multiple individual nuclear spins with a single electron spin quantum sensor, relevant for nanoscale nuclear magnetic resonance and for the implementation of multiqubit quantum networking nodes. With the number of environmental spins unknown a priori, our Bayesian approach is able to correctly identify the model, i.e., the number of spins and their couplings. We benchmark the algorithm on both simulated and experimental data, using standard figures of merit, and demonstrate that we can estimate dozens of parameters within minutes. Our methodology is compatible with the implementation of real-time adaptive choice of experimental settings in multiparameter estimation, which has the potential to greatly decrease the estimation time in a variety of metrology, tomography, and characterization tasks by prioritizing the measurements with the highest content of information.
Identifying and calibrating quantitative dynamical models for physical quantum systems is important for a variety of applications. Here we present a closed-loop Bayesian learning algorithm for estimating multiple unknown parameters in a dynamical model, using optimized experimental ``probe'' controls and measurement. The estimation algorithm is based on a Bayesian particle filter, and is designed to autonomously choose informationally optimized probe experiments with which to compare to model predictions. We demonstrate the performance of the algorithm in both simulated calibration tasks and in an experimental single-qubit ion-trap system. Experimentally, we find that, with 606060 times fewer samples, we exceed the precision of conventional calibration methods, delivering an approximately 939393 times improvement in efficiency (as quantified by the reduction of measurements and resets required to achieve a target residual uncertainty and multiplied by the increase in accuracy). In simulated and experimental demonstrations, we see that successively longer pulses are selected as the posterior uncertainty iteratively decreases, leading to an exponential improvement in the accuracy of model parameters with the number of experimental queries, and a commensurate increase in the per-query time.