This research focuses on the theoretical asymptotic stability and long-time decay of the zero solution for a system of time-fractional nonlinear Schrödinger delay equations (NSDEs) in the context of the Caputo fractional derivative. Using the fractional Halanay inequality, we demonstrate theoretically when the considered system decays and behaves asymptotically, employing an energy function in the sense of the L2 norm. Together with utilizing the finite difference method for the spatial variables, we investigate the long-time stability for the semi-discrete system. Furthermore, we operate the L1 scheme to approximate the Caputo fractional derivative and analyze the long-time stability of the fully discrete system through the discrete energy of the system. Moreover, we demonstrate that the proposed numerical technique energetically captures the long-time behavior of the original system of NSDEs. Finally, we provide numerical examples to validate the theoretical results.