Showing 1–43 of 43 results for author: Sawangwit, U

$\textit{Kilonova Seekers}$: the GOTO project for real-time citizen science in time-domain astrophysics

Authors: T. L. Killestein, L. Kelsey, E. Wickens, L. Nuttall, J. Lyman, C. Krawczyk, K. Ackley, M. J. Dyer, F. Jim�nez-Ibarra, K. Ulaczyk, D. O'Neill, A. Kumar, D. Steeghs, D. K. Galloway, V. S. Dhillon, P. O'Brien, G. Ramsay, K. Noysena, R. Kotak, R. P. Breton, E. Pall�, D. Pollacco, S. Awiphan, S. Belkin, P. Chote , et al. (29 additional authors not shown)

Abstract: Time-domain astrophysics continues to grow rapidly, with the inception of new surveys drastically increasing data volumes. Democratised, distributed approaches to training sets for machine learning classifiers are crucial to make the most of this torrent of discovery -- with citizen science approaches proving effective at meeting these requirements. In this paper, we describe the creation of and t… ▽ More Time-domain astrophysics continues to grow rapidly, with the inception of new surveys drastically increasing data volumes. Democratised, distributed approaches to training sets for machine learning classifiers are crucial to make the most of this torrent of discovery -- with citizen science approaches proving effective at meeting these requirements. In this paper, we describe the creation of and the initial results from the $\textit{Kilonova Seekers}$ citizen science project, built to find transient phenomena from the GOTO telescopes in near real-time. $\textit{Kilonova Seekers}$ launched in July 2023 and received over 600,000 classifications from approximately 2,000 volunteers over the course of the LIGO-Virgo-KAGRA O4a observing run. During this time, the project has yielded 20 discoveries, generated a `gold-standard' training set of 17,682 detections for augmenting deep-learned classifiers, and measured the performance and biases of Zooniverse volunteers on real-bogus classification. This project will continue throughout the lifetime of GOTO, pushing candidates at ever-greater cadence, and directly facilitate the next-generation classification algorithms currently in development. △ Less

Submitted 24 July, 2024; v1 submitted 4 June, 2024; originally announced June 2024.

Comments: 20 pages, 15 figures. Accepted in MNRAS

Real-time Monitoring for the Next Core-Collapse Supernova in JUNO

Authors: Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, Jo�o Pedro Athayde Marcondes de Andr�, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato, Marco Beretta, Antonio Bergnoli , et al. (606 additional authors not shown)

Abstract: The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neu… ▽ More The core-collapse supernova (CCSN) is considered one of the most energetic astrophysical events in the universe. The early and prompt detection of neutrinos before (pre-SN) and during the supernova (SN) burst presents a unique opportunity for multi-messenger observations of CCSN events. In this study, we describe the monitoring concept and present the sensitivity of the system to pre-SN and SN neutrinos at the Jiangmen Underground Neutrino Observatory (JUNO), a 20 kton liquid scintillator detector currently under construction in South China. The real-time monitoring system is designed to ensure both prompt alert speed and comprehensive coverage of progenitor stars. It incorporates prompt monitors on the electronic board as well as online monitors at the data acquisition stage. Assuming a false alert rate of 1 per year, this monitoring system exhibits sensitivity to pre-SN neutrinos up to a distance of approximately 1.6 (0.9) kiloparsecs and SN neutrinos up to about 370 (360) kiloparsecs for a progenitor mass of 30 solar masses, considering both normal and inverted mass ordering scenarios. The pointing ability of the CCSN is evaluated by analyzing the accumulated event anisotropy of inverse beta decay interactions from pre-SN or SN neutrinos. This, along with the early alert, can play a crucial role in facilitating follow-up multi-messenger observations of the next galactic or nearby extragalactic CCSN. △ Less

Submitted 4 December, 2023; v1 submitted 13 September, 2023; originally announced September 2023.

Comments: 24 pages, 9 figures, accepted for the publication at JCAP

JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo

Authors: JUNO Collaboration, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, Jo�o Pedro Athayde Marcondes de Andr�, Didier Auguste, Weidong Bai, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato , et al. (581 additional authors not shown)

Abstract: We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon… ▽ More We discuss JUNO sensitivity to the annihilation of MeV dark matter in the galactic halo via detecting inverse beta decay reactions of electron anti-neutrinos resulting from the annihilation. We study possible backgrounds to the signature, including the reactor neutrinos, diffuse supernova neutrino background, charged- and neutral-current interactions of atmospheric neutrinos, backgrounds from muon-induced fast neutrons and cosmogenic isotopes. A fiducial volume cut, as well as the pulse shape discrimination and the muon veto are applied to suppress the above backgrounds. It is shown that JUNO sensitivity to the thermally averaged dark matter annihilation rate in 10 years of exposure would be significantly better than the present-day best limit set by Super-Kamiokande and would be comparable to that expected by Hyper-Kamiokande. △ Less

Submitted 13 September, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: 25 pages, 9 figures, matches the publised version

Journal ref: JCAP 09 (2023) 001

Model Independent Approach of the JUNO $^8$B Solar Neutrino Program

Authors: JUNO Collaboration, Jie Zhao, Baobiao Yue, Haoqi Lu, Yufeng Li, Jiajie Ling, Zeyuan Yu, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Abid Aleem, Tsagkarakis Alexandros, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Weidong Bai , et al. (579 additional authors not shown)

Abstract: The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low backg… ▽ More The physics potential of detecting $^8$B solar neutrinos will be exploited at the Jiangmen Underground Neutrino Observatory (JUNO), in a model independent manner by using three distinct channels of the charged-current (CC), neutral-current (NC) and elastic scattering (ES) interactions. Due to the largest-ever mass of $^{13}$C nuclei in the liquid-scintillator detectors and the {expected} low background level, $^8$B solar neutrinos would be observable in the CC and NC interactions on $^{13}$C for the first time. By virtue of optimized event selections and muon veto strategies, backgrounds from the accidental coincidence, muon-induced isotopes, and external backgrounds can be greatly suppressed. Excellent signal-to-background ratios can be achieved in the CC, NC and ES channels to guarantee the $^8$B solar neutrino observation. From the sensitivity studies performed in this work, we show that JUNO, with ten years of data, can reach the {1$σ$} precision levels of 5%, 8% and 20% for the $^8$B neutrino flux, $\sin^2θ_{12}$, and $Δm^2_{21}$, respectively. It would be unique and helpful to probe the details of both solar physics and neutrino physics. In addition, when combined with SNO, the world-best precision of 3% is expected for the $^8$B neutrino flux measurement. △ Less

Submitted 6 March, 2024; v1 submitted 15 October, 2022; originally announced October 2022.

Comments: 19 pages, 7 figures, accepted version to appear in The Astrophysical Journal. Yufeng Li and Jiajie Ling are corresponding authors

Journal ref: Astrophysical Journal 965 (2024) 122

Self-Supervised Clustering on Image-Subtracted Data with Deep-Embedded Self-Organizing Map

Authors: Y. -L. Mong, K. Ackley, T. L. Killestein, D. K. Galloway, M. Dyer, R. Cutter, M. J. I. Brown, J. Lyman, K. Ulaczyk, D. Steeghs, V. Dhillon, P. O'Brien, G. Ramsay, K. Noysena, R. Kotak, R. Breton, L. Nuttall, E. Palle, D. Pollacco, E. Thrane, S. Awiphan, U. Burhanudin, P. Chote, A. Chrimes, E. Daw , et al. (23 additional authors not shown)

Abstract: Developing an effective automatic classifier to separate genuine sources from artifacts is essential for transient follow-ups in wide-field optical surveys. The identification of transient detections from the subtraction artifacts after the image differencing process is a key step in such classifiers, known as real-bogus classification problem. We apply a self-supervised machine learning model, th… ▽ More Developing an effective automatic classifier to separate genuine sources from artifacts is essential for transient follow-ups in wide-field optical surveys. The identification of transient detections from the subtraction artifacts after the image differencing process is a key step in such classifiers, known as real-bogus classification problem. We apply a self-supervised machine learning model, the deep-embedded self-organizing map (DESOM) to this "real-bogus" classification problem. DESOM combines an autoencoder and a self-organizing map to perform clustering in order to distinguish between real and bogus detections, based on their dimensionality-reduced representations. We use 32x32 normalized detection thumbnails as the input of DESOM. We demonstrate different model training approaches, and find that our best DESOM classifier shows a missed detection rate of 6.6% with a false positive rate of 1.5%. DESOM offers a more nuanced way to fine-tune the decision boundary identifying likely real detections when used in combination with other types of classifiers, for example built on neural networks or decision trees. We also discuss other potential usages of DESOM and its limitations. △ Less

Submitted 13 September, 2022; originally announced September 2022.

Prospects for Detecting the Diffuse Supernova Neutrino Background with JUNO

Authors: JUNO Collaboration, Angel Abusleme, Thomas Adam, Shakeel Ahmad, Rizwan Ahmed, Sebastiano Aiello, Muhammad Akram, Fengpeng An, Qi An, Giuseppe Andronico, Nikolay Anfimov, Vito Antonelli, Tatiana Antoshkina, Burin Asavapibhop, João Pedro Athayde Marcondes de André, Didier Auguste, Nikita Balashov, Wander Baldini, Andrea Barresi, Davide Basilico, Eric Baussan, Marco Bellato, Antonio Bergnoli, Thilo Birkenfeld, Sylvie Blin , et al. (577 additional authors not shown)

Abstract: We present the detection potential for the diffuse supernova neutrino background (DSNB) at the Jiangmen Underground Neutrino Observatory (JUNO), using the inverse-beta-decay (IBD) detection channel on free protons. We employ the latest information on the DSNB flux predictions, and investigate in detail the background and its reduction for the DSNB search at JUNO. The atmospheric neutrino induced n… ▽ More We present the detection potential for the diffuse supernova neutrino background (DSNB) at the Jiangmen Underground Neutrino Observatory (JUNO), using the inverse-beta-decay (IBD) detection channel on free protons. We employ the latest information on the DSNB flux predictions, and investigate in detail the background and its reduction for the DSNB search at JUNO. The atmospheric neutrino induced neutral current (NC) background turns out to be the most critical background, whose uncertainty is carefully evaluated from both the spread of model predictions and an envisaged \textit{in situ} measurement. We also make a careful study on the background suppression with the pulse shape discrimination (PSD) and triple coincidence (TC) cuts. With latest DSNB signal predictions, more realistic background evaluation and PSD efficiency optimization, and additional TC cut, JUNO can reach the significance of 3$σ$ for 3 years of data taking, and achieve better than 5$σ$ after 10 years for a reference DSNB model. In the pessimistic scenario of non-observation, JUNO would strongly improve the limits and exclude a significant region of the model parameter space. △ Less

Submitted 13 October, 2022; v1 submitted 18 May, 2022; originally announced May 2022.

Comments: 29 pages, 11 figures, final published version in JCAP

Journal ref: JCAP 10 (2022) 033

CTA Sensitivity on TeV scale Dark Matter Models with Complementary Limits from Direct Detection

Authors: C. Duangchan, C. Pongkitivanichkul, P. Uttayarat, A. Jardin-Blicq, M. Wechakama, T. Klangburam, W. Treesukrat, D. Samart, U. Sawangwit, A. Aguirre-Santaella, M. A. Sánchez-Conde

Abstract: With ever increasing pressure from collider physics and direct detection experiments, particle physics models of TeV scale dark matter are gaining more attention. In this work, we consider two realizations of the class of scalar portal dark matter scenarios -- the inverse seesaw model and the inert doublet model. Observations by the Cherenkov Telescope Array (CTA) of very-high-energy $γ$ rays from… ▽ More With ever increasing pressure from collider physics and direct detection experiments, particle physics models of TeV scale dark matter are gaining more attention. In this work, we consider two realizations of the class of scalar portal dark matter scenarios -- the inverse seesaw model and the inert doublet model. Observations by the Cherenkov Telescope Array (CTA) of very-high-energy $γ$ rays from dark matter annihilation in the context of these models are simulated for the Draco and Sculptor dwarf spheroidal galaxies, and later analyzed using ctools. We study the potential of CTA for the 5$σ$ detection of a dark matter annihilation signal. In the absence of a signal, we also derive the 2$σ$ upper limits on the annihilation cross-section. We compare our projected CTA sensitivity against the projected sensitivity of the next generation of direct detection experiment, i.e. XENONnT. Although the limits from CTA are significantly improved compared with the previous generations of $γ$-ray experiments, they are still $\sim2$ orders of magnitude above the thermal relic cross-section for the considered targets. In the case of the inverse seesaw model, the constraint from the future direct detection experiment XENONnT is much weaker than the CTA sensitivity, whereas for the inert doublet model, XENONnT gives a bound an order of magnitude stronger compared to the CTA limits. △ Less

Submitted 27 May, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 18 pages, 6 figures, 3 tables

Journal ref: JCAP 05 (2022) 038

The Gravitational-wave Optical Transient Observer (GOTO): prototype performance and prospects for transient science

Authors: D. Steeghs, D. K. Galloway, K. Ackley, M. J. Dyer, J. Lyman, K. Ulaczyk, R. Cutter, Y. L. Mong, V. Dhillon, P. O'Brien, G. Ramsay, S. Poshyachinda, R. Kotak, L. K. Nuttall, E. Palle, R. P. Breton, D. Pollacco, E. Thrane, S. Aukkaravittayapun, S. Awiphan, U. Burhanudin, P. Chote, A. Chrimes, E. Daw, C. Duffy , et al. (28 additional authors not shown)

Abstract: The Gravitational-wave Optical Transient Observer (GOTO) is an array of wide-field optical telescopes, designed to exploit new discoveries from the next generation of gravitational wave detectors (LIGO, Virgo, KAGRA), study rapidly evolving transients, and exploit multi-messenger opportunities arising from neutrino and very high energy gamma-ray triggers. In addition to a rapid response mode, the… ▽ More The Gravitational-wave Optical Transient Observer (GOTO) is an array of wide-field optical telescopes, designed to exploit new discoveries from the next generation of gravitational wave detectors (LIGO, Virgo, KAGRA), study rapidly evolving transients, and exploit multi-messenger opportunities arising from neutrino and very high energy gamma-ray triggers. In addition to a rapid response mode, the array will also perform a sensitive, all-sky transient survey with few day cadence. The facility features a novel, modular design with multiple 40-cm wide-field reflectors on a single mount. In June 2017 the GOTO collaboration deployed the initial project prototype, with 4 telescope units, at the Roque de los Muchachos Observatory (ORM), La Palma, Canary Islands. Here we describe the deployment, commissioning, and performance of the prototype hardware, and discuss the impact of these findings on the final GOTO design. We also offer an initial assessment of the science prospects for the full GOTO facility that employs 32 telescope units across two sites. △ Less

Submitted 11 October, 2021; originally announced October 2021.

Comments: 19 pages, 16 Figures, accepted for publication in Monthly Notices of the Royal Astronomical Society

Eclipse timing variations in the WD+dM eclipsing binary RR Cae

Authors: R. Rattanamala, S. Awiphan, S. Komonjinda, A. Phriksee, P. Sappankum, N. A-thano, S. Chitchak, P. Rittipruk, U. Sawangwit, S. Poshyachinda, D. E. Reichart, J. B. Haislip

Abstract: We present the binary model and the eclipse timing variations of the eclipsing binary RR Cae, which consists of a white dwarf eclipsed by an M-type dwarf companion. The multi-wavelength optical photometry from the Transiting Exoplanet Survey Satellite (TESS), the 0.6-m PROMPT-8 telescope, and the 0.7-m Thai Robotic Telescope at Spring Brook Observatory, combined with archive H-alpha radial velocit… ▽ More We present the binary model and the eclipse timing variations of the eclipsing binary RR Cae, which consists of a white dwarf eclipsed by an M-type dwarf companion. The multi-wavelength optical photometry from the Transiting Exoplanet Survey Satellite (TESS), the 0.6-m PROMPT-8 telescope, and the 0.7-m Thai Robotic Telescope at Spring Brook Observatory, combined with archive H-alpha radial velocities from the Very Large Telescope (VLT) are analysed. From the data, the physical parameters of the system are obtained along with 430 new times of minima. The TESS light curves in 2018 and 2020 show out-of-eclipse variations, which might be caused by a large spot on the secondary component. The light travel time effect models due to the gravitational interaction of one or two circumbinary objects are adopted to fit the cyclic variations in the RR Cae's O-C curve. The fitting solution of the O-C curve with one circumbinary object model shows a periodic variation with a period of $16.6\pm0.2$ yr and an amplitude of $14\pm1$ s, which can be caused by a planet with a minimum mass of $3.4\pm0.2$ M$_{Jup}$. When we consider the model with two circumbinary objects, the O-C curve shows cyclic variations with periods of $15.0\pm0.5$ yr and $39\pm5$ yr and amplitudes of $12\pm1$ s and $20\pm5$ s, respectively, corresponding to minimum masses of $3.0\pm0.3$ M$_{Jup}$ and $2.7\pm0.7$ M$_{Jup}$. △ Less

Submitted 7 June, 2023; v1 submitted 20 September, 2021; originally announced September 2021.

Comments: 16 pages, 17 figures, 8 tables

Searching For Fermi GRB Optical Counterparts With The Prototype Gravitational-Wave Optical Transient Observer (GOTO)

Authors: Y. -L. Mong, K. Ackley, D. K. Galloway, M. Dyer, R. Cutter, M. J. I. Brown, J. Lyman, K. Ulaczyk, D. Steeghs, V. Dhillon, P. OBrien, G. Ramsay, K. Noysena, R. Kotak, R. Breton, L. Nuttall, E. Palle, D. Pollacco, E. Thrane, S. Awiphan, U. Burhanudin, P. Chote, A. Chrimes, E. Daw, C. Duffy , et al. (23 additional authors not shown)

Abstract: The typical detection rate of $\sim1$ gamma-ray burst (GRB) per day by the \emph{Fermi} Gamma-ray Burst Monitor (GBM) provides a valuable opportunity to further our understanding of GRB physics. However, the large uncertainty of the \emph{Fermi} localization typically prevents rapid identification of multi-wavelength counterparts. We report the follow-up of 93 \emph{Fermi} GRBs with the Gravitatio… ▽ More The typical detection rate of $\sim1$ gamma-ray burst (GRB) per day by the \emph{Fermi} Gamma-ray Burst Monitor (GBM) provides a valuable opportunity to further our understanding of GRB physics. However, the large uncertainty of the \emph{Fermi} localization typically prevents rapid identification of multi-wavelength counterparts. We report the follow-up of 93 \emph{Fermi} GRBs with the Gravitational-wave Optical Transient Observer (GOTO) prototype on La Palma. We selected 53 events (based on favourable observing conditions) for detailed analysis, and to demonstrate our strategy of searching for optical counterparts. We apply a filtering process consisting of both automated and manual steps to 60\,085 candidates initially, rejecting all but 29, arising from 15 events. With $\approx3$ GRB afterglows expected to be detectable with GOTO from our sample, most of the candidates are unlikely to be related to the GRBs. Since we did not have multiple observations for those candidates, we cannot confidently confirm the association between the transients and the GRBs. Our results show that GOTO can effectively search for GRB optical counterparts thanks to its large field of view of $\approx40$ square degrees and its depth of $\approx20$ mag. We also detail several methods to improve our overall performance for future follow-up programs of \emph{Fermi} GRBs. △ Less

Submitted 26 August, 2021; originally announced August 2021.

Probing extreme environments with the Cherenkov Telescope Array

Authors: C. Boisson, A. M. Brown, A. Burtovoi, M. Cerruti, M. Chernyakova, T. Hassan, J. -P. Lenain, M. Manganaro, P. Romano, H. Sol, F. Tavecchio, S. Vercellone, L. Zampieri, R. Zanin, A. Zech, I. Agudo, R. Alves Batista, E. O. Anguner, L. A. Antonelli, M. Backes, C. Balazs, J. Becerra González, C. Bigongiari, E. Bissaldi, J. Bolmont , et al. (105 additional authors not shown)

Abstract: The physics of the non-thermal Universe provides information on the acceleration mechanisms in extreme environments, such as black holes and relativistic jets, neutron stars, supernovae or clusters of galaxies. In the presence of magnetic fields, particles can be accelerated towards relativistic energies. As a consequence, radiation along the entire electromagnetic spectrum can be observed, and ex… ▽ More The physics of the non-thermal Universe provides information on the acceleration mechanisms in extreme environments, such as black holes and relativistic jets, neutron stars, supernovae or clusters of galaxies. In the presence of magnetic fields, particles can be accelerated towards relativistic energies. As a consequence, radiation along the entire electromagnetic spectrum can be observed, and extreme environments are also the most likely sources of multi-messenger emission. The most energetic part of the electromagnetic spectrum corresponds to the very-high-energy (VHE, E>100 GeV) gamma-ray regime, which can be extensively studied with ground based Imaging Atmospheric Cherenkov Telescopes (IACTs). The results obtained by the current generation of IACTs, such as H.E.S.S., MAGIC, and VERITAS, demonstrate the crucial importance of the VHE band in understanding the non-thermal emission of extreme environments in our Universe. In some objects, the energy output in gamma rays can even outshine the rest of the broadband spectrum. The Cherenkov Telescope Array (CTA) is the next generation of IACTs, which, with cutting edge technology and a strategic configuration of ~100 telescopes distributed in two observing sites, in the northern and southern hemispheres, will reach better sensitivity, angular and energy resolution, and broader energy coverage than currently operational IACTs. With CTA we can probe the most extreme environments and considerably boost our knowledge of the non-thermal Universe. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted as input to ASTRONET Science Vision and Infrastructure roadmap on behalf of the CTA consortium

Multi-messenger and transient astrophysics with the Cherenkov Telescope Array

Authors: Ž. Bošnjak, A. M. Brown, A. Carosi, M. Chernyakova, P. Cristofari, F. Longo, A. López-Oramas, M. Santander, K. Satalecka, F. Schüssler, O. Sergijenko, A. Stamerra, I. Agudo, R. Alves Batista, E. Amato, E. O. Anguner, L. A. Antonelli, M. Backes, Csaba Balazs, L. Baroncelli, J. Becker Tjus, C. Bigongiari, E. Bissaldi, C. Boisson, J. Bolmont , et al. (120 additional authors not shown)

Abstract: The discovery of gravitational waves, high-energy neutrinos or the very-high-energy counterpart of gamma-ray bursts has revolutionized the high-energy and transient astrophysics community. The development of new instruments and analysis techniques will allow the discovery and/or follow-up of new transient sources. We describe the prospects for the Cherenkov Telescope Array (CTA), the next-generati… ▽ More The discovery of gravitational waves, high-energy neutrinos or the very-high-energy counterpart of gamma-ray bursts has revolutionized the high-energy and transient astrophysics community. The development of new instruments and analysis techniques will allow the discovery and/or follow-up of new transient sources. We describe the prospects for the Cherenkov Telescope Array (CTA), the next-generation ground-based gamma-ray observatory, for multi-messenger and transient astrophysics in the decade ahead. CTA will explore the most extreme environments via very-high-energy observations of compact objects, stellar collapse events, mergers and cosmic-ray accelerators. △ Less

Submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted to ASTRONET roadmap on behalf of the CTA consortium

Origin and role of relativistic cosmic particles

Authors: A. Araudo, G. Morlino, B. Olmi, F. Acero, I. Agudo, R. Adam, R. Alves Batista, E. Amato, E. O. Anguner, L. A. Antonelli, Y. Ascasibar, C. Balazs, J. Becker Tjus, C. Bigongiari, E. Bissaldi, J. Bolmont, C. Boisson, P. Bordas, Ž. Bošnjak, A. M. Brown, M. Burton, N. Bucciantini, F. Cangemi, P. Caraveo, M. Cardillo , et al. (99 additional authors not shown)

Abstract: This white paper briefly summarizes the importance of the study of relativistic cosmic rays, both as a constituent of our Universe, and through their impact on stellar and galactic evolution. The focus is on what can be learned over the coming decade through ground-based gamma-ray observations over the 20 GeV to 300 TeV range. The majority of the material is drawn directly from "Science with the C… ▽ More This white paper briefly summarizes the importance of the study of relativistic cosmic rays, both as a constituent of our Universe, and through their impact on stellar and galactic evolution. The focus is on what can be learned over the coming decade through ground-based gamma-ray observations over the 20 GeV to 300 TeV range. The majority of the material is drawn directly from "Science with the Cherenkov Telescope Array", which describes the overall science case for CTA. We request that authors wishing to cite results contained in this white paper cite the original work. △ Less

Submitted 15 June, 2021; v1 submitted 7 June, 2021; originally announced June 2021.

Comments: Submitted as input to ASTRONET Science Vision and Infrastructure roadmap on behalf of the CTA consortium. arXiv admin note: text overlap with arXiv:1709.07997

Light curve classification with recurrent neural networks for GOTO: dealing with imbalanced data

Authors: U. F. Burhanudin, J. R. Maund, T. Killestein, K. Ackley, M. J. Dyer, J. Lyman, K. Ulaczyk, R. Cutter, Y. -L. Mong, D. Steeghs, D. K. Galloway, V. Dhillon, P. O'Brien, G. Ramsay, K. Noysena, R. Kotak, R. P. Breton, L. Nuttall, E. Pallé, D. Pollacco, E. Thrane, S. Awiphan, P. Chote, A. Chrimes, E. Daw , et al. (20 additional authors not shown)

Abstract: The advent of wide-field sky surveys has led to the growth of transient and variable source discoveries. The data deluge produced by these surveys has necessitated the use of machine learning (ML) and deep learning (DL) algorithms to sift through the vast incoming data stream. A problem that arises in real-world applications of learning algorithms for classification is imbalanced data, where a cla… ▽ More The advent of wide-field sky surveys has led to the growth of transient and variable source discoveries. The data deluge produced by these surveys has necessitated the use of machine learning (ML) and deep learning (DL) algorithms to sift through the vast incoming data stream. A problem that arises in real-world applications of learning algorithms for classification is imbalanced data, where a class of objects within the data is underrepresented, leading to a bias for over-represented classes in the ML and DL classifiers. We present a recurrent neural network (RNN) classifier that takes in photometric time-series data and additional contextual information (such as distance to nearby galaxies and on-sky position) to produce real-time classification of objects observed by the Gravitational-wave Optical Transient Observer (GOTO), and use an algorithm-level approach for handling imbalance with a focal loss function. The classifier is able to achieve an Area Under the Curve (AUC) score of 0.972 when using all available photometric observations to classify variable stars, supernovae, and active galactic nuclei. The RNN architecture allows us to classify incomplete light curves, and measure how performance improves as more observations are included. We also investigate the role that contextual information plays in producing reliable object classification. △ Less

Submitted 25 May, 2021; v1 submitted 24 May, 2021; originally announced May 2021.

Comments: 16 pages, 12 figures, to be published in Monthly Notices of the Royal Astronomical Society

Processing GOTO data with the Rubin Observatory LSST Science Pipelines II: Forced Photometry and light curves

Authors: L. Makrygianni, J. Mullaney, V. Dhillon, S. Littlefair, K. Ackley, M. J. Dyer, J. Lyman, K. Ulaczyk, R. Cutter, Y. -L. Mong, D. Steeghs, D. K. Galloway, P. O'Brien, G. Ramsay, S. Poshyachinda, R. Kotak, L. Nuttall, E. Pallé, D. Pollacco, E. Thrane, S. Aukkaravittayapun, S. Awiphan, R. Breton, U. Burhanudin, P. Chote , et al. (23 additional authors not shown)

Abstract: We have adapted the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines to process data from the Gravitational-Wave Optical Transient Observer (GOTO) prototype. In this paper, we describe how we used the Rubin Observatory LSST Science Pipelines to conduct forced photometry measurements on nightly GOTO data. By comparing the photometry measurements of sources taken on… ▽ More We have adapted the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST) Science Pipelines to process data from the Gravitational-Wave Optical Transient Observer (GOTO) prototype. In this paper, we describe how we used the Rubin Observatory LSST Science Pipelines to conduct forced photometry measurements on nightly GOTO data. By comparing the photometry measurements of sources taken on multiple nights, we find that the precision of our photometry is typically better than 20~mmag for sources brighter than 16 mag. We also compare our photometry measurements against colour-corrected PanSTARRS photometry, and find that the two agree to within 10~mmag (1$σ$) for bright (i.e., $\sim14^{\rm th}$~mag) sources to 200~mmag for faint (i.e., $\sim18^{\rm th}$~mag) sources. Additionally, we compare our results to those obtained by GOTO's own in-house pipeline, {\sc gotophoto}, and obtain similar results. Based on repeatability measurements, we measure a $5σ$ L-band survey depth of between 19 and 20 magnitudes, depending on observing conditions. We assess, using repeated observations of non-varying standard SDSS stars, the accuracy of our uncertainties, which we find are typically overestimated by roughly a factor of two for bright sources (i.e., $<15^{\rm th}$~mag), but slightly underestimated (by roughly a factor of 1.25) for fainter sources ($>17^{\rm th}$~mag). Finally, we present lightcurves for a selection of variable sources, and compare them to those obtained with the Zwicky Transient Factory and GAIA. Despite the Rubin Observatory LSST Science Pipelines still undergoing active development, our results show that they are already delivering robust forced photometry measurements from GOTO data. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: Accepted for publication in PASA

Transient-optimised real-bogus classification with Bayesian Convolutional Neural Networks -- sifting the GOTO candidate stream

Authors: T. L. Killestein, J. Lyman, D. Steeghs, K. Ackley, M. J. Dyer, K. Ulaczyk, R. Cutter, Y. -L. Mong, D. K. Galloway, V. Dhillon, P. O'Brien, G. Ramsay, S. Poshyachinda, R. Kotak, R. P. Breton, L. K. Nuttall, E. Pall�, D. Pollacco, E. Thrane, S. Aukkaravittayapun, S. Awiphan, U. Burhanudin, P. Chote, A. Chrimes, E. Daw , et al. (23 additional authors not shown)

Abstract: Large-scale sky surveys have played a transformative role in our understanding of astrophysical transients, only made possible by increasingly powerful machine learning-based filtering to accurately sift through the vast quantities of incoming data generated. In this paper, we present a new real-bogus classifier based on a Bayesian convolutional neural network that provides nuanced, uncertainty-aw… ▽ More Large-scale sky surveys have played a transformative role in our understanding of astrophysical transients, only made possible by increasingly powerful machine learning-based filtering to accurately sift through the vast quantities of incoming data generated. In this paper, we present a new real-bogus classifier based on a Bayesian convolutional neural network that provides nuanced, uncertainty-aware classification of transient candidates in difference imaging, and demonstrate its application to the datastream from the GOTO wide-field optical survey. Not only are candidates assigned a well-calibrated probability of being real, but also an associated confidence that can be used to prioritise human vetting efforts and inform future model optimisation via active learning. To fully realise the potential of this architecture, we present a fully-automated training set generation method which requires no human labelling, incorporating a novel data-driven augmentation method to significantly improve the recovery of faint and nuclear transient sources. We achieve competitive classification accuracy (FPR and FNR both below 1%) compared against classifiers trained with fully human-labelled datasets, whilst being significantly quicker and less labour-intensive to build. This data-driven approach is uniquely scalable to the upcoming challenges and data needs of next-generation transient surveys. We make our data generation and model training codes available to the community. △ Less

Submitted 19 February, 2021; originally announced February 2021.

Comments: 17 pages, 12 figures, resubmitted to MNRAS following reviewer comments

Processing GOTO data with the Rubin Observatory LSST Science Pipelines I : Production of coadded frames

Authors: J. R. Mullaney, L. Makrygianni, V. Dhillon, S. Littlefair, K. Ackley, M. Dyer, J. Lyman, K. Ulaczyk, R. Cutter, Y. L. Mong, D. Steeghs, D. K. Galloway, P. O'Brien, G. Ramsay, S. Poshyachinda, R. Kotak, L. Nuttall, E. Pall�, D. Pollacco, E. Thrane, S. Aukkaravittayapun, S. Awiphan, R. Breton, U. Burhanudin, P. Chote , et al. (22 additional authors not shown)

Abstract: The past few decades have seen the burgeoning of wide field, high cadence surveys, the most formidable of which will be the Legacy Survey of Space and Time (LSST) to be conducted by the Vera C. Rubin Observatory. So new is the field of systematic time-domain survey astronomy, however, that major scientific insights will continue to be obtained using smaller, more flexible systems than the LSST. On… ▽ More The past few decades have seen the burgeoning of wide field, high cadence surveys, the most formidable of which will be the Legacy Survey of Space and Time (LSST) to be conducted by the Vera C. Rubin Observatory. So new is the field of systematic time-domain survey astronomy, however, that major scientific insights will continue to be obtained using smaller, more flexible systems than the LSST. One such example is the Gravitational-wave Optical Transient Observer (GOTO), whose primary science objective is the optical follow-up of Gravitational Wave events. The amount and rate of data production by GOTO and other wide-area, high-cadence surveys presents a significant challenge to data processing pipelines which need to operate in near real-time to fully exploit the time-domain. In this study, we adapt the Rubin Observatory LSST Science Pipelines to process GOTO data, thereby exploring the feasibility of using this "off-the-shelf" pipeline to process data from other wide-area, high-cadence surveys. In this paper, we describe how we use the LSST Science Pipelines to process raw GOTO frames to ultimately produce calibrated coadded images and photometric source catalogues. After comparing the measured astrometry and photometry to those of matched sources from PanSTARRS DR1, we find that measured source positions are typically accurate to sub-pixel levels, and that measured L-band photometries are accurate to $\sim50$ mmag at $m_L\sim16$ and $\sim200$ mmag at $m_L\sim18$. These values compare favourably to those obtained using GOTO's primary, in-house pipeline, GOTOPHOTO, in spite of both pipelines having undergone further development and improvement beyond the implementations used in this study. Finally, we release a generic "obs package" that others can build-upon should they wish to use the LSST Science Pipelines to process data from other facilities. △ Less

Submitted 28 October, 2020; originally announced October 2020.

Comments: 17 pages, 13 figures. Accepted for publication in PASA

Sensitivity of the Cherenkov Telescope Array for probing cosmology and fundamental physics with gamma-ray propagation

Authors: The Cherenkov Telescope Array Consortium, :, H. Abdalla, H. Abe, F. Acero, A. Acharyya, R. Adam, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves B, L. Amati, E. Amato, G. Ambrosi, E. O. Angüner, A. Araudo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, Y. Ascasíbar, M. Ashley , et al. (474 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nucle… ▽ More The Cherenkov Telescope Array (CTA), the new-generation ground-based observatory for $γ$-ray astronomy, provides unique capabilities to address significant open questions in astrophysics, cosmology, and fundamental physics. We study some of the salient areas of $γ$-ray cosmology that can be explored as part of the Key Science Projects of CTA, through simulated observations of active galactic nuclei (AGN) and of their relativistic jets. Observations of AGN with CTA will enable a measurement of $γ$-ray absorption on the extragalactic background light with a statistical uncertainty below 15% up to a redshift $z=2$ and to constrain or detect $γ$-ray halos up to intergalactic-magnetic-field strengths of at least 0.3pG. Extragalactic observations with CTA also show promising potential to probe physics beyond the Standard Model. The best limits on Lorentz invariance violation from $γ$-ray astronomy will be improved by a factor of at least two to three. CTA will also probe the parameter space in which axion-like particles could constitute a significant fraction, if not all, of dark matter. We conclude on the synergies between CTA and other upcoming facilities that will foster the growth of $γ$-ray cosmology. △ Less

Submitted 26 February, 2021; v1 submitted 3 October, 2020; originally announced October 2020.

Comments: 71 pages (including affiliations and references), 13 figures, 6 tables. Accepted in JCAP; matches published version. Corresponding authors: Jonathan Biteau, Julien Lefaucheur, Humberto Martinez-Huerta, Manuel Meyer, Santiago Pita, Ievgen Vovk

Journal ref: JCAP 02 (2021) 048

Machine Learning for Transient Recognition in Difference Imaging With Minimum Sampling Effort

Authors: Yik-Lun Mong, Kendall Ackley, Duncan Galloway, Tom Killestein, Joe Lyman, Danny Steeghs, Vik Dhillon, Paul O'Brien, Gavin Ramsay, Saran Poshyachinda, Rubina Kotak, Laura Nuttall, Enric Pall'e, Don Pollacco, Eric Thrane, Martin Dyer, Krzysztof Ulaczyk, Ryan Cutter, James McCormac, Paul Chote, Andrew Levan, Tom Marsh, Elizabeth Stanway, Ben Gompertz, Klaas Wiersema , et al. (23 additional authors not shown)

Abstract: The amount of observational data produced by time-domain astronomy is exponentially in-creasing. Human inspection alone is not an effective way to identify genuine transients fromthe data. An automatic real-bogus classifier is needed and machine learning techniques are commonly used to achieve this goal. Building a training set with a sufficiently large number of verified transients is challenging… ▽ More The amount of observational data produced by time-domain astronomy is exponentially in-creasing. Human inspection alone is not an effective way to identify genuine transients fromthe data. An automatic real-bogus classifier is needed and machine learning techniques are commonly used to achieve this goal. Building a training set with a sufficiently large number of verified transients is challenging, due to the requirement of human verification. We presentan approach for creating a training set by using all detections in the science images to be thesample of real detections and all detections in the difference images, which are generated by the process of difference imaging to detect transients, to be the samples of bogus detections. This strategy effectively minimizes the labour involved in the data labelling for supervised machine learning methods. We demonstrate the utility of the training set by using it to train several classifiers utilizing as the feature representation the normalized pixel values in 21-by-21pixel stamps centered at the detection position, observed with the Gravitational-wave Optical Transient Observer (GOTO) prototype. The real-bogus classifier trained with this strategy can provide up to 95% prediction accuracy on the real detections at a false alarm rate of 1%. △ Less

Submitted 13 October, 2020; v1 submitted 23 August, 2020; originally announced August 2020.

Comments: 9 pages, 8 figures

Sensitivity of the Cherenkov Telescope Array to a dark matter signal from the Galactic centre

Authors: The Cherenkov Telescope Array Consortium, :, A. Acharyya, R. Adam, C. Adams, I. Agudo, A. Aguirre-Santaella, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, L. Amati, G. Ambrosi, E. O. Ang�ner, L. A. Antonelli, C. Aramo, A. Araudo, T. Armstrong, F. Arqueros, K. Asano, Y. Ascas�bar, M. Ashley, C. Balazs, O. Ballester , et al. (427 additional authors not shown)

Abstract: We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models giv… ▽ More We provide an updated assessment of the power of the Cherenkov Telescope Array (CTA) to search for thermally produced dark matter at the TeV scale, via the associated gamma-ray signal from pair-annihilating dark matter particles in the region around the Galactic centre. We find that CTA will open a new window of discovery potential, significantly extending the range of robustly testable models given a standard cuspy profile of the dark matter density distribution. Importantly, even for a cored profile, the projected sensitivity of CTA will be sufficient to probe various well-motivated models of thermally produced dark matter at the TeV scale. This is due to CTA's unprecedented sensitivity, angular and energy resolutions, and the planned observational strategy. The survey of the inner Galaxy will cover a much larger region than corresponding previous observational campaigns with imaging atmospheric Cherenkov telescopes. CTA will map with unprecedented precision the large-scale diffuse emission in high-energy gamma rays, constituting a background for dark matter searches for which we adopt state-of-the-art models based on current data. Throughout our analysis, we use up-to-date event reconstruction Monte Carlo tools developed by the CTA consortium, and pay special attention to quantifying the level of instrumental systematic uncertainties, as well as background template systematic errors, required to probe thermally produced dark matter at these energies. "Full likelihood tables complementing our analysis are provided here [ https://doi.org/10.5281/zenodo.4057987 ]" △ Less

Submitted 30 January, 2021; v1 submitted 31 July, 2020; originally announced July 2020.

Comments: 68 pages (including references) and 26 figures; text identical to the version published in JCAP

Journal ref: JCAP01(2021)057

Searching for Electromagnetic Counterparts to Gravitational-wave Merger Events with the Prototype Gravitational-wave Optical Transient Observer (GOTO-4)

Authors: B. P. Gompertz, R. Cutter, D. Steeghs, D. K. Galloway, J. Lyman, K. Ulaczyk, M. J. Dyer, K. Ackley, V. S. Dhillon, P. T. O'Brien, G. Ramsay, S. Poshyachinda, R. Kotak, L. Nuttall, R. P. Breton, E. Pallé, D. Pollacco, E. Thrane, S. Aukkaravittayapun, S. Awiphan, M. J. I. Brown, U. Burhanudin, P. Chote, A. A. Chrimes, E. Daw , et al. (27 additional authors not shown)

Abstract: We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light… ▽ More We report the results of optical follow-up observations of 29 gravitational-wave triggers during the first half of the LIGO-Virgo Collaboration (LVC) O3 run with the Gravitational-wave Optical Transient Observer (GOTO) in its prototype 4-telescope configuration (GOTO-4). While no viable electromagnetic counterpart candidate was identified, we estimate our 3D (volumetric) coverage using test light curves of on- and off-axis gamma-ray bursts and kilonovae. In cases where the source region was observable immediately, GOTO-4 was able to respond to a GW alert in less than a minute. The average time of first observation was $8.79$ hours after receiving an alert ($9.90$ hours after trigger). A mean of $732.3$ square degrees were tiled per event, representing on average $45.3$ per cent of the LVC probability map, or $70.3$ per cent of the observable probability. This coverage will further improve as the facility scales up alongside the localisation performance of the evolving gravitational-wave detector network. Even in its 4-telescope prototype configuration, GOTO is capable of detecting AT2017gfo-like kilonovae beyond 200~Mpc in favourable observing conditions. We cannot currently place meaningful electromagnetic limits on the population of distant ($\hat{D}_L = 1.3$~Gpc) binary black hole mergers because our test models are too faint to recover at this distance. However, as GOTO is upgraded towards its full 32-telescope, 2 node (La Palma \& Australia) configuration, it is expected to be sufficiently sensitive to cover the predicted O4 binary neutron star merger volume, and will be able to respond to both northern and southern triggers. △ Less

Submitted 28 July, 2020; v1 submitted 31 March, 2020; originally announced April 2020.

Comments: 15 pages, 7 figures, 3 tables. Accepted for publication in MNRAS. Author's final submitted version

Weak Lensing Analysis of CODEX Clusters using Dark Energy Camera Legacy Survey : Mass-Richness Relation

Authors: Anirut Phriksee, Eric Jullo, Marceau Limousin, HuanYuan Shan, Alexis Finoguenov, Siramas Komonjinda, Suwicha Wannawichian, Utane Sawangwit

Abstract: We present the weak lensing analysis of 279 CODEX clusters using imaging data from 4200 $\text{deg}^{2}$ of the DECam Legacy Survey (DECaLS) Data Release 3. The cluster sample results from a joint selection in X-ray, optical richness in the range $20 \leq λ< 110$, and redshift in the range $0.1 \leq z \leq 0.2$. We model the cluster mass ($M_{\rm 200c}$) and the richness relation with the expressi… ▽ More We present the weak lensing analysis of 279 CODEX clusters using imaging data from 4200 $\text{deg}^{2}$ of the DECam Legacy Survey (DECaLS) Data Release 3. The cluster sample results from a joint selection in X-ray, optical richness in the range $20 \leq λ< 110$, and redshift in the range $0.1 \leq z \leq 0.2$. We model the cluster mass ($M_{\rm 200c}$) and the richness relation with the expression $\left\langle M_{\rm 200c} | λ\right\rangle \propto M_{0} \, (λ/ 40)^{F_λ}$. By measuring the CODEX cluster sample as an individual cluster, we obtain the best-fit values, $M_{0} = 3.24^{+0.29}_{-0.27} \times 10^{14} \text{M}_{\odot}$, and $F_λ = 1.00 ^{+0.22}_{-0.22}$ for the richness scaling index, consistent with a power law relation. Moreover, we separate the cluster sample into three richness groups; $λ= 20 - 30, 30 - 50$ and $50 - 110$, and measure the stacked excess surface mass density profile in each group. The results show that both methods are consistent. In addition, we find an excellent agreement between our weak lensing based scaling relation and the relation obtained with dynamical masses estimated from cluster member velocity dispersions measured by the SDSS-IV/SPIDERS team. This suggests that the cluster dynamical equilibrium assumption involved in the dynamical mass estimates is statistically robust for a large sample of clusters. △ Less

Submitted 24 October, 2019; originally announced October 2019.

Comments: 14 pages, 14 figures, accepted for publication in MNRAS

Signature of stochastic acceleration and cooling processes in an outburst phase of the TeV blazar ON 231

Authors: Nibedita Kalita, Utane Sawangwit, Alok C. Gupta, Paul J. Wiita

Abstract: We present a detailed spectral and temporal study of the intermediate-type blazar ON 231 during the TeV outburst phase in 2008 June with observations performed by Swift and XMM-Newton. The X-ray flux of the source, which was significantly dominated by the soft photons (below $3-4$ keV), varies between 27$\%$ and 38$\%$ on day timescales, while mild variations were observed in the optical/UV emissi… ▽ More We present a detailed spectral and temporal study of the intermediate-type blazar ON 231 during the TeV outburst phase in 2008 June with observations performed by Swift and XMM-Newton. The X-ray flux of the source, which was significantly dominated by the soft photons (below $3-4$ keV), varies between 27$\%$ and 38$\%$ on day timescales, while mild variations were observed in the optical/UV emissions. We found a maximum soft lag of $\sim 1$ hr between the UV and soft X-ray bands, which can be understood if the magnetic field of the emitting region is $\sim 5.6~ δ^{-1/3}$ G. The $0.6-10$ keV spectra can be well represented by a broken power-law model, which indicates the presence of both synchrotron and inverse Compton components in the studied X-ray regime. The synchrotron part of the SEDs constructed with simultaneous optical/UV and X-ray data follows a log-parabolic shape. A time-resolved spectral analysis shows that the break energy varies significantly between 2.4 and 7.3 keV with the changing flux state of the source, and the similar variations of the spectral slopes of the two components support the SSC scenario. The synchrotron tail, following a log-parabolic function, shows that the peak frequency ($ν_{p}$) varies by two orders of magnitude ($\sim 10^{14}-10^{16}$ Hz) during the event. A significantly positive $E_{p}-β$ relation is observed from both SED and time-resolved spectral analyses. The most feasible scenario for the observed trend during the flaring event could be associated with a magnetic-field-driven stochastic process evolving toward an equilibrium energy level. △ Less

Submitted 8 July, 2019; originally announced July 2019.

Comments: 16 pages, 13 figures, 7 tables, accepted for publication in ApJ

Swift and XMM-Newton monitoring of the ISP blazar ON 231 in outburst state

Authors: Nibedita Kalita, Utane Sawangwit

Abstract: We present a detailed temporal and spectral study of the intermediate type blazar, ON 231 with observations taken by Swift and XMM-Newton satellites during an outburst phase in June 2008. The source shows multiple X-ray flares in that period when the flux amplitude varies in the range of 27-38%. The X-ray spectra are well fitted with a broken power-law model indicating the presence of both synchro… ▽ More We present a detailed temporal and spectral study of the intermediate type blazar, ON 231 with observations taken by Swift and XMM-Newton satellites during an outburst phase in June 2008. The source shows multiple X-ray flares in that period when the flux amplitude varies in the range of 27-38%. The X-ray spectra are well fitted with a broken power-law model indicating the presence of both synchrotron and inverse Compton components. We find that the source shows strong and variable emission in the soft energy bands (below 3-4 keV), where the hard band emission is weak and stable. All the soft X-ray bands show correlated variability with zero lag, while a soft lag of -600 s between the 0.3-0.5 and 4-10 keV bands is observed with DCF analysis. A time-resolved spectral study of the flares gives a positive relation between the total fluxes and the break energies of the two emission components. △ Less

Submitted 28 May, 2019; originally announced June 2019.

Comments: Proceedings IAU XXX General Assembly, Vienna 2018; Focus Meeting FM03, Radio Galaxies: Resolving the AGN phenomenon

Journal ref: Astronomy in Focus, Proceedings of the IAU, Volume 1, 2018

Monte Carlo studies for the optimisation of the Cherenkov Telescope Array layout

Authors: A. Acharyya, I. Agudo, E. O. Angüner, R. Alfaro, J. Alfaro, C. Alispach, R. Aloisio, R. Alves Batista, J. -P. Amans, L. Amati, E. Amato, G. Ambrosi, L. A. Antonelli, C. Aramo, T. Armstrong, F. Arqueros, L. Arrabito, K. Asano, H. Ashkar, C. Balazs, M. Balbo, B. Balmaverde, P. Barai, A. Barbano, M. Barkov , et al. (445 additional authors not shown)

Abstract: The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possi… ▽ More The Cherenkov Telescope Array (CTA) is the major next-generation observatory for ground-based very-high-energy gamma-ray astronomy. It will improve the sensitivity of current ground-based instruments by a factor of five to twenty, depending on the energy, greatly improving both their angular and energy resolutions over four decades in energy (from 20 GeV to 300 TeV). This achievement will be possible by using tens of imaging Cherenkov telescopes of three successive sizes. They will be arranged into two arrays, one per hemisphere, located on the La Palma island (Spain) and in Paranal (Chile). We present here the optimised and final telescope arrays for both CTA sites, as well as their foreseen performance, resulting from the analysis of three different large-scale Monte Carlo productions. △ Less

Submitted 2 April, 2019; originally announced April 2019.

Comments: 48 pages, 16 figures, accepted for publication in Astroparticle Physics

The relationship between the radio core-dominance parameter and spectral index in different classes of extragalactic radio sources (II)

Authors: Zhi-Yuan Pei, Jun-Hui Fan, Denis Bastieri, Utane Sawangwit, Jiang-He Yang

Abstract: Active galactic nuclei (AGNs) can be divided into two major classes, namely radio loud and radio quiet AGNs. A small subset of the radio-loud AGNs is called blazars, which are believed to be unified with Fanaroff & Riley type I/II (FRI/II) radio galaxies. Following our previous work (Fan et al. 2011), we present a sample of 2400 sources with measured radio flux densities of the core and extended c… ▽ More Active galactic nuclei (AGNs) can be divided into two major classes, namely radio loud and radio quiet AGNs. A small subset of the radio-loud AGNs is called blazars, which are believed to be unified with Fanaroff & Riley type I/II (FRI/II) radio galaxies. Following our previous work (Fan et al. 2011), we present a sample of 2400 sources with measured radio flux densities of the core and extended components. The sample contains 250 BL Lacs, 520 quasars, 175 Seyferts, 1178 galaxies, 153 Fanaroff & Riley type I and type II galaxies and 104 unidentified sources. We then calculate the radio core-dominance parameters and spectral indices and study their relationship. Our analysis shows that the core-dominance parameters and spectral indices are quite different for different types of sources. We also confirm that the correlation between core-dominance parameter and spectral index exists for a large sample presented in this work. △ Less

Submitted 21 November, 2018; originally announced November 2018.

Comments: 24 pages, 10 figures

Laying the Groundwork for the Development of the Data Archive of the New Robotic Telescope

Authors: Marco C. Lam, Bovornpratch Vijarnwannaluk, Pathompong Butpan, Christopher M. Copperwheat, Andrzej S. Piascik, Utane Sawangwit, Robert J. Smith, Iain A. Steele

Abstract: The Liverpool Telescope has been in fully autonomous operation since 2004. The supporting data archive facility has largely been untouched. The data provision service has not been an issue although some modernisation of the system is desirable. This project is timely. Not only does it suit the upgrade of the current LT data archive, it is in line with the design phase of the New Robotic Telescope… ▽ More The Liverpool Telescope has been in fully autonomous operation since 2004. The supporting data archive facility has largely been untouched. The data provision service has not been an issue although some modernisation of the system is desirable. This project is timely. Not only does it suit the upgrade of the current LT data archive, it is in line with the design phase of the New Robotic Telescope which will be online in the early-2020s; and with the development of a new data archive facility for a range of telescopes at the National Astronomical Research Institute of Thailand. The Newton Fund enabled us to collaborate in designing a new versatile generic system that serves all purposes. In the end, we conclude that a single system would not meet the needs of all parties and only adopt similar front-ends while the back-ends are bespoke to our respective systems and data-flows. △ Less

Submitted 3 October, 2018; originally announced October 2018.

Comments: 11 pages, 3 tables

Journal ref: Proc. SPIE 10707, Software and Cyberinfrastructure for Astronomy V, 1070721 (6 July 2018)

The SUrvey for Pulsars and Extragalactic Radio Bursts II: New FRB discoveries and their follow-up

Authors: S. Bhandari, E. F. Keane, E. D. Barr, A. Jameson, E. Petroff, S. Johnston, M. Bailes, N. D. R. Bhat, M. Burgay, S. Burke-Spolaor, M. Caleb, R. P. Eatough, C. Flynn, J. A. Green, F. Jankowski, M. Kramer, V. Venkatraman Krishnan, V. Morello, A. Possenti, B. Stappers, C. Tiburzi, W. van Straten, I. Andreoni, T. Butterley, P. Chandra , et al. (25 additional authors not shown)

Abstract: We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time sca… ▽ More We report the discovery of four Fast Radio Bursts (FRBs) in the ongoing SUrvey for Pulsars and Extragalactic Radio Bursts (SUPERB) at the Parkes Radio Telescope: FRBs 150610, 151206, 151230 and 160102. Our real-time discoveries have enabled us to conduct extensive, rapid multi-messenger follow-up at 12 major facilities sensitive to radio, optical, X-ray, gamma-ray photons and neutrinos on time scales ranging from an hour to a few months post-burst. No counterparts to the FRBs were found and we provide upper limits on afterglow luminosities. None of the FRBs were seen to repeat. Formal fits to all FRBs show hints of scattering while their intrinsic widths are unresolved in time. FRB 151206 is at low Galactic latitude, FRB 151230 shows a sharp spectral cutoff, and FRB 160102 has the highest dispersion measure (DM = $2596.1\pm0.3$ pc cm$^{-3}$) detected to date. Three of the FRBs have high dispersion measures (DM >$1500$ pc cm$^{-3}$), favouring a scenario where the DM is dominated by contributions from the Intergalactic Medium. The slope of the Parkes FRB source counts distribution with fluences $>2$ Jyms is $α=-2.2^{+0.6}_{-1.2}$ and still consistent with a Euclidean distribution ($α=-3/2$). We also find that the all-sky rate is $1.7^{+1.5}_{-0.9}\times10^3$FRBs/($4π$ sr)/day above $\sim2$ Jyms and there is currently no strong evidence for a latitude-dependent FRB sky-rate. △ Less

Submitted 21 November, 2017; originally announced November 2017.

Comments: 21 pages, 8 figures and accepted for publication in MNRAS

A polarized fast radio burst at low Galactic latitude

Authors: E. Petroff, S. Burke-Spolaor, E. F. Keane, M. A. McLaughlin, R. Miller, I. Andreoni, M. Bailes, E. D. Barr, S. R. Bernard, S. Bhandari, N. D. R. Bhat, M. Burgay, M. Caleb, D. Champion, P. Chandra, J. Cooke, V. S. Dhillon, J. S. Farnes, L. K. Hardy, P. Jaroenjittichai, S. Johnston, M. Kasliwal, M. Kramer, S. P. Littlefair, J. P. Macquart , et al. (13 additional authors not shown)

Abstract: We report on the discovery of a new fast radio burst, FRB 150215, with the Parkes radio telescope on 2015 February 15. The burst was detected in real time with a dispersion measure (DM) of 1105.6$\pm$0.8 pc cm^{-3}, a pulse duration of 2.8^{+1.2}_{-0.5} ms, and a measured peak flux density assuming the burst was at beam center of 0.7^{+0.2}_{-0.1} Jy. The FRB originated at a Galactic longitude and… ▽ More We report on the discovery of a new fast radio burst, FRB 150215, with the Parkes radio telescope on 2015 February 15. The burst was detected in real time with a dispersion measure (DM) of 1105.6$\pm$0.8 pc cm^{-3}, a pulse duration of 2.8^{+1.2}_{-0.5} ms, and a measured peak flux density assuming the burst was at beam center of 0.7^{+0.2}_{-0.1} Jy. The FRB originated at a Galactic longitude and latitude of 24.66^{\circ}, 5.28^{\circ}, 25 degrees away from the Galactic Center. The burst was found to be 43$\pm$5% linearly polarized with a rotation measure (RM) in the range -9 < RM < 12 rad m^{-2} (95% confidence level), consistent with zero. The burst was followed-up with 11 telescopes to search for radio, optical, X-ray, gamma-ray and neutrino emission. Neither transient nor variable emission was found to be associated with the burst and no repeat pulses have been observed in 17.25 hours of observing. The sightline to the burst is close to the Galactic plane and the observed physical properties of FRB 150215 demonstrate the existence of sight lines of anomalously low RM for a given electron column density. The Galactic RM foreground may approach a null value due to magnetic field reversals along the line of sight, a decreased total electron column density from the Milky Way, or some combination of these effects. A lower Galactic DM contribution might explain why this burst was detectable whereas previous searches at low latitude have had lower detection rates than those out of the plane. △ Less

Submitted 8 May, 2017; originally announced May 2017.

Comments: 19 pages, 9 figures, accepted for publication in MNRAS

Void profile from Planck lensing potential map

Authors: Teeraparb Chantavat, Utane Sawangwit, Benjamin D. Wandelt

Abstract: We use the lensing potential map from Planck CMB lensing reconstruction analysis and the "Public Cosmic Void Catalog" to measure the stacked void lensing potential. In this profile, four parameters are needed to describe the shape of voids with different characteristic radii $R_V$. However, we have found that after reducing the background noise by subtracting the average background, there is a res… ▽ More We use the lensing potential map from Planck CMB lensing reconstruction analysis and the "Public Cosmic Void Catalog" to measure the stacked void lensing potential. In this profile, four parameters are needed to describe the shape of voids with different characteristic radii $R_V$. However, we have found that after reducing the background noise by subtracting the average background, there is a residue lensing power left in the data. The inclusion of the environment shifting parameter, $γ_V$, is necessary to get a better fit to the data with the residue lensing power. We divide the voids into two redshift bins: cmass1 ($0.45 < z < 0.5$) and cmass2 ($0.5 < z < 0.6$). Our best-fit parameters are $α= 1.989\pm0.149$, $β= 12.61\pm0.56$, $δ_c=-0.697\pm0.025$, $R_S/R_V=1.039\pm0.030$, $γ_v=(-7.034\pm0.150) \times 10^{-2}$ for the cmass1 sample with 123 voids and $α= 1.956\pm0.165$, $β= 12.91\pm0.60$, $δ_c=-0.673\pm0.027$, $R_S/R_V=1.115\pm0.032$, $γ_v=(-4.512\pm0.114) \times 10^{-2}$ for the cmass2 sample with 393 voids at 68% C.L. The addition of the environment parameter is consistent with the conjecture that the Sloan Digital Sky Survey voids reside in an underdense region. △ Less

Submitted 3 February, 2017; originally announced February 2017.

Comments: 8 pages, 5 figures, 1 table, accepted for publication in ApJ

The 2QDES Pilot : The luminosity and redshift dependence of quasar clustering

Authors: Ben Chehade, T. Shanks, J. Findlay, N. Metcalfe, U. Sawangwit, M. Irwin, E. González-Solares, S. Fine, M. J. Drinkwater, S. Croom, R. J. Jurek, D. Parkinson, R. Bielby

Abstract: We present a new redshift survey, the 2dF Quasar Dark Energy Survey pilot (2QDESp), which consists of ${\approx}10000$ quasars from ${\approx}150$ deg$^2$ of the southern sky, based on VST-ATLAS imaging and 2dF/AAOmega spectroscopy. Combining our optical photometry with the WISE (W1,W2) bands we can select essentially contamination free quasar samples with $0.8{<}z{<}2.5$ and $g{<}20.5$. At fainte… ▽ More We present a new redshift survey, the 2dF Quasar Dark Energy Survey pilot (2QDESp), which consists of ${\approx}10000$ quasars from ${\approx}150$ deg$^2$ of the southern sky, based on VST-ATLAS imaging and 2dF/AAOmega spectroscopy. Combining our optical photometry with the WISE (W1,W2) bands we can select essentially contamination free quasar samples with $0.8{<}z{<}2.5$ and $g{<}20.5$. At fainter magnitudes, optical UVX selection is still required to reach our $g{\approx}22.5$ limit. Using both these techniques we observed quasar redshifts at sky densities up to $90$ deg$^{-2}$. By comparing 2QDESp with other surveys (SDSS, 2QZ and 2SLAQ) we find that quasar clustering is approximately luminosity independent, with results for all four surveys consistent with a correlation scale of $r_{0}{=}6.1{\pm}0.1 \: h^{-1}$Mpc, despite their decade range in luminosity. We find a significant redshift dependence of clustering, particularly when BOSS data with $r_{0}{=}7.3{\pm}0.1 \: h^{-1}$Mpc are included at $z{\approx}2.4$. All quasars remain consistent with having a single host halo mass of ${\approx}2{\pm}1{\times}10^{12} \: h^{-1}M_\odot$. This result implies that either quasars do not radiate at a fixed fraction of the Eddington luminosity or AGN black hole and dark matter halo masses are weakly correlated. No significant evidence is found to support fainter, X-ray selected quasars at low redshift having larger halo masses as predicted by the `hot halo' mode AGN model of Fanidakis et al. 2013. Finally, although the combined quasar sample reaches an effective volume as large as that of the original SDSS LRG sample, we do not detect the BAO feature in these data. △ Less

