Origin of the broadband emission from the transition blazar B2 1308+326

Transition blazars exhibit a shift from one subclass to the next during different flux states. It is therefore crucial to study them to understand the underlying physics of blazars. We probe the origin of the multi-wavelength emission from the transition blazar B2 1308+326 using 14-year-long gamma-ray light curve from Fermi and the quasi-simultaneous data from Swift. We used the Bayesian block algorithm to identify epochs of flaring and quiescent flux states and modelled the broadband SEDs for these epochs. We employed the one-zone leptonic model in which the synchrotron emission causes the low-energy part of the SED and the high-energy part is produced by the IC emission of external seed photons. We also investigated its multi-band variability properties and gamma-ray flux distribution, and the correlation between optical and gamma-ray emissions. We observed a historically bright flare from B2 1308+326 across the optical to gamma-ray bands in June and July 2022. The highest daily averaged gamma-ray flux was (14.24$\pm$2.36) $\times$ 10$^{-7}$ ph cm$^{-2}$ s$^{-1}$ and was detected on 1 July 2022. The gamma-ray flux distribution was found to be log-normal. The optical and gamma-ray emissions are well correlated with zero time lag. The synchrotron peak frequency changes from $\sim 8 \times$ 10$^{12}$ Hz (in the quiescent state) to $\sim 6 \times$ 10$^{14}$ Hz (in the flaring state), together with a decrease in the Compton dominance providing a hint that the source transitions from a LSP to an ISP. The SEDs for these two states are well-fitted by one-zone leptonic models. The parameters in the model fits are essentially consistent between both SEDs, except for the Doppler-beaming factor, which changes from $\sim$15.6 to $\sim$27 during the transition. An increase in the Doppler factor might cause both the flare and the transition of B2 1308+326 from an LSP to an ISP blazar.


Introduction
Blazars are active galactic nuclei with relativistic jets, whose trajectories are closely aligned to the observer's line of sight (Urry & Padovani 1995;Padovani 2017).The main characteristics of blazars are the Doppler-boosted non-thermal emission from the jets, the high-amplitude flux variability over the entire electromagnetic spectrum, and radio-to-X-ray polarizations (e.g.Wagner & Witzel 1995;Pandey et al. 2022;Rajput et al. 2022;Liodakis et al. 2022).Blazars are categorized into flatspectrum radio quasars (FSRQs; EW 1 rest > 5Å) and BL Lacertae objects (BLLs; EW rest < 5Å) on the basis of their optical/ultraviolet (UV) emission line properties (e.g.Stocke et al. 1991;Marcha et al. 1996).The observed non-thermal broadband spectral energy distributions (SEDs) of blazars show doublehump structures.The peak of the low-energy hump lies in infrared to X-ray energies, and that of the high-energy part ranges from GeV to TeV energies (e.g.Fossati et al. 1998).The low-1 equivalent width of the emission lines in rest frame frequency hump is interpreted as the synchrotron emission produced by relativistic electrons within the jet.On the other hand, two different models, leptonic and hadronic models, have been proposed to explain the high-frequency hump.In the leptonic scenario, the high-energy part of the blazar SED is caused by the inverse-Compton (IC) scattering of either the synchrotron photons (synchrotron-self Compton, SSC; e.g.Ghisellini et al. (1985); Bloom & Marscher (1996)) or by the external photons (external Compton, EC; e.g.(Begelman et al. 1987;Sikora et al. 1994)) by the same electrons producing synchrotron emission.In hadronic models, in contrast, the high-energy γ−ray radiation of a blazar is attributed to hadronic processes such as proton and muon synchrotron emissions (e.g.Mücke & Protheroe 2001;Böttcher et al. 2013;Petropoulou & Mastichiadis 2015).
Based on the low-energy component peak (ν p syn ) of their SEDs, blazars are further classified into low-synchrotron peaked (LSPs; ν p syn ≤ 10 14 Hz ), intermediate-synchrotron peaked (ISPs; 10 14 < ν p syn < 10 15 Hz), and high-synchrotron peaked (HSPs; ν p syn ≥ 10 15 Hz) blazars (Abdo et al. 2010).FSRQs belong to the LSP class, and BL Lacs can be LSPs, ISPs, or HSPs based on their ν p syn .It has been found that the blazar population generally follows an empirical trend, known as the blazar sequence, such that ν p syn correlates with the γ−ray peak frequency and is anticorrelated with the Compton dominance (CD; the ratio of the γ−ray peak luminosity to the low-energy peak luminosity (Fossati et al. 1998;Ghisellini et al. 2017;Prandini & Ghisellini 2022)).A possible explanation for the blazar sequence is the difference in the electron cooling efficiency, as proposed by Ghisellini et al. (1998).The blazar sequence can also be interpreted as an artefact of the differences in Doppler boosting (e.g.Nieppola et al. 2008;Fan et al. 2017).
In addition to these conventional classifications, certain blazars showed characteristics of both FSRQs and BL Lacs during their different flux states (e.g.Ghisellini et al. 2011Ghisellini et al. , 2013;;Ruan et al. 2014;Mishra et al. 2021).Blazars that exhibit a transition from FSRQ to BL Lac or vice versa are known as transition blazars.These blazars can be identified by investigating the shape of their broadband SEDs (e.g.Ghisellini et al. 2011Ghisellini et al. , 2013) ) and/or by estimating the EW of the broad emission lines in their optical/UV spectra (e.g.Ruan et al. 2014;Mishra et al. 2021).

