Implications For First-Order Cosmological Phase Transitions From The Third Ligo-Virgo Observing Run

We place constraints on the nornialized energy density in gravitational waves from first-order strong phase transitions using data from Advanced LIGO and Virgo's first, second, and third observing runs. First, adopting a broken power law model, we place 95% confidence level upper limits simultaneously on the gravitational-wave energy density at 25 Hz from unresolved compact binary mergers, Omega(CBC) < 6.1 x 10(-9), and strong first-order phase transitions, Omega(BPL) < 4.4 x 10(-9). The inclusion of the former is necessary since we expect this astrophysical signal to be the foreground of any detected spectrum. We then consider two more complex phenomenological models, limiting at 25 Hz the gravitational-wave background due to bubble collisions to a Omega(pt) < 5.0 x 10(-9) and the background due to sound waves to Omega(pt) < 5.8 x 10(-9) at 95% confidence level for phase transitions occurring at temperatures above 10(8) GeV.