Sky Position Research Articles

State-space analysis of a continuous gravitational wave source with a pulsar timing array: inclusion of the pulsar terms

ABSTRACT Pulsar timing arrays (PTA) can detect continuous nanohertz gravitational waves (GW) emitted by individual supermassive black hole binaries. The data analysis procedure can be formulated within a time-domain, state-space framework, in which the radio timing observations are related to a temporal sequence of latent states, namely the intrinsic pulsar spin frequency. The achromatic wandering of the pulsar spin frequency is tracked using a Kalman filter concurrently with the pulse frequency modulation induced by a GW from a single source. The modulation is the sum of terms proportional to the GW strain at the Earth and at every pulsar in the array. Here, we generalize previous state-space formulations of the PTA problem to include the pulsar terms; that is, we copy the pulsar terms from traditional, non-state-space analyses over to the state-space framework. The performance of the generalized Kalman filter is tested using astrophysically representative software injections in Gaussian measurement noise. It is shown that including the pulsar terms corrects for previously identified biases in the parameter estimates (especially the sky position of the source) which also arise in traditional matched-filter analyses that exclude the pulsar terms. Additionally, including the pulsar terms decreases the minimum detectable strain by 14 per cent. Overall, the study verifies that the pulsar terms do not raise any special extra impediments for the state-space framework, beyond those studied in traditional analyses. The inspiral-driven evolution of the wave frequency at the Earth and at the retarded time at every pulsar in the array is also investigated.

Kalman tracking and parameter estimation of continuous gravitational waves with a pulsar timing array

Abstract Continuous nanohertz gravitational waves from individual supermassive black hole binaries may be detectable with pulsar timing arrays. A novel search strategy is developed, wherein intrinsic achromatic spin wandering is tracked simultaneously with the modulation induced by a single gravitational wave source in the pulse times of arrival. A two-step inference procedure is applied within a state-space framework, such that the modulation is tracked with a Kalman filter, which then provides a likelihood for nested sampling. The procedure estimates the static parameters in the problem, such as the sky position of the source, without fitting for ensemble-averaged statistics such as the power spectral density of the timing noise, and therefore complements traditional parameter estimation methods. It also returns the Bayes factor relating a model with a single gravitational wave source to one without, complementing traditional detection methods. It is shown via astrophysically representative software injections in Gaussian measurement noise that the procedure distinguishes a gravitational wave from pure noise down to a characteristic wave strain of h0 ≈ 2 × 10−15. Full posterior distributions of model parameters are recovered and tested for accuracy. There is a bias of ≈0.3 rad in the marginalised one-dimensional posterior for the orbital inclination ι, introduced by dropping the so-called ‘pulsar terms’. Smaller biases $\lesssim 10~{{\%}}$ are also observed in other static parameters.

The Chandra Source Catalog Release 2 Series

The Chandra Source Catalog (CSC) is a virtual X-ray astrophysics facility that enables both detailed individual source studies and statistical studies of large samples of X-ray sources detected in Advanced CCD Imaging Spectrometer and High Resolution Camera-I imaging observations obtained by the Chandra X-ray Observatory. The catalog provides carefully curated, high-quality, and uniformly calibrated and analyzed tabulated positional, spatial, photometric, spectral, and temporal source properties, as well as science-ready X-ray data products. The latter includes multiple types of source- and field-based FITS format products that can be used as a basis for further research, significantly simplifying follow-up analysis of scientifically meaningful source samples. We discuss in detail the algorithms used for the CSC Release 2 Series, including CSC 2.0, which includes 317,167 unique X-ray sources on the sky identified in observations released publicly through the end of 2014, and CSC 2.1, which adds Chandra data released through the end of 2021 and expands the catalog to 407,806 sources. Besides adding more recent observations, the CSC Release 2 Series includes multiple algorithmic enhancements that provide significant improvements over earlier releases. The compact source sensitivity limit for most observations is ∼5 photons over most of the field of view, which is ∼2× fainter than Release 1, achieved by coadding observations and using an optimized source detection approach. A Bayesian X-ray aperture photometry code produces robust fluxes even in crowded fields and for low-count sources. The current release, CSC 2.1, is tied to the Gaia-CRF3 astrometric reference frame for the best sky positions for catalog sources.

