X-Ray Polarization of the BL Lacertae Type Blazar 1ES 0229+200

We present polarization measurements in the 2–8 keV band from blazar 1ES 0229+200, the first extreme high synchrotron peaked source to be observed by the Imaging X-ray Polarimetry Explorer (IXPE). Combining two exposures separated by about two weeks, we find the degree of polarization to be ΠX = 17.9% ± 2.8% at an electric-vector position angle ψ X = 25.°0 ± 4.°6 using a spectro-polarimetric fit from joint IXPE and XMM-Newton observations. There is no evidence for the polarization degree or angle varying significantly with energy or time on both short timescales (hours) or longer timescales (days). The contemporaneous polarization degree at optical wavelengths was >7× lower, making 1ES 0229+200 the most strongly chromatic blazar yet observed. This high X-ray polarization compared to the optical provides further support that X-ray emission in high-peaked blazars originates in shock-accelerated, energy-stratified electron populations, but is in tension with many recent modeling efforts attempting to reproduce the spectral energy distribution of 1ES 0229+200, which attribute the extremely high energy synchrotron and Compton peaks to Fermi acceleration in the vicinity of strongly turbulent magnetic fields.


Introduction
The jets of active galactic nuclei (AGNs) have been shown to be luminous sources at energies spanning from the lowestfrequency radio waves to the highest-energy γ-rays.These observations clearly show that AGNs are powerful sites of nonthermal radiation originating from the acceleration of particles to highly relativistic energies.One particularly important subclass of AGNs for studying the mechanisms of particle acceleration and their subsequent radiation are blazars.Blazars are AGNs whose relativistic, highly energetic plasma jets are pointed toward our line of sight (e.g., Blandford et al. 2019;Hovatta & Lindfors 2019).Blazar emission is usually characterized by two broad spectral humps from radio to X-rays and X-rays to TeV γ-rays.The low-energy hump is interpreted as synchrotron radiation from energetic electrons, which for the most extreme blazars can peak at energies of ∼1-50 keV (Costamante et al. 2001;Di Gesu et al. 2022;Liodakis et al. 2022).
One of the most extreme blazars, at least from the standpoint of its panchromatic spectral energy distribution (SED) observed to date is 1ES 0229+200 (R.A. = 02 h 32 m 48 6, decl.= +20°1 7′17 4, z = 0.14).1ES 0229+200 is a BL Lacertae object which belongs to the rare class of extreme high synchrotron peaked (HSP) objects (Costamante et al. 2002;Biteau et al. 2020) with a synchrotron peak frequency of ν syn ∼ 10 19 Hz or ∼40 keV (Ajello et al. 2020).Due to its extreme spectrum, 1ES 0229+200 has been used to study the extragalactic background light (Aharonian et al. 2007), extragalactic magnetic fields (Tavecchio et al. 2010;Acciari et al. 2023), and Lorentz invariance violations (Tavecchio & Bonnoli 2016).The origin of the second emission hump is Compton scattering of photons by higher-energy particles, but questions about the origin of the particles responsible for the scattering and the source of the seed photons remain unanswered.Recent results from the Imaging X-ray Polarimetry Explorer (IXPE) collaboration for blazars and AGNs with synchrotron peaks well below the IXPE bandpass (e.g., Ehlert et al. 2022;Middei et al. 2023;Peirson et al. 2023) suggest relativistic electrons are the dominant scattering particles, but the results are not yet sensitive enough to distinguish between different seed photon populations.
The nature of the extreme HSP sources and the physical processes that lead to such high ν syn are still unknown.Previous attempts to model the extremely high-energy synchrotron and Compton peaks (at ∼9-20 keV and ∼12 TeV, respectively; Costamante et al. 2018;Ajello et al. 2020) in conjunction with single-zone synchrotron self-Compton (SSC) models have resulted in extreme, finely tuned parameters (e.g., Kaufmann et al. 2011).In particular, the simplest single-zone SSC models predict that 1ES 0229+200 has the following unusual properties: (1) an electron energy distribution that spans an unusually small range of energies; (2) a very high jet bulk Doppler factor of δ ∼ 40 or larger; and (3) extremely weak magnetic fields with B ∼ 30 μG to 2 mG (Kaufmann et al. 2011;Costamante et al. 2018).Models with these parameters are also dramatically out of equipartition, with electron energy densities orders of magnitude higher than the magnetic field energy density.For these reasons, more sophisticated particle acceleration models for these extreme HSP sources have been proposed details of these different proposed models can differ significantly, they all result in fits that reproduce the panchromatic SED of 1ES 0229+200 with somewhat stronger magnetic fields (B ∼ 10 −3 G).The X-ray synchrotron photons in these models originate from electrons that have recently undergone stochastic acceleration (SA) and/or diffuse shock acceleration (DSA) in the turbulent, magnetized plasma of the jet.
The IXPE satellite, launched in 2021 December (Weisskopf et al. 2022), is the first X-ray polarization mission, offering a new way to study high-energy and extreme phenomena in the Universe.It is particularly important for the study of extragalactic jets, as it can directly test particle acceleration and emission mechanisms in blazars by providing information on the geometry of the magnetic fields involved (Zhang & Böttcher 2013;Zhang et al. 2016;Liodakis et al. 2019;Peirson et al. 2022).During the first year (2022) of IXPE, several observations of Mrk 501 (Liodakis et al. 2022) and Mrk 421 (Di Gesu et al. 2022) were obtained, with more HSPs scheduled to be observed in the second year.Here we present the first X-ray polarization observations of 1ES 0229+200.All of the prior X-ray and multiwavelength polarization observations of HSPs point to a shock-accelerated, energy-stratified electron population for the origin of the synchrotron X-rays in blazar jets (Marscher & Gear 1985;Angelakis et al. 2016;Tavecchio 2021).However, this is the first time IXPE has observed an extreme HSP.In Sections 2 and 3 we present our X-ray observations and modeling, in Sections 4 and 5 the X-ray polarization results, and in Section 6 our contemporaneous radio and optical observations.In Section 7 we discuss our results and present our conclusions.Unless otherwise noted, all uncertainties described and error bars plotted correspond to 68.3% (1σ) confidence intervals for the measurement in question.

