We present a fractional order model for the spread of malware in wireless sensor networks that builds memory directly into the dynamics through the proportional Hadamard–Caputo operator. The network population is organized into six groups, namely susceptible, exposed, infectious, quarantined, recovered, and vaccinated devices. We recast the system as an integral equation using a logarithmic change of time and we prove two fixed point results, where existence follows from the nonlinear Leray–Schauder alternative and uniqueness is obtained by Banach’s contraction principle. We then establish stability in the sense of Ulam–Hyers and in its extended form, showing that small modeling or data errors lead to proportionally small changes in the solutions. For computation, we build a predictor and corrector scheme in the modified Adams Bashforth Moulton framework adapted to the proportional Hadamard Caputo kernel with exponential memory in logarithmic time. Simulations show that stronger memory or a lower fractional order slows decay and extends spread while values near the classical case bring rapid stabilization.
Barakat et al. (Sat,) studied this question.