Testing Z boson rare decays Z → H1γ , A1γ with (g-2 )μ, MW , and B R (hSM → Zγ ) in the NMSSM

We study the rare decay process of Z boson into photon, accompanied by a CP-even or CP-odd scalar. We present the analytical delineation of the processes through the model-independent parametrizations of the new physics couplings and, finally, consider the Next-to-Minimal Supersymmetric Standard Model to mark out the parameter space where the branching fraction can have the maximum value. As a part of the necessary phenomenological and experimental cross-checks, we aim to fit the anomalous magnetic moment of the muon and W boson mass anomaly through the supersymmetric contributions. We also find that the decays Z -> H-1 gamma , A(1)gamma can serve as an excellent complementary test to BR ( h(SM) -> Z gamma ). In fact, to facilitate future searches, we unveil a few benchmark points that additionally satisfy the deviation of BR ( h(SM) -> Z gamma ) from the SM value based on the recent measurements of ATLAS and CMS. Future proposals such as ILC, CEPC, and FCC-ee are anticipated to operate for multiple years, focusing on center-of-mass energy near the Z pole. Consequently, these projects will be capable of conducting experiments at the Giga-Z (10 (9) of Z bosons) and Tera-Z (10(12) of Z bosons) phases, which may probe the aforesaid rare decay processes, thus the model as well. These unconventional yet complementary searches offer different routes to explore the supersymmetric models with extended Higgs sectors like NMSSM.