Submitted 15 March, 2016; originally announced March 2016.

Comments: 20 pages, 15 figures, accepted by MNRAS

Cosmological parameter constraints from CMB lensing with cosmic voids

Authors: Teeraparb Chantavat, Utane Sawangwit, P. M. Sutter, Benjamin D. Wandelt

Abstract: We investigate the potential of using cosmic voids as a probe to constrain cosmological parameters through the gravitational lensing effect of the cosmic microwave background (CMB) and make predictions for the next generation surveys. By assuming the detection of a series of $\approx 5 - 10$ voids along a line of sight within a square-degree patch of the sky, we found that they can be used to brea… ▽ More We investigate the potential of using cosmic voids as a probe to constrain cosmological parameters through the gravitational lensing effect of the cosmic microwave background (CMB) and make predictions for the next generation surveys. By assuming the detection of a series of $\approx 5 - 10$ voids along a line of sight within a square-degree patch of the sky, we found that they can be used to break the degeneracy direction of some of the cosmological parameter constraints (for example $ω_b$ and $Ω_Λ$) in comparison with the constraints from random CMB skies with the same size area for a survey with extensive integration time. This analysis is based on our current knowledge of the average void profile and analytical estimates of the void number function. We also provide combined cosmological parameter constraints between a sky patch where series of voids are detected and a patch without voids (a randomly selected patch). The full potential of this technique relies on an accurate determination of the void profile to $\approx 10$% level. For a small-area CMB observation with extensive integration time and a high signal-to-noise ratio, CMB lensing with such series of voids will provide a complementary route to cosmological parameter constraints to the CMB observations. Example of parameter constraints with a series of five voids on a $1.0^{\circ} \times 1.0^{\circ}$ patch of the sky are $100ω_b = 2.20 \pm 0.27$, $ω_c = 0.120 \pm 0.022$, $Ω_Λ= 0.682 \pm 0.078$, $Δ_{\mathcal{R}}^2 = \left(2.22 \pm 7.79\right) \times 10^{-9}$, $n_s = 0.962 \pm 0.097$ and $τ= 0.925 \pm 1.747$ at 68% C.L. △ Less