B2
1308+326 (OP 313) is a high-redshift (z=0.9980±0.0005;Hewett & Wild (2010)) blazar.It has been observed several times at different wavelengths because of its variable emissions and uncertain classification (e.g.Gabuzda et al. 1993;Watson et al. 2000;Britzen et al. 2017).It was initially classified as a BLL (Stickel et al. 1991) due to its nearly featureless optical spectra (Wills & Wills 1979), high optical polarization, and extreme optical variability (Angel & Stockman 1980).However, using very long baseline interferometry (VLBI) polarization images, Gabuzda et al. (1993) observed the polarized flux from the inner part of its jet with a position angle perpendicular to the jet.In addition, they also detected tentative superluminal motion in its VLBI jet.Because these are characteristics of a quasar, they classified B2 1308+326 as an FSRQ with unusually weak emission lines.Using their VLBI data together with the optical data of Stickel et al. (1991), they suggested that B2 1308+326 might be a gravitationally microlensed quasar.However, Urry et al. (1999) did not detect any spatially extended emission from B2 1308+326 in the high-resolution imaging observations carried out with HST WFPC2.
The blazar B2 1308+326 was observed simultaneously at X-ray, optical, and radio wavelengths by Watson et al. (2000).They found that it could be a radio-selected BLL based on its optical variability and synchrotron peak power, but because of its high bolometric luminosity, high Doppler factor, and variable line emission, it appears more likely to be an FSRQ.They concluded that it may be an intermediate or a transitional blazar or a gravitationally microlensed quasar (due to excess absorption at X-rays).In the γ−regime, B2 1308+326 was detected by the Large Area Telescope (LAT) on board the Fermi telescope (Ackermann et al. 2013).It is designated as an FSRQ in the Fermi-LAT Fourth Source Catalog (4FGL; Abdollahi et al. (2020)).
In this work, we investigate the physical processes that cause the broadband emission of the blazar B2 1308+326 using optical-to-γ−ray data for a period of ∼ 14 yr.The paper is organized as follows.In Section 2 we discuss the steps we used to reduce the multi-wavelength data.Section 3 presents the results of this study.The details of the SED modelling are given in Section 4. A discussion of our results and our conclusions is given in section 5. Our findings are summarized in Section 6.

