---
title: 'QLIF-CAST: Quantum-Classical Recurrent Forecasting'
url: https://www.emergentmind.com/topics/qlif-cast
type: topic
---

# QLIF-CAST: Quantum-Classical Recurrent Forecasting

Searching arXiv for the exact QLIF-CAST paper and closely related context papers.
QLIF-CAST is a hybrid quantum–classical recurrent forecasting architecture that adapts the Quantum Leaky Integrate-and-Fire (QLIF) neuron from classification to continuous-valued time-series regression, with a specific emphasis on short-term multivariate weather prediction. In this formulation, neuron excitation states are encoded as single-qubit quantum superpositions updated by \(R_x\) rotation gates and a \(T_1\) relaxation leak, and the resulting excitation probabilities are embedded within a seven-layer recurrent network for next-step forecasting [2605.18333]. The model was evaluated both against a parameter-matched classical LIF baseline on multivariate weather data and against state-of-the-art quantum LSTM and quantum neural network baselines on air-quality and wind-speed benchmarks, where it achieved lower error than the matched classical LIF model and substantially shorter training times than deeper variational quantum alternatives [2605.18333].

## 1. Definition and conceptual role

QLIF-CAST is defined as a hybrid quantum–classical recurrent forecaster in which the conventional leaky-integrate-and-fire neuronal update is replaced by a single-qubit quantum LIF element. The central design choice is to represent excitation probability through a qubit state whose dynamics are driven by shallow quantum evolution rather than by a purely classical membrane-potential recursion [2605.18333].

The motivating task is short-term multivariate time-series forecasting of environmental variables, including temperature, humidity, wind speed, and pressure, as well as single-target regression tasks in air quality and wind speed [2605.18333]. The paper positions the model between two existing regimes: classical LSTM-based architectures, which offer strong predictive accuracy at higher computational cost, and classical LIF-based spiking models, which are computationally efficient but have more limited state dynamics. QLIF-CAST is presented as an attempt to retain shallow, efficient neuron-level computation while introducing richer temporal state evolution through superposition, interference, and probabilistic measurement [2605.18333].

A plausible implication is that QLIF-CAST is best understood not as a deep variational quantum model, but as a shallow quantum-neuronal replacement inside an otherwise standard recurrent regression stack. This interpretation is consistent with the fact that the architecture contains no trainable quantum angles and relies on classical backpropagation through surrounding layers rather than parameter-shift-based quantum optimization [2605.18333].

## 2. Quantum neuron formulation

The quantum component is a single-qubit representation of excitation. The qubit state is written as
\[
|\psi\rangle = \alpha_0 |0\rangle + \beta_1 |1\rangle,
\]
with excitation probability
\[
p_1 = |\beta_1|^2.
\]
In QLIF-CAST, the excitation probability \(\alpha\) is represented through a rotation angle \(\phi\) such that
\[
\alpha = \sin^2(\phi/2), \qquad \phi = 2 \arcsin(\sqrt{\alpha}).
\]
The per-timestep circuit is
\[
|0\rangle \rightarrow R_x(\phi_t) \rightarrow R_x(\theta_{\text{input},t}) \rightarrow \text{Measure},
\]
yielding
\[
\alpha_{t+1} = \sin^2\!\left(\frac{\phi_t + \theta_{\text{input},t}}{2}\right).
\]
A spike is emitted if \(\alpha_{t+1} \ge 0.75\), after which the state is reset [2605.18333].

The gate dynamics use the standard Pauli-\(X\) rotation
\[
R_x(\theta) = \exp(-i\theta X/2) = \cos(\theta/2) I - i \sin(\theta/2) X.
\]
This is a depth-2, single-qubit construction, and each QLIF neuron executes two \(R_x\) gates per timestep [2605.18333].

