In this work, we calculate the merger rate of primordial black hole–neutron star (PBH–NS) binaries within the framework of ellipsoidal-collapse dark matter models and compare it with that obtained from spherical-collapse dark matter halo models. Our results exhibit that ellipsoidal-collapse dark matter halo models can potentially amplify the merger rate of PBH–NS binaries in such a way that it is very close to the range estimated by the LIGO–Virgo observations. In contrast, spherical-collapse dark matter halo models cannot justify PBH–NS merger events as consistent results with the latest gravitational wave data reported by the LIGO–Virgo collaborations. In addition, we calculate the merger rate of PBH–NS binaries as a function of PBH mass and fraction within the context of ellipsoidal-collapse dark matter halo models. The results indicate that PBH–NS merger events with masses of M
PBH ≤ 5M
⊙,
M
NS ≃ 1.4M
⊙ will be consistent with the LIGO–Virgo observations if f
PBH ≃ 1.