Black Hole Stars Across the Universe: Identifying Central Engine Dominated Little Red Dots at $z\sim1.5-9.5$

Photometric selections of Little Red Dots (LRDs) largely rely on identifying their ``V-shaped'' spectral energy distribution (SED). Recent work suggests this V-shape stems from a combination of a central engine -- also referred to as a Black Hole Star (BH*) -- and a star-forming host galaxy. We present a new and highly complementary photometric selection that is based on incorporating BH* templates in the \texttt{eazy} redshift fitting code. Selecting compact sources where a BH* template contributes $>80$\% to the best fitting SED in the rest-optical, we compile a sample of 241 BH*-dominated candidates from $\sim1000\,{\rm arcmin}^2$ of legacy and pure parallel JWST imaging. Our selection does not require a blue UV-component, and it successfully identifies objects that resemble the paradigmatic sources ``MoM-BH*-1'' and ``The Cliff''. We find that BH*-dominated sources exist across a wide range of redshifts ($z\sim1.7-9.3$) and optical luminosities (log$(L_{5100}/{\rm erg}\,{\rm s}^{-1})\sim42-44.5$), and we measure a median Balmer break strength of $\sim3$, with some breaks reaching values $>10$. We estimate bolometric luminosities in the range log$(L_{\rm bol}/{\rm erg}\,{\rm s}^{-1})\sim42-45$, which, assuming accretion at the Eddington-limit, would translate to black hole masses of $M_{\rm BH}\sim10^4-10^7{\rm M_\odot}$, spanning the intermediate mass black hole to the quasar regime. The number density of BH*-dominated candidates peaks at $z\sim5-6$ ($\sim10^{-5}\,{\rm Mpc}^{-3}$) and it declines by an order of magnitude down to $z\sim2$. Tentatively, comparing to V-shaped LRD samples suggests that the fraction of BH*-dominated sources among the broader LRD population does not decrease towards lower redshift. Crucially, our work demonstrates that BH*-dominated sources are not merely an early-Universe phenomenon but rather persist at least until cosmic noon.