Showing 1–43 of 43 results for author: Reeves, A

Association of the IceCube neutrinos with blazars in the CGRaBS sample

Authors: Pouya M. Kouch, Elina Lindfors, Talvikki Hovatta, Ioannis Liodakis, Karri I. I. Koljonen, Kari Nilsson, Sebastian Kiehlmann, Walter Max-Moerbeck, Anthony C. S. Readhead, Rodrigo A. Reeves, Timothy J. Pearson, Jenni Jormanainen, Vandad Fallah Ramazani, Matthew J. Graham

Abstract: The origin of high-energy (HE) astrophysical neutrinos has remained an elusive hot topic in the field of HE astrophysics for the past decade. Apart from a handful of individual associations, the vast majority of HE neutrinos arise from unknown sources. While there are theoretically-motivated candidate populations, such as blazars -- a subclass of AGN with jets pointed towards our line-of-sight --… ▽ More The origin of high-energy (HE) astrophysical neutrinos has remained an elusive hot topic in the field of HE astrophysics for the past decade. Apart from a handful of individual associations, the vast majority of HE neutrinos arise from unknown sources. While there are theoretically-motivated candidate populations, such as blazars -- a subclass of AGN with jets pointed towards our line-of-sight -- they have not yet been convincingly linked to HE neutrino production. Here, we perform a spatio-temporal association analysis between a sample of blazars (from CGRaBS catalog) in the radio and optical bands and the most up-to-date IceCube HE neutrino catalog. We find that if the IceCube error regions are enlarged by 1$^\circ$ in quadrature, to account for unknown systematic errors at maximal level, a spatio-temporal correlation between the multiwavelength light curves of the CGRaBS blazars and the IceCube HE neutrinos is hinted at least at a 2.17$σ$ significance level. On the other hand, when the IceCube error regions are taken as their published values, we do not find any significant correlations. A discrepancy in the blazar-neutrino correlation strengths, when using such minimal and enlarged error region scenarios, was also obtained in a recent study by the IceCube collaboration. In our study, this difference arises because several flaring blazars -- coinciding with a neutrino arrival time -- happen to narrowly miss the published 90\%-likelihood error region of the nearest neutrino event. For all of the associations driving our most significant correlations, the flaring blazar is much less than 1$^\circ$ away from the published error regions. Therefore, our results indicate that the question of the blazar-neutrino connection is highly sensitive to the reconstruction of the neutrino error regions, whose reliability is expected to improve with the next generation of neutrino observatories. △ Less

Submitted 9 July, 2024; originally announced July 2024.

Comments: Accepted for publication in A&A, 24 pages, 10 figures, 2 tables, 3 electronic tables

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

Unprecedented extreme high-frequency radio variability in early-stage active galactic nuclei

Authors: E. Järvelä, T. Savolainen, M. Berton, A. Lähteenmäki, S. Kiehlmann, T. Hovatta, I. Varglund, A. C. S. Readhead, M. Tornikoski, W. Max-Moerback, R. A. Reeves, S. Suutarinen

