In this paper, we study a non-minimal gauged $$ \textrm{U}{(1)}_{L_{\mu }-{L}_{\tau }} $$
U
1
L
μ
−
L
τ
model, where we add two complex singlet scalars, three right-handed Majorana neutrinos (RHN), and a vector-like dark fermion to the Standard Model (SM), all non-trivially charged under the extra gauge symmetry. The model offers an easy resolution to the muon (g – 2) anomaly, which fixes the scale of spontaneous symmetry breaking. Furthermore, the two-zero minor structure in the RHN mass matrix provides successful predictions for neutrino oscillation parameters, including the Dirac phase. The extended scalar sector can easily induce first-order phase transitions. We identify all possible phase transition patterns in the three-dimensional field space. We quantify the associated gravitational waves from the sound wave source and demonstrate that the signatures can be observed in future space-based experiments. We find that strong first-order phase transitions require large values of scalar quartic couplings which constrain the scalar dark matter (DM) relic density to a maximum of 10−2 and 10−5 when we consider the DM direct detection bound. Nonetheless, the model successfully explains the DM relic density via contribution from the vector-like dark fermion. We show the allowed range of the model parameters that can address all the beyond SM issues targeted in this study.