We report the first observations of a long-duration very-high-energy (VHE; E > 100 GeV) flare from BL Lacertae (VER J2202+422), taken with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). On October 15, 2022, the Fermi-Large Area Telescope (LAT) detected elevated GeV activity originating from this blazar. This triggered a multiwavelength campaign, which includes observations from VERITAS, Swift, NuSTAR, and select optical and radio observatories. VERITAS observed the source for a total of ∼ 9.8 hours between September 1, 2022 and December 1, 2022. An analysis of these data yields a ∼ 28 σ detection of the source. While previously observed VHE flares from BL Lacertae have lasted on time-scales of minutes to days, VERITAS continued to detect flaring activity from the source for over a month (∼ 40 days) after the original flaring activity was detected with Fermi-LAT. Broadband spectral modeling shows that a synchrotron self-Compton (SSC) model with an external inverse-Compton (EC) component is preferred over a one-zone SSC model.

V1674 Her (Nova Her 2021) is the fastest classical nova ever recorded, with an optical decline time of t(2) similar to 1 day, typically interpreted as evidence for a white dwarf mass close to the Chandrasekhar limit. We present a broadband X-ray study of V1674 Her combining contemporaneous XMM-Newton and NuSTAR observations in quiescence to directly constrain the white dwarf's mass and magnetic field strength. The hard X-ray emission is modeled using a physically motivated postshock accretion column model that accounts for the temperature gradient in the flow and reflection from the white dwarf's surface. Under the assumption that the accretion disk is truncated at the corotation radius, we obtain a white dwarf mass of M = 1.09( -0.06 )(+ 0.07) M-circle dot . An independent constraint derived from timing analysis of the X-ray power spectrum yields a consistent value of M = 1.12 +/- 0.06 M-circle dot. These values are significantly lower than those inferred from empirical decline-time relations, suggesting that such relations may overestimate white dwarf masses in extreme fast novae. From the inferred accretion rate and magnetospheric radius, we estimate a surface magnetic field strength of B = 21.3 (- 5.7) (+ 6.6) ( stat )( -8.1) (+ 12.9) (sys ) MG , placing V1674 Her at the high end of the magnetic field distribution for intermediate polars. Our results demonstrate that even the fastest novae do not necessarily host near-Chandrasekhar white dwarfs, highlighting the importance of direct X-ray constraints and suggesting that additional parameters beyond white dwarf mass play a key role in setting nova timescales.

The origin of intergalactic magnetic fields (IGMFs) remains one of the key open questions in cosmology. Gamma-ray pair halos produced by electromagnetic cascades from TeV-emitting blazars provide a powerful indirect probe of these fields. In this work, we present a novel search for pair halos that explicitly exploits their expected anisotropic morphology, aligning with the projected orientation of blazar jets on the sky. Using a Monte Carlo framework to model the spatial distribution of cascade emission, we identify an optimal sample of 21 high-synchrotron-peaked BL Lac objects with well-constrained jet position angles from radio interferometry. By rotating and stacking Fermi-LAT observations of these sources along their jet directions, we enhance sensitivity to anisotropic extended emission that would be diluted in traditional orientation-agnostic analyses. Applying a likelihood analysis to the combined dataset, we find evidence for a non-zero IGMF, excluding the null hypothesis at 3.8σ level and obtaining a best-fit field strength of B_0 = 2.8 × 10^-16 G, with a 99% confidence interval of 0.9 × 10^-16 G < B_0 < 8.9 × 10^-16 G. Our result is consistent with previous constraints from spectral, spatial, and temporal studies, while demonstrating that incorporating anisotropic information provides a significant gain in sensitivity. This approach opens a new avenue for probing intergalactic magnetism and highlights the potential of future high-angular-resolution gamma-ray observations to directly image pair halos and map magnetic fields in cosmic voids.