IXPE
IXPE is a NASA mission in partnership with the Italian Space Agency (ASI).As described in detail in Weisskopf et al. (2022, and references therein), the IXPE observatory includes three identical X-ray telescopes, each comprising an X-ray mirror assembly (NASA furnished) and a polarization-sensitive pixelated detector (ASI furnished), to provide imaging polarimetry over a nominal 2-8 keV band.IXPE data telemetered to ground stations in Malindi (primary) and in Singapore (secondary) are transmitted to the Mission Operations Center (MOC; at the Laboratory for Atmospheric and Space Physics, University of Colorado) and then to the Science Operations Center (SOC; at NASA Marshall Space Flight Center).Using software developed jointly by ASI and NASA, the SOC processes science and relevant engineering and ancillary data, to produce data products that are archived at the High-Energy Astrophysics Science Archive Research Center (HEASARC; at NASA Goddard Space Flight Center), for use by the international astrophysics community.IXPE targeted 1ES 0229+200 starting on 2023 January 15 with its three detector units (DUs).IXPE observations were taken in two separate segments.The first exposure took place from 2023 January 15 at 10:02 (UT) to 2023 January 18 16:44 (UT), while the second occurred from 2023 January 27 00:53 (UT) to 2023 February 01 14:44 (UT).A total of 401 ks of total exposure on source was collected, with 37% and 63% of the total exposure time taking place in the first and second time segments, respectively.Cleaned level 2 event files were computed and calibrated using standard filtering criteria with the dedicated FTOOLS tasks and latest IXPE calibration database (version 20211118).Likely background events were removed from the event lists using the selection criteria of Di Marco et al. (2023).Stokes Q and U background spectra were derived from source-free circular regions with a radius of 102″ .Extraction radii for the I Stokes spectra of 1ES 0229+200 were selected via an iterative process aimed at maximizing the signal-to-noise ratio (S/N) in the 2-8 keV energy band.This method is similar to the approach described in Piconcelli et al. (2004).We thus adopted circular regions centered on the target with radii of 42″ for DU1 and 47″ for both DU2 and DU3. 55 The same radii were also used for the Q and U Stokes spectra.We then binned the I Stokes spectra, requiring an S/N higher than 5 in each spectral channel, while a constant binning (0.2 keV) was adopted for the Q and U Stokes spectra.