Abstract: We report on the discovery of one of the most extreme cases of high-frequency radio variability ever measured in active galactic nuclei (AGN), observed on timescales of days and exhibiting variability amplitudes of three to four orders of magnitude. These sources, all radio-weak narrow-line Seyfert 1 (NLS1) galaxies, were discovered some years ago at Aalto University Metsähovi Radio Observatory (M… ▽ More We report on the discovery of one of the most extreme cases of high-frequency radio variability ever measured in active galactic nuclei (AGN), observed on timescales of days and exhibiting variability amplitudes of three to four orders of magnitude. These sources, all radio-weak narrow-line Seyfert 1 (NLS1) galaxies, were discovered some years ago at Aalto University Metsähovi Radio Observatory (MRO) based on recurring flaring at 37 GHz, strongly indicating the presence of relativistic jets. In subsequent observations with the Karl G. Jansky Very Large Array (JVLA) at 1.6, 5.2, and 9.0~GHz no signs of jets were seen. To determine the cause of their extraordinary behaviour, we observed them with the JVLA at 10, 15, 22, 33, and 45 GHz, and with the Very Long Baseline Array (VLBA) at 15 GHz. These observations were complemented with single-dish monitoring at 37 GHz at MRO, and at 15 GHz at Owens Valley Radio Observatory (OVRO). Intriguingly, all but one source either have a steep radio spectrum up to 45 GHz, or were not detected at all. Based on the 37 GHz data the timescales of the radio flares are a few days, and the derived variability brightness temperatures and variability Doppler factors comparable to those seen in blazars. We discuss alternative explanations for their extreme behaviour, but so far no definite conclusions can be made. These sources exhibit radio variability at a level rarely, if ever, seen in AGN. They might represent a new type of jetted AGN, or a new variability phenomenon, and thus deserve our continued attention. △ Less

Submitted 4 December, 2023; originally announced December 2023.

Comments: 31 pages, 37 figures, submitted to MNRAS

How many stars form in galaxy mergers?

Authors: Andrew M. M. Reeves, Michael J. Hudson

Abstract: We forward model the difference in stellar age between post-coalescence mergers and a control sample with the same stellar mass, environmental density, and redshift. In particular, we use a pure sample of 445 post-coalescence mergers from the recent visually-confirmed post-coalescence merger sample identified by Bickley et al. and find that post-coalescence mergers are on average younger than cont… ▽ More We forward model the difference in stellar age between post-coalescence mergers and a control sample with the same stellar mass, environmental density, and redshift. In particular, we use a pure sample of 445 post-coalescence mergers from the recent visually-confirmed post-coalescence merger sample identified by Bickley et al. and find that post-coalescence mergers are on average younger than control galaxies for $10<\log (M_\star/\mathrm{M}_\odot)<11$. The difference in age from matched controls is up to 1.5 Gyr, highest for lower stellar mass galaxies. We forward model this difference using parametric star formation histories, accounting for the pre-coalescence inspiral phase of enhanced star formation using close pair data, and a final additive burst of star formation at coalescence. We find a best-fitting stellar mass burst fraction of $f_\mathrm{burst}=ΔM_\star/M_{\star,\mathrm{merger}}=0.18 \pm 0.02$ for $10<\log (M_\star/\mathrm{M}_\odot)<11$ galaxies, with no evidence of a trend in stellar mass. The modeled burst fraction is robust to choice of parametric star formation history, as well as differences in burst duration. The result appears consistent with some prior observationally-derived values, but is significantly higher than that found in hydrodynamical simulations. Using published Luminous InfraRed Galaxy (LIRG) star formation rates, we find a burst duration increasing with stellar mass, from $120-250$ Myr. A comparison to published cold gas measurements indicates there is enough molecular gas available in very close pairs to fuel the burst. Additionally, given our stellar mass burst estimate, the predicted cold gas fraction remaining after the burst is consistent with observed post-coalescence mergers. △ Less

Submitted 17 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS; 12 pages, 9 figures

CLASS Observations of Atmospheric Cloud Polarization at Millimeter Wavelengths

Authors: Yunyang Li, John W. Appel, Charles L. Bennett, Ricardo Bustos, David T. Chuss, Joseph Cleary, Jullianna Denes Couto, Sumit Dahal, Rahul Datta, Rolando Dünner, Joseph R. Eimer, Thomas Essinger-Hileman, Kathleen Harrington, Jeffrey Iuliano, Tobias A. Marriage, Matthew A. Petroff, Rodrigo A. Reeves, Karwan Rostem, Rui Shi, Deniz A. N. Valle, Duncan J. Watts, Oliver F. Wolff, Edward J. Wollack, Zhilei Xu

Abstract: The dynamic atmosphere imposes challenges to ground-based cosmic microwave background observation, especially for measurements on large angular scales. The hydrometeors in the atmosphere, mostly in the form of clouds, scatter the ambient thermal radiation and are known to be the main linearly polarized source in the atmosphere. This scattering-induced polarization is significantly enhanced for ice… ▽ More The dynamic atmosphere imposes challenges to ground-based cosmic microwave background observation, especially for measurements on large angular scales. The hydrometeors in the atmosphere, mostly in the form of clouds, scatter the ambient thermal radiation and are known to be the main linearly polarized source in the atmosphere. This scattering-induced polarization is significantly enhanced for ice clouds due to the alignment of ice crystals under gravity, which are also the most common clouds seen at the millimeter-astronomy sites at high altitudes. This work presents a multifrequency study of cloud polarization observed by the Cosmology Large Angular Scale Surveyor (CLASS) experiment on Cerro Toco in the Atacama Desert of northern Chile, from 2016 to 2022, at the frequency bands centered around 40, 90, 150, and 220 GHz. Using a machine-learning-assisted cloud classifier, we made connections between the transient polarized emission found in all four frequencies with the clouds imaged by monitoring cameras at the observing site. The polarization angles of the cloud events are found to be mostly $90^\circ$ from the local meridian, which is consistent with the presence of horizontally aligned ice crystals. The 90 and 150 GHz polarization data are consistent with a power law with a spectral index of $3.90\pm0.06$, while an excess/deficit of polarization amplitude is found at 40/220 GHz compared with a Rayleigh scattering spectrum. These results are consistent with Rayleigh-scattering-dominated cloud polarization, with possible effects from supercooled water absorption and/or Mie scattering from a population of large cloud particles that contribute to the 220 GHz polarization. △ Less

Submitted 13 September, 2023; originally announced September 2023.

Comments: 16 pages, 14 figures, submitted to ApJ

$\mathbf{12\times2}$pt combined probes: pipeline, neutrino mass, and data compression

Authors: Alexander Reeves, Andrina Nicola, Alexandre Refregier, Tomasz Kacprzak, Luis Fernando Machado Poletti Valle

Abstract: With the rapid advance of wide-field surveys it is increasingly important to perform combined cosmological probe analyses. We present a new pipeline for simulation-based multi-probe analyses, which combines tomographic large-scale structure (LSS) probes (weak lensing and galaxy clustering) with cosmic microwave background (CMB) primary and lensing data. These are combined at the $C_\ell$-level, yi… ▽ More With the rapid advance of wide-field surveys it is increasingly important to perform combined cosmological probe analyses. We present a new pipeline for simulation-based multi-probe analyses, which combines tomographic large-scale structure (LSS) probes (weak lensing and galaxy clustering) with cosmic microwave background (CMB) primary and lensing data. These are combined at the $C_\ell$-level, yielding 12 distinct auto- and cross-correlations. The pipeline is based on $\texttt{UFalconv2}$, a framework to generate fast, self-consistent map-level realizations of cosmological probes from input lightcones, which is applied to the $\texttt{CosmoGridV1}$ N-body simulation suite. It includes a non-Gaussian simulation-based covariance for the LSS tracers, several data compression schemes, and a neural network emulator for accelerated theoretical predictions. We validate our framework, apply it to a simulated $12\times2$pt tomographic analysis of KiDS, BOSS, and $\textit{Planck}$, and forecast constraints for a $Λ$CDM model with a variable neutrino mass. We find that, while the neutrino mass constraints are driven by the CMB data, the addition of LSS data helps to break degeneracies and improves the constraint by up to 35%. For a fiducial $M_ν=0.15\mathrm{eV}$, a full combination of the above CMB+LSS data would enable a $3σ$ constraint on the neutrino mass. We explore data compression schemes and find that MOPED outperforms PCA. We also study the impact of an internal lensing tension in the CMB data, parametrized by $A_L$, on the neutrino mass constraint, finding that the addition of LSS to CMB data including all cross-correlations is able to mitigate the impact of this systematic. $\texttt{UFalconv2}$ and a MOPED compressed $\textit{Planck}$ CMB primary + CMB lensing likelihood are made publicly available. [abridged] △ Less

Submitted 6 December, 2023; v1 submitted 6 September, 2023; originally announced September 2023.

Comments: Accepted for publication in JCAP. 51 pages including 7 Tables and 17 Figures. Associated $\texttt{UFalconv2}$ code available at https://cosmology.ethz.ch/research/software-lab/UFalcon.html, compressed $\textit{Planck}$ likelihood at https://github.com/alexreevesy/planck_compressed

CLASS Angular Power Spectra and Map-Component Analysis for 40 GHz Observations through 2022

Authors: Joseph R. Eimer, Yunyang Li, Michael K. Brewer, Rui Shi, Aamir Ali, John W. Appel, Charles L. Bennett, Sarah Marie Bruno, Ricardo Bustos, David T. Chuss, Joseph Cleary, Sumit Dahal, Rahul Datta, Jullianna Denes Couto, Kevin L. Denis, Rolando Dünner, Thomas Essinger-Hileman, Pedro Fluxá, Johannes Hubmayer, Kathleen Harrington, Jeffrey Iuliano, John Karakla, Tobias A. Marriage, Carolina Núñez, Lucas Parker , et al. (9 additional authors not shown)

Abstract: Measurement of the largest angular scale ($\ell < 30$) features of the cosmic microwave background (CMB) polarization is a powerful way to constrain the optical depth to reionization and search for the signature of inflation through the detection of primordial $B$-modes. We present an analysis of maps covering 73.6\% of the sky made from the $40\,\mathrm{GHz}$ channel of the Cosmology Large Angula… ▽ More Measurement of the largest angular scale ($\ell < 30$) features of the cosmic microwave background (CMB) polarization is a powerful way to constrain the optical depth to reionization and search for the signature of inflation through the detection of primordial $B$-modes. We present an analysis of maps covering 73.6\% of the sky made from the $40\,\mathrm{GHz}$ channel of the Cosmology Large Angular Scale Surveyor (CLASS) from 2016 August to 2022 May. Taking advantage of the measurement stability enabled by front-end polarization modulation and excellent conditions from the Atacama Desert, we show this channel achieves higher sensitivity than the analogous frequencies from satellite measurements in the range $10 < \ell < 100$. Simulations show the CLASS linear (circular) polarization maps have a white noise level of $125 \,(130)\,\mathrm{μK\, arcmin}$. We measure the Galaxy-masked $EE$ and $BB$ spectra of diffuse synchrotron radiation and compare to space-based measurements at similar frequencies. In combination with external data, we expand measurements of the spatial variations of the synchrotron spectral energy density (SED) to include new sky regions and measure the diffuse SED in the harmonic domain. We place a new upper limit on a background of circular polarization in the range $5 < \ell < 125$ with the first bin showing $D_\ell < 0.023$ $\mathrm{μK^2_{CMB}}$ at 95\% confidence. These results establish a new standard for recovery of the largest-scale CMB polarization from the ground and signal exciting possibilities when the higher sensitivity and higher-frequency CLASS channels are included in the analysis. △ Less

Submitted 14 February, 2024; v1 submitted 1 September, 2023; originally announced September 2023.

Comments: 38 pages, 25 figures, 6 tables. Accepted to The Astrophysical Journal

When the Well Runs Dry: Modeling Environmental Quenching of High-mass Satellites in Massive Clusters at \boldmath$z \gtrsim 1$

Authors: Devontae C. Baxter, Michael C. Cooper, Michael L. Balogh, Gregory H. Rudnick, Gabriella De Lucia, Ricardo Demarco, Alexis Finoguenov, Ben Forrest, Adam Muzzin, Andrew Reeves, Florian Sarron, Benedetta Vulcani, Gillian Wilson, Dennis Zaritsky

Abstract: We explore models of massive ($\gt 10^{10}~{\rm M}_{\odot}$) satellite quenching in massive clusters at $z\gtrsim1$ using an MCMC framework, focusing on two primary parameters: $R_{\rm quench}$ (the host-centric radius at which quenching begins) and $τ_{\rm quench}$ (the timescale upon which a satellite quenches after crossing $R_{\rm quench}$). Our MCMC analysis shows two local maxima in the 1D p… ▽ More We explore models of massive ($\gt 10^{10}~{\rm M}_{\odot}$) satellite quenching in massive clusters at $z\gtrsim1$ using an MCMC framework, focusing on two primary parameters: $R_{\rm quench}$ (the host-centric radius at which quenching begins) and $τ_{\rm quench}$ (the timescale upon which a satellite quenches after crossing $R_{\rm quench}$). Our MCMC analysis shows two local maxima in the 1D posterior probability distribution of $R_{\rm quench}$ at approximately $0.25$ and $1.0~R_{\rm{200}}$. Analyzing four distinct solutions in the $τ_{\rm quench}$-$R_{\rm quench}$ parameter space, nearly all of which yield quiescent fractions consistent with observational data from the GOGREEN survey, we investigate whether these solutions represent distinct quenching pathways and find that they can be separated between \textquote{starvation} and \textquote{core quenching} scenarios. The starvation pathway is characterized by quenching timescales that are roughly consistent with the total cold gas (H$_{2}$+H{\scriptsize I}) depletion timescale at intermediate $z$, while core quenching is characterized by satellites with relatively high line-of-sight velocities that quench on short timescales ($\sim 0.25$ Gyr) after reaching the inner region of the cluster ($\lt 0.30~R_{\rm{200}}$). Lastly, we break the degeneracy between these solutions by comparing the observed properties of transition galaxies from the GOGREEN survey. We conclude that only the \textquote{starvation} pathway is consistent with the projected phase-space distribution and relative abundance of transition galaxies at $z \sim 1$. However, we acknowledge that ram pressure might contribute as a secondary quenching mechanism. △ Less

Submitted 17 October, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 15 pages; 8 figures; Accepted for publication in Monthly Notices of the Royal Astronomical Society

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 526, Issue 3, December 2023

Long-term multi-wavelength study of 1ES 0647+250

Authors: MAGIC Collaboration, V. A. Acciari, 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, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch , et al. (195 additional authors not shown)

Abstract: The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-in… ▽ More The BL Lac object 1ES 0647+250 is one of the few distant $γ$-ray emitting blazars detected at very high energies (VHE, $\gtrsim$100 GeV) during a non-flaring state. It was detected with the MAGIC telescopes during its low activity in the years 2009-2011, as well as during three flaring activities in the years 2014, 2019 and 2020, with the highest VHE flux in the latter epoch. An extensive multi-instrument data set was collected within several coordinated observing campaigns throughout these years. We aim to characterise the long-term multi-band flux variability of 1ES 0647+250, as well as its broadband spectral energy distribution (SED) during four distinct activity states selected in four different epochs, in order to constrain the physical parameters of the blazar emission region under certain assumptions. We evaluate the variability and correlation of the emission in the different energy bands with the fractional variability and the Z-transformed Discrete Correlation Function, as well as its spectral evolution in X-rays and $γ$ rays. Owing to the controversy in the redshift measurements of 1ES 0647+250 reported in the literature, we also estimate its distance in an indirect manner through the comparison of the GeV and TeV spectra from simultaneous observations with Fermi-LAT and MAGIC during the strongest flaring activity detected to date. Moreover, we interpret the SEDs from the four distinct activity states within the framework of one-component and two-component leptonic models, proposing specific scenarios that are able to reproduce the available multi-instrument data. △ Less

Submitted 23 November, 2022; originally announced November 2022.

Comments: 20 pages, 7 figures. Accepted in A&A. Corresponding authors: Jorge Otero-Santos; Daniel Morcuende; Vandad Fallah Ramazani; Daniela Dorner; David Paneque (mailto: contact.magic@mpp.mpg.de)

Journal ref: A&A 670, A49 (2023)

Constraining quenching timescales in galaxy clusters by forward-modelling stellar ages and quiescent fractions in projected phase space

Authors: Andrew M. M. Reeves, Michael J. Hudson, Kyle A. Oman

Abstract: We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in $\log(M_{\mathrm{vir}}/\mathrm{M}_{\odot})>14$ galaxy clusters at $z\sim 0$. We find the average deviation in MWA from the MWA-$M_\star$ relation depends on position in PPS, with a max… ▽ More We forward-model mass-weighted stellar ages (MWAs) and quiescent fractions in projected phase space (PPS), using data from the Sloan Digital Sky Survey, to jointly constrain an infall quenching model for galaxies in $\log(M_{\mathrm{vir}}/\mathrm{M}_{\odot})>14$ galaxy clusters at $z\sim 0$. We find the average deviation in MWA from the MWA-$M_\star$ relation depends on position in PPS, with a maximum difference between the inner cluster and infalling interloper galaxies of $\sim 1$ Gyr. Our model employs infall information from N-body simulations and stochastic star-formation histories from the UniverseMachine model. We find total quenching times of $t_\mathrm{Q}=3.7\pm 0.4$ Gyr and $t_\mathrm{Q}=4.0\pm 0.2$ Gyr after first pericentre, for $9<\log(M_{\star}/\mathrm{M}_{\odot})<10$ and $10<\log(M_{\star}/\mathrm{M}_{\odot})<10.5$ galaxies, respectively. By using MWAs, we break the degeneracy in time of quenching onset and timescale of star formation rate (SFR) decline. We find that time of quenching onset relative to pericentre is $t_{\mathrm{delay}}=3.5^{+0.6}_{-0.9}$ Gyr and $t_{\mathrm{delay}}=-0.3^{+0.8}_{-1.0}$ Gyr for our lower and higher stellar mass bins, respectively, and exponential SFR suppression timescales are $τ_{\mathrm{env}}\leq 1.0$ Gyr and $τ_{\mathrm{env}}\sim 2.3$ Gyr for our lower and higher stellar mass bins, respectively. Stochastic star formation histories remove the need for rapid infall quenching to maintain the bimodality in the SFR of cluster galaxies; the depth of the green valley prefers quenching onsets close to first pericentre and a longer quenching envelope, in slight tension with the MWA-driven results. Taken together these results suggest that quenching begins close to, or just after pericentre, but the timescale for quenching to be fully complete is much longer and therefore ram-pressure stripping is not complete on first pericentric passage. △ Less

Submitted 7 June, 2024; v1 submitted 16 November, 2022; originally announced November 2022.

Comments: Accepted for publication in MNRAS (April 1 2023), 21 pages, 13 figures

GOGREEN: a critical assessment of environmental trends in cosmological hydrodynamical simulations at z ~ 1

Authors: Egidijus Kukstas, Michael L. Balogh, Ian G. McCarthy, Yannick M. Bahe, Gabriella De Lucia, Pascale Jablonka, Benedetta Vulcani, Devontae C. Baxter, Andrea Biviano, Pierluigi Cerulo, Jeffrey C. Chan, M. C. Cooper, Ricardo Demarco, Alexis Finoguenov, Andreea S. Font, Chris Lidman, Justin Marchioni, Sean McGee, Adam Muzzin, Julie Nantais, Lyndsay Old, Irene Pintos-Castro, Bianca Poggianti, Andrew M. M. Reeves, Gregory Rudnick , et al. (6 additional authors not shown)

Abstract: Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been… ▽ More Recent observations have shown that the environmental quenching of galaxies at z ~ 1 is qualitatively different to that in the local Universe. However, the physical origin of these differences has not yet been elucidated. In addition, while low-redshift comparisons between observed environmental trends and the predictions of cosmological hydrodynamical simulations are now routine, there have been relatively few comparisons at higher redshifts to date. Here we confront three state-of-the-art suites of simulations (BAHAMAS+MACSIS, EAGLE+Hydrangea, IllustrisTNG) with state-of-the-art observations of the field and cluster environments from the COSMOS/UltraVISTA and GOGREEN surveys, respectively, at z ~ 1 to assess the realism of the simulations and gain insight into the evolution of environmental quenching. We show that while the simulations generally reproduce the stellar content and the stellar mass functions of quiescent and star-forming galaxies in the field, all the simulations struggle to capture the observed quenching of satellites in the cluster environment, in that they are overly efficient at quenching low-mass satellites. Furthermore, two of the suites do not sufficiently quench the highest-mass galaxies in clusters, perhaps a result of insufficient feedback from AGN. The origin of the discrepancy at low stellar masses (Mstar <~ 1E10 Msun), which is present in all the simulations in spite of large differences in resolution, feedback implementations, and hydrodynamical solvers, is unclear. The next generation of simulations, which will push to significantly higher resolution and also include explicit modelling of the cold interstellar medium, may help to shed light on the low-mass tension. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 20 pages, 10 figures, resubmitted to MNRAS following minor revisions

Long-Timescale Stability in CMB Observations at Multiple Frequencies using Front-End Polarization Modulation

Authors: Joseph Cleary, Rahul Datta, John W. Appel, Charles L. Bennett, David T. Chuss, Jullianna Denes Couto, Sumit Dahal, Francisco Espinoza, Thomas Essinger-Hileman, Kathleen Harrington, Jeffrey Iuliano, Yunyang Li, Tobias A. Marriage, Carolina Nunez, Matthew A. Petroff, Rodrigo A. Reeves, Rui Shi, Duncan J. Watts, Edward J. Wollack, Zhilei Xu

Abstract: The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array observing the Cosmic Microwave Background (CMB) at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the CMB polarization on the largest angular scales to constrain the inflationary tensor-to-scalar ratio and the optical depth due to reionization. To achieve the long time-scale stability necessary for this… ▽ More The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array observing the Cosmic Microwave Background (CMB) at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the CMB polarization on the largest angular scales to constrain the inflationary tensor-to-scalar ratio and the optical depth due to reionization. To achieve the long time-scale stability necessary for this measurement from the ground, CLASS utilizes a front-end, variable-delay polarization modulator on each telescope. Here we report on the improvements in stability afforded by front-end modulation using data across all four CLASS frequencies. Across one month of modulated linear polarization data in 2021, CLASS achieved median knee frequencies of 9.1, 29.1, 20.4, and 36.4 mHz for the 40, 90, 150, and 220 GHz observing bands. The knee frequencies are approximately an order of magnitude lower than achieved via CLASS pair-differencing orthogonal detector pairs without modulation. △ Less

Submitted 9 August, 2022; originally announced August 2022.

Comments: Submitted to SPIE Astronomical Telescopes + Instrumentation 2022 Conference (AS22)

The GOGREEN Survey: Constraining the Satellite Quenching Timescale in Massive Clusters at $\boldsymbol{z} \gtrsim 1$

Authors: Devontae Baxter, Michael Cooper, Michael Balogh, Tim Carleton, Pierluigi Cerulo, Gabriella De Lucia, Ricardo Demarco, Sean McGee, Adam Muzzin, Julie Nantais, Irene Pintos Castro, Andrew Reeves, Gregory Rudnick, Florian Sarron, Remco van der Burg, Benedetta Vulcani, Gillian Wilson, Dennis Zaritsky

Abstract: We model satellite quenching at $z \sim 1$ by combining $14$ massive ($10^{13.8} < M_{\mathrm{halo}}/\mathrm{M}_{\odot} < 10^{15}$) clusters at $0.8 < z < 1.3$ from the GOGREEN and GCLASS surveys with accretion histories of $56$ redshift-matched analogs from the IllustrisTNG simulation. Our fiducial model, which is parameterized by the satellite quenching timescale ($τ_{\rm quench}$), accounts for… ▽ More We model satellite quenching at $z \sim 1$ by combining $14$ massive ($10^{13.8} < M_{\mathrm{halo}}/\mathrm{M}_{\odot} < 10^{15}$) clusters at $0.8 < z < 1.3$ from the GOGREEN and GCLASS surveys with accretion histories of $56$ redshift-matched analogs from the IllustrisTNG simulation. Our fiducial model, which is parameterized by the satellite quenching timescale ($τ_{\rm quench}$), accounts for quenching in our simulated satellite population both at the time of infall by using the observed coeval field quenched fraction and after infall by tuning $τ_{\rm quench}$ to reproduce the observed satellite quenched fraction versus stellar mass trend. This model successfully reproduces the observed satellite quenched fraction as a function of stellar mass (by construction), projected cluster-centric radius, and redshift and is consistent with the observed field and cluster stellar mass functions at $z \sim 1$. We find that the satellite quenching timescale is mass dependent, in conflict with some previous studies at low and intermediate redshift. Over the stellar mass range probed ($M_{\star}> 10^{10}~\mathrm{M}_{\odot}$), we find that the satellite quenching timescale decreases with increasing satellite stellar mass from $\sim1.6~{\rm Gyr}$ at $10^{10}~\mathrm{M}_{\odot}$ to $\sim 0.6 - 1~{\rm Gyr}$ at $10^{11}~\mathrm{M}_{\odot}$ and is roughly consistent with the total cold gas (H{\scriptsize I}+H$_{2}$) depletion timescales at intermediate $z$, suggesting that starvation may be the dominant driver of environmental quenching at $z < 2$. Finally, while environmental mechanisms are relatively efficient at quenching massive satellites, we find that the majority ($\sim65-80\%$) of ultra-massive satellites ($M_{\star} > 10^{11}~\mathrm{M}_{\odot}$) are quenched prior to infall. △ Less

Submitted 28 July, 2022; originally announced July 2022.

Comments: Key Figures: 4, 6 and 9; 16 pages; 11 figures; accepted for publication in MNRAS 26 July 2022;

Restoring cosmological concordance with early dark energy and massive neutrinos?

Authors: Alexander Reeves, Laura Herold, Sunny Vagnozzi, Blake D. Sherwin, Elisa G. M. Ferreira

Abstract: The early dark energy (EDE) solution to the Hubble tension comes at the cost of an increased clustering amplitude that has been argued to worsen the fit to galaxy clustering data. We explore whether freeing the total neutrino mass $M_ν$, which can suppress small-scale structure growth, improves EDE's fit to galaxy clustering. Using Planck Cosmic Microwave Background and BOSS galaxy clustering data… ▽ More The early dark energy (EDE) solution to the Hubble tension comes at the cost of an increased clustering amplitude that has been argued to worsen the fit to galaxy clustering data. We explore whether freeing the total neutrino mass $M_ν$, which can suppress small-scale structure growth, improves EDE's fit to galaxy clustering. Using Planck Cosmic Microwave Background and BOSS galaxy clustering data, a Bayesian analysis shows that freeing $M_ν$ does not appreciably increase the inferred EDE fraction $f_{\rm EDE}$: we find the 95% C.L. upper limits $f_{\rm EDE}<0.092$ and $M_ν<0.15\,{\rm eV}$. Similarly, in a frequentist profile likelihood setting (where our results support previous findings that prior volume effects are important), we find that the baseline EDE model (with $M_ν=0.06\,{\rm eV}$) provides the overall best fit. For instance, compared to baseline EDE, a model with $M_ν=0.24\,{\rm eV}$ maintains the same $H_0$(km/s/Mpc)=(70.08, 70.11, respectively) whilst decreasing $S_8$=(0.837, 0.826) to the $Λ$CDM level, but worsening the fit significantly by $Δχ^2=7.5$. For the datasets used, these results are driven not by the clustering amplitude, but by background modifications to the late-time expansion rate due to massive neutrinos, which worsen the fit to measurements of the BAO scale. △ Less

Submitted 26 January, 2023; v1 submitted 4 July, 2022; originally announced July 2022.

Comments: 9 pages, 6 figures. v3: included additional discussions and 3 figures which further strengthen our results and their interpretation, added further references, slight format change. Accepted for publication in MNRAS. and associated Jupyter notebook available at https://github.com/alexreevesy/ede_mnu

Journal ref: Mon. Not. Roy. Astron. Soc. 520 (2023) 3688

Investigating the blazar TXS 0506+056 through sharp multi-wavelength eyes during 2017-2019

Authors: MAGIC Collaboration, V. A. Acciari, T. Aniello, S. Ansoldi, L. A. Antonelli, A. Arbet Engels, M. Artero, K. Asano, D. Baack, A. Babić, A. Baquero, U. Barres de Almeida, J. A. Barrio, I. Batković, J. Becerra González, W. Bednarek, E. Bernardini, M. Bernardos, A. Berti, J. Besenrieder, W. Bhattacharyya, C. Bigongiari, A. Biland, O. Blanch, H. Bökenkamp , et al. (192 additional authors not shown)

Abstract: The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the Ic… ▽ More The blazar TXS 0506+056 got into the spotlight of the astrophysical community in September 2017, when a high-energy neutrino detected by IceCube (IceCube-170922A) was associated at the 3 $σ$ level to a $γ$-ray flare from this source. This multi-messenger photon-neutrino association remains, as per today, the most significant one ever observed. TXS 0506+056 was a poorly studied object before the IceCube-170922A event. To better characterize its broad-band emission, we organized a multi-wavelength campaign lasting 16 months (November 2017 to February 2019), covering the radio-band (Metsähovi, OVRO), the optical/UV (ASAS-SN, KVA, REM, Swift/UVOT), the X-rays (Swift/XRT, NuSTAR), the high-energy $γ$ rays (Fermi/LAT) and the very-high-energy (VHE) $γ$ rays (MAGIC). In $γ$ rays, the behaviour of the source was significantly different from the 2017 one: MAGIC observations show the presence of flaring activity during December 2018, while the source only shows an excess at the 4$σ$ level during the rest of the campaign (74 hours of accumulated exposure); Fermi/LAT observations show several short (days-to-week timescale) flares, different from the long-term brightening of 2017. No significant flares are detected at lower energies. The radio light curve shows an increasing flux trend, not seen in other wavelengths. We model the multi-wavelength spectral energy distributions in a lepto-hadronic scenario, in which the hadronic emission emerges as Bethe-Heitler and pion-decay cascade in the X-rays and VHE $γ$ rays. According to the model presented here, the December 2018 $γ$-ray flare was connected to a neutrino emission that was too brief and not bright enough to be detected by current neutrino instruments. △ Less

Submitted 1 May, 2022; v1 submitted 5 February, 2022; originally announced February 2022.

Comments: 18 pages, 6 figures; V2 to match published version, typo fixed

New Tests of Millilensing in the Blazar PKS 1413+135

Authors: A. L. Peirson, I. Liodakis, A. C. S. Readhead, M. L. Lister, E. S. Perlman, M. F. Aller, R. D. Blandford, K. J. B. Grainge, D. A. Green, M. A. Gurwell, M. W. Hodges, T. Hovatta, S. Kiehlmann, A. Lähteenmäki, W. Max-Moerbeck, T. Mcaloone, S. O'Neill, V. Pavlidou, T. J. Pearson, V. Ravi, R. A. Reeves, P. F. Scott, G. B. Taylor, D. J. Titterington, M. Tornikoski , et al. (4 additional authors not shown)

Abstract: Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If… ▽ More Symmetric Achromatic Variability (SAV) is a rare form of radio variability in blazars that has been attributed to gravitational millilensing by a ~$10^2 - 10^5$ $M_\odot$ mass condensate. Four SAVs have been identified between 1980 and 2020 in the long-term radio monitoring data of the blazar PKS 1413+135. We show that all four can be fitted with the same, unchanging, gravitational lens model. If SAV is due to gravitational millilensing, PKS 1413+135 provides a unique system for studying active galactic nuclei with unprecedented microarcsecond resolution, as well as for studying the nature of the millilens itself. We discuss two possible candidates for the putative millilens: a giant molecular cloud hosted in the intervening edge-on spiral galaxy, and an undetected dwarf galaxy with a massive black hole. We find a significant dependence of SAV crossing time on frequency, which could indicate a fast shock moving in a slower underlying flow. We also find tentative evidence for a 989-day periodicity in the SAVs, which, if real, makes possible the prediction of future SAVs: the next three windows for possible SAVs begin in August 2022, May 2025, and February 2028. △ Less

Submitted 8 January, 2022; v1 submitted 4 January, 2022; originally announced January 2022.

Comments: Accepted to ApJ

The Unanticipated Phenomenology of the Blazar PKS~2131$-$021: A Unique Super-Massive Black Hole Binary Candidate

Authors: S. O'Neill, S. Kiehlmann, A. C. S. Readhead, M. F. Aller, R. D. Blandford, I. Liodakis, M. L. Lister, P. Mróz, C. P. O'Dea, T. J. Pearson, V. Ravi, M. Vallisneri, K. A. Cleary, M. J. Graham, K. J. B. Grainge, M. W. Hodges, T. Hovatta, A. Lähteenmäki, J. W. Lamb, T. J. W. Lazio, W. Max-Moerbeck, V. Pavlidou, T. A. Prince, R. A. Reeves, M. Tornikoski , et al. (2 additional authors not shown)

Abstract: Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio lig… ▽ More Most large galaxies host supermassive black holes in their nuclei and are subject to mergers, which can produce a supermassive black hole binary (SMBHB), and hence periodic signatures due to orbital motion. We report unique periodic radio flux density variations in the blazar PKS~2131$-$021, which strongly suggest an SMBHB with an orbital separation of $\sim 0.001-0.01$ pc. Our 45.1-year radio light curve shows two epochs of strong sinusoidal variation with the same period and phase to within $<2\%$ and $\sim 10\%$, respectively, straddling a 20-year period when this variation was absent. Our simulated light curves accurately reproduce the ``red noise'' of this object, and Lomb-Scargle, weighted wavelet Z-transform, and least-squares sine wave analyses demonstrate conclusively, at the $4.6σ$ significance level, that the periodicity in this object is not due to random fluctuations in flux density. The observed period translates to $2.082\pm 0.003$ years in the rest frame at the $z=1.285$ redshift of PKS~2131$-$021. The periodic variation in PKS~2131$-$021 is remarkably sinusoidal. We present a model in which orbital motion, combined with the strong Doppler boosting of the approaching relativistic jet, produces a sine-wave modulation in the flux density which easily fits the observations. Given the rapidly-developing field of gravitational wave experiments with pulsar timing arrays, closer counterparts to PKS~2131$-$021 and searches using the techniques we have developed are strongly motivated. These results constitute a compelling demonstration that the phenomenology, not the theory, must provide the lead in this field. △ Less

Submitted 24 January, 2022; v1 submitted 3 November, 2021; originally announced November 2021.

Comments: 24 pages, 13 figure, 3 Tables, accepted for publication in APJL

The GOGREEN survey: Transition Galaxies and The Evolution of Environmental Quenching

Authors: Karen McNab, Michael L. Balogh, Remco F. J. van der Burg, Anya Forestell, Kristi Webb, Benedetta Vulcani, Gregory Rudnick, Adam Muzzin, M. C. Cooper, Sean McGee, Andrea Biviano, Pierluigi Cerulo, Jeffrey C. C. Chan, Gabriella De Lucia, Ricardo Demarco, Alexis Finoguenov, Ben Forrest, Caelan Golledge, Pascale Jablonka, Chris Lidman, Julie Nantais, Lyndsay Old, Irene Pintos-Castro, Bianca Poggianti, Andrew M. M. Reeves , et al. (3 additional authors not shown)

Abstract: We measure the rate of environmentally-driven star formation quenching in galaxies at $z\sim 1$, using eleven massive ($M\approx 2\times10^{14}\,\mathrm{M}_\odot$) galaxy clusters spanning a redshift range $1.0<z<1.4$ from the GOGREEN sample. We identify three different types of transition galaxies: "green valley" (GV) galaxies identified from their rest-frame $(NUV-V)$ and $(V-J)$ colours; "blue… ▽ More We measure the rate of environmentally-driven star formation quenching in galaxies at $z\sim 1$, using eleven massive ($M\approx 2\times10^{14}\,\mathrm{M}_\odot$) galaxy clusters spanning a redshift range $1.0<z<1.4$ from the GOGREEN sample. We identify three different types of transition galaxies: "green valley" (GV) galaxies identified from their rest-frame $(NUV-V)$ and $(V-J)$ colours; "blue quiescent" (BQ) galaxies, found at the blue end of the quiescent sequence in $(U-V)$ and $(V-J)$ colour; and spectroscopic post-starburst (PSB) galaxies. We measure the abundance of these galaxies as a function of stellar mass and environment. For high stellar mass galaxies ($\log{M/\mathrm{M}_\odot}>10.5$) we do not find any significant excess of transition galaxies in clusters, relative to a comparison field sample at the same redshift. It is likely that such galaxies were quenched prior to their accretion in the cluster, in group, filament or protocluster environments. For lower stellar mass galaxies ($9.5<\log{M/\mathrm{M}_\odot}<10.5$) there is a small but significant excess of transition galaxies in clusters, accounting for an additional $\sim 5-10$ per cent of the population compared with the field. We show that our data are consistent with a scenario in which 20--30 per cent of low-mass, star-forming galaxies in clusters are environmentally quenched every Gyr, and that this rate slowly declines from $z=1$ to $z=0$. While environmental quenching of these galaxies may include a long delay time during which star formation declines slowly, in most cases this must end with a rapid ($τ<1$ Gyr) decline in star formation rate. △ Less

Submitted 7 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS Sept 6, 2021

The GOGREEN survey: Dependence of galaxy properties on halo mass at z > 1 and implications for environmental quenching

Authors: Andrew M. M. Reeves, Michael L. Balogh, Remco F. J. van der Burg, Alexis Finoguenov, Egidijus Kukstas, Ian G. McCarthy, Kristi Webb, Adam Muzzin, Sean McGee, Gregory Rudnick, Andrea Biviano, Pierluigi Cerulo, Jeffrey C. C. Chan, M. C. Cooper, Ricardo Demarco, Pascale Jablonka, Gabriella De Lucia, Benedetta Vulcani, Gillian Wilson, Howard K. C. Yee, Dennis Zaritsky

Abstract: We use photometric redshifts and statistical background subtraction to measure stellar mass functions in galaxy group-mass ($4.5-8\times10^{13}~\mathrm{M}_\odot$) haloes at $1<z<1.5$. Groups are selected from COSMOS and SXDF, based on X-ray imaging and sparse spectroscopy. Stellar mass ($M_{\mathrm{stellar}}$) functions are computed for quiescent and star-forming galaxies separately, based on thei… ▽ More We use photometric redshifts and statistical background subtraction to measure stellar mass functions in galaxy group-mass ($4.5-8\times10^{13}~\mathrm{M}_\odot$) haloes at $1<z<1.5$. Groups are selected from COSMOS and SXDF, based on X-ray imaging and sparse spectroscopy. Stellar mass ($M_{\mathrm{stellar}}$) functions are computed for quiescent and star-forming galaxies separately, based on their rest-frame $UVJ$ colours. From these we compute the quiescent fraction and quiescent fraction excess (QFE) relative to the field as a function of $M_{\mathrm{stellar}}$. QFE increases with $M_{\mathrm{stellar}}$, similar to more massive clusters at $1<z<1.5$. This contrasts with the apparent separability of $M_{\mathrm{stellar}}$ and environmental factors on galaxy quiescent fractions at $z\sim 0$. We then compare our results with higher mass clusters at $1<z<1.5$ and lower redshifts. We find a strong QFE dependence on halo mass at fixed $M_{\mathrm{stellar}}$; well fit by a logarithmic slope of $\mathrm{d}(\mathrm{QFE})/\mathrm{d}\log (M_{\mathrm{halo}}) \sim 0.24 \pm 0.04$ for all $M_{\mathrm{stellar}}$ and redshift bins. This dependence is in remarkably good qualitative agreement with the hydrodynamic simulation BAHAMAS, but contradicts the observed dependence of QFE on $M_{\mathrm{stellar}}$. We interpret the results using two toy models: one where a time delay until rapid (instantaneous) quenching begins upon accretion to the main progenitor ("no pre-processing") and one where it starts upon first becoming a satellite ("pre-processing"). Delay times appear to be halo mass dependent, with a significantly stronger dependence required without pre-processing. We conclude that our results support models in which environmental quenching begins in low-mass ($<10^{14}M_\odot$) haloes at $z>1$. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: Accepted July 6, 2021, MNRAS

The GOGREEN Survey: Evidence of an excess of quiescent disks in clusters at $1.0<z<1.4$

Authors: Jeffrey C. C. Chan, Gillian Wilson, Michael Balogh, Gregory Rudnick, Remco F. J. van der Burg, Adam Muzzin, Kristi A. Webb, Andrea Biviano, Pierluigi Cerulo, M. C. Cooper, Gabriella De Lucia, Ricardo Demarco, Ben Forrest, Pascale Jablonka, Chris Lidman, Sean L. McGee, Julie Nantais, Lyndsay Old, Irene Pintos-Castro, Bianca Poggianti$, Andrew M. M. Reeves, Benedetta Vulcani, Howard K. C. Yee, Dennis Zaritsky

Abstract: We present results on the measured shapes of 832 galaxies in 11 galaxy clusters at 1.0 < z <1.4 from the GOGREEN survey. We measure the axis ratio ($q$), the ratio of the minor to the major axis, of the cluster galaxies from near-infrared Hubble Space Telescope imaging using Sérsic profile fitting and compare them with a field sample. We find that the median $q$ of both star-forming and quiescent… ▽ More We present results on the measured shapes of 832 galaxies in 11 galaxy clusters at 1.0 < z <1.4 from the GOGREEN survey. We measure the axis ratio ($q$), the ratio of the minor to the major axis, of the cluster galaxies from near-infrared Hubble Space Telescope imaging using Sérsic profile fitting and compare them with a field sample. We find that the median $q$ of both star-forming and quiescent galaxies in clusters increases with stellar mass, similar to the field. Comparing the axis ratio distributions between clusters and the field in four mass bins, the distributions for star-forming galaxies in clusters are consistent with those in the field. Conversely, the distributions for quiescent galaxies in the two environments are distinct, most remarkably in $10.1\leq\log(M/{\rm M}_{\odot})<10.5$ where clusters show a flatter distribution, with an excess at low $q$. Modelling the distribution with oblate and triaxial components, we find that the cluster and field sample difference is consistent with an excess of flattened oblate quiescent galaxies in clusters. The oblate population contribution drops at high masses, resulting in a narrower $q$ distribution in the massive population than at lower masses. Using a simple accretion model, we show that the observed $q$ distributions and quenched fractions are consistent with a scenario where no morphological transformation occurs for the environmentally quenched population in the two intermediate mass bins. Our results suggest that environmental quenching mechanism(s) likely produce a population that has a different morphological mix than those resulting from the dominant quenching mechanism in the field. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: Accepted for publication in ApJ. 25 pages, 15 figures

The Relativistic Jet Orientation and Host Galaxy of the Peculiar Blazar PKS 1413+135

Authors: A. C. S. Readhead, V. Ravi, I. Liodakis, M. L. Lister, V. Singh, M. F. Aller, R. D. Blandford, I. W. A. Browne, V. Gorjian, K. J. B. Grainge, M. A. Gurwell, M. W. Hodges, T. Hovatta, S. Kiehlmann, A. Lähteenmäki, T. McAloone, W. Max-Moerbeck, V. Pavlidou, T. J. Pearson, A. L. Peirson, E. S. Perlman, R. A. Reeves, B. T. Soifer, G. B. Taylor, M. Tornikoski , et al. (4 additional authors not shown)

Abstract: PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and i… ▽ More PKS 1413+135 is one of the most peculiar blazars known. Its strange properties led to the hypothesis almost four decades ago that it is gravitationally lensed by a mass concentration associated with an intervening galaxy. It exhibits symmetric achromatic variability, a rare form of variability that has been attributed to gravitational milli-lensing. It has been classified as a BL Lac object, and is one of the rare objects in this class with a visible counterjet. BL Lac objects have jet axes aligned close to the line of sight. It has also been classified as a compact symmetric object, which have jet axes not aligned close to the line of sight. Intensive efforts to understand this blazar have hitherto failed to resolve even the questions of the orientation of the relativistic jet, and the host galaxy. Answering these two questions is important as they challenge our understanding of jets in active galactic nuclei and the classification schemes we use to describe them. We show that the jet axis is aligned close to the line of sight and PKS 1413+135 is almost certainly not located in the apparent host galaxy, but is a background object in the redshift range $0.247 < z < 0.5$. The intervening spiral galaxy at $z = 0.247$ provides a natural host for the putative lens responsible for symmetric achromatic variability and is shown to be a Seyfert 2 galaxy. We also show that, as for the radio emission, a "multizone" model is needed to account for the high-energy emission. △ Less

Submitted 7 December, 2020; originally announced December 2020.

Comments: Accepted for publication in The Astrophysical Journal. 28 pages, 15 figures

Filament rotation in the California L1482 cloud

Authors: Rodrigo H. Álvarez-Gutiérrez, Amelia M. Stutz, Chi Y. Law, Stefan Reissl, Ralf S. Klessen, Nathan W. C. Leigh, Hong-Li Liu, Rodrigo A. Reeves

Abstract: We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C… ▽ More We analyze the gas mass distribution, the gas kinematics, and the young stellar object (YSO) content of the California Molecular Cloud (CMC) L1482 filament. We derive a Gaia DR2 YSO distance of 511$^{+17}_{-16}$ pc. We derive scale-free power-laws for the mean gas line-mass (M/L) profiles; we calculate the gravitational potential and field profiles consistent with these. We present IRAM 30 m C$^{18}$O (1-0) (and other tracers) position-velocity (PV) diagrams that exhibit complex velocity twisting and turning structures. We find a rotational profile in C$^{18}$O perpendicular to the southern filament ridgeline. The profile is regular, confined ($r\lesssim0.4$ pc), anti-symmetric, and to first order linear with a break at $r\sim0.25$ pc. The timescales of the inner (outer) gradients are $\sim$0.7 (6.0) Myr. We show that the centripetal force, compared to gravity, increases toward the break; when the ratio of forces approaches unity, the profile turns over, just before filament breakup is achieved. The timescales and relative roles of gravity to rotation indicate that the structure is stable, long lived ($\sim$ a few times 6 Myr), and undergoing outside-in evolution. Moreover, this filament has practically no star formation, a perpendicular Planck plane-of-the-sky (POS) magnetic field morphology, and POS "zig-zag" morphology, which together with the rotation profile lead to the suggestion that the 3D shape is a corkscrew filament with a helical magnetic field. These results, combined with results in Orion and G035.39-00.33, suggest evolution toward higher densities as rotating filaments shed angular momentum. Thus, magnetic fields may be an essential feature of high-mass (M $\sim10^5$ M$_{\odot}$) cloud filament evolution toward cluster formation. △ Less

Submitted 17 December, 2020; v1 submitted 21 October, 2020; originally announced October 2020.

Comments: 23 pages, 17 figures. Accepted for publication in the Astrophysical Journal

The GOGREEN and GCLASS Surveys: First Data Release

Authors: Michael L. Balogh, Remco F. J. van der Burg, Adam Muzzin, Gregory Rudnick, Gillian Wilson, Kristi Webb, Andrea Biviano, Kevin Boak, Pierluigi Cerulo, Jeffrey Chan, M. C. Cooper, David G. Gilbank, Stephen Gwyn, Chris Lidman, Jasleen Matharu, Sean L. McGee, Lyndsay Old, Irene Pintos-Castro, Andrew M. M. Reeves, Heath Shipley, Benedetta Vulcani, Howard K. C. Yee, M. Victoria Alonso, Callum Bellhouse, Kevin C. Cooke , et al. (20 additional authors not shown)

Abstract: We present the first public data release of the GOGREEN and GCLASS surveys of galaxies in dense environments, spanning a redshift range $0.8<z<1.5$. The surveys consist of deep, multiwavelength photometry and extensive Gemini GMOS spectroscopy of galaxies in 26 overdense systems ranging in halo mass from small groups to the most massive clusters. The objective of both projects was primarily to und… ▽ More We present the first public data release of the GOGREEN and GCLASS surveys of galaxies in dense environments, spanning a redshift range $0.8<z<1.5$. The surveys consist of deep, multiwavelength photometry and extensive Gemini GMOS spectroscopy of galaxies in 26 overdense systems ranging in halo mass from small groups to the most massive clusters. The objective of both projects was primarily to understand how the evolution of galaxies is affected by their environment, and to determine the physical processes that lead to the quenching of star formation. There was an emphasis on obtaining unbiased spectroscopy over a wide stellar mass range ($M\gtrsim 2\times 10^{10}~\mathrm{M}_\odot$), throughout and beyond the cluster virialized regions. The final spectroscopic sample includes 2771 unique objects, of which 2257 have reliable spectroscopic redshifts. Of these, 1704 have redshifts in the range $0.8<z<1.5$, and nearly 800 are confirmed cluster members. Imaging spans the full optical and near-infrared wavelength range, at depths comparable to the UltraVISTA survey, and includes \textit{HST}/WFC3 F160W (GOGREEN) and F140W (GCLASS). This data release includes fully reduced images and spectra, with catalogues of advanced data products including redshifts, line strengths, star formation rates, stellar masses and rest-frame colours. Here we present an overview of the data, including an analysis of the spectroscopic completeness and redshift quality. △ Less

Submitted 28 September, 2020; originally announced September 2020.

Comments: Accepted for publication in MNRAS Sept 28, 2020

Association of IceCube neutrinos with radio sources observed at Owens Valley and Metsähovi Radio Observatories

Authors: T. Hovatta, E. Lindfors, S. Kiehlmann, W. Max-Moerbeck, M. Hodges, I. Liodakis, A. Lahteenmaki, T. J. Pearson, A. C. S. Readhead, R. A. Reeves, S. Suutarinen, J. Tammi, M. Tornikoski

Abstract: Identifying the most likely sources for high-energy neutrino emission has been one of the main topics in high-energy astrophysics ever since the first observation of high-energy neutrinos by the IceCube Neutrino Observatory. Active galactic nuclei with relativistic jets, also known as blazars, have been considered to be one of the main candidates because of their ability to accelerate particles to… ▽ More Identifying the most likely sources for high-energy neutrino emission has been one of the main topics in high-energy astrophysics ever since the first observation of high-energy neutrinos by the IceCube Neutrino Observatory. Active galactic nuclei with relativistic jets, also known as blazars, have been considered to be one of the main candidates because of their ability to accelerate particles to high energies. We study the connection between radio emission and IceCube neutrino events using data from the Owens Valley Radio Observatory (OVRO) and Metsähovi Radio Observatory blazar monitoring programs. We identify sources in our radio monitoring sample that are positionally consistent with IceCube high-energy neutrino events. We estimate their mean flux density and variability amplitudes around the neutrino arrival time, and compare these with values from random samples to establish the significance of our results. We find radio source associations within our samples with 15 high-energy neutrino events detected by IceCube. Nearly half of the associated sources are not detected in the $γ$-ray energies, but their radio variability properties and Doppler boosting factors are similar to the $γ$-ray detected objects in our sample, meaning that they could still be potential neutrino emitters. We find that the number of strongly flaring objects in our statistically complete OVRO samples is unlikely to be a random coincidence (at $2σ$ level). Based on our results, we conclude that although it is clear that not all neutrino events are associated with strong radio flaring blazars, observations of large-amplitude radio flares in a blazar at the same time as a neutrino event are unlikely to be a random coincidence. △ Less

Submitted 15 April, 2021; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: Accepted for publication in A&A (18 pages, 6 figures). Accepted version (v2) includes changes to methodology and associations, but the main conclusions are unaffected. V3 minor fix to numbers in Sect. 3.3.1. and some language edits

Journal ref: A&A 650, A83 (2021)

The GOGREEN survey: Post-infall environmental quenching fails to predict the observed age difference between quiescent field and cluster galaxies at z>1

Authors: Kristi Webb, Michael L. Balogh, Joel Leja, Remco F. J. van der Burg, Gregory Rudnick, Adam Muzzin, Kevin Boak, Pierluigi Cerulo, David Gilbank, Chris Lidman, Lyndsay J. Old, Irene Pintos-Castro, Sean McGee, Heath Shipley, Andrea Biviano, Jeffrey C. C. Chan, Michael Cooper, Gabriella De Lucia, Ricardo Demarco, Ben Forrest, Pascale Jablonka, Egidijus Kukstas, Ian G. McCarthy, Karen McNab, Julie Nantais , et al. (7 additional authors not shown)

Abstract: We study the star formation histories (SFHs) and mass-weighted ages of 331 UVJ-selected quiescent galaxies in 11 galaxy clusters and in the field at 1<z<1.5 from the Gemini Observations of Galaxies in Rich Early ENvironments (GOGREEN) survey. We determine the SFHs of individual galaxies by simultaneously fitting rest-frame optical spectroscopy and broadband photometry to stellar population models.… ▽ More We study the star formation histories (SFHs) and mass-weighted ages of 331 UVJ-selected quiescent galaxies in 11 galaxy clusters and in the field at 1<z<1.5 from the Gemini Observations of Galaxies in Rich Early ENvironments (GOGREEN) survey. We determine the SFHs of individual galaxies by simultaneously fitting rest-frame optical spectroscopy and broadband photometry to stellar population models. We confirm that the SFHs are consistent with more massive galaxies having on average earlier formation times. Comparing galaxies found in massive clusters with those in the field, we find galaxies with $M_\ast<10^{11.3}$ M$_{\odot}$ in the field have more extended SFHs. From the SFHs we calculate the mass-weighted ages, and compare age distributions of galaxies between the two environments, at fixed mass. We constrain the difference in mass-weighted ages between field and cluster galaxies to $0.31_{^{-0.33}}^{_{+0.51}}$ Gyr, in the sense that cluster galaxies are older. We place this result in the context of two simple quenching models and show that neither environmental quenching based on time since infall (without pre-processing) nor a difference in formation times alone can reproduce both the average age difference and relative quenched fractions. This is distinctly different from local clusters, for which the majority of the quenched population is consistent with having been environmentally quenched upon infall. Our results suggest that quenched population in galaxy clusters at z>1 has been driven by different physical processes than those at play at z=0. △ Less

Submitted 8 September, 2020; originally announced September 2020.

Comments: accepted Sept 7 2020, MNRAS

The GOGREEN Survey: A deep stellar mass function of cluster galaxies at 1.0<z<1.4 and the complex nature of satellite quenching

Authors: Remco F. J. van der Burg, Gregory Rudnick, Michael L. Balogh, Adam Muzzin, Chris Lidman, Lyndsay J. Old, Heath Shipley, David Gilbank, Sean McGee, Andrea Biviano, Pierluigi Cerulo, Jeffrey C. C. Chan, Michael Cooper, Gabriella De Lucia, Ricardo Demarco, Ben Forrest, Stephen Gwyn, Pascale Jablonka, Egidijus Kukstas, Danilo Marchesini, Julie Nantais, Allison Noble, Irene Pintos-Castro, Bianca Poggianti, Andrew M. M. Reeves , et al. (6 additional authors not shown)

Abstract: We study the stellar mass functions (SMFs) of star-forming and quiescent galaxies in 11 galaxy clusters at 1.0<z<1.4, drawn from the Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) survey. Based on more than 500 hours of Gemini/GMOS spectroscopy, and deep multi-band photometry taken with a range of observatories, we probe the SMFs down to a stellar mass limit of 10^9.7 Msun (1… ▽ More We study the stellar mass functions (SMFs) of star-forming and quiescent galaxies in 11 galaxy clusters at 1.0<z<1.4, drawn from the Gemini Observations of Galaxies in Rich Early Environments (GOGREEN) survey. Based on more than 500 hours of Gemini/GMOS spectroscopy, and deep multi-band photometry taken with a range of observatories, we probe the SMFs down to a stellar mass limit of 10^9.7 Msun (10^9.5 Msun for star-forming galaxies). At this early epoch, the fraction of quiescent galaxies is already highly elevated in the clusters compared to the field at the same redshift. The quenched fraction excess (QFE) represents the fraction of galaxies that would be star-forming in the field, but are quenched due to their environment. The QFE is strongly mass dependent, and increases from ~30% at Mstar=10^9.7 Msun, to ~80% at Mstar=10^11.0 Msun. Nonetheless, the shapes of the SMFs of the two individual galaxy types, star-forming and quiescent galaxies, are identical between the clusters and the field - to high statistical precision. Yet, along with the different quiescent fractions is the total galaxy SMF environmentally dependent, with a relative deficit of low-mass galaxies in the clusters. These results are in stark contrast with findings in the local Universe, and thus require a substantially different quenching mode to operate at early times. We discuss these results in the light of several popular quenching models. △ Less

Submitted 22 April, 2020; originally announced April 2020.

Comments: 19 pages, 9 figures (excluding appendices). Accepted for publication in A&A

Journal ref: A&A 638, A112 (2020)

The GOGREEN survey: The environmental dependence of the star-forming galaxy main sequence at $1.0<z<1.5$

Authors: Lyndsay J. Old, Michael. L. Balogh, Remco F. J van der Burg, Andrea Biviano, Howard K. C. Yee, Irene Pintos-Castro, Kristi Webb, Adam Muzzin, Gregory Rudnick, Benedetta Vulcani, Bianca Poggianti, Michael Cooper, Dennis Zaritsky, Pierluigi Cerulo, Gillian Wilson, Jeffrey C. C. Chan, Chris Lidman, Sean McGee, Ricardo Demarco, Ben Forrest, Gabriella De Lucia, David Gilbank, Egidijus Kukstas, Ian G. McCarthy, Pascale Jablonka , et al. (4 additional authors not shown)

Abstract: We present results on the environmental dependence of the star-forming galaxy main sequence in 11 galaxy cluster fields at $1.0 < z < 1.5$ from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN) survey. We use a homogeneously selected sample of field and cluster galaxies whose membership is derived from dynamical analysis. Using [OII]-derived star formation rates (SFRs… ▽ More We present results on the environmental dependence of the star-forming galaxy main sequence in 11 galaxy cluster fields at $1.0 < z < 1.5$ from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN) survey. We use a homogeneously selected sample of field and cluster galaxies whose membership is derived from dynamical analysis. Using [OII]-derived star formation rates (SFRs), we find that cluster galaxies have suppressed SFRs at fixed stellar mass in comparison to their field counterparts by a factor of 1.4 $\pm$ 0.1 ($\sim3.3σ$) across the stellar mass range: $9.0 < \log(M_{*} /M_{\odot}) < 11.2$. We also find that this modest suppression in the cluster galaxy star-forming main sequence is mass and redshift dependent: the difference between cluster and field increases towards lower stellar masses and lower redshift. When comparing the distribution of cluster and field galaxy SFRs to the star-forming main sequence, we find an overall shift towards lower SFRs in the cluster population, and note the absence of a tail of high SFR galaxies as seen in the field. Given this observed suppression in the cluster galaxy star-forming main sequence, we explore the implications for several scenarios such as formation time differences between cluster and field galaxies, and environmentally-induced star formation quenching and associated timescales. △ Less