Identifying the physical mechanism driving blazar flares remains a central challenge in high-energy astrophysics. We show that the energy dependence of the standard deviation of the polarization angle variability (sigma PA) provides a powerful and robust discriminator of blazar flaring mechanisms. Using particle-in-cell-integrated polarized radiative transfer simulations, we perform the most rigorous statistical analyses of polarization variability to date. We demonstrate that magnetic reconnection and magnetized turbulence imprint qualitatively distinct energy dependence of sigma PA that directly reflects their different magnetic field evolution and particle transport. Reconnection predicts higher sigma PA with higher photon energy until the synchrotron spectral peak, whereas turbulence produces nearly flat sigma PA across the synchrotron spectral component. These trends are resilient to realistic observational limitations. Applying our results to optical and IXPE data of Mrk 421 and 1ES 1959+650, we find strong evidence for reconnection-driven flares embedded in a turbulent blazar zone. Energy-dependent sigma PA emerges as a decisive new probe of particle acceleration in relativistic jets.

The polarization of X-ray synchrotron emission in blazars directly probes the magnetic field geometry and particle acceleration processes in relativistic jets. We use particle-in-cell simulations of magnetic reconnection and magnetized turbulence, coupled to polarization-sensitive radiative transfer code, to interpret Imaging X-ray Polarimetry Explorer (IXPE) observations of Mrk 421 during a high flux state recorded in December of 2023. To evaluate the fitness of the two theoretical scenarios, we rely on a quantitative comparison of the statistical properties of simulated and observed X-ray flux and polarization light curves using five evaluation metrics, rather than attempting to fit individual data points. We propose a turbulence-driven multi-zone model where jet emission is represented as the sum of the radiative output of N independent cells, each described by a particle-in-cell simulation. Comparison of ensembles of simulated Stokes-parameter light curves with IXPE data shows that magnetic-reconnection-dominated models provide the best match to the observed X-ray flux and polarization dynamics. The optimal configuration corresponds to N = 15 emitting cells, which reproduces the observed amplitudes and timescales of the X-ray flux and polarization variations. Magnetized turbulence models underpredict both the flux and polarization variability. Our results indicate that a multi-zone, reconnection-powered emission scenario can describe the X-ray polarization behavior of Mrk 421 and establish a quantitative framework for testing theoretical models against IXPE observations of other high-synchrotron-peaked blazars.

While the sources of the diffuse astrophysical neutrino flux detected by the IceCube Neutrino Observatory are still largely unknown, one of the promising methods to improve our understanding of them is investigating the potential temporal and spatial correlations between neutrino alerts and the electromagnetic radiation from blazars. We report on the multiwavelength target-of-opportunity observations of the blazar B3 2247+381, taken in response to an IceCube multiplet alert for a cluster of muon neutrino events compatible with the source location between 2022 May 20 and 2022 November 10. B3 2247+381 was not detected with VERITAS during this time period. The source was found to be in a low-flux state in the optical, ultraviolet, and gamma-ray bands for the time interval corresponding to the neutrino event, but was detected in the hard X-ray band with NuSTAR during this period. We find the multiwavelength spectral energy distribution is described well using a simple one-zone leptonic synchrotron self-Compton radiation model. Moreover, assuming the neutrinos originate from hadronic processes within the jet, the neutrino flux would be accompanied by a photon flux from the cascade emission, and the integrated photon flux required in such a case would significantly exceed the total multiwavelength fluxes and the VERITAS upper limits presented here. The lack of flaring activity observed with VERITAS, combined with the low multiwavelength flux levels, as well as the significance of the neutrino excess being at a 3 σ level (uncorrected for trials), makes B3 2247+381 an unlikely source of the IceCube multiplet. We conclude that the neutrino excess is likely a background fluctuation.

We report the Imaging X-ray Polarimetry Explorer (IXPE) polarimetric and simultaneous multiwavelength observations of the high-energy-peaked BL Lacertae object (HBL) 1ES 1959+650, performed in 2022 October and 2023 August. In 2022 October, IXPE measured an average polarization degree Π _X = 9.4% ± 1.6% and an electric-vector position angle ψ _X = 53° ± 5°. The polarized X-ray emission can be decomposed into a constant component, plus a rotating component, with the rotation velocity ω _EVPA = (−117 ± 12) deg day ^−1 . In 2023 August, during a period of pronounced activity of the source, IXPE measured an average Π _X = 12.4% ± 0.7% and ψ _X = 20° ± 2°, with evidence (∼0.4% chance probability) for a rapidly rotating component with ω _EVPA = 1864 ± 34 deg day ^−1 . These findings suggest the presence of a helical magnetic field in the jet of 1ES 1959+650 or stochastic processes governing the field in turbulent plasma. Our multiwavelength campaigns from radio to X-ray reveal variability in both polarization and flux from optical to X-rays. We interpret the results in terms of a relatively slowly varying component dominating the radio and optical emission, while rapidly variable polarized components dominate the X-ray and provide minor contribution at optical wavelengths. The radio and optical data indicate that on parsec scales the magnetic field is primarily orthogonal to the jet direction. On the contrary, X-ray measurements show a magnetic field almost aligned with the parsec jet direction. Confronting with other IXPE observations, we guess that the magnetic field of HBLs on subparsec scale should be rather unstable, often changing its direction with respect to the Very Long Baseline Array jet.

