Showing 1–50 of 125 results for author: Smith, P S

A gamma-ray flare from TXS 1508+572: characterizing the jet of a $z=4.31$ blazar in the early Universe

Authors: Andrea Gokus, Markus Böttcher, Manel Errando, Michael Kreter, Jonas Heßdörfer, Florian Eppel, Matthias Kadler, Paul S. Smith, Petra Benke, Leonid I. Gurvits, Alex Kraus, Mikhail Lisakov, Felicia McBride, Eduardo Ros, Florian Rösch, Jörn Wilms

Abstract: Blazars can be detected from very large distances due to their high luminosity. However, the detection of $γ$-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… ▽ More Blazars can be detected from very large distances due to their high luminosity. However, the detection of $γ$-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 Fermi-LAT, we present the first detection of a $γ$-ray flare from the $z=4.31$ blazar TXS 1508+572. This source showed high $γ$-ray activity from February to August 2022, reaching a peak luminosity comparable to the most luminous flares ever detected with Fermi -LAT. We conducted a multiwavelength observing campaign involving XMM-Newton, Swift, 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 NEOWISE at optical and infrared wavelengths, respectively. We find that the source is particularly variable in the infrared band on daily time scales. The spectral energy distribution collected during our campaign is well described by a one-zone leptonic model, with the $γ$-ray flare originating from an increase of external Compton emission as a result of a fresh injection of accelerated electrons. △ Less

Submitted 1 August, 2024; v1 submitted 11 June, 2024; originally announced June 2024.

Comments: 19 pages, 8 figures; accepted by ApJ on July 30 - Initial version on arXiv is the submitted one

The variability patterns of the TeV blazar PG 1553+113 from a decade of MAGIC and multi-band observations

Authors: MAGIC Collaboration, H. Abe, S. Abe, J. Abhir, V. A. Acciari, I. Agudo, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, C. Arcaro, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, I. Batković, J. Baxter, J. Becerra González, E. Bernardini, J. Bernete, A. Berti, J. Besenrieder, C. Bigongiari , et al. (242 additional authors not shown)

Abstract: PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m a… ▽ More PG 1553+113 is one of the few blazars with a convincing quasi-periodic emission in the gamma-ray band. The source is also a very high-energy (VHE; >100 GeV) gamma-ray emitter. To better understand its properties and identify the underlying physical processes driving its variability, the MAGIC Collaboration initiated a multiyear, multiwavelength monitoring campaign in 2015 involving the OVRO 40-m and Medicina radio telescopes, REM, KVA, and the MAGIC telescopes, Swift and Fermi satellites, and the WEBT network. The analysis presented in this paper uses data until 2017 and focuses on the characterization of the variability. The gamma-ray data show a (hint of a) periodic signal compatible with literature, but the X-ray and VHE gamma-ray data do not show statistical evidence for a periodic signal. In other bands, the data are compatible with the gamma-ray period, but with a relatively high p-value. The complex connection between the low and high-energy emission and the non-monochromatic modulation and changes in flux suggests that a simple one-zone model is unable to explain all the variability. Instead, a model including a periodic component along with multiple emission zones is required. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Monthly Notices of the Royal Astronomical Society. 19 pages, 9 figures. Corresponding authors: Elisa Prandini, Antonio Stamerra, Talvikki Hovatta

VERITAS discovery of very high energy gamma-ray emission from S3 1227+25 and multiwavelength observations

Authors: Atreya Acharyya, Colin Adams, Avery Archer, Priyadarshini Bangale, Wystan Benbow, Aryeh Brill, Jodi Christiansen, Alisha Chromey, Manel Errando, Abe Falcone, Qi Feng, John Finley, Gregory Foote, Lucy Fortson, Amy Furniss, Greg Gallagher, William Hanlon, David Hanna, Olivier Hervet, Claire Hinrichs, John Hoang, Jamie Holder, Weidong Jin, Madalyn Johnson, Philip Kaaret , et al. (46 additional authors not shown)

Abstract: We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th… ▽ More We report the detection of very high energy gamma-ray emission from the blazar S3 1227+25 (VER J1230+253) with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). VERITAS observations of the source were triggered by the detection of a hard-spectrum GeV flare on May 15, 2015 with the Fermi-Large Area Telescope (LAT). A combined five-hour VERITAS exposure on May 16th and May 18th resulted in a strong 13$σ$ detection with a differential photon spectral index, $Γ$ = 3.8 $\pm$ 0.4, and a flux level at 9% of the Crab Nebula above 120 GeV. This also triggered target of opportunity observations with Swift, optical photometry, polarimetry and radio measurements, also presented in this work, in addition to the VERITAS and Fermi-LAT data. A temporal analysis of the gamma-ray flux during this period finds evidence of a shortest variability timescale of $τ_{obs}$ = 6.2 $\pm$ 0.9 hours, indicating emission from compact regions within the jet, and the combined gamma-ray spectrum shows no strong evidence of a spectral cut-off. An investigation into correlations between the multiwavelength observations found evidence of optical and gamma-ray correlations, suggesting a single-zone model of emission. Finally, the multiwavelength spectral energy distribution is well described by a simple one-zone leptonic synchrotron self-Compton radiation model. △ Less

Submitted 4 May, 2023; originally announced May 2023.

Comments: 18 pages, 6 figures. Accepted for publication in the Astrophysical Journal (ApJ)

The Sloan Digital Sky Survey Reverberation Mapping Project: Key Results

Authors: Yue Shen, Catherine J. Grier, Keith Horne, Zachary Stone, Jennifer I. Li, Qian Yang, Yasaman Homayouni, Jonathan R. Trump, Scott F. Anderson, W. N. Brandt, Patrick B. Hall, Luis C. Ho, Linhua Jiang, Patrick Petitjean, Donald P. Schneider, Charling Tao, Fergus. R. Donnan, Yusra AlSayyad, Matthew A. Bershady, Michael R. Blanton, Dmitry Bizyaev, Kevin Bundy, Yuguang Chen, Megan C. Davis, Kyle Dawson , et al. (22 additional authors not shown)

Abstract: We present the final data from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a precursor to the SDSS-V Black Hole Mapper Reverberation Mapping program. This data set includes 11-year photometric and 7-year spectroscopic light curves for 849 broad-line quasars over a redshift range of 0.1<z<4.5 and a luminosity range of Lbol=1E44-47.5 erg/s, along with spectral and variabili… ▽ More We present the final data from the Sloan Digital Sky Survey Reverberation Mapping (SDSS-RM) project, a precursor to the SDSS-V Black Hole Mapper Reverberation Mapping program. This data set includes 11-year photometric and 7-year spectroscopic light curves for 849 broad-line quasars over a redshift range of 0.1<z<4.5 and a luminosity range of Lbol=1E44-47.5 erg/s, along with spectral and variability measurements. We report 23, 81, 125, and 110 reverberation mapping lags (relative to optical continuum variability) for broad Halpha, Hbeta, MgII and CIV using the SDSS-RM sample, spanning much of the luminosity and redshift ranges of the sample. Using 30 low-redshift RM AGNs with dynamical-modeling black hole masses, we derive a new estimate of the average virial factor of <log f>=0.62+-0.07 for the line dispersion measured from the RMS spectrum. The intrinsic scatter of individual virial factors is 0.31+-0.07 dex, indicating a factor of two systematic uncertainty in RM black hole masses. Our lag measurements reveal significant R-L relations for Hbeta and MgII at high redshift, consistent with the latest measurements based on heterogeneous samples. While we are unable to robustly constrain the slope of the R-L relation for CIV given the limited dynamical range in luminosity, we found substantially larger scatter in CIV lags at fixed L1350. Using the SDSS-RM lag sample, we derive improved single-epoch (SE) mass recipes for Hbeta, MgII and CIV, which are consistent with their respective RM masses as well as between the SE recipes from two different lines, over the luminosity range probed by our sample. The new Hbeta and MgII recipes are approximately unbiased estimators at given RM masses, but there are systematic biases in the CIV recipe. The intrinsic scatter of SE masses around RM masses is ~0.45 dex for Hbeta and MgII, increasing to ~0.58 dex for CIV. △ Less

Submitted 1 April, 2024; v1 submitted 1 May, 2023; originally announced May 2023.

Comments: Replaced with accepted version (ApJS in press). All measurements remain unchanged from the previous version. 38 pages. Data products available at https://ariel.astro.illinois.edu/sdssrm/final result/ and ftp://quasar.astro.illinois.edu/public/sdssrm/final_result/

Multi-Epoch Spectropolarimetry for a Sample of Type IIn Supernovae: Persistent Asymmetry in Dusty Circumstellar Material

Authors: Christopher Bilinski, Nathan Smith, G. Grant Williams, Paul S. Smith, Douglas C. Leonard, Jennifer L. Hoffman, Jennifer Andrews, Peter Milne

Abstract: We present multi-epoch spectropolarimetry and spectra for a sample of 14 Type IIn supernovae (SNe IIn). We find that after correcting for likely interstellar polarization, SNe IIn commonly show intrinsic continuum polarization of 1--3% at the time of peak optical luminosity, although a few show weaker or negligible polarization. While some SNe IIn have even stronger polarization at early times, th… ▽ More We present multi-epoch spectropolarimetry and spectra for a sample of 14 Type IIn supernovae (SNe IIn). We find that after correcting for likely interstellar polarization, SNe IIn commonly show intrinsic continuum polarization of 1--3% at the time of peak optical luminosity, although a few show weaker or negligible polarization. While some SNe IIn have even stronger polarization at early times, their polarization tends to drop smoothly over several hundred days after peak. We find a tendency for the intrinsic polarization to be stronger at bluer wavelengths, especially at early times. While polarization from an electron scattering region is expected to be grey, scattering of SN light by dusty circumstellar material (CSM) may induce such a wavelength-dependent polarization. For most SNe IIn, changes in polarization degree and wavelength dependence are not accompanied by changes in the position angle, requiring that asymmetric pre-SN mass loss had a persistent geometry. While 2--3% polarization is typical, about 30% of SNe IIn have very low or undetected polarization. Under the simplifying assumption that all SN IIn progenitors have axisymmetric CSM (i.e. disk/torus/bipolar), then the distribution of polarization values we observe is consistent with similarly asymmetric CSM seen from a distribution of random viewing angles. This asymmetry has very important implications for understanding the origin of pre-SN mass loss in SNe IIn, suggesting that it was shaped by binary interaction. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: 76 pages, 54 figures (13 in main text, 41 in appendix A)

Record-breaking polarization from the interacting superluminous supernova 2017hcc

Authors: Jon C. Mauerhan, Nathan Smith, G. Grant Williams, Paul S. Smith, Alexei V. Filippenko, Christopher Bilinski, WeiKang Zheng, Thomas G. Brink, Jennifer L. Hoffman, Douglas C. Leonard, Peter Milne, Benjamin Jeffers, Shaunak Modak, Samantha Stegman, Keto D. Zhang

Abstract: We present multiepoch spectropolarimetry of the superluminous interacting Type IIn supernova SN2017hcc, covering 16 to 391 days after explosion. In our first epoch we measure continuum polarization as high as 6%, making SN 2017hcc the most intrinsically polarized SN ever reported. During the first 29 days of coverage, when the polarization is strongest, the continuum polarization has a wavelength… ▽ More We present multiepoch spectropolarimetry of the superluminous interacting Type IIn supernova SN2017hcc, covering 16 to 391 days after explosion. In our first epoch we measure continuum polarization as high as 6%, making SN 2017hcc the most intrinsically polarized SN ever reported. During the first 29 days of coverage, when the polarization is strongest, the continuum polarization has a wavelength dependence that rises toward blue wavelengths, and becomes wavelength independent by day 45. The polarization strength drops rapidly during the first month, even as the SN flux is still climbing to peak brightness. Nonetheless, record-high polarization is maintained until day 68, at which point the source polarization declines to 1.9%, comparable to peak levels in previous well-studied SNe IIn. Thereafter the SN continues in polarization decline, while exhibiting only minor changes in position angle on the sky. The blue slope of the polarized continuum during the first month, accompanied by short-lived polarized flux for Balmer emission, suggests that an aspherical distribution of dust grains in pre-shock circumstellar material (CSM) is echoing the SN IIn spectrum and strongly influencing the polarization, while the subsequent decline during the wavelength-independent phase appears broadly consistent with electron scattering near the SN/CSM interface. The persistence of the polarization position angle between these two phases suggests that the pre-existing CSM responsible for the dust scattering at early times is part of the same geometric structure as the electron-scattering region that dominates the polarization at later times. SN2017hcc appears to be yet another, but much more extreme, case of aspherical yet well-ordered CSM in Type IIn SNe, possibly resulting from pre-SN mass loss shaped by a binary progenitor system. △ Less

Submitted 24 April, 2023; originally announced April 2023.

Comments: Submitted to MNRAS April 21, 2023

A High-Velocity Scatterer Revealed in the Thinning Ejecta of a Type II Supernova

Authors: Douglas C. Leonard, Luc Dessart, D. John Hillier, Giuliano Pignata, G. Grant Williams, Jennifer L. Hoffman, Peter Milne, Nathan Smith, Paul S. Smith, Harish G. Khandrika

Abstract: We present deep, nebular-phase spectropolarimetry of the Type II-P/L SN 2013ej, obtained 167 days after explosion with the European Southern Observatory's Very Large Telescope. The polarized flux spectrum appears as a nearly perfect (92% correlation), redshifted (by ~4,000 km/sec) replica of the total flux spectrum. Such a striking correspondence has never been observed before in nebular-phase sup… ▽ More We present deep, nebular-phase spectropolarimetry of the Type II-P/L SN 2013ej, obtained 167 days after explosion with the European Southern Observatory's Very Large Telescope. The polarized flux spectrum appears as a nearly perfect (92% correlation), redshifted (by ~4,000 km/sec) replica of the total flux spectrum. Such a striking correspondence has never been observed before in nebular-phase supernova spectropolarimetry, although data capable of revealing it have heretofore been only rarely obtained. Through comparison with 2D polarized radiative transfer simulations of stellar explosions, we demonstrate that localized ionization produced by the decay of a high-velocity, spatially confined clump of radioactive 56-Ni -- synthesized by and launched as part of the explosion with final radial velocity exceeding 4,500 km/sec -- can reproduce the observations through enhanced electron scattering. Additional data taken earlier in the nebular phase (day 134) yield a similarly strong correlation (84%) and redshift, whereas photospheric-phase epochs that sample days 8 through 97, do not. This suggests that the primary polarization signatures of the high-velocity scattering source only come to dominate once the thick, initially opaque hydrogen envelope has turned sufficiently transparent. This detection in an otherwise fairly typical core-collapse supernova adds to the growing body of evidence supporting strong asymmetries across Nature's most common types of stellar explosions, and establishes the power of polarized flux -- and the specific information encoded by it in line photons at nebular epochs -- as a vital tool in such investigations going forward. △ Less

