In this paper, the ABC fractional derivative is used to provide a mathematical model for the dynamic systems of substance addiction. The basic reproduction number is investigated, as well as the equilibrium points' stability. Using fixed point theory and nonlinear analytic techniques, we verify the theoretical results of solution existence and uniqueness for the proposed model. A numerical technique for getting the approximate solution of the suggested model is established by using the Adams type predictor‐corrector rule for the ABC‐fractional integral operator. There are several numerical graphs that correspond to different fractional orders. Furthermore, we present a numerical simulation for the transmission of substance addiction in two scenarios with fundamental reproduction numbers greater than and fewer than one.