  We study the effect of dark matter scalar and vector-mediated interactions on
dark matter admixed neutron stars. In particular, we exploit the two-fluid
formalism of Tolman-Oppenheimer-Volkoff equations as generic framework for dark
matter admixed neutron stars. The fluids couple to each other only by gravity.
In particular, the baryonic sector is described by the BSk22 equation of state,
whereas we employ a relativistic mean field model for the dark matter, where we
include the interaction of dark massive fermions with light dark scalar and
vector mediators. We consider both the linear and the quadratic scalar
interactions with the dark fermion. In the quadratic scalar scenario, we take
into account a quartic self-interaction that significantly affects the stellar
properties. Interestingly, in both cases, the effect of the scalar coupling is
smaller than that of the vector, though quantitative differences arise. We also
compute the sound speed of DM, finding that the scalar and quadratic couplings
have an important influence on it. We compare our results with GW1708017,
GW190425 and NICER data and constrain DM couplings and mass.