Fermi-LAT data
We used the γ−ray data of B2 1308+326 measured with the Fermi-LAT between 4 August 2008, and 12 December 2022 covering ∼ 14 yr of Fermi operations.We adopted the standard LAT data analysis procedures2 to perform the data analysis.The data were reduced using the FermiTool version 2.2.0 and the Fermipy version 1.2 (Wood et al. 2017).We selected the Pass 8 Data (P8R3) in the energy range of 0.1−500 GeV and considered all the SOURCE class events (ev-class=128 and evtype=3) in the region of interest (ROI) of 10 • × 10 • centred on the target source position (RA: 197.619,Dec: 32.3455).We applied the standard cuts (zenith angle, z max < 90 • and "(DATA_QUAL>0)&&(LAT_CONFIG==1)") to obtain the good time intervals.Because at low energy (E<100 MeV), the effective area of Fermi-LAT rapidly decreases and the point-spread function (PSF) increases (e.g. at 100 MeV, the 68% containment angle of the acceptance-weighted PSF is ∼5 degrees3 ) we included all the point sources from the 4FGL within 15 • of the ROI centre, along with the Galactic (gll_iem_v07 ) and extragalactic isotropic diffuse emission ( iso_P8R3_SOURCE_V3_v1 ) components in our initial model file.To optimize the spectral parameters, we performed a binned likelihood analysis with eight energy bins per decade and a spatial binning of 0.1 • per pixel.
After initial optimization, we removed sources with test statistics (TS) < 1 and allowed the spectral normalization of sources with TS > 10 to vary in the model file.We also set the spectral shape parameters free to vary for sources within 3 • of the ROI center.The normalizations of the Galactic and extragalactic diffuse components were left free, together with the spectral index of the Galactic diffuse component.We fitted the ROI again, and when the fit converged successfully, we used the best-fitting model file to generate the γ−ray light curve and spectra.For the SED analyses, we used a larger ROI of 15 degrees to avoid the systematics.

Swift-XRT data
The X-ray light curve and spectra of B2 1308+326 were generated using the online Swift-XRT data products generator tool 4 .The details of the data reduction process followed by this online tool can be found in Evans et al. (2007Evans et al. ( , 2009)).
The 0.3-10 keV X-ray spectra of B2 1308+326 were first grouped to a minimum of 20 counts per bin using the GRP-PHA task of FTOOLS and were then fitted with an absorbed power-law (PL) model (tbabs*powerlaw) in the XSPEC version 12.12.0.While fitting, we fixed the Galactic hydrogen column density, n H , to 1.22 × 10 20 cm −2 (Willingale et al. 2013).

Swift-UVOT data
During the Swift monitoring, the UVOT instrument observed the FSRQ B2 1308+326 in its three optical (v, b, and u) and three ultraviolet (uvw1, uvm2, and uvw2) filters.We downloaded the UVOT data from the HEASARC Data archive 5 and performed the data analysis using the HEASoft package version 6.29 and the CALDB version 20211108.
We first summed the multiple observations taken in the same filter over a given epoch using the task uvotimsum and then extracted the source magnitudes using the task uvotsource Ṫhe source magnitudes were extracted from a circular region of radius 5 arcsec centred on the source, while the background magnitudes were derived from a source-free circular region of radius 20 arcsec.The magnitudes were corrected for the Galactic extinction using the E(B − V) value of 0.0115 taken from Schlafly & Finkbeiner (2011) and the extinction laws from Cardelli et al. (1989).The reddening-corrected magnitudes were then converted into flux densities using the zero points given in Breeveld et al. (2011).

Archival data
We also used the publicly available g-band optical data of B2 1308+326 from the ASAS-SN6 (All Sky Automated Survey for Supernovae) data archive (Shappee et al. 2014;Kochanek et al. 2017) to investigate the correlation between optical and γ−ray emissions.

