Determining the location of $\gamma$-ray emission in blazar jets is a challenging task. Pinpointing the exact location of $\gamma$-ray production within a relativistic jet can place strong constraints on our understanding of high-energy astrophysics and astroparticle physics. We present a study of the radio- and $\gamma$-bright flat-spectrum radio quasar (FSRQ) 4C +01.28 (PKS B1055+018) in which we try to pinpoint the emission site of several prominent GeV flares. This source shows prominent high-amplitude broadband variability on time scales ranging from days to years. We combine high-resolution VLBI observations and multi-band radio light curves over a period of around nine years. We can associate two bright and compact newly ejected jet components with bright flares observed by the Fermi/LAT $\gamma$-ray telescope and at various radio frequencies. A cross-correlation analysis reveals the radio light curves systematically lag behind the $\gamma$-rays. In combination with the jet kinematics as measured by the VLBA, we use these cross-correlations to constrain a model in which the flares become observable at a given frequency when a plasma component passes through the region at which the bulk energy dissipation takes place at that frequency. We derive a lower limit of the location of the $\gamma$-ray emitting region in 4C +01.28 of several parsecs from the jet base, well beyond the expected extent of the broad-line region. This observational limit challenges blazar-emission models that rely on the broad-line region as a source of seed photons for inverse-Compton scattering.

Active Galactic Nuclei (AGN) are some of the brightest and most variable objects in the universe. Those with relativistic jets observed at small viewing angles are blazars. Due to Doppler-boosting, blazars exhibit extreme stochastic variability. While the origin of this variability is thought to be changes in the accretion flow and jet dynamics, much about blazar variability remains unknown. In this paper we use several blazar-dominated AGN samples to form a catalog of 7918 blazars and candidates -- the largest to date. We also collect source types, redshifts, spectral energy distribution (SED) peak frequencies, radio variability Doppler factors, and X-ray flux densities for as many sources as possible. We use all-sky surveys (CRTS, ATLAS, and ZTF, abbreviated as ``CAZ'') to extract their optical multiband flux density on a nightly basis between 2007 and 2023, and construct as long and as high cadence light curves as possible for as many sources as attainable. The catalog and its light curves are provided in the accompanying electronic tables, enabling many analyses involving AGN variability with unprecedented sample-sizes. We quantify the variability of the light curves, and apply the Bayesian blocks algorithm to determine their flaring periods. We find: (1) optical flares generally have a faster rise than decay; (2) optical brightness and variability are strongly dependent on the synchrotron peak frequency; (3) flat spectrum radio quasars and BL Lac objects have comparable optical variability and flare characteristics at the same synchrotron peak frequency; and (4) optical flare times tend to decrease and amplitudes increase with increasing radio variability Doppler factor.

Over the past decade, the IceCube Neutrino Observatory has detected a few hundreds of high-energy (HE) neutrinos from cosmic sources. Despite numerous studies searching for their origin, it is still not known which source populations emit them. A few confident individual associations exist with active galactic nuclei (AGN), mostly with blazars which are jetted AGN whose jet points in our direction. Nonetheless, on a population level, blazar-neutrino correlation strengths are rather weak. This could mean that blazars as a population do not emit HE neutrinos, or that the detection power of the tests is insufficient due to the strong atmospheric neutrino background. By assuming an increase in HE neutrino emission during major blazar flares, in our previous studies we leveraged the arrival time of the neutrinos to boost the detection power. In this paper we utilize the same principle while substantially increasing the number of blazars. We search for the spatio-temporal correlation of 356 IceCube HE neutrinos with major optical flares of 3225 radio- and 3814 $\gamma$-ray-selected blazars. We find that, despite the increase in data size, the number of confident spatio-temporal associations remains low and the overall correlation strengths weak. Two individual associations drive our strongest and the only $>$2$\sigma$ post-trial spatio-temporal correlation, occurring with the BL Lac objects of the radio-selected blazar sample. We estimate that $\lesssim$8\% of the detected cosmic neutrinos were emitted by blazars during major optical flares. As a complementary analysis, we compare the synchrotron peak frequency, redshift, Doppler factor, X-ray brightness, and optical variability of spatially neutrino-associated blazars to those of the general blazar population. We find that spatially neutrino-associated blazars of the tested samples have higher than average Doppler factor and X-ray brightness.

