Accretion and heating rates of PBHs with varying masses in the density vs. velocity plane. The color-coded panels illustrate the behavior of PBHs with varying masses (1 M_⊙, 33 M_⊙, and 100 M_⊙ columns from left to right) in the density vs. velocity plane, depicting the accretion rate (upper row) and the resulting heating rate (bottom row). The dashed white lines delineate the upper limits of specific accretion regimes (eADAF, ADAF, and LHAF) as shown in Fig. 1. The red rectangles enclose the density and characteristic velocity ranges around each PBH within the halos, employing a rescaled version of the gas density profiles from Liu et al. (2022) and the CIELO particle velocities. The blue rectangles encapsulate weighted values around PBHs in the IGM. The dashed yellow line marks the peak density of a gas particle in the CIELO simulation. These results are representative of gas conditions at z ∼ 23.