The X-ray polarization observations, made possible with the Imaging X-ray Polarimetry Explorer (IXPE), offer new ways of probing high-energy emission processes in astrophysical jets from blazars. Here, we report the first X-ray polarization observation of the blazar S4 0954 +65 in a high optical and X-ray state. During our multi-wavelength (MWL) campaign of the source, we detected an optical flare whose peak coincided with the peak of an X-ray flare. This optical-X-ray flare most likely took place in a feature moving along the parsec-scale jet, imaged at 43 GHz by the Very Long Baseline Array (VLBA). The 43 GHz polarization angle of the moving component underwent a rotation near the time of the flare. In the optical band, prior to the IXPE observation, we measured the polarization angle to be aligned with the jet axis. In contrast, during the optical flare, the optical polarization angle was perpendicular to the jet axis; after the flare, it reverted to being parallel to the jet axis. Due to the smooth behavior of the optical polarization angle during the flare, we favor shocks as the main acceleration mechanism. We also infer that the ambient magnetic field lines in the jet were parallel to the jet position angle. The average degree of optical polarization during the IXPE observation was (14.3 +/- 4.1)%. Despite the flare, we only detected an upper limit of 14% (at 3 sigma level) on the X-ray polarization degree; however, a reasonable assumption on the X-ray polarization angle results in an upper limit of 8.8% (3 sigma). We modeled the spectral energy distribution (SED) and spectral polarization distribution (SPD) of S4 0954 +65 with leptonic (synchrotron self-Compton) and hadronic (proton and pair synchrotron) models. Our combined MWL polarization observations and SED modeling tentatively disfavor the use of hadronic models for the X-ray emission in S4 0954 +65.

The Moon is the closest celestial gamma-ray emitting object. Its gamma-ray emission arises from interactions between Galactic cosmic rays (CRs) and the lunar surface. While the lunar GeV gamma-ray spectrum is dominated by a continuum from hadronic decay processes, the MeV emission exhibits both continuum and distinctive spectral lines from nuclear de-excitation and radioactive decay processes. Using Geant4 Monte Carlo particle simulations, we model the lunar gamma-ray spectrum. Our results demonstrate its consistency with Fermi-Large Area Telescope observations, and predict that next-generation MeV gamma-ray instruments will detect both the lunar MeV continuum and several key spectral line features, notably the 1.779 MeV line from ^28 Si de-excitation enhanced by the lunar surface composition, the e ^+ e ^− annihilation line, and radioactive decay lines from ^22 Na ( τ ≈ 3.75 yr) and long-lived ^26 Al ( τ ≈ 1 Myr). These gamma-ray lines are sensitive to CRs with energies ≲1 GeV nuc ^−1 , offering unique temporal probes of CR activity over different timescales. Observations of the lunar MeV gamma-ray spectrum will therefore open a new window to study the current irradiation of the solar-terrestrial environment by low-energy CRs and its long-term temporal evolution.