Submitted 20 October, 2021; originally announced October 2021.

Comments: Accepted for publication in The Astrophysical Journal Letters

Mrk 1239: a Type-2 Counterpart of Narrow-line Seyfert-1?

Authors: Xiang Pan, Hongyan Zhou, Chenwei Yang, Luming Sun, Paul S. Smith, Tuo Ji, Ning Jiang, Peng Jiang, Wenjuan Liu, Honglin Lu, Xiheng Shi, Xuejie Dai, Shaohua Zhang

Abstract: We present new spectrophotometric and spectropolarimetric observations of Mrk 1239, one of the 8 prototypes that defines type-1 narrow-line Seyfert galaxies (NLS1s). Unlike the other typical NLS1s though, a high degree of polarization ($P\sim$5.6%) and red optical-IR ($g-W_4$ = 12.35) colors suggest that Mrk 1239 is more similar to type-2 active galactic nuclei like NGC 1068. Detailed analysis of… ▽ More We present new spectrophotometric and spectropolarimetric observations of Mrk 1239, one of the 8 prototypes that defines type-1 narrow-line Seyfert galaxies (NLS1s). Unlike the other typical NLS1s though, a high degree of polarization ($P\sim$5.6%) and red optical-IR ($g-W_4$ = 12.35) colors suggest that Mrk 1239 is more similar to type-2 active galactic nuclei like NGC 1068. Detailed analysis of spectral energy distribution in the UV-optical-IR yields two components from the nucleus: a direct and transmitted component that is heavily obscured ($E_{B-V} \approx 1.6$), and another indirect and scattered one with mild extinction ($E_{B-V} \sim$ 0.5). Such a two-light-paths scenario is also found in previous reports based on the X-ray data. Comparison of emission lines and the detection of He\,{\footnotesize I}*$λ$10830 BAL at [-3000,-1000] km s$^{-1}$ indicates that the obscuring clouds are at physical scale between the sublimation radius and that of the narrow emission line regions. The potential existence of powerful outflows is found as both the obscurer and scatterer are outflowing. Similar to many other type-2s, jet-like structure in the radio band is found in Mrk 1239, perpendicular to the polarization angle, suggesting polar scattering. We argue that Mrk 1239 is very probably a type-2 counterpart of NLS1s. The identification of 1 out of 8 prototype NLS1s as a type-2 counterpart implies that there can be a substantial amount of analogs of Mrk 1239 misidentified as type-1s in the optical band. Properties of these misidentified objects are going to be explored in our future works. △ Less

Submitted 28 May, 2021; originally announced May 2021.

Comments: 14 pages, 13 figures

Journal ref: The Astrophysical Journal, Volume 912 (2021), Issue 2, id.118, 13 pp

The electron-capture origin of supernova 2018zd

Authors: Daichi Hiramatsu, D. Andrew Howell, Schuyler D. Van Dyk, Jared A. Goldberg, Keiichi Maeda, Takashi J. Moriya, Nozomu Tominaga, Ken'ichi Nomoto, Griffin Hosseinzadeh, Iair Arcavi, Curtis McCully, Jamison Burke, K. Azalee Bostroem, Stefano Valenti, Yize Dong, Peter J. Brown, Jennifer E. Andrews, Christopher Bilinski, G. Grant Williams, Paul S. Smith, Nathan Smith, David J. Sand, Gagandeep S. Anand, Chengyuan Xu, Alexei V. Filippenko , et al. (5 additional authors not shown)

Abstract: In the transitional mass range ($\sim$ 8-10 solar masses) between white dwarf formation and iron core-collapse supernovae, stars are expected to produce an electron-capture supernova. Theoretically, these progenitors are thought to be super-asymptotic giant branch stars with a degenerate O+Ne+Mg core, and electron capture onto Ne and Mg nuclei should initiate core collapse. However, no supernovae… ▽ More In the transitional mass range ($\sim$ 8-10 solar masses) between white dwarf formation and iron core-collapse supernovae, stars are expected to produce an electron-capture supernova. Theoretically, these progenitors are thought to be super-asymptotic giant branch stars with a degenerate O+Ne+Mg core, and electron capture onto Ne and Mg nuclei should initiate core collapse. However, no supernovae have unequivocally been identified from an electron-capture origin, partly because of uncertainty in theoretical predictions. Here we present six indicators of electron-capture supernovae and show that supernova 2018zd is the only known supernova having strong evidence for or consistent with all six: progenitor identification, circumstellar material, chemical composition, explosion energy, light curve, and nucleosynthesis. For supernova 2018zd, we infer a super-asymptotic giant branch progenitor based on the faint candidate in the pre-explosion images and the chemically-enriched circumstellar material revealed by the early ultraviolet colours and flash spectroscopy. The light-curve morphology and nebular emission lines can be explained with the low explosion energy and neutron-rich nucleosynthesis produced in an electron-capture supernova. This identification provides insights into the complex stellar evolution, supernova physics, cosmic nucleosynthesis, and remnant populations in the transitional mass range. △ Less

Submitted 28 June, 2021; v1 submitted 4 November, 2020; originally announced November 2020.

Comments: Author version of the published letter in Nature Astronomy, 28 June 2021

Mid-infrared Studies of HD 113766 and HD 172555: Assessing Variability in the Terrestrial Zone of Young Exoplanetary Systems

Authors: Kate Y. L. Su, George H. Rieke, Carl Melis, Alan P. Jackson, Paul S. Smith, Huan Y. A. Meng, Andras Gaspar

Abstract: We present multi-epoch infrared photometry and spectroscopy obtained with warm Spitzer, Subaru and SOFIA to assess variability for the young ($\sim$20 Myr) and dusty debris systems around HD 172555 and HD 113766A. No variations (within 0.5%) were found for the former at either 3.6 or 4.5 $μ$m, while significant non-periodic variations (peak-to-peak of $\sim$10-15% relative to the primary star) wer… ▽ More We present multi-epoch infrared photometry and spectroscopy obtained with warm Spitzer, Subaru and SOFIA to assess variability for the young ($\sim$20 Myr) and dusty debris systems around HD 172555 and HD 113766A. No variations (within 0.5%) were found for the former at either 3.6 or 4.5 $μ$m, while significant non-periodic variations (peak-to-peak of $\sim$10-15% relative to the primary star) were detected for the latter. Relative to the Spitzer IRS spectra taken in 2004, multi-epoch mid-infrared spectra reveal no change in either the shape of the prominent 10 $μ$m solid-state features or the overall flux levels (no more than 20%) for both systems, corroborating that the population of sub-$μ$m-sized grains that produce the pronounced solid-state features is stable over a decadal timescale. We suggest that these sub-$μ$m-sized grains were initially generated in an optically thick clump of debris of mm-sized vapor condensates resulting from a recent violent impact between large asteroidal or planetary bodies. Because of the shielding from the stellar photons provided by this clump, intense collisions led to an over-production of fine grains that would otherwise be ejected from the system by radiation pressure. As the clump is sheared by its orbital motion and becomes optically thin, a population of very fine grains could remain in stable orbits until Poynting-Robertson drag slowly spirals them into the star. We further suggest that the 3-5 $μ$m disk variation around HD 113766A is consistent with a clump/arc of such fine grains on a modestly eccentric orbit in its terrestrial zone. △ Less

Submitted 11 June, 2020; originally announced June 2020.

Comments: to be published in ApJ

The Great Markarian 421 Flare of February 2010: Multiwavelength variability and correlation studies

Authors: A. U. Abeysekara, W. Benbow, R. Bird, A. Brill, R. Brose, M. Buchovecky, J. H. Buckley, J. L. Christiansen, A. J. Chromey, M. K. Daniel, J. Dumm, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, N. Galante, A. Gent, G. H. Gillanders, C. Giuri, O. Gueta, T. Hassan, O. Hervet, J. Holder, G. Hughes , et al. (234 additional authors not shown)

Abstract: We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a co… ▽ More We report on variability and correlation studies using multiwavelength observations of the blazar Mrk 421 during the month of February, 2010 when an extraordinary flare reaching a level of $\sim$27~Crab Units above 1~TeV was measured in very-high-energy (VHE) $γ$-rays with the VERITAS observatory. This is the highest flux state for Mrk 421 ever observed in VHE $γ$-rays. Data are analyzed from a coordinated campaign across multiple instruments including VHE $γ$-ray (VERITAS, MAGIC), high-energy (HE) $γ$-ray (Fermi-LAT), X-ray (Swift}, RXTE, MAXI), optical (including the GASP-WEBT collaboration and polarization data) and radio (Metsähovi, OVRO, UMRAO). Light curves are produced spanning multiple days before and after the peak of the VHE flare, including over several flare `decline' epochs. The main flare statistics allow 2-minute time bins to be constructed in both the VHE and optical bands enabling a cross-correlation analysis that shows evidence for an optical lag of $\sim$25-55 minutes, the first time-lagged correlation between these bands reported on such short timescales. Limits on the Doppler factor ($δ\gtrsim 33$) and the size of the emission region ($ δ^{-1}R_B \lesssim 3.8\times 10^{13}\,\,\mbox{cm}$) are obtained from the fast variability observed by VERITAS during the main flare. Analysis of 10-minute-binned VHE and X-ray data over the decline epochs shows an extraordinary range of behavior in the flux-flux relationship: from linear to quadratic to lack of correlation to anti-correlation. Taken together, these detailed observations of an unprecedented flare seen in Mrk 421 are difficult to explain by the classic single-zone synchrotron self-Compton model. △ Less

Submitted 10 February, 2020; originally announced February 2020.

Comments: 41 pages including 3 appendices, 13 figures; version accepted to Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 890, Number 2 (2020)

Unravelling the complex behavior of Mrk 421 with simultaneous X-ray and VHE observations during an extreme flaring activity in April 2013

Authors: MAGIC collaboration, V. A. Acciari, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, D. Baack, A. Babic, B. Banerjee, U. Barres de Almeida, J. A. Barrio, J. Becerra Gonzalez, W. Bednarek, L. Bellizzi, E. Bernardini, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, G. Bonnoli, Z. Bosnjak, G. Busetto, R. Carosi, G. Ceribella , et al. (215 additional authors not shown)

Abstract: We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on… ▽ More We report on a multi-band variability and correlation study of the TeV blazar Mrk 421 during an exceptional flaring activity observed from 2013 April 11 to 2013 April 19. The study uses, among others, data from GASP-WEBT, Swift, NuSTAR, Fermi-LAT, VERITAS, and MAGIC. The large blazar activity, and the 43 hours of simultaneous NuSTAR and MAGIC/VERITAS observations, permitted variability studies on 15 minute time bins, and over three X-ray bands (3-7 keV, 7-30 keV and 30-80 keV) and three very-high-energy (>0.1 TeV, hereafter VHE) gamma-ray bands (0.2-0.4 TeV, 0.4-0.8 TeV and >0.8 TeV). We detected substantial flux variations on multi-hour and sub-hour timescales in all the X-ray and VHE gamma-ray bands. The characteristics of the sub-hour flux variations are essentially energy-independent, while the multi-hour flux variations can have a strong dependence on the energy of the X-ray and the VHE gamma rays. The three VHE bands and the three X-ray bands are positively correlated with no time-lag, but the strength and the characteristics of the correlation changes substantially over time and across energy bands. Our findings favour multi-zone scenarios for explaining the achromatic/chromatic variability of the fast/slow components of the light curves, as well as the changes in the flux-flux correlation on day-long timescales. We interpret these results within a magnetic reconnection scenario, where the multi-hour flux variations are dominated by the combined emission from various plasmoids of different sizes and velocities, while the sub-hour flux variations are dominated by the emission from a single small plasmoid moving across the magnetic reconnection layer. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 46 pages, 17 figures, 9 tables, accepted for publication in ApJ Supplements

Multiwavelength behaviour of the blazar 3C279: decade-long study from $γ$-ray to radio

Authors: V. M. Larionov, S. G. Jorstad, A. P. Marscher, M. Villata, C. M. Raiteri, P. S. Smith, I. Agudo, S. S. Savchenko, D. A. Morozova, J. A. Acosta-Pulido, M. F. Aller, H. D. Aller, T. S. Andreeva, A. A. Arkharov, R. Bachev, G. Bonnoli, G. A. Borman, V. Bozhilov, P. Calcidese, M. I. Carnerero, D. Carosati, C. Casadio, W. -P. Chen, G. Damljanovic, A. V. Dementyev , et al. (62 additional authors not shown)

Abstract: We report the results of decade-long (2008-2018) $γ$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $γ$-ray-optical flux-flux relation changes with activi… ▽ More We report the results of decade-long (2008-2018) $γ$-ray to 1 GHz radio monitoring of the blazar 3C 279, including GASP/WEBT, $\it{Fermi}$ and $\it{Swift}$ data, as well as polarimetric and spectroscopic data. The X-ray and $γ$-ray light curves correlate well, with no delay > 3 hours, implying general co-spatiality of the emission regions. The $γ$-ray-optical flux-flux relation changes with activity state, ranging from a linear to a more complex dependence. The behaviour of the Stokes parameters at optical and radio wavelengths, including 43 GHz VLBA images, supports either a predominantly helical magnetic field or motion of the radiating plasma along a spiral path. Apparent speeds of emission knots range from 10 to 37c, with the highest values requiring bulk Lorentz factors close to those needed to explain $γ$-ray variability on very short time scales. The Mg II emission line flux in the `blue' and `red' wings correlates with the optical synchrotron continuum flux density, possibly providing a variable source of seed photons for inverse Compton scattering. In the radio bands we find progressive delays of the most prominent light curve maxima with decreasing frequency, as expected from the frequency dependence of the $τ=1$ surface of synchrotron self-absorption. The global maximum in the 86 GHz light curve becomes less prominent at lower frequencies, while a local maximum, appearing in 2014, strengthens toward decreasing frequencies, becoming pronounced at $\sim5$ GHz. These tendencies suggest different Doppler boosting of stratified radio-emitting zones in the jet. △ Less