Flare selection
The weekly binned γ−ray light curve of B2 1308+326 is shown in Figure 1.We show upper limits at the 95% confidence limit for the time bins within which either the source is detected with TS < 10 or the flux is lower than or equal to its error F i ≤ σ i .The average γ−ray flux F = (0.57± 0.05)× 10 −7 ph cm −2 s −1 , estimated using the maximum likelihood analysis over the entire ∼ 14 yr monitoring, is shown as a horizontal red band.The source exhibits multiple episodes of high γ−ray activity, with the highest weekly averaged γ−ray flux of 10.67 ×10 −7 ph cm −2 s −1 recorded recently on MJD 59761.
There is no generally accepted way of defining a flare in a light curve.However, the Bayesian blocks (BB) algorithm (Scargle et al. 2013) has recently been used by several authors to delineate flares in blazar light curves (e.g.Meyer et al. 2019;Sahakyan et al. 2022).The algorithm provides a block-wise constant representation of a light curve by identifying the statistically significant variations.We adopted the astropy 7 implementation of the BB algorithm with a false-alarm probability p 0 = 0.05 to identify the flares in the γ−ray light curves of B2 1308+326.Because the BB algorithm can be applied to nonuniform light curves, we ignored the upper limits while applying the BB method to the light curves.The BB representation of the weekly binned γ−ray light curve of B2 1308+326 is shown as solid blue lines in Figure 1, where several flares exceeding the mean flux level are easily recognized.To estimate the duration of the flares and to treat the overlapping flares, we used the HOP algorithm (Eisenstein & Hut 1998;Meyer et al. 2019).The algorithm includes two steps: (1) It identifies a block with a value F BB higher than both the preceding and following blocks as a peak, and (2) it moves downward in both directions from the peak and includes blocks with a value F BB ≥ F into the peak to form a HOP group.Using the algorithm, we iteratively searched the γ−ray light curve for HOP groups and arbitrar-ily selected only those whose peak value F BB ≥ 5 × F for the weekly binned light curve (e.g.Meyer et al. 2019).In this way, we found only one HOP group (outburst phase), which spans MJD 59460.16 to 59922.05 (∼ 462 days).For this duration, we reoptimized the spectral model following the steps described in Section 2.1 and extracted the one-day-binned light curve that is shown in Figure 2. We searched the one-day-binned light curve using the BB and HOP algorithm for the flare identification and found one flare that satisfied the arbitrary condition F BB ≥ 2 × F, which covers the period from MJD 59733.66 to MJD 59786.16(∼ 52.5 days).During the flare, the highest daily averaged γ−ray flux of (14.24±2.36)× 10 −7 ph cm −2 s −1 was observed on MJD 59761.65 (1 July 2022).

Multi-band flux variability
The γ−ray to optical energy band light curves of B2 1308+326 are plotted in Figure 3.The figure shows that B2 1308+326 exhibits strong flux variations at all these different energy bands.To quantify these flux variations, we estimated the fractional variability amplitude (e.g.Vaughan et al. 2003;Pandey et al. 2017, and references therein), which is defined as follows: The uncertainty in F var can be determined as where S 2 is the sample variance of the light curve, x is the mean flux, and σ 2 err is the mean square error.The values of F var together with its uncertainty for the different energy band light curves are given in Table 1.At low energies (optical/UV), the blazar B2 1308+326 showed more flux variations than at high energies (X-ray/γ−ray).
We also calculated the shortest flux doubling/halving timescales during the flaring period, as follows: where F(t 1 ) and F(t 2 ) denote the fluxes at times t 1 and t 2 , respectively, ∆t = t 2 − t 1 , and τ represents the flux doubling or halving timescale.The shortest γ−ray flux halving timescale is found to be (140.86±41.15)hr between MJD 59761.66 and 59770.66,while at the optical (g) band, we found a shorter flux doubling timescale of (80.71±2.53)hr between MJD 59757.25 and 59759.32.

Highest-energy photon
We also extracted the arrival time, energy, and probability of the highest-energy photons coming from the source using the tool gtsrcprob on the ULTRACLEAN (evclass=512) event class and a 0.5 • ROI.We plot the energies of the photons with a probability >99% against their arrival times in the bottom panel of

Gamma-ray flux distribution
We investigated the distribution of weekly averaged γ−ray fluxes of B2 1308+326 for bins with TS ≥10 and F i > σ i .In this way, we neglected the upper limits of the γ−ray flux, which are mostly below 0.5 ×10 −7 ph cm −2 s −1 and may form a disjoint population in the distribution.We obtained upper limits instead of a flux value in ∼ 64% time bins for the weekly binned light curve.We first performed the Anderson-Darling (AD) test.The null hypothesis for the AD test is that the sample follows a normal distribution.We applied the AD test on the flux and log8 flux distributions and estimated the p−values for the flux and log-flux distributions, which are 1.3e-51 and 0.16, respectively.A p−value < 0.01 indicates a clear deviation from a nominally normal distribution.Therefore, the γ−ray flux distribution of B2 1308+326 does not follow a normal distribution on a linear scale.However, on a logarithmic scale, the flux distribution is consistent with normal, implying that the γ−ray flux distribution of B2 1308+326 is essentially log-normal.
We also constructed the normalized histograms of flux values on linear as well as logarithmic scales, shown in Figure 4, and fitted them with a Gaussian function using a χ 2 -fit.The reduced χ 2 values for the Gaussian fit to the flux distribution on linear and logarithmic scales are 2.98 and 1.05, respectively, which confirms that the γ−ray flux distribution of B2 1308+326 is log-normal.A similar log-normal distribution was observed by Meyer et al. (2019) for weekly binned γ−ray light curves of all the six blazars they studied.In Figure 4, we also show the distribution of the upper limits as dashed bars.

