  Recent experiments have observed large anisotropic collective flows in high
multiplicity proton-lead collisions at the Large Hadron Collider (LHC), which
indicates the possible formation of mini quark-gluon plasma (QGP) in small
collision systems. However, no jet quenching has been confirmed in such small
systems so far. To understand this intriguing result, the system size scan
experiments have been proposed to bridge the gap between large and small
systems. In this work, we perform a systematic study on both heavy and light
flavor jet quenching in different collision systems at the LHC energies. Using
our state-of-the-art jet quenching model, which combines the
next-to-leading-order perturbative QCD framework, a linear Boltzmann transport
model and the (3+1)-dimensional viscous hydrodynamics simulation, we provide a
good description of nuclear modification factor $R_{\rm AA}$ for charged
hadrons and $D$ mesons in central and mid-central Pb+Pb and Xe+Xe collisions
measured by CMS collaboration. We further predict the transverse momentum and
centrality dependences of $R_{AA}$ for charged hadrons, $D$ and $B$ mesons in
Pb+Pb, Xe+Xe, Ar+Ar and O+O collisions at the LHC energies. Our numerical
results show a clear system size dependence for both light and heavy flavor
hadron $R_{AA}$ across different collision systems. Sizable jet quenching
effect is obtained for both heavy and light flavor hadrons in central O+O
collisions at the LHC energies. Our study provides a significant bridge for jet
quenching from large to small systems, and should be helpful for finding the
smallest QGP droplet and the disappearance of QGP in relativistic nuclear
collisions.
