-
SKA Science Data Challenge 2: analysis and results
Authors:
P. Hartley,
A. Bonaldi,
R. Braun,
J. N. H. S. Aditya,
S. Aicardi,
L. Alegre,
A. Chakraborty,
X. Chen,
S. Choudhuri,
A. O. Clarke,
J. Coles,
J. S. Collinson,
D. Cornu,
L. Darriba,
M. Delli Veneri,
J. Forbrich,
B. Fraga,
A. Galan,
J. Garrido,
F. Gubanov,
H. Håkansson,
M. J. Hardcastle,
C. Heneka,
D. Herranz,
K. M. Hess
, et al. (83 additional authors not shown)
Abstract:
The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed t…
▽ More
The Square Kilometre Array Observatory (SKAO) will explore the radio sky to new depths in order to conduct transformational science. SKAO data products made available to astronomers will be correspondingly large and complex, requiring the application of advanced analysis techniques to extract key science findings. To this end, SKAO is conducting a series of Science Data Challenges, each designed to familiarise the scientific community with SKAO data and to drive the development of new analysis techniques. We present the results from Science Data Challenge 2 (SDC2), which invited participants to find and characterise 233245 neutral hydrogen (Hi) sources in a simulated data product representing a 2000~h SKA MID spectral line observation from redshifts 0.25 to 0.5. Through the generous support of eight international supercomputing facilities, participants were able to undertake the Challenge using dedicated computational resources. Alongside the main challenge, `reproducibility awards' were made in recognition of those pipelines which demonstrated Open Science best practice. The Challenge saw over 100 participants develop a range of new and existing techniques, with results that highlight the strengths of multidisciplinary and collaborative effort. The winning strategy -- which combined predictions from two independent machine learning techniques to yield a 20 percent improvement in overall performance -- underscores one of the main Challenge outcomes: that of method complementarity. It is likely that the combination of methods in a so-called ensemble approach will be key to exploiting very large astronomical datasets.
△ Less
Submitted 14 March, 2023;
originally announced March 2023.
-
Integration of storage endpoints into a Rucio data lake, as an activity to prototype a SKA Regional Centres Network
Authors:
Manuel Parra-Royón,
Jesús Sánchez-Castañeda,
Julián Garrido,
Susana Sánchez-Expósito,
Rohini Joshi,
James Collinson,
Rob Barnsley,
Jesús Salgado,
Lourdes Verdes-Montenegro
Abstract:
The Square Kilometre Array (SKA) infrastructure will consist of two radio telescopes that will be the most sensitive telescopes on Earth. The SKA community will have to process and manage near exascale data, which will be a technical challenge for the coming years. In this respect, the SKA Global Network of Regional Centres plays a key role in data distribution and management. The SRCNet will prov…
▽ More
The Square Kilometre Array (SKA) infrastructure will consist of two radio telescopes that will be the most sensitive telescopes on Earth. The SKA community will have to process and manage near exascale data, which will be a technical challenge for the coming years. In this respect, the SKA Global Network of Regional Centres plays a key role in data distribution and management. The SRCNet will provide distributed computing and data storage capacity, as well as other important services for the network. Within the SRCNet, several teams have been set up for the research, design and development of 5 prototypes. One of these prototypes is related to data management and distribution, where a data lake has been deployed using Rucio. In this paper we focus on the tasks performed by several of the teams to deploy new storage endpoints within the SKAO data lake. In particular, we will describe the steps and deployment instructions for the services required to provide the Rucio data lake with a new Rucio Storage Element based on StoRM and WebDAV within the Spanish SRC prototype.
△ Less
Submitted 13 March, 2023;
originally announced March 2023.
-
The Changing-Look Quasar Mrk 590 is Awakening
Authors:
S. Mathur,
K. D. Denney,
A. Gupta,
M. Vestergaard,
G. De Rosa,
Y. Krongold,
F. Nicastro,
J. Collinson,
M. Goad,
K. Korista,
R. W. Pogge,
B. M. Peterson
Abstract:
Mrk 590 was originally classified as a Seyfert 1 galaxy, but then it underwent dramatic changes: the nuclear luminosity dropped by over two orders of magnitude and the broad emission lines all but disappeared from the optical spectrum. Here we present followup observations to the original discovery and characterization of this "changing look" active galactic nucleus (AGN). The new Chandra and HST…
▽ More
Mrk 590 was originally classified as a Seyfert 1 galaxy, but then it underwent dramatic changes: the nuclear luminosity dropped by over two orders of magnitude and the broad emission lines all but disappeared from the optical spectrum. Here we present followup observations to the original discovery and characterization of this "changing look" active galactic nucleus (AGN). The new Chandra and HST observations from 2014 show that Mrk 590 is awakening, changing its appearance again. While the source continues to be in a low state, its soft excess has re-emerged, though not to the previous level. The UV continuum is brighter by more than a factor of two and the broad MgII emission line is present, indicating that the ionizing continuum is also brightening. These observations suggest that the soft excess is not due to reprocessed hard X-ray emission. Instead, it is connected to the UV continuum through warm Comptonization. Variability of the Fe K-alpha emission lines suggests that the reprocessing region is within about 10 light years or 3 pc of the central source. The AGN type change is neither due to obscuration, nor due to one-way evolution from type-1 to type-2, as suggested in literature, but may be related to episodic accretion events.
△ Less
Submitted 15 October, 2018;
originally announced October 2018.
-
The broad band SEDs of four 'hypervariable' AGN
Authors:
James S. Collinson,
Martin J. Ward,
Andy Lawrence,
Alastair Bruce,
Chelsea L. MacLeod,
Martin Elvis,
Suvi Gezari,
Philip J. Marshall,
Chris Done
Abstract:
We present an optical to X-ray spectral analysis of four 'hypervariable' AGN (HVAs) discovered by comparing Pan-STARRS data to that from SDSS over a 10 year baseline (Lawrence et al 2016). There is some evidence that these objects are X-ray loud for their corresponding UV luminosities, but given that we measured them in a historic high state, it is not clear whether to take the high-state or low-s…
▽ More
We present an optical to X-ray spectral analysis of four 'hypervariable' AGN (HVAs) discovered by comparing Pan-STARRS data to that from SDSS over a 10 year baseline (Lawrence et al 2016). There is some evidence that these objects are X-ray loud for their corresponding UV luminosities, but given that we measured them in a historic high state, it is not clear whether to take the high-state or low-state as typical of the properties of these HVAs. We estimate black hole masses based on MgII and H$α$ emission line profiles, and either the high and low state luminosities, finding mass ranges $\log(M_{\rm BH}/M_{\odot}) = 8.2-8.8$ and $\log(M_{\rm BH}/M_{\odot}) = 7.9-8.3$ respectively. We then fit energy conserving models to the SEDs, obtaining strong constraints on the bolometric luminosity and $α_{\rm OX}$. We compare the SED properties with a larger, X-ray selected AGN sample for both of these scenarios, and observe distinct groupings in spectral shape versus luminosity parameter space. In general, the SED properties are closer to normal if we assume that the low-state is representative. This supports the idea that the large slow outbursts may be due to extrinsic effects (for example microlensing) as opposed to accretion rate changes, but a larger sample of HVAs is needed to be confident of this conclusion.
△ Less
Submitted 20 November, 2017;
originally announced November 2017.
-
Reaching the Peak of the quasar spectral energy distribution - II. Exploring the accretion disc, dusty torus and host galaxy
Authors:
James S. Collinson,
Martin J. Ward,
Hermine Landt,
Chris Done,
Martin Elvis,
Jonathan C. McDowell
Abstract:
We continue our study of the spectral energy distributions (SEDs) of 11 AGN at 1.5 < z < 2.2, with optical-NIR spectra, X-ray data and mid-IR photometry. In a previous paper we presented the observations and models; in this paper we explore the parameter space of these models. We first quantify uncertainties on the black hole masses (M$_{\rm BH}$) and degeneracies between SED parameters. The effec…
▽ More
We continue our study of the spectral energy distributions (SEDs) of 11 AGN at 1.5 < z < 2.2, with optical-NIR spectra, X-ray data and mid-IR photometry. In a previous paper we presented the observations and models; in this paper we explore the parameter space of these models. We first quantify uncertainties on the black hole masses (M$_{\rm BH}$) and degeneracies between SED parameters. The effect of BH spin is tested, and we find that while low to moderate spin values (a$_*$ $\leq$ 0.9) are compatible with the data in all cases, maximal spin (a$_*$ = 0.998) can only describe the data if the accretion disc is face-on. The outer accretion disc radii are well constrained in 8/11 objects, and are found to be a factor ~5 smaller than the self-gravity radii. We then extend our modelling campaign into the mid-IR regime with WISE photometry, adding components for the host galaxy and dusty torus. Our estimates of the host galaxy luminosities are consistent with the M$_{\rm BH}$-bulge relationship, and the measured torus properties (covering factor and temperatures) are in agreement with earlier work, suggesting a predominantly silicate-based grain composition. Finally, we deconvolve the optical-NIR spectra using our SED continuum model. We claim that this is a more physically motivated approach than using empirical descriptions of the continuum such as broken power-laws. For our small sample, we verify previously noted correlations between emission linewidths and luminosities commonly used for single-epoch M$_{\rm BH}$ estimates, and observe a statistically significant anti-correlation between [O III] equivalent width and AGN luminosity.
△ Less
Submitted 13 October, 2016;
originally announced October 2016.
-
Spectral analysis of four 'hypervariable' AGN: a microneedle in the haystack?
Authors:
A. Bruce,
A. Lawrence,
C. MacLeod,
M. Elvis,
M. J. Ward,
J. S. Collinson,
S. Gezari,
P. J. Marshall,
M. C. Lam,
R. Kotak,
C. Inserra,
J. Polshaw,
N. Kaiser,
R-P. Kudritzki,
E. A. Magnier,
C. Waters
Abstract:
We analyze four extreme active galactic nuclei (AGN) transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 mag in the g-band on a time-scale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differentia…
▽ More
We analyze four extreme active galactic nuclei (AGN) transients to explore the possibility that they are caused by rare, high-amplitude microlensing events. These previously unknown type-I AGN are located in the redshift range 0.6-1.1 and show changes of > 1.5 mag in the g-band on a time-scale of ~years. Multi-epoch optical spectroscopy, from the William Herschel Telescope, shows clear differential variability in the broad line fluxes with respect to the continuum changes and also evolution in the line profiles. In two cases, a simple point-source, point-lens microlensing model provides an excellent match to the long-term variability seen in these objects. For both models, the parameter constraints are consistent with the microlensing being due to an intervening stellar mass object but as yet there is no confirmation of the presence of an intervening galaxy. The models predict a peak amplification of 10.3/13.5 and an Einstein time-scale of 7.5/10.8 yr, respectively. In one case, the data also allow constraints on the size of the CIII] emitting region, with some simplifying assumptions, to to be ~1.0-6.5 light-days and a lower limit on the size of the MgII emitting region to be > 9 light-days (half-light radii). This CIII] radius is perhaps surprisingly small. In the remaining two objects, there is spectroscopic evidence for an intervening absorber but the extra structure seen in the light curves requires a more complex lensing scenario to adequately explain.
△ Less
Submitted 3 March, 2017; v1 submitted 22 August, 2016;
originally announced August 2016.
-
Physical properties and the variability mechanism of the He I outflow in NGC 4151
Authors:
Conor Wildy,
Hermine Landt,
Michael R. Goad,
Martin J. Ward,
James S. Collinson
Abstract:
We report on variable helium absorption lines in NGC 4151 observed across six epochs of quasi-simultaneous near-infrared and optical data. These observations cover the transitions from the metastable 2^3S state at 3889 A and 10830 A, and from the 2^1S state at 20587 A. This is the first AGN absorption line variability study to include measurements of the 20587 A line. The physical properties of th…
▽ More
We report on variable helium absorption lines in NGC 4151 observed across six epochs of quasi-simultaneous near-infrared and optical data. These observations cover the transitions from the metastable 2^3S state at 3889 A and 10830 A, and from the 2^1S state at 20587 A. This is the first AGN absorption line variability study to include measurements of the 20587 A line. The physical properties of the absorber recorded at the fifth observational epoch are relatively well constrained by the presence of absorption in both the optical and near-infrared components, with the 10830 A line likely saturated. The observations suggest variations in this absorber's strength are best explained by ionization changes in response to a variable incident continuum. Photoionization simulations constrain the total hydrogen number density of the epoch 5 absorber to 7.1<log(n_H/cm^-3)<8.8, the hydrogen column density to 21.2<log(N_H/cm^-2)<23.3 and the ionization parameter range to -1.9<logU<0.4. The simulations also suggest the absorber is located between 0.03 and 0.49 pc from the continuum emission region. This range in physical properties is consistent with an absorber of similar velocity seen in NGC 4151 from previous ultraviolet and optical studies, but with high column density X-ray absorbing components not present. The mass outflow rate due to the fifth epoch absorber is in the range 0.008 to 0.38 M_sun/yr, too low to contribute to galaxy feedback effects.
△ Less
Submitted 10 June, 2016;
originally announced June 2016.
-
Reaching the Peak of the Quasar Spectral Energy Distribution I: Observations and Models
Authors:
James S. Collinson,
Martin J. Ward,
Chris Done,
Hermine Landt,
Martin Elvis,
Jonathan C. McDowell
Abstract:
We perform a spectral analysis of a sample of 11 medium redshift (1.5 < z < 2.2) quasars. Our sample all have optical spectra from the SDSS, infrared spectra from GNIRS and TSPEC, and X-ray spectra from XMM-Newton. We first analyse the Balmer broad emission line profiles which are shifted into the IR spectra to constrain black hole masses. Then we fit an energy-conserving, three component accretio…
▽ More
We perform a spectral analysis of a sample of 11 medium redshift (1.5 < z < 2.2) quasars. Our sample all have optical spectra from the SDSS, infrared spectra from GNIRS and TSPEC, and X-ray spectra from XMM-Newton. We first analyse the Balmer broad emission line profiles which are shifted into the IR spectra to constrain black hole masses. Then we fit an energy-conserving, three component accretion model of the broadband spectral energy distribution (SED) to our multi-wavelength data. Five out of the 11 quasars show evidence of an SED peak, allowing us to constrain their bolometric luminosity from these models and estimate their mass accretion rates. Based on our limited sample, we suggest that estimating bolometric luminosities from L_5100A and L_2-10keV may be unreliable, as has been also noted for a low-redshift, X-ray selected AGN sample.
△ Less
Submitted 20 February, 2015;
originally announced February 2015.