Correlation analysis
We investigated the correlations between the γ−ray and optical emissions of B2 1308+326 using the discrete correlation function (DCF; Edelson & Krolik (1988)).The DCF method has been extensively used to search for correlations between nonuniform light curves at different wavebands (e.g.Cohen et al. 2014;Pandey et al. 2017;Raiteri et al. 2021).The peaks in the DCF plot denote the correlations whose significance increases as the peak value increases towards unity.The result of DCF analysis between the one-day-binned γ−ray and optical g−band light curves is plotted in Figure 5, where a positive lag would indicate that the γ−rays lead the optical emission.The strong positive DCF peak at about zero time lag indicates that the γ−ray and optical emissions are strongly correlated.

Origin of the multi-wavelength emission
The multiwavelength data analysed in Section 2 allow us to generate the broadband SEDs of B2 1308+326 in different flux states.We generated the SEDs for the following epochs: (b) Flaring state (from MJD 59754 -59770 ): When the source was in the bright state in all the bands.This period corresponds to a Bayesian block within the flaring duration given by the BB algorithm.
The broadband SEDs for these epochs are shown in Figure 6, which depicts two highly intriguing characteristics.First is a change in the Compton dominance; the γ−ray peak dominates the synchrotron peak in the low state, while in the flaring state, both peaks are similar.Second is an indication of a shift of the optical-UV synchrotron peak, where the lower state SED is like that of an FSRQ, while in the flaring state, it is closer to that of an ISP.
To investigate and comprehend the emission mechanism and potential cause of the peak shift, we modelled the broadband SEDs of B2 1308+326 using a one-zone leptonic scenario.In this model, the emission originates from a spherical blob of particles with a radius, R, within the jet.The blob is filled with a uniform magnetic field, B, and moves down the jet with a bulk Lorentz factor of Γ at a small angle θ to the observer, such that the Doppler factor δ ≃ Γ.The blob is assumed to be filled with relativistic electrons with a broken power-law distribution of the form where p and p 1 are the indices below and above the break energy γ br .γ min and γ max are the minimum and maximum electron Lorentz factors, respectively, and N e is the normalization.The electrons interact with the magnetic field and produce synchrotron emission, which causes the low-energy component of the SED.The HE part of the SED is generated by the IC scattering of seed photons by the same population of electrons.When the seed photons are the synchrotron photons, the HE component of the SED is explained by the SSC process (e.g.Bloom & Marscher 1996).On the other hand, if the seed photons come from the external fields such as an accretion disk, BLR, and/or dusty torus (DT) the HE component is caused by ECdisk, EC-BLR, and/or EC-DT, respectively (e.g.Dermer et al. 1992;Sikora et al. 1994;Błażejowski et al. 2000).
The exact location of the γ−ray emitting region in blazar jets is still unclear.However, we can constrain the location of the emission region, R diss , using the observed minimum variability timescale and the energy of the highest-energy photon (e.g.Prince 2020).We detected a minimum variability timescale of ∼ 81 hr, which corresponds to a region size of The R diss is then given by R diss = R/ψ, where ψ is the semiaperture angle of the jet, whose values generally range between 0.1 and 0.25 (e.g.Dermer et al. 2009;Ghisellini & Tavecchio 2015).Radio observations have also reported values of ψ < 0.1, but radio features do not occupy the whole section of the jet (e.g.Jorstad et al. 2005).Here, we assumed a typical value of ψ = 0.1 (Ghisellini & Tavecchio 2015).It is important to note that the BLR is opaque for the high-energy (E > 20 GeV/(1+z)) γ−rays because they are absorbed via γγ pair-production (Liu & Bai 2006;Stern & Poutanen 2014).The maximum energy of a γ−ray photon detected during the flaring state is 27 GeV, which implies that the R diss location is probably outside the BLR.
For the broadband SED fitting of B2 1308+326, we assumed that the emission region is either at the outer edge of the BLR or outside the BLR and the HE part of the SED is represented by the SSC + EC-BLR + EC-DT.We considered BLR as a spherical shell with inner and outer radii of R BLR and 1.1 × R BLR (Kaspi et al. 2007), respectively, where R BLR = 10 17 (L disk /10 45 ) 0.5 cm and L disk = 9 × 10 45 erg/s (Ghisellini & Tavecchio 2015).The distance to the dusty torus was assumed to be R DT = 2 × 10 18 (L disk /10 45 ) 0.5 cm, and its temperature was fixed at a typical value of 800 K (Ghisellini & Tavecchio 2015).
We performed the SED modelling using the publicly available code JetSeT (Massaro et al. 2006;Tramacere et al. 2009Tramacere et al. , 2011;;Tramacere 2020).The free parameters were p, p 1 , γ min , γ max , γ br , N e , B, R, and δ, which were constrained during the fitting.The Minuit optimizer was used to limit the parameters initially, and subsequently, MCMC sampling of their distributions was used to enhance them.The values of the best-fitted model parameters for both quiescent and flaring SED modelling are given in Table 2.The optimal model SEDs for both states are plotted in Figure 7.
In the quiescent state, the synchrotron radiation from the jet dominates the optical/UV band, the SSC contributes most in the X-ray band, and the EC-DT produces the γ−rays (see the left panel of Fig. 7).The SED modelling shows that the slope of the electron distribution changes from 1.95 to 3.24 at the break energy of 1.01 ×10 3 .The minimum energy of electrons γ min is ∼ 49.61, which is a typical value for LSPs (Kang et al. 2014).The magnetic field is found to be 0.14±0.01G.The Doppler factor is 15.61, which is consistent with the value found by (Ghisellini et al. 2010).The emitting region size is 2.32 ×10 17 cm.The synchrotron and IC peak frequencies are ∼ 8×10 12 and ∼ 5×10 22 Hz, respectively.The value of the CD is 2.93.
In the flaring state, the synchrotron radiation still dominates the optical/UV band, but the X-ray band receives contributions from both the synchrotron and SSC emissions.The SSC emission extends up to the low-energy Fermi-LAT band, while the high-energy γ−ray data are fitted by EC-DT (see the right panel of Fig. 7).The power-law slope of the electron distribution at low energy is slightly harder than that in the quiescent state, but the high-energy slope has increased to 5.14 with a higher break energy of 6.82 ×10 3 .The minimum electron Lorentz factor γ min is only ∼ 3, indicating that electrons with low energy are also effectively accelerated.The strength of the magnetic field is marginally (0.29±0.02G) higher than that in the quiescent state, but the Doppler factor has now increased to 26.97.The radius of the emitting region is 7.49 ×10 16 cm, which corresponds to a variability timescale of t var = R(1 + z)/δc ∼ 2 days.The peak frequencies for the synchrotron and IC components are ∼ 6×10 14 and ∼ 2×10 24 Hz, respectively.The value of the CD is 1.39.