XMM-Newton
XMM-Newton (Jansen et al. 2001) observed 1ES 0229+200 on 2023 January 15, quasi simultaneously with IXPE, for about 18 ks.We extracted the event lists of the European Photon Imaging Camera (EPIC-pn; Strüder et al. 2001) with the standard System Analysis Software (SAS; version 19.0.0) and the calibration database corresponding to this release.We extracted the source spectrum by selecting a region in the CCD image of 40″ radius centered on the source, and the background by extracting a source-free larger region (radius = 50″).The response matrices and auxiliary response files were generated with the SAS commands rmfgen and arfgen, respectively.Spectra were then grouped by allowing 30 counts for each spectral bin in order not to oversample the instrumental resolution by a factor larger than 3.

Swift
Before, during and after the IXPE pointing, we monitored 1ES 0229+200 with the Neil Gehrels Swift X-Ray Telescope (Swift-XRT).The Swift-XRT observations each had an exposure time of ∼1 ks and were performed in Photon Counting (PC) mode.We used the XRT Data Analysis Software (XRTDAS; 56 v. 3.6.1)to reduce, clean, and process the data.In the analysis, we used the latest calibration files available in the Swift-XRT CALDB (version 20220331).The source spectrum was extracted from the cleaned event file, adopting a circular region with a radius of 47″ .A concentric annulus with inner (outer) radii of 120 (150) arcseconds was then adopted to determine the background.The background was computed using long exposures available in the Swift archive.Finally, the spectra were binned to achieve at least 25 counts in each energy bin.

Spectral Modeling
We fit the joint XMM-Newton and IXPE Stokes I spectral data to an absorbed log-parabolic model of the form const * tbabs * logpar within XSPEC.The log-parabola model is a straightforward extension of a simple power-law model, and 55 DU1 has a slightly sharper point-spread function (PSF) than DU2 and DU3. 56https://swift.gsfc.nasa.gov/analysis/xrt_swguide_v1_2.pdfparameterized as This model is commonly fit to X-ray spectra of blazars that show evidence of curvature beyond a single power law.We find best-fit parameters of Γ = 1.82 ± 0.01 and β = 0.13 ± 0.02, indicating that the spectrum is steepening with increasing energy.The Galactic absorption column density was frozen at N H = 7.81 × 10 20 cm −2 , as determined by the HI4PI survey (HI4PI Collaboration et al. 2016) for our fiducial model, and the pivot energy was fixed to E piv = 1 keV.Allowing the column density to be freely fit by the model results in a best-fit value of N H = 8.01 ± 0.01 × 10 20 cm −2 and a negligible improvement to the overall fit.
The overall χ 2 value of this fit is χ 2 = 409 with 370 degrees of freedom.We compare these values to a simple power-law model, which has a best-fit photon index value of Γ = 1.894 ± 0.006 with χ 2 = 465 and 371 degrees of freedom.When E piv = 1 keV, the log-parabola model is a simple extension of a power-law model, and we can use an F-test to compare the statistical significance of the additional parameter to the fit improvement.Under a null hypothesis where the true model is a power law, the probability of such an improvement is P null ∼ 6 × 10 −12 .We can also use the Akaike and Bayesian information criteria (AIC and BIC, respectively) to compare their overall fit quality.For this total change in the fit statistic (Δχ 2 = −56.12),both criteria strongly favor the log-parabolic model over the simple power-law model after accounting for the additional parameters (ΔAIC = −52 and ΔBIC = −52).We therefore consider this log-parabolic model as our fiducial spectrum for the remainder of this work.The const term, which accounts for cross-calibration terms in the effective areas of the four different detectors, was fixed to unity for IXPE DU1.The best-fit constant normalization offsets for DU2, DU3, and our XMM-Newton spectra are 0.96 ± 0.02, 0.89 ± 0.01, and 0.92 ± 0.01, respectively.These factors are consistent with previous results reported by spectral fits using IXPE (Ehlert et al. 2022).The spectra from all four detectors (IXPE DU1, IXPE DU2, IXPE DU3, and XMM-Newton), along with the best-fit log-parabolic model, are displayed in Figure 1.
Although the log-parabola model deployed here represents a "good" fit to the data, the best-fit parameters differ significantly from the results of Costamante et al. (2018), who used an identical model to fit Swift and NuSTAR observations of 1ES 0229+200.In their best-fit model,57 Γ = 1.49± 0.04 and β = 0.27 ± 0.02.These different values suggest that 1ES 0229 +200 during the IXPE observations had a much softer photon index at 1 keV, but did not further soften as rapidly with energy as it did in 2013 when the NuSTAR observations were performed.During our observations, the X-ray flux is at the 1 × 10 −11 erg cm −2 s −1 level in the 2-10 keV band, whereas during the Costamante et al. (2018) observations it is about a factor of two brighter (1.95 × 10 −11 erg cm −2 s −1 ).1ES0229 +200 becomes harder when brighter as was found in Acciari et al. (2020), hence the difference in spectral shape could potentially arise from the difference in flux.However, comparing these results directly is subject to the caveat that our analysis considers a much lower energy band than the result of Costamante et al. (2018).
We present best-fit log-parabolic models for each of the individual Swift observations of 1ES 0229+200 along with their fluxes in the 2-10 keV band in  second segment are ∼2× larger than for the first segment.We note that there is a small mismatch between the spectral parameters derived from the Swift+IXPE and the XMM-Newton+IXPE fits.This could be due to the observations not being strictly simultaneous, as well as the fact that XMM-Newton's significantly higher effective area above ∼5 keV can better constrain the shape of the spectrum in the ∼5-10 keV band.Nevertheless, the parameters are within 1σ and the polarimetric results (see Section 4) are not sensitive to small variations of the spectral shape.

