Paper deep dive
QARIMA: A Quantum Approach To Classical Time Series Analysis
Nishikanta Mohanty, Bikash K. Behera, Badshah Mukherjee, Pravat Dash
Intelligence
Status: succeeded | Model: google/gemini-3.1-flash-lite-preview | Prompt: intel-v1 | Confidence: 95%
Last extracted: 4/14/2026, 1:42:46 AM
Summary
The paper introduces QARIMA, a quantum-inspired ARIMA methodology that replaces classical time series analysis components with quantum-assisted algorithms. It utilizes compact swap-test-driven quantum autocorrelation (QACF) and quantum partial autocorrelation (QPACF) for lag discovery, and employs fixed-configuration Variational Quantum Circuits (VQCs) for autoregressive (AR) and moving-average (MA) parameter estimation, alongside a weak-lag refinement mechanism.
Entities (6)
Relation Signals (3)
QARIMA → integrates → Variational Quantum Circuits
confidence 95% · integrates quantum-assisted lag discovery with fixed-configuration variational quantum circuits (VQCs)
QARIMA → uses → Compact Swap Test
confidence 95% · Differencing and candidate lags are identified via swap-test-driven quantum autocorrelation
QARIMA → improves → ARIMA
confidence 90% · We introduce a hybrid, quantum–inspired ARIMA framework that replaces key classical steps with quantum counterparts
Cypher Suggestions (0)
No Cypher suggestions yet.
Abstract
Abstract:We present a quantum-inspired ARIMA methodology that integrates quantum-assisted lag discovery with fixed-configuration variational quantum circuits (VQCs) for parameter estimation and weak-lag refinement. Differencing and candidate lags are identified via swap-test-driven quantum autocorrelation (QACF) and quantum partial autocorrelation (QPACF), with a delayed-matrix construction that aligns quantum projections to time-domain regressors, followed by standard information-criterion parsimony. Given the screened orders (p,d,q), we retain a fixed VQC ansatz, optimizer, and training budget, preventing hyperparameter leakage, and deploy the circuit in two estimation roles: VQC-AR for autoregressive coefficients and VQC-MA for moving-average coefficients. Between screening and estimation, a lightweight VQC weak-lag refinement re-weights or prunes screened AR lags without altering (p,d,q). Across environmental and industrial datasets, we perform rolling-origin evaluations against automated classical ARIMA, reporting out-of-sample mean squared error (MSE), mean absolute percentage error (MAPE), and Diebold-Mariano tests on MSE and MAE. Empirically, the seven quantum contributions (1) differencing selection, (2) QACF, (3) QPACF, (4) swap-test primitives with delayed-matrix construction, (5) VQC-AR, (6) VQC weak-lag refinement, and (7) VQC-MA collectively reduce meta-optimization overhead and make explicit where quantum effects enter order discovery, lag refinement, and AR/MA parameter estimation.
Tags
Links
- Source: https://arxiv.org/abs/2604.08277v2
- Canonical: https://arxiv.org/abs/2604.08277v2
Trouble viewing inline? Open PDF directly →
Full Text
109,806 characters extracted from source content.
Expand or collapse full text
[1] 1] for Quantum Software and Information, of Technology Sydney, 15 Broadway, Ultimo, , 2007, , 2] ’s Quantum (OPC) Pvt. Ltd., , , 741246, , 3] à degli Studi di Cagliari, Is Mirrions, , 09123, 4] Middle East Fz-Llc., , 9262, QARIMA: A Quantum Approach To Classical Time Series Analysis nishikanta.mohanty@uts.edu.au K. bikas.riki@gmail.com badshah.mukherjee@outlook.com pravat.dash@outlook.com [ [ [ [ Abstract We present a quantum-inspired ARIMA methodology that integrates quantum-assisted lag discovery with fixed-configuration variational quantum circuits (VQCs) for parameter estimation and weak-lag refinement. Differencing and candidate lags are identified via swap-test–driven quantum autocorrelation (QACF) and quantum partial autocorrelation (QPACF), with a delayed-matrix construction that aligns quantum projections to time-domain regressors, followed by standard information-criterion parsimony. Given the screened orders (p,d,q)(p,d,q), we retain a fixed VQC ansatz, optimizer, and training budget preventing hyperparameter leakage and deploy the circuit in two estimation roles: VQC-AR for autoregressive coefficients and VQC-MA for moving-average coefficients. Between screening and estimation, a lightweight VQC weak-lag refinement re-weights or prunes screened AR lags without altering (p,d,q)(p,d,q). Across environmental and industrial datasets, we perform rolling-origin evaluations against automated classical ARIMA, reporting out-of-sample mean squared error (MSE), mean absolute percentage error (MAPE), and Diebold–Mariano tests on MSE and MAE. Empirically, the seven quantum contributions (1) differencing selection, (2) QACF, (3) QPACF, (4) swap-test primitives with delayed-matrix construction, (5) VQC-AR, (6) VQC weak-lag refinement, and (7) VQC-MA collectively reduce meta-optimization overhead and make explicit where quantum effects enter order discovery, lag refinement, and AR/MA parameter estimation. keywords: Compact Swaptest, Variational Quantum Circuits, Quantum Optimisation, ARIMA, Quantum-ACF, Quantum-PACF, Forecast Error Metrics, Diebold-Mariano test 1 Introduction Time series data forecasting is an essential topic in statistical modelling and machine learning, with applications across economics, energy, healthcare, climate research, and more. The AutoRegressive Integrated Moving Average (ARIMA) framework [1] is one of the most prevalent models for both stationary and non-stationary univariate sequences, since it decomposes signals into Auto-Regressive (AR), Differencing (I), and Moving-Average (MA) components. Although ARIMA models are interpretable and computationally efficient, they often fall short in capturing long-range dependencies, nonlinear dynamics, or stochastic effects that are inherent in real-world situations. Furthermore, ARIMA significantly depends on traditional heuristics, including partial autocorrelation functions (PACF), information criteria (e.g., AIC, BIC), and residual-based diagnostics for order estimation and parameter optimization [2]. These approaches may become fragile or computationally intensive in high-dimensional or noisy environments. Recent advancements in quantum computing and quantum-inspired algorithms have created opportunities to overcome these limitations by reconsidering classical learning objectives via quantum principles [3]. Quantum similarity measures (e.g., swap tests, quantum cosine projections), variational quantum circuits (VQCs), and hybrid optimisation loops have shown potential in enhancing learning for structured data [4, 5]. Nonetheless, a comprehensive integration of these methods into time series modelling, especially within the ARIMA framework, remains little investigated. We introduce a hybrid, quantum–inspired ARIMA framework that replaces key classical steps with quantum counterparts while preserving ARIMA’s interpretability and modularity. The pipeline comprises seven coordinated components: (1) differencing selection, (2) quantum autocorrelation (Q–ACF), and (3) quantum partial autocorrelation (Q–PACF) obtained from swap–test projections, together with (4) a delayed–matrix construction aligning quantum measurements to time–domain regressors. Given the screened orders (p,d,q)(p,d,q), estimation proceeds with a fixed–configuration VQC (fixed ansatz, optimizer, and budget) used in two roles: (5) VQC–AR for autoregressive coefficients b and (7) VQC–MA for moving–average coefficients θ, each optimized under a norm–constrained, quantum–augmented loss. Between screening and estimation, (6) a lightweight VQC weak–lag refinement re–weights or prunes screened AR lags without altering (p,d,q)(p,d,q). Instead of selecting lags via classical PACF/ACF heuristics or sweeping AIC/BIC grids, we estimate QACF and QPACF using compact-swap–test projections and phase-corrected cosine similarities; candidate (p,q)(p,q) are gated by these quantum diagnostics, and only then do we apply information-criterion parsimony. This keeps order discovery consistent with the projection logic used later in VQC-refined AR/MA estimation and weak-lag pruning. Classical ARIMA treats model selection and estimation as largely separate, our system integrates order discovery, weak–lag refinement, and parameter estimation within a unified quantum–classical loop that maintains an interpretable ARIMA form. The resulting methodology is theoretically grounded, scalable across datasets, and implementable on near-term hardware through simulation backends, providing a principled bridge between statistical forecasting and quantum machine-learning techniques. 2 Background ARIMA models are among the most widely used approaches for modeling and forecasting univariate time series. The method was originally formalized by Box and Jenkins [6] and has since become foundational in econometrics and applied statistics. Given a time series ytt=1N\y_t\_t=1^N, an ARIMA (p,d,q)(p,d,q) process is defined by differencing the original series d times to induce stationarity and then modeling the result as a combination of AR and MA components. Formally, the model is expressed as: Bp(L)(1−L)dyt=Θq(L)εt, B_p(L)\,(1-L)^d\,y_t\;=\; _q(L)\, _t, (1) where: • L is the lag operator, Lyt=yt−1Ly_t=y_t-1. • Bp(L)=1−b1L−⋯−bpLpB_p(L)=1-b_1L-·s-b_pL^p denotes the AR polynomial. • Θq(L)=1+θ1L+⋯+θqLq _q(L)=1+ _1L+·s+ _qL^q denotes the MA polynomial. • εt _t is white noise, typically assumed Gaussian, εt∼(0,σ2) _t (0,σ^2). After applying (1−L)d(1-L)^d differencing, the ARIMA model specifies that the current observation depends linearly on its past p values and on q past forecast errors: yt=∑i=1pbiyt−i+εt+∑j=1qθjεt−j. y_t= _i=1^pb_i\,y_t-i\;+\; _t\;+\; _j=1^q _j\, _t-j. (2) ARIMA models are commonly applied for: • Forecasting economic indicators (e.g., GDP, inflation). • Modeling financial time series (e.g., stock prices, returns). • Predicting demand in supply chain applications. Parameters (p,d,q)(p,d,q) are typically identified via criteria such as AIC, BIC, or cross-validation. In our preceding section of 3, we will detail the process of calculating the parameters (p,d,q)(p,d,q) using a combination of quantum processes and classical optimization. These parameters will be further tuned by a VQC for improved performance in case-to-case scenarios. 3 Proposed Methodology Our methodology proposes a modular approach to construct Quantum-ARIMA that consists of several distinct phases, each complementing the other. The workflow begins by preparing a lagged representation of the series and selecting the differencing order required for approximate stationarity. Autocorrelation and partial–autocorrelation structure are then inferred by quantum projections: compact swap–test primitives provide phase–corrected cosine similarities that we aggregate into quantum analogues of ACF and PACF. These diagnostics are computed from compact swap-test projections with phase-corrected cosine similarities and feed directly into the delayed-matrix construction used for estimation. A delayed–matrix construction then aligns the quantum measurements with time–domain regressors, and standard information criteria are subsequently applied to retain a concise order triple (p⋆,d⋆,q⋆)(p ,d ,q ). Given these screened orders, estimation proceeds with the same fixed VQC configuration at every rolling origin. Autoregressive coefficients ϕ φ are obtained by minimizing a quantum–augmented loss that combines a prediction–error term with cosine–alignment and projection–entropy penalties derived from the swap–test outputs; an ℓ2 _2 constraint stabilizes the solution. Between screening and estimation we introduce a lightweight weak–lag refinement that re–weights or prunes the screened AR lags without altering (p⋆,d⋆,q⋆)(p ,d ,q ), allowing the model to discourage marginal lags while preserving the selected order. Moving–average coefficients θ are then estimated from the delayed innovation matrix using the same loss family; initialization follows a conditional–least–squares (OLS) fit on the innovations with coefficients clipped for numerical stability. Because the circuit, optimizer, and budget are held constant, the effect of the variational layer is comparable across datasets and time. This section provides the architectural overview; the next section details the concrete procedures. The differencing routine, quantum ACF/PACF construction, and delayed–matrix alignment appear as Algorithms 4–8. The AR estimator, weak–lag refinement, and MA estimator are specified in Algorithms 10–17, and the final ARMA composition used for forecasting is summarized in Algorithm 18. Together these components deliver an interpretable ARIMA model in which quantum effects enter precisely at order discovery, weak–lag refinement, and AR/MA estimation, with all other choices fixed for reproducibility. 4 Processes and Algorithms In this section, we will detail the algorithms that are used to calculate Quantum ARIMA components for parameter estimation. Along with these, we will describe some supporting mechanisms that are used to aid the major algorithms. 4.1 Quantum-Inspired Cosine Similarity via Compact Swap Test To quantify the similarity between the input vector tx_t and a parameter vector θ, we employ a compact swap test formulation [7, 8]. This test estimates the cosine similarity using quantum state overlaps, encoded as probability amplitudes. Given two real-valued vectors ,∈ℝnx, θ ^n, we define normalized states: ϕ=[‖Z,−‖Z],Z=‖2+‖2 φ= [ \|x\| Z,- \| θ\| Z ], Z=\|x\|^2+\| θ\|^2 (3) ψ=12[x1‖,θ1‖,…,xn‖,θn‖] ψ= 1 2 [ x_1\|x\|, _1\| θ\|,…, x_n\|x\|, _n\| θ\| ] (4) These vectors are then padded (if necessary) to match quantum register lengths as powers of 2. The cosine similarity is estimated using the outcome probability of measuring 0 from the control qubit of a controlled-swap (Fredkin) circuit: cos(θ,x)≈2p0−1,p0=Pr(measuring 0) (θ,x)≈ 2p_0-1, p_0=Pr(measuring 0) (5) This estimation is executed using repeated quantum circuit simulations over a specified number of shots. The resulting cosine similarity acts as a multifaceted quantum proxy for partial autocorrelation, AR and MA in time series analysis. Algorithm 1 outlines state preparation, and Algorithm 2 details the overall compact swap test computation. Algorithm 1 State Preparation for Compact Swap Test 1:,∈ℝnx, θ ^n 2:Compute norms: ‖\|x\|, ‖\| θ\| 3:if any norm is zero then 4: return default states [1,0],[1,0][1,0],[1,0] 5:end if 6:Compute Z←‖2+‖2Z←\|x\|^2+\| θ\|^2 7:Compute ϕ←[‖Z,−‖Z]φ← [ \|x\| Z,- \| θ\| Z ] 8:Initialize ψ←[]ψ←[] 9:for i=1i=1 to n do 10: Append to ψ: xi‖2,θi‖2 x_i\|x\| 2, _i\| θ\| 2 11:end for 12:return ϕ,ψφ,ψ Algorithm 2 Compact Swap Test Projection 1:Input ,∈ℝnx, θ ^n, number of shots S 2:Compute ϕ,ψ←φ,ψ← prep-swaptest(,x, θ) 3:Normalize ϕ,ψφ,ψ to unit ℓ2 _2 norm 4:Pad ψ to length 2k2^k for minimal k such that 2k≥len(ψ)2^k (ψ) 5:Construct quantum registers: one control qubit, one ancilla for ϕφ, and k ancillas for ψ 6:Initialize ϕφ into ancilla qubit, ψ into multi-qubit register 7:Apply Hadamard to control qubit 8:Apply controlled-swap between ϕφ and ψ 9:Apply Hadamard to control qubit 10:Measure control qubit S times to estimate p0p_0 11:return cos(θ,x)≈2p0−1 (θ,x)≈ 2p_0-1 Algorithm 3 Compact Swap Test for Scaled Dot Product 1:Vectors ,∈ℝnx, θ ^n, shots S 2:Compute norms: ‖\|x\|, ‖\| θ\| 3:if ‖=0\|x\|=0 or ‖=0\| θ\|=0 then 4: return 0 5:end if 6:Compute state encodings: ϕ,ψ←φ,ψ← prep-swaptest(,x, θ) 7:Normalize: ϕ←ϕ/‖ϕ‖φ←φ/\|φ\|, ψ←ψ/‖ψ‖ψ←ψ/\|ψ\| 8:Compute target length: L←2⌈log2(len(ψ))⌉L← 2 _2(len(ψ)) 9:if len(ψ)<L(ψ)<L then 10: Pad ψ with zeros to length L 11:end if 12:k←log2(L)k← _2(L) ⊳ Number of qubits for ψ 13:Allocate registers: one control qubit q1q_1, 1+k1+k ancilla qubits q2q_2, and 1 classical bit c 14:Initialize ϕφ onto q2[0]q_2[0] and ψ onto q2[1],…,q2[k]q_2[1],…,q_2[k] 15:Apply: X and H to control qubit q1q_1 16:Apply: Controlled-swap between q2[0]q_2[0] and q2[1]q_2[1] 17:Apply: H to q1q_1 and measure into c 18:Simulate circuit for S shots using QASM backend 19:Obtain measurement probabilities: p0←Pr(0)p_0← (0), p1←Pr(1)p_1← (1) 20:Compute swap test projection: Pswap←1−2p0+p1P_swap← 1-2p_0+p_1 21:return ‖⋅‖⋅max(Pswap,0)\|x\|·\| θ\|· (P_swap,0) Figure 1: Compact swap test estimating similarity between encoded states |ϕ⟩ φ and |ψ⟩ ψ. The control qubit is prepared in |0⟩|0 , Hadamard → CSWAP → Hadamard, then measured to yield p0p_0, from which cosθswap _swap and the entropy term H(p0)H(p_0) in ℒARL_AR / ℒMAL_MA are computed. 4.2 Differentiability To determine the optimal differencing order d required to induce stationarity in the time series, we employ a quantum-inspired projection framework in place of traditional unit root tests such as the Augmented Dickey-Fuller (ADF) test [9]. The process begins by constructing a delay matrix from the raw series (Algorithm 4) and generating differenced variants up to a pre-defined maximum order dmaxd_ (Algorithm 5). For each differencing level d, a pair of parameters (α,γ)(α,γ) are estimated to minimize the mean squared error (MSE) between the differenced target DeldDel_d and its predicted value y^t y_t, computed via swap-test–based quantum projection. Specifically, we define the projection-based prediction as y^t=⟨ψ(t)|ψ()⟩, y_t= ψ(x_t)|ψ( θ) , (6) where t=[1,xt]and _t=[1,x_t]\ and\ θ = = [α,γ] [α,γ] (7) are encoded into quantum states using amplitude encoding. The objective at each differencing level is to minimize: ℒd(α,γ)=1N∑t(Deld(t)−y^t)2. _d(α,γ)= 1N _t (Del_d(t)- y_t )^2. (8) The optimization is performed using a classical method such as COBYLA. A convergence mechanism is employed: if the learned γ parameter stabilizes over multiple iterations (indicating no further drift), the process halts early. The selected order d∗d^* is the smallest d for which the loss is minimized and γ≈0γ≈ 0, ensuring that the series is stationary while avoiding unnecessary over-differencing. The complete estimation routine is summarized in Algorithm 6. Algorithm 4 Build Delay Matrix 1:Series y, maximum lag p 2:for i=1i=1 to p do 3: Create column lagi←yt−ilag_i← y_t-i 4:end for 5:Append target yty_t return DataFrame of p lag columns and target yty_t Algorithm 5 Generate Differenced Series 1:Series y, max differencing order dmaxd_ 2:Set prev←yprev← y 3:for i=1i=1 to dmaxd_ do 4: Deli←prev−prevt−1Del_i -prev_t-1 5: prev←Deliprev _i 6:end forreturn DataFrame with Del1,…,DeldmaxDel_1,…,Del_d_ Algorithm 6 Quantum-Inspired Estimation of Differencing Order d 1:Time series y, maximum differencing depth dmaxd_ , lag order p, convergence threshold ϵε, patience T, max optimizer iterations ImaxI_ , primary loss metric ℒL, random seed s 2:Construct delay matrix from y using p lags (Algorithm 4) 3:Generate differenced series up to dmaxd_ (Algorithm 5) 4:Initialize d←0d← 0, metrics_log ←∅← 5:Evaluate d=0d=0: 6:Let xt←lag1(t)x_t _1(t), yt←y(t)y_t← y(t) 7:Optimize (α0,γ0)( _0, _0) to minimize: ℒ0(α,γ)=1N∑t(yt−y^t)2L_0(α,γ)= 1N _t (y_t- y_t )^2 8:Compute predictions via swap test projection: y^t=⟨ψ([1,xt])∣ψ([α,γ])⟩ y_t= ψ([1,x_t]) ψ([α,γ]) 9:Store α0 _0, γ0 _0, ℒ0L_0 in metrics_log 10:for d=1d=1 to dmaxd_ do 11: if DeldDel_d or Deld−1Del_d-1 not defined then 12: Break 13: end if 14: Let xt←Deld−1(t)x_t _d-1(t), yt←Deld(t)y_t _d(t) 15: Initialize optimizer with seed s+ds+d 16: Optimize (αd,γd)( _d, _d) to minimize: ℒd(α,γ)=1N∑t(yt−y^t)2L_d(α,γ)= 1N _t (y_t- y_t )^2 17: Compute predictions via swap test projection: y^t=⟨ψ([1,xt])∣ψ([α,γ])⟩ y_t= ψ([1,x_t]) ψ([α,γ]) 18: Store evaluation metrics in metrics_log 19: if last T values of γd _d satisfy max−min<ϵ - <ε then 20: Break ⊳ Early stopping if γ stabilizes 21: end if 22:end for 23:Select d∗d^* with lowest ℒdL_d and γd≈0 _d≈ 0 return d∗,αd∗,γd∗d^*, _d^*, _d^* The estimation process is governed by several hyperparameters that control convergence behavior, optimization stability, and evaluation fidelity. These include the maximum differencing depth dmaxd_ , maximum number of optimization iterations, convergence threshold for γ, a patience parameter that determines early stopping, and a tunable random seed for initialization. Additionally, the primary evaluation metric such as MSE or mean absolute error (MAE) can be selected to align with specific forecasting goals. Table 1: Hyperparameters for Estimating Differencing Order d Symbol Name Description dmaxd_ Max differencing depth Maximum number of differencing levels tested. p Lag order Number of past lags used to construct the delay matrix. ϵε Convergence threshold Threshold for detecting convergence of γ over recent iterations. T Patience Number of iterations required to confirm γ stability. ImaxI_ Max optimizer iterations Maximum steps allowed per optimization of (α,γ)(α,γ). ℒL Loss metric Primary evaluation metric (e.g., MSE, MAE) used for model selection. s Random seed Seed for reproducibility across differencing levels. 4.3 Quantum-Inspired Autocorrelation Function (ACF) The ACF quantifies linear dependence between yty_t and its lagged values yt−ky_t-k without conditioning on intermediate lags, and is a primary diagnostic for identifying MA structure in classical ARIMA [1, 10]. We introduce a quantum-inspired estimator that replaces classical inner products by overlap amplitudes computed via a compact swap test [7]. Given a univariate series ytt=1N\y_t\_t=1^N (optionally centered and/or differenced), define aligned pairs (yt,yt−k)t=k+1N\(y_t,y_t-k)\_t=k+1^N for lag k. Let ψ(⋅)ψ(·) denote amplitude encoding of a real scalar into a single-qubit state (or a short register when batching); the swap-test estimate of cosine similarity between encodings is obtained from the control-qubit outcome probability p0p_0: cosθ^t,k= 2p0(t,k)−1, θ_t,k\;=\; \,2p_0^(t,k)-1\,\,, (9) which serves as a quantum proxy for normalized correlation at lag k. Averaging over aligned pairs yields the quantum-inspired ACF value ρ^kQ=1N−k∑t=k+1Ncosθ^t,k. ρ^\,Q_k\;=\; 1N-k _t=k+1^N θ_t,k. (10) To enhance robustness, we optionally apply a phase-corrected cosine (as in our PACF module) that blends classical and swap-test angles: cos(ϕt,kcorr)=cos(θswap(t,k)+ω(θdot(t,k)−θswap(t,k))), (φ^corr_t,k )\;=\; \! ( _swap^(t,k)\;+\;ω\, ( _dot^(t,k)- _swap^(t,k) ) ), (11) and set ρ^kQ=1N−k∑tcos(ϕt,kcorr) ρ^\,Q_k= 1N-k _t (φ^corr_t,k). Significant lags are then selected by thresholding |ρ^kQ|| ρ^\,Q_k| using either (i) classical asymptotic bounds ±Z/N± Z/ N (with Z=1.96Z\!=\!1.96 for 95% bands), (i) a percentile rule on magnitudes, or (i) a mean++σ cutoff; a fallback band at βτβτ retains near-significant lags to stabilize downstream MA identification. In classical diagnostics, slowly decaying ACF suggests MA order q [1]; here, the quantum-inspired ACF provides the same role but with overlap-based similarity that is naturally compatible with our swap-test losses and VQC refinement. Algorithm 7 Quantum-Inspired ACF Estimation via Compact Swap Test 1:Series y1:Ny_1:N (optionally centered/differenced), max lag K, shots S, threshold mode Tm∈static,percentile,stdT_m∈\static,percentile,std\, phase weight ω, fallback ratio β 2:ACF values ρ^kQk=1K\ ρ^\,Q_k\_k=1^K, selected lags ℒL 3:Initialize list ACF←[]ACF←[\ ] 4:for k=1k=1 to K do 5: if N≤k+1N≤ k+1 then 6: Append 0 to ACF; continue 7: end if 8: Form aligned pairs (yt,yt−k)t=k+1N\(y_t,y_t-k)\_t=k+1^N 9: Initialize accumulator s←0s← 0, count m←0m← 0 10: for t=k+1t=k+1 to N do 11: Encode yt,yt−ky_t,y_t-k; run compact swap test with S shots to get cosθ^t,k θ_t,k 12: Optionally compute classical cosine and apply phase correction to get cos(ϕt,kcorr) (φ^corr_t,k) 13: s←s+(cos(ϕt,kcorr)orcosθ^t,k)s← s+ ( (φ^corr_t,k)\;or\; θ_t,k ); m←m+1m← m+1 14: end for 15: ρ^kQ←s/m ρ^\,Q_k← s/m; Append ρ^kQ ρ^\,Q_k to ACF 16:end for 17:Compute magnitudes ak←|ρ^kQ|a_k← | ρ^\,Q_k | 18:if Tm=T_m= static then 19: τ←Z/Nτ← Z/ N 20:else if Tm=T_m= percentile then 21: τ←Percentile(a,q)τ (a,q) 22:else if Tm=T_m= std then 23: τ←mean(a)+σstd(a)τ (a)+σ\,std(a) 24:else 25: error 26:end if 27:τf←β⋅τ _f←β·τ; ℒsig←k:ak≥τL_sig←\k:\,a_k≥τ\; ℒfb←k:τf≤ak<τL_fb←\k:\, _f≤ a_k<τ\ 28:return ACF=ρ^kQACF=\ ρ^\,Q_k\, ℒ←ℒsig∪ℒfbL _sig _fb Table 2: Hyperparameters for Quantum-Inspired ACF Estimation Symbol Name Description K max_lag Maximum lag evaluated for ACF. S shots Number of swap-test measurements per pair (yt,yt−k)(y_t,y_t-k). TmT_m threshold_mode Thresholding strategy: static (±Z/N± Z/ N), percentile, or std. Z threshold_z Z-score for static bounds (typically 1.961.96 for 95%). q dynamic_percentile Percentile for magnitude cutoff (e.g., 90th). σ std_multiplier Multiplier for std-based bound (mean ++ σ⋅σ·std). ω phase_weight Weight for phase-corrected cosine blend (optional). β fallback_ratio Fallback fraction of the primary threshold (0<β<10<β<1). In the classical Box–Jenkins workflow, the ACF guides MA order selection by inspecting damped sine waves and cutoffs [1]. Our quantum-inspired ACF preserves this role but replaces sample correlations with overlap-based similarities that (i) align natively with our swap-test losses used later in MA training, and (i) integrate a principled uncertainty control through shot noise and optional entropy regularization in downstream objectives. Together with the quantum-inspired PACF, this yields a consistent diagnostic pair for proposing candidate (p,q)(p,q) before VQC-based refinement. 4.4 Quantum-Inspired PACF Estimation. The PACF is a fundamental diagnostic tool in time series analysis used to determine the order p of the autoregressive (AR) model by measuring the correlation between yty_t and yt−ky_t-k, after eliminating the influence of intermediate lags [1]. In this work, we extend PACF estimation through a quantum-inspired approach that leverages projection-based inner products via the compact swap test, a method to quantify the similarity between the input vector tx_t and a parameter vector θ in a quantum-inspired PACF estimation [7, 8] Algorithms 1, 2. For a given lag k, the projection amplitude is computed as ρ^k=1N−k∑t=k+1N⟨ψ(yt)|ψ(yt−k)⟩, ρ_k= 1N-k _t=k+1^N ψ(y_t) |ψ(y_t-k) , (12) where ψ(⋅)ψ(·) denotes a quantum state encoding of the input scalar into qubit amplitude representations [11]. The swap test is repeated over N−kN-k aligned pairs to generate a mean partial correlation value at lag k. After computing PACF values across all lags up to a maximum K, a thresholding mechanism is applied to select significant lags. We implement three strategies: • a static threshold τ=1.96Nτ= 1.96 N based on classical asymptotic confidence intervals [10] • a dynamic threshold based on the q-th percentile of PACF magnitudes • a standard deviation threshold defined as μ+σμ+σ, where μ and σ are the mean and standard deviation of the PACF magnitudes To enhance robustness, a fallback mechanism selects lags with projections between a fallback threshold τf=β⋅τ _f=β·τ and the primary threshold τ, ensuring the retention of moderately informative lags. This framework preserves the interpretability of classical PACF while embedding quantum-aligned projection logic into lag selection, a technique not commonly adopted in the time series literature, thereby introducing a novel avenue for hybrid quantum-classical model design. Algorithm 8 Quantum-Inspired PACF Estimation Using Swap Test 1:Time series y, maximum lag K, swap test function SwapTest, number of shots S, threshold mode TmT_m, fallback flag f 2:Initialize PACF←[] PACF←[] 3:for k=1k=1 to K do 4: if length(y)≤k+1length(y)≤ k+1 then 5: Append 0 to PACF and continue 6: end if 7: Extract yt←yk+1,…,yNy_t←\y_k+1,…,y_N\, yt−k←y1,…,yN−ky_t-k←\y_1,…,y_N-k\ 8: for each t do 9: Compute projection: ρt(k)=⟨ψ(yt)|ψ(yt−k)⟩≈SwapTest(yt,yt−k,S) _t^(k)= ψ(y_t)|ψ(y_t-k) ≈ SwapTest(y_t,y_t-k,S) 10: end for 11: Compute ρ^k=1N−k∑tρt(k) ρ_k= 1N-k _t _t^(k) 12: Append ρ^k ρ_k to PACF 13:end for 14:Compute threshold τ: 15:if TmT_m = “static” then 16: τ←ZNτ← Z N 17:else if TmT_m = “percentile” then 18: τ←τ← top-q percentile of |PACF|| PACF| 19:else if TmT_m = “std” then 20: τ←μ+στ←μ+σ 21:end if 22:Define fallback threshold τf←β⋅τ _f←β·τ 23:Identify significant lags: ℒs=k∣|ρ^k|≥τ,ℒf=k∣τf≤|ρ^k|<τL_s=\k | ρ_k|≥τ\, _f=\k _f≤| ρ_k|<τ\ 24:Combine: ℒfinal←ℒs∪ℒfL_final _s _f if f=Truef=True return PACF, ℒfinalL_final, τ Table 3: Hyperparameters for Quantum-Inspired PACF Estimation Symbol Name Description K max_lag Maximum lag to compute PACF for. S shots Number of measurement shots used in the swap test circuit. TmT_m threshold_mode Strategy to determine the threshold: “static”, “percentile”, or “std”. Z threshold_numerator Z-score for static thresholding (typically 1.96 for 95%). q dynamic_percentile Percentile threshold (e.g., 90th percentile) for projection magnitude. σ std_multiplier Number of standard deviations above mean (used in “std” mode). β fallback_ratio Proportion of τ used as the fallback threshold. f enable_fallback Boolean flag to include fallback lags in final selection. 4.5 AutoRegressive Order p In ARIMA modeling, the autoregressive (AR) order p governs how many past values of a time series yty_t are used to predict its current value. Classical AR models rely on linear regression, using dot-product formulations over lagged vectors. In our quantum-inspired extension, we replace the classical regression-based fitting with a loss function that combines classical and quantum elements: phase-corrected cosine similarity from compact swap tests, entropy-based uncertainty, and misalignment penalties. The process is formulated as follows. For a lag vector tx_t and coefficient vector b at time t, the classical prediction is: y^t(dot)=t⊤ y_t^(dot)=x_t b (13) We compute the classical cosine similarity as: cosθdot=t⊤^‖t‖,where ^=‖ _dot= x_t b\|x_t\|, b= b\|b\| (14) The quantum-inspired variant estimates the cosine angle using a simulated swap test: cosθswap=SWAP(t,) _swap=SWAP(x_t,b) (15) This estimate is phase-corrected based on the discrepancy between θdot=arccos(cosθdot) _dot= ( _dot) and θswap=arccos(cosθswap) _swap= ( _swap), resulting in the corrected prediction: y^t(quantum)=‖t‖⋅‖⋅cos(θswap+λphase(θdot−θswap)) y_t^(quantum)=\|x_t\|·\|b\|· ( _swap+ _phase( _dot- _swap) ) (16) The total AR loss integrates prediction error, cosine similarity misalignment, and binary entropy uncertainty from the swap test probability: ℒAR=∑t(yt−y^t(quantum))2+λcos∑t(cosθdot−cosθswap)2+λent∑tH(cosθswap) _AR= _t (y_t- y_t^(quantum) )^2+ _cos _t( _dot- _swap)^2+ _ent _tH( _swap) where the binary entropy is computed as: H(p)=−plog2p−(1−p)log2(1−p) H(p)=-p _2p-(1-p) _2(1-p) (18) The following algorithms compute the total loss ℒARL_AR for a given p and coefficient vector b: Algorithm 9 Phase-Corrected Cosine Calculation 1:Classical cosine cosθdot _dot, Quantum cosine cosθswap _swap, Phase correction weight ω 2:Corrected cosine value and intermediate angles 3:Clip classical cosine: c~dot←clip(cosθdot,−1.0,1.0) c_dot ( _dot,-1.0,1.0) 4:Clip quantum cosine: c~swap←clip(cosθswap,−1.0,1.0) c_swap ( _swap,-1.0,1.0) 5:Compute angles: θdot←arccos(c~dot) _dot← ( c_dot), θswap←arccos(c~swap) _swap← ( c_swap) 6:Phase difference: Δθ←θdot−θswap θ← _dot- _swap 7:Corrected angle: θcorr←θswap+ω⋅Δθ _corr← _swap+ω· θ 8:Corrected cosine: cosθcorr←cos(θcorr) _corr← ( _corr) 9:return cosθcorr,θdot,θswap,Δθ,θcorr _corr, _dot, _swap, θ, _corr Algorithm 10 Quantum-Inspired AR Loss Estimation 1:Lagged features t\x_t\, targets yt\y_t\, candidate AR coefficients b, swap test function, hyperparameters λcos,λent,λphase _cos, _ent, _phase 2:Total quantum-inspired AR loss ℒARL_AR 3:Normalize b to obtain unit vector unitb_unit 4:Pad unitb_unit and tx_t to power-of-2 length 5:for each time step t do 6: Compute classical dot product: y^t(dot)←t⊤ y_t^(dot) _t b 7: Compute classical cosine: cosθdot←t⊤unit‖t‖ _dot← x_t b_unit\|x_t\| 8: Estimate cosθswap _swap via compact swap test 9: Compute angular correction: Δθ←θdot−θswap θ← _dot- _swap 10: Corrected angle: θcorr←θswap+λphase⋅Δθ _corr← _swap+ _phase· θ 11: Quantum prediction: y^t(quantum)←‖t‖⋅‖⋅cos(θcorr) y_t^(quantum)←\|x_t\|·\|b\|· ( _corr) 12: Entropy: Ht←−p0log2p0−(1−p0)log2(1−p0)H_t←-p_0 _2p_0-(1-p_0) _2(1-p_0) where p0=1−cos2θswapp_0=1- ^2 _swap 13: Accumulate: ℒ+=(yt−y^t(quantum))2+λcos(cosθdot−cosθswap)2+λentHtL +=(y_t- y_t^(quantum))^2+ _cos( _dot- _swap)^2+ _entH_t 14:end for 15:return ℒL Table 4: Hyperparameters for Quantum-Inspired AR Loss Estimation Symbol Name Description S shots Number of swap test measurement shots used for cosine similarity estimation. λcos _cos lambda_cosine_align Penalty weight applied to alignment difference between classical dot and swap-based cosine. λent _ent lambda_entropy Penalty weight for quantum entropy derived from swap test probability. ω phase_weight Scaling factor to correct the phase discrepancy between classical and quantum cosine. σ noise_std Standard deviation of added Gaussian noise to simulate quantum fluctuations. P padded_dim Size of padded vector, adjusted to nearest power-of-two for swap test register. f save_csv Boolean flag to save log metrics and loss values for each AR trial. 4.5.1 AR Order Estimation via VQC Refinement In classical time series analysis, the autoregressive order p is typically selected by scanning information criteria such as AIC or BIC, or by grid search guided by PACF [1]. In our quantum-inspired framework, we preserve the PACF-guided candidate set P, but replace global search methods with a VQC refinement of the AR coefficients for each candidate p. For each p∈p , we proceed as follows: 1. Difference the series using the identified order d⋆d ; 2. Build the AR design matrix (p)X^(p) from lagged values; 3. Initialize the coefficient vector (0)b^(0) via ordinary least squares (OLS); 4. Refine b by minimizing the quantum-inspired AR loss ℒARL_AR (cf. Eq. (4.5)) using a depth-controlled VQC ansatz and a classical derivative-free optimizer such as COBYLA [4, 12]. The VQC offers an expressive, hardware-compatible parameterization while retaining PACF interpretability. OLS initialization provides a warm start that reduces susceptibility to barren plateaus and poor local minima [13]. Let ∈ℝTy ^T denote the differenced series. For each p∈p , we form (p)∈ℝ(T−p)×pX^(p) ^(T-p)× p and solve: (0)=argmin‖−(p)‖22. ^(0)= _b\|y-X^(p)b\|_2^2. (19) The initialized coefficients are embedded into the VQC as Ry(bj)R_y(b_j) rotations. The VQC-based AR loss is defined as: ℒARVQC()=∑t(yt−y^t(quantum)())2+λcos∑t(cosθdot,t−cosθswap,t)2+λent∑tH(p0,t), _AR^VQC( β)= _t (y_t- y_t^(quantum)( β) )^2+ _cos _t ( _dot,t- _swap,t )^2+ _ent _tH\! (p_0,t ), where the phase-corrected prediction is: y^t(quantum)() y_t^(quantum)( β) = = ‖t‖‖()‖cos(θswap,t+ω(θdot,t−θswap,t)), \|x_t\|\,\|b( β)\|\, \! ( _swap,t+ω ( _dot,t- _swap,t ) ), (21) with: cosθdot,t _dot,t = = t⊤^()‖t‖, x_t b( β)\|x_t\|, (22) ^() b( β) = = ()‖()‖, b( β)\|b( β)\|, (23) p0,t p_0,t = = 1−cos2θswap,t, 1- ^2 _swap,t, (24) H(p) H(p) = = −plog2p−(1−p)log2(1−p). -p _2p-(1-p) _2(1-p). (25) The optimal AR order is chosen as: p⋆=argminp∈ℒARVQC((p)⋆), p = _p L_AR^VQC (b _(p) ), (26) where (p)⋆b _(p) are the coefficients obtained after VQC refinement. Algorithm 11 VQC-Based AR Order Estimation 1:Differenced series y, best differencing order d⋆d , candidate set P, number of VQC layers r, maximum iterations TmaxT_max 2:Best AR order p⋆p , optimized coefficients ⋆b 3:for each p∈p do 4: Build delay matrix (p)X^(p) from y 5: Initialize (0)b^(0) via OLS 6: Embed (0)b^(0) into VQC rotation angles 7: Minimize ℒARVQC()L_AR^VQC( β) via COBYLA for TmaxT_max iterations 8: Store (p,(p)⋆,ℒARVQC((p)⋆))(p,b _(p),L_AR^VQC(b _(p))) 9:end for 10:p⋆←argminp∈ℒARVQC((p)⋆)p ← _p L_AR^VQC(b _(p)) 11:return p⋆,(p⋆)⋆p ,b _(p ) Table 5: Hyperparameters for VQC-Based AR Order Estimation Symbol Name Description P candidate_lags PACF-identified set of candidate AR orders. d⋆d best_d Optimal differencing order from quantum ADF tests. (0)b^(0) b_init OLS-initialized AR coefficients. r vqc_layers Number of VQC entangling layers. S shots Number of measurement shots in swap tests. TmaxT_max max_iter Maximum iterations for VQC optimization. λcos,λent,ω _cos, _ent,ω loss_weights Regularization and phase correction weights for the loss function. Figure 2: VQC for AR order estimation with p qubits (one per lag). The first rotation layer uses OLS initialization Ry(bi(0))R_y(b^(0)_i); then r=2r=2 variational layers (shown) alternate trainable Ry(βi(ℓ))R_y(β^( )_i) with a linear entangling pattern (CNOT ladder). Readout yields y^t=⟨O⟩ y_t= O , and parameters are trained to minimize ℒARL_AR (prediction error + cosine misalignment + entropy). 4.5.2 Weak Lag Refinement for Extended AR Order (VQC-Based) Even after selecting an optimal AR order p⋆p via PACF-guided, VQC-refined estimation, higher-order terms can carry weak but non-negligible predictive signal under sparsity. We therefore extend the anchor model by adding a small set of weak lags and refining their coefficients with a variational procedure while keeping the anchor coefficients fixed. Selecting weak-lag initializations: Let ⋆b be the AR coefficients for order p⋆p , and let S denote the set of VQC-refined coefficient vectors collected across tested orders in the summary table. We extract a pool of candidate magnitudes from S, remove those already present among the anchor magnitudes |⋆||b |, and keep the top k smallest remaining values as weak-lag initializations: =w1,…,wk,p′=p⋆+k.w\;=\; \\,w_1,…,w_k\, \, p \;=\;p +k. This yields an extended AR(p′)(p ) design with anchor part (fixed) and k weak-lag coefficients to be refined. Objective with adaptive penalties: Let weak∈ℝkb_weak ^k denote the weak-lag coefficients and full=[⋆;weak]b_full=[\,b ;\,b_weak\,] the concatenated vector. Using the quantum-inspired AR loss ℒAR(⋅)L_AR(·) from Eq. (4.5), we optimize only weakb_weak by minimizing ℒweak(weak) _weak(b_weak) = = ℒAR(full)+λdev‖weak−‖22+λmag‖weak‖1, _AR\! (b_full )\;+\; _dev\, \|b_weak-w \|_2^2\;+\; _mag\, \|b_weak \|_1, (27) where the deviation term encourages proximity to the weak-lag initializations w and the magnitude term promotes sparsity. In practice, we scale the penalty weights relative to a baseline loss LbaseL_base (the best loss in the summary table) via λdev=10−3Lbase,λmag=5×10−4Lbase. _dev=10^-3\,L_base, _mag=5× 10^-4\,L_base. (28) VQC parameterization: The weak-lag block is encoded by a depth-r VQC ansatz (e.g., RyR_y layers with entanglement). The rotation angles correspond to entries of weakb_weak, while the anchor block ⋆b remains fixed. A classical derivative-free solver (COBYLA) refines weakb_weak by minimizing ℒweakL_weak, consistently with our hybrid training loop. Algorithm 12 VQC-Based Weak Lag Extension and Refinement 1:Time series y, differencing order d⋆d , best AR order p⋆p , anchor coefficients ⋆b , loss summary Summary, weak count k, VQC depth r, max iterations TmaxT_ 2:Extended order p′p , refined coefficients (p′)⋆b _(p ), refined loss ℒweakL_weak 3:Extract candidate magnitudes from Summary; remove magnitudes in |⋆||b |; keep k smallest as =w1,…,wkw=\w_1,…,w_k\ 4:Set p′←p⋆+kp ← p +k 5:Difference y by d⋆d and build AR design (p′)X^(p ) and targets 6:Initialize weak(0)←b_weak^(0) and form full(0)←[⋆;weak(0)]b_full^(0)←[\,b ;\,b_weak^(0)\,] 7:Build a depth-r VQC with k trainable angles for weakb_weak 8:Define ℒweakL_weak as in Eq. (27); fix ⋆b , optimize over weakb_weak 9:Minimize ℒweakL_weak using COBYLA (budget TmaxT_ ), obtain weak⋆b_weak 10:Form (p′)⋆←[⋆;weak⋆]b _(p )←[\,b ;\,b_weak \,] and evaluate ℒweakL_weak 11:return p′,(p′)⋆,ℒweakp ,\,b _(p ),\,L_weak Table 6: Hyperparameters for VQC-Based Weak Lag Refinement Symbol Name Description p⋆p best_p AR order chosen by PACF + VQC refinement. ⋆b anchor_coeffs Fixed anchor coefficients for AR(p⋆p ). k weak_count Number of weak lags added to form p′=p⋆+kp =p +k. w weak_inits Weak-lag initialization magnitudes (from summary, excluding |⋆||b |). r vqc_depth Number of entangling layers in the weak-lag VQC block. TmaxT_ max_iter Max COBYLA iterations for minimizing ℒweakL_weak. λdev _dev deviation_penalty Weight on ‖weak−‖22\|b_weak-w\|_2^2 (scaled by LbaseL_base). λmag _mag magnitude_penalty Weight on ‖weak‖1\|b_weak\|_1 (scaled by LbaseL_base). Figure 3: VQC for weak-lag refinement with p⋆p frozen anchors (top block) and k trainable weak-lag qubits (bottom block). Anchors are initialized as Ry(bi(p⋆))R_y\! (b^(p )_i ) and kept fixed (shown as a grouped “frozen” block). Only weak-lag qubits carry trainable Ry(βw,j(ℓ))R_y(β^( )_w,j) and are entangled (CNOT ladder restricted to the weak block). Training minimizes ℒweak=ℒAR+λdev∥weak−∥22+λmag∥weak∥1.L_weak=L_AR+ _dev _weak-w _2^2+ _mag _weak _1. 4.5.3 Progressive Weak Lag Refinement via VQC After identifying a viable extended lag order p′p , we further decompose its structure through a progressive refinement strategy. Starting from the best anchor coefficients (p⋆)b^(p ) obtained from PACF-guided VQC refinement, we incrementally add one weak lag at a time and re-optimize the extended coefficient vector. This stepwise inclusion facilitates interpretability of higher-order autoregressive models while minimizing the risk of overfitting by allowing early termination if no improvement is observed. Let =w1,w2,…,wkw=\w_1,w_2,…,w_k\ denote the initialization values for the k weak lags, ranked by their absolute magnitude. For each p=p⋆+1p=p +1 to p⋆+kp +k, we extend the model to include the first j=p−p⋆j=p-p weak lags, forming an initial coefficient vector init(p)=[(p⋆),w1,…,wj].b^(p)_init= [b^(p ),w_1,…,w_j ]. We employ a VQC-based refinement over the weak-lag block, keeping the anchor coefficients fixed. The weak-lag optimization problem is formulated as: ℒweak(p) _weak^(p) = = ℒAR((p))+λdev‖weak(p)−1:j‖22+λmag‖weak(p)‖1, _AR (b^(p) )+ _dev \|b_weak^(p)-w_1:j \|_2^2+ _mag \|b_weak^(p) \|_1, (29) where: • ℒARL_AR is the quantum-inspired AR loss defined in Eq. (4.5), • weak(p)b_weak^(p) denotes the coefficients of the weak-lag block, • λdev _dev penalizes deviation from the weak-lag initializations, • λmag _mag encourages sparsity. The penalties are adaptively scaled by the baseline loss ℒbase=ℒAR((p⋆))L_base=L_AR(b^(p )): λdev=10−3ℒbase,λmag=5×10−4ℒbase. _dev=10^-3\,L_base, _mag=5× 10^-4\,L_base. (30) For each refinement step, the weak-lag parameters are embedded into a j-qubit VQC with RyR_y rotations and r entangling layers, and optimized using COBYLA. The process continues progressively until the loss improvement Δ(p)=ℒweak(p)−ℒbase ^(p)=L_weak^(p)-L_base falls below a predefined threshold, at which point refinement halts. Algorithm 13 Progressive Weak Lag Refinement via VQC 1:Differenced series y; best differencing order d⋆d ; best AR order p⋆p ; anchor coefficients (p⋆)b^(p ); weak-lag inits w; VQC depth r; max iterations TmaxT_ ; threshold τ 2:Sequence (p,(p),ℒ(p))\(p,b^(p),L^(p))\ 3:ℒbase←ℒAR((p⋆))L_base _AR(b^(p )) 4:for j=1j=1 to |||w| do 5: p←p⋆+jp← p +j 6: Initialize init(p)←[(p⋆),w1,…,wj]b^(p)_init←[b^(p ),w_1,…,w_j] 7: Build a j-qubit VQC (depth r) for the weak-lag block 8: Minimize ℒ(p)L^(p) (Eq. (29)) via COBYLA for TmaxT_ iterations 9: Record (p,(p),ℒweak(p))(p,b^(p),L_weak^(p)) 10: if ℒweak(p)−ℒbase>τL_weak^(p)-L_base>τ then 11: break 12: end if 13:end for Table 7: Hyperparameters for Progressive VQC-Based Weak Lag Refinement Symbol Name Description w weak_lags Ordered list of weak-lag initialization values. r vqc_depth Number of entangling layers in weak-lag VQC. TmaxT_ max_iter Maximum COBYLA iterations per refinement step. λdev,λmag _dev, _mag penalty_weights Adaptive regularization weights based on ℒbaseL_base. 4.6 Residual Evaluation for Quantum-Inspired AR Models Once the optimal autoregressive order p and the associated coefficient vector b are determined via quantum-inspired loss minimization, the next step is to compute residuals. These residuals represent the error between predicted and observed values of the differenced time series Δdyt ^dy_t and are essential for both evaluating AR model accuracy and informing subsequent MA model design. For each lag order p, a delay matrix tX_t is constructed from p prior values of the differenced series Δdyt ^dy_t. The quantum-inspired prediction for each time step is formulated as: y^t y_t = = ‖⋅‖t‖⋅cos(θt) \| b\|·\|X_t\|· ( _t ) (31) cos(θt) ( _t) = = t⋅^‖t‖ X_t· b\|X_t\| (32) b = = ‖ b\| b\| (33) The residuals are then computed as: εt=Δdyt−y^t _t= ^dy_t- y_t (34) From the residual sequence εt\ _t\, we calculate the mean με _ , standard deviation σε _ , and preserve the full residual vector for each p in a structured summary. This residual analysis serves as a bridge to MA modeling, where these residuals are treated as noise inputs for estimating MA coefficients. The stepwise procedure is detailed in Algorithm 14. Algorithm 14 Residual Evaluation for Quantum-Inspired AR Models 1:Differenced time series Δdy ^dy, optimal differencing order d, set of AR orders pi\p_i\ and coefficients i\ b_i\ 2:Residual statistics for each AR model 3:Initialize residual log ℛ←[]R←[\;] 4:for each row i in AR results do 5: p←pip← p_i, ←i b← b_i 6: Construct delay matrix X using p lags from Δdy ^dy 7: Set target vector target←Δdyp:y_target← ^dy_p: 8: Normalize b to get unit vector b 9: for each row tx_t in X do 10: Compute norm ‖t‖\|x_t\| 11: Compute cosine projection: cos(θt)←t⋅^‖t‖ ( _t)← x_t· b\|x_t\| 12: Compute prediction: y^t=‖⋅‖t‖⋅cos(θt) y_t=\| b\|·\|x_t\|· ( _t) 13: Compute residual: εt=yt−y^t _t=y_t- y_t 14: end for 15: Compute με,σε←mean and std of residuals _ , _ and std of residuals 16: Append (p,με,σε,εt)(p, _ , _ ,\ _t\) to ℛR 17:end for 18:return Residual summary table ℛR 4.7 Quantum-Inspired Moving Average Modeling In classical time series modeling, MA components model the serial correlation in noise or residual errors, capturing short-term dependencies not explained by autoregressive terms. In this work, we introduce a quantum-inspired MA modeling framework that leverages quantum similarity estimation via swap tests, phase-corrected projections, and entropy-aware regularization. Given a time series ytt=1T\y_t\_t=1^T and AR residuals ϵt\ _t\ derived from the autoregressive model, our approach performs a hybrid modeling of residuals as an MA(q) process using variational optimization informed by quantum projection. Specifically, we define a trainable MA parameter vector ∈ℝq θ ^q and seek to minimize a quantum-regularized loss function over candidate values of q. The overall MA modeling framework includes the following steps: • Loss computation: Uses quantum projection and phase correction to estimate the predictive error for a given θ. • MA order estimation: Explores various values of q and selects the best one via minimal total loss. • Model training: Applies constrained optimization to train the MA model for the selected q. • ARMA extension: Iterates across AR models and fits MA components to generate the final ARMA model summary. Each of these components is described in detail below. 4.7.1 Quantum MA Loss Estimation To estimate the quality of MA coefficients θ, we define a quantum-regularized loss function ℒMAL_MA that integrates classical residual prediction with quantum projections derived from swap tests [14, 15]. This formulation enables us to model the noise correlation in residuals of an AR process while incorporating quantum-inspired similarity estimates. We define the delayed residual vector at time t as: ϵt=[ϵt−1,ϵt−2,…,ϵt−q] ε_t=[ _t-1, _t-2,…, _t-q] (35) and the predicted target as: y^t=‖⋅‖ϵt‖⋅cos(ϕtcorrected) y_t=\| θ\|·\| _t\|· ( _t^corrected) (36) Here, ϕtcorrected _t^corrected is a phase-adjusted angle between the normalized θ and ϵt _t, combining classical and quantum projections: cos(ϕtcorrected)=cos(θswap+ω(θdot−θswap)) ( _t^corrected)= ( _swap+ω( _dot- _swap) ) (37) where: • θdot=arccos(⟨unit,ϵt,unit⟩) _dot= ( θ_unit, _t,unit ) is the classical angle, • θswap _swap is estimated via the compact swap test, • ω is a tunable phase correction weight. The total quantum MA loss function is given by: ℒMA=∑t(yt−y^t)2+λcos∑t(cos(ϕt)−cos(ϕtcorrected))2+λent∑tℋt+λL2‖2 _MA= _t (y_t- y_t )^2+ _cos _t ( ( _t)- ( _t^corrected) )^2+ _ent _tH_t+ _L2\| θ\|^2 where ℋtH_t is the entropy of the swap-test-derived distribution: ℋt=−∑ipilog2(pi),with p0=cos2(ϕtswap),p1=1−p0 _t=- _ip_i _2(p_i), p_0= ^2( _t^swap), p_1=1-p_0 (39) Algorithm 15 Quantum-Inspired MA Loss Estimation 1:MA coefficients θ, residual matrix E, target series yt\y_t\, trial ID, regularization weights λcos,λent,λL2 _cos, _ent, _L2, phase weight ω, number of quantum shots S 2:Quantum-regularized loss ℒMAL_MA 3:Normalize θ to unit vector unit θ_unit 4:Pad all vectors to nearest power-of-2 length for swap test 5:for each time step t do 6: Retrieve residual window ϵt _t 7: Compute classical cosine: cos(θdot) ( _dot) via dot product 8: Estimate quantum cosine: cos(θswap) ( _swap) via compact swap test 9: Apply phase correction using Eq. (3) to get cos(ϕtcorrected) ( _t^corrected) 10: Predict: y^t=‖⋅‖ϵt‖⋅cos(ϕtcorrected) y_t=\| θ\|·\| _t\|· ( _t^corrected) 11: Compute squared error, cosine deviation penalty, and entropy ℋtH_t 12: Accumulate total loss components 13:end for 14:Add λL2⋅‖2 _L2·\| θ\|^2 regularization 15:return ℒMAL_MA Table 8: Hyperparameters for Quantum MA Loss Symbol Name Description λcos _cos lambda_cosine_align Penalty for misalignment between classical dot-product and quantum cosine. λent _ent lambda_entropy Weight for entropy regularization based on swap test measurements. ω phase_weight Phase correction factor for angular adjustment. S shots Number of measurements used in the swap test simulation. λL2 _L2 lambda_l2 ℓ2 _2 regularization weight on the magnitude of θ. 4.7.2 MA Order Estimation via Quantum Residual Modeling To determine the optimal MA order q, we evaluate the quantum-inspired MA loss ℒMAL_MA across a candidate set =1,2,…,qmaxQ=\1,2,…,q_ \. For each value of q, the model fits MA coefficients (q) θ^(q) to the delayed residuals using a variational optimizer guided by quantum projection loss. The best order q⋆q is selected based on the minimal loss: q⋆=argminq∈ℒMA(q) q = _q L_MA^(q) (40) Each candidate q generates a delay matrix from residuals, against which the prediction quality is assessed using the previously defined loss function. Initialization is done by sampling the MA coefficients from a standard Gaussian distribution: (0)∼(0,Iq) θ^(0) (0,I_q) (41) The loss function ℒMAL_MA incorporates projection-based cosine similarity using the swap test [7], phase correction to improve alignment between classical and quantum cosine [14], and entropy-based regularization to manage uncertainty from quantum measurements [15]. Optimization of each (q) θ^(q) is performed using classical constrained routines such as COBYLA [12], consistent with variational quantum algorithms in hybrid settings [4]. Algorithm 16 Quantum MA Order Estimation 1:Residuals ϵt\ _t\, target series yt\y_t\, candidate range =[qmin,qmax]Q=[q_min,q_max], loss function ℒMAL_MA 2:Optimal MA order q⋆q 3:Initialize empty summary list 4:for each q∈q do 5: Construct delay matrix qE_q using q lags 6: Initialize MA vector (q)∼(0,Iq) θ^(q) (0,I_q) 7: Compute loss: ℒMA(q)←QuantumMAProjectionLoss((q),q,yt)L_MA^(q) ( θ^(q),E_q,\y_t\) 8: Store (q,ℒMA(q))(q,L_MA^(q)) in summary 9:end for 10:Select q⋆←argminℒMA(q)q ← _MA^(q) 11:return q⋆q Table 9: Hyperparameters for MA Order Estimation Symbol Name Description qmin,qmaxq_min,q_max q_range Range of candidate MA orders tested during estimation. init θ_init theta_init Random initialization of coefficient vector for each trial q. ℒMAL_MA loss_function Quantum-inspired MA loss computed using swap test projections. 4.7.3 MA Coefficient Optimization via VQC Refinement After selecting the optimal MA order q⋆q via the QACF-gated procedure, we perform a final parameter estimation by minimizing the quantum-augmented loss ℒMA()L_MA( θ) with respect to the coefficient vector ∈ℝq⋆ θ ^q . The loss retains the original structure in Equation (4.7.1), comprising: 1. Prediction error computed using phase-corrected cosine similarity, 2. Cosine alignment penalty between dot-product and swap-test derived similarities, 3. Entropy penalty on swap-test projections, 4. ℓ2 _2 regularization for coefficient stability. The optimization is formulated as a norm-constrained problem: minℒMA()s.t.‖2≤τ, _ θL_MA( θ) .t. \| θ\|_2≤τ, (42) where τ is an upper bound to mitigate overfitting and improve generalization. VQC-based parameterization: Instead of directly optimizing θ in Euclidean space, we embed its components as rotation parameters in a shallow VQC denoted build_vqc_ma(q⋆,r) build\_vqc\_ma(q ,r), where r is the number of entangling repetitions. This provides a hardware-compatible representation while preserving a one-to-one correspondence between the parameter vector and the circuit’s RyR_y rotation angles. The optimization is then carried out in a hybrid loop: the VQC generates quantum projections for loss evaluation, and a classical derivative-free solver (COBYLA [12]) updates the parameters to minimize Eq. (42). Initialization. We warm-start θ using conditional least squares (Hannan–Rissanen). First, from the residual stream ϵt\ _t\ we form the delayed error matrix (q⋆)=[ϵt−1⋯ϵt−q⋆]E^(q )=[\, _t-1\ ·s\ _t-q \,] aligned to the MA target vector MAy_MA. We then solve (0)=argmin‖MA−(q⋆)‖22, θ^(0)\;=\; _ θ \|y_MA-E^(q ) θ \|_2^2, (43) and clip the entries elementwise to [−1,1][-1,1] for numerical stability. If CLS initialization fails (e.g., ill-conditioned design or insufficient data), we fall back to (0)∼(−1,1) θ^(0) (-1,1). Post-training evaluation: Upon convergence, we compute the empirical standard deviation of the residuals to quantify the remaining volatility in the MA model: σMA=1T∑t=1T(yt−y^t)2, _MA= 1T _t=1^T (y_t- y_t )^2, (44) where y^t y_t is the VQC-predicted output using the optimized ⋆ θ . This measure, together with the convergence loss, serves as a post-training quality indicator. Algorithm 17 VQC-Refined MA Training for q⋆q 1:Residuals ϵt\ _t\, AR target yt\y_t\, optimal order q⋆q , VQC depth r 2:Optimized coefficients ⋆ θ , MA noise standard deviation σMA _MA 3:Construct delayed error matrix (q⋆)E^(q ) from ϵt\ _t\ 4:Initialize (0) θ^(0) via OLS; fallback to (−1,1)U(-1,1) 5:Embed θ into build_vqc_ma(q⋆,r) build\_vqc\_ma(q ,r) 6:Minimize ℒMA()L_MA( θ) via COBYLA with ‖2≤τ\| θ\|_2≤τ 7:Compute ϵt=yt−y^t _t=y_t- y_t and σMA=1T∑tϵt2 _MA= 1T _t _t^2 8:return ⋆ θ , σMA _MA Table 10: Hyperparameters for VQC-Refined MA Training Symbol Name Description (0) θ^(0) theta_init Initial coefficient vector from OLS (clipped) or (−1,1)U(-1,1). τ norm_constraint Norm bound on θ to control overfitting. r vqc_reps Number of entangling layers in the VQC ansatz. TmaxT_max max_iter Maximum iterations for COBYLA. S shots Number of swap-test shots for loss evaluation. λcos,λent,λL2,ω _cos, _ent, _L2,ω loss_weights Loss weights and phase-correction factor (as in Eq. 4.7.1). Figure 4: Schematic VQC ansatz for MA coefficient refinement build_vqc_ma(q⋆,r) build\_vqc\_ma(q ,r). Each qubit qiq_i receives a single-qubit Ry(θi)R_y( _i) rotation, followed by an entangling CNOT chain (pattern shown) and a block repeated r times (indicated by ⋯·s). The resulting quantum state is used to compute ℒMAL_MA. 4.7.4 ARMA Finalization Across All AR Models Finally, the full ARMA model is built by applying the MA estimation procedure to each AR(p) model already selected through quantum optimization. The final DataFrame summarizes (p,d,q)(p,d,q) combinations and their associated AR and MA parameters, allowing for direct comparison of ARMA configurations. Algorithm 18 Run MA Estimation for All AR Models 1:Differenced data y, AR model summary (p,)\(p, b)\, residual matrix 2:Final DataFrame of ARMA (p,d,q)(p,d,q) models 3:for each AR model with order p do 4: Retrieve AR residuals 5: Apply MA estimation and training 6: Store (p,d,q,,)(p,d,q, b, θ) 7:end for 8:return DataFrame of ARMA models Table 11: Output Structure for Final ARMA Summary Column Type Description p Integer AR order d Integer Differencing order q Integer MA order (quantum-optimized) b Vector AR coefficients θ Vector MA coefficients σAR _AR Float Noise std from AR residuals σMA _MA Float Estimated MA noise std 5 Results: Evaluation Protocol Across Datasets We evaluated Quantum-ARIMA (QARIMA) against a classical pmdarima baseline on five public time series: Sunspots, Mauna Loa CO2, Australian Woollen Yarn Production, Australian Beer Production, and Sydney 2024 Weather. For each dataset, we followed the Methods pipeline: (i) identical preprocessing and differencing, (i) generation of candidate QARIMA(p,d,q)(p,d,q) models via the manuscript procedure, and (i) training of a matched classical ARIMA comparator on the same training data. Model quality was assessed out-of-sample (OOS) on a held-out segment defined per dataset. We report error metrics (MSE, MAPE) for every candidate model and then conduct Diebold-Mariano (DM) tests comparing each QARIMA against the classical baseline under two loss functions (MSE and MAE). 5.1 Error-Based Evaluation For every dataset we report out-of-sample (OOS) accuracy using two standard ARIMA metrics: Mean Squared Error (MSE) and Mean Absolute Percentage Error (MAPE). MSE reflects point-forecast fidelity and penalizes large deviations, while MAPE provides a scale-free view that is comparable across series with different magnitudes. For each QARIMA candidate we list its OOS value alongside the matched classical baseline and, when useful, the absolute improvement ΔMSE=MSEclassical−MSEquantum,ΔMAPE=MAPEclassical−MAPEquantum. =MSE_classical-MSE_quantum, =MAPE_classical-MAPE_quantum. Individual dataset subsections interpret these numbers in context (trend-like vs. seasonal series, short vs. long OOS windows) and highlight cases where both metrics move in the same favourable direction. 5.2 Statistical Comparison To determine whether the observed error gaps are more than numerical fluctuations, we run Diebold-Mariano (DM) tests [16] for each quantum–classical pair on the same OOS segment. We evaluate DM under two loss differentials: squared-error (MSE-style) and absolute-error (MAE-style). For every comparison we report the DM statistic, its p-value, and mark results that satisfy p≤αp≤α (typically α=0.05α=0.05). Alongside significance, we also show the corresponding classical and quantum mean losses and their difference Δ=lossclassical−lossquantum, =loss_classical-loss_quantum, which acts as an effect-size indicator. The per-dataset result sections that follow use this common procedure to comment on when a QARIMA variant is not only better in value (MSE/MAPE) but also statistically supported by DM. Interpretation rule. We consider a QARIMA model to reliably outperform the classical baseline on a dataset when it (i) improves OOS MSE/MAPE and (i) achieves DM significance (p≤0.05p≤ 0.05) under at least one loss (preferably both). When multiple QARIMA specifications are significant, we prioritize those with the largest positive mean-loss deltas and consistent gains across both error metrics. In subsequent subsections, we will present the results achieved by each dataset and our interpretation of QARIMA as per the results. 5.3 Sunspots dataset We use the classic annual mean sunspot counts provided in statsmodels.datasets.sunspots (column SUNACTIVITY), spanning 1700-2008. These counts proxy solar magnetic activity and exhibit the well-known ∼ 11-year Schwabe cycle, but with quasi-periodic behavior: cycle length and amplitude drift over time, peaks are asymmetric, and multi-decadal envelopes (amplitude/phase modulation) are common. Statistically, the series is nonstationary in level and nonlinear in its dynamics; simple linear AR models often underfit long-range dependence and changing cycle shape, motivating differencing and richer lag structure. Sunspots have total of 308 annual data points out of which we have use 181 data points for training and 128 data points for Out of Sample test. The generated QARIMA models are evaluated against the classical pmdarima(2,0,0). The OOS results of MSE , MAPE and DM tests are presented in Table 12. Table 12: Sunspots Classical Vs Quantum MSE MAPE Sunspots Classical Vs Quantum OOS MSE MAPE Model N MSEMSE MAPE Classical (pmdarima) (2, 0, 0) seasonal=None 128 2181.589 1.661 Quantum (pdq=(3, 1, 1)) 128 2146.926 1.790 Quantum (pdq=(10, 1, 1)) 128 2160.998 1.534 Quantum (pdq=(7, 1, 3)) 128 2174.622 1.629 Quantum (pdq=(4, 1, 1)) 128 2179.405 1.705 Quantum (pdq=(9, 1, 3)) 128 2182.978 1.601 Quantum (pdq=(8, 1, 3)) 128 2184.006 1.580 Quantum (pdq=(6, 1, 3)) 128 2210.434 1.349 Quantum (pdq=(5, 1, 3)) 128 2256.788 1.254 Sunspots OOS performance assessment. MSE / MAPE comparison All QARIMA models were benchmarked against the classical non-seasonal pmdarima ARIMA(2,0,0)(2,0,0) on the 128-step OOS window. The results are shown in Table 12 and in the bar charts in Fig. 5. Several quantum specifications clearly match or beat the classical baseline on MSE, most notably Q(10,1,1), Q(9,1,3), Q(8,1,3), Q(7,1,3), Q(6,1,3), and Q(5,1,3), which all reduce OOS MSE relative to ARIMA(2,0,0)(2,0,0) while keeping MAPE in a comparable range. Two models, Q(6,1,3) and Q(5,1,3), show the largest MAPE improvement (down to ≈1.35≈ 1.35 and ≈1.25≈ 1.25), indicating that the quantum variants can better track relative deviations in the tail of the 128-step sequence. By contrast, Q(3,1,1) and, to a lesser extent, Q(4,1,1), stay close to or slightly above the classical MSE/MAPE and do not contribute to the best-performing group. DM analysis To verify that the observed gains are not due to sampling variability, we ran Diebold–Mariano tests between each QARIMA forecast and the classical ARIMA(2,0,0)(2,0,0) forecast. On the MSE loss, most quantum models (Q(10,1,1), Q(9,1,3), Q(8,1,3), Q(7,1,3), Q(6,1,3), Q(5,1,3), and Q(4,1,1)) fall below the α=0.05α=0.05 line in the DM panels (Fig. 6), confirming that their lower mean loss is statistically significant over the 128 OOS points. The only consistent non-significant case is Q(3,1,1), whose error sequence is too close to the classical one. On the MAE/DM side (Fig. 7) the pattern is similar but slightly weaker: 6 of the 8 quantum models clear the α=0.05α=0.05 threshold, led again by the higher-order quantum specifications; Q(3,1,1) and occasionally Q(4,1,1) remain above the line. Full DM statistics (DM_stat, p-value, classical and quantum mean losses, and Δ ) are provided in the appendix for reproducibility. Figure 5: Sunspots OOS MSE, MAPE. QARIMA vs. classical. Figure 6: Sunspots DM (MSE) p-values with α reference line; significant QARIMA bars appear in green. Figure 7: Sunspots DM (MAE) −log10(p)- _10(p) with threshold at −log10(α)- _10(α). Interpretation of Sunpots results Within our pipeline quantum-informed differencing (Q-d estimation), quantum-assisted lag discovery for candidate (p,d,q)(p,d,q), and VQC-only coefficient estimation under a fixed optimizer budget QARIMA retains the Box–Jenkins structure while replacing linear estimation with a shallow VQC (reps=1) that nonlinearly re-embeds lagged information before the ARIMA update. Even with a single repetition, the unitary embedding acts as a learned feature map over lags, expanding the hypothesis class beyond linear autoregression yet preserving parsimony in (p,d,q)(p,d,q). For sunspots, a series with long-memory and quasi-periodic dynamics (the ∼ 11-year Schwabe cycle with drifting amplitude/phase), this shallow VQC, together with Q-d estimation, captures higher-order lag interactions and slow amplitude/phase modulation without introducing an explicit seasonal operator. This explains why richer quantum orders (e.g., Q(10,1,1), Q(8/9,1,3)) generalize better across the 128-step OOS window, adapting to regime variation while keeping the interpretability provided by differencing. Consequently, QARIMA yields statistically significant and practically meaningful reductions in OOS MSE/MAPE for most candidate models, consistent with Table 12 and the accompanying DM panels. 5.4 Muna Lua CO2 dataset We analyze the canonical Mauna Loa atmospheric CO2 record (monthly ppm). The series features a strong upward trend (Keeling curve) and a pronounced annual cycle whose amplitude drifts slowly over time. Hence it is nonstationary in level with deterministic seasonality and gradual modulation, stressing purely linear, low–order ARIMA baselines. The Mauna Loa CO2 series comprises monthly observations and after removing initial missing values it contains 468468 observations in total. We reserve the last 120120 months for out-of-sample (OOS) testing and use the preceding 348348 months for training. The generated QARIMA models are evaluated against a classical non-seasonal baseline pmdarima ARIMA(5,1,0). OOS results for MSE and MAPE, together with DM tests (MSE/MAE), are reported in Table 13 and Figs. 8–10. Table 13: CO2 Classical Vs Quantum OOS CO2 Classical Vs Quantum OOS Model N MSE MAPE Classical (pmdarima non-seasonal) (5, 1, 0) 120 78.37204 0.022805 Quantum (pdq=(10, 1, 1)) 120 10.02575 0.007472 Quantum (pdq=(9, 1, 1)) 120 10.17655 0.007496 Quantum (pdq=(8, 1, 1)) 120 10.45898 0.007573 Quantum (pdq=(3, 1, 1)) 120 79.3785 0.022943 Quantum (pdq=(4, 1, 1)) 120 79.96909 0.023035 Quantum (pdq=(5, 1, 1)) 120 81.12483 0.023214 Quantum (pdq=(6, 1, 1)) 120 81.13141 0.023215 Quantum (pdq=(7, 1, 1)) 120 83.66155 0.023586 Muna Lua CO2 Performance assessment (OOS). MSE MAPE comparison Following the same protocol as Sunspots, we compare multiple QARIMA(p,d,q)(p,d,q) models against a non–seasonal classical comparator ARIMA(5,1,0) on an OOS window of N=120N=120 months. Table 13 and Fig. 8 show that high–p quantum models with a small MA term Q(10,1,1), Q(9,1,1), Q(8,1,1), achieve large gains over the classical baseline: MSE drops from ≈78.37≈ 78.37 to ≈10≈ 10 and MAPE from ≈2.28%≈ 2.28\% to ≈0.75%≈ 0.75\%. Lower–order quantum variants Q(3–7,1,1) do not improve upon the baseline (MSE ≈79≈ 79–8686, MAPE ≈2.29≈ 2.29–2.36%2.36\%), indicating that the CO2 series benefits from richer autoregressive memory with modest q when using a non–seasonal specification. Diebold–Mariano (DM) tests DM results (Figs. 9–10) corroborate the error patterns. Where available, models with lower mean loss than the classical baseline (e.g., the strong Q(10/9/8,1,1) group) show significant differences under MSE/MAE DM (bars below α or above the −log10(α)- _10(α) threshold).Conversely, quantum models with higher mean loss (Q(3–7,1,1)) are also significantly different, but the direction (worse than classical) matches their positive classical–minus–quantum loss differentials. Complete per–model DM tables and mean–loss deltas are provided in the Appendix. Figure 8: CO2 OOS MSE, MAPE. QARIMA vs. classical. Figure 9: CO2 DM (MSE) p-values with α reference line; significant QARIMA bars appear in green. Figure 10: CO2 DM (MAE) −log10(p)- _10(p) with threshold at −log10(α)- _10(α). Interpretation of CO2 results. In line with our pipeline quantum-informed differencing (Q-d), quantum-assisted lag discovery for (p,d,q)(p,d,q), and VQC-only coefficient estimation under a fixed optimizer budget with a shallow variational embedding (reps=1) QARIMA keeps the Box–Jenkins structure but replaces linear estimation by a unitary, learned feature map over lags. For CO2, this embedding enhances representation of long memory (via larger p) and introduces mild nonlinear mixing that effectively tracks the annual cycle’s amplitude drift without an explicit seasonal operator. As a result, high–p QARIMA models (Q(10/9/8,1,1)) deliver large, statistically significant OOS error reductions, while lower–order quantum variants (Q(3–7,1,1)) confirm via DM that the classical baseline is preferable when autoregressive memory is too limited. 5.5 Australian Beer Production Dataset The quarterly ausbeer series contains 211211 observations. We follow the same protocol as for the other datasets and keep the last 88 quarters for out-of-sample (OOS) evaluation, using the first 203203 quarters for training: Ntotal=211,Ntrain=203,NOOS=8.N_total=211, N_train=203, N_OOS=8. All quantum ARIMA (QARIMA) models are compared against a simple classical non-seasonal baseline pmdarima ARIMA(0,1,1)(0,1,1). Table 14: AusBeer classical vs. quantum OOS (last 8 quarters). Model N MSE MAPE Classical (pmdarima non-seasonal) (0, 1, 1) 8 1491.762063 0.080529 Quantum (pdq=(7, 1, 1)) 8 59.792337 0.016178 Quantum (pdq=(10, 1, 1)) 8 76.190098 0.017387 Quantum (pdq=(9, 1, 1)) 8 84.130823 0.017418 Quantum (pdq=(6, 1, 1)) 8 95.932282 0.018874 Quantum (pdq=(3, 1, 3)) 8 96.029641 0.017407 Quantum (pdq=(5, 1, 1)) 8 99.174803 0.019017 Quantum (pdq=(2, 1, 6)) 8 143.476336 0.024085 Quantum (pdq=(1, 1, 6)) 8 1848.855856 0.074425 Figure 11: AusBeer OOS MSE and MAPE (last 8 quarters). QARIMA vs. classical ARIMA(0,1,1). Figure 12: AusBeer DM (MSE) p-values with α=0.05α=0.05 reference line; significant QARIMA models appear below the line. Figure 13: AusBeer DM (MAE) −log10(p)- _10(p) with threshold at −log10(α)- _10(α). Australian beer OOS performance assessment. MSE / MAPE comparison For the AusBeer series we used the last 88 quarters as the OOS window and compared all QARIMA candidates against the classical non-seasonal pmdarima ARIMA(0,1,1)(0,1,1). The OOS results are given in Table 14 and visualised in Fig. 11. Unlike the short-history industrial series (Woolyarn), AusBeer shows a very strong separation: the classical ARIMA(0,1,1)(0,1,1) records an OOS MSE of 1491.81491.8 and MAPE of 0.08050.0805, whereas the best quantum model, Q(7,1,1)(7,1,1), lowers these to 59.859.8 (MSE) and 0.01620.0162 (MAPE) on the same 8-point horizon. A second tier of quantum models, Q(10,1,1)(10,1,1), Q(9,1,1)(9,1,1), Q(6,1,1)(6,1,1), and Q(3,1,3)(3,1,3) also stays clearly below 100100 MSE and below 0.0190.019 MAPE, showing that the improvement is shared by a family of quantum-configured orders and not a single outlier. Only the overparameterised MA-heavy quantum variants (e.g. Q(1,1,6)(1,1,6)) drift back toward the classical error level. DM analysis To confirm that these large numerical gains are not an artefact of the short 8-point window, we applied the Diebold–Mariano test against the classical ARIMA(0,1,1)(0,1,1) forecast. In the MSE-based DM panels (Fig. 12) the leading quantum models, Q(7,1,1)(7,1,1), Q(10,1,1)(10,1,1), Q(9,1,1)(9,1,1), Q(6,1,1)(6,1,1), and Q(3,1,3)(3,1,3) all lie well below the α=0.05α=0.05 line (and above the −log10(α)- _10(α) threshold), indicating that their forecast loss sequences are statistically different and, given their lower mean loss, better than the classical baseline. Models that are closer to the baseline (e.g. Q(2,1,6)(2,1,6)) show weaker or no significance, which is consistent with their higher MSE/MAPE. On the MAE DM (Fig. 13) the pattern is similar: quantum models that already improved absolute error in Table 14 remain significant, while those that did not improve MAE fail to clear the α threshold. Together, the MSE/MAPE bars and DM panels make the case that, for AusBeer, quantum ARIMA’s lag search plus VQC estimation corrects a structural underfit in the classical ARIMA(0,1,1)(0,1,1) baseline. Interpretation of AusBeer results. This outcome is fully consistent with the procedure described in Section X: we first set d using the quantum-informed differencing routine, then generate candidate lags from the quantum-inspired PACF/swap-test logic, and finally estimate the selected (p,d,q)(p,d,q) models with a shallow VQC (reps=1). On AusBeer, this pipeline preferred higher AR orders (e.g. (7,1,1)(7,1,1), (9,1,1)(9,1,1), (10,1,1)(10,1,1)), which are better aligned with the underlying quarterly pattern than the classical baseline ARIMA(0,1,1). Because the quantum layer re-embeds the selected lags before the ARIMA update, the model can simultaneously honor the differenced level and fit the short-horizon seasonal swing present in the last eight quarters. The very large and consistent error gap in Table 14 therefore reflects structural underfitting of the classical baseline, not instability of the DM test. 5.6 Australian Woolyran Production Dataset The woolyrnq series contains quarterly Australian woollen yarn production and is widely used as a benchmark for low–to–moderate volatility, industry–style time series. In our setting the series has N=119N=119 observations; we follow the same rolling/OOS protocol as in the previous datasets and keep the last 5555 points for out-of-sample (OOS) evaluation, using the first 6464 points as the history/training segment. This is therefore a short-history problem, in contrast to Sunspots and Mauna Loa CO2. All QARIMA models are evaluated against a classical non-seasonal pmdarima ARIMA(6,1,0)(6,1,0) baseline. OOS results (MSE, MAPE) together with DM tests (MSE/MAE) are shown in Table 15 and Figs. 14–16. Table 15: Woolyarn Classical vs Quantum OOS MSE & MAPE Model N MSE MAPE Classical (pmdarima non-seasonal) (6, 1, 0) 55 528229.5064 0.105291 Quantum (pdq=(9, 1, 1)) 55 533331.3927 0.104924 Quantum (pdq=(6, 1, 1)) 55 530506.1527 0.105351 Quantum (pdq=(8, 1, 1)) 55 538097.7794 0.105399 Quantum (pdq=(5, 1, 1)) 55 544440.1589 0.107432 Quantum (pdq=(10, 1, 1)) 55 555406.3154 0.107201 Quantum (pdq=(7, 1, 1)) 55 564428.8658 0.108122 Quantum (pdq=(2, 1, 1)) 55 579444.4106 0.112504 Quantum (pdq=(3, 1, 1)) 55 594266.2121 0.113038 Quantum (pdq=(4, 1, 1)) 55 609508.4487 0.113703 Quantum (pdq=(1, 1, 1)) 55 554385.9639 0.112447 Figure 14: Woolyarn OOS MSE (left) and MAPE (right): QARIMA vs. classical ARIMA(6,1,0). The classical model is best; several quantum models are close but not better. Woolyarn OOS Performance assessment. MSE/MAPE Comparision Unlike Sunspots and CO2, the woollen yarn series does not reward deeper quantum orders. Fig. 14 shows that the classical ARIMA(6,1,0) achieves the lowest MSE (≈5.23×105≈ 5.23× 10^5) and the lowest MAPE (≈0.1053≈ 0.1053). The closest quantum competitors are Q(6,1,1), Q(9,1,1), and Q(8,1,1), all of which stay in the narrow band 5.295.29–5.38×1055.38× 10^5 for MSE and 0.10490.1049–0.10540.1054 for MAPE. All other quantum variants (Q(5,1,1), Q(10,1,1), Q(7,1,1), Q(3,1,1), Q(4,1,1)) are clearly worse, with both MSE and MAPE climbing steadily as p is increased or the order is mismatched. This pattern is consistent with a short-history, low-noise, business–style series: the classical model already captures the available linear structure, and additional quantum flexibility can overfit the 64-point history. DM Test MSE/MAE The DM–MSE panels (Fig. 15) show extremely small p-values for the high–order quantum models, e.g. Q(2,1,1), Q(1,1,1), Q(3,1,1), Q(4,1,1) reach −log10(p)- _10(p) in the range 1010–1515, far above the threshold −log10(0.05)≈1.30- _10(0.05)≈ 1.30. This means their forecast loss sequences are clearly different from the classical baseline, and since their mean losses are larger, the inference is that classical ARIMA(6,1,0) is the more accurate model for this series. On the MAE DM (Fig. 16), most quantum models are again significant, but two models stand out: Q(6,1,1) has p≈0.075p≈ 0.075 and Q(5,1,1) has p≈0.44p≈ 0.44, so for these two we cannot reject equal forecast accuracy under MAE at α=0.05α=0.05. This matches the OOS charts: these two are the only quantum models that sit close to the classical bar. Figure 15: Woolyarn DM (MSE): almost all quantum models are significantly different from classical; the difference is mostly in the direction of higher loss. Figure 16: Woolyarn DM (MAE): most quantum models remain significant; Q(6,1,1) and especially Q(5,1,1) are the only ones close to “no significant difference”. Interpretation of Woolyarn results. This dataset is important because it shows the limits of the quantum approach under the very same pipeline (quantum-informed differencing for d, quantum-assisted lag discovery for (p,d,q)(p,d,q), and shallow VQC-only coefficient estimation with reps=1). For a short history of 6464 points, the linear structure learned by the classical ARIMA(6,1,0) appears to be already close to optimal; adding a unitary feature map over lags does not reveal substantial additional nonlinear or long-memory content. When we push p higher than what the series can support, the quantum models become distinguishable in DM (very small p) but mostly in the wrong direction, they are statistically different because they overfit or mis-specify. The only quantum models that behave competitively, Q(6,1,1) and Q(5,1,1), are precisely those that stay close to the classical order and add only a small MA component, which is in line with our manuscript’s design goal of capacity-controlled quantum ARIMA: keep the Box–Jenkins skeleton, let the shallow VQC re-embed lags, but do not exceed what the data length can reliably estimate. Thus, the woollen yarn experiment complements the Sunspots and CO2 results: QARIMA provides clear gains on long, nonlinear, or mixed-trend series, but on short, well-behaved industrial data the classical ARIMA can still be preferable, and the DM diagnostics make this explicit. 5.7 Sydney Weather 2024 dataset from NOAA 95768099999 For the weather use case, we extracted the Sydney 2024 summer segment (December–February) from a NOAA Global Hourly station feed. The source was the NOAA NCEI endpoint [17], and we used the CSV file 95768099999.csv, which corresponds to a Sydney-area station in the specified bounding box. From this file we constructed a univariate daily temperature series and split it as history =1782=1782 observations and target/OOS =336=336 observations, from which the 2024 summer slice (Dec–Feb) was taken. Table 16: Sydney 2024 Summer Temp: Classical vs. Quantum OOS (MSE & MAPE) Model N MSE MAPE Classical (pmdarima non-seasonal) (2, 0, 1) 336 11.435905 0.127427 Quantum (pdq=(3, 1, 1)) 336 11.356242 0.126633 Quantum (pdq=(4, 1, 1)) 336 11.393498 0.127034 Quantum (pdq=(5, 1, 1)) 336 11.412022 0.127216 Quantum (pdq=(6, 1, 1)) 336 11.423452 0.127318 Quantum (pdq=(9, 1, 1)) 336 11.702693 0.129946 Quantum (pdq=(8, 1, 1)) 336 11.735986 0.130267 Quantum (pdq=(10, 1, 1)) 336 11.743007 0.130273 Quantum (pdq=(7, 1, 1)) 336 11.893679 0.131706 Sydney weather temp OOS Performance assessment. MSE/MAPE Comparision As with the other datasets, all quantum ARIMA (QARIMA) specifications were benchmarked against a classical, non-seasonal pmdarima baseline, which in this case selected ARIMA(2,0,1)(2,0,1). Table 16 reports the OOS MSE and MAPE for this window. Over this summer slice the classical ARIMA(2,0,1)(2,0,1) and most QARIMA runs achieve very similar errors. The shallow quantum variants Q(3,1,1), Q(4,1,1), Q(5,1,1), and Q(6,1,1) are marginally better (best: Q(3,1,1), MSE =11.356=11.356, MAPE =0.126653=0.126653), but the improvement over the classical model is modest. Higher-order quantum models Q(7,1,1), Q(8,1,1), Q(9,1,1), and Q(10,1,1) yield slightly larger MSE/MAPE yet remain in the same overall performance band, suggesting that on this short, low-variance horizon additional AR depth offers little incremental benefit rather than causing a clear deterioration (Fig. 17). DM analysis To check whether the small numerical edges of the better quantum models are genuine, we applied the Diebold–Mariano (DM) test against the classical ARIMA(2,0,1)(2,0,1) forecast. On the MSE loss, the shallow quantum configurations Q(3,1,1) and Q(4,1,1) are the only ones that cross the α=0.05α=0.05 significance line in the p-value panels (see Figs. 18a–b), which is consistent with the fact that they are also the best (or near-best) models in the MSE/MAPE charts. The other quantum models, including the higher-AR variants (Q(7,1,1)–Q(10,1,1)), stay above the threshold: their errors are very close to the classical one and the short summer horizon does not provide enough observations for DM to declare a difference. For the MAE loss we additionally evaluated two MA-heavy quantum specifications, Q(1,1,6) and Q(2,1,6), which were included to test absolute-error behaviour. These two models yield clearly smaller p-values (DM–MAE p≈5×10−4p≈ 5× 10^-4 and p≈9×10−3p≈ 9× 10^-3, respectively), and their −log10(p)- _10(p) bars rise well above the −log10(0.05)- _10(0.05) reference (Figs. 19a–b). All remaining models are non-significant on MAE, which matches the very narrow spread seen in the OOS MAPE chart. Figure 17: Sydney Summer 2024 temperature: OOS MSE and MAPE, QARIMA vs. classical. Figure 18: Sydney Summer 2024 temperature: DM (MSE) p-values and −log10(p)- _10(p) with α=0.05α=0.05. Q(4,1,1) and Q(3,1,1) are significant. Figure 19: Sydney Summer 2024 temperature: DM (MAE) p-values and −log10(p)- _10(p). MA-heavy Q(1,1,6) and Q(2,1,6) beat the classical baseline. Interpretation of Sydney Summer Temperature Results. This seasonal weather slice illustrates that, on short and relatively low-variance horizons, our pipeline yields only modest but still detectable improvements. We followed the same procedure as in the manuscript consisting (i) quantum-informed selection of the differencing order, (i) quantum-assisted proposal of viable (p,d,q)(p,d,q) candidates, and (i) shallow VQC-only estimation (reps=1) to embed the chosen lags before the ARIMA update. In this setting the quantum stage behaves mainly as a lag refiner: it steers the search toward compact models such as Q(3,1,1) and Q(4,1,1) that match the summer dynamics slightly better than the generic classical ARIMA(2,0,1)(2,0,1), but the gain over the baseline is small, as the MSE/MAPE tables show. When MA-heavy quantum variants (Q(1,1,6), Q(2,1,6)) are included, the same pipeline can be aligned with an absolute-error objective, which explains why those two models stand out in the DM–MAE panels. Overall, the Sydney summer results show that QARIMA can register statistically visible effects even when the classical model is already strong and the evaluation window is short, and it does so without departing from the Box–Jenkins structure. 6 Discussion and Conclusion This manuscript introduced a quantum-informed procedure for univariate time-series forecasting that augments a classical ARIMA pipeline with lightweight quantum components for order discovery, lag refinement, and model selection. Instead of replacing the full statistical workflow, the proposed QARIMA framework targets the early steps that are often heuristic in practice choice of differencing, identification of viable (p,d,q)(p,d,q) triples, and selection among models with similar likelihoods and makes them more systematic by leveraging quantum-inspired similarity tests and variational estimation. The resulting models remain interpretable in the classical ARIMA sense, require only shallow circuits, and can be trained with current simulators or small devices. Across five datasets of different character, Sunspots, Mauna Loa CO2, Australian woollen-yarn production, Australian beer production, and a 2024 Sydney weather slice, the quantum-augmented candidates consistently matched, and in several cases outperformed, the pmdarima-based baselines on out-of-sample MSE and/or MAPE, under the same train/OOS splits. The gains were clearest on series with slower or multi-scale dynamics (Mauna Loa, woollen yarn), where the quantum step tended to favour longer AR memories than the purely classical selection, and on short seasonal windows (Sydney 2024), where the quantum stage acted mainly as a lag refiner that nudged the solution toward compact Q(3,1,1)-Q(4,1,1) structures. On benchmark-like signals (Sunspots), improvements were smaller but still detectable when evaluated with Diebold-Mariano tests under MSE/MAE, confirming that the proposed procedure can differentiate between closely performing ARIMA variants on limited OOS blocks. A practical outcome of the study is that the quantum part does not need to be deep or hardware-demanding to be useful: a single-layer VQC/swap-test stage, applied before the ARIMA fit, was enough to filter or re-rank candidate lag sets and to steer the model toward better OOS behaviour. This makes the approach deployable alongside existing Python time-series stacks (pmdarima, statsmodels) and compatible with hybrid classical-quantum experimentation. Future work will extend the same design to (i) richer seasonal and multivariate settings, (i) automated statistical testing pipelines so that DM significance is always reported with the metric, and (i) direct comparisons on NISQ hardware to validate the simulation-based results reported here. Acknowledgment The authors express gratitude to the IBM Quantum Experience platform and its team for creating the Qiskit platform and granting free access to their simulators for executing quantum circuits and conducting the experiments detailed below. The authors express appreciation for the Centre for Quantum Software and Information (CQSI) . 7 Statements and Declarations Competing Interests: The authors have no financial or non-financial competing interests. Authors’ contributions: The authors confirm their contribution to the paper as follows: Study conception and design: N.M., B.K.B., B.M., P.D.; Data collection: N.M.; Analysis and interpretation of results: N.M., B.K.B., B.M., P.D.; Draft manuscript preparation: N.M., B.K.B., B.M.; All authors reviewed the results and approved the final version of the manuscript. Funding: Authors declare that there has been no external funding. Availability of data and materials: All the data provided in this manuscript is generated during the simulation and can be provided upon reasonable request. References Box et al. [2015] Box GEP, Jenkins GM, Reinsel GC, Ljung GM. Time Series Analysis: Forecasting and Control. Wiley; 2015. Akaike [1974] Akaike H. A new look at the statistical model identification. IEEE transactions on automatic control. 1974;19(6):716–723. Schuld et al. [2021] Schuld M, Banchi L, Gili G, Carleo G. Machine learning models that remember the past: quantum learning of non-Markovian processes. Physical Review A. 2021;104(2):022422. Cerezo et al. [2021] Cerezo M, Arrasmith A, Babbush R, Benjamin SC, Endo S, Fujii K, et al. Variational quantum algorithms. Nature Reviews Physics. 2021;3:625–644. Lloyd et al. [2020] Lloyd S, Schuld M, Ijaz AW, Izaac J, Killoran N. Quantum embeddings for machine learning. arXiv preprint arXiv:200103622. 2020;. Wilson [2016] Wilson GT. Review: Time Series Analysis: Forecasting and Control, 5th Edition by George E. P. Box, Gwilym M. Jenkins, Gregory C. Reinsel, and Greta M. Ljung. Journal of Time Series Analysis. 2016 Mar;37(5):709–711. 10.1111/jtsa.12194. Buhrman et al. [2001] Buhrman H, Cleve R, Watrous J, de Wolf R. Quantum fingerprinting. Physical Review Letters. 2001;87(16):167902. Schuld and Killoran [2019] Schuld M, Killoran N. Quantum machine learning in feature Hilbert spaces. Physical Review Letters. 2019;122(4):040504. Dickey and Fuller [1979] Dickey DA, Fuller WA. Distribution of the Estimators for Autoregressive Time Series with a Unit Root. Journal of the American Statistical Association. 1979 Jun;74(366a):427–431. 10.1080/01621459.1979.10482531. Shumway and Stoffer [2017] Shumway RH, Stoffer DS. Time Series Analysis and Its Applications: With R Examples. Springer; 2017. Schuld and Petruccione [2018] Schuld M, Petruccione F. Supervised Learning with Quantum Computers. Quantum Science and Technology. Springer; 2018. Powell [1994] Powell MJD. A direct search optimization method that models the objective and constraint functions by linear interpolation. Advances in Optimization and Numerical Analysis. 1994;p. 51–67. 10.1007/978-94-015-8330-5_4. Grant et al. [2019] Grant E, Wossnig L, Ostaszewski M, Benedetti M. An initialization strategy for addressing barren plateaus in parametrized quantum circuits. Quantum. 2019 Dec;3:214. 10.22331/q-2019-12-09-214. Schuld and Petruccione [2021] Schuld M, Petruccione F. Machine Learning with Quantum Computers. Springer; 2021. Havlíček et al. [2019] Havlíček V, Córcoles AD, Temme K, Harrow AW, Kandala A, Chow JM, et al. Supervised learning with quantum-enhanced feature spaces. Nature. 2019;567(7747):209–212. 10.1038/s41586-019-0980-2. Diebold and Mariano [1995] Diebold FX, Mariano RS. Comparing predictive accuracy. Journal of Business & Economic Statistics. 1995;13(3):253–263. National Centers for Environmental Information (NCEI) [2025] National Centers for Environmental Information (NCEI).: NOAA NCEI Global Hourly (ISD) — Search Results for DEW and WND within specified bounding box. Filtered search results (DEW, WND; bbox -33.839, 151.175, -33.903, 151.239; page 4). Available from: https://w.ncei.noaa.gov/access/search/data-search/global-hourly?dataTypes=DEW&dataTypes=WND&bbox=-33.839,151.175,-33.903,151.239&pageNum=4. 8 Appendix 8.1 Sunspots Tables Table 17: Sunspots Classical VS Quantum - DM Stats MSE Sunspots Classical VS Quantum DM Stats MSE Quantum (pdq) blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (10, 1, 1) 10 MSE 4.941708 0.0000007744 2411.531786 1563.023210 848.508576 (9, 1, 3) 10 MSE 4.831355 0.0000013561 2411.531786 1608.609316 802.922470 (8, 1, 3) 10 MSE 4.812065 0.0000014938 2411.531786 1613.363812 798.167975 (6, 1, 3) 10 MSE 4.680626 0.0000028600 2411.531786 1603.130658 808.401129 (7, 1, 3) 10 MSE 4.608440 0.0000040570 2411.531786 1629.559531 781.972255 (5, 1, 3) 10 MSE 4.404979 0.0000105794 2411.531786 1682.786868 728.744918 (4, 1, 1) 10 MSE 2.754791 0.0058729550 2411.531786 2056.829346 354.702441 (3, 1, 1) 10 MSE 1.917276 0.0552029000 2411.531786 2179.575810 231.955977 Table 18: Sunspots Classical VS Quantum - DM Stats MAE Sunspots Classical VS Quantum DM Stats MAE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (10, 1, 1) 10 MAE 4.908669 0.0000009170 35.734117 28.387620 7.346497 (8, 1, 3) 10 MAE 4.609235 0.0000040415 35.734117 29.106819 6.627298 (9, 1, 3) 10 MAE 4.602818 0.0000041681 35.734117 29.066595 6.667522 (6, 1, 3) 10 MAE 4.539380 0.0000056420 35.734117 28.866935 6.867182 (7, 1, 3) 10 MAE 4.324537 0.0000152853 35.734117 29.306929 6.427188 (5, 1, 3) 10 MAE 4.242686 0.0000220860 35.734117 29.837574 5.896543 (4, 1, 1) 10 MAE 1.155393 0.2479295000 35.734117 34.423927 1.310189 (3, 1, 1) 10 MAE -0.005323 0.9957526000 35.734117 35.739900 -0.005784 8.2 Woolyarn Tables Table 19: Woolyarn Classical VS Quantum - DM Stats MSE Woolyarn Classical VS Quantum DM Stats MSE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (2, 1, 1) 1 MSE -7.968394 0.0000000000 934889.823 1136324 -201434.0761 (1, 1, 1) 1 MSE -7.952058 0.0000000000 934889.823 1185477 -250586.7362 (3, 1, 1) 1 MSE -7.245462 0.0000000000 934889.823 1131699 -196809.3958 (4, 1, 1) 1 MSE -6.427643 0.0000000001 934889.823 1042081 -107190.8507 (7, 1, 1) 1 MSE -4.500098 0.0000067922 934889.823 1343679 -408789.1725 (10, 1, 1) 1 MSE -4.071137 0.0000467843 934889.823 1026709 -91818.85786 (8, 1, 1) 1 MSE -4.023652 0.0000573025 934889.823 1024751 -89861.61058 (9, 1, 1) 1 MSE -3.937144 0.0000824572 934889.823 1009798 -74908.25687 (6, 1, 1) 1 MSE -2.992495 0.0027670760 934889.823 942247.1 -7357.235698 (5, 1, 1) 1 MSE -2.548627 0.0108148000 934889.823 948237.4 -13347.5792 Table 20: Woolyarn Classical VS Quantum - DM Stats MAE Woolyarn Classical VS Quantum DM Stats MAE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (3, 1, 1) 1 MAE -7.36933 0.0000000000 826.718879 904.906818 -78.187938 (4, 1, 1) 1 MAE -7.148529 0.0000000000 826.718879 865.442915 -38.724036 (2, 1, 1) 1 MAE -6.434141 0.0000000001 826.718879 902.821988 -76.103109 (1, 1, 1) 1 MAE -5.220583 0.0000001784 826.718879 926.307162 -99.588282 (7, 1, 1) 1 MAE -4.272585 0.0000193220 826.718879 994.530279 -167.811399 (10, 1, 1) 1 MAE -3.251013 0.0011499450 826.718879 869.822437 -43.103558 (8, 1, 1) 1 MAE -3.194489 0.0014007860 826.718879 869.053863 -42.334983 (9, 1, 1) 1 MAE -3.039556 0.0023693500 826.718879 862.313859 -35.594980 (6, 1, 1) 1 MAE -1.778833 0.0752672000 826.718879 829.852377 -3.133498 (5, 1, 1) 1 MAE 0.767760 0.4426297000 826.718879 823.707948 3.010931 8.3 Muna Lua CO2 Tables Table 21: CO2 Classical vs Quantum - DM Stats (MSE) CO2 Classical vs Quantum - DM Stats (MSE) Quantum (pdq) blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (8, 1, 1) 174 MSE 8.681193 0.00E+00 9.080382 3.599858 5.480524 (9, 1, 1) 172 MSE 8.991140 0.00E+00 9.099370 3.282792 5.816578 (10, 1, 1) 174 MSE 9.194297 0.00E+00 9.080382 3.192490 5.887893 (5, 1, 1) 174 MSE -6.729980 1.70E-11 9.080382 9.236257 -0.155875 (6, 1, 1) 174 MSE -6.673577 2.50E-11 9.080382 9.247513 -0.167131 (4, 1, 1) 174 MSE -5.947968 2.71E-09 9.080382 9.217262 -0.136880 (3, 1, 1) 174 MSE -5.288334 1.23E-07 9.080382 9.234463 -0.154080 (7, 1, 1) 174 MSE 5.230577 1.69E-07 9.080382 6.563634 2.516748 Table 22: CO2 Classical VS Quantum - DM Stats MAE CO2 Classical VS Quantum DM Stats MAE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (8, 1, 1) 174 MAE 9.252065 0.00E+00 2.315727 1.406047 0.909681 (9, 1, 1) 172 MAE 9.561353 0.00E+00 2.314338 1.348458 0.965881 (10, 1, 1) 174 MAE 9.857719 0.00E+00 2.315727 1.324859 0.990868 (3, 1, 1) 174 MAE -5.639472 1.71E-08 2.315727 2.338754 -0.023027 (7, 1, 1) 174 MAE 5.486284 4.10E-08 2.315727 1.887741 0.427986 (4, 1, 1) 174 MAE -5.417399 6.05E-08 2.315727 2.332440 -0.016713 (5, 1, 1) 174 MAE -4.536355 5.72E-06 2.315727 2.328392 -0.012664 (6, 1, 1) 174 MAE -4.474334 7.66E-06 2.315727 2.329375 -0.013647 8.4 Australian Beer Production Tables Table 23: AusBeer Classical VS Quantum - DM Stats (MSE) AusBeer Classical VS Quantum DM Stats MSE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (1, 1, 6) 84 MSE 11.878028 0.00E+00 2559.328583 1296.194617 1263.133966 (2, 1, 6) 84 MSE 15.353045 0.00E+00 2559.328583 564.606641 1994.721942 (7, 1, 1) 84 MSE 15.626542 0.00E+00 2559.328583 388.120011 2171.208572 (6, 1, 1) 84 MSE 15.626973 0.00E+00 2559.328583 412.991120 2146.337463 (5, 1, 1) 84 MSE 15.708202 0.00E+00 2559.328583 408.258903 2151.069680 (10, 1, 1) 84 MSE 15.730669 0.00E+00 2559.328583 381.998708 2177.329875 (9, 1, 1) 84 MSE 15.775209 0.00E+00 2559.328583 378.526593 2180.801990 (3, 1, 3) 84 MSE 16.062654 0.00E+00 2559.328583 410.665446 2148.663137 Table 24: AusBeer Classical VS Quantum - DM Stats (MAE) AusBeer Classical VS Quantum DM Stats MAE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (1, 1, 6) 84 MAE 11.772301 0.00E+00 41.666725 29.306057 12.360668 (2, 1, 6) 84 MAE 18.632163 0.00E+00 41.666725 18.651992 23.014733 (3, 1, 3) 84 MAE 19.375578 0.00E+00 41.666725 16.011428 25.655297 (6, 1, 1) 84 MAE 19.859260 0.00E+00 41.666725 16.104594 25.562131 (7, 1, 1) 84 MAE 20.302789 0.00E+00 41.666725 15.588372 26.078353 (5, 1, 1) 84 MAE 20.327477 0.00E+00 41.666725 16.109335 25.557390 (9, 1, 1) 84 MAE 20.668863 0.00E+00 41.666725 15.037503 26.629222 (10, 1, 1) 84 MAE 20.854476 0.00E+00 41.666725 15.106466 26.560259 8.5 Sydney weather 2024 tables Table 25: Sydney 2024 Summer Temp - Classical vs. Quantum DM Stats (MSE) Sydney 2024 Summer Temp Classical VS Quantum DM Stats MSE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (4, 1, 1) 160 MSE -3.41981 6.27E-04 8.1677 9.023377 -0.855678 (3, 1, 1) 161 MSE -2.943498 3.25E-03 8.177383 19.809067 -11.631684 (5, 1, 1) 161 MSE -2.456801 1.40E-02 8.177383 8.713587 -0.536204 (6, 1, 1) 161 MSE -1.544826 1.22E-01 8.177383 8.668534 -0.491151 (10, 1, 1) 155 MSE 0.610649 5.41E-01 8.230546 7.906464 0.324082 (7, 1, 1) 161 MSE -0.577497 5.64E-01 8.177383 8.416235 -0.238851 (9, 1, 1) 159 MSE 0.555303 5.79E-01 8.258326 7.989777 0.268549 (8, 1, 1) 161 MSE 0.175505 8.61E-01 8.177383 8.099677 0.077707 Table 26: Sydney 2024 Summer Classical VS Quantum DM Stats (MAE) Sydney 2024 Summer Classical VS Quantum DM Stats MAE Quantum pdq blocks_used loss DM_stat p_value classical_mean_loss quantum_mean_loss delta_mean_loss (3, 1, 1) 161 MAE -4.038151 5.40E-05 2.183719 2.875055 -0.691336 (4, 1, 1) 160 MAE -2.975823 2.92E-03 2.180798 2.280019 -0.099221 (5, 1, 1) 161 MAE -1.775833 7.58E-02 2.183719 2.245571 -0.061852 (10, 1, 1) 155 MAE 1.205682 2.28E-01 2.193426 2.106318 0.087108 (6, 1, 1) 161 MAE -1.158140 2.47E-01 2.183719 2.237213 -0.053494 (9, 1, 1) 159 MAE 1.088950 2.76E-01 2.198211 2.124441 0.073770 (8, 1, 1) 161 MAE 0.601969 5.47E-01 2.183719 2.145549 0.038169 (7, 1, 1) 161 MAE -0.164042 8.70E-01 2.183719 2.193769 -0.010050