Submitted 17 January, 2020; originally announced January 2020.

Comments: 21 pages, 18 figures, accepted for publication in MNRAS

Journal ref: MNRAS 492 (2020) 3829

Quasi-periodic behaviour in the optical and γ-ray light curves of blazars 3C 66A and B2 1633+38

Authors: J. Otero-Santos, J. A. Acosta-Pulido, J. Becerra González, C. M. Raiteri, V. M. Larionov, P. Peñil, P. S. Smith, C. Ballester Niebla, G. A. Borman, M. I. Carnerero, N. Castro Segura, T. S. Grishina, E. N. Kopatskaya, E. G. Larionova, D. A. Morozova, A. A. Nikiforova, S. S. Savchenko, Yu. V. Troitskaya, I. S. Troitsky, A. A. Vasilyev, M. Villata

Abstract: We report on quasi-periodic variability found in two blazars included in the Steward Observatory Blazar Monitoring data sample: the BL Lac object 3C 66A and the Flat Spectrum Radio Quasar B2 1633+38. We collect optical photometric and polarimetric data in V and R bands of these sources from different observatories: St. Petersburg University, Crimean Astrophysical Observatory, WEBT-GASP, Catalina R… ▽ More We report on quasi-periodic variability found in two blazars included in the Steward Observatory Blazar Monitoring data sample: the BL Lac object 3C 66A and the Flat Spectrum Radio Quasar B2 1633+38. We collect optical photometric and polarimetric data in V and R bands of these sources from different observatories: St. Petersburg University, Crimean Astrophysical Observatory, WEBT-GASP, Catalina Real-Time Transient Survey, Steward Observatory, STELLA Robotic Observatory and Katzman Automatic Imaging Telescope. In addition, an analysis of the $γ$-ray light curves from $\textit{Fermi}$-LAT is included. Three methods are used to search for any periodic behaviour in the data: the Z-transform Discrete Correlation Function, the Lomb-Scargle periodogram and the Weighted Wavelet Z-transform. We find evidences of possible quasi-periodic variability in the optical photometric data of both sources with periods of $\sim$3 years for 3C 66A and $\sim$1.9 years for B2 1633+38, with significances between 3$σ$ and 5$σ$. Only B2 1633+38 shows evidence of this behaviour in the optical polarized data set at a confidence level of 2$σ$-4$σ$. This is the first reported evidence of quasi-periodic behaviour in the optical light curve of B2 1633+38. Also a hint of quasi-periodic behaviour is found in the $γ$-ray light curve of B2 1633+38 with a confidence level $\geqslant$2$σ$, while no periodicity is observed for 3C 66A in this energy range. We propose different jet emission models that could explain the quasi-periodic variability and the differences found between these two sources. △ Less

Submitted 5 February, 2020; v1 submitted 16 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS

H.E.S.S. detection of very-high-energy gamma-ray emission from the quasar PKS 0736+017

Authors: H. E. S. S. Collaboration, :, H. Abdalla, R. Adam, F. Aharonian, F. Ait Benkhali, E. O. Angüner, M. Arakawa, C. Arcaro, C. Armand, H. Ashkar, M. Backes, V. Barbosa Martins, M. Barnard, Y. Becherini, D. Berge, K. Bernlöhr, R. Blackwell, M. Böttcher, C. Boisson, J. Bolmont, S. Bonnefoy, J. Bregeon, M. Breuhaus, F. Brun , et al. (203 additional authors not shown)

Abstract: Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes… ▽ More Flat-spectrum radio-quasars (FSRQs) are rarely detected at very-high-energies (VHE; E>100 GeV) due to their low-frequency-peaked SEDs. At present, only 6 FSRQs are known to emit VHE photons, representing only 7% of the VHE extragalactic catalog. Following the detection of MeV-GeV gamma-ray flaring activity from the FSRQ PKS 0736+017 (z=0.189) with Fermi, the H.E.S.S. array of Cherenkov telescopes triggered ToO observations on February 18, 2015, with the goal of studying the gamma-ray emission in the VHE band. H.E.S.S. ToO observations were carried out during the nights of February 18, 19, 21, and 24, 2015. Together with Fermi-LAT, the multi-wavelength coverage of the flare includes Swift observations in soft-X-rays and optical/UV, and optical monitoring (photometry and spectro-polarimetry) by the Steward Observatory, the ATOM, the KAIT and the ASAS-SN telescope. VHE emission from PKS 0736+017 was detected with H.E.S.S. during the night of February 19, 2015, only. Fermi data indicate the presence of a gamma-ray flare, peaking at the time of the H.E.S.S. detection, with a flux doubling time-scale of around six hours. The gamma-ray flare was accompanied by at least a 1 mag brightening of the non-thermal optical continuum. No simultaneous observations at longer wavelengths are available for the night of the H.E.S.S. detection. The gamma-ray observations with H.E.S.S. and Fermi are used to put constraints on the location of the gamma-ray emitting region during the flare: it is constrained to be just outside the radius of the broad-line-region with a bulk Lorentz factor $\simeq 20$, or at the level of the radius of the dusty torus with Gamma > 60. PKS 0736+017 is the seventh FSRQ known to emit VHE photons and, at z=0.189, is the nearest so far. The location of the gamma-ray emitting region during the flare can be tightly constrained thanks to opacity, variability, and collimation arguments. △ Less

Submitted 12 November, 2019; originally announced November 2019.

Comments: In press in A&A. 13 pages, 5 figures. Abstract abridged

Journal ref: A&A 633, A162 (2020)

Investigating the multiwavelength behaviour of the flat spectrum radio quasar CTA 102 during 2013-2017

Authors: F. D'Ammando, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, I. Agudo, A. A. Arkharov, R. Bachev, G. V. Baida, E. Benitez, G. A. Borman, W. Boschin, V. Bozhilov, M. S. Butuzova, P. Calcidese, M. I. Carnerero, D. Carosati, C. Casadio, N. Castro-Segura, W. -P. Chen, G. Damljanovic, A. Di Paola, J. Echevarria, N. V. Efimova, Sh. A. Ehgamberdiev, C. Espinosa , et al. (72 additional authors not shown)

Abstract: We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and… ▽ More We present a multiwavelength study of the flat-spectrum radio quasar CTA 102 during 2013-2017. We use radio-to-optical data obtained by the Whole Earth Blazar Telescope, 15 GHz data from the Owens Valley Radio Observatory, 91 and 103 GHz data from the Atacama Large Millimeter Array, near-infrared data from the Rapid Eye Monitor telescope, as well as data from the Swift (optical-UV and X-rays) and Fermi ($γ$ rays) satellites to study flux and spectral variability and the correlation between flux changes at different wavelengths. Unprecedented $γ$-ray flaring activity was observed during 2016 November-2017 February, with four major outbursts. A peak flux of (2158 $\pm$ 63)$\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, corresponding to a luminosity of (2.2 $\pm$ 0.1)$\times$10$^{50}$ erg s$^{-1}$, was reached on 2016 December 28. These four $γ$-ray outbursts have corresponding events in the near-infrared, optical, and UV bands, with the peaks observed at the same time. A general agreement between X-ray and $γ$-ray activity is found. The $γ$-ray flux variations show a general, strong correlation with the optical ones with no time lag between the two bands and a comparable variability amplitude. This $γ$-ray/optical relationship is in agreement with the geometrical model that has successfully explained the low-energy flux and spectral behaviour, suggesting that the long-term flux variations are mainly due to changes in the Doppler factor produced by variations of the viewing angle of the emitting regions. The difference in behaviour between radio and higher energy emission would be ascribed to different viewing angles of the jet regions producing their emission. △ Less

Submitted 11 November, 2019; v1 submitted 8 October, 2019; originally announced October 2019.

Comments: 24 pages, 18 figures, and 4 tables. Accepted for publication in Monthly Notices of the Royal Astronomical Society

Mid-IR Variability and Dust Reverberation Mapping of Low-$z$ Quasars. I. Data, Methods and Basic Results

Authors: Jianwei Lyu, George H. Rieke, Paul S. Smith

Abstract: The continued operation of the Wide-field Infrared Survey Explorer (WISE) combined with several ground-based optical transient surveys (e.g., CRTS, ASAS-SN and PTF) offer an unprecedented opportunity to explore the dust structures in luminous AGNs. We use these data for a mid-IR dust reverberation mapping (RM) study of 87 archetypal Palomar-Green quasars at $z\lesssim0.5$. To cope with various con… ▽ More The continued operation of the Wide-field Infrared Survey Explorer (WISE) combined with several ground-based optical transient surveys (e.g., CRTS, ASAS-SN and PTF) offer an unprecedented opportunity to explore the dust structures in luminous AGNs. We use these data for a mid-IR dust reverberation mapping (RM) study of 87 archetypal Palomar-Green quasars at $z\lesssim0.5$. To cope with various contaminations of the photometry data and the sparse time sampling of the light curves, procedures to combine these datasets and retrieve the dust RM signals have been developed. We find $\sim$70% of the sample (with a completeness correction, up to 95%) have convincing mid-IR time-lags in the WISE W1 ($\sim3.4~μm$) and W2 ($\sim4.5~μm$) bands and they are proportional to the square root of the AGN luminosity. Combined with previous K-band ($\sim2.2~μm$) RM results in the literature, the inferred dust emission size ratios are $R_{K}:R_{W1}:R_{W2}=0.6:1:1.2$. Under simple assumptions, we put preliminary constraints on the projected dust surface density at these bands and reveal the possibly different torus structures among hot-dust-deficient, warm-dust-deficient and normal quasars from the reverberation signals. With multi-epoch Spitzer data and later WISE photometry, we also explore AGN IR variability at 10--24 $μm$ over a 5-year time-scale. Except for blazars and flat-spectrum radio sources, the majority of AGNs have typical variation amplitudes at 24 $μm$ no more than 10% of that in the W1 band, indicating the dust reverberation signals damp out quickly at longer wavelengths. In particular, steep-spectrum radio quasars also lack strong 24 $μm$ variability, consistent with the unification picture of radio-loud AGN. △ Less

Submitted 17 November, 2019; v1 submitted 24 September, 2019; originally announced September 2019.

Comments: ApJ accepted (35 pages, 21 figures), revised to match the published paper

Journal ref: ApJ, 886, 33 (2019)

Discovery and Rapid Follow-up Observations of the Unusual Type II SN 2018ivc in NGC 1068

Authors: K. A. Bostroem, S. Valenti, D. J. Sand, J. E. Andrews, S. D. Van Dyk, L. Galbany, D. Pooley, R. C. Amaro, N. Smith, S. Yang, G. C. Anupama, I. Arcavi, E. Baron, P. J. Brown, J. Burke, R. Cartier, D. Hiramatsu, Y. Dong, E. Egami, S. Ertel, A. V. Filippenko, O. D. Fox, J. Haislip, G. Hosseinzadeh, D. A. Howell , et al. (22 additional authors not shown)

Abstract: We present the discovery and high-cadence follow-up observations of SN 2018ivc, an unusual Type II supernova that exploded in NGC 1068 (D=10.1 Mpc). The light curve of SN 2018ivc declines piecewise-linearly, changing slope frequently, with four clear slope changes in the first 30 days of evolution. This rapidly changing light curve indicates that interaction between the circumstellar material and… ▽ More We present the discovery and high-cadence follow-up observations of SN 2018ivc, an unusual Type II supernova that exploded in NGC 1068 (D=10.1 Mpc). The light curve of SN 2018ivc declines piecewise-linearly, changing slope frequently, with four clear slope changes in the first 30 days of evolution. This rapidly changing light curve indicates that interaction between the circumstellar material and ejecta plays a significant role in the evolution. Circumstellar interaction is further supported by a strong X-ray detection. The spectra are rapidly evolving and dominated by hydrogen, helium, and calcium emission lines. We identify a rare high-velocity emission-line feature blueshifted at ~7800 km/s (in Ha, Hb, Pb, Pg, HeI, CaII), which is visible from day 18 until at least day 78 and could be evidence of an asymmetric progenitor or explosion. From the overall similarity between SN 2018ivc and SN 1996al, the \Ha{} equivalent width of its parent HII region, and constraints from pre-explosion archival Hubble Space Telescope images, we find that the progenitor of SN 2018ivc could be as massive as 52 Msun but is more likely <12 Msun. SN 2018ivc demonstrates the importance of the early discovery and rapid follow-up observations of nearby supernovae to study the physics and progenitors of these cosmic explosions. △ Less

Submitted 27 May, 2020; v1 submitted 16 September, 2019; originally announced September 2019.

Comments: Accepted by ApJ. Revised version includes more extensive progenitor analysis

Journal ref: The Astrophysical Journal, 2020, Volume 895, Number 1

Searches After Gravitational-waves Using ARizona Observatories (SAGUARO): System Overview and First Results from Advanced LIGO/Virgo's Third Observing Run

Authors: M. J. Lundquist, K. Paterson, W. Fong, D. J. Sand, J. E. Andrews, I. Shivaei, P. N. Daly, S. Valenti, S. Yang, E. Christensen, A. R. Gibbs, F. Shelly, S. Wyatt, O. Kuhn, R. C. Amaro, I. Arcavi, P. Behroozi, N. Butler, L. Chomiuk, A. Corsi, M. R. Drout, E. Egami, X. Fan, R. J. Foley, B. Frye , et al. (31 additional authors not shown)