Submitted 18 February, 2016; v1 submitted 11 September, 2014; originally announced September 2014.

Comments: 8 pages, 5 figures, accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 93, 043523, 2016

Clustering analysis of high-redshift Luminous Red Galaxies in Stripe 82

Authors: Nikolaos Nikoloudakis, Tom Shanks, Utane Sawangwit

Abstract: We present a clustering analysis of Luminous Red Galaxies in SDSS Stripe 82. We study the angular 2-point correlation function, of 130,000 LRG candidates via colour-cut selections in izK with the K band coverage coming from UKIDSS LAS. We have used the cross-correlation technique of Newman (2008) to establish the LRG redshift distribution. Cross-correlating with SDSS QSOs, MegaZ-LRGs and DEEP2 gal… ▽ More We present a clustering analysis of Luminous Red Galaxies in SDSS Stripe 82. We study the angular 2-point correlation function, of 130,000 LRG candidates via colour-cut selections in izK with the K band coverage coming from UKIDSS LAS. We have used the cross-correlation technique of Newman (2008) to establish the LRG redshift distribution. Cross-correlating with SDSS QSOs, MegaZ-LRGs and DEEP2 galaxies implies an average LRG redshift of z~1 with space density, n_g~3.2 +/-0.16 x10^-4 h^3 Mpc^-3. For theta<10', w(theta) significantly deviates from a single power-law. A double power-law with a break at r_b~2.4 h^-1 Mpc fits the data better, with best-fit scale length, r_0,1=7.63+/-0.27 h^-1Mpc and slope gamma_1=2.01 +/-0.02 at small scales and r_0,2=9.92 +/-0.40 h^-1 Mpc and gamma_2=1.64 +/-0.04 at large scales. Due to the flat slope at large scales, we find that a standard LambdaCDM linear model is accepted only at 2-3sigma, with the best-fit bias factor, b=2.74+/-0.07. We also fitted HOD models and estimate an effective halo mass of M_eff=3.3 +/-0.6x10^13 h^-1 M_sun. But at large scales, the current HOD models did not help explain the power excess in the clustering signal. We then compare the w(theta) results to those of Sawangwit et al. (2011) from 3 samples of LRGs at lower redshifts to measure clustering evolution. We find that a long-lived model may be a poorer fit than at lower redshifts, although this assumes that the Stripe 82 LRGs are luminosity-matched to the AAOmega LRGs. We find stronger evidence for evolution in the form of the z~1 LRG correlation function, with the above flat 2-halo slope maintaining to r>50 h^-1 Mpc. Applying the cross-correlation test of Ross et al. (2011), we find little evidence that the result is due to systematics. Otherwise it may provide evidence for primordial non-Gaussianity in the matter distribution, with f^local_NL=90+/-30.[Abridged] △ Less