Leak is implemented through a \(T_1\)-style amplitude damping mechanism. The Kraus operators are
\[
E_0 =
\begin{bmatrix}
1 & 0\\
0 & \sqrt{1-\gamma}
\end{bmatrix},
\qquad
E_1 =
\begin{bmatrix}
0 & \sqrt{\gamma}\\
0 & 0
\end{bmatrix},
\]
with
\[
\gamma = 1 - e^{-\Delta t / T_1}.
\]
In the operational model, leak acts in probability space as \(\alpha \rightarrow \alpha e^{-\tau/T_1}\), and the corresponding decay angle is
\[
\phi_{\text{decay},t} = 2\arcsin\!\left(\sqrt{\alpha_t e^{-\tau/T_1}}\right).
\]
The effective input rotation is then
\[
\theta_{\text{input},t} = X_t \theta_t + (1-X_t)(-\phi_t + \phi_{\text{decay},t}),
\]
so that a previous spike causes the input-driven angle \(\theta_t\) to dominate, while otherwise the evolution is governed by decay [2605.18333].

The recurrent update for neuron \(j\) is therefore specified by
\[
\phi_{j,t} = 2\arcsin(\sqrt{\alpha_{j,t}}),
\]
\[
\theta_{j,t} = w_j^\top f_t + b_j,
\]
\[
\theta_{\text{input},j,t} = X_{j,t}\theta_{j,t} + (1-X_{j,t})\gamma_{j,t},
\]
\[
\alpha_{j,t+1} = \sin^2\!\left(\frac{\phi_{j,t} + \theta_{\text{input},j,t}}{2}\right).
\]
If \(\alpha_{j,t+1} \ge 0.75\), then \(s_{j,t+1}=1\) and \(\alpha_{j,t+1}\leftarrow 0\); otherwise \(s_{j,t+1}=0\) [2605.18333].

For comparison, the classical reference neuron follows the discrete LIF update
\[
U_{t+1} = \beta U_t + (1-\beta) I_t,
\]
with \(\beta = e^{-\Delta t/\tau}\), and emits a spike when \(U_{t+1}\) crosses threshold [2605.18333].

## 3. Network architecture and learning pipeline

QLIF-CAST is embedded in a seven-layer recurrent network whose total parameter count is approximately \(11{,}140\), and this total is identical for the quantum and classical baselines [2605.18333]. The layer stack is fixed except for the neuronal update in the second layer.

| Layer | Function | Parameters |
|---|---|---:|
| L1 | TimeDistributed Dense(48, ReLU) + Dropout(0.1) | 240 |
| L2 | 48-neuron QLIF layer or classical LIF layer | 2,400 |
| L3 | BatchNormalization + Dropout(0.2) | 96 |
| L4 | LSTM(24) | 7,008 |
| L5 | Dense(32, ReLU) + Dropout(0.2) | 800 |
| L6 | Dense(16, ReLU) | 528 |
| L7 | Linear Dense(\(N\)) output | \(\le 68\) |

Each input window has length \(T=12\). For weather and wind forecasting, this corresponds to 12 hours; for air quality, 12 days [2605.18333]. Each window is standardized to zero mean and unit variance. At each timestep \(t\), the input \(x_t\) is first projected through
\[
f_t = \mathrm{ReLU}(W_1 x_t + b_1),
\]
after which the QLIF layer computes 48 excitation trajectories in parallel, producing a sequence of shape \((T,48)\). This sequence is normalized, passed to an LSTM with hidden size 24, and then processed by two dense regression heads to produce the next-step forecast \(\hat{z}(t+1)\) [2605.18333].

Training uses mean squared error,
\[
\mathrm{MSE}(\hat{z},y)=\frac{1}{N}\sum_i (\hat{z}_i-y_i)^2,
\]
with mean absolute error, RMSE, and \(R^2\) reported as evaluation metrics on inverse-standardized predictions [2605.18333]. Optimization uses Adam with initial learning rate \(1\times 10^{-3}\), exponential decay, batch size 64, and early stopping with patience 5. L2 regularization is set to \(1\times 10^{-4}\), and dropout ranges from 0.1 to 0.2 [2605.18333].