Comparing Gaia, NED and SIMBAD source classifications in nearby galaxies

Abstract Gaia Data Release 3 (DR3) provides the first classifications for the sources in Gaia’s all-sky database. Most Gaia sources are stars in the Milky Way, but DR3 also contains many sources that belong to nearby galaxies, as well as background galaxies and quasars. In this work, we compare the Gaia classifications from the Discrete Source Classifier (CU8-DSC) module to the more detailed and heterogeneous classifications in NED and/or SIMBAD for sources with sky positions within twice the Holmberg radius of nearby galaxies. Matching these catalogues gives approximately 3.2 × 105 unique Gaia matches for 4 × 105 sources over 1040 galaxies (excluding some large Local Group galaxies) in the Local Volume Galaxy catalogue. Matched sources contain a lower fraction of Gaia-classified stars and higher fractions of galaxies and quasars (∼95, 2 and 2 per cent, respectively) than DR3 overall. Considering NED (SIMBAD) classifications as truth values, the balanced accuracy of Gaia classification is 0.80 (0.83): the most common disagreements are literature-classified galaxies Gaia-classified as stars and literature-classified stars Gaia-classified as quasars. Purity (P) and completeness (C) metrics show that agreement between Gaia classification and NED/SIMBAD classification is best for stars (P, C ∼ 0.9), and decreases for quasars (P < 0.3, 0.7 < C < 0.8), galaxies (0.7 < P < 0.8, 0.3 < C < 0.6), white dwarfs (0.04 < P < 0.6, C ∼ 0.6), and binary stars (P, C < 0.1). NED or SIMBAD sources classified only by detection wavelength are most often Gaia-classified as stars, while non-stellar components of galaxies appear in all Gaia classes.

Positive and unlabelled machine learning reveals new fast radio burst repeater candidates

ABSTRACT Fast radio bursts (FRBs) are astronomical radio transients of unknown origin. A minority of FRBs have been observed to originate from repeating sources, and it is unknown which apparent one-off bursts are hidden repeaters. Recent studies increasingly suggest that there are intrinsic physical differences between repeating and non-repeating FRBs. Previous research has used machine learning classification techniques to identify apparent non-repeaters with repeater characteristics, whose sky positions would be ideal targets for future observation campaigns. However, these methods have not sufficiently accounted for the positive and unlabelled (PU) nature of the data, wherein true labels are only available for repeaters. Modified techniques that do not inadvertently learn properties of hidden repeaters as characteristic of non-repeaters are likely to identify additional repeater candidates with greater accuracy. We present in this paper the first known attempt at applying PU-specific machine learning techniques to study FRBs. We train an ensemble of five PU-specific classifiers on the available data and use them to identify 66 repeater candidates in burst data from the CHIME/FRB collaboration, 18 of which were not identified with the use of machine learning classifiers in past research. Our results additionally support repeaters and non-repeaters having intrinsically different physical properties, particularly spectral index, frequency width, and burst width. This work additionally opens new possibilities to study repeating and non-repeating FRBs using the framework of PU learning.

Tessilator: a one-stop shop for measuring TESS rotation periods

ABSTRACT We present a software package designed to produce photometric light curves and measure rotation periods from full-frame images taken by the Transiting Exoplanet Survey Satellite (TESS), which we name ‘tessilator’. tessilator is the only publicly available code that will run a full light curve and rotation period ($P_{\rm rot}$) analysis based on just a (list of) target identifier(s) or sky position(s) via a simple command-line prompt. This paper sets out to introduce the rationale for developing tessilator, and then describes the methods, considerations, and assumptions for: extracting photometry; dealing with potential contamination; accounting for natural and instrumental systematic effects; light-curve normalization and detrending; removing outliers and unreliable data; and finally, measuring the $P_{\rm rot}$ value and several periodogram attributes. Our methods have been tuned specifically to optimize TESS light curves and are independent from the pipelines developed by the TESS Science Processing Operations Center (SPOC), meaning tessilator can, in principle, analyse any target across the entire celestial sphere. We compare tessilator$P_{\rm rot}$ measurements with TESS-SPOC-derived light curves of 1560 (mainly solar-type and low-mass) stars across four benchmark open clusters (Pisces–Eridanus, the Pleiades, the Hyades, and Praesepe) and a sample of nearby field M-dwarfs. From a vetted subsample of 864 targets we find an excellent return of $P_{\rm rot}$ matches for the first three open clusters ($\gt 85$ per cent) and a moderate ($\sim 60$ per cent) match for the 700 Myr Praesepe and MEarth sample, which validates tessilator as a tool for measuring $P_{\rm rot}$. The tessilator code is available at https://github.com/alexbinks/tessilator.