Submitted 22 November, 2012; v1 submitted 16 April, 2012; originally announced April 2012.

Comments: 20 pages, 25 figures, accepted for publication in MNRAS, more tests for systematics errors checked (Fig. 22-23), conclusions remained the same

Measuring BAO and non-Gaussianity via QSO clustering

Authors: U. Sawangwit, T. Shanks, S. M. Croom, M. J. Drinkwater, S. Fine, D. Parkinson, Nicholas P. Ross

Abstract: Our goals are (i) to search for BAO and large-scale structure in current QSO survey data and (ii) to use these and simulation/forecast results to assess the science case for a new, >10x larger, QSO survey. We first combine the SDSS, 2QZ and 2SLAQ surveys to form a survey of ~60000 QSOs. We find a hint of a peak in the QSO 2-point correlation function, xi(s), at the same scale (~105h^-1 Mpc) as det… ▽ More Our goals are (i) to search for BAO and large-scale structure in current QSO survey data and (ii) to use these and simulation/forecast results to assess the science case for a new, >10x larger, QSO survey. We first combine the SDSS, 2QZ and 2SLAQ surveys to form a survey of ~60000 QSOs. We find a hint of a peak in the QSO 2-point correlation function, xi(s), at the same scale (~105h^-1 Mpc) as detected by Eisenstein et al (2005) in their sample of DR5 LRGs but only at low statistical significance. We then compare these data with QSO mock catalogues from the Hubble Volume simulation used by Hoyle et al (2002) and find that both routes give statistical error estimates that are consistent at ~100h^-1 Mpc scales. Mock catalogues are then used to estimate the nominal survey size needed for a 3-4 sigma detection of the BAO peak. We find that a redshift survey of ~250000 z<2.2 QSOs is required over ~3000 deg^2. This is further confirmed by static log-normal simulations where the BAO are clearly detectable in the QSO power spectrum and correlation function. The nominal survey would on its own produce the first detection of, for example, discontinuous dark energy evolution in the so far uncharted 1<z<2.2 redshift range. A survey with ~50% higher QSO sky densities and 50% bigger area will give an ~6sigma BAO detection, leading to an error ~60% of the size of the BOSS error on the dark energy evolution parameter, w_a. Another important aim for a QSO survey is to place new limits on primordial non-Gaussianity at large scales, testing tentative evidence we have found for the evolution of the linear form of the combined QSO xi(s) at z~1.6. Such a QSO survey will also determine the gravitational growth rate at z~1.6 via z-space distortions, allow lensing tomography via QSO magnification bias while also measuring the exact luminosity dependence of small-scale QSO clustering. △ Less

