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

Journal of Instrumentation

Optimization-based unfolding in high-energy physics

Simone Gasperini, G. Bianco, Marco Lorusso, Carla Rieger, Michele Grossi

Published Aug 1, 2026• 0 citations
Publisher

Abstract

Abstract In experimental High-Energy Physics, unfolding refers to the problem of estimating the underlying distribution of a physical observable from detector-level data, in the presence of statistical fluctuations and systematic uncertainties. Starting from its reformulation as a regularized quadratic optimization problem, we develop a framework to address unfolding using both classical and quantum-compatible methods. In particular, we derive a Quadratic Unconstrained Binary Optimization (QUBO) representation of the unfolding objective, allowing direct implementation on quantum annealing and hybrid quantum-classical solvers. The proposed approach is implemented in QUnfold, an open-source Python package integrating classical mixed-integer solvers and D-Wave's hybrid quantum solver. We benchmark the method against widely used unfolding techniques in RooUnfold, including response Matrix Inversion, Iterative Bayesian Unfolding, and Singular Value Decomposition unfolding, using synthetic datasets with controlled distortion effects. Our results demonstrate that the optimization-based approach achieves competitive reconstruction accuracy across multiple distributions while naturally accommodating regularization within the objective function. This work establishes a unified optimization perspective on unfolding and provides a practical pathway for exploring quantum-enhanced methods in experimental HEP data analysis.

DOI
10.1088/1748-0221/21/08/t08004

Related papers

Optimization-based Unfolding in High-Energy Physics

Simone Gasperini, Gianluca Bianco, Marco Lorusso +2 more

In experimental High-Energy Physics, unfolding refers to the problem of estimating the underlying distribution of a physical observable from detector-level data, in the presence of statistical fluctuations and systematic uncertainties. Starting from its reformulation as a regularized quadratic optimization problem, we develop a framework to address unfolding using both classical and quantum-compatible methods. In particular, we derive a Quadratic Unconstrained Binary Optimization (QUBO) representation of the unfolding objective, allowing direct implementation on quantum annealing and hybrid quantum-classical solvers. The proposed approach is implemented in QUnfold, an open-source Python package integrating classical mixed-integer solvers and D-Wave's hybrid quantum solver. We benchmark the method against widely used unfolding techniques in RooUnfold, including response Matrix Inversion, Iterative Bayesian Unfolding, and Singular Value Decomposition unfolding, using synthetic datasets with controlled distortion effects. Our results demonstrate that the optimization-based approach achieves competitive reconstruction accuracy across multiple distributions while naturally accommodating regularization within the objective function. This work establishes a unified optimization perspective on unfolding and provides a practical pathway for exploring quantum-enhanced methods in experimental HEP data analysis.

quant-pharXiv:2602.22776•Feb 27, 2026•0 citations

Unfolding measurement distributions via quantum annealing

K. J. R. Cormier, Riccardo Di Sipio, Peter Wittek

A bstract High-energy physics is replete with hard computational problems and it is one of the areas where quantum computing could be used to speed up calculations. We present an implementation of likelihood-based regularized unfolding on a quantum computer. The inverse problem is recast in terms of quadratic unconstrained binary optimization (QUBO), which has the same form of the Ising hamiltonian and hence it is solvable on a programmable quantum annealer. We tested the method using a model that captures the essence of the problem, and compared the results with a baseline method commonly used in precision measurements at the Large Hadron Collider (LHC) at CERN. The unfolded distribution is in very good agreement with the original one. We also show how the method can be extended to include the effect of nuisance parameters representing sources of systematic uncertainties affecting the measurement.

arXiv:1908.08519•Journal of High Energy Physics•Nov 1, 2019•13 citations

A Review of Quantum Supremacy in QUBO-Based Applications

Arim Ryou, Kiwoong Kim, Kyungtaek Jun