Polarization Measurements
We have determined the broadband polarization of 1ES 0229 +200 using two different analysis methods: one by measuring the average normalized Stokes parameters over various energy bands to determine the average polarization degree and angle using IXPEOBSSIM (Pesce-Rollins et al. 2019; Baldini et al. 2022), and the other by a simultaneous spectro-polarimetric fit of the Stokes Q and U spectra along with the XMM-Newton and IXPE Stokes I spectra.We discuss the results of these two analysis methods separately, since they measure slightly different, albeit related, polarization quantities.

Polarization Cube
Our polarization analysis determines the average polarization using the statistical framework of Kislat et al. (2015) without any event-specific weights.Over the entire nominal IXPE bandpass of 2-8 keV with 33,502 net counts, the average values of the Stokes parameters are Q = 0.063 ± 0.035 and U = 0.110 ± 0.035.Under the null hypothesis of zero true polarization, the χ 2 for these Stokes parameter values is 12.80 with two degrees of freedom.This value corresponds to a confidence level of 99.8%, indicating detection.The corresponding polarization degree of this measurement is Π X = 12.5% ± 3.2% and the electric-vector polarization angle is ψ = 30°.0± 8°.0, east of north.The significance of the detection depends strongly on the choice of energy band.At lower energies of 2-4 keV with 29,561 net counts, the normalized Stokes parameters differ from 0 at the 99.94% confidence level.On the other hand, the 3941 net counts in the 4-8 keV band show no statistically significant evidence of polarization at 99% confidence.The 99% confidence upper limit in the 4-8 keV band is Π X < 29.8%.Although no statistically significant polarization is detected in the 4-8 keV band, the upper limit is consistent with the polarization degree observed at lower energies.The significance of the detection is highest in the 2-6 keV energy band, for which the significance of the detection is securely above 99.99% confidence: Π X = 14.9% ± 3.0% and ψ = 33°.9± 5°.7.
The polarization behavior of this source is consistent between the two observation segments.Restricting the data to only include events gathered during each of the individual segments gives nearly identical results (see Table 1).We therefore conclude that we can combine the results from both  time segments without loss of any polarization information.
The Stokes parameters for each segment, as well as the timeintegrated averages, are shown in Figure 3.