Submitted 23 September, 2011; v1 submitted 4 August, 2011; originally announced August 2011.

Comments: 11 pages, 9 figures, typo corrected, references added, accepted by MNRAS

Evolution in the clustering strength of radio galaxies

Authors: S. Fine, T. Shanks, N. Nikoloudakis, U. Sawangwit

Abstract: We cross match the NVSS and FIRST surveys with three large photometric catalogues of luminous red galaxies (LRGs) to define radio-loud samples. These have median redshifts 0.35, 0.55 and 0.68 and, by matching rest-frame optical and radio properties, we construct uniform samples across the three surveys. This paper is concerned with the clustering properties of these samples derived from the angula… ▽ More We cross match the NVSS and FIRST surveys with three large photometric catalogues of luminous red galaxies (LRGs) to define radio-loud samples. These have median redshifts 0.35, 0.55 and 0.68 and, by matching rest-frame optical and radio properties, we construct uniform samples across the three surveys. This paper is concerned with the clustering properties of these samples derived from the angular correlation function. The primary aim is to characterise any evolution in the clustering amplitude of radio galaxies bellow z~0.68. We find no evidence for evolution in the large-scale (~1-50h^{-1}Mpc) clustering amplitude. Our radio galaxy autocorrelations are consistent with previous findings indicating little-to-no evolution in the redshift range 0.68 to 0 (~6Gyr of time). We also cross correlate radio galaxies with the parent LRG samples to increase the precision of our results and again find no evidence for evolution. Our results are inconsistent with a long-lived model for the clustering evolution that assumes radio sources randomly sample the LRG population. A model where the halo mass is constant with redshift is consistent with the data. This is similar to QSOs that have clustering amplitudes consistent with a single halo mass at all redshifts. Given that the brightest radio sources show stronger evolution in space density compared to fainter radio sources we restrict our samples to include only objects with L>10^{26}W/Hz and repeat the analysis. Again we find no evidence for evolution in the comoving correlation amplitude. These radio sources appear to inhabit the same mass halos as fainter radio galaxies (~9x10^{13}h^{-1}M_\odot). These halos are ~twice as massive as those of the general LRG population and ~30 times as massive as optical AGN/QSOs. △ Less

