Showing 1–50 of 267 results for author: Kong, H

The Construction of Large-scale Structure Catalogs for the Dark Energy Spectroscopic Instrument

Authors: A. J. Ross, J. Aguilar, S. Ahlen, S. Alam, A. Anand, S. Bailey, D. Bianchi, S. Brieden, D. Brooks, E. Burtin, A. Carnero Rosell, E. Chaussidon, T. Claybaugh, S. Cole, K. Dawson, A. de la Macorra, A. de Mattia, Arjun Dey, Biprateep Dey, P. Doel, K. Fanning, S. Ferraro, J. Ereza, A. Font-Ribera, J. E. Forero-Romero , et al. (59 additional authors not shown)

Abstract: We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produ… ▽ More We present the technical details on how large-scale structure (LSS) catalogs are constructed from redshifts measured from spectra observed by the Dark Energy Spectroscopic Instrument (DESI). The LSS catalogs provide the information needed to determine the relative number density of DESI tracers as a function of redshift and celestial coordinates and, e.g., determine clustering statistics. We produce catalogs that are weighted subsamples of the observed data, each matched to a weighted `random' catalog that forms an unclustered sampling of the probability density that DESI could have observed those data at each location. Precise knowledge of the DESI observing history and associated hardware performance allows for a determination of the DESI footprint and the number of times DESI has covered it at sub-arcsecond level precision. This enables the completeness of any DESI sample to be modeled at this same resolution. The pipeline developed to create LSS catalogs has been designed to easily allow robustness tests and enable future improvements. We describe how it allows ongoing work improving the match between galaxy and random catalogs, such as including further information when assigning redshifts to randoms, accounting for fluctuations in target density, accounting for variation in the redshift success rate, and accommodating blinding schemes. △ Less

Submitted 26 May, 2024; originally announced May 2024.

Comments: Supporting publication of DESI 2024II: Sample definitions, characteristics, and two-point clustering statistics

Forward modeling fluctuations in the DESI LRGs target sample using image simulations

Authors: Hui Kong, Ashley J. Ross, Klaus Honscheid, Dustin Lang, Anna Porredon, Arnaud de Mattia, Mehdi Rezaie, Rongpu Zhou, Edward Schlafly, John Moustakas, Alberto Rosado-Marin, Jessica Nicole Aguilar, Steven Ahlen, David Brooks, Edmond Chaussidon, Todd Claybaugh, Shaun Cole, Axel de la Macorra, Arjun Dey, Biprateep Dey, Peter Doel, Kevin Fanning, Jaime E. Forero-Romero, Enrique Gaztanaga, Satya Gontcho A Gontcho , et al. (28 additional authors not shown)

Abstract: We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). We update the Obiwan pipeline, which had previously been developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition makes it possible to simulate the DES… ▽ More We use the forward modeling pipeline, Obiwan, to study the imaging systematics of the Luminous Red Galaxies (LRGs) targeted by the Dark Energy Spectroscopic Instrument (DESI). We update the Obiwan pipeline, which had previously been developed to simulate the optical images used to target DESI data, to further simulate WISE images in the infrared. This addition makes it possible to simulate the DESI LRGs sample, which utilizes WISE data in the target selection. Deep DESI imaging data combined with a method to account for biases in their shapes is used to define a truth sample of potential LRG targets. We simulate a total of 15 million galaxies to obtain a simulated LRG sample (Obiwan LRGs) that predicts the variations in target density due to imaging properties. We find that the simulations predict the trends with depth observed in the data, including how they depend on the intrinsic brightness of the galaxies. We observe that faint LRGs are the main contributing power of the imaging systematics trend induced by depth. We also find significant trends in the data against Galactic extinction that are not predicted by Obiwan. These trends depend strongly on the particular map of Galactic extinction chosen to test against, implying Large-Scale Structure systematic contamination (e.g. Cosmic-Infrared Background) in the Galactic extinction maps is a likely root cause. We additionally observe that the DESI LRGs sample exhibits a complex dependency on a combination of seeing, depth, and intrinsic galaxy brightness, which is not replicated by Obiwan, suggesting discrepancies between the current simulation settings and the actual observations. The detailed findings we present should be used to guide any observational systematics mitigation treatment for the clustering of the DESI LRG sample. △ Less

Submitted 25 May, 2024; originally announced May 2024.

Comments: 46 pages, 26 figures

Autoregressive Search of Gravitational Waves: Denoising

Authors: Sangin Kim, C. Y. Hui, Jianqi Yan, Alex P. Leung, Kwangmin Oh, A. K. H. Kong, L. C. -C. Lin, Kwan-Lok Li

Abstract: Because of the small strain amplitudes of gravitational-wave (GW) signals, unveiling them in the presence of detector/environmental noise is challenging. For visualizing the signals and extracting its waveform for a comparison with theoretical prediction, a frequency-domain whitening process is commonly adopted for filtering the data. In this work, we propose an alternative template-free framework… ▽ More Because of the small strain amplitudes of gravitational-wave (GW) signals, unveiling them in the presence of detector/environmental noise is challenging. For visualizing the signals and extracting its waveform for a comparison with theoretical prediction, a frequency-domain whitening process is commonly adopted for filtering the data. In this work, we propose an alternative template-free framework based on autoregressive modeling for denoising the GW data and extracting the waveform. We have tested our framework on extracting the injected signals from the simulated data as well as a series of known compact binary coalescence (CBC) events from the LIGO data. Comparing with the conventional whitening procedure, our methodology generally yields improved cross-correlation and reduced root mean square errors with respect to the signal model. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: Phys. Rev. D in press, 16 pages, 11 figures, 1 table

Extract non-Gaussian Features in Gravitational Wave Observation Data Using Self-Supervised Learning

Authors: Yu-Chiung Lin, Albert K. H. Kong

Abstract: We propose a self-supervised learning model to denoise gravitational wave (GW) signals in the time series strain data without relying on waveform information. Denoising GW data is a crucial intermediate process for machine-learning-based data analysis techniques, as it can simplify the model for downstream tasks such as detections and parameter estimations. We use the blind-spot neural network and… ▽ More We propose a self-supervised learning model to denoise gravitational wave (GW) signals in the time series strain data without relying on waveform information. Denoising GW data is a crucial intermediate process for machine-learning-based data analysis techniques, as it can simplify the model for downstream tasks such as detections and parameter estimations. We use the blind-spot neural network and train it with whitened strain data with GW signals injected as both input data and target. Under the assumption of a Gaussian noise model, our model successfully denoises 38% of GW signals from binary black hole mergers in H1 data and 49% of signals in L1 data detected in the O1, O2, and O3 observation runs with an overlap greater than 0.5. We also test the model's potential to extract glitch features and loud inspiral compact binary coalescence signals a few seconds before the merger. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: 39 pages, 15 figures in the main article, and 43 figures in the appendix

Ultralight vector dark matter search using data from the KAGRA O3GK run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi , et al. (1778 additional authors not shown)

Abstract: Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we prese… ▽ More Among the various candidates for dark matter (DM), ultralight vector DM can be probed by laser interferometric gravitational wave detectors through the measurement of oscillating length changes in the arm cavities. In this context, KAGRA has a unique feature due to differing compositions of its mirrors, enhancing the signal of vector DM in the length change in the auxiliary channels. Here we present the result of a search for $U(1)_{B-L}$ gauge boson DM using the KAGRA data from auxiliary length channels during the first joint observation run together with GEO600. By applying our search pipeline, which takes into account the stochastic nature of ultralight DM, upper bounds on the coupling strength between the $U(1)_{B-L}$ gauge boson and ordinary matter are obtained for a range of DM masses. While our constraints are less stringent than those derived from previous experiments, this study demonstrates the applicability of our method to the lower-mass vector DM search, which is made difficult in this measurement by the short observation time compared to the auto-correlation time scale of DM. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: 20 pages, 5 figures

Report number: LIGO-P2300250

Rapid Generation of Kilonova Light Curves Using Conditional Variational Autoencoder

Authors: Surojit Saha, Michael J. Williams, Laurence Datrier, Fergus Hayes, Matt Nicholl, Albert K. H. Kong, Martin Hendry, IK Siong Heng, Gavin P. Lamb, En-Tzu Lin, Daniel Williams

Abstract: The discovery of the optical counterpart, along with the gravitational waves from GW170817, of the first binary neutron star merger, opened up a new era for multi-messenger astrophysics. Combining the GW data with the optical counterpart, also known as AT2017gfo, classified as a kilonova, has revealed the nature of compact binary merging systems by extracting enriched information about the total b… ▽ More The discovery of the optical counterpart, along with the gravitational waves from GW170817, of the first binary neutron star merger, opened up a new era for multi-messenger astrophysics. Combining the GW data with the optical counterpart, also known as AT2017gfo, classified as a kilonova, has revealed the nature of compact binary merging systems by extracting enriched information about the total binary mass, the mass ratio, the system geometry, and the equation of state. Even though the detection of kilonova brought about a revolution in the domain of multi-messenger astronomy, since there has been only one kilonova from a gravitational wave detected binary neutron star merger event so far, this limits the exact understanding of the origin and propagation of the kilonova. Here, we use a conditional variational autoencoder trained on light curve data from two kilonova models having different temporal lengths, and consequently, generate kilonova light curves rapidly based on physical parameters of our choice with good accuracy. Once trained, the time scale for light curve generation is of the order of a few milliseconds, thus speeding up generating light curves by $1000$ times compared to the simulation. The mean squared error between the generated and original light curves is typically $0.015$ with a maximum of $0.08$ for each set of considered physical parameter; while having a maximum of $\approx0.6$ error across the whole parameter space. Hence, implementing this technique provides fast and reliably accurate results. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 19 pages, 7 figures (3 additional figures in appendix), accepted to ApJ

A Joint Fermi-GBM and Swift-BAT Analysis of Gravitational-Wave Candidates from the Third Gravitational-wave Observing Run

Authors: C. Fletcher, J. Wood, R. Hamburg, P. Veres, C. M. Hui, E. Bissaldi, M. S. Briggs, E. Burns, W. H. Cleveland, M. M. Giles, A. Goldstein, B. A. Hristov, D. Kocevski, S. Lesage, B. Mailyan, C. Malacaria, S. Poolakkil, A. von Kienlin, C. A. Wilson-Hodge, The Fermi Gamma-ray Burst Monitor Team, M. Crnogorčević, J. DeLaunay, A. Tohuvavohu, R. Caputo, S. B. Cenko , et al. (1674 additional authors not shown)

Abstract: We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses,… ▽ More We present Fermi Gamma-ray Burst Monitor (Fermi-GBM) and Swift Burst Alert Telescope (Swift-BAT) searches for gamma-ray/X-ray counterparts to gravitational wave (GW) candidate events identified during the third observing run of the Advanced LIGO and Advanced Virgo detectors. Using Fermi-GBM on-board triggers and sub-threshold gamma-ray burst (GRB) candidates found in the Fermi-GBM ground analyses, the Targeted Search and the Untargeted Search, we investigate whether there are any coincident GRBs associated with the GWs. We also search the Swift-BAT rate data around the GW times to determine whether a GRB counterpart is present. No counterparts are found. Using both the Fermi-GBM Targeted Search and the Swift-BAT search, we calculate flux upper limits and present joint upper limits on the gamma-ray luminosity of each GW. Given these limits, we constrain theoretical models for the emission of gamma-rays from binary black hole mergers. △ Less

Submitted 25 August, 2023; originally announced August 2023.

Chandra Observation of NGC 1559: Eight Ultraluminous X-ray Sources Including a Compact Binary Candidate

Authors: Chen-Hsun Ma, Kwan-Lok Li, You-Hua Chu, Albert K. H. Kong

Abstract: Despite the 30-year history of ultra-luminous X-ray sources (ULXs) studies, issues like the majority of their physical natures (i.e., neutron stars, stellar-mass black holes, or intermediate black holes) as well as the accretion mechanisms are still under debate. Expanding the ULX sample size in the literature is clearly a way to help. To this end, we investigated the X-ray source population, ULXs… ▽ More Despite the 30-year history of ultra-luminous X-ray sources (ULXs) studies, issues like the majority of their physical natures (i.e., neutron stars, stellar-mass black holes, or intermediate black holes) as well as the accretion mechanisms are still under debate. Expanding the ULX sample size in the literature is clearly a way to help. To this end, we investigated the X-ray source population, ULXs in particular, in the barred spiral galaxy NGC 1559 using a Chandra observation made in 2016. In this 45-ks exposure, 33 X-ray point sources were detected within the 2.'7 isophotal radius of the galaxy. Among them, 8 ULXs were identified with the criterion of the X-ray luminosity $L_x>10^{39}$ erg s$^{-1}$ (0.3-7~keV). Both X-ray light curves and spectra of all the sources were examined. Except for some low-count spectra that only provide ambiguous spectral fitting results, all the X-ray sources were basically spectrally hard and therefore likely have non-thermal origins. While no strong X-ray variability was present in most of the sources owing to the relatively short exposure of the observation, we found an intriguing ULX, named X-24, exhibiting a periodicity of $\sim$7500s with a detection significance of 2.7$σ$. We speculate that it is the orbital period of the system. Roche-lobe over flow and Roche limit are consistent with the speculation. Thus, we suggest that X-24 may be the one of the rare compact binary ULXs, and hence, a good candidate as a stellar-mass black hole. △ Less

Submitted 10 August, 2023; originally announced August 2023.

Comments: Accepted for publication in ApJ

Influences of dynamical disruptions on the evolution of pulsars in globular clusters

Authors: Kwangmin Oh, C. Y. Hui, Jongsuk Hong, J. Takata, A. K. H. Kong, Pak-Hin Thomas Tam, Kwan-Lok Li, K. S. Cheng

Abstract: By comparing the physical properties of pulsars hosted by core-collapsed (CCed) and non-core-collapsed (Non-CCed) globular clusters (GCs), we find that pulsars in CCed GCs rotate significantly slower than their counterparts in Non-CCed GCs. Additionally, radio luminosities at 1.4 GHz in CCed GCs are higher. These findings are consistent with the scenario that dynamical interactions in GCs can inte… ▽ More By comparing the physical properties of pulsars hosted by core-collapsed (CCed) and non-core-collapsed (Non-CCed) globular clusters (GCs), we find that pulsars in CCed GCs rotate significantly slower than their counterparts in Non-CCed GCs. Additionally, radio luminosities at 1.4 GHz in CCed GCs are higher. These findings are consistent with the scenario that dynamical interactions in GCs can interrupt angular momentum transfer processes and surface magnetic field decay during the recycling phase. Our results suggest that such effects in CCed GCs are stronger due to more frequent disruptions of compact binaries. This is further supported by the observation that both estimated disruption rates and the fraction of isolated pulsars are predominantly higher in CCed GCs. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 9 pages, 8 figures, 3 tables, Accepted in MNRAS

Cosmic ray calorimetry in star-forming galaxy populations and implications for their contribution to the extra-galactic $γ$-ray background

Authors: Ellis R. Owen, Albert K. H. Kong, Kuo-Chuan Pan

Abstract: Star-forming galaxies (SFGs) have been established as an important source population in the extra-galactic $γ$-ray background (EGB). Their intensive star-formation creates an abundance of environments able to accelerate particles, and these build-up a rich sea of cosmic rays (CRs). Above GeV energies, CR protons can undergo hadronic interactions with their environment to produce $γ$-rays. SFGs can… ▽ More Star-forming galaxies (SFGs) have been established as an important source population in the extra-galactic $γ$-ray background (EGB). Their intensive star-formation creates an abundance of environments able to accelerate particles, and these build-up a rich sea of cosmic rays (CRs). Above GeV energies, CR protons can undergo hadronic interactions with their environment to produce $γ$-rays. SFGs can operate as CR proton "calorimeters", where a large fraction of the CR energy is converted to $γ$-rays. However, CRs also deposit energy and momentum to modify the thermal and hydrodynamic conditions of the gas in SFGs, and can become a powerful driver of outflows. Such outflows are ubiquitous among some types of SFGs, and have the potential to severely degrade their CR proton calorimetry. This diminishes their contribution to the EGB. In this work, we adopt a self-consistent treatment of particle transport in outflows from SFGs to assess their calorimetry. We use 1D numerical treatments of galactic outflows driven by CRs and thermal gas pressure, accounting for the dynamical effects and interactions of CRs. We show the impact CR-driven flows have on the relative contribution of SFG populations to the EGB, and investigate the properties of SFGs that contribute most strongly. △ Less

Submitted 9 August, 2023; originally announced August 2023.

Comments: 8 pages, 4 figures, 1 table. Presented at the 38th International Cosmic Ray Conference (ICRC2023)

Journal ref: PoS (ICRC2023), 554

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Evidence of stellar oscillations in the post-common envelop binary candidate ASASSN-V J205543.90+240033.5

Authors: J. Takata, A. K. H. Kong, X. F. Wang, F. F. Song, J. Mao, X. Hou, C. -P. Hu, L. C. -C. Lin, K. L. Li, C. Y. Hui

Abstract: ASASSN-V J205543.90+240033.5 (ASJ2055) is a possible post-common envelope binary system. Its optical photometric data shows an orbital variation about $0.52$~days and a fast period modulation of $P_0\sim 9.77$~minute, whose origin is unknown. In this {\it Letter}, we report an evidence of the stellar oscillation of the companion star as the origin of the fast period modulation. We analyze the phot… ▽ More ASASSN-V J205543.90+240033.5 (ASJ2055) is a possible post-common envelope binary system. Its optical photometric data shows an orbital variation about $0.52$~days and a fast period modulation of $P_0\sim 9.77$~minute, whose origin is unknown. In this {\it Letter}, we report an evidence of the stellar oscillation of the companion star as the origin of the fast period modulation. We analyze the photometric data taken by TESS, Liverpool telescope, and Lulin One-meter Telescope. It is found that the period of the 9.77-minute signal measured in 2022 August is significantly shorter than that in 2021 July/August, and the magnitude of the change is of the order of $|\triangle P_0|/P_0\sim 0.0008(4)$. Such a large variation will be incompatible with the scenario of the white dwarf spin as the origin of the 9.77-minute periodic modulation. We suggest that the fast periodic signal is related to the emission from the irradiated companion star rather than that of the white dwarf. Using existing photometric data covering a wide wavelength range, we estimate that the hot white dwarf in ASJ2055 has a temperature of $T_{eff}\sim 80000$~K and is heating the oscillating M-type main-sequence star with $T_{eff}\sim 3500$~K on its un-irradiated surface. The stellar oscillation of M-type main-sequence star has been predicted in theoretical studies, but no observational confirmation has been done. ASJ2055, therefore, has a potential to be a unique laboratory to investigate the stellar oscillation of a M-type main-sequence star and the heating effect on the stellar oscillation. △ Less

Submitted 3 August, 2023; v1 submitted 1 August, 2023; originally announced August 2023.

Comments: 9 pages, 5 figures, 2 tables. Accepted for publication in ApJ Letter

Local primordial non-Gaussianity from the large-scale clustering of photometric DESI luminous red galaxies

Authors: Mehdi Rezaie, Ashley J. Ross, Hee-Jong Seo, Hui Kong, Anna Porredon, Lado Samushia, Edmond Chaussidon, Alex Krolewski, Arnaud de Mattia, Florian Beutler, Jessica Nicole Aguilar, Steven Ahlen, Shadab Alam, Santiago Avila, Benedict Bahr-Kalus, Jose Bermejo-Climent, David Brooks, Todd Claybaugh, Shaun Cole, Kyle Dawson, Axel de la Macorra, Peter Doel, Andreu Font-Ribera, Jaime E. Forero-Romero, Satya Gontcho A Gontcho , et al. (24 additional authors not shown)

Abstract: We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter $\fnl$. Our sample comprises over 12 million targets, covering 14,000 square degrees of the sky, with redshifts in the range $0.2< z < 1.35$. We identify Galactic extinction, survey depth, and astronomical seeing as the… ▽ More We use angular clustering of luminous red galaxies from the Dark Energy Spectroscopic Instrument (DESI) imaging surveys to constrain the local primordial non-Gaussianity parameter $\fnl$. Our sample comprises over 12 million targets, covering 14,000 square degrees of the sky, with redshifts in the range $0.2< z < 1.35$. We identify Galactic extinction, survey depth, and astronomical seeing as the primary sources of systematic error, and employ linear regression and artificial neural networks to alleviate non-cosmological excess clustering on large scales. Our methods are tested against simulations with and without $\fnl$ and systematics, showing superior performance of the neural network treatment. The neural network with a set of nine imaging property maps passes our systematic null test criteria, and is chosen as the fiducial treatment. Assuming the universality relation, we find $\fnl = 34^{+24(+50)}_{-44(-73)}$ at 68\%(95\%) confidence. We apply a series of robustness tests (e.g., cuts on imaging, declination, or scales used) that show consistency in the obtained constraints. We study how the regression method biases the measured angular power-spectrum and degrades the $\fnl$ constraining power. The use of the nine maps more than doubles the uncertainty compared to using only the three primary maps in the regression. Our results thus motivate the development of more efficient methods that avoid over-correction, protect large-scale clustering information, and preserve constraining power. Additionally, our results encourage further studies of $\fnl$ with DESI spectroscopic samples, where the inclusion of 3D clustering modes should help separate imaging systematics and lessen the degradation in the $\fnl$ uncertainty. △ Less

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

Comments: 21 pages, 17 figures, 7 tables (Appendix excluded). Published in MNRAS

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

A journey from the hard to the soft state: How do QPOs evolve in the 2021 outburst of GX 339-4?

Authors: H. Stiele, A. K. H. Kong

Abstract: We investigated the snapshots of five NICER observations of the black hole transient GX 339-4 when the source transited from the hard state into the soft state during its outburst in 2021. In this paper, we focused our study on the evolution of quasi-periodic oscillations (QPOs) and noise components using power-density spectra. In addition, we derived hardness ratios comparing count rates above an… ▽ More We investigated the snapshots of five NICER observations of the black hole transient GX 339-4 when the source transited from the hard state into the soft state during its outburst in 2021. In this paper, we focused our study on the evolution of quasi-periodic oscillations (QPOs) and noise components using power-density spectra. In addition, we derived hardness ratios comparing count rates above and below 2 keV. The evolution from the hard to the soft state was a somewhat erratic process showing several transitions between states that are dominated by top-flat noise and can show type-C QPOs; those that are dominated by red noise and can show type-B QPOs. From the parameters that we studied, we only found a strong correlation between the hardness ratio and the type of QPO observed. This implies that the appearance of type-B QPOs is related to a change in the accretion geometry of the system that also reflects in altered spectral properties. We also observed that the type-B QPO forms from or disintegrates into a broad peaked feature when the source comes out of or goes to the hard-intermediate state, respectively. This implies some strong decoherence in the process that creates this feature. △ Less

Submitted 12 April, 2023; originally announced April 2023.

Comments: 6 pages, 5 figures, supplementary online material as appendices (13 pages), accepted for publication in MNRAS

A comparison of millisecond pulsar populations between globular clusters and the Galactic field

Authors: Jongsu Lee, C. Y. Hui, J. Takata, A. K. H. Kong, Pak-Hin Thomas Tam, Kwan-Lok Li, K. S. Cheng

Abstract: We have performed a systematic study of the rotational, orbital and X-ray properties of millisecond pulsars (MSPs) in globular clusters (GCs) and compared their nature with those of the MSPs in the Galactic field (GF). We found that GC MSPs generally rotate slower than their counterparts in the GF. Different from the expectation of a simple recycling scenario, no evidence for the correlation betwe… ▽ More We have performed a systematic study of the rotational, orbital and X-ray properties of millisecond pulsars (MSPs) in globular clusters (GCs) and compared their nature with those of the MSPs in the Galactic field (GF). We found that GC MSPs generally rotate slower than their counterparts in the GF. Different from the expectation of a simple recycling scenario, no evidence for the correlation between the orbital period and the rotation period can be found from the MSP binaries in GCs. There is also an indication that the surface magnetic field of GC MSPs are stronger than those in the GF. All these suggest dynamical interactions in GCs can alter the evolution of MSPs/their progenitors which can leave an imprint on their X-ray emission properties. While the MSPs in both GF and GCs have similar distributions of X-ray luminosity and hardness, our sample supports the notion that these two populations follow different relation between the X-ray luminosity and spin-down power. We discuss this in terms of both pulsar emission model and the observational bias. △ Less

Submitted 17 February, 2023; originally announced February 2023.

Comments: 33 pages, 13 figures, 10 tables, Accepted for publication in ApJ

The optical light curve of GRB 221009A: the afterglow and the emerging supernova

Authors: M. D. Fulton, S. J. Smartt, L. Rhodes, M. E. Huber, A. V. Villar, T. Moore, S. Srivastav, A. S. B. Schultz, K. C. Chambers, L. Izzo, J. Hjorth, T. -W. Chen, M. Nicholl, R. J. Foley, A. Rest, K. W. Smith, D. R. Young, S. A. Sim, J. Bright, Y. Zenati, T. de Boer, J. Bulger, J. Fairlamb, H. Gao, C. -C. Lin , et al. (24 additional authors not shown)

Abstract: We present extensive optical photometry of the afterglow of GRB~221009A. Our data cover $0.9 - 59.9$\,days from the time of \textit{Swift} and \textit{Fermi} GRB detections. Photometry in $rizy$-band filters was collected primarily with Pan-STARRS and supplemented by multiple 1- to 4-meter imaging facilities. We analyzed the Swift X-ray data of the afterglow and found a single decline rate power-l… ▽ More We present extensive optical photometry of the afterglow of GRB~221009A. Our data cover $0.9 - 59.9$\,days from the time of \textit{Swift} and \textit{Fermi} GRB detections. Photometry in $rizy$-band filters was collected primarily with Pan-STARRS and supplemented by multiple 1- to 4-meter imaging facilities. We analyzed the Swift X-ray data of the afterglow and found a single decline rate power-law $f(t) \propto t^{-1.556\pm0.002}$ best describes the light curve. In addition to the high foreground Milky Way dust extinction along this line of sight, the data favour additional extinction to consistently model the optical to X-ray flux with optically thin synchrotron emission. We fit the X-ray-derived power-law to the optical light curve and find good agreement with the measured data up to $5-6$\,days. Thereafter we find a flux excess in the $riy$ bands which peaks in the observer frame at $\sim20$\,days. This excess shares similar light curve profiles to the type Ic broad-lined supernovae SN~2016jca and SN~2017iuk once corrected for the GRB redshift of $z=0.151$ and arbitrarily scaled. This may be representative of a supernova emerging from the declining afterglow. We measure rest-frame absolute peak AB magnitudes of $M_g=-19.8\pm0.6$ and $M_r=-19.4\pm0.3$ and $M_z=-20.1\pm0.3$. If this is an SN component, then Bayesian modelling of the excess flux would imply explosion parameters of $M_{\rm ej}=7.1^{+2.4}_{-1.7}$ M$_{\odot}$, $M_{\rm Ni}=1.0^{+0.6}_{-0.4}$ M$_{\odot}$, and $v_{\rm ej}=33,900^{+5,900}_{-5,700} kms^{-1}$, for the ejecta mass, nickel mass and ejecta velocity respectively, inferring an explosion energy of $E_{\rm kin}\simeq 2.6-9.0\times10^{52}$ ergs. △ Less

Submitted 23 March, 2023; v1 submitted 25 January, 2023; originally announced January 2023.

Comments: Submitted to ApJL on 23rd January 2023, for consideration for publication in the special issue on GRB 221009A. Accepted on 2nd March 2023. The results of this paper are under press embargo until 28th March 2023. 15 pages, 7 figures, 1 table

Multi-Wavelength Observations Of A New Redback Millisecond Pulsar 4FGL J1910.7-5320

Authors: Ka-Yui Au, Jay Strader, Samuel J. Swihart, Lupin C. C. Lin, Albert K. H. Kong, Jumpei Takata, Chung-Yue Hui, Teresa Panurach, Isabella Molina, Elias Aydi, Kirill Sokolovsky, Kwan-Lok Li

Abstract: We present the study of multi-wavelength observations of an unidentified Fermi Large Area Telescope (LAT) source, 4FGL J1910.7-5320, a new candidate redback millisecond pulsar binary. In the 4FGL 95% error region of 4FGL J1910.7-5320, we find a possible binary with a 8.36-hr orbital period from the Catalina Real-Time Transient Survey (CRTS), confirmed by optical spectroscopy using the SOAR telesco… ▽ More We present the study of multi-wavelength observations of an unidentified Fermi Large Area Telescope (LAT) source, 4FGL J1910.7-5320, a new candidate redback millisecond pulsar binary. In the 4FGL 95% error region of 4FGL J1910.7-5320, we find a possible binary with a 8.36-hr orbital period from the Catalina Real-Time Transient Survey (CRTS), confirmed by optical spectroscopy using the SOAR telescope. This optical source was recently independently discovered as a redback pulsar by the TRAPUM project, confirming our prediction. We fit the optical spectral energy distributions of 4FGL J1910.7-5320 with a blackbody model, inferring a maximum distance of 4.1 kpc by assuming that the companion fills its Roche-lobe with a radius of R = 0.7R_sun. Using a 12.6 ks Chandra X-ray observation, we identified an X-ray counterpart for 4FGL J1910.7-5320, with a spectrum that can be described by an absorbed power-law with a photon index of 1.0+/-0.4. The spectrally hard X-ray emission shows tentative evidence for orbital variability. Using more than 12 years of Fermi-LAT data, we refined the position of the γ-ray source, and the optical candidate still lies within the 68% positional error circle. In addition to 4FGL J1910.7-5320, we find a variable optical source with a periodic signal of 4.28-hr inside the 4FGL catalog 95% error region of another unidentified Fermi source, 4FGL J2029.5-4237. However, the γ-ray source does not have a significant X-ray counterpart in a 11.7 ks Chandra observation, with a 3-σ flux upper limit of 2.4*10^-14 erg cm^-2 s^-1 (0.3-7 keV). Moreover, the optical source is outside our updated Fermi-LAT 95% error circle. These observational facts all suggest that this new redback millisecond pulsar powers the γ-ray source 4FGL J1910.7-5320 while 4FGL J2029.5-4237 is unlikely the γ-ray counterpart to the 4.28-hr variable. △ Less

Submitted 26 December, 2022; v1 submitted 22 December, 2022; originally announced December 2022.

Comments: Accepted for publication in ApJ

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

Search for gravitational-wave transients associated with magnetar bursts in Advanced LIGO and Advanced Virgo data from the third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1645 additional authors not shown)

Abstract: Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bu… ▽ More Gravitational waves are expected to be produced from neutron star oscillations associated with magnetar giant flares and short bursts. We present the results of a search for short-duration (milliseconds to seconds) and long-duration ($\sim$ 100 s) transient gravitational waves from 13 magnetar short bursts observed during Advanced LIGO, Advanced Virgo and KAGRA's third observation run. These 13 bursts come from two magnetars, SGR 1935$+$2154 and Swift J1818.0$-$1607. We also include three other electromagnetic burst events detected by Fermi GBM which were identified as likely coming from one or more magnetars, but they have no association with a known magnetar. No magnetar giant flares were detected during the analysis period. We find no evidence of gravitational waves associated with any of these 16 bursts. We place upper bounds on the root-sum-square of the integrated gravitational-wave strain that reach $2.2 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at 100 Hz for the short-duration search and $8.7 \times 10^{-23}$ $/\sqrt{\text{Hz}}$ at $450$ Hz for the long-duration search, given a detection efficiency of 50%. For a ringdown signal at 1590 Hz targeted by the short-duration search the limit is set to $1.8 \times 10^{-22}$ $/\sqrt{\text{Hz}}$. Using the estimated distance to each magnetar, we derive upper bounds on the emitted gravitational-wave energy of $3.2 \times 10^{43}$ erg ($7.3 \times 10^{43}$ erg) for SGR 1935$+$2154 and $8.2 \times 10^{42}$ erg ($2.8 \times 10^{43}$ erg) for Swift J1818.0$-$1607, for the short-duration (long-duration) search. Assuming isotropic emission of electromagnetic radiation of the burst fluences, we constrain the ratio of gravitational-wave energy to electromagnetic energy for bursts from SGR 1935$+$2154 with available fluence information. The lowest of these ratios is $3 \times 10^3$. △ Less

Submitted 19 October, 2022; originally announced October 2022.

Comments: 30 pages with appendices, 5 figures, 10 tables

Report number: LIGO-P2100387

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

Searching for cataclysmic variable stars in unidentified X-ray sources

Authors: J. Takata, X. F. Wang, A. K. H. Kong, J. Mao, X. Hou, C. -P. Hu, L. C. -C. Lin, K. L. Li, C. Y. Hui

Abstract: We carry out a photometric search for new cataclysmic variable stars (CVs), with the goal of identification for candidates of AR~Scorpii-type binary systems. We select GAIA sources that are likely associated with unidentified X-ray sources, and analyze the light curves taken by the Zwicky Transient Facility, Transiting Exoplanet Survey Satellite, and Lulin One-meter Telescope in Taiwan. We investi… ▽ More We carry out a photometric search for new cataclysmic variable stars (CVs), with the goal of identification for candidates of AR~Scorpii-type binary systems. We select GAIA sources that are likely associated with unidentified X-ray sources, and analyze the light curves taken by the Zwicky Transient Facility, Transiting Exoplanet Survey Satellite, and Lulin One-meter Telescope in Taiwan. We investigate eight sources as candidates for CVs, among which six sources are new identifications. Another two sources have been recognized as CVs in previous studies, but no detailed investigations have been done. We identify two eclipsing systems that are associated with an unidentified XMM-Newton or Swift source, and one promising candidate for polar associated with an unidentified ASKA source. Two polar candidates may locate in the so-called period gap of a CV, and the other six candidates have an orbital period shorter than that of the period gap. Although we do not identify a promising candidate for AR~Scorpii-type binary systems, our study suggests that CV systems that have X-ray emission and do not show frequent outbursts may have been missed in previous surveys. △ Less

Submitted 24 August, 2022; v1 submitted 2 August, 2022; originally announced August 2022.

Comments: 21 pages, 24 figures, accepted for publication in ApJ. "Note added in proof" and "Appendix C" are added in the revised version

A Comprehensive Analysis of the Gravitational Wave Events with the Hilbert-Huang Transform: From Compact Binary Coalescence to Supernova

Authors: Chin-Ping Hu, Lupin Chun-Che Lin, Kuo-Chuan Pan, Kwan-Lok Li, Chien-Chang Yen, Albert K. H. Kong, C. Y. Hui

Abstract: We analyze the gravitational wave signals with a model-independent time-frequency analysis, which is improved from the Hilbert-Huang transform (HHT) and optimized for characterizing the frequency variability on the time-frequency map. Except for the regular HHT algorithm, i.e., obtaining intrinsic mode functions with ensemble empirical mode decomposition and yielding the instantaneous frequencies,… ▽ More We analyze the gravitational wave signals with a model-independent time-frequency analysis, which is improved from the Hilbert-Huang transform (HHT) and optimized for characterizing the frequency variability on the time-frequency map. Except for the regular HHT algorithm, i.e., obtaining intrinsic mode functions with ensemble empirical mode decomposition and yielding the instantaneous frequencies, we propose an alternative algorithm that operates the ensemble mean on the time-frequency map. We systematically analyze the known gravitational wave events of the compact binary coalescence observed in LIGO O1 and O2, and in the simulated gravitational wave signals from core-collapse supernovae (CCSNe) with our method. The time-frequency maps of the binary black hole coalescence cases show much better details compared to those wavelet spectra. Moreover, the oscillation in the instantaneous frequency caused by mode-mixing could be reduced with our algorithm. For the CCSNe data, the oscillation from the proto-neutron star and the radiation from the standing accretion shock instability can be precisely determined with the HHT in great detail. More importantly, the initial stage of different modes of oscillations can be clearly separated. These results provide new hints for further establishment of the detecting algorithm, and new probes to investigate the underlying physical mechanisms. △ Less

Submitted 14 July, 2022; originally announced July 2022.

Comments: 24 pages, 26 figures, 1 table, accepted for publication in ApJ

Noise subtraction from KAGRA O3GK data using Independent Component Analysis

Authors: KAGRA collaboration, H. Abe, T. Akutsu, M. Ando, A. Araya, N. Aritomi, H. Asada, Y. Aso, S. Bae, Y. Bae, R. Bajpai, K. Cannon, Z. Cao, E. Capocasa, M. Chan, C. Chen, D. Chen, K. Chen, Y. Chen, C-Y. Chiang, Y-K. Chu, S. Eguchi, M. Eisenmann, Y. Enomoto, R. Flaminio , et al. (178 additional authors not shown)

Abstract: In April 2020, KAGRA conducted its first science observation in combination with the GEO~600 detector (O3GK) for two weeks. According to the noise budget estimation, suspension control noise in the low frequency band and acoustic noise in the middle frequency band are identified as the dominant contribution. In this study, we show that such noise can be reduced in offline data analysis by utilizin… ▽ More In April 2020, KAGRA conducted its first science observation in combination with the GEO~600 detector (O3GK) for two weeks. According to the noise budget estimation, suspension control noise in the low frequency band and acoustic noise in the middle frequency band are identified as the dominant contribution. In this study, we show that such noise can be reduced in offline data analysis by utilizing a method called Independent Component Analysis (ICA). Here the ICA model is extended from the one studied in iKAGRA data analysis by incorporating frequency dependence while linearity and stationarity of the couplings are still assumed. By using optimal witness sensors, those two dominant contributions are mitigated in the real observational data. We also analyze the stability of the transfer functions for whole two weeks data in order to investigate how the current subtraction method can be practically used in gravitational wave search. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 26 pages, 10 figures

Report number: JGW-P2214018, RESCEU-5/22

XMM-Newton and NuSTAR observations of the compact millisecond pulsar binary PSR~J1653$-$0158

Authors: Jane SiNan Long, Albert K. H. kong, Kinwah Wu, Jumpei Takata, Qin Han, David C. Y. Hui, Kwan Lok Li

Abstract: We have presented the first joint {\it XMM-Newton} and {\it NuSTAR} analysis of the millisecond pulsar (MSP) binary PSR J1653$-$0158. The 75-minute orbital period inferred from optical and gamma-ray observations together with the 1.97-ms pulsation in the gamma-rays indicate that this system is the most compact Black Widow MSP system known to date. The orbital period was not detected in the {\it XM… ▽ More We have presented the first joint {\it XMM-Newton} and {\it NuSTAR} analysis of the millisecond pulsar (MSP) binary PSR J1653$-$0158. The 75-minute orbital period inferred from optical and gamma-ray observations together with the 1.97-ms pulsation in the gamma-rays indicate that this system is the most compact Black Widow MSP system known to date. The orbital period was not detected in the {\it XMM-Newton} and {\it NuSTAR} data, probably due to insufficient photon counts obtained in the observations. Fitting the joint X-ray spectrum of PSR J1653$-$0158 with a power law gives a photon index $Γ= 1.71 \pm 0.09$. The X-ray luminosity of the source in the ($0.2-40$)~keV band is deduced to be $1.18 \times 10^{31}\;\!{\rm erg~s}^{-1}$, for an adopted distance of 0.84~kpc. We have shown that the broad-band X-ray spectrum can be explained by synchrotron radiation from electrons accelerated in the intra-binary shock, and the gamma-rays detected in the {\it Fermi} data are curvature radiations from electrons and positrons in the pulsar magnetosphere. Our kinematic analysis of the Tidarren systems PSR~J1653--0158 and PSR~J1311--3430 indicates that the two Tidarren systems are likely to have originated in the Galactic Disk. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: Accepted by APJ

A 62-minute orbital period black widow binary in a wide hierarchical triple

Authors: Kevin B. Burdge, Thomas R. Marsh, Jim Fuller, Eric C. Bellm, Ilaria Caiazzo, Deepto Chakrabarty, Michael W. Coughlin, Kishalay De, V. S. Dhillon, Matthew J. Graham, Pablo Rodrí guez-Gil, Amruta D. Jaodand, David L. Kaplan, Erin Kara, Albert K. H. Kong, S. R. Kulkarni, Kwan-Lok Li, S. P. Littlefair, Walid A. Majid, Przemek Mróz, Aaron B. Pearlman, E. S. Phinney, Jan van Roestel, Robert A. Simcoe, Igor Andreoni , et al. (8 additional authors not shown)

Abstract: Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposin… ▽ More Over a dozen millisecond pulsars are ablating low-mass companions in close binary systems. In the original "black widow", the 8-hour orbital period eclipsing pulsar PSR J1959+2048 (PSR B1957+20), high energy emission originating from the pulsar is irradiating and may eventually destroy a low-mass companion. These systems are not only physical laboratories that reveal the dramatic result of exposing a close companion star to the relativistic energy output of a pulsar, but are also believed to harbour some of the most massive neutron stars, allowing for robust tests of the neutron star equation of state. Here, we report observations of ZTF J1406+1222, a wide hierarchical triple hosting a 62-minute orbital period black widow candidate whose optical flux varies by a factor of more than 10. ZTF J1406+1222 pushes the boundaries of evolutionary models, falling below the 80 minute minimum orbital period of hydrogen-rich systems. The wide tertiary companion is a rare low metallicity cool subdwarf star, and the system has a Galactic halo orbit consistent with passing near the Galactic center, making it a probe of formation channels, neutron star kick physics, and binary evolution. △ Less

Submitted 4 May, 2022; originally announced May 2022.

Comments: 63 pages, 15 figures, 2 tables, published in Nature on May 4, 2022

Journal ref: Nature 605, 41-45 (2022)

Multi-wavelength properties of 850-$μ$m selected sources from the North Ecliptic Pole SCUBA-2 survey

Authors: H. Shim, D. Lee, Y. Kim, D. Scott, S. Serjeant, Y. Ao, L. Barrufet, S. C. Chapman, D. Clements, C. J. Conselice, T. Goto, T. R. Greve, H. S. Hwang, M. Im, W. -S. Jeong, H. K. Kim, M. Kim, S. J. Kim, A. K. H. Kong, M. P. Koprowski, M. A. Malkan, M. Michalowski, C. Pearson, H. Seo, T. Takagi , et al. (3 additional authors not shown)

Abstract: We present the multi-wavelength counterparts of 850-$μ$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$μ$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we ide… ▽ More We present the multi-wavelength counterparts of 850-$μ$m selected submillimetre sources over a 2-deg$^2$ field centred on the North Ecliptic Pole. In order to overcome the large beam size (15 arcsec) of the 850-$μ$m images, deep optical to near-infrared (NIR) photometric data and arcsecond-resolution 20-cm images are used to identify counterparts of submillimetre sources. Among 647 sources, we identify 514 reliable counterparts for 449 sources (69 per cent in number), based either on probabilities of chance associations calculated from positional offsets or offsets combined with the optical-to-NIR colours. In the radio imaging, the fraction of 850-$μ$m sources having multiple counterparts is 7 per cent. The photometric redshift, infrared luminosity, stellar mass, star-formation rate (SFR), and the AGN contribution to the total infrared luminosity of the identified counterparts are investigated through spectral energy distribution fitting. The SMGs are infrared-luminous galaxies at an average $\langle z\rangle=2.5$ with $\mathrm{log}_{10} (L_\mathrm{IR}/\mathrm{L}_\odot)=11.5-13.5$, with a mean stellar mass of $\mathrm{log}_{10} (M_\mathrm{star}/\mathrm{M}_\odot)=10.90$ and SFR of $\mathrm{log}_{10} (\mathrm{SFR/M_\odot\,yr^{-1}})=2.34$. The SMGs show twice as large SFR as galaxies on the star-forming main sequence, and about 40 per cent of the SMGs are classified as objects with bursty star formation. At $z\ge4$, the contribution of AGN luminosity to total luminosity for most SMGs is larger than 30 per cent. The FIR-to-radio correlation coefficient of SMGs is consistent with that of main-sequence galaxies at $z\simeq2$. △ Less

Submitted 24 April, 2022; originally announced April 2022.

Comments: 24 pages, 16 figures. Accepted for publication in MNRAS

Search for continuous gravitational wave emission from the Milky Way center in O3 LIGO--Virgo data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1645 additional authors not shown)

Abstract: We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo… ▽ More We present a directed search for continuous gravitational wave (CW) signals emitted by spinning neutron stars located in the inner parsecs of the Galactic Center (GC). Compelling evidence for the presence of a numerous population of neutron stars has been reported in the literature, turning this region into a very interesting place to look for CWs. In this search, data from the full O3 LIGO--Virgo run in the detector frequency band $[10,2000]\rm~Hz$ have been used. No significant detection was found and 95$\%$ confidence level upper limits on the signal strain amplitude were computed, over the full search band, with the deepest limit of about $7.6\times 10^{-26}$ at $\simeq 142\rm~Hz$. These results are significantly more constraining than those reported in previous searches. We use these limits to put constraints on the fiducial neutron star ellipticity and r-mode amplitude. These limits can be also translated into constraints in the black hole mass -- boson mass plane for a hypothetical population of boson clouds around spinning black holes located in the GC. △ Less

Submitted 9 April, 2022; originally announced April 2022.

Comments: 25 pages, 5 figures

Neutrino follow-up with the Zwicky Transient Facility: Results from the first 24 campaigns

Authors: Robert Stein, Simeon Reusch, Anna Franckowiak, Marek Kowalski, Jannis Necker, Sven Weimann, Mansi M. Kasliwal, Jesper Sollerman, Tomas Ahumada, Pau Amaro-Seoane, Shreya Anand, Igor Andreoni, Eric C. Bellm, Joshua S. Bloom, Michael Coughlin, Kishalay De, Christoffer Fremling, Suvi Gezari, Matthew Graham, Steven L. Groom, George Helou, David L. Kaplan, Viraj Karambelkar, Albert K. H. Kong, Erik C. Kool , et al. (11 additional authors not shown)

Abstract: The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutri… ▽ More The Zwicky Transient Facility (ZTF) performs a systematic neutrino follow-up program, searching for optical counterparts to high-energy neutrinos with dedicated Target-of-Opportunity (ToO) observations. Since first light in March 2018, ZTF has taken prompt observations for 24 high-quality neutrino alerts from the IceCube Neutrino Observatory, with a median latency of 12.2 hours from initial neutrino detection. From two of these campaigns, we have already reported tidal disruption event (TDE) AT 2019dsg and likely TDE AT 2019fdr as probable counterparts, suggesting that TDEs contribute >7.8% of the astrophysical neutrino flux. We here present the full results of our program through to December 2021. No additional candidate neutrino sources were identified by our program, allowing us to place the first constraints on the underlying optical luminosity function of astrophysical neutrino sources. Transients with optical absolutes magnitudes brighter that $-21$ can contribute no more than 87% of the total, while transients brighter than $-22$ can contribute no more than 58% of the total, neglecting the effect of extinction and assuming they follow the star formation rate. These are the first observational constraints on the neutrino emission of bright populations such as superluminous supernovae. None of the neutrinos were coincident with bright optical AGN flares comparable to that observed for TXS 0506+056/IC170922A, with such optical blazar flares producing no more than 26% of the total neutrino flux. We highlight the outlook for electromagnetic neutrino follow-up programs, including the expected potential for the Rubin Observatory. △ Less

Submitted 4 April, 2024; v1 submitted 31 March, 2022; originally announced March 2022.

Journal ref: MNRAS, 521, 4, 2023, 5046-5063

Search for Gravitational Waves Associated with Fast Radio Bursts Detected by CHIME/FRB During the LIGO--Virgo Observing Run O3a

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, the CHIME/FRB Collaboration, :, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca , et al. (1633 additional authors not shown)

Abstract: We search for gravitational-wave transients associated with fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB), during the first part of the third observing run of Advanced LIGO and Advanced Virgo (1 April 2019 15:00 UTC-1 Oct 2019 15:00 UTC). Triggers from 22 FRBs were analyzed with a search that targets compact binary coal… ▽ More We search for gravitational-wave transients associated with fast radio bursts (FRBs) detected by the Canadian Hydrogen Intensity Mapping Experiment Fast Radio Burst Project (CHIME/FRB), during the first part of the third observing run of Advanced LIGO and Advanced Virgo (1 April 2019 15:00 UTC-1 Oct 2019 15:00 UTC). Triggers from 22 FRBs were analyzed with a search that targets compact binary coalescences with at least one neutron star component. A targeted search for generic gravitational-wave transients was conducted on 40 FRBs. We find no significant evidence for a gravitational-wave association in either search. Given the large uncertainties in the distances of the FRBs inferred from the dispersion measures in our sample, however, this does not conclusively exclude any progenitor models that include emission of a gravitational wave of the types searched for from any of these FRB events. We report $90\%$ confidence lower bounds on the distance to each FRB for a range of gravitational-wave progenitor models. By combining the inferred maximum distance information for each FRB with the sensitivity of the gravitational-wave searches, we set upper limits on the energy emitted through gravitational waves for a range of emission scenarios. We find values of order $10^{51}$-$10^{57}$ erg for a range of different emission models with central gravitational wave frequencies in the range 70-3560 Hz. Finally, we also found no significant coincident detection of gravitational waves with the repeater, FRB 20200120E, which is the closest known extragalactic FRB. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 35 pages, 6 figures, 8 tables

Report number: P2100124

In search of short gamma-ray burst optical counterpart with the Zwicky Transient Facility

Authors: Tomás Ahumada, Shreya Anand, Michael W. Coughlin, Igor Andreoni, Erik C. Kool, Harsh Kumar, Simeon Reusch, Ana Sagués-Carracedo, Robert Stein, S. Bradley Cenko, Mansi M. Kasliwal, Leo P. Singer, Rachel Dunwoody, Joseph Mangan, Varun Bhalerao, Mattia Bulla, Eric Burns, Matthew J. Graham, David L. Kaplan, Daniel Perley, Mouza Almualla, Joshua S. Bloom, Virginia Cunningham, Kishalay De, Pradip Gatkine , et al. (24 additional authors not shown)

Abstract: The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschi… ▽ More The Fermi Gamma-ray Burst Monitor (GBM) triggers on-board in response to $\sim$ 40 short gamma-ray bursts (SGRBs) per year; however, their large localization regions have made the search for optical counterparts a challenging endeavour. We have developed and executed an extensive program with the wide field of view of the Zwicky Transient Facility (ZTF) camera, mounted on the Palomar 48 inch Oschin telescope (P48), to perform target-of-opportunity (ToO) observations on 10 Fermi-GBM SGRBs during 2018 and 2020-2021. Bridging the large sky areas with small field of view optical telescopes in order to track the evolution of potential candidates, we look for the elusive SGRB afterglows and kilonovae (KNe) associated with these high-energy events. No counterpart has yet been found, even though more than 10 ground based telescopes, part of the Global Relay of Observatories Watching Transients Happen (GROWTH) network, have taken part in these efforts. The candidate selection procedure and the follow-up strategy have shown that ZTF is an efficient instrument for searching for poorly localized SGRBs, retrieving a reasonable number of candidates to follow-up and showing promising capabilities as the community approaches the multi-messenger era. Based on the median limiting magnitude of ZTF, our searches would have been able to retrieve a GW170817-like event up to $\sim$ 200 Mpc and SGRB afterglows to z = 0.16 or 0.4, depending on the assumed underlying energy model. Future ToOs will expand the horizon to z = 0.2 and 0.7 respectively. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Performance of the KAGRA detector during the first joint observation with GEO 600 (O3GK)

Authors: KAGRA Collaboration, H. Abe, R. X. Adhikari, T. Akutsu, M. Ando, A. Araya, N. Aritomi, H. Asada, Y. Aso, S. Bae, Y. Bae, R. Bajpai, S. W. Ballmer, K. Cannon, Z. Cao, E. Capocasa, M. Chan, C. Chen, D. Chen, K. Chen, Y. Chen, C-Y. Chiang, Y-K. Chu, J. C. Driggers, S. E. Dwyer , et al. (193 additional authors not shown)

Abstract: KAGRA, the kilometer-scale underground gravitational-wave detector, is located at Kamioka, Japan. In April 2020, an astrophysics observation was performed at the KAGRA detector in combination with the GEO 600 detector; this observation operation is called O3GK. The optical configuration in O3GK is based on a power recycled Fabry-Pérot Michelson interferometer; all the mirrors were set at room temp… ▽ More KAGRA, the kilometer-scale underground gravitational-wave detector, is located at Kamioka, Japan. In April 2020, an astrophysics observation was performed at the KAGRA detector in combination with the GEO 600 detector; this observation operation is called O3GK. The optical configuration in O3GK is based on a power recycled Fabry-Pérot Michelson interferometer; all the mirrors were set at room temperature. The duty factor of the operation was approximately 53%, and the strain sensitivity was $3\times10^{-22}~/\sqrt{\rm{Hz}}$ at 250 Hz. In addition, the binary-neutron-star (BNS) inspiral range was approximately 0.6 Mpc. The contributions of various noise sources to the sensitivity of O3GK were investigated to understand how the observation range could be improved; this study is called a "noise budget". According to our noise budget, the measured sensitivity could be approximated by adding up the effect of each noise. The sensitivity was dominated by noise from the sensors used for local controls of the vibration isolation systems, acoustic noise, shot noise, and laser frequency noise. Further, other noise sources that did not limit the sensitivity were investigated. This paper provides a detailed account of the KAGRA detector in O3GK including interferometer configuration, status, and noise budget. In addition, strategies for future sensitivity improvements such as hardware upgrades, are discussed. △ Less

Submitted 14 March, 2022; originally announced March 2022.

Report number: JGW-P2113405

First joint observation by the underground gravitational-wave detector, KAGRA, with GEO600

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1647 additional authors not shown)

Abstract: We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing… ▽ More We report the results of the first joint observation of the KAGRA detector with GEO600. KAGRA is a cryogenic and underground gravitational-wave detector consisting of a laser interferometer with three-kilometer arms, and located in Kamioka, Gifu, Japan. GEO600 is a British--German laser interferometer with 600 m arms, and located near Hannover, Germany. GEO600 and KAGRA performed a joint observing run from April 7 to 20, 2020. We present the results of the joint analysis of the GEO--KAGRA data for transient gravitational-wave signals, including the coalescence of neutron-star binaries and generic unmodeled transients. We also perform dedicated searches for binary coalescence signals and generic transients associated with gamma-ray burst events observed during the joint run. No gravitational-wave events were identified. We evaluate the minimum detectable amplitude for various types of transient signals and the spacetime volume for which the network is sensitive to binary neutron-star coalescences. We also place lower limits on the distances to the gamma-ray bursts analysed based on the non-detection of an associated gravitational-wave signal for several signal models, including binary coalescences. These analyses demonstrate the feasibility and utility of KAGRA as a member of the global gravitational-wave detector network. △ Less

Submitted 19 August, 2022; v1 submitted 2 March, 2022; originally announced March 2022.

Comments: Matches with published version

Report number: LIGO-P2100286

Journal ref: Progress of Theoretical and Experimental Physics, Volume 2022, Issue 6, 063F01 (2022)

All-sky search for continuous gravitational waves from isolated neutron stars using Advanced LIGO and Advanced Virgo O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1645 additional authors not shown)

Abstract: We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivativ… ▽ More We present results of an all-sky search for continuous gravitational waves which can be produced by spinning neutron stars with an asymmetry around their rotation axis, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. Four different analysis methods are used to search in a gravitational-wave frequency band from 10 to 2048 Hz and a first frequency derivative from $-10^{-8}$ to $10^{-9}$ Hz/s. No statistically-significant periodic gravitational-wave signal is observed by any of the four searches. As a result, upper limits on the gravitational-wave strain amplitude $h_0$ are calculated. The best upper limits are obtained in the frequency range of 100 to 200 Hz and they are ${\sim}1.1\times10^{-25}$ at 95\% confidence-level. The minimum upper limit of $1.10\times10^{-25}$ is achieved at a frequency 111.5 Hz. We also place constraints on the rates and abundances of nearby planetary- and asteroid-mass primordial black holes that could give rise to continuous gravitational-wave signals. △ Less

Submitted 3 January, 2022; originally announced January 2022.

Comments: 23 main text pages, 17 figures

Report number: LIGO-P2100367

Narrowband searches for continuous and long-duration transient gravitational waves from known pulsars in the LIGO-Virgo third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato , et al. (1636 additional authors not shown)

Abstract: Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational… ▽ More Isolated neutron stars that are asymmetric with respect to their spin axis are possible sources of detectable continuous gravitational waves. This paper presents a fully-coherent search for such signals from eighteen pulsars in data from LIGO and Virgo's third observing run (O3). For known pulsars, efficient and sensitive matched-filter searches can be carried out if one assumes the gravitational radiation is phase-locked to the electromagnetic emission. In the search presented here, we relax this assumption and allow the frequency and frequency time-derivative of the gravitational waves to vary in a small range around those inferred from electromagnetic observations. We find no evidence for continuous gravitational waves, and set upper limits on the strain amplitude for each target. These limits are more constraining for seven of the targets than the spin-down limit defined by ascribing all rotational energy loss to gravitational radiation. In an additional search we look in O3 data for long-duration (hours-months) transient gravitational waves in the aftermath of pulsar glitches for six targets with a total of nine glitches. We report two marginal outliers from this search, but find no clear evidence for such emission either. The resulting duration-dependent strain upper limits do not surpass indirect energy constraints for any of these targets. △ Less

Submitted 27 June, 2022; v1 submitted 21 December, 2021; originally announced December 2021.

Comments: 37 pages, 9 figures, submitted to ApJ

Report number: LIGO-P2100267

Journal ref: ApJ, 932, 133 (2022)

The extragalactic $γ$-ray background: imprints from the physical properties and evolution of star-forming galaxy populations

Authors: Ellis R. Owen, Albert K. H. Kong, Khee-Gan Lee

Abstract: Star-forming galaxies (SFGs) harbour an abundant reservoir of cosmic rays (CRs). At GeV energies, these CRs undergo interactions with their environment to produce $γ$-rays, and the unresolved $γ$-ray emission from populations of SFGs forms a component of the isotropic extragalactic $γ$-ray background (EGB). In this work, we investigate the contribution to the 0.01 - 50 GeV EGB from SFG populations… ▽ More Star-forming galaxies (SFGs) harbour an abundant reservoir of cosmic rays (CRs). At GeV energies, these CRs undergo interactions with their environment to produce $γ$-rays, and the unresolved $γ$-ray emission from populations of SFGs forms a component of the isotropic extragalactic $γ$-ray background (EGB). In this work, we investigate the contribution to the 0.01 - 50 GeV EGB from SFG populations located up to redshift $z=3$. We find this is dominated by starbursts, while the contribution from main sequence SFGs is marginal at all energies. We also demonstrate that most of the $γ$-ray contribution from SFGs emanates from low mass galaxies, with over 80 per cent of the emission originating from galaxies with stellar masses below $10^8 \;\!{\rm M}_{\odot}$. Many of these galaxies are located at relatively high redshift, with their peak EGB contribution arising $\sim$ 700 Myr before the noon of cosmic star-formation. We find that the precise redshift distributions of EGB sources at different energies imprint intensity signatures at different angular scales, which may allow their contribution to be distinguished using analyses of small-scale EGB intensity anisotropies, particularly if the diffuse EGB is dominated by hadronic CR-driven $γ$-ray emission from SFGs. We show that the EGB is sensitive to the evolution of low mass populations of galaxies, particularly around $z\sim2.5$, and that it provides a new means to probe the engagement of CRs in these galaxies before the high noon of cosmic star-formation. △ Less

Submitted 14 April, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: 14 pages, 8 figures, 2 tables. Accepted for publication in MNRAS (v2: updated to match accepted version)

Tests of General Relativity with GWTC-3

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, P. F. de Alarcón, S. Albanesi, R. A. Alfaidi, A. Allocca , et al. (1657 additional authors not shown)

Abstract: The ever-increasing number of detections of gravitational waves (GWs) from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of th… ▽ More The ever-increasing number of detections of gravitational waves (GWs) from compact binaries by the Advanced LIGO and Advanced Virgo detectors allows us to perform ever-more sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We perform a suite of tests of GR using the compact binary signals observed during the second half of the third observing run of those detectors. We restrict our analysis to the 15 confident signals that have false alarm rates $\leq 10^{-3}\, {\rm yr}^{-1}$. In addition to signals consistent with binary black hole (BH) mergers, the new events include GW200115_042309, a signal consistent with a neutron star--BH merger. We find the residual power, after subtracting the best fit waveform from the data for each event, to be consistent with the detector noise. Additionally, we find all the post-Newtonian deformation coefficients to be consistent with the predictions from GR, with an improvement by a factor of ~2 in the -1PN parameter. We also find that the spin-induced quadrupole moments of the binary BH constituents are consistent with those of Kerr BHs in GR. We find no evidence for dispersion of GWs, non-GR modes of polarization, or post-merger echoes in the events that were analyzed. We update the bound on the mass of the graviton, at 90% credibility, to $m_g \leq 1.27 \times 10^{-23} \mathrm{eV}/c^2$. The final mass and final spin as inferred from the pre-merger and post-merger parts of the waveform are consistent with each other. The studies of the properties of the remnant BHs, including deviations of the quasi-normal mode frequencies and damping times, show consistency with the predictions of GR. In addition to considering signals individually, we also combine results from the catalog of GW signals to calculate more precise population constraints. We find no evidence in support of physics beyond GR. △ Less

Submitted 13 December, 2021; originally announced December 2021.

Report number: LIGO-P2100275

All-sky search for gravitational wave emission from scalar boson clouds around spinning black holes in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1647 additional authors not shown)

Abstract: This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust to… ▽ More This paper describes the first all-sky search for long-duration, quasi-monochromatic gravitational-wave signals emitted by ultralight scalar boson clouds around spinning black holes using data from the third observing run of Advanced LIGO. We analyze the frequency range from 20~Hz to 610~Hz, over a small frequency derivative range around zero, and use multiple frequency resolutions to be robust towards possible signal frequency wanderings. Outliers from this search are followed up using two different methods, one more suitable for nearly monochromatic signals, and the other more robust towards frequency fluctuations. We do not find any evidence for such signals and set upper limits on the signal strain amplitude, the most stringent being $\approx10^{-25}$ at around 130~Hz. We interpret these upper limits as both an "exclusion region" in the boson mass/black hole mass plane and the maximum detectable distance for a given boson mass, based on an assumption of the age of the black hole/boson cloud system. △ Less

Submitted 9 May, 2022; v1 submitted 30 November, 2021; originally announced November 2021.

Comments: 28 pages, 16 figures

Report number: P2100343

Journal ref: Phys. Rev. D 105, 102001, 2022

Searches for Gravitational Waves from Known Pulsars at Two Harmonics in the Second and Third LIGO-Virgo Observing Runs

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1672 additional authors not shown)

Abstract: We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both… ▽ More We present a targeted search for continuous gravitational waves (GWs) from 236 pulsars using data from the third observing run of LIGO and Virgo (O3) combined with data from the second observing run (O2). Searches were for emission from the $l=m=2$ mass quadrupole mode with a frequency at only twice the pulsar rotation frequency (single harmonic) and the $l=2, m=1,2$ modes with a frequency of both once and twice the rotation frequency (dual harmonic). No evidence of GWs was found so we present 95\% credible upper limits on the strain amplitudes $h_0$ for the single harmonic search along with limits on the pulsars' mass quadrupole moments $Q_{22}$ and ellipticities $\varepsilon$. Of the pulsars studied, 23 have strain amplitudes that are lower than the limits calculated from their electromagnetically measured spin-down rates. These pulsars include the millisecond pulsars J0437\textminus4715 and J0711\textminus6830 which have spin-down ratios of 0.87 and 0.57 respectively. For nine pulsars, their spin-down limits have been surpassed for the first time. For the Crab and Vela pulsars our limits are factors of $\sim 100$ and $\sim 20$ more constraining than their spin-down limits, respectively. For the dual harmonic searches, new limits are placed on the strain amplitudes $C_{21}$ and $C_{22}$. For 23 pulsars we also present limits on the emission amplitude assuming dipole radiation as predicted by Brans-Dicke theory. △ Less

Submitted 20 July, 2022; v1 submitted 25 November, 2021; originally announced November 2021.

Comments: 37 pages

Report number: LIGO-P2100049

Candidate Tidal Disruption Event AT2019fdr Coincident with a High-Energy Neutrino

Authors: Simeon Reusch, Robert Stein, Marek Kowalski, Sjoert van Velzen, Anna Franckowiak, Cecilia Lunardini, Kohta Murase, Walter Winter, James C. A. Miller-Jones, Mansi M. Kasliwal, Marat Gilfanov, Simone Garrappa, Vaidehi S. Paliya, Tomas Ahumada, Shreya Anand, Cristina Barbarino, Eric C. Bellm, Valery Brinnel, Sara Buson, S. Bradley Cenko, Michael W. Coughlin, Kishalay De, Richard Dekany, Sara Frederick, Avishay Gal-Yam , et al. (24 additional authors not shown)

Abstract: The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flar… ▽ More The origins of the high-energy cosmic neutrino flux remain largely unknown. Recently, one high-energy neutrino was associated with a tidal disruption event (TDE). Here we present AT2019fdr, an exceptionally luminous TDE candidate, coincident with another high-energy neutrino. Our observations, including a bright dust echo and soft late-time X-ray emission, further support a TDE origin of this flare. The probability of finding two such bright events by chance is just 0.034%. We evaluate several models for neutrino production and show that AT2019fdr is capable of producing the observed high-energy neutrino, reinforcing the case for TDEs as neutrino sources. △ Less

Submitted 10 June, 2022; v1 submitted 17 November, 2021; originally announced November 2021.

Comments: 20 pages, 6 figures, 6 tables

Journal ref: Phys. Rev. Lett. 128, 221101 (2022)

The population of merging compact binaries inferred using gravitational waves through GWTC-3

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato , et al. (1612 additional authors not shown)

Abstract: We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8… ▽ More We report on the population properties of 76 compact binary mergers detected with gravitational waves below a false alarm rate of 1 per year through GWTC-3. The catalog contains three classes of binary mergers: BBH, BNS, and NSBH mergers. We infer the BNS merger rate to be between 10 $\rm{Gpc^{-3} yr^{-1}}$ and 1700 $\rm{Gpc^{-3} yr^{-1}}$ and the NSBH merger rate to be between 7.8 $\rm{Gpc^{-3}\, yr^{-1}}$ and 140 $\rm{Gpc^{-3} yr^{-1}}$ , assuming a constant rate density versus comoving volume and taking the union of 90% credible intervals for methods used in this work. Accounting for the BBH merger rate to evolve with redshift, we find the BBH merger rate to be between 17.9 $\rm{Gpc^{-3}\, yr^{-1}}$ and 44 $\rm{Gpc^{-3}\, yr^{-1}}$ at a fiducial redshift (z=0.2). We obtain a broad neutron star mass distribution extending from $1.2^{+0.1}_{-0.2} M_\odot$ to $2.0^{+0.3}_{-0.3} M_\odot$. We can confidently identify a rapid decrease in merger rate versus component mass between neutron star-like masses and black-hole-like masses, but there is no evidence that the merger rate increases again before 10 $M_\odot$. We also find the BBH mass distribution has localized over- and under-densities relative to a power law distribution. While we continue to find the mass distribution of a binary's more massive component strongly decreases as a function of primary mass, we observe no evidence of a strongly suppressed merger rate above $\sim 60 M_\odot$. The rate of BBH mergers is observed to increase with redshift at a rate proportional to $(1+z)^κ$ with $κ= 2.9^{+1.7}_{-1.8}$ for $z\lesssim 1$. Observed black hole spins are small, with half of spin magnitudes below $χ_i \simeq 0.25$. We observe evidence of negative aligned spins in the population, and an increase in spin magnitude for systems with more unequal mass ratio. △ Less

Submitted 23 February, 2022; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: v2: minor edits, most to Table 1 and caption; v3: rerun with public data; Data release: https://zenodo.org/record/5655785; v4: update Fig 14

Report number: LIGO-P2100239

Search for Gravitational Waves Associated with Gamma-Ray Bursts Detected by Fermi and Swift During the LIGO-Virgo Run O3b

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato , et al. (1610 additional authors not shown)

Abstract: We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target bina… ▽ More We search for gravitational-wave signals associated with gamma-ray bursts detected by the Fermi and Swift satellites during the second half of the third observing run of Advanced LIGO and Advanced Virgo (1 November 2019 15:00 UTC-27 March 2020 17:00 UTC).We conduct two independent searches: a generic gravitational-wave transients search to analyze 86 gamma-ray bursts and an analysis to target binary mergers with at least one neutron star as short gamma-ray burst progenitors for 17 events. We find no significant evidence for gravitational-wave signals associated with any of these gamma-ray bursts. A weighted binomial test of the combined results finds no evidence for sub-threshold gravitational wave signals associated with this GRB ensemble either. We use several source types and signal morphologies during the searches, resulting in lower bounds on the estimated distance to each gamma-ray burst. Finally, we constrain the population of low luminosity short gamma-ray bursts using results from the first to the third observing runs of Advanced LIGO and Advanced Virgo. The resulting population is in accordance with the local binary neutron star merger rate. △ Less

Submitted 5 November, 2021; originally announced November 2021.

Comments: 26 pages, 6 figures, 2 tables

Report number: P2100091

GWTC-3: Compact Binary Coalescences Observed by LIGO and Virgo During the Second Part of the Third Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, S. Akcay, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin , et al. (1637 additional authors not shown)

Abstract: The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There ar… ▽ More The third Gravitational-Wave Transient Catalog (GWTC-3) describes signals detected with Advanced LIGO and Advanced Virgo up to the end of their third observing run. Updating the previous GWTC-2.1, we present candidate gravitational waves from compact binary coalescences during the second half of the third observing run (O3b) between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. There are 35 compact binary coalescence candidates identified by at least one of our search algorithms with a probability of astrophysical origin $p_\mathrm{astro} > 0.5$. Of these, 18 were previously reported as low-latency public alerts, and 17 are reported here for the first time. Based upon estimates for the component masses, our O3b candidates with $p_\mathrm{astro} > 0.5$ are consistent with gravitational-wave signals from binary black holes or neutron star-black hole binaries, and we identify none from binary neutron stars. However, from the gravitational-wave data alone, we are not able to measure matter effects that distinguish whether the binary components are neutron stars or black holes. The range of inferred component masses is similar to that found with previous catalogs, but the O3b candidates include the first confident observations of neutron star-black hole binaries. Including the 35 candidates from O3b in addition to those from GWTC-2.1, GWTC-3 contains 90 candidates found by our analysis with $p_\mathrm{astro} > 0.5$ across the first three observing runs. These observations of compact binary coalescences present an unprecedented view of the properties of black holes and neutron stars. △ Less

Submitted 23 October, 2023; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: 88 pages (10 pages author list, 31 pages main text, 1 page acknowledgements, 24 pages appendices, 22 pages bibliography), 17 figures, 16 tables. Update to match version to be published in Physical Review X. Data products available from https://gwosc.org/GWTC-3/

Report number: LIGO-P2000318

Journal ref: Phys. Rev. X; 13(4):041039; 2023

Constraints on the cosmic expansion history from GWTC-3

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, A. Allocca, P. A. Altin , et al. (1654 additional authors not shown)

Abstract: We use 47 gravitational-wave sources from the Third LIGO-Virgo-KAGRA Gravitational-Wave Transient Catalog (GWTC-3) to estimate the Hubble parameter $H(z)$, including its current value, the Hubble constant $H_0$. Each gravitational-wave (GW) signal provides the luminosity distance to the source and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog.… ▽ More We use 47 gravitational-wave sources from the Third LIGO-Virgo-KAGRA Gravitational-Wave Transient Catalog (GWTC-3) to estimate the Hubble parameter $H(z)$, including its current value, the Hubble constant $H_0$. Each gravitational-wave (GW) signal provides the luminosity distance to the source and we estimate the corresponding redshift using two methods: the redshifted masses and a galaxy catalog. Using the binary black hole (BBH) redshifted masses, we simultaneously infer the source mass distribution and $H(z)$. The source mass distribution displays a peak around $34\, {\rm M_\odot}$, followed by a drop-off. Assuming this mass scale does not evolve with redshift results in a $H(z)$ measurement, yielding $H_0=68^{+12}_{-7} {\rm km\,s^{-1}\,Mpc^{-1}}$ ($68\%$ credible interval) when combined with the $H_0$ measurement from GW170817 and its electromagnetic counterpart. This represents an improvement of 17% with respect to the $H_0$ estimate from GWTC-1. The second method associates each GW event with its probable host galaxy in the catalog GLADE+, statistically marginalizing over the redshifts of each event's potential hosts. Assuming a fixed BBH population, we estimate a value of $H_0=68^{+8}_{-6} {\rm km\,s^{-1}\,Mpc^{-1}}$ with the galaxy catalog method, an improvement of 42% with respect to our GWTC-1 result and 20% with respect to recent $H_0$ studies using GWTC-2 events. However, we show that this result is strongly impacted by assumptions about the BBH source mass distribution; the only event which is not strongly impacted by such assumptions (and is thus informative about $H_0$) is the well-localized event GW190814. △ Less

Submitted 19 November, 2021; v1 submitted 5 November, 2021; originally announced November 2021.

Comments: Main paper: 30 pages, 15 figure, 7 tables

Report number: LIGO-P2100185-v6

Investigation of the timing and spectral properties of an Ultra-luminous X-ray pulsar NGC 7793 P13

Authors: Lupin Chun-Che Lin, Chin-Ping Hu, Jumpei Takata, Kwan-Lok Li, C. Y. Hui, A. K. H. Kong

Abstract: We perform both timing and spectral analyses using the archival X-ray data taken with Swift, XMM-Newton, NICER, and NuSTAR from 2016 to 2020 to study an ultra-luminous pulsar, NGC 7793 P13, that showed a long period of super-Eddington accretion. We use the Rayleigh test to investigate the pulsation at different epochs, and confirm the variation of the pulse profile with the finite Gaussian mixture… ▽ More We perform both timing and spectral analyses using the archival X-ray data taken with Swift, XMM-Newton, NICER, and NuSTAR from 2016 to 2020 to study an ultra-luminous pulsar, NGC 7793 P13, that showed a long period of super-Eddington accretion. We use the Rayleigh test to investigate the pulsation at different epochs, and confirm the variation of the pulse profile with the finite Gaussian mixture modelling and two-sample Kuiper test. Taking into account the periodic variation of the spin periods caused by the orbital Doppler effect, we further determine an orbital period of ~65 d and show that no significant correlation can be detected between the orbital phase and the pulsed fraction. The pulsed spectrum of NGC7793 P13 in 0.5-20keV can be simply described using a power-law with a high energy exponential cutoff, while the broad-band phase-averaged spectrum of the same energy range requires two additional components to account for the contribution of a thermal accretion disk and the Comptonization photons scattered into the hard X-rays. We find that NGC 7793 P13 stayed in the hard ultra-luminous state and the pulsed spectrum was relatively soft when the source was faint in the end of 2019. Moreover, an absorption feature close to 1.3 keV is marginally detected from the pulsed spectra and it is possibly associated with a cyclotron resonant scattering feature. △ Less

Submitted 3 November, 2021; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: 13 pages, 6 figures and 5 tables, accepted by ApJ

The HASHTAG project: The First Submillimeter Images of the Andromeda Galaxy from the Ground

Authors: Matthew W. L. Smith, Stephen A. Eales, Thomas G. Williams, Bumhyun Lee, Zongnan Li, Pauline Barmby, Martin Bureau, Scott Chapman, Brian S. Cho, Aeree Chung, Eun Jung Chung, Hui-Hsuan Chung, Christopher J. R. Clark, David L. Clements, Timothy A. Davis, Ilse De Looze, David J. Eden, Gayathri Athikkat-Eknath, George P. Ford, Yu Gao, Walter Gear, Haley L. Gomez, Richard de Grijs, Jinhua He, Luis C. Ho , et al. (24 additional authors not shown)

Abstract: Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing… ▽ More Observing nearby galaxies with submillimeter telescopes on the ground has two major challenges. First, the brightness is significantly reduced at long submillimeter wavelengths compared to the brightness at the peak of the dust emission. Second, it is necessary to use a high-pass spatial filter to remove atmospheric noise on large angular scales, which has the unwelcome by-product of also removing the galaxy's large-scale structure. We have developed a technique for producing high-resolution submillimeter images of galaxies of large angular size by using the telescope on the ground to determine the small-scale structure (the large Fourier components) and a space telescope (Herschel or Planck) to determine the large-scale structure (the small Fourier components). Using this technique, we are carrying out the HARP and SCUBA-2 High Resolution Terahertz Andromeda Galaxy Survey (HASHTAG), an international Large Program on the James Clerk Maxwell Telescope, with one aim being to produce the first high-fidelity high-resolution submillimeter images of Andromeda. In this paper, we describe the survey, the method we have developed for combining the space-based and ground-based data, and present the first HASHTAG images of Andromeda at 450 and 850um. We also have created a method to predict the CO(J=3-2) line flux across M31, which contaminates the 850um band. We find that while normally the contamination is below our sensitivity limit, the contamination can be significant (up to 28%) in a few of the brightest regions of the 10 kpc ring. We therefore also provide images with the predicted line emission removed. △ Less

Submitted 30 September, 2021; originally announced October 2021.

Comments: 26 pages, 19 figures. Submitted to ApJS June 2021, Accepted September 2021

A Bayesian Inference Framework for Gamma-Ray Burst Afterglow Properties

Authors: En-Tzu Lin, Fergus Hayes, Gavin P. Lamb, Ik Siong Heng, Albert K. H. Kong, Michael J. Williams, Surojit Saha, John Veitch

Abstract: In the field of multi-messenger astronomy, Bayesian inference is commonly adopted to compare the compatibility of models given the observed data. However, to describe a physical system like neutron star mergers and their associated gamma-ray burst (GRB) events, usually more than ten physical parameters are incorporated in the model. With such a complex model, likelihood evaluation for each Monte C… ▽ More In the field of multi-messenger astronomy, Bayesian inference is commonly adopted to compare the compatibility of models given the observed data. However, to describe a physical system like neutron star mergers and their associated gamma-ray burst (GRB) events, usually more than ten physical parameters are incorporated in the model. With such a complex model, likelihood evaluation for each Monte Carlo sampling point becomes a massive task and requires a significant amount of computational power. In this work, we perform quick parameter estimation on simulated GRB X-ray light curves using an interpolated physical GRB model. This is achieved by generating a grid of GRB afterglow light curves across the parameter space and replacing the likelihood with a simple interpolation function in the high-dimensional grid that stores all light curves. This framework, compared to the original method, leads to a $\sim$90$\times$ speedup per likelihood estimation. It will allow us to explore different jet models and enable fast model comparison in the future. △ Less

Submitted 30 September, 2021; originally announced September 2021.

Comments: 9 pages, 4 figures, accepted to the special issue of Universe, "Waiting for GODOT -- Present and Future of Multi-Messenger Astronomy"

Search for subsolar-mass binaries in the first half of Advanced LIGO and Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato , et al. (1612 additional authors not shown)

Abstract: We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio… ▽ More We report on a search for compact binary coalescences where at least one binary component has a mass between 0.2 $M_\odot$ and 1.0 $M_\odot$ in Advanced LIGO and Advanced Virgo data collected between 1 April 2019 1500 UTC and 1 October 2019 1500 UTC. We extend previous analyses in two main ways: we include data from the Virgo detector and we allow for more unequal mass systems, with mass ratio $q \geq 0.1$. We do not report any gravitational-wave candidates. The most significant trigger has a false alarm rate of 0.14 $\mathrm{yr}^{-1}$. This implies an upper limit on the merger rate of subsolar binaries in the range $[220-24200] \mathrm{Gpc}^{-3} \mathrm{yr}^{-1}$, depending on the chirp mass of the binary. We use this upper limit to derive astrophysical constraints on two phenomenological models that could produce subsolar-mass compact objects. One is an isotropic distribution of equal-mass primordial black holes. Using this model, we find that the fraction of dark matter in primordial black holes is $f_\mathrm{PBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 6\%$. The other is a dissipative dark matter model, in which fermionic dark matter can collapse and form black holes. The upper limit on the fraction of dark matter black holes depends on the minimum mass of the black holes that can be formed: the most constraining result is obtained at $M_\mathrm{min}=1 M_\odot$, where $f_\mathrm{DBH} \equiv Ω_\mathrm{PBH} / Ω_\mathrm{DM} \lesssim 0.003\%$. These are the tightest limits on spinning subsolar-mass binaries to date. △ Less

Submitted 24 September, 2021; originally announced September 2021.

Report number: LIGO-P2100163-v8

Search for continuous gravitational waves from 20 accreting millisecond X-ray pulsars in O3 LIGO data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, T. D. Abbott, F. Acernese, K. Ackley, C. Adams, N. Adhikari, R. X. Adhikari, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, K. Agatsuma, N. Aggarwal, O. D. Aguiar, L. Aiello, A. Ain, T. Akutsu, S. Albanesi, A. Allocca, P. A. Altin, A. Amato, C. Anand , et al. (1612 additional authors not shown)

Abstract: Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an… ▽ More Results are presented of searches for continuous gravitational waves from 20 accreting millisecond X-ray pulsars with accurately measured spin frequencies and orbital parameters, using data from the third observing run of the Advanced LIGO and Advanced Virgo detectors. The search algorithm uses a hidden Markov model, where the transition probabilities allow the frequency to wander according to an unbiased random walk, while the $\mathcal{J}$-statistic maximum-likelihood matched filter tracks the binary orbital phase. Three narrow sub-bands are searched for each target, centered on harmonics of the measured spin frequency. The search yields 16 candidates, consistent with a false alarm probability of 30% per sub-band and target searched. These candidates, along with one candidate from an additional target-of-opportunity search done for SAX J1808.4$-$3658, which was in outburst during one month of the observing run, cannot be confidently associated with a known noise source. Additional follow-up does not provide convincing evidence that any are a true astrophysical signal. When all candidates are assumed non-astrophysical, upper limits are set on the maximum wave strain detectable at 95% confidence, $h_0^{95\%}$. The strictest constraint is $h_0^{95\%} = 4.7\times 10^{-26}$ from IGR J17062$-$6143. Constraints on the detectable wave strain from each target lead to constraints on neutron star ellipticity and $r$-mode amplitude, the strictest of which are $ε^{95\%} = 3.1\times 10^{-7}$ and $α^{95\%} = 1.8\times 10^{-5}$ respectively. This analysis is the most comprehensive and sensitive search of continuous gravitational waves from accreting millisecond X-ray pulsars to date. △ Less

Submitted 21 January, 2022; v1 submitted 19 September, 2021; originally announced September 2021.

Comments: 40 pages, 6 figures. This version contains minor typographical revisions to match published article

Report number: LIGO-P2100221

Journal ref: Phys. Rev. D 105, 022002 (2022)

Inclination estimates from off-axis GRB afterglow modelling

Authors: Gavin P Lamb, Joseph J Fernández, Fergus Hayes, Albert K H Kong, En-Tzu Lin, Nial R Tanvir, Martin Hendry, Ik Siong Heng, Surojit Saha, John Veitch

Abstract: For gravitational wave (GW) detected neutron star mergers, one of the leading candidates for electromagnetic (EM) counterparts is the afterglow from an ultra-relativistic jet. Where this afterglow is observed, it will likely be viewed off-axis, such as the afterglow following GW170817/GRB 170817A. The temporal behaviour of an off-axis observed GRB afterglow can be used to reveal the lateral jet st… ▽ More For gravitational wave (GW) detected neutron star mergers, one of the leading candidates for electromagnetic (EM) counterparts is the afterglow from an ultra-relativistic jet. Where this afterglow is observed, it will likely be viewed off-axis, such as the afterglow following GW170817/GRB 170817A. The temporal behaviour of an off-axis observed GRB afterglow can be used to reveal the lateral jet structure, and statistical model fits can put constraints on the various model free-parameters. Amongst these parameters is the inclination of the system to the line of sight. Along with the GW detection, the afterglow modelling provides the best constraint on the inclination to the line-of-sight and can improve the estimates of cosmological parameters e.g. the Hubble constant, from GW-EM events. However, modelling of the afterglow depends on the assumed jet structure and, often overlooked, the effects of lateral spreading. Here we show how the inclusion of lateral spreading in the afterglow models can affect the estimated inclination of GW-EM events. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: 14 pages, 5 figures, Accepted to the special issue of Universe, "Waiting for GODOT -- Present and Future of Multi-Messenger Astronomy"