We study the general NMSSM with an emphasis on the parameter regions with a very strong first-order electroweak phase transition (EWPT). In the presence of heavy fields coupled to the Higgs sector, the analysis can be problematic due to the existence of sizable radiative corrections. In this paper we propose a subtraction scheme that helps to circumvent this problem. For simplicity we focus on a parameter region that is by construction hidden from the current collider searches. The analysis proves that (at least) in the identified parameter region the EWPT can be very strong and striking gravitational wave signals can be produced. The corresponding gravitational stochastic background can potentially be detected at the planned space-based gravitational wave observatory eLISA, depending on the specific experiment design that will be approved.

By means of a relativistic microscopic approach we calculate the expansion velocity of bubbles generated during a first-order electroweak phase transition. In particular, we use the gradient expansion of the Kadanoff-Baym equations to set up the fluid system. This turns out to be equivalent to the one found in the semi-classical approach in the non-relativistic limit. Finally, by including hydrodynamic deflagration effects and solving the Higgs equations of motion in the fluid, we determine velocity and thickness of the bubble walls. Our findings are compared with phenomenological models of wall velocities. As illustrative examples, we apply these results to three theories providing first-order phase transitions with a particle content in the thermal plasma that resembles the Standard Model.