Atmospheric neutrinos provide a unique avenue to study neutrino interactions in matter. In this work, the prospects of constraining non-standard neutrino interactions with atmospheric neutrino oscillations are investigated for the proposed ESSnuSB far detector. By analyzing atmospheric neutrino samples equivalent to 5.4 Mt·year exposure, it is found that ESSnuSB could be able to set the upper bounds \( \mid {\epsilon}_{e\mu}^m\mid <0.053 \) , \( \mid {\epsilon}_{e\tau}^m\mid <0.057 \) , \( \mid {\epsilon}_{\mu \tau}^m\mid <0.021 \) , \( {\epsilon}_{ee}^m-{\epsilon}_{\mu \mu}^m<0.075 \) and \( \mid {\epsilon}_{\tau \tau}^m-{\epsilon}_{\mu \mu}^m\mid <0.031 \) at 90% CL, when the results are minimized for \( {\phi}_{e\mu}^m \) , \( {\phi}_{e\tau}^m \) and \( {\phi}_{\mu \tau}^m \) and normal ordering is assumed for neutrino masses. It is also shown that the presence of non-standard interactions could affect the sensitivities to neutrino mass ordering and θ23 octant in comparison to the standard interaction scheme. The results of this work highlight the complementarity between atmospheric and accelerator neutrino programs in ESSnuSB.