Submitted 28 September, 2020; v1 submitted 26 February, 2020; originally announced February 2020.

Comments: Accepted for publication in MNRAS on February 25 2020, appended erratum accepted in MNRAS on September 17 2020

AOtools -- a Python package for adaptive optics modelling and analysis

Authors: M. J. Townson, O. J. D. Farley, G. Orban de Xivry, J. Osborn, A. P. Reeves

Abstract: AOtools is a Python package which is open-source and aimed at providing tools for adaptive optics users and researchers. We present version 1.0 which contains tools for adaptive optics processing, including analysing data in the pupil plane, images and point spread functions in the focal plane, wavefront sensors, modelling of atmospheric turbulence, physical optical propagation of wavefronts, and… ▽ More AOtools is a Python package which is open-source and aimed at providing tools for adaptive optics users and researchers. We present version 1.0 which contains tools for adaptive optics processing, including analysing data in the pupil plane, images and point spread functions in the focal plane, wavefront sensors, modelling of atmospheric turbulence, physical optical propagation of wavefronts, and conversion between frequently used adaptive optics and astronomical units. The main drivers behind AOtools is that it should be easy to install and use. To achieve this the project features extensive documentation, automated unit testing and is registered on the Python Package Index. AOtools is under continuous active development to expand the features available and we encourage everyone involved in adaptive optics to become involved and contribute to the project. △ Less

