Gravitational wave energy spectral density properties from BPASS Galactic binary population in the Milky Way galaxy

Abstract We analyse the energy spectral density properties of Gravitational waves from Galactic binary populations in the mHz band targeted by the Laser Interferometer Space Antenna mission. Our analysis is based on combining BPASS with a Milky Way analogue galaxy from the Feedback In Realistic Environment (FIRE) simulations and the GWs these populations emit. Our investigation compares different functional forms of gravitational wave (GW) ESDs, namely the single power-law, broken power-law, and single-peak models, revealing disparities within and among Galactic binary populations. We estimate the ESDs for six different Galactic binary populations and the ESD of the total Galactic binary population for LISA. Employing a single power-law model, we predict a total Galactic binary GW signal amplitude α = $2.0^{+0.2}_{-0.2} \times 10^{-8}$ and a slope β = $-2.64 ^{+0.03}_{-0.04}$ and the ESD $\rm h^2 \Omega _{GW}$ = $1.1 ^{+0.1}_{-0.1} \times 10^{-9}$ at 3 mHz. For the Galactic WDB binary GW signal $\alpha = 1^{+0.02}_{-0.02} \times 10^{-10}$, $\beta = -1.56 ^{+0.03}_{-0.03}$ and $\rm h^2 \Omega _{GW} = 18 ^{+1}_{-1} \times 10^{-12}$. Our analysis underscores the importance of accurate noise parameter estimation and highlights the complexities of modelling realistic observations, prompting future exploration into more flexible models.