In 2023 November, the Fermi Large Area Telescope detected a γ -ray flare from the high-redshift blazar GB6 B1428+4217 ( z = 4.715). We initiated a multiwavelength follow-up campaign involving Swift, NuSTAR, the Sierra Nevada and Perkins Observatories, and the Effelsberg 100 m radio telescope. This source, also known as 5BZQ J1430+4204, has shown an anomalous soft X-ray spectrum in previous observations, including possible ionized absorption features or signatures of bulk Comptonization of thermal electrons, which are also detected during the flaring episode. Simultaneous optical data revealed a polarization fraction of ∼8% in the R band, confirming that synchrotron emission dominated over thermal emission from the accretion disk. The hard X-ray flux was enhanced during the flare. Modeling of the broadband spectral energy distribution suggests that the high-energy component is dominated by Compton scattering by external seed photons from the accretion disk. The origin of the flare is consistent with the injection of a hard-spectrum electron population in the emission region. With a γ -ray luminosity among the top 5% of flaring events, GB6 B1428+4217 exemplifies a prototypical MeV blazar. Its Compton-dominated spectral energy distribution and extreme luminosity are in line with expectations from the blazar sequence. High-redshift flares like this are critical for understanding jet physics in the early Universe and may improve detection prospects with future missions such as the Compton Spectrometer and Imager.

Assuming Galactic cosmic rays originate in supernovae and the winds of massive stars, starburst galaxies should produce very-high-energy (VHE; E > 100 GeV) gamma-ray emission via the interaction of their copious quantities of cosmic rays with the large reservoirs of dense gas within the galaxies. Such VHE emission was detected by VERITAS from the starburst galaxy M82 in 2008–09. An extensive, multiyear campaign followed these initial observations, yielding a total of 254 hr of good-quality VERITAS data on M82. Leveraging modern analysis techniques and the larger exposure, these VERITAS data show a more statistically significant VHE signal (∼6.5 standard deviations, σ ). The corresponding photon spectrum is well fit by a power law (Γ = 2.3 ± 0.3 _stat ± 0.2 _sys ), and the observed integral flux is F (>450 GeV) = (3.2 ± 0.6 _stat ± 0.6 _sys ) × 10 ^−13 cm ^−2 s ^−1 , or ∼0.4% of the Crab Nebula flux above the same energy threshold. The improved VERITAS measurements, when combined with various multiwavelength data, enable modeling of the underlying emission and transport processes. A purely leptonic scenario is found to be a poor representation of the gamma-ray spectral energy distribution (SED). A lepto-hadronic scenario with cosmic rays following a power-law spectrum in momentum (index s ≃ 2.25) and with significant bremsstrahlung below 1 GeV provides a good match to the observed SED. The synchrotron emission from the secondary electrons indicates that efficient nonradiative losses of cosmic-ray electrons may be related to advective escape from the starburst core.

TeV-emitting BL Lac–type blazars represent the extreme end of the blazar population. They are characterized by relatively weak jets and radiatively inefficient accretion disks. Particles accelerated in these jets experience fewer radiative losses, allowing them to reach energies beyond the TeV scale and produce TeV gamma-ray emission. The study of TeV blazars is constrained by the limited number of known sources in this category. Currently, only 56 high-synchrotron-peaked BL Lacs have been detected at energies above 0.1 TeV. Searches for TeV emission from BL Lacs typically target sources with bright X-ray emission and a synchrotron peak at or above 1 keV. The recently released eRASS1 catalog by the eROSITA collaboration, which covers half of the sky, represents the deepest X-ray survey in the soft X-ray band to date. Utilizing the eROSITA survey, combined with infrared data from the Wide-field Infrared Survey Explorer and archival radio observations, we have identified 121 TeV-emitting blazar candidates. Our search introduces selection criteria based on the radio to infrared that remove quasar-like objects that have similar infrared spectra and X-ray fluxes as TeV-emitting BL Lacs. In our search, we find 23 objects that had not been detected in the ROSAT X-ray survey and 11 that have not been previously associated with blazars. The candidates resulting from our search are suitable for follow-up observations with currently operating imaging atmospheric Cherenkov telescopes, as well as future facilities like the CTAO Observatory.

