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.
@article{arxiv.2509.01677,
title = {Machine Learning in the 2HDM2S model for Dark Matter},
author = {Rafael Boto and Tiago P. Rebelo and Jorge C. Romão and João P. Silva},
journal= {arXiv preprint arXiv:2509.01677},
year = {2026}
}