Calibrating the global network of gravitational wave observatories via laser power calibration at NIST and PTB

Current gravitational wave (GW) observatories rely on photon calibrators that use laser radiation pressure to generate displacement fiducials used to calibrate detector output signals. Reducing calibration uncertainty enables optimal extraction of astrophysical information such as source distance and sky position from detected signals. For the ongoing O4 observing run that started on 24 May 2023, the global GW detector network is employing a new calibration scheme with transfer standards calibrated at both the National Institute of Standards and Technology (NIST) and the Physikalisch-Technische Bundesanstalt (PTB). These transfer standards will circulate between the observatories and the metrology institutes to provide laser power calibration traceable to the International System of Units (SI) and enable assessment and reduction of relative calibration errors for the observatory network. The Laser Interferometer Gravitational-Wave Observatory (LIGO) project and the Virgo project are currently participating in the new calibration scheme. The Large-scale Cryogenic Gravitational-wave Telescope project (KAGRA) is expected to join in 2024, with the LIGO Aundha Observatory in India joining later. Before implementing this new scheme, a NIST-PTB bilateral comparison was conducted. It validated the scale representation by both laboratories, with a degree of equivalence of −0.2% and an associated expanded uncertainty of 0.32% (k = 2) which is significantly lower than previous studies. We describe the transfer of power sensor calibration, including detailed uncertainty estimates, from the transfer standards calibrated by NIST and PTB to the sensors operating continuously at the interferometer end stations. Finally, we discuss the ongoing calibration of Pcal-induced displacement fiducials for the O4 observing run. Achieved combined standard uncertainty levels as low as 0.3% facilitate calibrating the interferometer output signals with sub-percent accuracy.

High-priority targets for transient gravitational waves from glitching pulsars

ABSTRACT Glitching pulsars are expected to be important sources of gravitational waves (GWs). In this paper, we explore six different models that propose the emission of transient continuous waves, lasting days to months, coincident with glitches. The maximal GW energy is calculated for each model, which is then used to determine whether associated GWs could be detectable with LIGO-Virgo-KAGRA’s O4 detectors. We provide an analytical approximation to calculate the signal-to-noise ratio (SNR) which includes information about the source’s sky position, improving on previous estimates that assume isotropic or sky and orientation averaged sensitivities. By analysing the entire glitching population, we find that certain models predict detectable signals in O4, whereas others do not. We also rank glitching pulsars by SNR, based on archival data, and we find that for all models, the Vela pulsar (PSR J0835$-$4510) would provide the strongest signal. Moreover, PSR J0537$-$6910 is not expected to yield a detectable signal in O4, but will start becoming relevant for next-generation detectors. Our analysis also extends to the entire pulsar population, regardless of whether they have glitched, and we provide a list of pulsars that would present a significant signal, if they were to glitch. Finally, we apply our analysis to the 2024 April Vela glitch and find that a signal should be detectable under certain models. The non-detection of a supposedly detectable signal would provide an efficiency factor that quantifies a model’s contribution to GW emission, eventually leading to a differentiation of models and independent constraints on physical parameters.

Improved source classification and performance analysis using Gaia DR3

