The Hyper Suprime-Cam extended Point Spread Functions and applications to measuring the intra-halo light

We present extended point spread function (PSF) models for the Hyper Suprime-Cam Subaru Strategic Program Public Data Release 3 (HSC-SSP PDR3) in all $\textit{g,r,i,Z}$ and $\textit{Y}$-bands. Due to its 8.2m primary mirror and long exposure periods, HSC combines deep images with wide-field coverage, making it one of the most suitable observing facilities for low surface brightness (LSB) studies. By applying a median stacking technique of point sources with different brightnesses, we show how to construct the HSC-SSP PDR3 PSF models to an extent of R $\sim$ 5.6 arcmin. These new PSFs provide the community with a crucial tool to characterise LSB properties at large angles. We apply our HSC PSFs and demonstrate that they behave reasonably in two cases: first, to generate a 2-D model of a bright star, and second, to remove the PSF-scattered light from an Ultra Deep image of the 400020 Galaxy And Mass Assembly (GAMA) group in the SXDS field. Our main focus in this second application is characterising the $\textit{r}$-band intra-halo light (IHL) component of 400020. Building on advanced source extraction techniques with careful consideration of PSF flux, we measure the IHL surface brightness (SB) group profile up to $\sim$ 31 mag arcsec$^{-2}$ and R = 300 kpc. We estimate the IHL fraction ($\mathrm{f_{IHL}}$) profile, with a mean of $\mathrm{f_{IHL}}$ $\sim$ 0.13. Our results show that not removing the PSF light can overestimate the IHL SB by $\sim$ 1.7 mag arcsec$^{-2}$ and the $\mathrm{f_{IHL}}$ by $\sim$ 30%.Comment: 21 pages, 19 figures, submitted to MNRA

Preparing for low surface brightness science with the Vera C. Rubin Observatory: A Comparison of Observable and Simulated Intracluster Light Fractions

Intracluster Light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (Surface Brightness Threshold-SB, Non-Parametric Measure-NP, Composite Models-CM, Multi-Galaxy Fitting-MGF) and new approaches under development (Wavelet Decomposition-WD) applied to mock images of 61 galaxy clusters (14<log10 M_200c/M_solar <14.5) from four cosmological hydrodynamical simulations. We compare the BCG+ICL and ICL fractions from observational measures with those using simulated measures (aperture and kinematic separations). The ICL fractions measured by kinematic separation are significantly larger than observed fractions. We find the measurements are related and provide equations to estimate kinematic ICL fractions from observed fractions. The different observational techniques give consistent BCG+ICL and ICL fractions but are biased to underestimating the BCG+ICL and ICL fractions when compared with aperture simulation measures. Comparing the different methods and algorithms we find that the MGF algorithm is most consistent with the simulations, and CM and SB methods show the smallest projection effects for the BCG+ICL and ICL fractions respectively. The Ahad (CM), MGF and WD algorithms are best set up to process larger samples, however, the WD algorithm in its current form is susceptible to projection effects. We recommend that new algorithms using these methods are explored to analyse the massive samples that Rubin Observatory's Legacy Survey of Space and Time will provide.Comment: Submitted for publication in MNRAS, posted to arXiv after responding to two positive rounds of referee comments. Key results in Figs 3, 5, 6 and 1

Preparing for low surface brightness science with the Vera C. Rubin Observatory: a comparison of observable and simulated intracluster light fractions

Intracluster light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this, we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (surface brightness threshold-SB, non-parametric measure-NP, composite models-CM, and multi-galaxy fitting-MGF) and new approaches under development (wavelet decomposition-WD) applied to mock images of 61 galaxy clusters (14 <log10M200c/M⊙ < 14.5) from four cosmological hydrodynamical simulations. We compare the BCG+ICL and ICL fractions from observational measures with those using simulated measures (aperture and kinematic separations). The ICL fractions measured by kinematic separation are significantly larger than observed fractions. We find the measurements are related and provide equations to estimate kinematic ICL fractions from observed fractions. The different observational techniques give consistent BCG+ICL and ICL fractions but are biased to underestimating the BCG+ICL and ICL fractions when compared with aperture simulation measures. Comparing the different methods and algorithms, we find that the MGF algorithm is most consistent with the simulations, and CM and SB methods show the smallest projection effects for the BCG+ICL and ICL fractions, respectively. The Ahad (CM), MGF, and WD algorithms are best set up to process larger samples; however, the WD algorithm in its current form is susceptible to projection effects. We recommend that new algorithms using these methods are explored to analyse the massive samples that Rubin Observatory’s Legacy Survey of Space and Time will provide

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