Discussion and conclusions
We extensively investigated the multi-wavelength emission from the transition blazar B2 1308+326 with a focus on its most recent outburst in 2022 using data from Fermi, Swift, and ASAS-SN.The blazar showed flux variations in each of the studied energy bands.The maximum daily averaged γ−ray flux of (14.24±2.36)× 10 −7 ph cm −2 s −1 was detected on MJD 59761.65 during the 2022 flare.The flux variability amplitudes are higher at lower energies.The derived minimum variability timescale in the optical band is ∼81 hr (3.4 days), which is shorter than the observed minimum variability timescale of ∼141 hr (5.9 days) in the γ−ray.The inability to identify faster variability in the γ−ray band may, however, be due to the relatively large error bars.The measured minimum variability timescales are consistent with those observed from the SED modelling.
We observed that the γ−ray flux of B2 1308+326 follows a lognormal distribution.The study of the flux distribution is a unique tool for probing the underlying physical processes.A Gaussian or normal flux distribution, which is expected for linear stochastic variations, suggests additive processes.For non-linear stochastic variations, the flux distribution is log-normal, which more naturally arises from multiplicative processes (Uttley et al. 2005).Log-normal flux distributions are often found in blazars (e.g.Giebels & Degrange 2009;Kushwaha et al. 2017).They can be interpreted as the result of multiplicative processes that initially originate in the accretion disk and then propagate to the jet (McHardy 2008), or as the sum of emission from the randomly oriented mini-jets within the jet (Biteau & Giebels 2012).Variations in the particle acceleration or escape timescales can also produce a log-normal flux distribution (Sinha et al. 2018).The strong correlation at zero lag found between optical and gammaray radiations suggests that they have a common spatial origin.
To understand the origin of the broadband emission of B2 1308+326, we generated SEDs in the quiescent and flaring states and fitted them with the one-zone leptonic model.Both SEDs are well represented by the model, defining the low-energy component as synchrotron emission and the high-energy compo-  nent as the sum of SSC and EC-DT.By comparing the bestfit model parameters for the two SEDs, we found that the electron distribution has a similar low-energy power-law index, but the high-energy indices for quiescent and flaring states are ∼ 3.24 and ∼5.14, respectively.The differences between the spectral indices, ∆p = p 1 − p in the two states are 1.29 and 3.39, respectively.The expected value of ∆p for a standard radiative cooling is 0.5 (Kardashev 1962).However, Reynolds (2009) suggested that such a large ∆p > 1 could be naturally produced due to inhomogeneities in the source.The break energy is higher, γ br ∼ 6823, during the flare than in the quiescent state, γ br ∼ 1009, indicating that the particles are accelerated to higher energies during the flare.In the two states, the strength of the magnetic field is comparable, but the Doppler factor has increased from ∼ 15.6 to ∼27 during the flare.
We observed that the peak frequencies for the synchrotron and IC components were pushed to higher values during the flare.In the quiescent stage, the synchrotron peak frequency was ∼ 8× 10 12 Hz, suggesting that the source was an LSP, but it shifts to ∼ 6× 10 14 Hz during the flare, indicating a source transition from LSP to ISP.We also noted a shift in IC peak frequency from ∼ 5× 10 22 Hz (quiescent) to ∼ 2× 10 24 Hz (flare).A similar shift in the synchrotron peak frequency of B2 1308+326 during a high flux state has also been reported by Watson et al. (2000).This shift in the synchrotron peak frequency during the flare has been observed in several blazars (e.g.Cavaliere & D'Elia 2002;D'Ammando et al. 2011;Giommi et al. 2012;Cutini et al. 2014) The transition from one blazar subclass to the next has also been reported using spectroscopic data (e.g.Shaw et al. 2012;Ruan et al. 2014;Peña-Herazo et al. 2021).The shift from LSP to ISP/HSP can be attributed to (i) variations in the Doppler factor alone (e.g.Bianchin et al. 2009); (ii) hiding of the broad emission lines by overwhelming synchrotron emission that peaks in the UV in the sources with radiatively weak cooling (Ghisellini et al. 2012); and (iii) swamping of broad emission lines by the variability in the jet continuum emission in the sources with radiatively efficient accretion flows and strongly beamed jets (e.g.Giommi et al. 2012;Ruan et al. 2014).The results of our broadband SED modelling suggest that the transition of B2 1308+326 from LSP to ISP is due to an increase in the Doppler factor, which also explains the shift in the IC peak and the decrease in the CD.The change in the Doppler factor can be explained by geometrical effects, such as the change in the viewing angle of the emitting region within the jet (e.g.Villata & Raiteri 1999;Raiteri et al. 2017).An inhomogeneous and helically curved jet can undergo orientation changes caused by the magnetohydrodynamic instabilities or rotation of the twisted jet.This leads to variations in the viewing angle toward jet-emitting regions, and hence, to variations in the Doppler factor.When the orientation of the emitting region is closely aligned to the observer, the emission from it is more strongly Doppler boosted.As a result, a flare can be seen.Using long-term Very Long Baseline Array (VLBA) observations, Britzen et al. (2017) found that the jet of B2 1308+326 has a helical structure.The 2022 flare of B2 1308+326 may therefore be caused by a compact emitting region viewed at smaller viewing angles as compared to the entire jet, which enhances the Doppler factor (e.g Raiteri et al. 2017).
A similar transitional behaviour was also observed in the blazar PMN J2345−1555 by Ghisellini et al. (2013).It is an FSRQ source with a synchrotron component that generally peaks in the far-IR region.However, during the January 2013 flare, the synchrotron peak moved to the optical-UV frequencies, which changed it to a BL Lac source.This transitional behaviour was interpreted as due to the change in the location of the dissipation region, R diss , from within the BLR (during the low state) to just outside the BLR (during the flare).The corresponding decrease in the radiative cooling allowed the relativistic particles to reach higher energies, resulting in a shift in the synchrotron and IC component peaks.In our case, the broadband SED modelling suggested that the flare and the shift in the peaks were due to the change in the Doppler factor, which is plausibly caused by the change in the orientation of the emitting region.Although the location of the dissipation region was shifted from 10.01 ×10 17 cm (in the low state) to 4.05 ×10 17 cm (during a flare), it was outside the BLR in both cases.
Modelling the SED also provides information about the jet luminosity.The total jet luminosity is defined as (Celotti & Ghisellini 2008) L total = L e + L p + L B + L r = πR 2 Γ 2 c(U e + U p + U B + U r ), (6) where L e (U e ), L p (U p ), L B (U B ), and L r (U r ) are the luminosities (energy densities) of the electron, proton, magnetic field, and radiation, respectively.On the assumption of one proton per relativistic electron and cold protons, the value of the luminosity for each component and the total jet luminosity are listed in Table 2 for both SEDs.The values of L e /L B in the quiescent and flaring states are ∼ 4.29 and ∼3.63, respectively, indicating that the system is close to equipartition.Compared to the quiescent state, the overall jet power is marginally higher when the source is flaring.Additionally, in the flaring state, the overall jet luminosity exceeds the disk luminosity, whereas in the quiescent state, it is comparable to the disk luminosity.Ghisellini et al. (2014) found that the total jet power of blazars can be higher than their disk luminosities.Our results concur with their findings.For a black hole mass of 5.25 × 10 8 M ⊙ (Zamaninasab et al. 2014), the Eddington luminosity for B2 1308+326 is L Edd ≃ 6.61 × 10 46 erg/s.The total jet power in the quiescent state is ∼ 14% of L Edd , while it is ∼ 35% of L Edd in the flaring state.