Observations of linear polarization in the 2–8 keV energy range with the Imaging X-ray Polarimetry Explorer (IXPE) explore the magnetic field geometry and dynamics of the regions generating nonthermal radiation in relativistic jets of blazars. These jets, particularly in blazars whose spectral energy distribution peaks at X-ray energies, emit X-rays via synchrotron radiation from high-energy particles within the jet. IXPE observations of the X-ray-selected BL Lac–type blazar 1ES 1959+650 on 2022 May 3–4 showed a significant linear polarization degree of Π x = 8.0% ± 2.3% at an electric-vector position angle ψ x = 123° ± 8°. However, on 2022 June 9–12, only an upper limit of Π x ≤ 5.1% could be derived (at the 99% confidence level). The degree of optical polarization at that time, Π O ∼ 5%, is comparable to the X-ray measurement. We investigate possible scenarios for these findings, including temporal and geometrical depolarization effects. Unlike some other X-ray-selected BL Lac objects, there is no significant chromatic dependence of the measured polarization in 1ES 1959+650, and its low X-ray polarization may be attributed to turbulence in the jet flow with dynamical timescales shorter than 1 day.

We report the X-ray polarization properties of the high-synchrotron-peaked (HSP) blazar PKS 2155$-$304 based on observations with the Imaging X-ray Polarimetry Explorer (IXPE). We observed the source between Oct 27 and Nov 7, 2023. We also conducted an extensive contemporaneous multiwavelength (MW) campaign. We find that during the first half ($T_1$) of the IXPE pointing, the source exhibited the highest X-ray polarization degree detected for an HSP blazar thus far, (30.7$\pm$2.0)%, which dropped to (15.3$\pm$2.1)% during the second half ($T_2$). The X-ray polarization angle remained stable during the IXPE pointing at 129.4$^\circ$$\pm$1.8$^\circ$ and 125.4$^\circ$$\pm$3.9$^\circ$ during $T_1$ and $T_2$, respectively. Meanwhile, the optical polarization degree remained stable during the IXPE pointing, with average host-galaxy-corrected values of (4.3$\pm$0.7)% and (3.8$\pm$0.9)% during the $T_1$ and $T_2$, respectively. During the IXPE pointing, the optical polarization angle changed achromatically from $\sim$140$^\circ$ to $\sim$90$^\circ$ and back to $\sim$130$^\circ$. Despite several attempts, we only detected (99.7% conf.) the radio polarization once (during $T_2$, at 225.5 GHz): with degree (1.7$\pm$0.4)% and angle 112.5$^\circ$$\pm$5.5$^\circ$. The direction of the broad pc-scale jet is rather ambiguous and has been found to point to the east and south at different epochs; however, on larger scales (> 1.5 pc) the jet points toward the southeast ($\sim$135$^\circ$), similar to all of the MW polarization angles. Moreover, the X-ray to optical polarization degree ratios of $\sim$7 and $\sim$4 during $T_1$ and $T_2$, respectively, are similar to previous IXPE results for several HSP blazars. These findings, combined with the lack of correlation of temporal variability between the MW polarization properties, agree with an energy-stratified shock-acceleration scenario in HSP blazars.

High-frequency-peaked BL Lacs (HBLs) dominate the extragalactic TeV sky, with more than 50 objects detected with the current generation of ground-based TeV gamma-ray observatories. In the last three years, the VERITAS telescope array has observed a flux-limited sample of 36 X-ray selected HBLs with the goal of producing the first unbiased census of TeV emission from HBL blazars. The VERITAS HBL sample contains known TeV sources as well as 15 objects for which TeV emission has not been reported before. The results of this VERITAS campaign include the detection of new TeV blazars as well as unbiased estimates of the TeV flux of HBLs that have previously been reported only during flaring states. The implications of our results in understanding the intrinsic properties of HBLs as a source population will be discussed.

A fraction of the active supermassive black holes at the centers of galaxies in our Universe are capable of launching extreme kiloparsec-long relativistic jets. These jets are known multiband (radio to γ-ray) and multimessenger (neutrino) emitters, and some of them have been monitored over decades at all accessible wavelengths. However, several open questions remain unanswered about the processes powering these highly energetic phenomena. These jets intrinsically produce soft-to-hard X-ray emission that extends from E>0.1keV up to E>100keV, and simultaneous broadband X-ray coverage, combined with excellent timing and imaging capabilities, is required to uncover the physics of jets. Indeed, truly simultaneous soft-to-hard X-ray coverage, in synergy with current and upcoming high-energy facilities (such as IXPE, COSI, CTAO, etc.) and neutrino detectors (e.g., IceCube), would enable us to disentangle the particle population responsible for the high-energy radiation from these jets. A sensitive hard X-ray survey (F20−80keV<10−15ergcm−2s−1) could unveil the bulk of their population in the early Universe. Acceleration and radiative processes responsible for the majority of their X-ray emission would be pinned down by microsecond timing capabilities at both soft and hard X-rays. Furthermore, imaging jet structures for the first time in the hard X-ray regime could unravel the origin of their high-energy emission. The proposed Probe-class mission concept High Energy X-ray Probe (HEX-P) combines all these required capabilities, making it the crucial next-generation X-ray telescope in the multi-messenger, time-domain era. HEX-P will be the ideal mission to unravel the science behind the most powerful accelerators in the Universe.

