Gamma-ray variability and multi-wavelength insights into the unprecedented outburst from 4C 31.03

The blazar 4C 31.03 recently underwent a major gamma-ray outburst at the beginning of 2023 after a prolonged quiescent phase. Fermi-LAT reported a daily average flux of 5x10^-6 phs cm^-2 s^-1, which is about 60 times its average value. We investigated this extraordinary outbreak through temporal and multi-wavelength analysis. From the statistical analysis of the gamma-ray lightcurves using Bayesian blocks, we identified 3 epochs of prominent flares. The fastest flux decay during this major outburst was observed within 5.5 +/- 0.7 hours. The highest energy of gamma-ray photons found from the source during the active phase is ~ 82 GeV. Using the transparency of gamma-rays against pair production and light crossing time argument, we could obtain the minimum jet Doppler factor as 17 corresponding to the flaring state.The broadband spectral energy distribution study performed using synchrotron, SSC and EC emission processes supports the external Compton scattering of IR photons as the likely mechanism for the gamma-ray emission from the source. The results of this study suggest the scenario of the emission region in 4C 31.03, being located beyond the Broad-line region from the central blackhole. Long-term gamma-ray flux distribution of 4C 31.03 depicts a double log-normal variability, indicating that two distinct flux states are active in this energy band. The index distribution also reveals a two distinct variability patterns, suggesting that the gamma-ray spectrum can be more precisely described by two photon indices.