Submitted 10 October, 2019; originally announced October 2019.

Comments: Accepted in Optics Express

The Rest-frame $H$-band Luminosity Function of Red Sequence Galaxies in Clusters at $1.0 < z < 1.3$

Authors: Jeffrey C. C. Chan, Gillian Wilson, Gregory Rudnick, Adam Muzzin, Michael Balogh, Julie Nantais, Remco F. J. van der Burg, Pierluigi Cerulo, Andrea Biviano, Michael C. Cooper, Ricardo Demarco, Ben Forrest, Chris Lidman, Allison Noble, Lyndsay Old, Irene Pintos-Castro, Andrew M. M. Reeves, Kristi A. Webb, Howard K. C. Yee, Mohamed H. Abdullah, Gabriella De Lucia, Danilo Marchesini, Sean L. McGee, Mauro Stefanon, Dennis Zaritsky

Abstract: We present results on the rest-frame $H$-band luminosity functions (LF) of red sequence galaxies in seven clusters at 1.0 < z < 1.3 from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN). Using deep GMOS-z' and IRAC $3.6 μ$m imaging, we identify red sequence galaxies and measure their LFs down to $M_{H} \sim M_{H}^{*} + (2.0 - 3.0)$. By stacking the entire sample, we… ▽ More We present results on the rest-frame $H$-band luminosity functions (LF) of red sequence galaxies in seven clusters at 1.0 < z < 1.3 from the Gemini Observations of Galaxies in Rich Early Environments Survey (GOGREEN). Using deep GMOS-z' and IRAC $3.6 μ$m imaging, we identify red sequence galaxies and measure their LFs down to $M_{H} \sim M_{H}^{*} + (2.0 - 3.0)$. By stacking the entire sample, we derive a shallow faint end slope of $ α\sim -0.35^{+0.15}_{-0.15} $ and $ M_{H}^{*} \sim -23.52^{+0.15}_{-0.17} $, suggesting that there is a deficit of faint red sequence galaxies in clusters at high redshift. By comparing the stacked red sequence LF of our sample with a sample of clusters at z~0.6, we find an evolution in the faint end of the red sequence over the ~2.6 Gyr between the two samples, with the mean faint end red sequence luminosity growing by more than a factor of two. The faint-to-luminous ratio of our sample ($0.78^{+0.19}_{-0.15}$) is consistent with the trend of decreasing ratio with increasing redshift as proposed in previous studies. A comparison with the field shows that the faint-to-luminous ratios in clusters are consistent with the field at z~1.15 and exhibit a stronger redshift dependence. Our results support the picture that the build up of the faint red sequence galaxies occurs gradually over time and suggest that faint cluster galaxies, similar to bright cluster galaxies, experience the quenching effect induced by environment already at z~1.15. △ Less