Because the quantum angles are not trainable, the model does not use parameter-shift rules. Instead, backpropagation flows through classical parameters, while the spike threshold is differentiated through a surrogate gradient based on
\[
\hat{S}(\alpha) = \frac{1}{\pi}\arctan(\pi \alpha)+0.5,
\qquad
\hat{S}'(\alpha)=\frac{1}{1+\pi^2\alpha^2}.
\]
This design is one of the main reasons the architecture trains substantially faster than variational quantum sequence models [2605.18333].

## 4. Datasets, implementation, and hardware execution

The primary controlled experiment uses the D1 Weather History dataset from Kaggle, containing 96,453 hourly samples from April 2006 to September 2016. Inputs and targets are Temperature (\(^\circ\)C), Relative Humidity (fraction), Wind Speed (km/h), and Pressure (mb). The forecasting target is the next hour’s 4-variable vector, and the chronological split uses 10,000 training samples and 2,000 test samples [2605.18333].

Two additional datasets are used for cross-domain comparison. D2 is the Bangkok Air Quality dataset from WAQI, containing approximately 3,290 daily records from July 2016 to January 2026, with inputs PM10, O\(_3\), and NO\(_2\), and target next-day PM2.5. D3 is an Open-Meteo wind-speed dataset with 35,064 hourly samples from January 2020 to December 2023 at 100 m altitude, with next-hour wind speed as the target [2605.18333].

Simulation is performed with PennyLane on the `lightning.qubit` state-vector backend, and the surrounding network is implemented in TensorFlow/Keras [2605.18333]. The reported training platform is an Intel Core i9-13900H CPU with 16 GB RAM. Crucially, the per-timestep quantum workload is shallow: for 48 neurons, there are 96 \(R_x\) gates and 48 measurements per timestep, with no entanglement and no multi-qubit gates [2605.18333].

The implementation uses vectorized execution by flattening angles across neurons, timesteps, and batch elements into a single parallel PennyLane call, producing approximately \(500\times\) speedup relative to naïve sequential circuit execution [2605.18333]. This is a central implementation detail: the practical viability of the model depends less on raw qubit count than on the fact that the quantum subroutine is embarrassingly parallel and extremely shallow.

Hardware verification was performed on IBM Marrakesh, a 156-qubit QPU, by running the core single-qubit sequence
\[
|0\rangle \rightarrow R_x(\phi) \rightarrow R_x(\theta) \rightarrow \text{Measure}
\]
with 1,000 shots per test case [2605.18333]. The average deviation between QPU-measured \(P(|1\rangle)\) and ideal simulation was \(1.2\%\), with representative comparisons \(0.1516\) vs. \(0.1590\), \(0.7081\) vs. \(0.6850\), and \(0.9682\) vs. \(0.9620\) [2605.18333]. This supports the claim that the neuron-level circuit is hardware-compatible in the NISQ setting.

## 5. Empirical performance

In the parameter-matched comparison on the multivariate weather task, QLIF-CAST achieved MSE \(17{,}897\), MAE \(35.54\), RMSE \(133.8\), and \(R^2=0.896\), whereas the classical LIF baseline achieved MSE \(21{,}152\), MAE \(37.18\), RMSE \(145.4\), and \(R^2=0.877\) [2605.18333]. These correspond to a \(15.4\%\) reduction in MSE, a \(4.4\%\) reduction in MAE, and an \(8.0\%\) reduction in RMSE for the quantum model under identical parameter count and training schedule [2605.18333].

The per-variable MAE breakdown is more mixed. QLIF-CAST improves Temperature from \(2.93\) to \(2.65\) and Pressure from \(142.6\) to \(135.9\), but performs worse on Humidity (\(0.068\) vs. \(0.056\)) and Wind Speed (\(3.51\) vs. \(3.11\)) [2605.18333]. The paper interprets this as indicating that gains are concentrated on variables with stronger temporal structure, especially temperature and pressure [2605.18333].