Abstract: We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pi… ▽ More We present Searches After Gravitational-waves Using ARizona Observatories (SAGUARO), a comprehensive effort dedicated to the discovery and characterization of optical counterparts to gravitational wave (GW) events. SAGUARO utilizes ground-based facilities ranging from 1.5m to 10m in diameter, located primarily in the Northern Hemisphere. We provide an overview of SAGUARO's telescopic resources, pipeline for transient detection, and database for candidate visualization. We describe SAGUARO's discovery component, which utilizes the $5$~deg$^2$ field-of-view optical imager on the Mt. Lemmon 1.5m telescope, reaching limits of $\approx 21.3$~AB mag while rapidly tiling large areas. We also describe the follow-up component of SAGUARO, used for rapid vetting and monitoring of optical candidates. With the onset of Advanced LIGO/Virgo's third observing run, we present results from the first three SAGUARO searches following the GW events S190408an, S190425z and S190426c, which serve as a valuable proof-of-concept of SAGUARO. We triggered and searched 15, 60 and 60 deg$^{2}$ respectively, 17.6, 1.4 and 41.8 hrs after the initial GW alerts. We covered 7.8, 3.0 and 5.1\% of the total probability within the GW event localizations, reaching 3$σ$ limits of 19.8, 21.3 and 20.8 AB mag, respectively. Although no viable counterparts associated with these events were found, we recovered 6 known transients and ruled out 5 potential candidates. We also present Large Binocular Telescope spectroscopy of PS19eq/SN2019ebq, a promising kilonova candidate that was later determined to be a supernova. With the ability to tile large areas and conduct detailed follow-up, SAGUARO represents a significant addition to GW counterpart searches. △ Less

Submitted 17 July, 2019; v1 submitted 14 June, 2019; originally announced June 2019.

Comments: 16 pages, 7 figures, 1 table. Accepted to ApJL

The June 2016 Optical and Gamma-Ray Outburst and Optical Micro-Variability of the Blazar 3C454.3

Authors: Zachary R. Weaver, Thomas J. Balonek, Svetlana G. Jorstad, Alan P. Marscher, Valeri M. Larionov, Paul S. Smith, Samantha J. Boni, George A. Borman, K. J. Chapman, Leah G. Jenks, Evgenia N. Kopatskaya, Elena G. Larionova, Daria A. Morozova, Anna A. Nikiforova, Alina Sabyr, Sergey S. Savchenko, Ryan W. Stahlin, Yulia V. Troitskaya, Ivan S. Troitsky, Saiyang Zhang

Abstract: The quasar 3C454.3 underwent a uniquely-structured multi-frequency outburst in June 2016. The blazar was observed in the optical $R$ band by several ground-based telescopes in photometric and polarimetric modes, at $γ$-ray frequencies by the \emph{Fermi}\ Large Area Telescope, and at 43 GHz with the Very Long Baseline Array. The maximum flux density was observed on 2016 June 24 at both optical and… ▽ More The quasar 3C454.3 underwent a uniquely-structured multi-frequency outburst in June 2016. The blazar was observed in the optical $R$ band by several ground-based telescopes in photometric and polarimetric modes, at $γ$-ray frequencies by the \emph{Fermi}\ Large Area Telescope, and at 43 GHz with the Very Long Baseline Array. The maximum flux density was observed on 2016 June 24 at both optical and $γ$-ray frequencies, reaching $S^\mathrm{max}_\mathrm{opt}=18.91\pm0.08$ mJy and $S_γ^\mathrm{max} =22.20\pm0.18\times10^{-6}$ ph cm$^{-2}$ s$^{-1}$, respectively. The June 2016 outburst possessed a precipitous decay at both $γ$-ray and optical frequencies, with the source decreasing in flux density by a factor of 4 over a 24-hour period in $R$ band. Intraday variability was observed throughout the outburst, with flux density changes between 1 and 5 mJy over the course of a night. The precipitous decay featured statistically significant quasi-periodic micro-variability oscillations with an amplitude of $\sim 2$-$3\%$ about the mean trend and a characteristic period of 36 minutes. The optical degree of polarization jumped from $\sim3\%$ to nearly 20\% during the outburst, while the position angle varied by $\sim120\degr$. A knot was ejected from the 43 GHz core on 2016 Feb 25, moving at an apparent speed $v_\mathrm{app}=20.3c\pm0.8c$. From the observed minimum timescale of variability $τ_\mathrm{opt}^\mathrm{min}\approx2$ hr and derived Doppler factor $δ=22.6$, we find a size of the emission region $r\lesssim2.6\times10^{15}$ cm. If the quasi-periodic micro-variability oscillations are caused by periodic variations of the Doppler factor of emission from a turbulent vortex, we derive a rotational speed of the vortex $\sim0.2c$. △ Less

Submitted 11 March, 2019; originally announced March 2019.

Comments: 19 pages, 13 figures, 3 tables, accepted to the Astrophysical Journal 2019 March 9

A Deeply Buried Narrow-Line Seyfert 1 Nucleus Uncovered in Scattered Light

Authors: Xiang Pan, Honglin Lu, S. Komossa, Dawei Xu, Weimin Yuan, Luming Sun, Paul S. Smith, Shaohua Zhang, Peng Jiang, Chenwei Yang, Wenjuan Liu, Ning Jiang, Y. E. Rashed, A. Eckart, Jens Dierkes, Hongyan Zhou

Abstract: We present spectropolarimetric and spectrophotometric observations of the peculiar active galactic nucleus (AGN) SDSS J120300.19+162443.7 (hereafter J1203+1624) at $z=0.1656$. Its optical total flux spectra clearly show broad emission lines (BELs) in H$α$ and H$β$. After removal of narrow emission lines (NELs), the full width at half maximum of the lorentzian BEL is $FWHM_{\rm BEL} \approx 1,270$… ▽ More We present spectropolarimetric and spectrophotometric observations of the peculiar active galactic nucleus (AGN) SDSS J120300.19+162443.7 (hereafter J1203+1624) at $z=0.1656$. Its optical total flux spectra clearly show broad emission lines (BELs) in H$α$ and H$β$. After removal of narrow emission lines (NELs), the full width at half maximum of the lorentzian BEL is $FWHM_{\rm BEL} \approx 1,270$ km s$^{-1}$, fulfilling the conventional definition of a narrow line Seyfert 1 (NLS1) galaxy. However, its NELs are unprecedentedly strong when compared to type 1AGNs. This, together with its large MIR excess ($g-W_4 = 13.172$ mag), implies that the observer and the narrow emission line region (NELR) might see a different ionization continuum. Our optical spectropolarimetry confirms its type 2 nature by detecting a polarized blue continuum and Balmer BELs ($FWHM_{\rm Polarized~BEL} \approx 1,183$ km s$^{-1}$), with a high polarization degree of $> 20$% in the blue wing. The mass and Eddington rate of the central black hole is estimated based on both transmitted and scattered AGN radiation, which is $M_{\bullet} < 2.9 \times 10^7 M_\odot$ and $L_{\rm bol}/L_{\rm Edd} > 1.5$. Severe extinction of the AGN emission also enables a clear view of the compact host galaxy. Discovery of J1203+1624 suggests that NLS1 follows AGN unfication scheme, studying of its analogs could blaze a new trail for exploring the connection between black hole growth and star formation in the host galaxy. Interesting features of J1203+1624 like the peculiar NELs, inflowing scattering clouds within sublimation radius are worth detailed follow-ups in the future. △ Less

Submitted 11 October, 2019; v1 submitted 15 November, 2018; originally announced November 2018.

Comments: 9 pages, 6 figures

Journal ref: 2019,ApJ,850,75

Kinematics of the Broad-line Region of 3C 273 from a Ten-year Reverberation Mapping Campaign

Authors: Zhi-Xiang Zhang, Pu Du, Paul S. Smith, Yulin Zhao, Chen Hu, Ming Xiao, Yan-Rong Li, Ying-Ke Huang, Kai Wang, Jin-Ming Bai, Luis C. Ho, Jian-Min Wang

Abstract: Despite many decades of study, the kinematics of the broad-line region of 3C~273 are still poorly understood. We report a new, high signal-to-noise, reverberation mapping campaign carried out from November 2008 to March 2018 that allows the determination of time lags between emission lines and the variable continuum with high precision. The time lag of variations in H$β$ relative to those of the 5… ▽ More Despite many decades of study, the kinematics of the broad-line region of 3C~273 are still poorly understood. We report a new, high signal-to-noise, reverberation mapping campaign carried out from November 2008 to March 2018 that allows the determination of time lags between emission lines and the variable continuum with high precision. The time lag of variations in H$β$ relative to those of the 5100 Angstrom continuum is $146.8_{-12.1}^{+8.3}$ days in the rest frame, which agrees very well with the Paschen-$α$ region measured by the GRAVITY at The Very Large Telescope Interferometer. The time lag of the H$γ$ emission line is found to be nearly the same as for H$β$. The lag of the Fe II emission is $322.0_{-57.9}^{+55.5}$ days, longer by a factor of $\sim$2 than that of the Balmer lines. The velocity-resolved lag measurements of the H$β$ line show a complex structure which can be possibly explained by a rotation-dominated disk with some inflowing radial velocity in the H$β$-emitting region. Taking the virial factor of $f_{\rm BLR} = 1.3$, we derive a BH mass of $M_{\bullet} = 4.1_{-0.4}^{+0.3} \times 10^8 M_{\odot}$ and an accretion rate of $9.3\,L_{\rm Edd}\,c^{-2}$ from the H$β$ line. The decomposition of its $HST$ images yields a host stellar mass of $M_* = 10^{11.3 \pm 0.7} M_\odot$, and a ratio of $M_{\bullet}/M_*\approx 2.0\times 10^{-3}$ in agreement with the Magorrian relation. In the near future, it is expected to compare the geometrically-thick BLR discovered by the GRAVITY in 3C 273 with its spatially-resolved torus in order to understand the potential connection between the BLR and the torus. △ Less

Submitted 12 March, 2019; v1 submitted 9 November, 2018; originally announced November 2018.

Comments: 17 pages, 12 figures, 6 tables, accepted for publication in The Astrophysical Journal

The extreme HBL behaviour of Markarian 501 during 2012

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, W. Bhattacharyya, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, A. Chatterjee, S. M. Colak, P. Colin, E. Colombo, J. L. Contreras, J. Cortina, S. Covino , et al. (254 additional authors not shown)

Abstract: A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0… ▽ More A multiwavelength campaign was organized to take place between March and July of 2012. Excellent temporal coverage was obtained with more than 25 instruments, including the MAGIC, FACT and VERITAS Cherenkov telescopes, the instruments on board the Swift and Fermi spacecraft, and the telescopes operated by the GASP-WEBT collaboration. Mrk 501 showed a very high energy (VHE) gamma-ray flux above 0.2 TeV of $\sim$0.5 times the Crab Nebula flux (CU) for most of the campaign. The highest activity occurred on 2012 June 9, when the VHE flux was $\sim$3 CU, and the peak of the high-energy spectral component was found to be at $\sim$2 TeV. This study reports very hard X-ray spectra, and the hardest VHE spectra measured to date for Mrk 501. The fractional variability was found to increase with energy, with the highest variability occurring at VHE, and a significant correlation between the X-ray and VHE bands. The unprecedentedly hard X-ray and VHE spectra measured imply that their low- and high-energy components peaked above 5 keV and 0.5 TeV, respectively, during a large fraction of the observing campaign, and hence that Mrk 501 behaved like an extreme high-frequency- peaked blazar (EHBL) throughout the 2012 observing season. This suggests that being an EHBL may not be a permanent characteristic of a blazar, but rather a state which may change over time. The one-zone synchrotron self-Compton (SSC) scenario can successfully describe the segments of the SED where most energy is emitted, with a significant correlation between the electron energy density and the VHE gamma-ray activity, suggesting that most of the variability may be explained by the injection of high-energy electrons. The one-zone SSC scenario used reproduces the behaviour seen between the measured X-ray and VHE gamma-ray fluxes, and predicts that the correlation becomes stronger with increasing energy of the X-rays. △ Less

Submitted 14 August, 2018; v1 submitted 13 August, 2018; originally announced August 2018.

Comments: 25 pages, 9 figures, accepted for publication in A&A. Corresponding authors: Gareth Hughes ([email protected]), David Paneque ([email protected]), Amit Shukla ([email protected])

Journal ref: A&A 620, A181 (2018)

Multiwavelength observations of the blazar BL Lacertae: a new fast TeV gamma-ray flare

Authors: A. U. Abeysekara, W. Benbow, R. Bird, T. Brantseg, R. Brose, M. Buchovecky, J. H. Buckley, V. Bugaev, M. P. Connolly, W. Cui, M. K. Daniel, A. Falcone, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, G. H. Gillanders, I. Gunawardhana, M. Hütten, D. Hanna, O. Hervet, J. Holder, G. Hughes, T. B. Humensky, C. A. Johnson , et al. (52 additional authors not shown)

Abstract: Combined with very-long-baseline interferometry measurements, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL~Lac objects and flat-spectrum radio quasars. We report on a fast TeV… ▽ More Combined with very-long-baseline interferometry measurements, the observations of fast TeV gamma-ray flares probe the structure and emission mechanism of blazar jets. However, only a handful of such flares have been detected to date, and only within the last few years have these flares been observed from lower-frequency-peaked BL~Lac objects and flat-spectrum radio quasars. We report on a fast TeV gamma-ray flare from the blazar BL~Lacertae observed by VERITAS, with a rise time of $\sim$2.3~hr and a decay time of $\sim$36~min. The peak flux above 200 GeV is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photon} \;\text{m}^{-2}\; \text{s}^{-1}$ measured with a 4-minute-binned light curve, corresponding to $\sim$180\% of the flux which is observed from the Crab Nebula above the same energy threshold. Variability contemporaneous with the TeV gamma-ray flare was observed in GeV gamma-ray, X-ray, and optical flux, as well as in optical and radio polarization. Additionally, a possible moving emission feature with superluminal apparent velocity was identified in VLBA observations at 43 GHz, potentially passing the radio core of the jet around the time of the gamma-ray flare. We discuss the constraints on the size, Lorentz factor, and location of the emitting region of the flare, and the interpretations with several theoretical models which invoke relativistic plasma passing stationary shocks. △ Less

Submitted 27 February, 2018; originally announced February 2018.