The Discrete Source Classifier (DSC) provides probabilistic classification of sources in Gaia Data Release 3 (GDR3) using a Bayesian framework and a global prior. The DSC Combmod classifier in GDR3 achieved for the extragalactic classes (quasars and galaxies) a high completeness of 92%, but a low purity of 22% (all sky, all magnitudes) due to contamination from the far larger star class. However, these single metrics mask significant variation in performance with magnitude and sky position. Furthermore, a better combination of the individual classifiers that comprise Combmod is possible. Here we compute two-dimensional (2D) representations of the completeness and the purity as a function of Galactic latitude and source brightness, and also exclude the Magellanic Clouds where stellar contamination significantly reduces the purity. Reevaluated on a cleaner validation set and without introducing changes to the published GDR3 DSC probabilities themselves, we here achieve for Combmod average 2D completenesses of 92% and 95% and average 2D purities of 55% and 89% for the quasar and galaxy classes, respectively. Since the relative proportions of extragalactic objects to stars in Gaia is expected to vary significantly with brightness and latitude, we then introduce a new prior that is a continuous function of brightness and latitude, and compute new class probabilities from the GDR3 DSC component classifiers, Specmod and Allosmod. Contrary to expectations, this variable prior only improves the performance by a few percentage points, mostly at the faint end. Significant improvement, however, is obtained by a new additive combination of Specmod and Allosmod. This classifier, Combmod-α, achieves average 2D completenesses of 82% and 93% and average 2D purities of 79% and 93% for the quasar and galaxy classes, respectively, when using the global prior. Thus, we achieve a significant improvement in purity for a small loss of completeness. The improvement is most significant for faint quasars (G≥20) where the purity rises from 20% to 62%.

X-Ray Quasi-periodic Eruptions and Tidal Disruption Events Prefer Similar Host Galaxies

In the past 5 yr, six X-ray quasi-periodic eruption (QPE) sources have been discovered in the nuclei of nearby galaxies. Their origin remains an open question. We present Multi Unit Spectroscopic Explorer integral field spectroscopy of five QPE host galaxies to characterize their properties. We find that 3/5 galaxies host extended emission-line regions (EELRs) up to 10 kpc in size. The EELRs are photoionized by a nonstellar continuum, but the current nuclear luminosity is insufficient to power the observed emission lines. The EELRs are decoupled from the stars both kinematically and in projected sky position, and the low velocities and velocity dispersions (<100 km s−1 and ≲75 km s−1, respectively) are inconsistent with being driven by active galactic nuclei (AGNs) or shocks. The origin of the EELRs is likely a previous phase of nuclear activity. QPE host galaxies share several similarities with tidal disruption event (TDE) hosts, including an overrepresentation of galaxies with strong Balmer absorption and little ongoing star formation, as well as a preference for a short-lived (the typical EELR lifetime is ∼15,000 yr), gas-rich phase where the nucleus has recently faded significantly. This suggests that QPEs and TDEs may share a common formation channel, disfavoring AGN accretion disk instabilities as the origin of QPEs. If QPEs are related to extreme mass ratio inspiral systems (EMRIs), e.g., stellar-mass objects on bound orbits about massive black holes, the high incidence of EELRs and recently faded nuclei could be used to localize the hosts of EMRIs discovered by low-frequency gravitational-wave observatories.

The Massive and Distant Clusters of WISE Survey 2: Equatorial First Data Release

The Massive and Distant Clusters of WISE Survey 2 (MaDCoWS2) is a new survey designed as the successor of the original MaDCoWS survey. MaDCoWS2 improves upon its predecessor by using deeper optical and infrared data and a more powerful detection algorithm (PZWav). As input to the search, we use grz photometry from the DECam Legacy Survey (DECaLS) in combination with W1 and W2 photometry from the CatWISE2020 catalog to derive the photometric redshifts with full redshift probability distribution functions for Wide-field Infrared Survey Explorer (WISE)-selected galaxies. Cluster candidates are then detected using the PZWav algorithm to find three-dimensional galaxy overdensities from the sky positions and photometric redshifts. This paper provides the first MaDCoWS2 data release, covering 1461 (1838 without masking) deg2 centered on the Hyper-SuprimeCam Subaru Strategic Program equatorial fields. Within this region, we derive a catalog of 22,970 galaxy cluster candidates detected at a signal-to-noise ratio (S/N) > 5. These clusters span the redshift range 0.1 < z < 2, including 1312 candidates at z > 1.5. We compare MaDCoWS2 to six existing catalogs in the area. We rediscover 60%–92% of the clusters in these surveys at S/N > 5. The medians of the absolute redshift offset are <0.02 relative to these surveys, while the standard deviations are less than 0.06. The median offsets between the detection position from MaDCoWS2 and other surveys are less than 0.25 Mpc. We quantify the relation between S/N and gas mass, total mass, luminosity, and richness from other surveys using a redshift-dependent power law relation. We find that the S/N-richness relation exhibits the lowest scatter.

