Insights from the first flaring activity of a high synchrotron peaked blazar with X-ray polarization and VHE gamma rays
Authors:   MAGIC Collaboration:, Abe K. , Abe S. , Abhir J. , Abhishek A. , Acciari V.A. , Aguasca-Cabot A. , Agudo I. , Aniello T. , Ansoldi S. , Antonelli L.A. , Arbet Engels A. , Arcaro C. , Asano K. , Babić A. , Barres de Almeida U. , Barrio J.A. , Barrios-Jiménez L. , Batković I. , Baxter J. , Becerra Gonz'alez J. , Bednarek W. , Bernardini E. , Bernete J. , Berti A. , Besenrieder J. , Bigongiari C. , Biland A. , Blanch O. , Bonnoli G. , Bošnjak Ž. , Bronzini E. , Burelli I. , Campoy-Ordaz A. , Carosi A. , Carosi R. , Carretero-Castrillo M. , Castro-Tirado A.J. , Cerasole D. , Ceribella G. , Chai Y. , Chilingarian A. , Cifuentes A. , Colombo E. , Contreras J.L. , Cortina J. , Covino S. , D'Ammando F. , D'Amico G. , Da Vela P. , Dazzi F. , De Angelis A. , De Lotto B. , de Menezes R. , Delfino M. , Delgado J. , Delgado Mendez C. , Di Pierro F. , Di Tria R. , Di Venere L. , Dinesh A. , Dominis Prester D. , Donini A. , Dorner D. , Doro M. , Eisenberger L. , Elsaesser D. , Escudero J. , Fariña L. , Foffano L. , Font L. , Fröse S. , Fukazawa Y. , García López R.J. , Garczarczyk M. , Gasparyan S. , Gaug M. , Giesbrecht Paiva J.G. , Giglietto N. , Giordano F. , Gliwny P. , Godinović N. , Gradetzke T. , Grau R. , Green D. , Green J.G. , Günther P. , Hadasch D. , Hahn A. , Hassan T. , Heckmann L. , Herrera Llorente J. , Hrupec D. , Imazawa R. , Israyelyan D. , Itokawa T. , Jiménez Martínez I. , Jiménez Quiles J. , Jormanainen J. , Kankkunen S. , Kayanoki T. , Kerszberg D. , Khachatryan M. , Kluge G.W. , Kobayashi Y. , Konrad J. , Kouch P.M. , Kubo H. , Kushida J. , L'ainez M. , Lamastra A. , Lindfors E. , Lombardi S. , Longo F. , López-Coto R. , López-Moya M. , López-Oramas A. , Loporchio S. , Lorini A. , Lyard E. , Majumdar P. , Makariev M. , Maneva G. , Manganaro M. , Mangano S. , Mannheim K. , Mariotti M. , Martínez M. , Maruševec P. , Mas-Aguilar A. , Mazin D. , Menchiari S. , Mender S. , Miceli D. , Miranda J.M. , Mirzoyan R. , Molero Gonz'alez M. , Molina E. , Mondal H.A. , Moralejo A. , Nakamori T. , Nanci C. , Neustroev V. , Nickel L. , Nievas Rosillo M. , Nigro C. , Nikolić L. , Nilsson K. , Nishijima K. , Njoh Ekoume T. , Noda K. , Nozaki S. , Okumura A. , Paiano S. , Paneque D. , Paoletti R. , Paredes J.M. , Peresano M. , Persic M. , Pihet M. , Pirola G. , Podobnik F. , Prada Moroni P.G. , Prandini E. , Principe G. , Rhode W. , Ribó M. , Rico J. , Righi C. , Sahakyan N. , Saito T. , Saturni F.G. , Schmuckermaier F. , Schubert J.L. , Sciaccaluga A. , Silvestri G. , Sitarek J. , Sliusar V. , Sobczynska D. , Stamerra A. , Striškovi'c J. , Strom D. , Strzys M. , Suda Y. , Tajima H. , Takahashi M. , Takeishi R. , Temnikov P. , Terauchi K. , Terzić T. , Teshima M. , Truzzi S. , Tutone A. , Ubach S. , van Scherpenberg J. , Ventura S. , Verna G. , Viale I. , Vigliano A. , Vigorito C.F. , Vitale V. , Vovk I. , Walter R. , Wersig F. , Will M. , Yamamoto T. , Yeung P.K.H. ,   Other groups and collaborations:, Liodakis I. , Middei R. , Kiehlmann S. , Gesu L.D. , Kim D.E. , Ehlert S.R. , Saade M.L. , Kaaret P. , Maksym W.P. , Chen C.T. , De La Calle Pérez I. , Perri M. , Verrecchia F. , Domann O. , Dürr S. , Feige M. , Heidemann M. , Koppitz O. , Manhalter G. , Reinhart D. , Steineke R. , Lorey C. , McCall C. , Jermak H.E. , Steele I.A. , Fallah Ramazani V. , Otero-Santos J. , Morcuende D. , Aceituno F.J. , Casanova V. , Sota A. , Jorstad S.G. , Marscher A.P. , Pauley C. , Sasada M. , Kawabata K.S. , Uemura M. , Mizuno T. , Nakaoka T. , Akitaya H. , Myserlis I. , Gurwell M. , Keating G.K. , Rao R. , Angelakis E. , Kraus A. 
Published: February 2025