Submitted 11 July, 2019; v1 submitted 25 June, 2019; originally announced June 2019.

Comments: Accepted for publication in ApJ. 24 pages, 13 figures

Optimizing the accuracy and efficiency of optical turbulence profiling using adaptive optics telemetry for extremely large telescopes

Authors: Douglas J Laidlaw, James Osborn, Timothy J Morris, Alastair G Basden, Olivier Beltramo-Martin, Timothy Butterley, Eric Gendron, Andrew P Reeves, Gérard Rousset, Matthew J Townson, Richard W Wilson

Abstract: Advanced adaptive optics (AO) instruments on ground-based telescopes require accurate knowledge of the atmospheric turbulence strength as a function of altitude. This information assists point spread function reconstruction, AO temporal control techniques and is required by wide-field AO systems to optimize the reconstruction of an observed wavefront. The variability of the atmosphere makes it imp… ▽ More Advanced adaptive optics (AO) instruments on ground-based telescopes require accurate knowledge of the atmospheric turbulence strength as a function of altitude. This information assists point spread function reconstruction, AO temporal control techniques and is required by wide-field AO systems to optimize the reconstruction of an observed wavefront. The variability of the atmosphere makes it important to have a measure of the optical turbulence profile in real time. This measurement can be performed by fitting an analytically generated covariance matrix to the cross-covariance of Shack-Hartmann wavefront sensor (SHWFS) centroids. In this study we explore the benefits of reducing cross-covariance data points to a covariance map region of interest (ROI). A technique for using the covariance map ROI to measure and compensate for SHWFS misalignments is also introduced. We compare the accuracy of covariance matrix and map ROI optical turbulence profiling using both simulated and on-sky data from CANARY, an AO demonstrator on the 4.2 m William Herschel telescope, La Palma. On-sky CANARY results are compared to contemporaneous profiles from Stereo-SCIDAR - a dedicated high-resolution optical turbulence profiler. It is shown that the covariance map ROI optimizes the accuracy of AO telemetry optical turbulence profiling. In addition, we show that the covariance map ROI reduces the fitting time for an extremely large telescope-scale system by a factor of 72. The software package we developed to collect all of the presented results is now open source. △ Less

