Are there BPS dyons in the generalized SU(2) Yang–Mills–Higgs model?

We use the well-known Bogomolny’s equations, in general coordinate system, for BPS monopoles and dyons in the SU (2) Yang–Mills–Higgs model to obtain an explicit form of BPS Lagrangian density under the BPS Lagrangian method. We then generalize this BPS Lagrangian density and use it to derive several possible generalized Bogomolny’s equations, with(out) additional constraint equations, for BPS monopoles and dyons in the generalized SU (2) Yang–Mills–Higgs model. We also compute the stress–energy–momentum tensor of the generalized model, and argue that the BPS monopole and dyon solutions are stable if all components of the stress-tensor density are zero in the BPS limit. This stability requirement implies the scalar fields-dependent couplings to be related to each other by an equation, which is different from the one obtained in Atmaja and Prasetyo (Adv High Energy Phys 2018:7376534, arXiv: 1803.06122 , 2018), and then picks particular generalized Bogomolny’s equations, with no additional constraint equation, out of those possible equations. We show that the computations in [1] are actually incomplete. Under the Julia–Zee ansatz, the generalized Bogomolny’s equations imply all scalar fields-dependent couplings must be constants, whose solutions are the BPS dyons of the SU (2) Yang–Mills–Higgs model (Prasad and Sommerfield in Phys Rev Lett 35:760, 1975), or in another words there are no generalized BPS dyon solutions under the Julia–Zee ansatz. We propose two possible ways for obtaining generalized BPS dyons, where at least one of the scalar fields-dependent couplings is not constant, that are by using different ansatze, such as axially symmetric ansatz for higher topological charge dyons; and/or by considering the most general BPS Lagrangian density.