Optimization of continuous-variable objectives on standard gate-based quantum computers via variational algorithms such as QAOA is typically approached by first discretizing each decision variable into a finite binary representation. This increases qubit requirements and restricts solution precision through fixed-resolution encodings. We propose a qubit-efficient variational framework for continuous optimization that instead encodes each decision variable into the squared amplitude or equivalently, the measurement probability of a single qubit state. This removes explicit discretization from the variable representation while remaining entirely within the standard qubit circuit model unlike methods like CV-QAOA employing qumode based hardware to achieve the same. To read out encoded variables, we propose using amplitude estimation rather than naive sampling or tomographic reconstruction, with the goal of improving precision scaling for continuous-value recovery. We outline how amplitude-estimation error propagates to decision-variable error and then to objective-value error under standard regularity assumptions, suggesting a distinct width-versus-precision tradeoff relative to discretized approaches. In particular, the framework replaces the logarithmic increase in qubits needed for finer binary precision with a constant cost of one qubit per decision variable, while shifting accuracy requirements into the estimation procedure. We position this approach relative to traditional discretized variational formulations, and argue that it provides a promising new direction for continuous optimization on standard qubit architectures.
Variational quantum algorithms offer a promising route to combinatorial optimization, but their applicability is limited by the challenge of encoding large-scale problems within restricted qubit resources. In this work, we introduce a qubit-efficient variational framework based on Pauli correlation encoding (PCE) and apply it to electric power demand portfolio optimization. Binary variables are represented through expectation values of Pauli correlation operators, which encode multi-body correlations of the quantum state and provide a continuous relaxation enabling compact representations with few qubits. We further propose a two-stage hybrid formulation, in which a time-averaged problem provides initialization for a time-resolved optimization. Numerical simulations demonstrate near-optimal performance across problem sizes ranging from
m=18 to
10,296, with normalized cost gaps on the order of
10−4 relative to solutions with certified optimality. We show that the performance is governed by the interplay between continuous relaxation and discretization: the effective resolution of the correlator representation determines how reliably improvements in the continuous loss translate into better discrete solutions, with larger systems exhibiting more consistent behavior. Finally, we demonstrate robustness on a trapped-ion quantum processor, where high-quality solutions are obtained despite noise and finite sampling. These results establish PCE as a physically motivated and qubit-efficient framework for large-scale combinatorial optimization.
We study a Continuous-Variable Quantum Approximate Optimization Algorithm (CV-QAOA) for high-dimensional continuous optimization. Our formulation extends an earlier CV-QAOA proposal with a variationally optimized initial state and recovers the convergence guarantees of Quantum Hamiltonian Descent (QHD) in the high-depth limit. We prove rigorous performance guarantees of CV-QAOA on several families of cost functions. First, we show
d-step CV-QAOA minimizes any
d-dimensional strictly convex quadratic function with
2d quantum queries to the cost function. We then analyze a family of nonconvex "Rotated Double Well" (RDW) functions with
2d local minima introduced by arXiv:2311.00811. While prior work showed QHD reaches its global minimum with
O~(d3) queries, we prove that 1-step CV-QAOA solves RDW with just two quantum queries. Although general-purpose classical solvers need superpolynomial time for RDW and structure-awareness can reduce the cost to polynomial time, we show that the 1-step CV-QAOA protocol can be efficiently dequantized, and that a gradient-aligned line search succeeds with
O(d) queries, nearly matching the information-theoretic
Ω(d/logd) query lower bound. To move beyond the dequantizable regime, we introduce a ``Rotated Square Well'' (RSW) problem, whose globally flat landscape suppresses useful local gradient information. For this family, we show that an adiabatic evolution simulated by CV-QAOA can reach the global minimum using
do(1) queries. On the other hand, any classical algorithm that learn the hidden rotation in RSW provably requires
Ω(d2/logd) queries, a bound we nearly match with an explicit
Θ(d2logd)-query classical algorithm.Numerical simulations on deflected corrugated spring and Easom functions illustrate the promising performance of CV-QAOA on more general problems.
A novel class of hybrid quantum-classical algorithms based on the variational approach have recently emerged from separate proposals addressing, for example, quantum chemistry and combinatorial problems. These algorithms provide an approximate solution to the problem at hand by encoding it in the state of a quantum computer. The operations used to prepare the state are not a priori fixed but, quite the opposite, are subjected to a classical optimization procedure that modifies the quantum gates and improves the quality of the approximate solution. While the quantum hardware determines the size of the problem and what states are achievable (limited, respectively, by the number of qubits and by the kind and number of possible quantum gates), it is the classical optimization procedure that determines the way in which the quantum states are explored and whether the best available solution is actually reached. In addition, the quantities required in the optimization, for example the objective function itself, have to be estimated with finite precision in any experimental implementation. While it is desirable to have very precise estimates, this comes at the cost of repeating the state preparation and measurement multiple times. Here we analyze the competing requirements of high precision and low number of repetitions and study how the overall performance of the variational algorithm is affected by the precision level and the choice of the optimization method. Finally, this study introduces quasi-Newton optimization methods in the general context of hybrid variational algorithms and presents quantitative results for the Quantum Approximate Optimization Algorithm.
We introduce a framework for simulating hybrid oscillator-qubit quantum processors on qubit-only systems through position encoding. By encoding continuous-variable position and momentum wave functions into qubit amplitudes, our method efficiently simulates all Gaussian and conditional Gaussian operations -- encompassing the phase-space instruction set (beam splitter, single-qubit rotation, conditional displacement) and extending to squeezing, conditional squeezing, conditional rotation, and conditional beam splitter -- using
O(log2(Γ+log(1/ε))) qubit gates per hybrid gate, where
Γ is the Fock-level bound and
ε is the target precision. This polylogarithmic per-gate complexity represents an exponential improvement over Fock basis encoding approaches, which require exponential quantum or classical resources in the number of qubits per mode. We provide rigorous numerical characterization of quantum Fourier transform errors for Fock-bounded states, enabling precise resource estimation for practical implementations. This work establishes that hybrid oscillator-qubit algorithms can be implemented on qubit processors with polynomial overhead, providing new insights into the computational power trade-offs between discrete-variable and hybrid continuous-discrete-variable quantum computing.
We investigate the use of amplitude amplification on the gate-based model of quantum computing as a means for solving combinatorial optimization problems. This study focuses primarily on quadratic unconstrained binary optimization (QUBO) problems, which are well-suited for qubit superposition states. Specifically, we demonstrate circuit designs which encode QUBOs as ‘cost oracle’ operations UC, which distribute phases across the basis states proportional to a cost function. We then show that when UC is combined with the standard Grover diffusion operator Us, one can achieve high probabilities of measurement for states corresponding to optimal and near optimal solutions while still only requiring O(π42N/M) iterations. In order to achieve these probabilities, a single scalar parameter ps is required, which we show can be found through a variational quantum–classical hybrid approach and can be used for heuristic solutions.
Quantum Amplitude Estimation (QAE) is a fundamental component of many quantum computing algorithms and applications. The practical implementation of QAE workflows on current Noisy Intermediate-Scale Quantum (NISQ) platforms requires accurate probability amplitude representation and reliable measurement recovery. Understanding the behavior of individual computational stages is therefore essential for improving the overall performance of practical quantum estimation procedures. This study presents a component-wise validation of the amplitude encoding process within a controlled quantum simulation environment. The experimental implementation was developed using the Qiskit framework, where classical probability values were encoded into single-qubit quantum states through parameterized rotation operations and subsequently recovered from measurement statistics. A systematic evaluation was performed across the probability interval from 0.1 to 0.9 using 8192 measurement shots for each experiment. The experimental results demonstrate that the recovered measurement probabilities closely match the corresponding target values throughout the investigated domain. Statistical evaluation produced a Mean Absolute Error (MAE) of 0.003521, a Mean Squared Error (MSE) of 1.67 × 10⁻⁵, and a Root Mean Squared Error (RMSE) of 0.004087, indicating high numerical consistency under ideal simulation conditions. The findings suggest that the amplitude encoding stage accurately represents classical probability values and provides a reliable baseline for studying more advanced quantum estimation workflows. The component-wise validation methodology adopted in this work offers a practical framework for investigating individual stages of quantum algorithms and supports future research on hybrid quantum-classical processing strategies for NISQ-era quantum computing.