We propose a unified approach for reformulating diverse computational challenges into Quadratic Unconstrained Binary Optimization (QUBO) models and solving them via quantum annealing. By applying this framework across multiple domains, we demonstrate its flexibility in handling structurally different problems within a single mathematical formulation. The study highlights QUBO’s capability to exploit quantum parallelism, integrate regularization seamlessly, and address ill-posed inverse problems. These advantages position QUBO-based quantum annealing as a promising tool for large-scale optimization and reconstruction tasks, with strong potential for impact in scientific computing, engineering, and data-driven applications.

IEEE International Conference on Consumer Electronics-Asia (ICCE-Asia)•Oct 27, 2025•0 citations

QVaR: a Quantum Variational Regularization method for Linear Inverse Problems

Siiri Rautio, Hjørdis Schlüter, Andreas Hauptmann +1 more

We present a tailored framework for solving regularized linear inverse problems using quantum optimization methods. By discretizing the solution space and encoding data fidelity and regularization terms into quadratic unconstrained binary optimization (QUBO) models, we formulate both Tikhonov- and sparsity-promoting regularized inverse problems within a unified quantum optimization framework. We further introduce a notion of quantum sensitivity that characterizes the effect of perturbations arising from approximate quantum evolution and discretization. We derive bounds relating these perturbations to stability properties of the underlying inverse problem, thereby establishing a theoretical connection between quantum solution variability and classical notions of ill-posedness. The framework is extended to variational inverse problems through wavelet-based representations and complemented by reduced-order modeling strategies in both parameter and Hamiltonian spaces to mitigate current hardware limitations. Numerical experiments on simulated and physical quantum hardware indicate that the resulting low-energy solution distributions retain information about the underlying inverse problem, while also revealing the limitations imposed by finite-time evolution, discretization, and hardware noise.

arXiv:2606.21392•Jun 19, 2026•0 citations

A Quantum Optimization Framework for Data-Assimilation-Augmented Parameter Estimation

Muhammad Jalil Ahmad, Mohammadhossein Mohammadisiahroudi, Animikh Biswas +1 more

Parameter estimation is a fundamental challenge in the calibration of ordinary differential equation (ODE) models, where repeated numerical integration can lead to high computational cost. In this work, we investigate whether quantum algorithms can be leveraged to assist parameter estimation in nonlinear dynamical systems. We develop a hybrid classical-quantum framework that reformulates a data-assimilation-augmented parameter estimation problem as a combinatorial optimization task. Model dynamics and data assimilation are enforced entirely on the classical side, while the resulting parameter estimation cost functional is discretized and approximated by a quadratic unconstrained binary optimization (QUBO) surrogate. This surrogate is mapped to an Ising Hamiltonian, and quantum optimizers are used to search for low-energy configurations corresponding to candidate parameter estimates. We apply the framework to SIS and SIR epidemic models, the chaotic Lorenz-63 system, and a high-dimensional two-layer Lorenz-96 system. In this setting, the method is used to recover classical system parameters from partial state observations across steady-state, chaotic, and high-dimensional multiscale dynamical systems. Numerical experiments with synthetic data show that the proposed approach accurately recovers parameters while requiring data-assimilation solves only on a prescribed coarse grid. The framework avoids quantum state tomography, illustrating a viable pathway for integrating quantum optimization into data-driven parameter estimation for nonlinear dynamical systems.

quant-pharXiv:2608.13614•Aug 14, 2026•0 citations

Quantum Annealing Methods and Experimental Evaluation to the Phase-Unwrapping Problem in Synthetic Aperture Radar Imaging

Khaled A. Helal Kelany, N.J. Dimopoulos, Clemens P. J. Adolphs +1 more

