 | \small \emph{Left panel}: 1D and 2D posterior distributions derived from the PL template reconstruction of the benchmark PL-BNK\_1 ($\{\log_{10}(h^2\,\Omega_*), n_T\}=\{-12.5,\,2.085\}$) in the presence of galactic and extragalactic foregrounds. The posteriors of foregrounds and noise parameters are omitted for clarity. In the corner plot, the gold and red dots, and corresponding vertical lines, show the injected parameters and their reconstructed mean values. The 68\% C.L.~(95\% C.L.) reconstruction region is displayed in dark (light) gray. The equivalent regions and 1D posteriors obtained with the Fisher approximation are in blue. The top-right inset visualises the injected and reconstructed signals, with 68 and 95 \% CL error bands. The error bands on the galactic foreground and instrumental noise are too small to be visible. The LISA PLS is plotted in solid black. \emph{Right panel}: Like the left panel but for PL-BNK\_2 ($\{\log_{10}(h^2\,\Omega_*), n_T\}=\{-11,\,0.77\}$). In the inset plots of both panels, the error bands on the galactic foreground and instrumental noise are too small to be visible.\small 1D and 2D posterior distributions derived from the LN template reconstruction of the benchmarks LN-BNK\_1 (left panel) and LN-BNK\_2 (right panel). Lines' styles and color codes are as in \cref{fig:PL_MCMC}. LN-BNK\_1 and LN-BNK\_2 are defined as in \cref{eq:master_modelIII_bum} with $\{\log_{10}(h^2\,\Omega_*),\,\log_{10}\rho,\,f_*/{\rm mHz}\}$ set to $\{-9.9,\,-0.13,\, 11\}$ and $\{-11.8,\,-0.37,\,1\}$, respectively.\small 1D and 2D posterior distributions derived from the BPL template reconstruction of the benchmark BPL-BNK\_1. Lines' styles and color codes are as in \cref{fig:PL_MCMC}. In addition, the purple contours display the 1D and 2D posteriors recovered by assuming the injected parameter $\delta=1$ to be known a priori. BPL-BNK\_1 is defined as in eq.~\eqref{eq:master_modelIV_stBBN} with $\log_{10}(h^2\,\Omega_*), f_* /{\rm mHz},\, n_{t,1},\, n_{t,2},\,\delta$ equal to $-10.5, 1,\, 4, 0,\,1$, respectively.\small 1D and 2D posterior distributions derived from the BPL template reconstruction of the benchmark BPL-BNK\_2. Lines' styles and color codes are as in \cref{fig:PL_MCMC}. BPL-BNK\_2 is defined as in eq.~\eqref{eq:master_modelIV_stBBN} with $\log_{10}(h^2\,\Omega_*), f_* /{\rm mHz},\, n_{t,1},\, n_{t,2},\,\delta$ equal to $-9.3, 1,\, 2.65, -2.1,\,5.3$, respectively.\small 1D and 2D posterior distributions derived from the DP template reconstruction of the benchmarks DP-BNK\_1 (top left), DP-BNK\_2 (top right), DP-BNK\_3 (bottom left) and DP-BNK\_4 (bottom right) defined as in \cref{eq:templateDP} with $\{\log_{10}h^2\Omega_* ,\,f_*/\textrm{Hz}\}$ equal to $\{-9.35, 5\}$, $\{-11.35, 5\}$, $\{-9.35, 1\}$ and $\{-11.35, 1\}$, respectively, and $\{\beta , \,\kappa_1,\, \kappa_2,\, \rho, \,\gamma \} = \{0.242, \,0.456,\, 1.234,\, 0.08, \,6.91 \}$ fixed for all benchmarks. Lines' styles and color codes are as in \cref{fig:PL_MCMC}.\small 1D and 2D posterior distributions derived from the ES template reconstruction of the benchmarks ES-BNK\_1 (top left), ES-BNK\_2 (top right) and ES-BNK\_3 (bottom). Lines' styles and color codes are as in \cref{fig:PL_MCMC}. ES-BNK\_1, ES-BNK\_2 and ES-BNK\_3 are defined as in \cref{templateexcited} with $\{\log_{10}h^2\Omega_*,\gamma_{\mathrm{ES}},\omega_{\mathrm{ES}}\mathrm{Hz}\}$ set at $\{-11,10,6 \cdot 10^{3} \}$, $\{-11.5,10,6 \cdot 10^{3} \}$ and $\{-11,10,1.2 \cdot 10^{3}\}$, respectively.\small 1D and 2D posterior distributions derived from the LO template reconstruction of the benchmarks LO-BNK\_1 (left) and LO-BNK\_2 (right). Lines' styles and color codes are as in \cref{fig:PL_MCMC}. LO-BNK\_1 and LO-BNK\_2 are defined as in \cref{eq:sharp-template} with $\{f_*/ \mathrm{mHz},\,\log_{10}(\omega_{\mathrm{lin}} \,\mathrm{Hz}) \}= \{5,\,4\} $ and $\{1,\,4.7\}$ respectively, and fixed $\{\rho, {\cal A}_{\mathrm {lin}}\} = \{ 0.08, 0.3\}$. |