Submitted 28 July, 2011; originally announced July 2011.

Comments: 9 pages, 8 figures, accepted MNRAS

Testing the WMAP cosmology via Planck radio catalogues

Authors: J. R. Whitbourn, T. Shanks, U. Sawangwit

Abstract: The prime evidence underpinning the standard LCDM cosmological model is the CMB power spectrum as observed by WMAP and other microwave experiments. But Sawangwit & Shanks (2010) have recently shown that the WMAP CMB power spectrum is highly sensitive to the beam profile of the WMAP telescope. Here, we use the source catalogue from the Planck Early Data Release to test further the WMAP beam profile… ▽ More The prime evidence underpinning the standard LCDM cosmological model is the CMB power spectrum as observed by WMAP and other microwave experiments. But Sawangwit & Shanks (2010) have recently shown that the WMAP CMB power spectrum is highly sensitive to the beam profile of the WMAP telescope. Here, we use the source catalogue from the Planck Early Data Release to test further the WMAP beam profiles. We confirm that stacked beam profiles at Q, V and particularly at W appear wider than expected when compared to the Jupiter beam, normalised either directly to the radio source profiles or using Planck fluxes. The same result is also found based on WMAP-CMBfree source catalogues and NVSS sources. The accuracy of our beam profile measurements is supported by analysis of CMB sky simulations. However the beam profiles from WMAP7 at the W band are narrower than previously found in WMAP5 data and the rejection of the WMAP beam is now only at the ~3sigma level. We also find that the WMAP source fluxes demonstrate possible non-linearity with Planck fluxes. But including ground-based and Planck data for the bright Weiland et al (2011) sources may suggest that the discrepancy is a linear offset rather than a non-linearity. Additionally, we find that the stacked Sunyaev-Zel'dovich (SZ) decrements of ~151 galaxy clusters observed by Planck are in agreement with the WMAP data. We find that there is no evidence for a WMAP SZ deficit as has previously been reported. In the particular case of Coma we find evidence for the presence of an ~0.1mK downwards CMB fluctuation. We conclude that beam profile systematics can have significant effects on both the amplitude and position of the acoustic peaks, with potentially important implications for cosmology parameter fitting. △ Less

