Total: 1
Microquasars (MQs) are Galactic binary systems, consisting of a star and a compact object, a neutron star or a stellar mass black hole, which accretes matter from its companion star and gives rise to relativistic jets. Recent detection of very-high-energy (VHE; $E \gtrsim 100,\text{GeV}$) and ultra-high-energy (UHE; $E \gtrsim 100,\text{TeV}$) gamma-rays by LHAASO, HAWC and HESS from the MQ V4641 Sgr and SS 433 suggests them as Galactic PeVatrons. In this work, we studied a hadronic origin of the observed TeV-PeV gamma-ray emission from these MQs. We considered the hadronic scenario where the gamma-rays are produced by the interaction of relativistic protons in the MQ jet with the stellar wind. We fitted our model with observed data and constrained physical parameters like the hadronic jet power fraction, the proton spectral index, the maximum proton energy and the jet bulk Lorentz factor. Our best-fit model shows hard proton spectra ($1.84-2.44$) and maximum proton energies between 1 and 5 PeV. We also estimated the all-flavor neutrino fluxes corresponding to the gamma-ray fluxes from the hadronic model and found that V4641 sgr can be detected by next-generation neutrino telescopes like KM3NeT-ARCA and TRIDENT. Furthermore, we modeled a synthetic population of Galactic MQs and estimated their contribution to the diffuse TeV-PeV gamma-ray flux. For the inner Galaxy PSR contribution dominates in the range 10-100 TeV, and above 100 TeV diffused cosmic ray interactions with the molecular clouds is most dominant. We find that a population $\sim 14$ MQs is required to explain the LHAASO data above 100 TeV. For the outer Galaxy, we show that MQs are the dominant class of sources, and we constrain their population $\sim$14. Our findings strongly suggest that MQs are efficient particle accelerators, contributing to Galactic PeVatrons and potential multimessenger sources in our Galaxy.