Submitted 17 January, 2019; originally announced January 2019.

Journal ref: Published in MNRAS, 483, 4, 4341-4353, 2018

Optimizing astrophotonic spatial reformatters using simulated on-sky performance

Authors: Theodoros Anagnos, Robert J. Harris, Mark K. Corrigan, Andrew P. Reeves, Matthew J. Townson, David G. MacLachlan, Robert R. Thomson, Tim J. Morris, Christian Schwab, Andreas Quirrenbach

Abstract: One of the most useful techniques in astronomical instrumentation is image slicing. It enables a spectrograph to have a more compact angular slit, whilst retaining throughput and increasing resolving power. Astrophotonic components like the photonic lanterns and photonic reformatters can be used to replace bulk optics used so far. This study investigates the performance of such devices using end-t… ▽ More One of the most useful techniques in astronomical instrumentation is image slicing. It enables a spectrograph to have a more compact angular slit, whilst retaining throughput and increasing resolving power. Astrophotonic components like the photonic lanterns and photonic reformatters can be used to replace bulk optics used so far. This study investigates the performance of such devices using end-to-end simulations to approximate realistic on-sky conditions. It investigates existing components, tries to optimize their performance and aims to understand better how best to design instruments to maximize their performance. This work complements the recent work in the field and provides an estimation for the performance of the new components. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: Conference proceedings in SPIE 2018 Austin Texas