In this study, we demonstrate some of the caveats in common statistical methods used for analysing astronomical variability timescales. We consider these issues specifically in the context of active galactic nuclei (AGNs) and use a more practical approach compared to mathematics literature, where the number of formulae may sometimes be overwhelming. We conducted a thorough literature review both on the statistical properties of light-curve data, specifically in the context of sampling effects, as well as on the methods used to analyse them. We simulated a wide range of data to test some of the known issues in AGN variability analysis as well as to investigate previously unknown or undocumented caveats. We discovered problems with some commonly used methods and confirmed how challenging it is to identify timescales from observed data. We find that interpolation of a light curve with biased sampling, specifically with bias towards flaring events, affects its measured power spectral density in a different manner than those of simulated light curves. We also find that an algorithm aiming to match the probability density function of a light curve has often been used incorrectly. These new issues appear to have been mostly overlooked and not necessarily addressed before, especially in astronomy literature.

The Event Horizon Telescope (EHT) observations carried out in 2018 April at 1.3 mm wavelengths included 9 stations in the array, comprising 7 single-dish telescopes and 2 phased arrays. The metadata package for the 2018 EHT observing campaign contains calibration tables required for the a-priori amplitude calibration of the 2018 April visibility data. This memo is the official documentation accompanying the release of the 2018 EHT metadata package, providing an overview of the contents of the package. We describe how telescope sensitivities, gain curves and other relevant parameters for each station in the EHT array were collected, processed, and validated to produce the calibration tables.

We present observations of the blazar 3C 279 at 22 GHz using the space VLBI mission RadioAstron on 2018 January 15. Images in both total intensity and fractional polarization are reconstructed using RML method implemented in the eht-imaging library. The electric vector position angles are found to be mostly aligned with the general jet direction, suggesting a predominantly toroidal magnetic field, in agreement with the presence of a helical magnetic field. Ground-space fringes were detected up to a projected baseline length of $\sim 8$G$\lambda$, achieving the angular resolution of around 26$\mu$as. The fine-scale structure of the relativistic jet is found in our study extending to a projected distance of $\sim 180$ parsec from the radio core. However, the filamentary structure reported by previous RadioAstron observations of 2014 is not detected in our current study. We discuss potential causes for this phenomenon, together with a comparison using public 43 GHz data from the BEAM-ME program, showing a significant drop in the jet's total intensity. The optically thick core is observed with a brightness temperature of $1.6 \times 10^{12}$ K, consistent with equipartition between the energy densities of the relativistic particles and the magnetic field. This yields an estimated magnetic field strength of 0.2 G.

We examine lengthy radio light curves of the flat spectrum radio galaxy 3C 454.3 for possible quasi-periodic oscillations (QPOs). The data used in this work were collected at five radio frequencies, 4.8, 8.0, 14.5, 22.0, and 37.0 GHz between 1979--2013 as observed at the University of Michigan Radio Astronomical Observatory, Crimean Astrophysical Observatory, and Aalto University Mets{ä}hovi Radio Observatory. We employ generalized Lomb-Scargle periodogram and weighted wavelet transform analyses to search for periodicities in these light curves. We confirm a QPO period of $\sim$ 2000 day to be at least 4$\sigma$ significant using both methods at all five radio frequencies between 1979 and 2007, after which a strong flare changed the character of the light curve. We also find a $\sim$~600 day period which is at least 4$\sigma$ significant, but only in the 22.0 and 37.0 GHz light curves. We briefly discuss physical mechanisms capable of producing such variations.

AT2022rze is a luminous, ambiguous transient located South-East of the geometric center of its host galaxy at redshift z = 0.08. The host appears to be formed by a merging galaxy system. The observed characteristics of AT2022rze are reminiscent of active galactic nuclei (AGN), tidal disruption events (TDEs), and superluminous supernovae (SLSNe). The transient reached a peak absolute magnitude of -20.2 +- 0.2 mag, showing a sharp rise (trise,1/e = 27.5 +- 0.6 days) followed by a slow decline (tdec,1/e = 382.9 +- 0.6). Its bumpy light curve and narrow Balmer lines indicate the presence of gas (and dust). Its light curve shows rather red colors, indicating that the transient could be affected by significant host extinction. The spectra reveal coronal lines, indicative of high-energy (X-ray/UV) emission. Archival data reveal no prior activity at this location, disfavoring a steady-state AGN, although an optical spectrum obtained prior to the transient is consistent with an AGN classification of the host. Based on this, we conclude that the transient most likely represents a Changing-look AGN at the center of the smallest component of the merging system.