Spectro-polarimetric Fits
To investigate further the extent to which we can detect and measure polarization in 1ES 0229+200, we have added to our spectral fit the Q and U spectra for all three IXPE DUs.We have performed the spectro-polarimetric model fit using XSPEC (Strohmayer 2017) by including an energy-independent polarization model component to our fiducial log-parabolic model.Unlike the polarization cube analysis, the spectra used for these fits were weighted using the method of Di Marco et al. (2022).In XSPEC terms, this corresponds to a model of the form const * tbabs * polconst * logpar.Our best-fit model for the polarization degree and angle from this model gives Π X = 17.9% ± 2.8% and ψ = 25°.0± 4°. 6.The total χ 2 of this fit is 586.56 with 542 degrees of freedom.
The spectro-polarimetric fit enables another hypothesis test for the presence of polarization.Assuming a polarization degree of Π = 0 and fixing = 0 results in a fit with χ 2 = 626.36 with 544 degrees of freedom.The probability of obtaining a χ 2 value equal to or exceeding this under the assumption of zero polarization is P null = 0.0028, providing strong evidence that these data are an improbable realization of the model assuming zero polarization.Furthermore, allowing the polarization degree and angle to be free parameters provides a statistically significant improvement to the overall fit (Δχ 2 = −39.8with two fewer degrees of freedom), which for these data correspond to a null hypothesis improvement (as determined by an F-test) of P null = 1.9 × 10 −8 .The corresponding changes in the AIC and BIC are ΔAIC = −36 and ΔBIC = −27, respectively.The mean Stokes Q and U spectra, along with the best-fit spectro-polarimetric model, are shown in Figure 4.

Energy-dependent Polarization
The nondetection of polarization in the 4-8 keV band, as compared to detection when the energy range is extended down to 2 keV, leads to the question of whether this is caused by insufficient photon statistics in the 4-8 keV band or by an energy-dependent polarization.While it is clear that the vast majority of the photons are observed at lower energies, we further test the null hypothesis of constant polarization as a function of energy using our spectro-polarimetric model.We replace the constant polarization model component with a constant + linear energy dependence of the polarization degree and angle (const * tbabs * pollin * logpar in XSPEC).We find that the 90% confidence interval for the polarization degree's linear term is consistent with zero, and that the total improvement in χ 2 with respect to the constant polarization model is Δχ 2 = −1.1 for two fewer degrees of freedom, entirely consistent with the expected improvement to the fit arising from arbitrary additional parameters.Both of these calculations indicate that the evidence for any dependence of the polarization on energy is marginal.We therefore conclude that a constant polarization degree across the entire 2-8 keV band is consistent with the observations.

Time-dependent Variations of the Stokes Parameters
Recent observations of other blazars, in particular Mrk 421 (Di Gesu et al. 2023) show clear evidence of intraobservation variability.In the case of Mrk 421, this variability is consistent with the electric-vector polarization angle rotating at a constant rate of ∼80 degrees day −1 .We test for variability in the Stokes parameters of 1ES 0229+200 using a χ 2 -based hypothesis test.For this test, the null hypothesis is the assumption that the true Stokes parameters are equal to the mean value in each of N = 10 time bins.These time bins are explicitly assigned to result in four time bins associated with the first segment and six with the second.For this test, we use polarization cubes in the 2-6 keV band in order to maximize the S/N of the polarization measurement in each time bin.Comparing the as-measured Stokes parameters in each time bin with the time-integrated mean, we find that the total χ 2 of this model is χ 2 = 14.30 with 18 degrees of freedom.Similarly "good" values of χ 2 are obtained for all values of N in the range of 8-20, suggesting that the lack of variability is not an artifact of our choice of the number of time bins.We therefore conclude there is insufficient evidence to reject the null hypothesis that the Stokes parameters in all of the time bins are statistically consistent with the time-integrated average.The variations of the Stokes parameters with time are visualized in Figure 5.