Multi-wavelength characterization of the blazar S5~0716+714 during an unprecedented outburst phase

Authors: MAGIC Collaboration, M. L. Ahnen, S. Ansoldi, L. A. Antonelli, C. Arcaro, D. Baack, A. Babić, B. Banerjee, P. Bangale, U. Barres de Almeida, J. A. Barrio, J. Becerra González, W. Bednarek, E. Bernardini, R. Ch. Berse, A. Berti, W. Bhattacharyya, A. Biland, O. Blanch, G. Bonnoli, R. Carosi, A. Carosi, G. Ceribella, A. Chatterjee, S. M. Colak , et al. (165 additional authors not shown)

Abstract: The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light… ▽ More The BL Lac object S5~0716+714, a highly variable blazar, underwent an impressive outburst in January 2015 (Phase A), followed by minor activity in February (Phase B). The MAGIC observations were triggered by the optical flux observed in Phase A, corresponding to the brightest ever reported state of the source in the R-band. The comprehensive dataset collected is investigated in order to shed light on the mechanism of the broadband emission. Multi-wavelength light curves have been studied together with the broadband Spectral Energy Distributions (SEDs). The data set collected spans from radio, optical photometry and polarimetry, X-ray, high-energy (HE, 0.1 GeV < E < 100 GeV) with \textit{Fermi}-LAT to the very-high-energy (VHE, E>100 GeV) with MAGIC. The flaring state of Phase A was detected in all the energy bands, providing for the first time a multi-wavelength sample of simultaneous data from the radio band to the VHE. In the constructed SED the \textit{Swift}-XRT+\textit{NuSTAR} data constrain the transition between the synchrotron and inverse Compton components very accurately, while the second peak is constrained from 0.1~GeV to 600~GeV by \textit{Fermi}+MAGIC data. The broadband SED cannot be described with a one-zone synchrotron self-Compton model as it severely underestimates the optical flux in order to reproduce the X-ray to $γ$-ray data. Instead we use a two-zone model. The EVPA shows an unprecedented fast rotation. An estimation of the redshift of the source by combined HE and VHE data provides a value of $z = 0.31 \pm 0.02_{stats} \pm 0.05_{sys}$, confirming the literature value. The data show the VHE emission originating in the entrance and exit of a superluminal knot in and out a recollimation shock in the inner jet. A shock-shock interaction in the jet seems responsible for the observed flares and EVPA swing. This scenario is also consistent with the SED modelling. △ Less

Submitted 1 July, 2018; originally announced July 2018.

Journal ref: Astronomy & Astrophysics 619, A45 (2018)

Simulation and Optimization of an Astrophotonic Reformatter

Authors: Th. Anagnos, R. J. Harris, M. K. Corrigan, A. P. Reeves, M. J. Townson, D. G. MacLachlan, R. R. Thomson, T. J. Morris, C. Schwab, A. Quirrenbach

Abstract: Image slicing is a powerful technique in astronomy. It allows the instrument designer to reduce the slit width of the spectrograph, increasing spectral resolving power whilst retaining throughput. Conventionally this is done using bulk optics, such as mirrors and prisms, however more recently astrophotonic components known as photonic lanterns (PLs) and photonic reformatters have also been used. T… ▽ More Image slicing is a powerful technique in astronomy. It allows the instrument designer to reduce the slit width of the spectrograph, increasing spectral resolving power whilst retaining throughput. Conventionally this is done using bulk optics, such as mirrors and prisms, however more recently astrophotonic components known as photonic lanterns (PLs) and photonic reformatters have also been used. These devices reformat the multi-mode (MM) input light from a telescope into single-mode (SM) outputs, which can then be re-arranged to suit the spectrograph. The photonic dicer (PD) is one such device, designed to reduce the dependence of spectrograph size on telescope aperture and eliminate modal noise. We simulate the PD, by optimising the throughput and geometrical design using Soapy and BeamProp. The simulated device shows a transmission between 8 and 20 %, depending upon the type of adaptive optics (AO) correction applied, matching the experimental results well. We also investigate our idealised model of the PD and show that the barycentre of the slit varies only slightly with time, meaning that the modal noise contribution is very low when compared to conventional fibre systems. We further optimise our model device for both higher throughput and reduced modal noise. This device improves throughput by 6.4 % and reduces the movement of the slit output by 50%, further improving stability. This shows the importance of properly simulating such devices, including atmospheric effects. Our work complements recent work in the field and is essential for optimising future photonic reformatters. △ Less

Submitted 2 May, 2018; originally announced May 2018.

37 GHz observations of narrow-line Seyfert 1 galaxies

Authors: A. Lähteenmäki, E. Järvelä, T. Hovatta, M. Tornikoski, D. L. Harrison, M. López-Caniego, W. Max-Moerbeck, M. Mingaliev, T. J. Pearson, V. Ramakrishnan, A. C. S. Readhead, R. A. Reeves, J. L. Richards, Y. Sotnikova, J. Tammi

Abstract: Observations at 37 GHz, performed at Metsähovi Radio Observatory, are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even… ▽ More Observations at 37 GHz, performed at Metsähovi Radio Observatory, are presented for a sample of 78 radio-loud and radio-quiet narrow-line Seyfert 1 (NLS1) galaxies, together with additional lower and higher frequency radio data from RATAN-600, Owens Valley Radio Observatory, and the Planck satellite. Most of the data have been gathered between February 2012 and April 2015 but for some sources even longer lightcurves exist. The detection rate at 37 GHz is around 19%, comparable to other populations of active galactic nuclei presumed to be faint at radio frequencies, such as BL Lac objects. Variability and spectral indices are determined for sources with enough detections. Based on the radio data, many NLS1 galaxies show a blazar-like radio spectra exhibiting significant variability. The spectra at a given time are often inverted or convex. The source of the high-frequency radio emission in NLS1 galaxies, detected at 37 GHz, is most probably a relativistic jet rather than star formation. Jets in NLS1 galaxies are therefore expected to be a much more common phenomenon than earlier assumed. △ Less

Submitted 30 March, 2017; originally announced March 2017.

Comments: Accepted for publication in A&A. Table of 37 GHz data will be available at the CDS soon

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

Symmetric Achromatic Variability in Active Galaxies -- A Powerful New Gravitational Lensing Probe?

Authors: H. K. Vedantham, A. C. S. Readhead, T. Hovatta, T. J. Pearson, R. D. Blandford, M. A. Gurwell, A. Lähteenmäki, W. Max-Moerbeck, V. Pavlidou, V. Ravi, R. A. Reeves, J. L. Richards, M. Tornikoski, J. A. Zensus

Abstract: We report the discovery of a rare new form of long-term radio variability in the light-curves of active galaxies (AG) --- Symmetric Achromatic Variability (SAV) --- a pair of opposed and strongly skewed peaks in the radio flux density observed over a broad frequency range. We propose that SAV arises through gravitational milli-lensing when relativistically moving features in AG jets move through g… ▽ More We report the discovery of a rare new form of long-term radio variability in the light-curves of active galaxies (AG) --- Symmetric Achromatic Variability (SAV) --- a pair of opposed and strongly skewed peaks in the radio flux density observed over a broad frequency range. We propose that SAV arises through gravitational milli-lensing when relativistically moving features in AG jets move through gravitational lensing caustics created by $10^3-10^6 \;{\rm M}_{\odot}$ subhalo condensates or black holes located within intervening galaxies. The lower end of this mass range has been inaccessible with previous gravitational lensing techniques. This new interpretation of some AG variability can easily be tested and if it passes these tests, will enable a new and powerful probe of cosmological matter distribution on these intermediate mass scales, as well as provide, for the first time, micro-arcsecond resolution of the nuclei of AG --- a factor of 30--100 greater resolution than is possible with ground-based millimeter VLBI. △ Less

Submitted 17 February, 2017; originally announced February 2017.

Comments: Submitted version

The peculiar light-curve of J1415+1320: A case study in extreme scattering events

Authors: H. K. Vedantham, A. C. S. Readhead, T. Hovatta, L. V. E. Koopmans, T. J. Pearson, R. D. Blandford, M. A. Gurwell, A. Lähteenmäki, W. Max-Moerbeck, V. Pavlidou, V. Ravi, R. A. Reeves, J. L. Richards, M. Tornikoski, J. A. Zensus

Abstract: The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing… ▽ More The radio light-curve of J1415+1320 (PKS 1413+135) shows time-symmetric and recurring U-shaped features across the cm-wave and mm-wave bands. The symmetry of these features points to lensing by an intervening object as the cause. U-shaped events in radio light curves in the cm-wave band have previously been attributed to Extreme scattering events (ESE). ESEs are thought to be the result of lensing by compact plasma structures in the Galactic interstellar medium, but the precise nature of these plasma structures remains unknown. Since the strength of a plasma lens evolves with wavelength $λ$ as $λ^2$, the presence of correlated variations at over a wide wavelength range casts doubt on the canonical ESE interpretation for J1415+1320. In this paper, we critically examine the evidence for plasma lensing in J1415+1320. We compute limits on the lensing strength, and the associated free-free opacity of the putative plasma lenses. We compare the observed and model ESE light curves, and also derive a lower limit on the lens distance based on the effects of parallax due to the Earth's orbit around the Sun. We conclude that plasma lensing is not a viable interpretation for J1415+1320's light curves and that symmetric U-shaped features in the radio light curves of extragalactic sources do not present {\em prima facie} evidence for ESEs. The methodology presented here is generic enough to be applicable to any plasma lensing candidate. △ Less

Submitted 17 February, 2017; originally announced February 2017.

Comments: Submitted version

On-sky demonstration of matched filters for wavefront measurements using ELT-scale elongated laser guide stars

