An inâdepth understanding of the mechanism concerning the parasitic disease called visceral leishmaniasis (VL) remains challenging. Thus, we modeled the dynamics of this illness using two fractionalâorder operators, including CaputoâFabrizio and AtanganaâBaleanu. In the proposed dynamical model, the endemic and diseaseâfree equilibrium points were considered the symmetrical components. The fractional Euler method was applied to simulate the developed model, thus determining the equilibrium points' stability. The numerical simulation results were compared with the measured data to validate the model. The results obtained from the optimum fractional operator disclosed the minimum absolute and relative errors. The primary outcome of our study is the successful application of fractionalâorder operators, specifically the AtanganaâBaleanu operators with
= 0.98, in modeling the dynamics of VL. Notably, the numerical simulation results, validated against real data from Sudan, demonstrated that the AtanganaâBaleanu operators with
= 0.98 yielded the best performance, with minimum absolute and relative errors. This underscores the precision of our fractional calculusâbased dynamical model in predicting VL dynamics compared to the classical framework, particularly for fractionalâorder values of
= 0.99 and
= 0.98.