A Stress-Adaptive Variable-Order Fractional Model for Motivational Dynamics with Memory Effects

Human motivation is governed by a long-memory cognitive process in which the depth of temporal integration—how far into the past the system draws upon accumulated experience—is not fixed, but dynamically compressed under cognitive stress. Despite extensive empirical evidence that acute stress impairs working memory and narrows temporal integration in decision-making, no existing mathematical framework has formally coupled the memory depth of the governing operator to a physiologically grounded stress indicator. To address this gap, we propose a stress-adaptive variable-order fractional model for motivational intensity M(t), in which the Caputo fractional order α(t) varies inversely with an aggregated stress indicator σ(t) through the Hill-type coupling α(t)=αmin+(αmax−αmin)C/(C+σ(t)), thereby encoding the empirically documented shift from deep integrative to shallow heuristic processing as cognitive load increases. Rather than deriving the model by algebraic manipulation of a differential equation, we formulate it directly as a causally consistent type-III Volterra integral equation, in which the memory kernel is evaluated at the history time s, ensuring that the weight assigned to each past state reflects the memory depth that was physiologically active when that state was experienced. Well-posedness is established rigorously via the Banach fixed-point theorem with explicit contraction constants, uniform boundedness and non-negativity of solutions are derived through the fractional Gronwall inequality, and numerical solutions are computed using an Adams–Bashforth–Moulton predictor–corrector scheme adapted to the variable-order kernel. Five numerical experiments demonstrate that stress-induced variation in α(t) produces qualitatively richer dynamics compared with the tested constant-order baselines: the proposed model achieves a steeper peak decline rate (0.48 versus 0.19–0.45), a larger burnout gap (3.15 versus 1.92–2.81), and faster recovery to ninety percent of peak motivation (4.2 versus 3.9–7.3 time units), while the empirically observed numerical convergence approaches O(h2) for sufficiently small step sizes. The framework offers a principled phenomenological substrate for memory-adaptive cognitive modelling, with direct implications for stress-aware intelligent tutoring systems that are capable of inferring α(t) in real time from biometric signals such as heart rate variability or galvanic skin response, and adjusting instructional complexity accordingly. Empirical calibration against learning-analytics and psychophysiological datasets, together with stochastic extensions for probabilistic burnout-risk prediction, are identified as immediate priorities for future research.