The focus of this work is to explore the use of quantum annealing solvers for the problem of phase unwrapping of synthetic aperture radar (SAR) images. Although solutions to this problem exist based on network programming, these techniques do not scale well to larger sized images. Our approach involves formulating the problem as a quadratic unconstrained binary optimization (QUBO) problem, which can be solved on a quantum annealer. Given that present embodiments of quantum annealers remain limited in the number of qubits they possess, we decompose the problem into a set of subproblems that can be solved individually. These individual solutions are close to optimal up to an integer constant, with one constant per subimage. In a second phase, these integer constants are determined as a solution to yet another QUBO problem. This basic idea is extended to several passes, where each pass results in an image which is subsequently decomposed to yet another set of subproblems until the resulting image can be accommodated by the annealer at hand. Additionally, we explore improvements to the method by decomposing the original image into overlapping subimages and ignoring the results on the overlapped (marginal) pixels. We test our approach with a variety of software-based QUBO solvers and on a variety of images, both synthetic and real. Additionally, we experiment using D-wave systems’ quantum annealer, the D-wave 2000Q_6 and developed an embedding method which, for our problem, yielded improved results. Our method resulted in high quality solutions, comparable to state-of-the-art phase-unwrapping solvers.

IEEE Transactions on Quantum Engineering•Jan 1, 2022•10 citations
Abstract In experimental High-Energy Physics, unfolding refers to the problem of estimating the underlying distribution of a physical observable from detector-level data, in the presence of statistical fluctuations and systematic uncertainties. Starting from its reformulation as a regularized quadratic optimization problem, we develop a framework to address unfolding using both classical and quantum-compatible methods. In particular, we derive a Quadratic Unconstrained Binary Optimization (QUBO) representation of the unfolding objective, allowing direct implementation on quantum annealing and hybrid quantum-classical solvers. The proposed approach is implemented in QUnfold, an open-source Python package integrating classical mixed-integer solvers and D-Wave's hybrid quantum solver. We benchmark the method against widely used unfolding techniques in RooUnfold, including response Matrix Inversion, Iterative Bayesian Unfolding, and Singular Value Decomposition unfolding, using synthetic datasets with controlled distortion effects. Our results demonstrate that the optimization-based approach achieves competitive reconstruction accuracy across multiple distributions while naturally accommodating regularization within the objective function. This work establishes a unified optimization perspective on unfolding and provides a practical pathway for exploring quantum-enhanced methods in experimental HEP data analysis.
In experimental High-Energy Physics, unfolding refers to the problem of estimating the underlying distribution of a physical observable from detector-level data, in the presence of statistical fluctuations and systematic uncertainties. Starting from its reformulation as a regularized quadratic optimization problem, we develop a framework to address unfolding using both classical and quantum-compatible methods. In particular, we derive a Quadratic Unconstrained Binary Optimization (QUBO) representation of the unfolding objective, allowing direct implementation on quantum annealing and hybrid quantum-classical solvers. The proposed approach is implemented in QUnfold, an open-source Python package integrating classical mixed-integer solvers and D-Wave's hybrid quantum solver. We benchmark the method against widely used unfolding techniques in RooUnfold, including response Matrix Inversion, Iterative Bayesian Unfolding, and Singular Value Decomposition unfolding, using synthetic datasets with controlled distortion effects. Our results demonstrate that the optimization-based approach achieves competitive reconstruction accuracy across multiple distributions while naturally accommodating regularization within the objective function. This work establishes a unified optimization perspective on unfolding and provides a practical pathway for exploring quantum-enhanced methods in experimental HEP data analysis.
A bstract High-energy physics is replete with hard computational problems and it is one of the areas where quantum computing could be used to speed up calculations. We present an implementation of likelihood-based regularized unfolding on a quantum computer. The inverse problem is recast in terms of quadratic unconstrained binary optimization (QUBO), which has the same form of the Ising hamiltonian and hence it is solvable on a programmable quantum annealer. We tested the method using a model that captures the essence of the problem, and compared the results with a baseline method commonly used in precision measurements at the Large Hadron Collider (LHC) at CERN. The unfolded distribution is in very good agreement with the original one. We also show how the method can be extended to include the effect of nuisance parameters representing sources of systematic uncertainties affecting the measurement.
We propose a unified approach for reformulating diverse computational challenges into Quadratic Unconstrained Binary Optimization (QUBO) models and solving them via quantum annealing. By applying this framework across multiple domains, we demonstrate its flexibility in handling structurally different problems within a single mathematical formulation. The study highlights QUBO’s capability to exploit quantum parallelism, integrate regularization seamlessly, and address ill-posed inverse problems. These advantages position QUBO-based quantum annealing as a promising tool for large-scale optimization and reconstruction tasks, with strong potential for impact in scientific computing, engineering, and data-driven applications.
We present a tailored framework for solving regularized linear inverse problems using quantum optimization methods. By discretizing the solution space and encoding data fidelity and regularization terms into quadratic unconstrained binary optimization (QUBO) models, we formulate both Tikhonov- and sparsity-promoting regularized inverse problems within a unified quantum optimization framework. We further introduce a notion of quantum sensitivity that characterizes the effect of perturbations arising from approximate quantum evolution and discretization. We derive bounds relating these perturbations to stability properties of the underlying inverse problem, thereby establishing a theoretical connection between quantum solution variability and classical notions of ill-posedness. The framework is extended to variational inverse problems through wavelet-based representations and complemented by reduced-order modeling strategies in both parameter and Hamiltonian spaces to mitigate current hardware limitations. Numerical experiments on simulated and physical quantum hardware indicate that the resulting low-energy solution distributions retain information about the underlying inverse problem, while also revealing the limitations imposed by finite-time evolution, discretization, and hardware noise.
Parameter estimation is a fundamental challenge in the calibration of ordinary differential equation (ODE) models, where repeated numerical integration can lead to high computational cost. In this work, we investigate whether quantum algorithms can be leveraged to assist parameter estimation in nonlinear dynamical systems. We develop a hybrid classical-quantum framework that reformulates a data-assimilation-augmented parameter estimation problem as a combinatorial optimization task. Model dynamics and data assimilation are enforced entirely on the classical side, while the resulting parameter estimation cost functional is discretized and approximated by a quadratic unconstrained binary optimization (QUBO) surrogate. This surrogate is mapped to an Ising Hamiltonian, and quantum optimizers are used to search for low-energy configurations corresponding to candidate parameter estimates. We apply the framework to SIS and SIR epidemic models, the chaotic Lorenz-63 system, and a high-dimensional two-layer Lorenz-96 system. In this setting, the method is used to recover classical system parameters from partial state observations across steady-state, chaotic, and high-dimensional multiscale dynamical systems. Numerical experiments with synthetic data show that the proposed approach accurately recovers parameters while requiring data-assimilation solves only on a prescribed coarse grid. The framework avoids quantum state tomography, illustrating a viable pathway for integrating quantum optimization into data-driven parameter estimation for nonlinear dynamical systems.
The focus of this work is to explore the use of quantum annealing solvers for the problem of phase unwrapping of synthetic aperture radar (SAR) images. Although solutions to this problem exist based on network programming, these techniques do not scale well to larger sized images. Our approach involves formulating the problem as a quadratic unconstrained binary optimization (QUBO) problem, which can be solved on a quantum annealer. Given that present embodiments of quantum annealers remain limited in the number of qubits they possess, we decompose the problem into a set of subproblems that can be solved individually. These individual solutions are close to optimal up to an integer constant, with one constant per subimage. In a second phase, these integer constants are determined as a solution to yet another QUBO problem. This basic idea is extended to several passes, where each pass results in an image which is subsequently decomposed to yet another set of subproblems until the resulting image can be accommodated by the annealer at hand. Additionally, we explore improvements to the method by decomposing the original image into overlapping subimages and ignoring the results on the overlapped (marginal) pixels. We test our approach with a variety of software-based QUBO solvers and on a variety of images, both synthetic and real. Additionally, we experiment using D-wave systems’ quantum annealer, the D-wave 2000Q_6 and developed an embedding method which, for our problem, yielded improved results. Our method resulted in high quality solutions, comparable to state-of-the-art phase-unwrapping solvers.