In a vacuum first order phase transition gravitational waves are generated from collision of bubbles of the true vacuum. The spectrum from such collisions takes the form of a broken power law. We consider a toy model for such a phase transition where the dynamics of the scalar field depends on a single parameter $\bar{\lambda}$, which controls how thin the bubble wall is at nucleation and how close to degenerate the vacua are relative to the barrier. We extend on our previous work by performing a series of simulations with a range of $\bar{\lambda}$. The peak of the gravitational wave power spectrum varies by up to a factor of $1.3$, which is probably an unobservable effect. We find that the ultraviolet (UV) power law in the gravitational wave spectrum becomes steeper as $\bar{\lambda} \rightarrow 0$, varying between $k^{-1.4}$ and $k^{-2.3}$ for the $\bar{\lambda}$ considered. This provides some evidence that the form of the underlying effective potential of a vacuum first order phase transition could be determined from the gravitational wave spectrum it produces.

We have performed the first 3-dimensional simulations of strong first-order thermal phase transitions in the early Universe. For deflagrations, we find that the rotational component of the fluid velocity increases as the transition strength is increased. For detonations, however, the rotational velocity component remains constant and small. We also find that the efficiency with which kinetic energy is transferred to the fluid falls below theoretical expectations as we increase the transition strength. The probable origin of the kinetic energy deficit is the formation of reheated droplets of the metastable phase during the collision, slowing the bubble walls. The rate of increase in the gravitational wave energy density for deflagrations in strong transitions is suppressed compared to that predicted in earlier work. This is largely accounted for by the reduction in kinetic energy. Current modelling therefore substantially overestimates the gravitational wave signal for strong transitions with deflagrations, in the most extreme case by a factor of $10^{3}$. Detonations, however, are less affected.

We conduct large scale numerical simulations of gravitational wave production at a first order vacuum phase transition. We find a power law for the gravitational wave power spectrum at high wavenumber which falls off as $k^{-1.5}$ rather than the $k^{-1}$ produced by the envelope approximation. The peak of the power spectrum is shifted to slightly lower wave numbers from that of the envelope approximation. The envelope approximation reproduces our results for the peak power less well, agreeing only to within an order of magnitude. After the bubbles finish colliding the scalar field oscillates around the true vacuum. An additional feature is produced in the UV of the gravitational wave power spectrum, and this continues to grow linearly until the end of our simulation. The additional feature peaks at a length scale close to the bubble wall thickness and is shown to have a negligible contribution to the energy in gravitational waves, providing the scalar field mass is much smaller than the Planck mass.