Radio and Optical Observations
During the IXPE observation we coordinated a radio, millimeter, and optical campaign with the Effelsberg 100 m radio telescope, the SubMillimeter Array (SMA), the Nordic Optical Telescope (NOT), the Perkins telescope, and the Observatorio de Sierra Nevada (OSN).Observations with the Effelsberg and SMA were obtained at 6 cm (4.85 GHz) and 1ES 0229+200 has a bright host galaxy, the starlight from which dominates the emission at optical wavelengths.This results in the source appearing to be less variable and less polarized than is the case for the active nucleus alone.To correct for the host galaxy contribution and estimate the intrinsic polarization degree ( intr P ) of the source, we have performed detailed modeling of the host galaxy flux distribution following the method of Nilsson et al. (2007).We estimate the contribution of the host galaxy to the R-band emission (I host ) to be 0.54 mJy for an aperture with a 5″ radius and 0.67 mJy for 7 5.We then use the estimated I host for the apertures of the respective telescopes to correct the observed polarization degree Π obs .This is achieved by subtracting the host contribution from the total flux density (I),

Discussion and Conclusions
The IXPE observations of 1ES 0229+200 have detected X-ray polarization of Π X ∼ 17.9% along ψ X ∼ 25°.The polarization degree we measure for this blazar is significantly higher than that measured during the first IXPE observation of Mrk 501 (∼10%; Liodakis et al. 2022) and similar to that measured during the first IXPE observation of Mrk 421 (∼15% ;Di Gesu et al. 2022).As for these two HSP blazars, the X-ray polarization degree of 1ES 0229+200 is significantly higher than the optical values of ∼2%.In this case, the ratio of X-ray to optical Π is >7, making it the most strongly chromatic source in polarization observed thus far.
1ES 0229+200 is not the first HSP blazar with a higher degree of X-ray polarization than at longer wavelengths, yet with similar (within a few degrees) polarization angles in X-ray, optical, and radio bands.As discussed at length in Di Gesu et al. (2022) and Liodakis et al. (2022), the increasing polarization degree with energy suggests that the electrons generating these photons are accelerated at a shock.In this model, the X-ray photons are generated by synchrotron radiation from electrons immediately downstream of the shock front, while lower-energy photons, because of their longer radiative cooling time, are generated by electrons further downstream where the magnetic fields are more turbulent (Marscher & Gear 1985).
There remains one crucial difference between the results of 1ES 0229+200 and those of Mrk 501 discussed in Liodakis et al. (2022), however.In Mrk 501 the X-ray polarization angle was parallel to the position angle of the radio jet.In 1ES 0229 +200 the radio jet has an apparent position angle of ∼160°east of north (Piner & Edwards 2018) while the X-ray polarization angle is ∼30°, hence there is no obvious relationship between the polarization angle and the radio jet position angle.A similar mismatch was observed in the first observation of Mrk 421 (Di Gesu et al. 2022), although later observations have also shown clear evidence of time variability in Mrk 421ʼs polarization angle (Di Gesu et al. 2023).The agreement between the X-ray polarization angle and the radio jet position angle in Mrk 501 was considered an important clue suggesting the presence of energy-stratified shock acceleration in the jet.It remains unclear how to reconcile the disagreement between the X-ray polarization angle and the radio jet direction for Mrk 421 and 1ES 0229+200 with the predictions of an energy-stratified shock when all three blazars show the same wavelengthdependent polarization degree.However, we note that the estimates for the position angle of the jet are not contemporaneous.There is now a wealth of evidence showing that jets can change their position angle over time, such as in NRAO 150 (Agudo et al. 2007), OJ 287 (Britzen et al. 2018), PG1553+113 (Lico et al. 2020), and others (Lister et al. 2013;Weaver et al. 2022).In some cases, much larger position angle variations than the 50°mismatch we observe have been seen over a few years (Lister et al. 2013;Weaver et al. 2022).Further monitoring of both the X-ray polarization and jet structure is needed to determine the extent to which such mismatches might be attributed to variability of the jet and polarization directions.The optical polarization angle of ∼0°also appears to have no obvious relationship to either the X-ray polarization angle or the radio jet, although the polarization angles are closer to each other than to the radio jet's position angle.Fully accounting for the apparent discrepancy between the two polarization angles is beyond the scope of this paper, but does provide additional evidence that the regions where the optical and X-ray photons are generated appear to be largely disconnected from one another.
The reasonably high polarization degree, consistent with less extreme HSP blazars such as Mrk 501 and Mrk 421, further complicates any effort to account for the extreme properties of 1ES 0229+200.Recent single-zone (e.g., Tavecchio et al. 2022) and two-zone models (e.g., Aguilar-Ruiz et al. 2022) often predict the electrons responsible for the X-rays originate in shock acceleration in highly turbulent regions of the magnetized jet plasma, where the magnetic fields are not expected to have any coherent direction.The relatively high polarization degree we observe is therefore in tension with such models.A possible way to reconcile these apparently disparate observations is to identify different length scales for the shock acceleration and synchrotron radiation-the magnetic fields may be turbulent on the small scales where electron acceleration occurs but more ordered and structured on the larger scales where these same highenergy electrons are emitting their X-rays.Other proposed models such as the multiple shock model proposed in, e.g., Zech & Lemoine (2021) do not explicitly require turbulence, but further testing and simulation work will be required to determine if such a model is able to reproduce the polarization results presented in this work.Magnetic reconnection (e.g., Matthews et al. 2020, and references therein) appears to be disfavored as a viable particle acceleration model based on the low magnetic field strengths estimated from simple SSC model fits in Kaufmann et al. (2011) and Costamante et al. (2018), but given the unusual fit parameters from these models there is reason to question whether or not this estimated magnetic field strength is an accurate measurement of what is physically present in the jet.Although it is beyond the scope of this paper to develop a possible particle acceleration model that fully explains the polarization signal and multiwavelength properties of 1ES 0229+200, the fact that it has a similar X-ray polarization degree to Mrk 501 and Mrk 421 (despite the very different properties of these three blazars at TeV energies) is new evidence that informs future theoretical work understanding the most extreme blazar jets.Single-zone acceleration models are more favored for less extreme HSP blazars due to correlations between the X-ray and γ-ray light curves (e.g., Katarzyński et al. 2005, and references therein), but such a hypothesis has not been fully tested with 1ES 0229+200.The main reason this test has not yet been performed is that the light curves of 1ES 0229+200 (in particular the γ-ray light curve) show less variability than the less extreme blazars.Further observations of 1ES 0229 +200 with X-ray and γ-ray telescopes58 including IXPE will help identify the extent to which this source is variable in either flux or spectral properties.We find different spectral parameters than previous observations (e.g., Costamante et al. 2018) despite no obvious changes in its detected X-ray flux.It remains clear that the IXPE results for this extremely HSP blazar will require further modeling efforts to reconcile the extreme photon energies with the IXPE measurements fully.Simultaneous IXPE and TeV observations may help further elucidate still unanswered questions about particle acceleration within the jet of this AGN.
(e.g., Zech & Lemoine 2021; Aguilar-Ruiz et al. 2022; Tavecchio et al. 2022).Although the Original content from this work may be used under the terms of the Creative Commons Attribution 4.0 licence.Any further distribution of this work must maintain attribution to the author(s) and the title of the work, journal citation and DOI.
Figure 2.This figure also includes the time-resolved IXPE light curve for each observation segment.Each of these spectral fits was performed over the 0.3-10 keV energy band.Both the Swift and IXPE data show evidence that 1ES 0229+200 became brighter during the second IXPE observation segment.We have fit the same logparabolic model to each IXPE observation segment to test for changes in the spectral parameters.No changes in the spectral parameters between the two segments were identified beyond the overall normalization of the model.The values of Γ for the first and second segments are Γ = 1.69 ± 0.28 and Γ = 1.68 ± 0.19, respectively.The best-fit values of β for the first and second segments are β = 0.41 ± 0.28 and β = 0.55 ± 0.19, respectively.The normalizations for the