Authors: A. G. Basden, L. Bardou, D. Bonaccini Calia, T. Buey, M. Centrone, F. Chemla, J. L. Gach, E. Gendron, D. Gratadour, I. Guidolin, D. R. Jenkins, E. Marchetti, T. J. Morris, R. M. Myers, J. Osborn, A. P. Reeves, M. Reyes, G. Rousset, G. Lombardi, M. J. Townson, F. Vidal

Abstract: The performance of adaptive optics systems is partially dependant on the algorithms used within the real-time control system to compute wavefront slope measurements. We demonstrate use of a matched filter algorithm for the processing of elongated laser guide star (LGS) Shack-Hartmann images, using the CANARY adaptive optics instrument on the 4.2m William Herschel Telescope and the European Souther… ▽ More The performance of adaptive optics systems is partially dependant on the algorithms used within the real-time control system to compute wavefront slope measurements. We demonstrate use of a matched filter algorithm for the processing of elongated laser guide star (LGS) Shack-Hartmann images, using the CANARY adaptive optics instrument on the 4.2m William Herschel Telescope and the European Southern Observatory Wendelstein LGS Unit placed 40m away. This algorithm has been selected for use with the forthcoming Thirty Meter Telescope, but until now had not been demonstrated on-sky. From the results of a first observing run, we show that the use of matched filtering improves our adaptive optics system performance, with increases in on-sky H-band Strehl measured up to about a factor of 1.1 with respect to a conventional centre of gravity approach. We describe the algorithm used, and the methods that we implemented to enable on-sky demonstration. △ Less

Submitted 13 February, 2017; originally announced February 2017.

Comments: Accepted for publication in MNRAS

Turbulence velocity profiling for high sensitivity and vertical-resolution atmospheric characterization with Stereo-SCIDAR

Authors: J. Osborn, T. Butterley, M. J. Townson, A. P. Reeves, T. J. Morris, R. W. Wilson

Abstract: As telescopes become larger, into the era of ~40 m Extremely Large Telescopes, the high- resolution vertical profile of the optical turbulence strength is critical for the validation, optimization and operation of optical systems. The velocity of atmospheric optical turbulence is an important parameter for several applications including astronomical adaptive optics systems. Here, we compare the ve… ▽ More As telescopes become larger, into the era of ~40 m Extremely Large Telescopes, the high- resolution vertical profile of the optical turbulence strength is critical for the validation, optimization and operation of optical systems. The velocity of atmospheric optical turbulence is an important parameter for several applications including astronomical adaptive optics systems. Here, we compare the vertical profile of the velocity of the atmospheric wind above La Palma by means of a comparison of Stereo-SCIntillation Detection And Ranging (Stereo- SCIDAR) with the Global Forecast System models and nearby balloon-borne radiosondes. We use these data to validate the automated optical turbulence velocity identification from the Stereo-SCIDAR instrument mounted on the 2.5 m Isaac Newton Telescope, La Palma. By comparing these data we infer that the turbulence velocity and the wind velocity are consistent and that the automated turbulence velocity identification of the Stereo-SCIDAR is precise. The turbulence velocities can be used to increase the sensitivity of the turbulence strength profiles, as weaker turbulence that may be misinterpreted as noise can be detected with a velocity vector. The turbulence velocities can also be used to increase the altitude resolution of a detected layer, as the altitude of the velocity vectors can be identified to a greater precision than the native resolution of the system. We also show examples of complex velocity structure within a turbulent layer caused by wind shear at the interface of atmospheric zones. △ Less

Submitted 10 January, 2017; originally announced January 2017.

Comments: Published in MNRAS, 464, 3998-4007, 2017

Journal ref: MNRAS 464, 3998-4007, 2017

The awakening of the gamma-ray narrow-Line Seyfert 1 galaxy PKS 1502+036

Authors: F. D'Ammando, M. Orienti, J. Finke, T. Hovatta, M. Giroletti, W. Max-Moerbeck, T. J. Pearson, A. C. S. Readhead, R. A. Reeves, J. L. Richards

Abstract: After a long low-activity period, a gamma-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z=0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20 the source reached a daily peak flux, in the 0.1-300 GeV band, of (93 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-h… ▽ More After a long low-activity period, a gamma-ray flare from the narrow-line Seyfert 1 PKS 1502+036 (z=0.4089) was detected by the Large Area Telescope (LAT) on board Fermi in 2015. On 2015 December 20 the source reached a daily peak flux, in the 0.1-300 GeV band, of (93 $\pm$ 19) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$, attaining a flux of (237 $\pm$ 71) $\times$10$^{-8}$ ph cm$^{-2}$ s$^{-1}$ on 3-hr time-scales, which corresponds to an isotropic luminosity of (7.3 $\pm$ 2.1) $\times$10$^{47}$ erg/s. The gamma-ray flare was not accompanied by significant spectral changes. We report on multi-wavelength radio-to-gamma-ray observations of PKS 1502+036 during 2008 August-2016 March by Fermi-LAT, Swift, XMM-Newton, Catalina Real-Time Transient Survey, and the Owens Valley Radio Observatory (OVRO). An increase in activity was observed on 2015 December 22 by Swift in optical, UV, and X-rays. The OVRO 15 GHz light curve reached the highest flux density observed from this source on 2016 January 12, indicating a delay of about three weeks between the gamma-ray and 15 GHz emission peaks. This suggests that the gamma-ray emitting region is located beyond the broad line region. We compared the spectral energy distribution (SED) of an average activity state with that of the flaring state. The two SED, with the high-energy bump modelled as an external Compton component with seed photons from a dust torus, could be fitted by changing the electron distribution parameters as well as the magnetic field. The fit of the disc emission during the average state constrains the black hole mass to values lower than 10$^8$ solar masses. The SED, high-energy emission mechanisms, and gamma-ray properties of the source resemble those of a flat spectrum radio quasar. △ Less

Submitted 15 September, 2016; originally announced September 2016.

Comments: 13 pages, 4 tables, 7 figures. Accepted for publication in the Monthly Notices of the Royal Astronomical Society. arXiv admin note: text overlap with arXiv:1305.0572

Experience with wavefront sensor and deformable mirror interfaces for wide-field adaptive optics systems

Authors: A. G. Basden, D. Atkinson, N. A. Bharmal, U. Bitenc, M. Brangier, T. Buey, T. Butterley, D. Cano, F. Chemla, P. Clark, M. Cohen, J. -M. Conan, F. J. de Cos, C. Dickson, N. A. Dipper, C. N. Dunlop, P. Feautrier, T. Fusco, J. L. Gach, E. Gendron, D. Geng, S. J. Goodsell, D. Gratadour, A. H. Greenaway, A. Guesalaga , et al. (34 additional authors not shown)

Abstract: Recent advances in adaptive optics (AO) have led to the implementation of wide field-of-view AO systems. A number of wide-field AO systems are also planned for the forthcoming Extremely Large Telescopes. Such systems have multiple wavefront sensors of different types, and usually multiple deformable mirrors (DMs). Here, we report on our experience integrating cameras and DMs with the real-time c… ▽ More Recent advances in adaptive optics (AO) have led to the implementation of wide field-of-view AO systems. A number of wide-field AO systems are also planned for the forthcoming Extremely Large Telescopes. Such systems have multiple wavefront sensors of different types, and usually multiple deformable mirrors (DMs). Here, we report on our experience integrating cameras and DMs with the real-time control systems of two wide-field AO systems. These are CANARY, which has been operating on-sky since 2010, and DRAGON, which is a laboratory adaptive optics real-time demonstrator instrument. We detail the issues and difficulties that arose, along with the solutions we developed. We also provide recommendations for consideration when developing future wide-field AO systems. △ Less

Submitted 24 March, 2016; originally announced March 2016.

Comments: Accepted for publication in MNRAS

A tomographic algorithm to determine tip-tilt information from laser guide stars

Authors: A. P. Reeves, T. J. Morris, R. M. Myers, N. A. Bharmal, J. Osborn

Abstract: Laser Guide Stars (LGS) have greatly increased the sky-coverage of Adaptive Optics (AO) systems. Due to the up-link turbulence experienced by LGSs, a Natural Guide Star (NGS) is still required, preventing full sky-coverage. We present a method of obtaining partial tip-tilt information from LGSs alone in multi-LGS tomographic LGS AO systems. The method of LGS up-link tip-tilt determination is deriv… ▽ More Laser Guide Stars (LGS) have greatly increased the sky-coverage of Adaptive Optics (AO) systems. Due to the up-link turbulence experienced by LGSs, a Natural Guide Star (NGS) is still required, preventing full sky-coverage. We present a method of obtaining partial tip-tilt information from LGSs alone in multi-LGS tomographic LGS AO systems. The method of LGS up-link tip-tilt determination is derived using a geometric approach, then an alteration to the Learn and Apply algorithm for tomographic AO is made to accommodate up-link tip-tilt. Simulation results are presented, verifying that the technique shows good performance in correcting high altitude tip-tilt, but not that from low altitudes. We suggest that the method is combined with multiple far off-axis tip-tilt NGSs to provide gains in performance and sky-coverage over current tomographic AO systems. △ Less

Submitted 15 March, 2016; v1 submitted 28 January, 2016; originally announced January 2016.

Comments: Accepted by MNRAS, 11th March 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)

The Q/U Imaging ExperimenT: Polarization Measurements of the Galactic Plane at 43 and 95 GHz

Authors: QUIET Collaboration, T. M. Ruud, U. Fuskeland, I. K. Wehus, M. Vidal, D. Araujo, C. Bischoff, I. Buder, Y. Chinone, K. Cleary, R. N. Dumoulin, A. Kusaka, R. Monsalve, S. K. Naess, L. B. Newburgh, R. A. Reeves, J. T. L. Zwart, L. Bronfman, R. D. Davies, R. Davis, C. Dickinson, H. K. Eriksen, T. Gaier, J. O. Gundersen, M. Hasegawa , et al. (17 additional authors not shown)

Abstract: We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations.… ▽ More We present polarization observations of two Galactic plane fields centered on Galactic coordinates (l,b)=(0 deg,0 deg) and (329 deg, 0 deg) at Q- (43 GHz) and W-band (95 GHz), covering between 301 and 539 square degrees depending on frequency and field. These measurements were made with the QUIET instrument between 2008 October and 2010 December, and include a total of 1263 hours of observations. The resulting maps represent the deepest large-area Galactic polarization observations published to date at the relevant frequencies with instrumental rms noise varying between 1.8 and 2.8 uK deg, 2.3-6 times deeper than corresponding WMAP and Planck maps. The angular resolution is 27.3' and 12.8' FWHM at Q- and W-band, respectively. We find excellent agreement between the QUIET and WMAP maps over the entire fields, and no compelling evidence for significant residual instrumental systematic errors in either experiment, whereas the Planck 44 GHz map deviates from these in a manner consistent with reported systematic uncertainties for this channel. We combine QUIET and WMAP data to compute inverse-variance-weighted average maps, effectively retaining small angular scales from QUIET and large angular scales from WMAP. From these combined maps, we derive constraints on several important astrophysical quantities, including a robust detection of polarized synchrotron spectral index steepening of ~0.2 off the plane, as well as the Faraday rotation measure toward the Galactic center (RM=-4000 +/- 200 rad m^-2), all of which are consistent with previously published results. Both the raw QUIET and the co-added QUIET+WMAP maps are made publicly available together with all necessary ancillary information. △ Less

Submitted 11 August, 2015; originally announced August 2015.

Comments: 22 pages, 18 figures; accepted for publication in ApJ; data are available on Lambda (http://lambda.gsfc.nasa.gov)

Journal ref: ApJ 811, 89 (2015)