A model for the acoustic production of gravitational waves at a first-order phase transition is presented. The source of gravitational radiation is the sound waves generated by the explosive growth of bubbles of the stable phase. The model assumes that the sound waves are linear and that their power spectrum is determined by the characteristic form of the sound shell around the expanding bubble. The predicted power spectrum has two length scales, the average bubble separation and the sound shell width when the bubbles collide. The peak of the power spectrum is at wave numbers set by the sound shell width. For a higher wave number k, the power spectrum decreases to k −3 . At wave numbers below the inverse bubble separation, the power spectrum goes to k 5 . For bubble wall speeds near the speed of sound where these two length scales are distinguished, there is an intermediate k 1 power law. The detailed dependence of the power spectrum on the wall speed and the other parameters of the phase transition raises the possibility of their constraint or measurement at a future space-based gravitational wave observatory such as LISA. DOI: 10.1103/PhysRevLett.120.071301 Interest in gravitational waves from a first-order electroweak phase transition in the early Universe [1-3] has greatly increased following the European Space Agency's approval of a space-based gravitational wave observatory [4] and the detection of gravitational waves from a merging black hole binary [5]. At the same time, it has been realized that early work on gravitational waves from a thermal phase transition [6,7] greatly underestimated their energy density [8]. The first three-dimensional hydrodynamic simulations [8,9] revealed that the dominant source of gravitational waves was acoustic production from sound waves generated by the explosive growth of bubbles of the stable phase. In fact, it had been pointed out long before that sound waves were a source of gravitational waves [10], but subsequent work had not appreciated that the sound wave source persisted for long after the phase transition completed, hence boosting the signal by orders of magnitude. The original model of gravitational radiation from the colliding bubble walls may still be relevant for nearvacuum transitions [11], and a semianalytic approach has recently been developed [12].The simulations in Refs. [8,9] revealed a power spectrum peaked at a wavelength around the average bubble separation R Ã , with a power law k −p at wave number k ≫ R −1 Ã . Where the power law was clear, the index was somewhere in the range −3 ≲ p ≲ −4. There was also evidence for some structure in the peak: where the bubble wall speed v w was closer to the speed of sound, the peak was broader. Understanding the gravitational wave power spectrum is of great importance for LISA's detection prospects [13], so it is vital to have a better physical understanding of the numerical simulations.In this Letter, I outline a model for the acoustic gravitational wave power spectrum based on the observation that the shel...