We report on long-term multiwavelengthmonitoring of blazar Mrk 421 by the GLAST-AGILE Support Program of the Whole Earth Blazar Telescope (GASP-WEBT) collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007–2015, characterized by several extreme flares. The ratio between the optical, X-ray and γ -ray fluxes is very variable. The γ -ray flux variations show a fair correlation with the optical ones starting from 2012.We analyse spectropolarimetric data and find wavelengthdependence of the polarization degree (P), which is compatible with the presence of the host galaxy, and no wavelength dependence of the electric vector polarization angle (EVPA). Optical polarimetry shows a lack of simple correlation between P and flux and wide rotations of the EVPA.We build broad-band spectral energy distributions with simultaneous near-infrared and optical data from the GASP-WEBT and ultraviolet and X-ray data from the Swift satellite. They show strong variability in both flux and X-ray spectral shape and suggest a shift of the synchrotron peak up to a factor of ∼50 in frequency. The interpretation of the flux and spectral variability is compatible with jet models including at least two emitting regions that can change their orientation with respect to the line of sight.http://10.0.4.69/mnras/stx2185Accepted manuscrip

Preparing for low surface brightness science with the Vera C. Rubin Observatory: a comparison of observable and simulated intracluster light fractions

Intracluster light (ICL) provides an important record of the interactions galaxy clusters have undergone. However, we are limited in our understanding by our measurement methods. To address this, we measure the fraction of cluster light that is held in the Brightest Cluster Galaxy and ICL (BCG+ICL fraction) and the ICL alone (ICL fraction) using observational methods (surface brightness threshold-SB, non-parametric measure-NP, composite models-CM, and multi-galaxy fitting-MGF) and new approaches under development (wavelet decomposition-WD) applied to mock images of 61 galaxy clusters (14 <log10M200c/M☉ < 14.5) from four cosmological hydrodynamical simulations. We compare the BCG+ICL and ICL fractions from observational measures with those using simulated measures (aperture and kinematic separations). The ICL fractions measured by kinematic separation are significantly larger than observed fractions. We find the measurements are related and provide equations to estimate kinematic ICL fractions from observed fractions. The different observational techniques give consistent BCG+ICL and ICL fractions but are biased to underestimating the BCG+ICL and ICL fractions when compared with aperture simulation measures. Comparing the different methods and algorithms, we find that the MGF algorithm is most consistent with the simulations, and CM and SB methods show the smallest projection effects for the BCG+ICL and ICL fractions, respectively. The Ahad (CM), MGF, and WD algorithms are best set up to process larger samples; however, the WD algorithm in its current form is susceptible to projection effects. We recommend that new algorithms using these methods are explored to analyse the massive samples that Rubin Observatory's Legacy Survey of Space and Time will provide

Preparing for low surface brightness science with the Vera C. Rubin Observatory:Characterization of tidal features from mock images

Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilizing automated techniques and human visual classification in conjunction with realistic mock images produced using the NewHorizon cosmological simulation, we investigate the nature, frequency, and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-yr depth of the Legacy Survey of Space and Time (30-31 mag arcsec-2), falling to 60 per cent assuming a shallower final depth of 29.5 mag arcsec-2. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M* ∼1011.5 M⊙). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterization of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimized, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z < 0.2)

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

© 2017 The Authors Published by Oxford University Press on behalf of the Royal Astronomical Society. We report on long-term multiwavelength monitoring of blazar Mrk 421 by the GLAST-AGILE Support Program of the Whole Earth Blazar Telescope (GASP-WEBT) collaboration and Steward Observatory, and by the Swift and Fermi satellites. We study the source behaviour in the period 2007-2015, characterized by several extreme flares. The ratio between the optical, X-ray and γ-ray fluxes is very variable. The γ-ray flux variations show a fair correlation with the optical ones starting from 2012. We analyse spectropolarimetric data and find wavelength-dependence of the polarization degree (P), which is compatible with the presence of the host galaxy, and no wavelength dependence of the electric vector polarization angle (EVPA). Optical polarimetry shows a lack of simple correlation between P and flux and wide rotations of the EVPA. We build broad-band spectral energy distributions with simultaneous near-infrared and optical data from the GASP-WEBT and ultraviolet and X-ray data from the Swift satellite. They show strong variability in both flux and X-ray spectr al shape and suggest a shift of the synchrotron peak up to a factor of ~50 in frequency. The interpretation of the flux and spectral variability is compatible with jet models including at least two emitting regions that can change their orientation with respect to the line of sight