Searching for late-time interaction signatures in Type Ia supernovae from the Zwicky Transient Facility

The nature of the progenitor systems and explosion mechanisms that give rise to Type Ia supernovae (SNe Ia) are still debated. The interaction signature of circumstellar material (CSM) being swept up by the expanding ejecta can constrain the type of system from which it was ejected. However, most previous studies have focussed on finding CSM ejected shortly before the SN Ia explosion, which still resides close to the explosion site resulting in short delay times until the interaction starts. We used a sample of 3\,627 SNe Ia from the Zwicky Transient Facility (ZTF) that were discovered between 2018 and 2020 and searched for interaction signatures greater than 100 days after peak brightness. By binning the late-time light curve data to push the detection limit as deep as possible, we identified potential late-time rebrightening in three SNe Ia (SN 2018grt, SN 2019dlf, and SN 2020tfc). The late-time optical detections occur between 550 and 1450\,d after peak brightness, have mean absolute r -band magnitudes of $-$16.4 to $-$16.8 mag, and last up to a few hundred days, which is significantly brighter than the late-time CSM interaction discovered in the prototype, SN 2015cp. The late-time detections in the three objects all occur within 0.8 kpc of the host nucleus and are not easily explained by nuclear activity, another transient at a similar sky position, or data quality issues. This is suggestive of environment or specific progenitor characteristics playing a role in the production of potential CSM signatures in these SNe Ia. Through simulating the ZTF survey, we estimate that $&lt;$0.5 per cent of normal SNe Ia display a late-time ($&gt;100$ d post peak) strong CSM interaction. This is equivalent to an absolute rate of $ to $54_ $ Gpc$^ $ yr$^ $ assuming a constant SN Ia rate of $2.4 $ Mpc$^ $ yr$^ $ for $z 0.1$. Weaker interaction signatures of emission, more similar to the strength seen in SN 2015cp, could be more common but are difficult to constrain with our survey depth.

The NANOGrav 12.5 yr Data Set: Search for Gravitational Wave Memory

We present the results of a Bayesian search for gravitational wave (GW) memory in the NANOGrav 12.5 yr data set. We find no convincing evidence for any gravitational wave memory signals in this data set. We find a Bayes factor of 2.8 in favor of a model that includes a memory signal and common spatially uncorrelated red noise (CURN) compared to a model including only a CURN. However, further investigation shows that a disproportionate amount of support for the memory signal comes from three dubious pulsars. Using a more flexible red-noise model in these pulsars reduces the Bayes factor to 1.3. Having found no compelling evidence, we go on to place upper limits on the strain amplitude of GW memory events as a function of sky location and event epoch. These upper limits are computed using a signal model that assumes the existence of a common, spatially uncorrelated red noise in addition to a GW memory signal. The median strain upper limit as a function of sky position is approximately 3.3 × 10−14. We also find that there are some differences in the upper limits as a function of sky position centered around PSR J0613−0200. This suggests that this pulsar has some excess noise that can be confounded with GW memory. Finally, the upper limits as a function of burst epoch continue to improve at later epochs. This improvement is attributable to the continued growth of the pulsar timing array.

A catalogue of Galactic GEMS: Globular cluster Extra-tidal Mock Stars

ABSTRACT This work presents the Globular cluster Extra-tidal Mock Star (GEMS) catalogue of extra-tidal stars and binaries created via three-body dynamical encounters in globular cluster cores. Using the particle-spray code Corespray, we sample $N=50\, 000$ extra-tidal stars and escaped recoil binaries for 159 Galactic globular clusters. Sky positions, kinematics, stellar properties, and escape information are provided for all simulated stars. Stellar orbits are integrated in seven different static and time-varying Milky Way gravitational potential models where the structure of the disc, perturbations from the Large Magellanic Cloud and the mass and sphericity of the Milky Way’s dark matter halo are all investigated. We find that the action coordinates of the mock extra-tidal stars are largely Galactic model independent, where minor offsets and broadening of the distributions between models are likely due to interactions with substructure. Importantly, we also report the first evidence for stellar stream contamination by globular cluster core stars and binaries for clusters with pericentre radii larger than five kiloparsecs. Finally, we provide a quantitative tool that uses action coordinates to match field stars to host clusters with probabilities. Ultimately, combining data from the GEMS catalogue with information of observed stars will allow for association of extra-tidal field stars with any Galactic globular cluster; a requisite tool for understanding population-level dynamics and evolution of clusters in the Milky Way.