Figure 1 .
Figure 1.The best-fit spectral model to the XMM-Newton and IXPE Stokes I data for 1ES 0229+200.Top: the spectra with the best-fit log-parabolic model overlaid, corresponding to Γ = 1.82 ± 0.01 and β = 0.13 ± 0.02.Bottom: the normalized residuals of the data with respect to the best-fit model.

Figure 2 .
Figure 2. The best-fit spectral model parameters and flux of 1ES 0229+200 as determined by Swift observations taken before and during the IXPE observations.The shaded gray regions correspond to the two segments where IXPE was observing 1ES 0229+200.Top panel: the best-fit values for Γ, the photon index at 1 keV, as a function of time.Upper-middle panel: the best-fit values for β, the curvature parameter, as a function of time.Lower-middle panel: the Swift 2-10 keV flux as a function of time, in units of 10 −11 erg cm −2 s −1 .We have also included the corresponding flux measurement from the XMM-Newton observation in blue, which is consistent with the measurements from Swift.Bottom panel: the IXPE count rate during the two segments.The parameters for each of the fits to the Swift data are consistent with each other, but slightly softer than the measurements of Costamante et al. (2018) when 1ES 0229+200 was a factor of two brighter.

Figure 3 .
Figure3.Time-averaged Stokes parameters in the 2-8 keV band for the two observation segments of 1ES 0229+200.See the main text for the exposure times and date ranges for each segment.It is clear that the mean Stokes parameters from the first segment (in red) are consistent with those from the second segment (in blue).We can therefore safely combine the events from both segments, the results of which are shown in black.For reference, silver circles with constant polarization degrees of Π = 5% and Π = 15% are also drawn.
et al. (2016).Only the R-band Π O estimates have been corrected.All of the multiwavelength observations during the IXPE pointings are displayed in Figure6.1ES 0229+200 is faint in radio (0.06 Jy at 4.85 GHz and 0.01 Jy at 225.5 GHz), which prevents us from detecting polarization at the 3σ level.A QUIVER observation during the first segment of the IXPE observation yields an upper limit of 7% (99% confidence).Similarly, two separate SMAPOL observations after the second IXPE segment (MJDs 59,980 and 59,981) yield <22% (99% confidence) and <7% (99% confidence), respectively.For the NOT observations during the first IXPE segment, Π O = 2.42% ± 0.72% along ψ O = −2°.4± 8°.5.We were unable to obtain contemporaneous optical observations during the second IXPE measurement.However, a few days later, the Perkins telescope measured Π O = 3.2% ± 0.7% along a ψ O = −5°.1 ± 8°.7, consistent with the NOT results.This suggests similar levels of optical polarization for both IXPE observations.