Comments: 23 pages, 9 figures

Journal ref: ApJ 856 (2018) 95

Dissecting the long-term emission behaviour of the BL Lac object Mrk 421

Authors: M. I. Carnerero, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, V. M. Larionov, P. S. Smith, F. D'Ammando

Abstract: We report on long-term multi-wavelength monitoring of blazar Mrk 421 by the GASP-WEBT collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007-2015, characterized by several extreme flares. The ratio between the optical, X-ray, and γ-ray fluxes is very variable. The γ-ray flux variations show fair correlation with the optical on… ▽ More We report on long-term multi-wavelength monitoring of blazar Mrk 421 by the GASP-WEBT collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007-2015, characterized by several extreme flares. The ratio between the optical, X-ray, and γ-ray fluxes is very variable. The γ-ray flux variations show fair correlation with the optical ones starting from 2012. We analyse spectropolarimetric data and find wavelength-dependence of the polarisation degree (P), which is compatible with the presence of the host galaxy, and no wavelength-dependence of the electric vector polarisation angle (EVPA). Optical polarimetry shows a lack of simple correlation between P and flux and wide rotations of the EVPA. We build broad-band spectral energy distributions with simultaneous near-infrared and optical data from the GASP-WEBT and UV and X-ray data from the Swift satellite. They show strong variability in both flux and X-ray spectral shape and suggest a shift of the synchrotron peak up to a factor $\sim 50$ in frequency. The interpretation of the flux and spectral variability is compatible with jet models including at least two emitting regions that can change their orientation with respect to the line of sight. △ Less

Submitted 7 September, 2017; originally announced September 2017.

Comments: 18 pages, 27 figures, 1 table, in press for MNRAS

Multiwavelength observations of the blazar BL Lacertae: a new fast TeV $γ$-ray flare

Authors: Q. Feng, S. G. Jorstad, A. P. Marscher, M. L. Lister, Y. Y. Kovalev, A. B. Pushkarev, T. Savolainen, I. Agudo, S. N. Molina, J. L. Gomez, V. M. Larionov, G. A. Borman, A. A. Mokrushina, P. S. Smith

Abstract: Observations of fast TeV $γ$-ray flares from blazars reveal the extreme compactness of emitting regions in blazar jets. Combined with very-long-baseline radio interferometry measurements, they probe the structure and emission mechanism of the jet. We report on a fast TeV $γ$-ray flare from BL Lacertae observed by VERITAS, with a rise time of about 2.3 hours and a decay time of about 36 minutes. Th… ▽ More Observations of fast TeV $γ$-ray flares from blazars reveal the extreme compactness of emitting regions in blazar jets. Combined with very-long-baseline radio interferometry measurements, they probe the structure and emission mechanism of the jet. We report on a fast TeV $γ$-ray flare from BL Lacertae observed by VERITAS, with a rise time of about 2.3 hours and a decay time of about 36 minutes. The peak flux at $>$200 GeV measured with the 4-minute binned light curve is $(4.2 \pm 0.6) \times 10^{-6} \;\text{photons} \;\text{m}^{-2}\, \text{s}^{-1}$, or $\sim$180% the Crab Nebula flux. Variability in GeV $γ$-ray, X-ray, and optical flux, as well as in optical and radio polarization was observed around the time of the TeV $γ$-ray flare. A possible superluminal knot was identified in the VLBA observations at 43 GHz. The flare constrains the size of the emitting region, and is consistent with several theoretical models with stationary shocks. △ Less

Submitted 21 August, 2017; originally announced August 2017.

Comments: 8 pages, 5 figures, in Proceedings of the 35th International Cosmic Ray Conference (ICRC 2017), Busan (South Korea)

Multiband variability studies and novel broadband SED modeling of Mrk 501 in 2009

Authors: M. L. Ahnen, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, A. Berti, B. Biasuzzi, A. Biland, O. Blanch, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, S. Buson, A. Carosi, A. Chatterjee, R. Clavero, P. Colin , et al. (268 additional authors not shown)

Abstract: We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux… ▽ More We present an extensive study of the BL Lac object Mrk 501 based on a data set collected during the multi-instrument campaign spanning from 2009 March 15 to 2009 August 1 which includes, among other instruments, MAGIC, VERITAS, Whipple 10-m, Fermi-LAT, RXTE, Swift, GASP-WEBT and VLBA. We find an increase in the fractional variability with energy, while no significant interband correlations of flux changes are found in the acquired data set. The higher variability in the very high energy (>100 GeV, VHE) gamma-ray emission and the lack of correlation with the X-ray emission indicate that the highest-energy electrons that are responsible for the VHE gamma-rays do not make a dominant contribution to the ~1 keV emission. Alternatively, there could be a very variable component contributing to the VHE gamma-ray emission in addition to that coming from the synchrotron self-Compton (SSC) scenarios. The space of SSC model parameters is probed following a dedicated grid-scan strategy, allowing for a wide range of models to be tested and offering a study of the degeneracy of model-to-data agreement in the individual model parameters. We find that there is some degeneracy in both the one-zone and the two-zone SSC scenarios that were probed, with several combinations of model parameters yielding a similar model-to-data agreement, and some parameters better constrained than others. The SSC model grid-scan shows that the flaring activity around 2009 May 22 cannot be modeled adequately with a one-zone SSC scenario, while it can be suitably described within a two-independent-zone SSC scenario. The observation of an electric vector polarization angle rotation coincident with the gamma-ray flare from 2009 May 1 resembles those reported previously for low frequency peaked blazars, hence suggesting that there are many similarities in the flaring mechanisms of blazars with different jet properties. △ Less

Submitted 30 December, 2016; originally announced December 2016.

Comments: 33 pages, 15 figures, accepted for publication in A&A. Corresponding authors: Marlene Doert ([email protected]) and David Paneque ([email protected])

Journal ref: A&A 603, A31 (2017)

Multi-epoch Spectropolarimetry of SN 2011fe

Authors: Peter A. Milne, G. Grant Williams, Amber Porter, Paul S. Smith, Nathan Smith, Mark D. Leising, Buell T. Jannuzi, E. M. Green

Abstract: We present multiple spectropolarimetric observations of the nearby Type Ia supernova (SN) 2011fe in M101, obtained before, during, and after the time of maximum apparent visual brightness. The excellent time coverage of our spectropolarimetry has allowed better monitoring of the evolution of polarization features than is typical, which has allowed us new insight into the nature of normal SNe Ia. S… ▽ More We present multiple spectropolarimetric observations of the nearby Type Ia supernova (SN) 2011fe in M101, obtained before, during, and after the time of maximum apparent visual brightness. The excellent time coverage of our spectropolarimetry has allowed better monitoring of the evolution of polarization features than is typical, which has allowed us new insight into the nature of normal SNe Ia. SN 2011fe exhibits time-dependent polarization in both the continuum and strong absorption lines. At early epochs, red wavelengths exhibit a degree of continuum polarization of up to 0.4\%, likely indicative of a mild asymmetry in the electron-scattering photosphere. This behavior is more common in sub-luminous SNe Ia than in normal events, such as SN2011fe. The degree of polarization across a collection of absorption lines varies dramatically from epoch to epoch. During the earliest epoch a $λ$4600-5000 Å complex of absorption lines shows enhanced polarization at a different position angle than the continuum. We explore the origin of these features, presenting a few possible interpretations, without arriving at a single favored ion. During two epochs near maximum, the dominant polarization feature is associated with the Si{\sc ii} $λ$6355 Å absorption line. This is common for SNeIa, but for SN2011fe the polarization of this feature increases after maximum light, whereas for other SNeIa, that polarization feature was strongest before maximum light. △ Less

Submitted 21 November, 2016; originally announced November 2016.

Comments: 15 pages, accepted by ApJ

A search for spectral hysteresis and energy-dependent time lags from X-ray and TeV gamma-ray observations of Mrk 421

Authors: A. U. Abeysekara, S. Archambault, A. Archer, W. Benbow, R. Bird, M. Buchovecky, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, X. Chen, L. Ciupik, M. P. Connolly, W. Cui, J. D. Eisch, A. Falcone, Q. Feng, J. P. Finley, H. Fleischhack, A. Flinders, L. Fortson, A. Furniss, S. Griffin, M. Hütten, N. Håkansson , et al. (198 additional authors not shown)

Abstract: Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421,… ▽ More Blazars are variable emitters across all wavelengths over a wide range of timescales, from months down to minutes. It is therefore essential to observe blazars simultaneously at different wavelengths, especially in the X-ray and gamma-ray bands, where the broadband spectral energy distributions usually peak. In this work, we report on three "target-of-opportunity" (ToO) observations of Mrk 421, one of the brightest TeV blazars, triggered by a strong flaring event at TeV energies in 2014. These observations feature long, continuous, and simultaneous exposures with XMM-Newton (covering X-ray and optical/ultraviolet bands) and VERITAS (covering TeV gamma-ray band), along with contemporaneous observations from other gamma-ray facilities (MAGIC and Fermi-LAT) and a number of radio and optical facilities. Although neither rapid flares nor significant X-ray/TeV correlation are detected, these observations reveal subtle changes in the X-ray spectrum of the source over the course of a few days. We search the simultaneous X-ray and TeV data for spectral hysteresis patterns and time delays, which could provide insight into the emission mechanisms and the source properties (e.g. the radius of the emitting region, the strength of the magnetic field, and related timescales). The observed broadband spectra are consistent with a one-zone synchrotron self-Compton model. We find that the power spectral density distribution at $\gtrsim 4\times 10^{-4}$ Hz from the X-ray data can be described by a power-law model with an index value between 1.2 and 1.8, and do not find evidence for a steepening of the power spectral index (often associated with a characteristic length scale) compared to the previously reported values at lower frequencies. △ Less

Submitted 14 November, 2016; originally announced November 2016.

Comments: 45 pages, 15 figures

Exceptional outburst of the blazar CTA 102 in 2012: The GASP-WEBT campaign and its extension

Authors: V. M. Larionov, M. Villata, C. M. Raiteri, S. G. Jorstad, A. P. Marscher, I. Agudo, P. S. Smith, J. A. Acosta-Pulido, M. J. Arévalo, A. A. Arkharov, R. Bachev, D. A. Blinov, G. Borisov, G. A. Borman, V. Bozhilov, A. Bueno, M. I. Carnerero, D. Carosati, C. Casadio, W. P. Chen, D. P. Clemens, A. Di Paola, Sh. A. Ehgamberdiev, J. L. Gómez, P. A González-Morales , et al. (27 additional authors not shown)

Abstract: After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished… ▽ More After several years of quiescence, the blazar CTA 102 underwent an exceptional outburst in 2012 September-October. The flare was tracked from gamma-ray to near-infrared frequencies, including Fermi and Swift data as well as photometric and polarimetric data from several observatories. An intensive GASP-WEBT collaboration campaign in optical and NIR bands, with an addition of previously unpublished archival data and extension through fall 2015, allows comparison of this outburst with the previous activity period of this blazar in 2004-2005. We find remarkable similarity between the optical and gamma-ray behaviour of CTA 102 during the outburst, with a time lag between the two light curves of ~1 hour, indicative of co-spatiality of the optical and gamma-ray emission regions. The relation between the gamma-ray and optical fluxes is consistent with the SSC mechanism, with a quadratic dependence of the SSC gamma-ray flux on the synchrotron optical flux evident in the post-outburst stage. However, the gamma-ray/optical relationship is linear during the outburst; we attribute this to changes in the Doppler factor. A strong harder-when-brighter spectral dependence is seen both the in gamma-ray and optical non-thermal emission. This hardening can be explained by convexity of the UV-NIR spectrum that moves to higher frequencies owing to an increased Doppler shift as the viewing angle decreases during the outburst stage. The overall pattern of Stokes parameter variations agrees with a model of a radiating blob or shock wave that moves along a helical path down the jet. △ Less

Submitted 21 June, 2016; originally announced June 2016.

Comments: 11 pages, 14 figures, accepted in MNRAS

Polarization angle swings in blazars: The case of 3C 279

Authors: S. Kiehlmann, T. Savolainen, S. G. Jorstad, K. V. Sokolovsky, F. K. Schinzel, A. P. Marscher, V. M. Larionov, I. Agudo, H. Akitaya, E. Benítez, A. Berdyugin, D. A. Blinov, N. G. Bochkarev, G. A. Borman, A. N. Burenkov, C. Casadio, V. T. Doroshenko, N. V. Efimova, Y. Fukazawa, J. L. Gómez, T. S. Grishina, V. A. Hagen-Thorn, J. Heidt, D. Hiriart, R. Itoh , et al. (39 additional authors not shown)

Abstract: Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting regio… ▽ More Over the past few years, several occasions of large, continuous rotations of the electric vector position angle (EVPA) of linearly polarized optical emission from blazars have been reported. These events are often coincident with high energy gamma-ray flares and they have attracted considerable attention, as they could allow one to probe the magnetic field structure in the gamma-ray emitting region of the jet. The flat-spectrum radio quasar 3C279 is one of the most prominent examples showing this behaviour. Our goal is to study the observed EVPA rotations and to distinguish between a stochastic and a deterministic origin of the polarization variability. We have combined multiple data sets of R-band photometry and optical polarimetry measurements of 3C279, yielding exceptionally well-sampled flux density and polarization curves that cover a period of 2008-2012. Several large EVPA rotations are identified in the data. We introduce a quantitative measure for the EVPA curve smoothness, which is then used to test a set of simple random walk polarization variability models against the data. 3C279 shows different polarization variation characteristics during an optical low-flux state and a flaring state. The polarization variation during the flaring state, especially the smooth approx. 360 degrees rotation of the EVPA in mid-2011, is not consistent with the tested stochastic processes. We conclude that during the two different optical flux states, two different processes govern the polarization variation, possibly a stochastic process during the low-brightness state and a deterministic process during the flaring activity. △ Less

Submitted 21 June, 2016; v1 submitted 1 March, 2016; originally announced March 2016.

Comments: 19 pages, 15 figures, 6 tables

Journal ref: A&A 592, C1 (2016)

Gamma rays from the quasar PKS 1441+25: story of an escape

