Recently the ATLAS collaboration has reported the first results of searches for heavy scalar resonances decaying into a Z boson and a lighter new scalar resonance, where the Z boson decays leptonically and the lighter scalar decays into a top-quark pair, giving rise to ℓ+ℓ−tt¯ final states. This had previously been identified as a smoking-gun signature at the LHC for a first-order electroweak phase transition (FOEWPT) within the framework of two Higgs doublet models (2HDMs). In addition, ATLAS also presented new limits where the Z boson decays into pairs of neutrinos and the lighter scalar resonance into bottom-quark pairs, giving rise to the ννbb¯ final state. We analyze the impact of these new searches on the 2HDM parameter space, with emphasis on their capability to probe currently allowed 2HDM regions featuring a strong FOEWPT. We also study the complementarity of these new searches with other LHC probes that could target the FOEWPT region of the 2HDM. Remarkably, the ATLAS search in the ℓ+ℓ−tt¯ final state shows a local 2.85 σ excess (for masses of about 650 GeV and 450 GeV for the heavy and light resonance) in the 2HDM parameter region that would yield a FOEWPT in the early universe, which could constitute the first experimental hint of baryogenesis at the electroweak scale. We analyze the implications of this excess, and discuss the detectability prospects for the associated gravitational wave signal from the FOEWPT. Furthermore, we project the sensitivity reach of the ℓ+ℓ−tt¯ signature for the upcoming runs of the LHC. Finally, we introduce the python package thdmTools, a state-of-art tool for the exploration of the 2HDM. © 2024, The Author(s).

The work of T.B. is supported by the German Bundesministerium für Bildung und Forschung (BMBF, Federal Ministry of Education and Research) — project 05H21VKCCA. The work of M.O.O.R is supported by the European Union’s Horizon 2020 research and innovation programme under grant agreement No 101002846, ERC CoG “CosmoChart”. J.M.N. was supported by the Ramón y Cajal Fellowship contract RYC-2017-22986, and acknowledges partial financial support by the grant PID2021-124704NB-I00, and by the European Union’s Horizon 2020 research and innovation programme under the Marie Skłodowska-Curie grant agreements No 860881-HIDDeN and 101086085-ASYMMETRY. S.H. and J.M.N. acknowledge partial financial support by the Spanish Research Agency (Agencia Estatal de Investigación) through the grant IFT Centro de Excelencia Severo Ochoa No CEX2020-001007-S funded by MCIN/AEI/10.13039/501100011033. The work of S.H. was supported in part by the grant PID2019-110058GB-C21 funded by MCIN/AEI/10.13039/501100011033 and by “ERDF A way of making Europe”, and in part by the Grant PID2022-142545NB-C21 funded by MCIN/AEI/10.13039/501100011033/FEDER, UE. K.R. and G.W acknowledge support by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC 2121 “Quantum Universe” – 390833306. This work has been partially funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) - 491245950.

Peer reviewed