Submitted 11 February, 2013; v1 submitted 13 July, 2011; originally announced July 2011.

Comments: v3; 21 pages, 15 figures, amended in line with referee comments, resubmitted to MNRAS

Do all QSOs have the same black hole mass?

Authors: T. Shanks, S. M. Croom, S. Fine, N. P. Ross, U. Sawangwit

Abstract: QSOs from SDSS, 2QZ and 2SLAQ covering an order of magnitude in luminosity at fixed redshift exhibit similar amplitudes of clustering. In addition, QSO clustering evolution at z>0.5 is well fitted by a model that assumes a fixed host halo mass, implying that QSOs may occur in a relatively narrow range of halo and BH mass. We argue that the slow evolution of early-type galaxies out to z~1-2 may als… ▽ More QSOs from SDSS, 2QZ and 2SLAQ covering an order of magnitude in luminosity at fixed redshift exhibit similar amplitudes of clustering. In addition, QSO clustering evolution at z>0.5 is well fitted by a model that assumes a fixed host halo mass, implying that QSOs may occur in a relatively narrow range of halo and BH mass. We argue that the slow evolution of early-type galaxies out to z~1-2 may also provide support for a slow evolution of QSO host BH masses. The result would mean that if high-z QSOs radiate at Eddington rates then low-z SyI must radiate at ~100x less than Eddington. We conclude that models where QSOs radiate at L_Edd require M_BH and M_halo to be decoupled to circumvent the clustering results. While single BH mass and flickering models fit the z>0.5 clustering results, they appear to be rejected by the z~0, M_BH-L relation from reverberation mapping. We find that the inclusion of z<0.5 QSO clustering data improves the fit of a long-lived QSO model and suggest that the predictions of a PLE model for QSO BH masses agree reasonably with UV-bump and reverberation estimates (abridged). △ Less

Submitted 29 May, 2011; v1 submitted 12 May, 2011; originally announced May 2011.

Comments: 10 pages, 5 figures, accepted by MNRAS, typos corrected, references added

Lambda-CDM and the WMAP power spectrum beam profile sensitivity

Authors: Utane Sawangwit, Tom Shanks

Abstract: We first discuss the sensitivity of the WMAP CMB power spectrum to systematic errors by calculating the raw CMB power spectrum from WMAP data. We find that the power spectrum is surprisingly sensitive to the WMAP radiometer beam profile even at the position of the first acoustic peak on ~1 degree scales. Although the WMAP beam profile core is only 12.6arcmin FWHM at W, there is a long power-law ta… ▽ More We first discuss the sensitivity of the WMAP CMB power spectrum to systematic errors by calculating the raw CMB power spectrum from WMAP data. We find that the power spectrum is surprisingly sensitive to the WMAP radiometer beam profile even at the position of the first acoustic peak on ~1 degree scales. Although the WMAP beam profile core is only 12.6arcmin FWHM at W, there is a long power-law tail to the beam due to side-lobes and this causes significant effects even at the first peak position. We then test the form of the beam-profile used by the WMAP team which is based on observations of Jupiter. We stacked radio source beam profiles as observed in each WMAP band and found that they showed a wider profile in Q, V, W than the Jupiter profile. We have now checked that this is not due to any Eddington or other bias in our sample by showing that the same results are obtained when radio sources are selected at 1.4GHz and that our methods retrieve the Jupiter beam when it is employed in simulations. Finally, we show that the uncertainty in the WMAP beam profile allows the position as well as the amplitude of the first peak to be changed and how this could allow simpler cosmologies than standard Lambda-CDM to fit the CMB data. △ Less

Submitted 7 June, 2010; originally announced June 2010.

Comments: 5 pages, 4 figures, Proceedings of the conference "45th Rencontres de Moriond" - Cosmology Session, La Thuile, Val d'Aosta, Italy, March 13 - 20, 2010

Beam profile sensitivity of the WMAP CMB power spectrum

Authors: U. Sawangwit, T. Shanks