Adapting the PyCBC pipeline to find and infer the properties of gravitational waves from massive black hole binaries in LISA

The laser interferometer space antenna (LISA), due for launch in the mid 2030s, is expected to observe gravitational waves (GWs) from merging massive black hole binaries (MBHBs). These signals can last from days to months, depending on the masses of the black holes, and are expected to be observed with high signal to noise ratios (SNRs) out to high redshifts. We have adapted the PyCBC software package to enable a template bank search and inference of GWs from MBHBs. The pipeline is tested on the LISA data challenge’s Challenge 2a (‘Sangria’), which contains MBHBs and thousands of galactic binaries (GBs) in simulated instrumental LISA noise. Our search identifies all six MBHB signals with more than 92% of the optimal SNR. The subsequent parameter inference step recovers the masses and spins within their 90% confidence interval. Sky position parameters have eight high likelihood modes which are recovered but often our posteriors favour the incorrect sky mode. We observe that the addition of GBs biases the parameter recovery of masses and spins away from the injected values, reinforcing the need for a global fit pipeline which will simultaneously fit the parameters of the GB signals before estimating the parameters of MBHBs.

Detecting Technosignatures from Earth-scale Civilizations

Assuming civilizations have similar interplanetary communications and radar capability to NASA’s Deep Space Network, what is the feasibility of intercepting their communications? Interplanetary conjunctions between Earth-like exoplanets, their stars, and other planets in their systems provide one of the most unique and pragmatic opportunities for detecting technosignatures. While eavesdropping on terrestrial communications becomes limited by a planet’s rotation, the beams of satellite communications and interplanetary radar transmissions are tracked, providing the most persistent and powerful opportunity for signal interception. In this study, we present a framework for assessing exoplanet habitability and establishing quantitative bounds for detecting Earth-scale technosignatures from Earth-like planets. These constraints for time, frequency, sky positions, and observatory sensitivity provide recommended observational guidelines for using state-of-the-art and future ground-based radio observatories toward technosignature detection. Applying this framework, 16 exoplanet targets are proposed for radio observation.

Estimating the selection function of Gaia DR3 subsamples

Context. Understanding the intricacies behind the presence and absence of sources in an astronomical catalogue is crucial for the accurate interpretation of astronomical data. In particular, for the multi-dimensional Gaia data, filters and cuts on different parameters or measurements introduce a selection function that may unintentionally alter scientific conclusions in subtle ways. Aims. We aim to develop a methodology to estimate the selection function for different subsamples of stars in the Gaia catalogue. Methods. Comparing the number of stars in a given subsample to that in the overall Gaia catalogue provides an estimate of the subsample membership probability as a function of sky position, magnitude, and colour. The method used to make this estimate must differentiate the stochastic absence of subsample stars from selection effects. When multiplied with the overall Gaia catalogue selection function, this provides the total selection function of the subsample. Results. We present our new method for estimating the selection function by applying it to the sources in Gaia DR3 with heliocentric radial velocity measurements. We also compute the selection function for the stars in the Gaia-Sausage/Enceladus sample, confirming that the apparent asymmetry of its debris across the sky is merely caused by selection effects. Conclusions. The method we have developed estimates the selection function of the stars present in a subsample of Gaia data, given that the subsample is completely contained in the Gaia parent catalogue (for which the selection function is known). This tool is made available in a GaiaUnlimited Python package.

Modeling Apsidal Motion in Eclipsing Binaries Using ELC

