The N2HDM is based on the CP-conserving 2HDM extended by a real scalar singlet field. Its enlarged parameter space and its fewer symmetry conditions as compared to supersymmetric models allow for an interesting phenomenology compatible with current experimental constraints, while adding to the 2HDM sector the possibility of Higgs-toHiggs decays with three different Higgs bosons. In this paper the N2HDM is subjected to detailed scrutiny. Regarding the theoretical constraints we implement tests of tree-level perturbativity and vacuum stability. Moreover, we present, for the first time, a thorough analysis of the global minimum of the N2HDM. The model and the theoretical constraints have been implemented in ScannerS, and we provide N2HDECAY, a code based on HDECAY, for the computation of the N2HDM branching ratios and total widths including the state-ofthe-art higher order QCD corrections and off-shell decays. We then perform an extensive parameter scan in the N2HDM parameter space, with all theoretical and experimental constraints applied, and analyse its allowed regions. We find that large singlet admixtures are still compatible with the Higgs data and investigate which observables will allow to restrict the singlet nature most effectively in the next runs of the LHC. Similarly to the 2HDM, the N2HDM exhibits a wrong-sign parameter regime, which will be constrained by future Higgs precision measurements.

Abstract:The 2-Higgs-Doublet Model (2HDM) belongs to the simplest extensions of the Standard Model (SM) Higgs sector that are in accordance with theoretical and experimental constraints. In order to be able to properly investigate the experimental Higgs data and, in the long term to distinguish between possible models beyond the SM, precise predictions for the Higgs boson observables have to be made available on the theory side. This requires the inclusion of the higher order corrections. In this work, we investigate in detail the renormalization of the 2HDM, a pre-requisite for the computation of higher order corrections. We pay particular attention to the renormalization of the mixing angles α and β, which diagonalize the Higgs mass matrices and which enter all Higgs observables. The implications of various renormalization schemes in next-to-leading order corrections to the sample processes H ± → W ± h/H and H → ZZ are investigated. Based on our findings, we will present a renormalization scheme that is at the same time process independent, gauge independent and numerically stable.

No abstract