Abstract: Using the published WMAP 5-year data, we first show how sensitive the WMAP power spectra are to the form of the WMAP beam. It is well known that the beam profile derived from observations of Jupiter is non-Gaussian and indeed extends, in the W band for example, well beyond its 12.'6 FWHM core out to more than 1 degree in radius. This means that even though the core width corresponds to wavenumber… ▽ More Using the published WMAP 5-year data, we first show how sensitive the WMAP power spectra are to the form of the WMAP beam. It is well known that the beam profile derived from observations of Jupiter is non-Gaussian and indeed extends, in the W band for example, well beyond its 12.'6 FWHM core out to more than 1 degree in radius. This means that even though the core width corresponds to wavenumber l\approx1800, the form of the beam still significantly affects the WMAP results even at l\approx200 which is the scale of the first acoustic peak. The difference between the beam convolved C_l and the final C_l is \approx70% at the scale of the first peak, rising to \approx400% at the scale of the second. New estimates of the Q, V and W-band beam profiles are then presented, based on a stacking analysis of the WMAP5 radio source catalogue and temperature maps. The radio sources show a significantly (3-4σ) broader beam profile on scales of 10'-30' than that found by the WMAP team whose beam analysis is based on measurements of Jupiter. Beyond these scales the beam profiles from the radio sources are too noisy to give useful information. Furthermore, we find tentative evidence for a non-linear relation between WMAP and ATCA/IRAM 95 GHz source fluxes. We discuss whether the wide beam profiles could be caused either by radio source extension or clustering and find that neither explanation is likely. We also argue against the possibility that Eddington bias is affecting our results. The reasons for the difference between the radio source and the Jupiter beam profiles are therefore still unclear. If the radio source profiles were then used to define the WMAP beam, there could be a significant change in the amplitude and position of even the first acoustic peak. It is therefore important to identify the reasons for the differences between these two beam profile estimates. △ Less

Submitted 8 June, 2010; v1 submitted 2 December, 2009; originally announced December 2009.

Comments: Accepted for publication in MNRAS Letters (5 pages, 4 figures, 1 table) New version includes stacking analysis checks against simulations plus point source detection and cut-sky corrected C(l)'s

Angular correlation function of 1.5 million LRGs: clustering evolution and a search for BAO

Authors: U. Sawangwit, T. Shanks, F. B. Abdalla, R. D. Cannon, S. M. Croom, A. C. Edge, Nicholas P. Ross, D. A. Wake

Abstract: We present the angular correlation function measured from photometric samples comprising 1,562,800 luminous red galaxies (LRGs). Three LRG samples were extracted from the Sloan Digital Sky Survey (SDSS) imaging data, based on colour-cut selections at redshifts, z ~ 0.35, 0.55 and 0.7 as calibrated by the spectroscopic surveys, SDSS-LRG, 2dF-SDSS LRG and QSO (2SLAQ), and the AAOmega LRG survey. The… ▽ More We present the angular correlation function measured from photometric samples comprising 1,562,800 luminous red galaxies (LRGs). Three LRG samples were extracted from the Sloan Digital Sky Survey (SDSS) imaging data, based on colour-cut selections at redshifts, z ~ 0.35, 0.55 and 0.7 as calibrated by the spectroscopic surveys, SDSS-LRG, 2dF-SDSS LRG and QSO (2SLAQ), and the AAOmega LRG survey. The galaxy samples cover ~7600 s.q.deg of the sky, probing a total cosmic volume of ~5.5 h^-3} Gpc^3. The small and intermediate scale correlation functions generally show significant deviations from a single power-law fit with a well-detected break at ~1 h^-1 Mpc, consistent with the transition scale between the 1- and 2-halo terms in halo occupation models. For galaxy separations 1-20 h^-1 Mpc and at fixed luminosity, we see virtually no evolution of the clustering with redshift and the data is consistent with a simple high peaks biasing model where the comoving LRG space density is constant with z. At fixed z, the LRG clustering amplitude increases with luminosity in accordance with the simple high peaks model, with a typical LRG dark matter halo mass 10^13-10^14 h^-1 M_sol. For r < 1 h^-1 Mpc, the evolution is slightly faster and the clustering decreases towards high redshift consistent with a virialised clustering model. However, assuming the HOD and LCDM halo merger frameworks, ~2-3 per cent per Gyr of the LRGs is required to merge in order to explain the small scales clustering evolution, consistent with previous results [abridge]. △ Less

Submitted 6 July, 2011; v1 submitted 2 December, 2009; originally announced December 2009.

Comments: Replaced to match the version accepted by MNRAS, Halo model analysis and additional systematic tests are included

Journal ref: Mon.Not.Roy.Astron.Soc.416:3033-3056,2011

Cross-correlating WMAP5 with 1.5 million LRGs: a new test for the ISW effect

Authors: U. Sawangwit, T. Shanks, R. D. Cannon, S. M. Croom, Nicholas P. Ross, D. A. Wake

Abstract: We present the cross-correlation of the density map of LRGs and the temperature fluctuation in the CMB as measured by the WMAP5 observations. The LRG samples were extracted from imaging data of the SDSS based on two previous spectroscopic redshift surveys, the SDSS-LRG and the 2SLAQ surveys at average redshifts z~0.35 and z~0.55. In addition we have added a higher-redshift photometric LRG sample… ▽ More We present the cross-correlation of the density map of LRGs and the temperature fluctuation in the CMB as measured by the WMAP5 observations. The LRG samples were extracted from imaging data of the SDSS based on two previous spectroscopic redshift surveys, the SDSS-LRG and the 2SLAQ surveys at average redshifts z~0.35 and z~0.55. In addition we have added a higher-redshift photometric LRG sample based on the selection of the AAOmega LRG redshift survey at z~0.7. The total LRG sample thus comprises 1.5 million galaxies, sampling a redshift range of 0.2 < z < 0.9 over ~7600 square degrees of the sky, probing a total cosmic volume of ~5.5 h^{-3} Gpc^3. We find that the new LRG sample at z~0.7 shows very little positive evidence for the ISW effect. Indeed, the cross-correlation is negative out to ~1 deg. The standard LCDM model is rejected at ~2-3% significance by the new LRG data. We then performed a new test on the robustness of the LRG ISW detections at z~0.35 and 0.55. We made 8 rotations through 360deg of the CMB maps with respect to the LRG samples around the galactic pole. We find that in both cases there are stronger effects at angles other than zero. This implies that the z~0.35 and 0.55 ISW detections may still be subject to systematic errors which combined with the known sizeable statistical errors may leave these ISW detections looking unreliable. We have further made the rotation test on several other samples where ISW detections have been claimed and find that they also show peaks when rotated. We conclude that in the samples we have tested the ISW effect may be absent and we argue that this result may not be in contradiction with previous results. △ Less

Submitted 2 December, 2009; v1 submitted 6 November, 2009; originally announced November 2009.

Comments: Accepted for publication in MNRAS (19 pages, 14 figures, 2 tables)

Journal ref: Mon.Not.Roy.Astron.Soc.402:2228,2010

Photometric Selection of Emission Line Galaxies, Clustering Analysis and a Search for the ISW effect

Authors: Rich Bielby, T. Shanks, U. Sawangwit, S. M. Croom, Nicholas P. Ross, D. A. Wake

Abstract: We investigate the use of simple colour cuts applied to the SDSS optical imaging to perform photometric selections of emission line galaxies out to z<1. From colour-cuts using the SDSS g, r and i bands, we obtain mean photometric redshifts of z=0.32+-0.08, z=0.44+-0.12 and z=0.65+-0.21. We further calibrate our high redshift selection using spectroscopic observations with the AAOmega spectrogr… ▽ More We investigate the use of simple colour cuts applied to the SDSS optical imaging to perform photometric selections of emission line galaxies out to z<1. From colour-cuts using the SDSS g, r and i bands, we obtain mean photometric redshifts of z=0.32+-0.08, z=0.44+-0.12 and z=0.65+-0.21. We further calibrate our high redshift selection using spectroscopic observations with the AAOmega spectrograph on the 4m Anglo-Australian Telescope (AAT), observing ~50-200 galaxy candidates in 4 separate fields. With just 1-hour of integration time and with seeing of ~1.6", we successfully determined redshifts for ~65% of the targeted candidates. We calculate the angular correlation functions of the samples and find correlation lengths of r0=2.64 h-1 Mpc, r0=3.62 h-1 Mpc and r0=5.88 h-1 Mpc for the low, mid and high redshift samples respectively. Comparing these results with predicted dark matter clustering, we estimate the bias parameter for each sample to be b=0.70, b=0.92 and b=1.46. We calculate the 2-point redshift-space correlation function at z~0.6 and find a clustering amplitude of s0=6.4 h-1 Mpc. Finally, we use our photometric sample to search for the Integrated Sachs-Wolfe signal in the WMAP 5yr data. We cross-correlate our three redshift samples with the WMAP W, V, Q and K bands and find an overall trend for a positive signal similar to that expected from models. However, the signal in each is relatively weak. Combining all three galaxy samples we find a signal of wTg(<100')=0.20+-0.12 microK in the WMAP W-band, a significance of 1.7sigma. △ Less

Submitted 6 February, 2009; originally announced February 2009.

Comments: 14 pages, 17 figures, submitted to MNRAS

QSO-LRG 2-Point Cross-Correlation Function and Redshift-Space Distorions

Authors: G. Mountrichas, T. Shanks, S. M. Croom, U. Sawangwit, D. P. Schneider, A. D. Myers, K. Pimbblet

Abstract: We have measured the bias of QSOs as a function of QSO luminosity at fixed redshift (z<1) by cross-correlating them with LRGs in the same spatial volume, hence breaking the degeneracy between QSO luminosity and redshift. We use three QSO samples from 2SLAQ, 2QZ and SDSS covering a QSO absolute magnitude range, -24.5<M_{b_J}<-21.5, and cross-correlate them with 2SLAQ (z~0.5) and AAOmega (z~0.7) p… ▽ More We have measured the bias of QSOs as a function of QSO luminosity at fixed redshift (z<1) by cross-correlating them with LRGs in the same spatial volume, hence breaking the degeneracy between QSO luminosity and redshift. We use three QSO samples from 2SLAQ, 2QZ and SDSS covering a QSO absolute magnitude range, -24.5<M_{b_J}<-21.5, and cross-correlate them with 2SLAQ (z~0.5) and AAOmega (z~0.7) photometric and spectroscopic LRGs in the same redshift ranges. The 2-D and 3-D cross-clustering measurements are generally in good agreement. Our (2SLAQ) QSO-LRG clustering amplitude (r_0=6.8_{-0.3}^{+0.1}h^{-1}Mpc) as measured from the semi-projected cross-correlation function appears similar to the (2SLAQ) LRG-LRG auto-correlation amplitude (r_0=7.45\pm0.35h^{-1}Mpc) and both are higher than the (2QZ+2SLAQ) QSO-QSO amplitude (r_0\simeq5.0h^{-1}Mpc). Our measurements show remarkably little QSO-LRG cross-clustering dependence on QSO luminosity. If anything, the results imply that brighter QSOs may be less highly biased than faint QSOs, the opposite direction expected from simple high peaks biasing models. Assuming a standard LCDM model and values for b_{LRG} measured from LRG autocorrelation analyses, we find b_Q=1.45\pm0.11 at M_{b_J}\approx-24 and b_Q=1.90\pm0.16 at M_{b_J}~-22. We also find consistent results for the QSO bias from a z-space distortion analysis of the QSO-LRG cross-clustering at z~0.55. The dynamical infall results give β_Q=0.55\pm0.10, implying b_Q=1.4\pm0.2. Thus both the z-space distortion and the amplitude analyses yield b_Q~1.5 at M_{b_J}~-23. The implied DM halo mass inhabited by QSOs at z~0.55 is \sim10^{13}h^{-1}M_{\sun}, again approximately independent of QSO luminosity. △ Less

Submitted 11 January, 2008; originally announced January 2008.

Comments: 22 pages, 25 figures, 8 tables