Blazars exhibit strong variability across the entire electromagnetic spectrum, including periods of high-flux states commonly known as flares. The physical mechanisms in blazar jets responsible for flares remain poorly understood to date. Our aim is to better understand the emission mechanisms during blazar flares using X-ray polarimetry and broadband observations from the archetypical TeV blazar Mrk 421, which can be studied with higher accuracy than other blazars that are dimmer and/or located farther away. We studied a flaring activity from December 2023 that was characterized from radio to very high-energy (VHE; E $>0.1$,TeV) gamma rays with MAGIC Fermi -LAT Swift XMM-Newton and several optical and radio telescopes. These observations included, for the first time for a gamma-ray flare of a blazar, simultaneous X-ray polarization measurements with IXPE, in addition to optical and radio polarimetry data. We quantify the variability and correlations among the multi-band flux and polarization measurements, and describe the varying broadband emission within a theoretical scenario constrained by the polarization data. We find substantial variability in both X-rays and VHE gamma rays throughout the campaign, with the highest VHE flux above 0.2,TeV occurring during the IXPE observing window, and exceeding twice the flux of the Crab Nebula. However, the VHE and X-ray spectra are on average softer, and the correlation between these two bands is weaker than those reported in the previous flares of Mrk 421. IXPE reveals an X-ray polarization degree significantly higher than that at radio and optical frequencies, similar to previous results for Mrk 421 and other high synchrotron peaked blazars. Differently to past observations, the X-ray polarization angle varies by ∼100^∘ on timescales of days, and the polarization degree changes by more than a factor of 4. The highest X-ray polarization degree, analyzed in 12,hr time intervals, reaches $26±2%$, around which an X-ray counter-clockwise hysteresis loop is measured with XMM-Newton . It suggests that the X-ray emission comes from particles close to the high-energy cutoff, hence possibly probing an extreme case of the Turbulent Extreme Multi-Zone model for which the chromatic trend in the polarization may be more pronounced than theoretically predicted. We model the broadband emission with a simplified stratified jet model throughout the flare. The polarization measurements imply an electron distribution in the X-ray emitting region with a very high minimum Lorentz factor (γ'_ min ≳10^4), which is expected in electron-ion plasma, as well as a variation of the emitting region size of up to a factor of 3 during the flaring activity. We find no correlation between the fluxes and the evolution of the model parameters, which indicates a stochastic nature of the underlying physical mechanism that likely explains the lack of a tight X-ray/VHE correlation during this flaring activity. Such behavior would be expected in a highly turbulent electron-ion plasma crossing a shock front.