Machine Learning in the 2HDM2S model for Dark Matter

We introduce a two real scalar singlet extension of the two Higgs doublet model. We study the vacuum structure, the bounded from below conditions, the restrictions from the oblique parameters S,T and U, as well as the unitarity constraints. We submit the model to collider and Dark Matter experimental constraints and explore its allowed parameter space. We compare randomly populated simulations, simulations starting near the alignment limit, and a Machine Learning based exploration. Using Evolutionary Strategies, we efficiently search for regions with a viable Dark Matter candidate.