Summary
We have carried out a multi-wavelength analysis of the highredshift blazar B2 1308+326 for around ∼14 yr.Our main findings are summarized as follows.
-B2 1308+326 exhibited an historically large flare in opticalto-gamma-ray frequencies in June and July 2022, which reached a maximum daily averaged γ−ray flux of 1.42 × 10 −6 ph cm −2 s −1 on MJD 59761.65 (1 July 2022).-The estimated fractional variability amplitude was larger at low frequencies (optical/UV) than at high frequencies (X-ray/γ−ray).The detected variability timescale was also shorter (∼3.4 days) at the optical frequency than the variability timescale (∼6 days) observed at γ−rays.-The γ−ray flux distribution of the source followed a lognormal distribution.-The optical and γ−ray emissions are positively correlated without any detectable time lag.-The synchrotron peak frequency increased from ∼ 8×10 12  Hz (in the quiescent state) to ∼ 6×10 14 Hz (in the flaring state) with a corresponding decrease in the CD indicating that the source changes from LSP to ISP during the flare.-The SED modelling suggested that the transition was most likely produced by an increase in the Doppler factor.

Fig. 4 .Fig. 5 .
Fig. 4. Gamma-ray flux distribution of B2 1308+326 for the entire ∼ 14 yr monitoring period.Left panel: Flux distribution fitted with a Gaussian function on a linear scale (normal distribution).Right panel: Flux distribution fitted with a Gaussian function on a logarithmic scale (log-normal distribution).The results of the fits are mentioned in each plot.The distribution of the upper limits is shown as dashed bars.

Fig. 6 .
Fig. 6.Broadband SEDs of B2 1308+326 in the quiescent and flaring states specified in the text.

Fig. 7 .
Fig. 7. Model-fit broadband SEDs of B2 1308+326 in the quiescent (left panel) and flaring (right panel) states.The solid blue curve represents the sum of all the components, while the grey shaded region denotes the uncertainty region from MCMC sampling of the model parameters.The legends provide the colour-coding for the other components.

Table 1 .
Fractional variability amplitudes (in percent) at different energy bands with a one-day binning.The value of F var for the ASAS-SN g band is estimated for the period of MJD 58458.65-59901.49.For the same duration, the F var for γ−rays is 67.42±2.31%.

Table 2 .
Best-fit model parameters for the broadband SEDs of B2 1308+326 in the quiescent and flaring states.