Studying the long-term radio variability (timescales of months to years) of blazars enables us to gain a better understanding of the physical structure of these objects on sub-parsec scales, and the physics of super massive black holes. In this study, we focus on the radio variability of 1157 blazars observed at 15 GHz through the Owens Valley Radio Observatory (OVRO) Blazar Monitoring Program. We investigate the dependence of the variability amplitudes and timescales, characterized based on model fitting to the structure functions, on the milliarcsecond core sizes measured by Very Long Baseline Interferometry. We find that the most compact sources at milliarcsecond scales exhibit larger variability amplitudes and shorter variability timescales than more extended sources. Additionally, for sources with measured redshifts and Doppler boosting factors, the correlation between linear core sizes against variability amplitudes and intrinsic timescales are also significant. The observed relationship between variability timescales and core sizes is expected, based on light travel-time arguments. This variability vs core size relation extends beyond the core sizes measured at 15 GHz; we see significant correlation between the 15 GHz variability amplitudes (as well as timescales) and core sizes measured at other frequencies, which can be attributed to a frequency-source size relationship arising from the intrinsic jet structure. At low frequencies of 1 GHz where the core sizes are dominated by interstellar scattering, we find that the variability amplitudes have significant correlation with the 1 GHz intrinsic core angular sizes, once the scatter broadening effects are deconvoluted from the intrinsic core sizes.

Narrow-line Seyfert 1 (NLS1) galaxies are a class of active galactic nuclei (AGN) that, in some cases, can harbor powerful relativistic jets. One of them, PKS 2004-447, shows gamma-ray emission, and underwent its first recorded multifrequency flare in 2019. However, past studies revealed that in radio this source can be classified as a compact steep-spectrum source (CSS), suggesting that, unlike other gamma-ray sources, the relativistic jets of PKS 2004-447 have a large inclination with respect to the line of sight. We present here a set of spectroscopic observations of this object, aimed at carefully measuring its black hole mass and Eddington ratio, determining the properties of its emission lines, and characterizing its long term variability. We find that the black hole mass is $(1.5\pm0.2)\times10^7$ M$_\odot$, and the Eddington ratio is 0.08. Both values are within the typical range of NLS1s. The spectra also suggest that the 2019 flare was caused mainly by the relativistic jet, while the accretion disk played a minor role during the event. In conclusion, we confirm that PKS 2004-447 is one of the rare examples of gamma-ray emitting CSS/NLS1s hybrid, and that these two classes of objects are likely connected in the framework of AGN evolution.