Authors: A. U. Abeysekara, S. Archambault, A. Archer, T. Aune, A. Barnacka, W. Benbow, R. Bird, J. Biteau, J. H. Buckley, V. Bugaev, J. V Cardenzana, M. Cerruti, X. Chen, J. L. Christiansen, L. Ciupik, M. P. Connolly, P. Coppi, W. Cui, H. J. Dickinson, J. Dumm, J. D. Eisch, M. Errando, A. Falcone, Q. Feng, J. P. Finley , et al. (81 additional authors not shown)

Abstract: Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EB… ▽ More Outbursts from gamma-ray quasars provide insights on the relativistic jets of active galactic nuclei and constraints on the diffuse radiation fields that fill the Universe. The detection of significant emission above 100 GeV from a distant quasar would show that some of the radiated gamma rays escape pair-production interactions with low-energy photons, be it the extragalactic background light (EBL), or the radiation near the supermassive black hole lying at the jet's base. VERITAS detected gamma-ray emission up to 200 GeV from PKS 1441+25 (z=0.939) during April 2015, a period of high activity across all wavelengths. This observation of PKS 1441+25 suggests that the emission region is located thousands of Schwarzschild radii away from the black hole. The gamma-ray detection also sets a stringent upper limit on the near-ultraviolet to near-infrared EBL intensity, suggesting that galaxy surveys have resolved most, if not all, of the sources of the EBL at these wavelengths. △ Less

Submitted 14 December, 2015; originally announced December 2015.

Comments: 7 pages, 3 figures, published in ApJ Letters 815, L22 (2015)

A multi-wavelength polarimetric study of the blazar CTA 102 during a Gamma-ray flare in 2012

Authors: Carolina Casadio, José L. Gómez, Svetlana G. Jorstad, Alan P. Marscher, Valeri M. Larionov, Paul S. Smith, Mark A. Gurwell, Anne Lähteenmäki, Iván Agudo, Sol N. Molina, Vishal Bala, Manasvita Joshi, Brian Taylor, Karen E. Williamson, Arkady A. Arkharov, Dmitry A. Blinov, George A. Borman, Andrea Di Paola, Tatiana S. Grishina, Vladimir A. Hagen-Thorn, Ryosuke Itoh, Evgenia N. Kopatskaya, Elena G. Larionova, Liudmila V. Larionova, Daria A. Morozova , et al. (6 additional authors not shown)

Abstract: We perform a multi-wavelength polarimetric study of the quasar CTA 102 during an extraordinarily bright $γ$-ray outburst detected by the {\it Fermi} Large Area Telescope in September-October 2012 when the source reached a flux of F$_{>100~\mathrm{MeV}} =5.2\pm0.4\times10^{-6}$ photons cm$^{-2}$ s$^{-1}$. At the same time the source displayed an unprecedented optical and NIR outburst. We study the… ▽ More We perform a multi-wavelength polarimetric study of the quasar CTA 102 during an extraordinarily bright $γ$-ray outburst detected by the {\it Fermi} Large Area Telescope in September-October 2012 when the source reached a flux of F$_{>100~\mathrm{MeV}} =5.2\pm0.4\times10^{-6}$ photons cm$^{-2}$ s$^{-1}$. At the same time the source displayed an unprecedented optical and NIR outburst. We study the evolution of the parsec scale jet with ultra-high angular resolution through a sequence of 80 total and polarized intensity Very Long Baseline Array images at 43 GHz, covering the observing period from June 2007 to June 2014. We find that the $γ$-ray outburst is coincident with flares at all the other frequencies and is related to the passage of a new superluminal knot through the radio core. The powerful $γ$-ray emission is associated with a change in direction of the jet, which became oriented more closely to our line of sight ($θ\sim$1.2$^{\circ}$) during the ejection of the knot and the $γ$-ray outburst. During the flare, the optical polarized emission displays intra-day variability and a clear clockwise rotation of EVPAs, which we associate with the path followed by the knot as it moves along helical magnetic field lines, although a random walk of the EVPA caused by a turbulent magnetic field cannot be ruled out. We locate the $γ$-ray outburst a short distance downstream of the radio core, parsecs from the black hole. This suggests that synchrotron self-Compton scattering of near-infrared to ultraviolet photons is the probable mechanism for the $γ$-ray production. △ Less

Submitted 10 September, 2015; v1 submitted 28 August, 2015; originally announced August 2015.

Comments: Accepted for publication in The Astrophysical Journal

Spectropolarimetry of SN 2011dh in M51: geometric insights on a Type IIb supernova progenitor and explosion

Authors: Jon C. Mauerhan, G. Grant Williams, Douglas C. Leonard, Paul S. Smith, Alexei V. Filippenko, Nathan Smith, Jennifer L. Hoffman, Leah Huk, Kelsey I. Clubb, Jeffrey M. Silverman, S. Bradley Cenko, Peter Milne, Avishay Gal-Yam, Sagi Ben-Ami

Abstract: We present seven epochs of spectropolarimetry of the Type IIb supernova (SN) 2011dh in M51, spanning 86 days of its evolution. The first epoch was obtained 9 days after the explosion, when the photosphere was still in the depleted hydrogen layer of the stripped-envelope progenitor. Continuum polarization is securely detected at the level of P~0.5% through day 14 and appears to diminish by day 30,… ▽ More We present seven epochs of spectropolarimetry of the Type IIb supernova (SN) 2011dh in M51, spanning 86 days of its evolution. The first epoch was obtained 9 days after the explosion, when the photosphere was still in the depleted hydrogen layer of the stripped-envelope progenitor. Continuum polarization is securely detected at the level of P~0.5% through day 14 and appears to diminish by day 30, which is different from the prevailing trends suggested by studies of other core-collapse SNe. Time-variable modulations in P and position angle are detected across P-Cygni line features. H-alpha and HeI polarization peak after 30 days and exhibit position angles roughly aligned with the earlier continuum, while OI and CaII appear to be geometrically distinct. We discuss several possibilities to explain the evolution of the continuum and line polarization, including the potential effects of a tidally deformed progenitor star, aspherical radioactive heating by fast-rising plumes of Ni-56 from the core, oblique shock breakout, or scattering by circumstellar material. While these possibilities are plausible and guided by theoretical expectations, they are not unique solutions to the data. The construction of more detailed hydrodynamic and radiative-transfer models that incorporate complex aspherical geometries will be required to further elucidate the nature of the polarized radiation from SN 2011dh and other Type IIb supernovae. △ Less

Submitted 4 September, 2015; v1 submitted 29 June, 2015; originally announced June 2015.

Comments: Post-proof edit. Accepted to MNRAS 2015 Aug 19

Multiwavelength behaviour of the blazar OJ 248 from radio to γ-rays

Authors: M. I. Carnerero, C. M. Raiteri, M. Villata, J. A. Acosta-Pulido, F. D'Ammando, P. S. Smith, V. M. Larionov

Abstract: We present an analysis of the multiwavelength behaviour of the blazar OJ 248 at z = 0.939 in the period 2006-2013. We use low-energy data (optical, near-infrared, and radio) obtained by 21 observatories participating in the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), as well as data from the Swift (optical-UV and X-rays) and Fermi (gamma-rays) satellites, to stud… ▽ More We present an analysis of the multiwavelength behaviour of the blazar OJ 248 at z = 0.939 in the period 2006-2013. We use low-energy data (optical, near-infrared, and radio) obtained by 21 observatories participating in the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT), as well as data from the Swift (optical-UV and X-rays) and Fermi (gamma-rays) satellites, to study flux and spectral variability and correlations among emissions in different bands. We take into account the effect of absorption by the Damped Lyman Alpha intervening system at z = 0.525. Two major outbursts were observed in 2006-2007 and in 2012-2013 at optical and near-IR wavelengths, while in the high-frequency radio light curves prominent radio outbursts are visible peaking at the end of 2010 and beginning of 2013, revealing a complex radio-optical correlation. Cross-correlation analysis suggests a delay of the optical variations after the gamma-ray ones of about a month, which is a peculiar behaviour in blazars. We also analyse optical polarimetric and spectroscopic data. The average polarization percentage P is less than 3 per cent, but it reaches about 19 per cent during the early stage of the 2012-2013 outburst. A vague correlation of P with brightness is observed. There is no preferred electric vector polarisation angle and during the outburst the linear polarization vector shows wide rotations in both directions, suggesting a complex behaviour or structure of the jet and possible turbulence. The analysis of 140 optical spectra acquired at the Steward Observatory reveals a strong Mg II broad emission line with an essentially stable flux of 6.2 e-15 erg cm-2 s-1 and a full width at half-maximum of 2053 km s-1. △ Less

Submitted 5 May, 2015; originally announced May 2015.

Comments: 16 pages, 20 figures, 2 tables, in press for MNRAS

The SMARTS Multi-epoch Optical Spectroscopy Atlas (SAMOSA): Using Emission Line Variability to Probe the Location of the Blazar Gamma-emitting Region

Authors: Jedidah C. Isler, C. M. Urry, C. Bailyn, P. S. Smith, P. Coppi, M. Brady, E. Macpherson, I. Hasan, M. Buxton

Abstract: We present multi-epoch optical spectroscopy of seven southern Fermi-monitored blazars from 2008 - 2013 using the Small and Medium Aperture Research Telescope System (SMARTS), with supplemental spectroscopy and polarization data from the Steward Observatory. We find that the emission lines are much less variable than the continuum; 4 of 7 blazars had no detectable emission line variability over the… ▽ More We present multi-epoch optical spectroscopy of seven southern Fermi-monitored blazars from 2008 - 2013 using the Small and Medium Aperture Research Telescope System (SMARTS), with supplemental spectroscopy and polarization data from the Steward Observatory. We find that the emission lines are much less variable than the continuum; 4 of 7 blazars had no detectable emission line variability over the 5 years. This is consistent with photoionization primarily by an accretion disk, allowing us to use the lines as a probe of disk activity. Comparing optical emission line flux with Fermi $γ$-ray flux and optical polarized flux, we investigate whether relativistic jet variability is related to the accretion flow. In general, we see no such dependence, suggesting the jet variability is likely caused by internal processes like turbulence or shock acceleration rather than a variable accretion rate. However, three sources showed statistically significant emission line flares in close temporal proximity to very large Fermi $γ$-ray flares. While we do not have sufficient emission line data to quantitatively assess their correlation with the $γ$-ray flux, it appears that in some cases, the jet might provide additional photoionizing flux to the broad line region, which implies some gamma-rays are produced within the broad line region, at least for these large flares. △ Less

Submitted 3 February, 2015; originally announced February 2015.

Comments: 15 pages, 15 figures, Accepted to the Astrophysical Journal

Searching for Planets in Holey Debris Disks with the Apodizing Phase Plate

Authors: Tiffany Meshkat, Vanessa P. Bailey, Kate Y. L. Su, Matthew A. Kenworthy, Eric E. Mamajek, Philip M. Hinz, Paul S. Smith

Abstract: We present our first results from a high-contrast imaging search for planetary mass companions around stars with gapped debris disks, as inferred from the stars' bright infrared excesses. For the six considered stars, we model the disks' unresolved infrared spectral energy distributions (SEDs) in order to derive the temperature and location of the disk components. With VLT/NaCo Apodizing Phase Pla… ▽ More We present our first results from a high-contrast imaging search for planetary mass companions around stars with gapped debris disks, as inferred from the stars' bright infrared excesses. For the six considered stars, we model the disks' unresolved infrared spectral energy distributions (SEDs) in order to derive the temperature and location of the disk components. With VLT/NaCo Apodizing Phase Plate coronagraphic L'-band imaging, we search for planetary mass companions that may be sculpting the disks. We detect neither disks nor companions in this sample, confirmed by comparing plausible point sources with archival data. In order to calculate our mass sensitivity limit, we revisit the stellar age estimates. One target, HD 17848, at 540$\pm$100 Myr old is significantly older than previously estimated. We then discuss our high-contrast imaging results with respect to the disk properties. △ Less

Submitted 16 December, 2014; originally announced December 2014.

Comments: Accepted for publication in ApJ, 10 pages, 3 figures

Unprecedented study of the broadband emission of Mrk 421 during flaring activity in March 2010

Authors: The MAGIC Collaboration, J. Aleksić, S. Ansoldi, L. A. Antonelli, P. Antoranz, A. Babic, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, B. Biasuzzi, A. Biland, O. Blanch, A. Boller, S. Bonnefoy, G. Bonnoli, F. Borracci, T. Bretz, E. Carmona, A. Carosi, P. Colin, E. Colombo, J. L. Contreras , et al. (230 additional authors not shown)

Abstract: A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameter… ▽ More A flare from the TeV blazar Mrk 421, occurring in March 2010, was observed for 13 consecutive days from radio to very high energy (VHE, E > 100 GeV) gamma-rays with MAGIC, VERITAS, Whipple, FermiLAT, MAXI, RXTE, Swift, GASP-WEBT, and several optical and radio telescopes. We model the day-scale SEDs with one-zone and two-zone synchrotron self-Compton (SSC) models, investigate the physical parameters, and evaluate whether the observed broadband SED variability can be associated to variations in the relativistic particle population. Flux variability was remarkable in the X-ray and VHE bands while it was minor or not significant in the other bands. The one-zone SSC model can describe reasonably well the SED of each day for the 13 consecutive days. This flaring activity is also very well described by a two-zone SSC model, where one zone is responsible for the quiescent emission while the other smaller zone, which is spatially separated from the first one, contributes to the daily-variable emission occurring in X-rays and VHE gamma-rays. Both the one-zone SSC and the two-zone SSC models can describe the daily SEDs via the variation of only four or five model parameters, under the hypothesis that the variability is associated mostly to the underlying particle population. This shows that the particle acceleration and cooling mechanism producing the radiating particles could be the main one responsible for the broadband SED variations during the flaring episodes in blazars. The two-zone SSC model provides a better agreement to the observed SED at the narrow peaks of the low- and high-energy bumps during the highest activity, although the reported one-zone SSC model could be further improved by the variation of the parameters related to the emitting region itself ($δ$, $B$ and $R$), in addition to the parameters related to the particle population. △ Less