Apsidal motion is the precession of the line of apsides in the orbit of a binary star due to perturbations from General Relativity (GR), tides, or third-body interactions. The rate of precession due to tidal effects depends on the interior structures of the stars, and as a result, binaries in which this precession occurs are of great interest. Apsidal motion is observed through the analysis of eclipse times, which reveal small changes in the average interval between successive primary and secondary eclipses, taking all available observed times of eclipse and yielding an estimate of the apsidal rate. Given that this is a single observed quantity, various degeneracies are unavoidably present. Ideally, one would have a model that predicts eclipse times given the orbital and stellar parameters. These parameters for a given binary could then be computed using least squares, provided a suitably large number of eclipse times. Here we use the eclipsing light curve (ELC) program as such a model. The Newtonian equations of motion with additional force terms accounting for GR contributions and tidal distortions are integrated, yielding precise sky positions as a function of time. Times of mid-eclipse and instantaneous orbital elements are computed as a function of time. In this paper, we outline the method and compare numerically computed apsidal rates with standard formulae using a set of 15 binaries based on real systems. For our simulated systems, the derived apsidal rates agree with the standard formula.

Panning for gold, but finding helium: Discovery of the ultra-stripped supernova SN 2019wxt from gravitational-wave follow-up observations

We present the results from multi-wavelength observations of a transient discovered during an intensive follow-up campaign of S191213g, a gravitational wave (GW) event reported by the LIGO-Virgo Collaboration as a possible binary neutron star merger in a low latency search. This search yielded SN 2019wxt, a young transient in a galaxy whose sky position (in the 80% GW contour) and distance (∼150 Mpc) were plausibly compatible with the localisation uncertainty of the GW event. Initially, the transient’s tightly constrained age, its relatively faint peak magnitude (Mi ∼ −16.7 mag), and the r-band decline rate of ∼1 mag per 5 days appeared suggestive of a compact binary merger. However, SN 2019wxt spectroscopically resembled a type Ib supernova, and analysis of the optical-near-infrared evolution rapidly led to the conclusion that while it could not be associated with S191213g, it nevertheless represented an extreme outcome of stellar evolution. By modelling the light curve, we estimated an ejecta mass of only ∼0.1 M⊙, with 56Ni comprising ∼20% of this. We were broadly able to reproduce its spectral evolution with a composition dominated by helium and oxygen, with trace amounts of calcium. We considered various progenitor channels that could give rise to the observed properties of SN 2019wxt and concluded that an ultra-stripped origin in a binary system is the most likely explanation. Disentangling genuine electromagnetic counterparts to GW events from transients such as SN 2019wxt soon after discovery is challenging: in a bid to characterise this level of contamination, we estimated the rate of events with a volumetric rate density comparable to that of SN 2019wxt and found that around one such event per week can occur within the typical GW localisation area of O4 alerts out to a luminosity distance of 500 Mpc, beyond which it would become fainter than the typical depth of current electromagnetic follow-up campaigns.

The probability of identifying the cosmic web environment of galaxies around clusters motivated by the Weave Wide Field Cluster Survey

ABSTRACT Upcoming wide-field spectroscopic surveys will observe galaxies in a range of cosmic web environments in and around galaxy clusters. In this paper, we test and quantify how successfully we will be able to identify the environment of individual galaxies in the vicinity of massive galaxy clusters, reaching out to ∼5R200 into the clusters’ infall region. We focus on the WEAVE Wide Field Cluster Survey (WWFCS), but the methods we develop can be easily generalized to any similar spectroscopic survey. Using numerical simulations of a large sample of massive galaxy clusters from thethreehundred project, we produce mock observations that take into account the selection effects and observational constraints imposed by the WWFCS. We then compare the ‘true’ environment of each galaxy derived from the simulations (cluster core, filament, and neither core nor filament, ‘NCF’) with the one derived from the observational data, where only galaxy sky positions and spectroscopic redshifts will be available. We find that, while cluster core galaxy samples can be built with a high level of completeness and moderate contamination, the filament and NCF galaxy samples will be significantly contaminated and incomplete due to projection effects exacerbated by the galaxies’ peculiar velocities. We conclude that, in the infall regions surrounding massive galaxy clusters, associating galaxies with the correct cosmic web environment is highly uncertain. However, with large enough spectroscopic samples like the ones the WWFCS will provide (thousands of galaxies per cluster, out to 5R200), and the correct statistical treatment that takes into account the probabilities we provide here, we expect we will be able to extract robust and well-quantified conclusions on the relationship between galaxy properties and their environment.