Perturbative unitarity for models with singlet and doublet scalars

We provide a complete description of perturbative unitarity bounds on the gauge-scalar sectors of models with extra $SU(2)$ doublet, neutral singlet, and charged singlet scalars. Such additions are very frequent in models beyond the Standard Model, and, in particular, they are almost universal in models explaining the dark matter problem. We propose a specific classification and minimal set of scattering matrices containing all the relevant information. We also developed a Mathematica implementation of our results, BounDS, and we use it to fully study a number of simple cases, comparing with the literature, when available. The Mathematica notebook BounDS is provided via a public GitHub repository.