QANTIS: A Hardware-Validated Quantum Platform for POMDP Planning and Multi-Target Data Association
Authors: Bayram Yüksel Eker, Suayb S. Arslan, Özgür Nazlı, Mustafa Serhat Demirgil, Furkan Deligöz
Organizations: Neura Parse Ltd., London, United Kingdom · Department of Computer Engineering, Boğaziçi University, Istanbul, Turkey · İzmir Ekonomi Üniversitesi, İzmir, Turkey · Faculty of Engineering, University of Toronto, Toronto, Canada · Istanbul Technical University, Istanbul, Turkey
Autonomous navigation under uncertainty requires solving partially observable Markov decision processes (POMDPs) for planning and assigning sensor measurements to tracked targets--a task known as multi-target data association (MTDA). Both problems become computationally demanding at scale: belief conditioning costs O(P(e)−1) per node under rare evidence, while MTDA is NP-hard. Quantum amplitude amplification can quadratically reduce the belief-update query cost to O(P(e)−1/2), while QUBO reformulations expose MTDA to quantum and quantum-inspired optimisation heuristics. We present QANTIS, a modular platform that integrates quantum belief update (Grover amplitude amplification and BIQAE), QUBO-based data association via FPC-QAOA, and composable error mitigation, and we report a 45-experiment hardware study on three IBM Heron backends. On hardware, a single Grover iterate applied to a Tiger belief oracle amplifies a rare observation probability from 0.179 to 0.907 (5.1×; ISA 18) while preserving the Bayesian posterior (Hellinger 0.0015), increasing usable-shot yield from 1,463 to 7,429. We interpret this as a hardware validation of the quadratic query-complexity mechanism at k=1 with posterior preservation, rather than a wall-clock advantage claim. We further demonstrate, to our knowledge, the first closed-loop hybrid quantum-classical Tiger POMDP on superconducting hardware (T=8, max Hellinger below 0.015), and empirically characterise NISQ feasibility boundaries: ZNE-based error mitigation is beneficial below ISA ≈100 and harmful above ISA ≳1,000; FPC-QAOA is meaningful at ≤15 QUBO variables (ISA ≲450). These results characterise practical operating regimes on current superconducting hardware rather than wall-clock quantum advantage at today's problem scales.
Figures & tables
Symbol
Description
Symbol
Description
S,A,Ω
State/action/observation spaces
P(e)
Evidence probability P(o∣b,a)
T,O,R,γ
Transition/obs/reward/discount
B
Bayesian belief update operator
b∈Δ∣S∣−1
Belief state
H
Planning horizon
N,M
Track/measurement count
xij∈{0,1}
Assignment variable
xi−,Pi−
Predicted state and covariance
λ
QUBO penalty coefficient
dij2
Mahalanobis distance
cij
Association cost (NLL)
Table 1: Principal notation.
System
Q-Plan
Q-DA
HW
EM
CL
Aaronson (2014) [ 10 ]
✓
–
–
–
–
Stollenwerk et al. (2021) [ 18 ]
–
✓
✓
–
–
Cunha et al. (2025) [ 4 ]
✓
–
–
–
–
Ribeiro (2026) [ 23 ]
–
✓
✓
✓
–
TrackHHL (2026) [ 20 ]
–
✓
–
–
–
Marshall et al. (2026) [ 25 ]
–
✓
✓
–
–
Table 2: Positioning of QANTIS relative to prior and concurrent work. Q-Plan: quantum POMDP planning. Q-DA: quantum data association. HW: real quantum hardware experiment. EM: error mitigation. CL: closed-loop action selection.
Figure 1: Quantum belief update circuit for a single POMDP step (schematic; register sizes ns=⌈log2∣S∣⌉ , na=⌈log2∣A∣⌉ , no=⌈log2∣Ω∣⌉ , nr= reward precision bits). Unitaries U1 , U2 , U3 encode the transition, observation, and reward models respectively. Measuring the observation register and conditioning on outcome o (via amplitude amplification) yields the updated belief in the St+1 register. See Appendix A for the full hardware qubit allocation.
P(e)
pomcp
QBRL root
Speedup
10−2
100
≈10
10×
10−3
1,000
≈32
31.6×
10−4
10,000
≈100
100×
Table 3: Per-belief-node query cost (theoretical, fault-tolerant regime; circuit queries ∝P(e)−1/2 require gate fidelities ≲10−6 for large k⋆ ). Speedup is QBRL over pomcp . Hardware demonstration at P(e)≈0.18 : 5.1× amplification at ISA 18 on ibm_marrakesh (Section 8.6 ).
N
M
nvar
ISA ( p=3 )
NISQ feasibility
2
3
11
∼ 435
Validated
3
4
19
∼ 640
Marginal ( 20.4% at p=1 )
5
8
53
> 2 000
Beyond current NISQ
10
15
175
≫ 10 000
FTQC target
Table 4: Estimated qubit resources for MTDA QUBO instances. Gate-model ISA depths use p=3 FPC-QAOA on Heron R2 heavy-hex routing. Our hardware experiments (Section 8.6 ) validate the N=2,M=3 ( n=11 ) instance.
Figure 2: Modular architecture of the QANTIS platform. The three packages maintain a strict acyclic dependency hierarchy: quantum-pomdp and quantum-mht depend on quantum-common but have no cross-dependency. The backend abstraction layer in quantum-common dispatches to three quantum hardware/simulator backends.
Backend
Qubits
Topology
2Q error
Use Case
Status
IBM Heron R1 (ibm_torino)
133
Heavy-hex
∼2.5×10−3
QAOA, BIQAE
Validated
IBM Heron R2 (ibm_fez)
156
Heavy-hex
∼1.5×10−3
QAOA, BIQAE
Validated
IBM Heron R2 (ibm_marrakesh)
156
Heavy-hex
∼1.5×10−3
Tiger, E2E loop
Validated
D-Wave Advantage2 †
4 400+
Zephyr
N/A
Annealing
Companion study
IonQ Aria-2 †
25
All-to-all
∼3×10−3
Small circuits
Planned
Table 5: Quantum hardware backends supported by QANTIS. Validated entries confirmed with real hardware experiments in this work (45-experiment pilot, Section 8.6 ). † Not part of the hardware campaign of this work.
Planner
Cumulative Reward
Belief KL
pomcp
18.3
±
12.7
—
despot
19.1
±
11.4
—
pbvi
17.6
±
13.2
0.000
QBRL (sim.) †
18.7
±
12.1
0.014
Table 6: Tiger problem: mean cumulative reward ( ± std. dev.) over 100 episodes of 50 steps (Qiskit Aer simulation; seed 42). Horizon H=5 . One-way ANOVA: F(3,396)=0.27 , p=0.85 ; all pairwise Welch t -tests p>0.54 ; Cohen ∣d∣<0.09 for all pairs. Tiger ( ∣S∣=2 ) is intentionally simple and does not discriminate among well-tuned planners—agreement of all methods is the expected result. Hardware-validated E2E QBRL results ( T=4 , ibm_marrakesh) are in Table 15 . †
Grid
pomcp
Exact Bayes (classical) / QBRL (sim.) †
Plan Time (ms)
4×4
0.081±0.024
0.043±0.015
12.4 / 34.7
6×6
0.127±0.038
0.068±0.021
41.3 / 112.5
8×8
0.194±0.052
0.091±0.029
108.7 / 347.2
Table 7: GPS-denied navigation: Hellinger distance between belief and true posterior (lower is better) and planning time. Mean ± std. dev., 50 episodes, 30 steps. QBRL (Aer sim.) calls classical_update() internally; it is numerically equivalent to Exact Bayes (classical) under the same action-selection policy. pomcp uses an approximate particle filter (1 000 particles). The Hellinger improvement reflects exact vs. approximate belief update, not quantum vs. classical.
N
M
Vars
Build (ms)
Hungarian (ms)
GNN (ms)
2
3
11
0.3
0.02
0.01
3
5
23
0.5
0.04
0.02
5
8
53
1.2
0.11
0.05
8
12
116
3.1
0.38
0.12
10
15
175
5.4
0.72
0.19
15
23
383
14.8
2.41
0.48
Table 8: QUBO scaling with target count ( M=⌈1.5N⌉ ). Times in milliseconds (single-threaded Python, Intel i7-12700K).
Solver
Obj.
Gap to opt.
Feas. (%)
Time (ms)
Hungarian (optimal)
34.2
0.0%
100.0
0.72
GNN
41.7
+21.9%
100.0
0.19
QAOA ( p=2 )
52.3
+52.9%
68.0
2140
FPC-QAOA ( k=3 , p=8 )
38.9
+13.7%
89.0
4580
Table 9: Solver comparison ( N=10 , M=15 , 175 variables, Aer simulation). Objective lower is better; gap = % above Hungarian optimum. FPC-QAOA outperforms GNN at this scale; hardware feasibility at 175 variables requires fault-tolerant gate fidelities. Real-hardware results at the ISA-feasible 11-variable scale are in Table 15 .
Scenario
CT
MD
FA
tˉ (ms)
Crossing ( N=5 )
5/5
2
4
8.3
Dense clutter ( N=3 )
3/3
1
12
6.1
Swarm patrol
5/5
3
7
14.7
Table 10: End-to-end tracking pipeline (30 steps; Hungarian solver; seed 42). CT = confirmed tracks at final step; MD/FA = events accumulated over 30 steps.
Figure 3: Log-log scaling of solve and build times with target count N (data from Table 8 ). The Hungarian algorithm exhibits the expected cubic scaling; the QUBO build time scales quadratically. GNN provides the fastest heuristic but yields suboptimal assignments. All measurements are from classical simulation (single-threaded Python, Intel i7-12700K).
∣S∣
State qubits
ISA depth
Circuit fidelity
Regime
2
1
12
>0.98
NISQ (confirmed)
4
2
162
>0.78
NISQ (confirmed)
8
3
∼800†
∼0.30
Marginal NISQ
16
4
∼5,000†
∼10−4
FTQC required
64
6
∼100,000†
≈0
FTQC only
Table 11: Belief-update circuit ISA depth vs. state-space size on Heron R2. Confirmed values from hardware; projected values ( † ) extrapolated from O(∣S∣⋅(log2∣S∣)2) CNOT scaling. Heron R2 two-qubit error ∼1.5×10−3 ; circuit fidelity ≈(1−ϵ)ISA .
N
Variables
Hungarian (ms)
Frame rate
Real-time?
20
650
5.9
170 Hz
✓
50
4,100
∼92
11 Hz
✓(slow UAV)
100
15,200
∼730
1.4 Hz
×
200
60,800
∼5,900
0.2 Hz
×
500
376,000
∼91,000
<0.1 Hz
×
Table 12: Hungarian solver runtime vs. target count (extrapolated from N=20 , M=30 measured at 5.87 ms; single-threaded Python, Intel i7-12700K; M=⌈1.5N⌉ ). Real-time budget: 33 ms at 30 Hz.
True a
R1 (torino) ∣a^−a∣
R2 (fez) ∣a^−a∣
R2 (marrakesh) ∣a^−a∣
0.10
0.346 / 0.008 §
—
0.006 §
0.20
0.144 / 0.002 §
0.140
—
0.30
0.040
0.039
—
0.40
0.006
—
—
0.50
0.008
0.0001
—
0.60
0.008
—
—
Table 13: BIQAE amplitude sweep on ibm_torino / ibm_fez / ibm_marrakesh . Uniform-prior runs: 6 iterations, 300 shots, Kt=3t . All 95% credible intervals cover truth. § Prior-matched: prior_mean=atrue , prior_std=0.05 , 10 iter., 500 shots. “—” denotes not tested on this backend.
Method
Wall-clock
Quality
Notes
Hungarian (optimal)
<0.1 ms
100%
O(n3) exact
GNN (greedy)
<0.1 ms
100%
O(NMlogNM) ; optimal at N=2
FPC-QAOA p=3 (Aer)
∼5 s
∼91%
COBYLA convergence
FPC-QAOA p=3 (Heron R2)
∼180 s
64.1%±3.3%
ISA 433, 3 runs
Table 14: Classical vs. quantum solver on the N=2 , M=3 benchmark (Hungarian optimal =−92.4 ; full QUBO objective including false-alarm terms).
Experiment
Metric
Simulator
Hardware
Bell ZNE, 2 q (4096 shots)
⟨ZZ⟩↑
1.000
0.938 raw / 0.953 ZNE
Tiger ZNE, 11 q (4096 shots)
Hellinger ↓
0.004
0.228 raw / 0.277 ZNE ↑
Tiger belief, 11 q (4096 shots)
Hellinger ↓
0.004
0.235
Tiger minimal, 2 q, obs =0 (8192 shots)
Hellinger ↓
0.005
0.025 ¶¶
Tiger minimal, 2 q, obs =1 (8192 shots)
Hellinger ↓
0.007
0.014 ¶¶
FPC-QAOA p=1 , 11 q (4096 shots)
ISA depth / hw obj
2 / −47.4
154 ¶ / −33.0
Table 15: Hardware pilot results on ibm_torino / ibm_fez / ibm_marrakesh (Heron R1/R2, 2026-02-25–26): forty-five experiments across four algorithm families (ZNE/error mitigation, BIQAE, FPC-QAOA, and quantum POMDP/E2E loop). FPC-QAOA hw obj decoded from top-10 bitstrings; Eλ[f] = avg QUBO expectation at scale λ ; optimal Hungarian =−92.4 .
Step
Prior bt
Obs
Action
atrue
BIQAE HW est.
BIQAE err.
CI?
Tiger HW post.
Hellinger
t=0
[0.500,0.500]
0
listen
0.500
0.4945
0.0055
✓
[0.846,0.154]
0.0043
t=1
[0.846,0.154]
0
listen
0.154
0.1715
0.0171
✓
[0.972,0.028]
0.0066
t=2
[0.972,0.028]
0
open-right
0.028
0.0268
0.0012
✓
[0.500,0.500]
0.0003
t=3
[0.500,0.500]
0
listen
0.500
0.4964
0.0033
✓
[0.867,0.133]
0.0169
Table 16: End-to-end quantum POMDP closed loop with QBRL action selection: per-step hardware results on ibm_marrakesh (Heron R2, 2026-02-26, obs =[0,0,0,0] ). BIQAE ISA depth 11; Tiger minimal ISA depth 12–13 (listen) or 4 (open-right at t=2 ). All four BIQAE credible intervals cover ground truth; max Hellinger 0.0169 (PASS). Three additional independent hardware replications of the same pipeline (different shot realisations) yield Tiger max Hellinger 0.0095 , 0.0141 , 0.0148 (all PASS), confirming stable operation.
Step
Obs
Action
atrue
BIQAE HW
Tiger HW post.
Hell.
ISA
t=0
0
listen
0.500
0.503
[0.847,0.153]
0.0031
12
t=1
0
listen
0.153
0.163
[0.973,0.027]
0.0091
13
t=2
0
open-right
0.027
0.029
[0.518,0.482]
0.0128
4
t=3
1
listen
0.482
0.479
[0.175,0.825]
0.0149
13
t=4
1
listen
0.825
0.824
[0.043,0.957]
0.0124
13
t=5
1
open-left
0.957
0.959
[0.505,0.495]
0.0036
4
Table 17: T=8 end-to-end quantum POMDP closed loop on ibm_marrakesh (Heron R2, obs =[0,0,0,1,1,1,0,0] ). Two action events at t=2 (open-right) and t=5 (open-left) trigger ISA 4 Hadamard-only circuits and near-exact belief resets. Max Hellinger 0.0149 (PASS). Single fixed-seed hardware trial; the observation sequence is chosen to elicit two symmetric action events and validate belief resets.
Table 19: Quantum register allocation for the belief update circuit. ns=⌈log2∣S∣⌉ state qubits, no=⌈log2∣Ω∣⌉ observation qubits, np amplitude precision bits.
k
p
Parameters ( 2k )
Objective
Feasibility (%)
2
4
4
48.1
72.0
2
8
4
44.6
78.0
3
4
6
43.2
80.0
3
8
6
38.9
89.0
3
16
6
37.4
91.0
4
8
8
39.5
87.0
Appendix
Table 20: FPC-QAOA sensitivity to parameter count k and circuit depth p . Objective value and feasibility rate are averaged over 50 independent runs. Schedule type: polynomial.
G
Hellinger ( ↓ )
KL Div. ( ↓ )
Amplified P(e)
0
0.312
0.487
0.050
1
0.184
0.201
0.192
2
0.097
0.064
0.485
3
0.038
0.009
0.912
4
0.071
0.031
0.721
5
0.156
0.143
0.298
Appendix
Table 21: Belief update accuracy vs. Grover iterations G for Grid- 4×4 ( ∣S∣=16 , P(e)=0.05 ). Theoretical optimal G∗=⌊4π1/P(e)⌋=3 . Hellinger distance and KL divergence to exact posterior, averaged over 200 belief updates.
Neura Parse Ltd. London, United Kingdom · Dept. of Computer Engineering Boğaziçi University Istanbul, Turkey · Dept. of Physics Izmir University of Economics Izmir, Turkey +2