Cross-domain comparisons place QLIF-CAST in a different region of the speed–accuracy trade-off from deeper quantum sequence models. On the Bangkok PM2.5 task, QLIF-CAST obtained MAE \(15.73\ \mu\mathrm{g}/\mathrm{m}^3\) and RMSE \(20.62\ \mu\mathrm{g}/\mathrm{m}^3\), compared with published QLSTM values of MAE \(11.24\) and RMSE \(15.06\). However, QLIF-CAST converged in 26 epochs versus 100 epochs for QLSTM [2605.18333]. On the wind-speed benchmark, QLIF-CAST achieved RMSE \(5.58\ \mathrm{km/h}\), MAE \(4.14\ \mathrm{km/h}\), and \(R^2=0.693\), while the published LSTM-QNN achieved RMSE \(3.92\), MAE \(2.87\), but required approximately \(65.3\) minutes of training, compared with \(3.88\) minutes for QLIF-CAST [2605.18333]. The reported training-time reduction is \(94\%\), or approximately \(16.8\times\) [2605.18333].

The paper therefore does not claim uniform dominance in absolute forecasting accuracy. Instead, it presents QLIF-CAST as occupying a distinct Pareto position in which shallow quantum neuronal dynamics improve over a matched classical LIF baseline, while training much faster than parameterized quantum recurrent baselines [2605.18333].

## 6. Interpretation, limitations, and research significance

The proposed explanation for QLIF-CAST’s behavior has three components. First, sequential \(R_x(\phi)\) and \(R_x(\theta)\) operations introduce constructive and destructive interference in amplitude space, which provides richer temporal modulation than scalar exponential decay [2605.18333]. Second, the amplitude-damping leak acts nonlinearly in probability space, so the effective decay differs qualitatively from the classical affine update \(U_{t+1}=\beta U_t + (1-\beta)I_t\) [2605.18333]. Third, the excitation probability \(\alpha \in [0,1]\) is interpreted as a probabilistic state variable that can naturally encode sensor uncertainty and may regularize optimization [2605.18333].

Several limitations are explicit. End-to-end QPU training was not demonstrated; hardware use was limited to circuit verification [2605.18333]. The cross-domain comparisons with QLSTM and LSTM-QNN rely on published numbers rather than retraining under identical conditions [2605.18333]. The architecture also remains shallow and single-qubit at the neuron level, so scaling to larger neuron counts or longer sequences may encounter hardware-throughput limits, and the effect of stronger decoherence remains an open question [2605.18333].

The paper suggests several future directions: multi-step forecasting horizons, ensembles of QLIF forecasters, hybrid QLIF–variational spiking-LSTM models, and multi-qubit neuron designs [2605.18333]. This suggests a broader research program in which quantum spiking elements are treated as modular replacements for classical neuronal dynamics rather than as end-to-end quantumized recurrent blocks.

Within the immediate literature supplied here, QLIF-CAST is distinct from work that studies QLIF as an information-flow diagnostic in many-body quantum chaos [2603.15315]. In that latter context, QLIF denotes Quantum Liang Information Flow, whereas in QLIF-CAST it denotes Quantum Leaky Integrate-and-Fire [2603.15315; 2605.18333]. The acronym overlap is purely nominal; the two uses refer to unrelated theoretical objects.

In summary, QLIF-CAST establishes that a shallow single-qubit spiking neuron can be integrated into a classical recurrent forecasting stack for regression, can outperform a parameter-matched classical LIF baseline on multivariate weather prediction, can train far faster than deeper variational quantum sequence models, and can be executed on present-day hardware with only \(1.2\%\) average deviation from simulation [2605.18333]. This suggests a design space in quantum machine learning where practical value may emerge less from deep trainable quantum circuits than from carefully engineered shallow quantum dynamical primitives embedded within classical sequence architectures.

Source: https://www.emergentmind.com/topics/qlif-cast