We present our observations and analysis of SN~2023gpw, a hydrogen-rich superluminous supernova (SLSN~II) with broad emission lines in its post-peak spectra. Unlike previously observed SLSNe~II, its light curve suggests an abrupt drop during a solar conjunction between $\sim$80 and $\sim$180~d after the light-curve peak, possibly analogous to a normal hydrogen-rich supernova (SN). Spectra taken at and before the peak show hydrogen and helium `flash' emission lines attributed to early interaction with a dense confined circumstellar medium (CSM). A well-observed ultraviolet excess appears as these lines disappear, also as a result of CSM interaction. The blackbody photosphere expands roughly at the same velocity throughout the observations, indicating little or no bulk deceleration. This velocity is much higher than what is seen in spectral lines, suggesting asymmetry in the ejecta. The high total radiated energy ($\gtrsim9\times10^{50}$~erg) and aforementioned lack of bulk deceleration in SN~2023gpw are difficult to reconcile with a neutrino-driven SN simply combined with efficient conversion from kinetic energy to emission through interaction. This suggests an additional energy source such as a central engine. While magnetar-powered models qualitatively similar to SN~2023gpw exist, more modeling work is required to determine if they can reproduce the observed properties in combination with early interaction. The required energy might alternatively be provided by accretion onto a black hole created in the collapse of a massive progenitor star.

We studied the polarization behavior of the quasar 3C 273 over the 1 mm wavelength band at ALMA with a total bandwidth of 7.5GHz across 223 to 243 GHz at 0.8 arcsec resolution, corresponding to 2.1 kpc at the distance of 3C 273. With these observations we were able to probe the optically thin polarized emission close to the jet base, and constrain the magnetic field structure. We computed the Faraday rotation measure using simple linear fitting and Faraday rotation measure synthesis. In addition, we modeled the broadband behavior of the fractional Stokes Q and U parameters (qu-fitting). The systematic uncertainties in the polarization observations at ALMA were assessed through Monte Carlo simulations. We find the unresolved core of 3C 273 to be 1.8% linearly polarized. We detect a very high rotation measure (RM) of ~5.0 x10^5 rad/m^2 over the 1 mm band when assuming a single polarized component and an external RM screen. This results in a rotation of >40 deg of the intrinsic electric vector position angle, which is significantly higher than typically assumed for millimeter wavelengths. The polarization fraction increases as a function of wavelength, which according to our qu-fitting could be due to multiple polarized components of different Faraday depth within our beam, or internal Faraday rotation. With our limited wavelength coverage, we cannot distinguish between the cases, and additional multifrequency and high angular resolution observations are needed to determine the location and structure of the magnetic field of the Faraday active region. Comparing our RM estimate with values obtained at lower frequencies, the RM increases as a function of observing frequency, following a power law with an index of ~2.0 consistent with a sheath surrounding a conically expanding jet. We also detect ~0.2% circular polarization, although further observations are needed to confirm this result.

PKS 1830$-$211 is a $\gamma$-ray emitting, high-redshift (z $= 2.507 \pm 0.002$), lensed flat-spectrum radio quasar. During the period mid-February to mid-April 2019, this source underwent a series of strong $\gamma$-ray flares that were detected by both AGILE-GRID and Fermi-LAT, reaching a maximum $\gamma$-ray flux of $F_{\rm E>100 MeV}\approx 2.3\times10^{-5}$ ph cm$^{-2}$ s$^{-1}$. Here we report on a coordinated campaign from both on-ground (Medicina, OVRO, REM, SRT) and orbiting facilities (AGILE, Fermi, INTEGRAL, NuSTAR, Swift, Chandra), with the aim of investigating the multi-wavelength properties of PKS 1830$-$211 through nearly simultaneous observations presented here for the first time. We find a possible break in the radio spectra in different epochs above 15 GHz, and a clear maximum of the 15 GHz data approximately 110 days after the $\gamma$-ray main activity periods. The spectral energy distribution shows a very pronounced Compton dominance (> 200) which challenges the canonical one-component emission model. Therefore we propose that the cooled electrons of the first component are re-accelerated to a second component by, e.g., kink or tearing instability during the $\gamma$-ray flaring periods. We also note that PKS 1830$-$211 could be a promising candidate for future observations with both Compton satellites (e.g., e-ASTROGAM) and Cherenkov arrays (CTAO) which will help, thanks to their improved sensitivity, in extending the data availability in energy bands currently uncovered.

The nature of Type Ia supernova (SN Ia) progenitor systems and the mechanisms that lead up to their explosions are still widely debated. In rare cases the SN ejecta interact with circumstellar material (CSM) that was ejected from the progenitor system prior to the SN. The unknown distance between the CSM and SN explosion site makes it impossible to predict when the interaction will start. If the time between the SN and start of CSM interaction is of the order of months to years the SN has generally faded and is not actively followed up anymore, making it even more difficult to detect the interaction while it happens. Here we report on a real-time monitoring program which ran between 13-11-2023 and 09-07-2024, monitoring 6914 SNe Ia for signs of late-time rebrightening using the Zwicky Transient Facility (ZTF). Flagged candidates were rapidly followed up with photometry and spectroscopy to confirm the late-time excess and its position. We report the discovery of a $\sim50$ day rebrightening event in SN 2020qxz around 1200 days after the peak of its light curve. SN 2020qxz had signs of early CSM interaction but faded from view over 2 years before its reappearance. Follow-up spectroscopy revealed 4 emission lines that faded shortly after the end of the ZTF detected rebrightening. Our best match for these emission lines are H$\beta$ (blue shifted by $\sim5900$ km s$^{-1}$) and CaII$_{\lambda8542}$, NI$_{\lambda8567}$, and KI$_{\lambda\lambda 8763, 8767}$, all blue shifted by 5100 km s$^{-1}$ (although we note that these identifications are uncertain). This shows that catching and following up on late-time interactions as they occur can give new clues about the nature of the progenitor systems that produce these SNe by putting constraints on the possible type of donor star, and the only way to do this systematically is to use large sky surveys such as ZTF to monitor a large sample of objects.

It has been a decade since the IceCube collaboration began detecting high-energy (HE) neutrinos originating from cosmic sources. Despite a few well-known individual associations and numerous phenomenological, observational, and statistical multiwavelength studies, the origin of astrophysical HE neutrinos largely remains a mystery. To date, the most convincing associations link HE neutrinos with active galactic nuclei (AGN). Consequently, many studies have attempted population-based correlation tests between HE neutrinos and specific AGN subpopulations (such as blazars). While some of associations are suggestive, no definitive population-based correlation has been established. This could result from either a lack of a population-based correlation or insufficient detection power, given the substantial atmospheric neutrino background. By leveraging blazar variability, we performed spatio-temporal blazar-neutrino correlation tests aimed at enhancing detection power by reliably incorporating temporal information into the statistical analysis. We used simulations to evaluate the detection power of our method under various test strategies. We find that: (1) with sufficiently large source samples, if 20% of astrophysical HE neutrinos originate from blazars, we should robustly observe $\sim$4$\sigma$ associations; (2) a counting-based test-statistic combined with a top-hat weighting scheme (rather than a Gaussian one) provides the greatest detection power; (3) applying neutrino sample cuts reduces detection power when a weighting scheme is used; (4) in top-hat-like weighting schemes, small p-values do not occur arbitrarily with an increase in HE neutrino error region size (any such occurrence is indicative of an underlying blazar-neutrino correlation).

Many active galaxies harbor powerful relativistic jets, however, the detailed mechanisms of their formation and acceleration remain poorly understood. To investigate the area of jet acceleration and collimation with the highest available angular resolution, we study the innermost region of the bipolar jet in the nearby low-ionization nuclear emission-line region (LINER) galaxy NGC 1052. We combined observations of NGC 1052 taken with VLBA, GMVA, and EHT over one week in the spring of 2017. For the first time, NGC 1052 was detected with the EHT, providing a size of the central region in-between both jet bases of 250 RS (Schwarzschild radii) perpendicular to the jet axes. This size estimate supports previous studies of the jets expansion profile which suggest two breaks of the profile at around 300 RS and 10000 RS distances to the core. Furthermore, we estimated the magnetic field to be 1.25 Gauss at a distance of 22 {\mu}as from the central engine by fitting a synchrotron-self absorption spectrum to the innermost emission feature, which shows a spectral turn-over at about 130 GHz. Assuming a purely poloidal magnetic field, this implies an upper limit on the magnetic field strength at the event horizon of 26000 Gauss, which is consistent with previous measurements. The complex, low-brightness, double-sided jet structure in NGC 1052 makes it a challenge to detect the source at millimeter (mm) wavelengths. However, our first EHT observations have demonstrated that detection is possible up to at least 230 GHz. This study offers a glimpse through the dense surrounding torus and into the innermost central region, where the jets are formed. This has enabled us to finally resolve this region and provide improved constraints on its expansion and magnetic field strength.

The FSRQ 4C+01.28 is a bright and highly variable radio and $\gamma$-ray emitter. We aim to pinpoint the location of the $\gamma$-ray emitting region within its jet in order to derive strong constraints on $\gamma$-ray emission models for blazar jets. We use radio and $\gamma$-ray data obtained with ALMA, OVRO, SMA and Fermi/LAT to study the cross-correlation between $\gamma$-ray and multi-frequency radio light curves. Moreover, we employ VLBA observations at 43 GHz over a period of around nine years to study the parsec-scale jet kinematics. To pinpoint the location of the $\gamma$-ray emitting region, we use a model in which outbursts shown in the $\gamma$-ray and radio light curves are produced when moving jet components pass through the $\gamma$-ray emitting and the radio core regions. We find two bright and compact newly ejected jet components that are likely associated with a high activity period visible in the $\gamma$-ray and radio light curves. The kinematic analysis of the VLBA observations leads to a maximum apparent jet speed of $\beta_{app}=19\pm10$ and an upper limit on the viewing angle of $\phi$ < 4 deg. We determine the power law indices that are characterizing the jet geometry, brightness temperature distribution, and core shift to be $l=0.974\pm0.098$, $s=-3.31\pm0.31$, and $k_r=1.09\pm0.17$, which are in agreement with a conical jet in equipartition. A cross-correlation analysis shows that the radio light curves follow the $\gamma$-ray light curve. We pinpoint the location of the $\gamma$-ray emitting region with respect to the jet base to the range of $2.6\,\mathrm{pc}\leq d_\gamma\leq20\,\mathrm{pc}$. Our derived observational limits places the location of $\gamma$-ray production in 4C+01.28 beyond the expected extent of the broad-line region (BLR) and therefore challenges blazar-emission models that rely on inverse Compton up-scattering of seed photons from the BLR.

CTA 102 is a $\gamma$-ray bright blazar that exhibited multiple flares in observations by the Large Area Telescope on board the Fermi Gamma-Ray Space Telescope during the period of 2016-2018. We present results from the analysis of multi-wavelength light curves aiming at revealing the nature of $\gamma$-ray flares from the relativistic jet in the blazar. We analyse radio, optical, X-ray, and $\gamma$-ray data obtained in a period from 2012 September 29 to 2018 October 8. We identify six flares in the $\gamma$-ray light curve, showing a harder-when-brighter-trend in the $\gamma$-ray spectra. We perform a cross-correlation analysis of the multi-wavelength light curves. We find nearly zero time lags between the $\gamma$-ray and optical and X-ray light curves, implying a common spatial origin for the emission in these bands. We find significant correlations between the $\gamma$-ray and radio light curves as well as negative/positive time lags with the $\gamma$-ray emission lagging/leading the radio during different flaring periods. The time lags between $\gamma$-ray and radio emission propose the presence of multiple $\gamma$-ray emission sites in the source. As seen in 43 GHz images from the Very Long Baseline Array, two moving disturbances (or shocks) were newly ejected from the radio core. The $\gamma$-ray flares from 2016 to 2017 are temporally coincident with the interaction between a traveling shock and a quasi-stationary one at $\sim$0.1 mas from the core. The other shock is found to emerge from the core nearly simultaneous with the $\gamma$-ray flare in 2018. Our results suggest that the $\gamma$-ray flares originated from shock-shock interactions.

Concomitant with the Imaging X-ray Polarimetry Explorer (IXPE) observation of the Perseus cluster, we obtained optical spectropolarimetry of its central active galactic nucleus, NGC 1275, using the Alhambra Faint Object Spectrograph and Camera (ALFOSC) on the Nordic Optical Telescope (NOT). While the total-light spectrum confirms its edge-on, core obscured (type-2) classification, the polarized spectrum shows a polarization angle aligned with the arcsecond radio jet axis -- an exceptional behavior for type-2 objects. Our polarization analysis also reveals wavelength-dependent linear polarization at level 2-3% in the continuum, likely rising from a combination of variable syn

The origin of the astrophysical neutrino flux discovered by IceCube remains largely unknown. Several individual neutrino source candidates were observed. Among them is the gamma-ray flaring blazar TXS 0506+056. A similar coincidence of a high-energy neutrino and a gamma-ray flare was found in blazar PKS 0735+178. By modeling the spectral energy distributions of PKS 0735+178, we expect to investigate the physical conditions for neutrino production during different stages of the source activity. We analyze the multi-wavelength data during the selected periods of time. Using numerical simulations of radiation processes in the source, we study the parameter space of one-zone leptonic and leptohadronic models and find the best-fit solutions that explain the observed photon fluxes. We show the impact of model parameter degeneracy on the prediction of the neutrino spectra. We show that the available mutli-wavelength data are not sufficient to predict the neutrino spectrum unambiguously. Still, under the condition of maximal neutrino flux, we propose a scenario in which 0.2 neutrino events are produced during the 50 days flare.