Synchrotron emission from the blazar PG 1553+113. An analysis of its flux and polarization variability

© 2016 The Authors.In 2015 July 29-September 1, the satellite XMM-Newton pointed at the BL Lac object PG 1553+133 six times, collecting data for 218 h. During one of these epochs, simultaneous observations by the Swift satellite were requested to compare the results of the X-ray and optical-UV instruments. Optical, near-infrared and radio monitoring was carried out by the Whole Earth Blazar Telescope (WEBT) collaboration for the whole observing season. We here present the results of the analysis of all these data, together with an investigation of the source photometric and polarimetric behaviour over the last 3 yr. The 2015 EPIC spectra show slight curvature and the corresponding light curves display fast X-ray variability with a time-scale of the order of 1 h. In contrast to previous results, during the brightest X-ray states detected in 2015 the simple log-parabolic model that best fits the XMM-Newton data also reproduces reasonably well the whole synchrotron bump, suggesting a peak in the near- UV band. We found evidence of a wide rotation of the polarization angle in 2014, when the polarization degree was variable, but the flux remained almost constant. This is difficult to interpret with deterministic jet emission models, while it can be easily reproduced by assuming some turbulence of the magnetic field

Preparing for low surface brightness science with the Vera C. Rubin Observatory: characterisation of tidal features from mock images

Tidal features in the outskirts of galaxies yield unique information about their past interactions and are a key prediction of the hierarchical structure formation paradigm. The Vera C. Rubin Observatory is poised to deliver deep observations for potentially of millions of objects with visible tidal features, but the inference of galaxy interaction histories from such features is not straightforward. Utilising automated techniques and human visual classification in conjunction with realistic mock images produced using the NEWHORIZON cosmological simulation, we investigate the nature, frequency and visibility of tidal features and debris across a range of environments and stellar masses. In our simulated sample, around 80 per cent of the flux in the tidal features around Milky Way or greater mass galaxies is detected at the 10-year depth of the Legacy Survey of Space and Time (30-31 mag / sq. arcsec), falling to 60 per cent assuming a shallower final depth of 29.5 mag / sq. arcsec. The fraction of total flux found in tidal features increases towards higher masses, rising to 10 per cent for the most massive objects in our sample (M*~10^{11.5} Msun). When observed at sufficient depth, such objects frequently exhibit many distinct tidal features with complex shapes. The interpretation and characterisation of such features varies significantly with image depth and object orientation, introducing significant biases in their classification. Assuming the data reduction pipeline is properly optimised, we expect the Rubin Observatory to be capable of recovering much of the flux found in the outskirts of Milky Way mass galaxies, even at intermediate redshifts (z<0.2)

Multiwavelength observations of a VHE gamma-ray flare from PKS 1510-089 in 2015

Context. PKS 1510-089 is one of only a few flat spectrum radio quasars detected in the VHE (very-high-energy, > 100 GeV) gamma-ray band. Aims. We study the broadband spectral and temporal properties of the PKS 1510-089 emission during a high gamma-ray state. Methods. We performed VHE gamma-ray observations of PKS 1510-089 with the MAGIC telescopes during a high gamma-ray state in May 2015. In order to perform broad-band modelling of the source, we have also gathered contemporaneous multiwavelength data in radio, IR, optical photometry and polarization, UV, X-ray and GeV gamma-ray ranges. We construct a broadband spectral energy distribution (SED) in two periods, selected according to VHE gamma-ray state. Results. PKS 1510-089 has been detected in a high optical and gamma-ray state, showing for the first time a significant VHE gamma-ray variability. Similarly to the optical and gamma-ray high state of the source detected in 2012, it was accompanied by a rotation of the optical polarization angle and the emission of a new jet component observed in radio. The spectral shape in the VHE band during the flare is similar to the ones obtained during previous measurements of the source. The observed flux variability sets for the first time constraints on the size of the region from which VHE gamma rays are emitted. The broadband SED can be explained in the External Compton scenario