Submitted 14 June, 2015; v1 submitted 11 December, 2014; originally announced December 2014.

Comments: Published in A&A, 30 pages, 15 figures, 6 tables. Online-data: multi-wavelength light curves (data in Fig. 1) and broadband spectral energy distributions (the data in Figs. 6, and B1-B4) are available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/578/A22. Corresponding authors: David Paneque ([email protected]), Shangyu Sun ([email protected]), Hajime Takami ([email protected])

Journal ref: A&A 578, A22 (2015)

Debris Distribution in HD 95086 - A Young Analog of HR 8799

Authors: Kate Y. L. Su, Sarah Morrison, Renu Malhotra, Paul S. Smith, Zoltan Balog, George H. Rieke

Abstract: HD 95086 is a young early-type star that hosts (1) a 5 MJ planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 um crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 um. Due to the low resolution of MIPS-SED mode, this feature is not spectrally resolved, but is consistent with… ▽ More HD 95086 is a young early-type star that hosts (1) a 5 MJ planet at the projected distance of 56 AU revealed by direct imaging, and (2) a prominent debris disk. Here we report the detection of 69 um crystalline olivine feature from the disk using the Spitzer/MIPS-SED data covering 55-95 um. Due to the low resolution of MIPS-SED mode, this feature is not spectrally resolved, but is consistent with the emission from crystalline forsterite contributing 5% of the total dust mass. We also present detailed analysis of the disk SED and re-analysis of resolved images obtained by Herschel. Our results suggest that the debris structure around HD 95086 consists of a warm (175 K) belt, a cold (55 K) disk, and an extended disk halo (up to 800 AU), and is very similar to that of HR 8799. We compare the properties of the three debris components, and suggest that HD 95086 is a young analog of HR 8799. We further investigate and constrain single-planet, two-planet, three-planet and four-planet architectures that can account for the observed debris structure and are compatible with dynamical stability constraints. We find that equal-mass four-planet configurations of geometrically spaced orbits, with each planet of mass 5 MJ, could maintain the gap between the warm and cold debris belts, and also be just marginally stable for timescales comparable to the age of the system. △ Less

Submitted 29 November, 2014; originally announced December 2014.

Comments: accepted for publication in ApJ

On the Morphology and Chemical Composition of the HR 4796A Debris Disk

Authors: Timothy J. Rodigas, Christopher C. Stark, Alycia Weinberger, John H. Debes, Philip M. Hinz, Laird Close, Christine Chen, Paul S. Smith, Jared R. Males, Andrew J. Skemer, Alfio Puglisi, Katherine B. Follette, Katie Morzinski, Ya-Lin Wu, Runa Briguglio, Simone Esposito, Enrico Pinna, Armando Riccardi, Glenn Schneider, Marco Xompero

Abstract: [abridged] We present resolved images of the HR 4796A debris disk using the Magellan adaptive optics system paired with Clio-2 and VisAO. We detect the disk at 0.77 \microns, 0.91 \microns, 0.99 \microns, 2.15 \microns, 3.1 \microns, 3.3 \microns, and 3.8 \microns. We find that the deprojected center of the ring is offset from the star by 4.76$\pm$1.6 AU and that the deprojected eccentricity is 0.… ▽ More [abridged] We present resolved images of the HR 4796A debris disk using the Magellan adaptive optics system paired with Clio-2 and VisAO. We detect the disk at 0.77 \microns, 0.91 \microns, 0.99 \microns, 2.15 \microns, 3.1 \microns, 3.3 \microns, and 3.8 \microns. We find that the deprojected center of the ring is offset from the star by 4.76$\pm$1.6 AU and that the deprojected eccentricity is 0.06$\pm$0.02, in general agreement with previous studies. We find that the average width of the ring is 14$^{+3}_{-2}%$, also comparable to previous measurements. Such a narrow ring precludes the existence of shepherding planets more massive than \about 4 \mj, comparable to hot-start planets we could have detected beyond \about 60 AU in projected separation. Combining our new scattered light data with archival HST/STIS and HST/NICMOS data at \about 0.5-2 \microns, along with previously unpublished Spitzer/MIPS thermal emission data and all other literature thermal data, we set out to constrain the chemical composition of the dust grains. After testing 19 individual root compositions and more than 8,400 unique mixtures of these compositions, we find that good fits to the scattered light alone and thermal emission alone are discrepant, suggesting that caution should be exercised if fitting to only one or the other. When we fit to both the scattered light and thermal emission simultaneously, we find mediocre fits (reduced chi-square \about 2). In general, however, we find that silicates and organics are the most favored, and that water ice is usually not favored. These results suggest that the common constituents of both interstellar dust and solar system comets also may reside around HR 4796A, though improved modeling is necessary to place better constraints on the exact chemical composition of the dust. △ Less

Submitted 28 October, 2014; originally announced October 2014.

Comments: Accepted to ApJ on October 27, 2014. 21 pages, 12 figures, 4 tables

Comprehensive Monitoring of Gamma-ray Bright Blazars. I. Statistical Study of Optical, X-ray, and Gamma-ray Spectral Slopes

Authors: Karen E. Williamson, Svetlana G. Jorstad, Alan P. Marscher, Valeri M. Larionov, Paul S. Smith, Iván Agudo, Arkady A. Arkharov, Dmitry A. Blinov, Carolina Casadio, Natalia V. Efimova, José L. Gómez, Vladimir A. Hagen-Thorn, Manasvita Joshi, Tatiana S. Konstantinova, Evgenia N. Kopatskaya, Elena G. Larionova, Liudmilla V. Larionova, Michael P. Malmrose, Ian M. McHardy, Sol N. Molina, Daria A. Morozova, Gary D. Schmidt, Brian W. Taylor, Ivan S. Troitsky

Abstract: We present $γ$-ray, X-ray, ultraviolet, optical, and near-infrared light curves of 33 $γ$-ray bright blazars over four years that we have been monitoring since 2008 August with multiple optical, ground-based telescopes and the Swift satellite, and augmented by data from the Fermi Gamma-ray Space Telescope and other publicly available data from Swift. The sample consists of 21 flat-spectrum radio q… ▽ More We present $γ$-ray, X-ray, ultraviolet, optical, and near-infrared light curves of 33 $γ$-ray bright blazars over four years that we have been monitoring since 2008 August with multiple optical, ground-based telescopes and the Swift satellite, and augmented by data from the Fermi Gamma-ray Space Telescope and other publicly available data from Swift. The sample consists of 21 flat-spectrum radio quasars (FSRQs) and 12 BL Lac objects (BL Lacs). We identify quiescent and active states of the sources based on their $γ$-ray behavior. We derive $γ$-ray, X-ray, and optical spectral indices, $α_γ$, $α_X$, and $α_o$, respectively ($F_ν\proptoν^α$), and construct spectral energy distributions (SEDs) during quiescent and active states. We analyze the relationships between different spectral indices, blazar classes, and activity states. We find (i) significantly steeper $γ$-ray spectra of FSRQs than for BL Lacs during quiescent states, but a flattening of the spectra for FSRQs during active states while the BL Lacs show no significant change; (ii) a small difference of $α_X$ within each class between states, with BL Lac X-ray spectra significantly steeper than in FSRQs; (iii) a highly peaked distribution of X-ray spectral slopes of FSRQs at $\sim-$0.60, but a very broad distribution of $α_X$ of BL Lacs during active states; (iv) flattening of the optical spectra of FSRQs during quiescent states, but no statistically significant change of $α_o$ of BL Lacs between states; and (v) a positive correlation between optical and $γ$-ray spectral slopes of BL Lacs, with similar values of the slopes. We discuss the findings with respect to the relative prominence of different components of high-energy and optical emission as the flux state changes. △ Less

Submitted 10 June, 2014; originally announced June 2014.

Comments: 40 pages, 28 figures, 10 Tables. Additional figures at http://www.bu.edu/blazars/VLBAproject.html. Accepted by ApJ

Multi-epoch spectropolarimetry of SN 2009ip: direct evidence for aspherical circumstellar material

Authors: Jon C. Mauerhan, G. Grant Williams, Nathan Smith, Paul S. Smith, Alexei V. Filippenko, Jennifer L. Hoffman, Douglas C. Leonard, Peter Milne, Kelsey I. Clubb, Ori D. Fox, Patrick L. Kelly

Abstract: We present spectropolarimetry of SN 2009ip throughout the evolution of its 2012 explosion. During the initial 2012a phase, when the source spectrum exhibits broad P-Cygni lines, we measure a V-band polarization of P~0.9% at a position angle of ~166 deg, indicating substantial asphericity for the 2012a outflow. Near the subsequent peak of the 2012b phase, when the spectrum shows signs of intense in… ▽ More We present spectropolarimetry of SN 2009ip throughout the evolution of its 2012 explosion. During the initial 2012a phase, when the source spectrum exhibits broad P-Cygni lines, we measure a V-band polarization of P~0.9% at a position angle of ~166 deg, indicating substantial asphericity for the 2012a outflow. Near the subsequent peak of the 2012b phase, when the spectrum shows signs of intense interaction with circumstellar material (CSM), we measure P~1.7% at a position angle of 72 deg, indicating a separate physical component of polarization, which has a higher degree of asphericity than the 2012a outflow and an orthogonal axis of symmetry on the sky. Around 30 days past peak, which is coincident with a bump in the declining light curve, we measure P~0.7% and significant variations in P across some broad lines, particularly HeI/NaI. By 60 days past peak the source appears to be approaching a low value of interstellar polarization (P<0.2%). The results are consistent with a scenario in which a potentially bipolar explosion during 2012a impacts a toroidal distribution of CSM, thereby causing the 2012b brightening. Orthogonality between the 2012a and 2012b geometries is inconsistent with the hypothesis that the 2012a event launched the CSM that was hit 40 days later by material from the 2012b event. Rather, the aspherical CSM probably has a separate origin associated with the outbursts of the stellar progenitor during prior years. Previous calculations that assumed spherical symmetry for the CSM have underestimated the required explosion energy, as the results suggest that <10% of the fast SN ejecta participated in strong CSM interaction during the 2012b phase. In light of the spectropolarimetric results, a kinetic energy of 1e51 erg for the ejecta is difficult to avoid, supporting the interpretation that the 2012 outburst was the result of a core-collapse supernova explosion (abridged) △ Less

Submitted 15 December, 2014; v1 submitted 17 March, 2014; originally announced March 2014.

Comments: Accepted to MNRAS 2014 Apr 7. Accepted version. 15 pages, 10 figures

Journal ref: 2014, MNRAS, 442, 1166

The awakening of BL Lacertae: observations by Fermi, Swift, and the GASP-WEBT

Authors: C. M. Raiteri, M. Villata, F. D'Ammando, V. M. Larionov, M. A. Gurwell, D. O. Mirzaqulov, P. S. Smith

Abstract: Since the launch of the Fermi satellite, BL Lacertae has been moderately active at gamma-rays and optical frequencies until May 2011, when the source started a series of strong flares. The exceptional optical sampling achieved by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) in collaboration with the Steward Observatory allows us to perform a detailed comparison… ▽ More Since the launch of the Fermi satellite, BL Lacertae has been moderately active at gamma-rays and optical frequencies until May 2011, when the source started a series of strong flares. The exceptional optical sampling achieved by the GLAST-AGILE Support Program (GASP) of the Whole Earth Blazar Telescope (WEBT) in collaboration with the Steward Observatory allows us to perform a detailed comparison with the daily gamma-ray observations by Fermi. Discrete correlation analysis between the optical and gamma-ray emission reveals correlation with a time lag of 0 +- 1 d, which suggests cospatiality of the corresponding jet emitting regions. A better definition of the time lag is hindered by the daily gaps in the sampling of the extremely fast flux variations. In general, optical flares present more structure and develop on longer time scales than corresponding gamma-ray flares. Observations at X-rays and at millimetre wavelengths reveal a common trend, which suggests that the region producing the mm and X-ray radiation is located downstream from the optical and gamma-ray-emitting zone in the jet. The mean optical degree of polarisation slightly decreases over the considered period and in general it is higher when the flux is lower. The optical electric vector polarisation angle (EVPA) shows a preferred orientation of about 15 deg, nearly aligned with the radio core EVPA and mean jet direction. Oscillations around it increase during the 2011-2012 outburst. We investigate the effects of a geometrical interpretation of the long-term flux variability on the polarisation. A helical magnetic field model predicts an evolution of the mean polarisation that is in reasonable agreement with the observations. These can be fully explained by introducing slight variations in the compression factor in a transverse shock waves model. △ Less

Submitted 5 September, 2013; originally announced September 2013.

Comments: 17 pages, 17 figures, 1 table, accepted for publication in MNRAS

Journal ref: MNRAS 436, 1530 (2013)

Discovery of a new TeV gamma-ray source: VER J0521+211

Authors: VERITAS Collaboration, S. Archambault, T. Arlen, T. Aune, M. Beilicke, W. Benbow, R. Bird, A. Bouvier, J. H. Buckley, V. Bugaev, K. Byrum, A. Cesarini, L. Ciupik, M. P. Connolly, W. Cui, M. Errando, A. Falcone, S. Federici, Q. Feng, J. P. Finley, L. Fortson, A. Furniss, N. Galante, D. Gall, G. H. Gillanders , et al. (61 additional authors not shown)

Abstract: We report the detection of a new TeV gamma-ray source, VER J0521+211, based on observations made with the VERITAS imaging atmospheric Cherenkov telescope array. These observations were motivated by the discovery of a cluster of >30GeV photons in the first year of Fermi-LAT observations. VER J0521+211 is relatively bright at TeV energies, with a mean photon flux of 1.93 +/- 0.13_stat +/- 0.78_sys 1… ▽ More We report the detection of a new TeV gamma-ray source, VER J0521+211, based on observations made with the VERITAS imaging atmospheric Cherenkov telescope array. These observations were motivated by the discovery of a cluster of >30GeV photons in the first year of Fermi-LAT observations. VER J0521+211 is relatively bright at TeV energies, with a mean photon flux of 1.93 +/- 0.13_stat +/- 0.78_sys 10^-11 cm-2 s-1 above 0.2 TeV during the period of the VERITAS observations. The source is strongly variable on a daily timescale across all wavebands, from optical to TeV, with a peak flux corresponding to ~0.3 times the steady Crab Nebula flux at TeV energies. Follow-up observations in the optical and X-ray bands classify the newly-discovered TeV source as a BL Lac-type blazar with uncertain redshift, although recent measurements suggest z=0.108. VER J0521+211 exhibits all the defining properties of blazars in radio, optical, X-ray, and gamma-ray wavelengths. △ Less

