Authors: Firuz Kamalov, Hana Sulieman, Ahmed El Sayed, Sivaraj R., Nitesh Sughnani
Organizations: Department of Computational Sciences, Canadian University Dubai, Dubai, UAE · Department of Mathematics and Statistics, American University of Sharjah, Sharjah, UAE · Department of Mathematics, Dr. B.R. Ambedkar National Institute of Technology, Jalandhar, India · Knowledge and Human Development Authority, Dubai, UAE
Generative artificial intelligence promises to democratize personalized tutoring, yet the longitudinal stability of adaptive systems remains under-theorized. Current ethical discourse focuses on static bias. This paper argues that a distinct threat lies in algorithmic drift: a feedback loop in which minor, demographically asymmetric observation errors compound over time into systemic inequities. We make three contributions. First, we formalize algorithmic drift as a coupled learner--tutor dynamical system. Second, we corroborate this theory with an agent-based Monte Carlo simulation in which either a transparent online tracker or an LSTM deep-knowledge-tracing network drives an IRT-based curriculum policy. Third, we introduce the temporal fairness stabilizer (TFS), a targeted intervention that periodically re-probes learner capability, and benchmark it against standard exploration baselines with explicit group-fairness metrics and replication-level inference. Within the modeled environment, TFS closes the achievement gap where generic exploration does not, and we characterize the conditions under which it succeeds. We separate what is established analytically, what is demonstrated in simulation, and what remains to be validated with real learners. Our results reframe algorithmic fairness in education as a temporal property of a closed loop rather than a deployment-time checkbox.
Figures & tables
Symbol
Meaning
Default
Status
N , T , K
learners, time steps, knowledge components
10,000 , 100 , 50
design choice
∣Sprot∣/N
protected-group share
0.20
assumed (stylized)
θi(0)
initial latent ability
N(0.5,0.1) , clipped
assumed; varied in Sec. 4.7
αi
learning rate
N(0.02,0.005) , >0
calibrated (reference ceiling ≈0.92 )
λ
forgetting rate
0.002
assumed
ω
ZPD kernel width
0.12
assumed; sets ϵ⋆
Table 1: Simulation parameters. “Calibrated” values were set so that the noiseless reference group reaches the documented mastery ceiling within the horizon; “assumed” values are modeling choices whose influence is examined in the sensitivity and robustness analyses.
Fig 1: Dynamical system representation of algorithmic drift in the adaptive learning environment. The recursive loop, illustrates how asymmetric observation noise ( ϵ ) generated by linguistic friction depresses DKT mastery estimations ( p^t ). This leads the generative policy to issue a remedial curriculum ( βt+1 ), which in turn stalls the learner’s actual cognitive state ( θt+1 ). The TFS serves as a regulatory layer that interrupts this structural trajectory decoupling by forcing periodic exploration probes.
Fig 2: Flowchart of the TFS algorithm. Operating as a targeted regularization layer over the standard generative policy, TFS tracks a longitudinal pessimism counter ( κ ) to detect sustained periods of remedial curriculum assignment ( β<μthreshold ). Once the consecutive threshold ( τ ) is reached, it overrides the depressed estimations of the DKT, forcing a diagnostic exploration probe to verify actual cognitive ceiling and break the algorithmic drift loop. For visual clarity the figure depicts the absolute-threshold form of the pessimism trigger; the deployed policy uses the cohort-relative criterion and capability ratchet of Algorithm 1 .
Grader
FNref
FNprot
ϵ^
McNemar p
Cohen’s h
Pˉ(C)ref
Pˉ(C)prot
ΔP ( p )
items ↓
FPref / FPprot
Qwen2.5-7B-Instruct
0.000
0.050
+0.050
0.5
+0.45
0.993
0.939
+0.055 (0.043)
95%
0.00 / 0.00
Phi-4-mini-instruct
0.000
0.000
+0.000
1
+0.00
0.995
0.993
+0.001 (0.64)
98%
0.10 / 0.05
Mistral-7B-Instruct-v0.3
0.000
0.000
+0.000
1
+0.00
1.000
0.999
+0.001 (0.31)
62%
0.30 / 0.25
Pooled (3 graders)
+0.017 (mean)
0.5
+0.019 (mean)
Table 2: Matched-pair LLM grader measurement (40 correct and 20 incorrect items, three instruction phrasings, modal verdicts). FN : false-negative rate on conceptually correct answers; ϵ^=FNprot−FNref ; McNemar: exact two-sided p on the paired verdicts; Pˉ(C) : mean grader probability of CORRECT; ΔP : paired difference Pˉ(C)ref−Pˉ(C)prot with its paired- tp -value; items ↓ : share of items on which P(C) is lower for the second-language variant; FP : false-positive rate on incorrect answers. The last row pools the discordant pairs over graders.
Fig 3: Algorithmic drift ( ϵ=0.05 , online tracker, mean trajectories over R=50 replications). Scenario A (left) illustrates the drift divergence, where the AI’s estimation of the protected group (dashed) decouples from their true latent ability (blue solid line), suppressing their learning trajectory compared to the reference group (green). Scenario B (right) shows that the TFS mitigates this estimation error, substantially reducing the achievement gap under modeled conditions.
Policy
Final ability ( Sprot )
Gap Δ
Drift D ( Sprot )
Learner-level d vs. control
Control (no exploration)
0.884 [0.883, 0.885]
0.038 [0.037, 0.039]
0.268
—
ϵ -greedy
0.888 [0.887, 0.889]
0.035 [0.034, 0.036]
0.270
+0.03
UCB
0.901 [0.900, 0.901]
0.018 [0.017, 0.019]
0.283
+0.13
Thompson sampling
0.836 [0.835, 0.837]
0.019 [0.019, 0.020]
0.263
−0.34
TFS (ours)
0.972 [0.971, 0.972]
0.003 [0.002, 0.004]
0.203
+0.83
Table 3: Comparison of TFS against standard exploration baselines under ϵ=0.05 . Values are means over R=50 seeds with replication-level 95% percentile-bootstrap CIs. The last column is Cohen’s d on protected-group final ability relative to the control policy, computed at the learner level and pooled over replications (positive favors the policy); the paired replication-level contrast for TFS is reported in Section 4.10 .
Policy
Gβ
GZ
Gb
FNprot / FNref
GFN
ρ
Control, ϵ=0 (parity reference)
+0.000
+0.000
+0.000
0.105 / 0.105
−0.001
1.000
Control (no exploration)
+0.070
+0.059
−0.047
0.138 / 0.105
+0.033
0.911
ϵ -greedy
+0.064
+0.052
−0.049
0.145 / 0.115
+0.030
0.916
UCB
+0.059
+0.046
−0.056
0.098 / 0.066
+0.032
0.957
Thompson sampling
+0.065
+0.022
−0.058
0.127 / 0.105
+0.022
0.945
TFS (ours)
+0.033
+0.045
−0.026
0.037 / 0.032
+0.006
0.994
Table 4: Group-fairness metrics at ϵ=0.05 (means over R=50 replications; time-averaged over t=1,…,T ). Gβ : exposure gap (reference minus protected mean served difficulty); GZ : ZPD-access gap; Gb : estimation-bias gap (protected minus reference signed calibration error; negative = protected group underestimated relative to reference); FNprot/FNref : remedial-misclassification rates with their gap GFN ; ρ : gain ratio. The first row is the control policy with a noiseless channel ( ϵ=0 ), the parity reference point. Replication-level 95% CI half-widths are at most 0.002 for every entry.
Table 5: Main experiment with the LSTM-DKT in the loop ( N=10,000 , T=100 , R=20 replications per cell; replication-level 95% CIs in brackets). Held-out next-outcome AUC of the network on the heterogeneous pre-training population and on the deployment population, against the oracle AUC of the true success probability. Δ : gap widening; Gb : estimation-bias gap; Gβ : exposure gap. The last column gives the online tracker’s gap at the same (a,ϵ) ( R=10 ).
Fig 4: Online tracker versus LSTM-DKT in the loop. Left: estimation-bias gap Gb under the control policy as a function of ϵ (negative: protected group under-estimated relative to the reference group). Right: ability gap Δ under control (solid) and TFS (dashed). The LSTM reproduces the estimation bias at every ϵ and the trapping transition at larger ϵ ; TFS holds the gap at or below zero for both estimators.
Fig 5: Algorithmic drift persists at a lower noise parameter ϵ=0.02 (online tracker, mean trajectories over R=50 replications; gap widening Δ=0.013 under control versus 0.001 under TFS).
Control
TFS
Noise ϵ
Final ability ( Sprot )
Δ
D ( Sprot )
Gb
Final ability ( Sprot )
Δ
Gb
0.00
0.922
0.000
0.235
+0.000
0.975
0.000
+0.000
0.01
0.916
0.006
0.242
−0.009
0.975
0.000
−0.005
0.02
0.909
0.013
0.248
−0.019
0.975
0.000
−0.011
0.03
0.901
0.020
0.255
−0.028
0.974
0.001
−0.016
0.05
0.884
0.037
0.268
−0.047
0.972
0.002
−0.026
Table 6: Sensitivity of learning outcomes to observation noise ( R=25 replications per cell; reference-group final ability 0.922 under control and 0.975 under TFS throughout). Gap widening Δ carries a replication-level 95% bootstrap CI of at most ±0.001 at every row; Gb is the estimation-bias gap.
Initial-ability distribution of both groups, θi(0)∼N(μ0,σ0)
μ0
σ0
Control Δ
Control D
TFS Δ
0.3
0.05
+0.042
0.246
+0.004
0.3
0.10
+0.040
0.247
+0.002
0.3
0.20
+0.034
0.252
−0.006
0.5
0.05
+0.039
0.267
+0.004
0.5
0.10
+0.037
0.268
+0.002
Table 7: Initial-state and group-heterogeneity sensitivity at ϵ=0.05 ( R=25 ). Top: initial-ability distribution of both groups (control gap Δ , control drift D , and TFS gap). Middle: initial estimate offset ξ of the tracker. Bottom: protected group with a different prior ability or learning rate; Δ0 is the gap with a noiseless channel, Δ0.05 the gap with the biased channel, and the excess Δ0.05−Δ0 is the drift-attributable component; the last column is the TFS gap at ϵ=0.05 .
Cohort margin m
τ=3
τ=5
τ=8
τ=12
0.05
−0.006
+0.000
+0.013
+0.023
0.10
−0.005
+0.002
+0.015
+0.025
0.15
−0.001
+0.007
+0.018
+0.026
Table 8: TFS hyperparameter ablation at ϵ=0.05 : residual protected-group gap Δ as a function of patience τ and cohort margin m ( R=25 ). The recovered gap is near zero across the grid.
ϵ=0.05
ϕ \ ψ
0.0
0.05
0.1
0.00
+0.002 ( 0.972 )
+0.004 ( 0.967 )
+0.004 ( 0.964 )
0.25
+0.003 ( 0.969 )
+0.004 ( 0.965 )
+0.005 ( 0.962 )
0.50
+0.003 ( 0.965 )
+0.005 ( 0.961 )
+0.005 ( 0.958 )
1.00
+0.006 ( 0.952 )
+0.007 ( 0.949 )
+0.008 ( 0.947 )
ϵ=0.15
Table 9: TFS under probe side effects ( R=25 ): residual gap Δ (and protected-group final ability in parentheses) as a function of the frustration penalty ϕ and the entrenchment penalty ψ . The unmitigated control gap is 0.037 at ϵ=0.05 and 0.147 at ϵ=0.15 .
Algorithmic fairness studies how automated decisions distribute errors, opportunities, and harms across groups, and how to measure and reduce disparities.
In the presence of feedback loops, algorithmic decisions at one time step can influence observed data in subsequent steps, potentially introducing selection bias or shifts in data distribution. Research has demonstrated that ignoring these effects can lead to unfairness, which may worsen over time as the model remains deployed within society. This raises important questions of how to measure and ensure fairness of algorithmic decision-making models over time. In this short paper, we present a taxonomy that distinguishes between different types of static and temporal fairness notions in recent advancements. We then suggest prospects for future research by identifying common implicit temporal assumptions, open questions to enhance our temporal understanding of fairness, and other open challenges.
Generative AI tutors provide real-time, personalized learning support, but also create a new education inequity: students with access to premium AI services may receive clearer explanations, more personalized guidance, and better scaffolding than students limited to free or low-cost services. To address this challenge, we propose FairTutor, an equity-aware model-routing framework that achieves cost-effective AI tutoring via pedagogically motivated multi-agent orchestration. FairTutor combines query analysis, pedagogical planning, low-cost model generation, evaluator-guided critique and revision, and selective escalation to premium AI models. We introduce access-tier AI Education (AIED) Advantage Gap to measure the quality difference between premium-access and budget-constrained tutoring, and TutorAccessEval, a benchmark spanning math, reading, writing, science, and language learning. Empirical evaluations show that FairTutor achieves 97.1% of premium pedagogical quality (in floor-adjusted Likert scale) while reducing serving cost by 71.6%. Sensitivity analysis reveals a tunable cost--quality Pareto frontier, enabling FairTutor to be tailored to the needs of diverse student populations.
We establish a theoretical link between the recently proposed "drifting" generative dynamics and gradient flows induced by the Sinkhorn divergence. In a particle discretization, the drift field admits a cross-minus-self decomposition: an attractive term toward the target distribution and a repulsive/self-correction term toward the current model, both expressed via one-sided normalized Gibbs kernels. We show that Sinkhorn divergence yields an analogous cross-minus-self structure, but with each term defined by entropic optimal-transport couplings obtained through two-sided Sinkhorn scaling (i.e., enforcing both marginals). This provides a precise sense in which drifting acts as a surrogate for a Sinkhorn-divergence gradient flow, interpolating between one-sided normalization and full two-sided Sinkhorn scaling. Crucially, this connection resolves an identifiability gap in prior drifting formulations: leveraging the definiteness of the Sinkhorn divergence, we show that zero drift (equilibrium of the dynamics) implies that the model and target measures match. Experiments show that Sinkhorn drifting reduces sensitivity to kernel temperature and improves one-step generative quality, trading off additional training time for a more stable optimization, without altering the inference procedure used by drift methods. These theoretical gains translate to strong low-temperature improvements in practice: on FFHQ-ALAE at the lowest temperature setting we evaluate, Sinkhorn drifting reduces mean FID from 187.7 to 37.1 and mean latent EMD from 453.3 to 144.4, while on MNIST it preserves full class coverage across the temperature sweep. Project page: https://mint-vu.github.io/SinkhornDrifting/
Ping He, Om Khangaonkar, Hamed Pirsiavash +2
Department of Computer Science, Vanderbilt University · Department of Computer Science, University of California, Davis · Department of Computer Science, Purdue University