Figure 4 .
Figure 4. Top: Stokes Q and U spectra as measured by IXPE.For presentation purposes only, we have grouped the spectra from all three IXPE detectors together and shifted the Stokes U spectra by 0.05 keV.The solid red and blue curves correspond to the best-fit models for the Q and U spectra, respectively.Bottom: the residuals between the spectra and the best-fit model as a function of energy.

Figure 5 .
Figure 5. Variations of the Stokes parameters as a function of time during the observations of 1ES 0229+200.The circles correspond to the mean Stokes parameters and their statistical uncertainties in the 2-6 keV band, with their colors corresponding to the time in kiloseconds since the beginning of the first observation segment.Both segments are included.The white cross denotes the time-averaged mean Stokes parameters and their uncertainties.As we describe in detail in the text, we are unable to reject the null hypothesis that the true Stokes parameters in each time bin are equal to the time-averaged values.

Figure 6 .
Figure 6.Multiwavelength optical observations of 1ES 0229+200.The panels show the optical brightness (top), polarization degree (middle), and polarization angle (bottom).The IXPE observing periods are marked with the gray shaded regions.A correction for the host galaxy contribution has been applied exclusively to the Rband magnitudes and polarization measurements.

Table 1
Normalized Stokes Parameters in the 2-8 keV Bandpass for the Two IXPE Observation Segments of 1ES 0229+200, as Determined by the Polarization