Submitted 22 August, 2013; originally announced August 2013.

Comments: 10 pages, 8 figures, 3 tables. Accepted by The Astrophysical Journal

Journal ref: The Astrophysical Journal, Volume 776, Issue 2, article id. 69, 10 pp. (2013)

Candidate Type II Quasars at 2 < z < 4.3 in the Sloan Digital Sky Survey III

Authors: Rachael Alexandroff, Michael A. Strauss, Jenny E. Greene, Nadia L. Zakamska, Nicholas P. Ross, W. N. Brandt, Guilin Liu, Paul S. Smith, Jian Ge, Fred Hamann, Adam D. Myers, Patrick Petitjean, Donald P. Schneider, Hassen Yesuf, Donald G. York

Abstract: At low redshifts, dust-obscured quasars often have strong yet narrow permitted lines in the rest-frame optical and ultraviolet, excited by the central active nucleus, earning the designation Type II quasars. We present a sample of 145 candidate Type II quasars at redshifts between 2 and 4.3, encompassing the epoch at which quasar activity peaked in the universe. These objects, selected from the qu… ▽ More At low redshifts, dust-obscured quasars often have strong yet narrow permitted lines in the rest-frame optical and ultraviolet, excited by the central active nucleus, earning the designation Type II quasars. We present a sample of 145 candidate Type II quasars at redshifts between 2 and 4.3, encompassing the epoch at which quasar activity peaked in the universe. These objects, selected from the quasar sample of the Baryon Oscillation Spectroscopic Survey of the Sloan Digital Sky Survey III, are characterized by weak continuum in the rest-frame ultraviolet (typical continuum magnitude of i \approx 22) and strong lines of CIV and Ly α, with Full Width at Half Maximum less than 2000 kms-1. The continuum magnitudes correspond to an absolute magnitude of -23 or brighter at redshift 3, too bright to be due exclusively to the host galaxies of these objects. Roughly one third of the objects are detected in the shorter-wavelength bands of the WISE survey; the spectral energy distributions (SEDs) of these objects appear to be intermediate between classic Type I and Type II quasars seen at lower redshift. Five objects are detected at rest frame 6μm by Spitzer, implying bolometric luminosities of several times 10^46 erg s-1. We have obtained polarization measurements for two objects; they are roughly 3% polarized. We suggest that these objects are luminous quasars, with modest dust extinction (A_V ~ 0.5 mag), whose ultraviolet continuum also includes a substantial scattering contribution. Alternatively, the line of sight to the central engines of these objects may be partially obscured by optically thick material. △ Less

Submitted 9 August, 2013; v1 submitted 27 July, 2013; originally announced July 2013.

Comments: 26 pages, 13 figures, 10 tables, 4 machine readable tables. Accepted for publication in MNRAS

A Tight Connection between Gamma-Ray Outbursts and Parsec-Scale Jet Activity in the Quasar 3C 454.3

Authors: Svetlana G. Jorstad, Alan P. Marscher, Paul S. Smith, Valeri M. Larionov, Iván Agudo, Mark Gurwell, Ann E. Wehrle, Anne Lähteenmäki, Maria G. Nikolashvili, Gary D. Schmidt, Arkady A. Arkharov, Dmitry A. Blinov, Kelly Blumenthal, Carolina Casadio, Revaz A. Chigladze, Natatia V. Efimova, Joseph R. Eggen, José L. Gómez, Dirk Grupe, Vladimir A. Hagen-Thorn, Manasvita Joshi, Givi N. Kimeridze, Tatiana S. Konstantinova, Evgenia N. Kopatskaya, Omar M. Kurtanidze , et al. (20 additional authors not shown)

Abstract: We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent γ-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each γ-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events… ▽ More We analyze the multifrequency behavior of the quasar 3C 454.3 during three prominent γ-ray outbursts: 2009 Autumn, 2010 Spring, and 2010 Autumn. The data reveal a repeating pattern, including a triple flare structure, in the properties of each γ-ray outburst, which implies similar mechanism(s) and location for all three events. The multi-frequency behavior indicates that the lower frequency events are co-spatial with the γ-ray outbursts, although the γ-ray emission varies on the shortest timescales. We determine that the variability from UV to IR wavelengths during an outburst results from a single synchrotron component whose properties do not change significantly over the different outbursts. Despite a general increase in the degree of optical linear polarization during an outburst, the polarization drops significantly at the peak of the γ-ray event, which suggests that both shocks and turbulent processes are involved. We detect two disturbances (knots) with superluminal apparent speeds in the parsec-scale jet associated with the outbursts in 2009 Autumn and 2010 Autumn. The kinematic properties of the knots can explain the difference in amplitudes of the γ-ray events, while their millimeter-wave polarization is related to the optical polarization during the outbursts. We interpret the multi-frequency behavior within models involving either a system of standing conical shocks or magnetic reconnection events located in the parsec-scale millimeter-wave core of the jet. We argue that γ-ray outbursts with variability timescales as short as ~ 3 hr can occur on parsec scales if flares take place in localized regions such as turbulent cells. △ Less

Submitted 9 July, 2013; originally announced July 2013.

Comments: 80 pages, 10 Tables, 31 figures, accepted for publication in ApJ

Locating the Gamma-ray Flaring Emission of Blazar AO 0235+164 in the Jet at Parsec Scales Through Multi Spectral Range Monitoring

Authors: Ivan Agudo, Alan P. Marscher, Svetlana G. Jorstad, Valeri M. Larionov, Jose L. Gomez, Anne Lahteenmaki, Paul S. Smith, Kari Nilsson, Anthony C. S. Readhead, Margo F. Aller, Jochen Heidt, Mark Gurwell, Clemens Thum, Ann E. Wehrle, Omar M. Kurtanidze

Abstract: We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164 in 2008. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with ~0.15 milliarcsecond resolution. The association of the events at different wav… ▽ More We present observations of a major outburst at centimeter, millimeter, optical, X-ray, and gamma-ray wavelengths of the BL Lacertae object AO 0235+164 in 2008. We analyze the timing of multi-waveband variations in the flux and linear polarization, as well as changes in Very Long Baseline Array (VLBA) images at 7mm with ~0.15 milliarcsecond resolution. The association of the events at different wavebands is confirmed at high statistical significance by probability arguments and Monte-Carlo simulations. A series of sharp peaks in optical linear polarization, as well as a pronounced maximum in the 7mm polarization of a superluminal jet knot, indicate rapid fluctuations in the degree of ordering of the magnetic field. These results lead us to conclude that the outburst occurred in the jet both in the quasi-stationary core and in the superluminal knot, both at >12 parsecs downstream of the supermassive black hole. We interpret the outburst as a consequence of the propagation of a disturbance, elongated along the line of sight by light-travel time delays, that passes through a standing recollimation shock in the core and propagates down the jet to create the superluminal knot. The multi-wavelength light curves vary together on long time-scales (months/years), but the correspondence is poorer on shorter time-scales. This, as well as the variability of the polarization and the dual location of the outburst, agrees with the expectations of a multi-zone emission model in which turbulence plays a major role in modulating the synchrotron and inverse Compton fluxes. △ Less

Submitted 8 March, 2013; originally announced March 2013.

Comments: 2012 Fermi Symposium proceedings - eConf C121028

Time-resolved Spectroscopy of the Polar EU Cancri in the Open Cluster Messier 67

Authors: Kurtis A. Williams, Steve B. Howell, James Liebert, Paul S. Smith, Andrea Bellini, Kate H. R. Rubin, Michael Bolte

Abstract: We present time-resolved spectroscopic and polarimetric observations of the AM Her system EU Cnc. EU Cnc is located near the core of the old open cluster Messier 67; new proper motion measurements indicate that EU Cnc is indeed a member of the star cluster, this system therefore is useful to constrain the formation and evolution of magnetic cataclysmic variables. The spectra exhibit two-component… ▽ More We present time-resolved spectroscopic and polarimetric observations of the AM Her system EU Cnc. EU Cnc is located near the core of the old open cluster Messier 67; new proper motion measurements indicate that EU Cnc is indeed a member of the star cluster, this system therefore is useful to constrain the formation and evolution of magnetic cataclysmic variables. The spectra exhibit two-component emission features with independent radial velocity variations as well as time-variable cyclotron emission indicating a magnetic field strength of 41 MG. The period of the radial velocity and cyclotron hump variations are consistent with the previously-known photometric period, and the spectroscopic flux variations are consistent in amplitude with previous photometric amplitude measurements. The secondary star is also detected in the spectrum. We also present polarimetric imaging measurements of EU Cnc that show a clear detection of polarization, and the degree of polarization drops below our detection threshold at phases when the cyclotron emission features are fading or not evident. The combined data are all consistent with the interpretation that EU Cnc is a low-state polar in the cluster Messier 67. The mass function of the system gives an estimate of the accretor mass of M_WD >= 0.68 M_sun with M_WD ~ 0.83 M_sun for an average inclination. We are thus able to place a lower limit on the progenitor mass of the accreting WD of >= 1.43 M_sun. △ Less

Submitted 24 January, 2013; originally announced January 2013.

Comments: Accepted for publication in the Astronomical Journal. 9 pages, 5 figures, uses emulateapj

Journal ref: The Astronomical Journal 145 (2013) 129

Rapid TeV Gamma-Ray Flaring of BL Lacertae

Authors: T. Arlen, T. Aune, M. Beilicke, W. Benbow, A. Bouvier, J. H. Buckley, V. Bugaev, A. Cesarini, L. Ciupik, M. P. Connolly, W. Cui, R. Dickherber, J. Dumm, M. Errando, A. Falcone, S. Federici, Q. Feng, J. P. Finley, G. Finnegan, L. Fortson, A. Furniss, N. Galante, D. Gall, S. Griffin, J. Grube , et al. (61 additional authors not shown)

Abstract: We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measur… ▽ More We report on the detection of a very rapid TeV gamma-ray flare from BL Lacertae on 2011 June 28 with the Very Energetic Radiation Imaging Telescope Array System (VERITAS). The flaring activity was observed during a 34.6-minute exposure, when the integral flux above 200 GeV reached $(3.4\pm0.6) \times 10^{-6} \;\text{photons}\;\text{m}^{-2}\text{s}^{-1}$, roughly 125% of the Crab Nebula flux measured by VERITAS. The light curve indicates that the observations missed the rising phase of the flare but covered a significant portion of the decaying phase. The exponential decay time was determined to be $13\pm4$ minutes, making it one of the most rapid gamma-ray flares seen from a TeV blazar. The gamma-ray spectrum of BL Lacertae during the flare was soft, with a photon index of $3.6\pm 0.4$, which is in agreement with the measurement made previously by MAGIC in a lower flaring state. Contemporaneous radio observations of the source with the Very Long Baseline Array (VLBA) revealed the emergence of a new, superluminal component from the core around the time of the TeV gamma-ray flare, accompanied by changes in the optical polarization angle. Changes in flux also appear to have occurred at optical, UV, and GeV gamma-ray wavelengths at the time of the flare, although they are difficult to quantify precisely due to sparse coverage. A strong flare was seen at radio wavelengths roughly four months later, which might be related to the gamma-ray flaring activities. We discuss the implications of these multiwavelength results. △ Less

Submitted 13 November, 2012; originally announced November 2012.

Comments: Accepted for publication in The Astrophysical Journal

Variability of the blazar 4C 38.41 (B3 1633+382) from GHz frequencies to GeV energies

Authors: C. M. Raiteri, M. Villata, P. S. Smith, V. M. Larionov, J. A. Acosta-Pulido, M. F. Aller, F. D'Ammando, M. A. Gurwell, S. G. Jorstad, M. Joshi, O. M. Kurtanidze, A. Lähteenmäki, D. O. Mirzaqulov, I. Agudo, H. D. Aller, M. J. Arévalo, A. A. Arkharov, U. Bach, E. Benítez, A. Berdyugin, D. A. Blinov, K. Blumenthal, C. S. Buemi, A. Bueno, T. M. Carleton , et al. (52 additional authors not shown)

Abstract: The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energ… ▽ More The quasar-type blazar 4C 38.41 (B3 1633+382) experienced a large outburst in 2011, which was detected throughout the entire electromagnetic spectrum. We present the results of low-energy multifrequency monitoring by the GASP project of the WEBT consortium and collaborators, as well as those of spectropolarimetric/spectrophotometric monitoring at the Steward Observatory. We also analyse high-energy observations of the Swift and Fermi satellites. In the optical-UV band, several results indicate that there is a contribution from a QSO-like emission component, in addition to both variable and polarised jet emission. The unpolarised emission component is likely thermal radiation from the accretion disc that dilutes the jet polarisation. We estimate its brightness to be R(QSO) ~ 17.85 - 18 and derive the intrinsic jet polarisation degree. We find no clear correlation between the optical and radio light curves, while the correlation between the optical and γ-ray flux apparently fades in time, likely because of an increasing optical to γ-ray flux ratio. As suggested for other blazars, the long-term variability of 4C 38.41 can be interpreted in terms of an inhomogeneous bent jet, where different emitting regions can change their alignment with respect to the line of sight, leading to variations in the Doppler factor δ. Under the hypothesis that in the period 2008-2011 all the γ-ray and optical variability on a one-week timescale were due to changes in δ, this would range between ~ 7 and ~ 21. If the variability were caused by changes in the viewing angle θ only, then θ would go from ~ 2.6 degr to ~ 5 degr. Variations in the viewing angle would also account for the dependence of the polarisation degree on the source brightness in the framework of a shock-in-jet model. △ Less

Submitted 17 July, 2012; originally announced July 2012.

Comments: 19 pages, 23 figures, in press for Astronomy and Astrophysics

Journal ref: A&A 545, A48 (2012)