Blazars can be detected from very large distances due to their high luminosity. However, the detection of gamma-ray emission of blazars beyond z = 3 has only been confirmed for a small number of sources. Such observations probe the growth of supermassive black holes close to the peak of star formation in the history of galaxy evolution. As a result from a continuous monitoring of a sample of 80 z > 3 blazars with the Fermi Large Area Telescope (Fermi-LAT), we present the first detection of a gamma-ray flare from the z = 4.31 blazar TXS 1508+572. This source showed high gamma-ray activity from 2022 February to August, reaching a peak luminosity comparable to the most luminous flares ever detected with Fermi-LAT. We conducted a multiwavelength observing campaign involving XMM-Newton, the Neil Gehrels Swift Observatory, the Effelsberg 100 m radio telescope, and the Very Long Baseline Array. In addition, we make use of the monitoring programs by the Zwicky Transient Facility and the Near-Earth Object Wide-field Infrared Survey Explorer at optical and infrared wavelengths, respectively. We find that the source is particularly variable in the infrared band on daily timescales. The spectral energy distribution collected during our campaign is well described by a one-zone leptonic model, with the gamma-ray flare originating from an increase of external Compton emission as a result of a fresh injection of accelerated electrons.

We present multiwavelength polarization measurements of the luminous blazar Mrk 501 over a 14 month period. The 2-8 keV X-ray polarization was measured with the Imaging X-ray Polarimetry Explorer (IXPE) with six 100 ks observations spanning from 2022 March to 2023 April. Each IXPE observation was accompanied by simultaneous X-ray data from NuSTAR, Swift/XRT, and/or XMM-Newton. Complementary optical-infrared polarization measurements were also available in the B, V, R, I, and J bands, as were radio polarization measurements from 4.85 GHz to 225.5 GHz. Among the first five IXPE observations, we did not find significant variability in the X-ray polarization degree and angle with IXPE. However, the most recent sixth observation found an elevated polarization degree at >3 sigma above the average of the other five observations. The optical and radio measurements show no apparent correlations with the X-ray polarization properties. Throughout the six IXPE observations, the X-ray polarization degree remained higher than, or similar to, the R-band optical polarization degree, which remained higher than the radio value. This is consistent with the energy-stratified shock scenario proposed to explain the first two IXPE observations, in which the polarized X-ray, optical, and radio emission arises from different regions.

The magnetic field conditions in astrophysical relativistic jets can be probed by multiwavelength polarimetry, which has been recently extended to X-rays. For example, one can track how the magnetic field changes in the flow of the radiating particles by observing rotations of the electric vector position angle Ψ. Here we report the discovery of a Ψ_x rotation in the X-ray band in the blazar Mrk 421 at an average flux state. Across the 5 days of Imaging X-ray Polarimetry Explorer (IXPE) observations of 4-6 and 7-9 June 2022, Ψ_x rotated in total by ≥360^∘. Over the two respective date ranges, we find constant, within uncertainties, rotation rates (80 ± 9 and 91 ± 8 ^∘/ day) and polarization degrees (Π_x=10%±1%). Simulations of a random walk of the polarization vector indicate that it is unlikely that such rotation(s) are produced by a stochastic process. The X-ray emitting site does not completely overlap the radio/infrared/optical emission sites, as no similar rotation of Ψ was observed in quasi-simultaneous data at longer wavelengths. We propose that the observed rotation was caused by a helical magnetic structure in the jet, illuminated in the X-rays by a localized shock propagating along this helix. The optically emitting region likely lies in a sheath surrounding an inner spine where the X-ray radiation is released.