Search

The stars of the Milky Way carry the chemical history of our Galaxy in their atmospheres as they journey through its vast expanse. Like barcodes, we can extract the chemical fingerprints of stars from high-resolution spectroscopy. The fourth data release (DR4) of the Galactic Archaeology with HERMES (GALAH) Survey, based on a decade of observations, provides the chemical abundances of up to 32 elements for 917 588 stars that also have exquisite astrometric data from the $Gaia$ satellite. For the first time, these elements include life-essential nitrogen to complement carbon, and oxygen as well as more measurements of rare-earth elements critical to modern-life electronics, offering unparalleled insights into the chemical composition of the Milky Way. For this release, we use neural networks to simultaneously fit stellar parameters and abundances across the full spectrum, leveraging synthetic grids computed with Spectroscopy Made Easy. These grids account for atomic line formation in non-local thermodynamic equilibrium for 14 elements. In a two-iteration process, we first fit stellar labels for all 1 085 520 spectra, then co-add repeated observations and refine these labels using astrometric data from $Gaia$ and 2MASS photometry, improving the accuracy and precision of stellar parameters and abundances. Our validation thoroughly assesses the reliability of spectroscopic measurements and highlights key caveats for catalogue users. GALAH DR4 represents yet another milestone in Galactic archaeology, combining detailed chemical compositions from multiple nucleosynthetic channels with kinematic information and age estimates. The resulting dataset, covering nearly a million stars, opens new avenues for understanding not only the chemical and dynamical history of the Milky Way, but also the broader questions of the origin of elements and the evolution of planets, stars, and galaxies., Comment: 43 pages, 38 figures to be submitted to PASA. Accompanying the GALAH Data Release 4, see https://www.galah-survey.org and https://cloud.datacentral.org.au/teamdata/GALAH/public/GALAH_DR4/. All code available on http://github.com/svenbuder/GALAH_DR4/ and https://github.com/svenbuder/galah_dr4_paper. Comments welcome

We study open-ended evolution by focusing on computational and information-processing dynamics underlying major evolutionary transitions. In doing so, we consider biological organisms as hierarchical dynamical systems that generate regularities in their phase-spaces through interactions with their environment. These emergent information patterns can then be encoded within the organism's components, leading to self-modelling "tangled hierarchies". Our main conjecture is that when macro-scale patterns are encoded within micro-scale components, it creates fundamental tensions (computational inconsistencies) between what is encodable at a particular evolutionary stage and what is potentially realisable in the environment. A resolution of these tensions triggers an evolutionary transition which expands the problem-space, at the cost of generating new tensions in the expanded space, in a continual process. We argue that biological complexification can be interpreted computation-theoretically, within the G\"odel--Turing--Post recursion-theoretic framework, as open-ended generation of computational novelty. In general, this process can be viewed as a meta-simulation performed by higher-order systems that successively simulate the computation carried out by lower-order systems. This computation-theoretic argument provides a basis for hypothesising the biological arrow of time., Comment: 30 pages, 13 figures

While quiescent galaxies have comparable amounts of cool gas in their outer circumgalactic medium (CGM) compared to star-forming galaxies, they have significantly less interstellar gas. However, open questions remain on the processes causing galaxies to stop forming stars and stay quiescent . Theories suggest dynamical interactions with the hot corona prevent cool gas from reaching the galaxy, therefore predicting the inner regions of quiescent galaxy CGMs are devoid of cool gas. However, there is a lack of understanding of the inner regions of CGMs due to the lack of spatial information in quasar-sightline methods. We present integral-field spectroscopy probing 10--20~kpc (2.4--4.8 R\textsubscript{e}) around a massive quiescent galaxy using a gravitationally lensed star-forming galaxy. We detect absorption from Magnesium (MgII) implying large amounts of cool atomic gas (10\textsuperscript{8.4} -- 10\textsuperscript{9.3} M\textsubscript{$\odot$} with T$\sim$10\textsuperscript{4} Kelvin), in comparable amounts to star-forming galaxies. Lens modeling of Hubble imaging also reveals a diffuse asymmetric component of significant mass consistent with the spatial extent of the MgII absorption, and offset from the galaxy light profile. This study demonstrates the power of galaxy-scale gravitational lenses to not only probe the gas around galaxies, but to also independently probe the mass of the CGM due to it's gravitational effect., Comment: accepted nature communications physics

We present $griz$ photometric light curves for the full 5 years of the Dark Energy Survey Supernova program (DES-SN), obtained with both forced Point Spread Function (PSF) photometry on Difference Images (DIFFIMG) performed during survey operations, and Scene Modelling Photometry (SMP) on search images processed after the survey. This release contains $31,636$ DIFFIMG and $19,706$ high-quality SMP light curves, the latter of which contains $1635$ photometrically-classified supernovae that pass cosmology quality cuts. This sample spans the largest redshift ($z$) range ever covered by a single SN survey ($0.1

Type Ia Supernovae (SNe Ia) are a critical tool in measuring the accelerating expansion of the universe. Recent efforts to improve these standard candles have focused on incorporating the effects of dust on distance measurements with SNe Ia. In this paper, we use the state-of-the-art Dark Energy Survey 5 year sample to evaluate two different families of dust models: empirical extinction models derived from SNe Ia data, and physical attenuation models from the spectra of galaxies. Among the SNe Ia-derived models, we find that a logistic function of the total-to-selective extinction RV best recreates the correlations between supernova distance measurements and host galaxy properties, though an additional 0.02 magnitudes of grey scatter are needed to fully explain the scatter in SNIa brightness in all cases. These empirically-derived extinction distributions are highly incompatible with the physical attenuation models from galactic spectral measurements. From these results, we conclude that SNe Ia must either preferentially select extreme ends of galactic dust distributions, or that the characterisation of dust along the SNe Ia line-of-sight is incompatible with that of galactic dust distributions.

We present a precise measurement of cosmological time dilation using the light curves of 1504 type Ia supernovae from the Dark Energy Survey spanning a redshift range $0.1\lesssim z\lesssim 1.2$. We find that the width of supernova light curves is proportional to $(1+z)$, as expected for time dilation due to the expansion of the Universe. Assuming type Ia supernovae light curves are emitted with a consistent duration $\Delta t_{\rm em}$, and parameterising the observed duration as $\Delta t_{\rm obs}=\Delta t_{\rm em}(1+z)^b$, we fit for the form of time dilation using two methods. Firstly, we find that a power of $b \approx 1$ minimises the flux scatter in stacked subsamples of light curves across different redshifts. Secondly, we fit each target supernova to a stacked light curve (stacking all supernovae with observed bandpasses matching that of the target light curve) and find $b=1.003\pm0.005$ (stat) $\pm\,0.010$ (sys). Thanks to the large number of supernovae and large redshift-range of the sample, this analysis gives the most precise measurement of cosmological time dilation to date, ruling out any non-time-dilating cosmological models at very high significance., Comment: 14 pages, 13 figures. Updated in response to reviewer feedback. Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS)

We measure the current expansion rate of the Universe, Hubble's constant $H_0$, by calibrating the absolute magnitudes of supernovae to distances measured by Baryon Acoustic Oscillations. This `inverse distance ladder' technique provides an alternative to calibrating supernovae using nearby absolute distance measurements, replacing the calibration with a high-redshift anchor. We use the recent release of 1829 supernovae from the Dark Energy Survey spanning $0.01\lt z \lt1.13$ anchored to the recent Baryon Acoustic Oscillation measurements from DESI spanning $0.30 \lt z_{\mathrm{eff}} \lt 2.33$. To trace cosmology to $z=0$, we use the third-, fourth- and fifth-order cosmographic models, which, by design, are agnostic about the energy content and expansion history of the universe. With the inclusion of the higher-redshift DESI-BAO data, the third-order model is a poor fit to both data sets, with the fourth-order model being preferred by the Akaike Information Criterion. Using the fourth-order cosmographic model, we find $H_0=67.19^{+0.66}_{-0.64}\mathrm{~km} \mathrm{~s}^{-1} \mathrm{~Mpc}^{-1}$, in agreement with the value found by Planck without the need to assume Flat-$\Lambda$CDM. However the best-fitting expansion history differs from that of Planck, providing continued motivation to investigate these tensions.

We examine the sky distribution of radio galaxies in the NRAO VLA Sky Survey (NVSS) and the Rapid ASKAP Continuum Survey (RACS). Analyses of these samples have reported tension between their inferred dipoles and the kinematic dipole of the Cosmic Microwave Background (CMB). This represents a challenge to the traditional assumption that the Universe is homogeneous and isotropic on large scales: the cosmological principle. We find that NVSS and RACS contain local radio sources which give a non-negligible contribution to the overall dipole signal. These need to be adequately accounted for since the aim is to probe the composition of the Universe at large scales. By appropriately considering these sources, the inferred dipole amplitude in either sample is reduced. Nonetheless, we find support for a dipole aligning with that of the CMB but larger in amplitude, especially in the joint analysis. However, the 'clustering dipole' - the contribution of local sources to the net inferred dipole - appears to align with the direction of the CMB dipole, and its magnitude increases as deeper nearby sources are considered up to a comoving distance of $\approx 130$ Mpc ($h=0.7$). The significance of this observation in the context of the cosmological principle is unclear and prompts further inquiry., Comment: 15 pages, 13 figures, accepted for publication in MNRAS

Observed evolution of the total mass distribution with redshift is crucial to testing galaxy evolution theories. To measure the total mass distribution, strong gravitational lenses complement the resolved dynamical observations currently limited to $z \lesssim 0.5$. Here we present the lens models for a pilot sample of seven galaxy-scale lenses from the ASTRO3D Galaxy Evolution with Lenses (AGEL) survey. The AGEL lenses, modeled using HST/WFC3-F140W images with Gravitational Lens Efficient Explorer (GLEE) software, have deflector redshifts between $0.3 < z_{\rm defl} < 0.9$. Assuming a power-law density profile with slope $\gamma$, we measure the total density profile for the deflector galaxies via lens modeling. We also measure the stellar velocity dispersions ($\sigma_{\rm obs}$) for four lenses and obtain $\sigma_{\rm obs}$ from SDSS-BOSS for the remaining lenses to test our lens models by comparing observed and model-predicted velocity dispersions. For the seven AGEL lenses, we measure an average density profile slope of $-1.95 \pm 0.09$ and a $\gamma$--$z$ relation that does not evolve with redshift at $z<1$. Although our result is consistent with some observations and simulations, it differs from other studies at $z<1$ that suggest the $\gamma$--$z$ relation evolves with redshift. The apparent conflicts among observations and simulations may be due to a combination of 1) systematics in the lensing and dynamical modeling; 2) challenges in comparing observations with simulations; and 3) assuming a simple power-law for the total mass distribution. By providing more lenses at $z_{\rm defl} > 0.5$, the AGEL survey will provide stronger constraints on whether the mass profiles evolve with redshift as predicted by current theoretical models., Comment: Accepted for publication in the Astrophysical Journal (ApJ) 24 Pages, 8 Figures, 3 Tables

We explore the globular cluster population of NGC 1052-DF4, a dark matter deficient galaxy, using Bayesian inference to search for the presence of rotation. The existence of such a rotating component is relevant to the estimation of the mass of the galaxy, and therefore the question of whether NGC 1052-DF4 is truly deficient of dark matter, similar to NGC 1052-DF2 another galaxy in the same group. The rotational characteristics of seven globular clusters in NGC 1052-DF4 were investigated, finding that a non-rotating kinematic model has a higher Bayesian evidence than a rotating model, by a factor of approximately 2.5. In addition, we find that under the assumption of rotation, its amplitude must be small. This distinct lack of rotation strengthens the case that, based on its intrinsic velocity dispersion, NGC 1052-DF4 is a truly dark matter deficient galaxy., Comment: 9 pages 6 figures. Accepted for publication in PASA

Reverberation mapping is the leading technique used to measure direct black hole masses outside of the local Universe. Additionally, reverberation measurements calibrate secondary mass-scaling relations used to estimate single-epoch virial black hole masses. The Australian Dark Energy Survey (OzDES) conducted one of the first multi-object reverberation mapping surveys, monitoring 735 AGN up to $z\sim4$, over 6 years. The limited temporal coverage of the OzDES data has hindered recovery of individual measurements for some classes of sources, particularly those with shorter reverberation lags or lags that fall within campaign season gaps. To alleviate this limitation, we perform a stacking analysis of the cross-correlation functions of sources with similar intrinsic properties to recover average composite reverberation lags. This analysis leads to the recovery of average lags in each redshift-luminosity bin across our sample. We present the average lags recovered for the H$\beta$, Mg II and C IV samples, as well as multi-line measurements for redshift bins where two lines are accessible. The stacking analysis is consistent with the Radius-Luminosity relations for each line. Our results for the H$\beta$ sample demonstrate that stacking has the potential to improve upon constraints on the $R-L$ relation, which have been derived only from individual source measurements until now., Comment: 20 pages, 15 figures. Accepted by MNRAS

Stellar mergers and accretion events have been crucial in shaping the evolution of the Milky Way (MW). These events have been dynamically identified and chemically characterised using red giants and main-sequence stars. RR Lyrae (RRL) variables can play a crucial role in tracing the early formation of the MW since they are ubiquitous, old (t$\ge$10 Gyr) low-mass stars and accurate distance indicators. We exploited Data Release 3 of the GALAH survey to identify 78 field RRLs suitable for chemical analysis. Using synthetic spectra calculations, we determined atmospheric parameters and abundances of Fe, Mg, Ca, Y, and Ba. Most of our stars exhibit halo-like chemical compositions, with an iron peak around [Fe/H]$\approx -$1.40, and enhanced Ca and Mg content. Notably, we discovered a metal-rich tail, with [Fe/H] values ranging from $-$1 to approximately solar metallicity. This sub-group includes almost ~1/4 of the sample, it is characterised by thin disc kinematics and displays sub-solar $\alpha$-element abundances, marginally consistent with the majority of the MW stars. Surprisingly, they differ distinctly from typical MW disc stars in terms of the s-process elements Y and Ba. We took advantage of similar data available in the literature and built a total sample of 535 field RRLs for which we estimated kinematical and dynamical properties. We found that metal-rich RRLs (1/3 of the sample) likely represent an old component of the MW thin disc. We also detected RRLs with retrograde orbits and provided preliminary associations with the Gaia-Sausage-Enceladus, Helmi, Sequoia, Sagittarius, and Thamnos stellar streams., Comment: Accepted for publication in MNRAS. 29 pages, 20 figures

In the cold dark matter paradigm, our Galaxy is predicted to contain >10000 dark matter subhaloes in the $10^5-10^8M_\odot$ range which should be completely devoid of stars. Stellar streams are sensitive to the presence of these subhaloes, which can create small-scale features in streams if they pass closely enough. Modelling these encounters can therefore, potentially recover the subhalo's properties. In this work, we demonstrate this for streams generated in numerical simulations, modelled on eccentric orbits in a realistic Milky Way potential, which includes the Large Magellanic Cloud and the subhalo itself. We focus on a mock model of the ATLAS-Aliqa Uma stream and inject a $10^7 M_\odot$ subhalo, creating a similar discontinuous morphology to current observations. We then explore how well subhalo properties are recovered using mock stream observations, consisting of no observational errors, as well as assuming realistic observational setups. These setups include present day style observations, and what will be possible with 4MOST and Gaia DR5 in the future. We show that we can recover all parameters describing the impact even with uncertainties matching existing data, including subhalo positions, velocities, mass and scale radius. Modelling the subhalo on an orbit instead of assuming an impulse approximation, we greatly reduce the degeneracy between subhalo mass and velocity seen in previous works. However, we find a slight bias in the subhalo mass (~0.1 dex). This demonstrates that we should be able to reliably extract the properties of subhaloes with stellar streams in the near future., Comment: 13 pages, 14 figures. Submitted to MMRAS. Comments welcome!

We report the discovery and characterisation of TIC 350842552 ("Zvrk"), an apparently isolated, rapidly-rotating ($P_\text{rot} \sim 99\ \mathrm{d}$) red giant observed by TESS in its Southern Continuous Viewing Zone. The star's fast surface rotation is independently verified by the use of p-mode asteroseismology, strong periodicity in TESS and ASAS-SN photometry, and measurements of spectroscopic rotational broadening. A two-component fit to APOGEE spectra indicates a coverage fraction of its surface features consistent with the amplitude of the photometric rotational signal. Variations in the amplitude of its photometric modulations over time suggest the evolution of its surface morphology, and therefore enhanced magnetic activity. We further develop and deploy new asteroseismic techniques to characterise radial differential rotation, and find weak evidence for rotational shear within Zvrk's convective envelope. This feature, in combination with such a high surface rotation rate, is incompatible with models of angular-momentum transport in single-star evolution. Spectroscopic abundance estimates also indicate a high lithium abundance, among other chemical anomalies. Taken together, all of these suggest a planet-ingestion scenario for the formation of this rotational configuration, various models for which we examine in detail., Comment: 31 pages, 17 figures. Accepted for publication in ApJ

The ability to measure precise and accurate stellar effective temperatures ($T_{\rm{eff}}$) and surface gravities ($\log(g)$) is essential in determining accurate and precise abundances of chemical elements in stars. Measuring $\log(g)$ from isochrones fitted to colour-magnitude diagrams of open clusters is significantly more accurate and precise compared to spectroscopic $\log(g)$. By determining the ranges of ages, metallicity, and extinction of isochrones that fit the colour-magnitude diagram, we constructed a joint probability distribution of $T_{\rm{eff}}$ and $\log(g)$. The joint photometric probability shows the complex correlations between $T_{\rm{eff}}$ and $\log(g)$, which depend on the evolutionary stage of the star. We show that by using this photometric prior while fitting spectra, we can acquire more precise spectroscopic stellar parameters and abundances of chemical elements. This reveals higher-order abundance trends in open clusters like traces of atomic diffusion. We used photometry and astrometry provided by the \textit{Gaia} DR3 catalogue, Padova isochrones, and Galactic Archaeology with HERMES (GALAH) DR4 spectra. We analysed the spectra of 1979 stars in nine open clusters, using MCMC to fit the spectroscopic abundances of 26 elements, $T_{\rm{eff}}$, $\log(g)$, $v_{\rm{mic}}$, and $v_{\rm{broad}}$. We found that using photometric priors improves the accuracy of abundances and $\log(g)$, which enables us to view higher-order trends of abundances caused by atomic diffusion in M67 and Ruprecht 147., Comment: 47 Pages. 14 figures in the main body. Accepted by MNRAS on the 8th of February 2024

Blue Horizontal Branch (BHB) stars, excellent distant tracers for probing the Milky Way's halo density profile, are distinguished in the $(g-r)_0$ vs $(i-z)_0$ color space from another class of stars, blue straggler stars (BSs). We develop a Bayesian mixture model to classify BHB stars using high-precision photometry data from the Dark Energy Survey Data Release 2 (DES DR2). We select $\sim2100$ highly-probable BHBs based on their $griz$ photometry and the associated uncertainties, and use these stars to map the stellar halo over the Galactocentric radial range $20 \lesssim R \lesssim 70$ kpc. After excluding known stellar overdensities, we find that the number density $n_\star$ of BHBs can be represented by a power law density profile $n_\star \propto R^{-\alpha}$ with an index of $\alpha=4.28_{-0.12}^{+0.13}$, consistent with existing literature values. In addition, we examine the impact of systematic errors and the spatial inhomogeneity on the fitted density profile. Our work demonstrates the effectiveness of high-precision $griz$ photometry in selecting BHB stars. The upcoming photometric survey from the Rubin Observatory, expected to reach depths 2-3 magnitudes greater than DES during its 10-year mission, will enable us to investigate the density profile of the Milky Way's halo out to the virial radius, unravelling the complex processes of formation and evolution in our Galaxy., Comment: 22 pages, 18 figures. Submitted AAS Journal. Comments Welcome Code used in this work can be found at: https://github.com/sazabi4/Yu2024_BHB/

We present cosmological constraints from the sample of Type Ia supernovae (SN Ia) discovered during the full five years of the Dark Energy Survey (DES) Supernova Program. In contrast to most previous cosmological samples, in which SN are classified based on their spectra, we classify the DES SNe using a machine learning algorithm applied to their light curves in four photometric bands. Spectroscopic redshifts are acquired from a dedicated follow-up survey of the host galaxies. After accounting for the likelihood of each SN being a SN Ia, we find 1635 DES SNe in the redshift range $0.100.5$ SNe compared to the previous leading compilation of Pantheon+, and results in the tightest cosmological constraints achieved by any SN data set to date. To derive cosmological constraints we combine the DES supernova data with a high-quality external low-redshift sample consisting of 194 SNe Ia spanning $0.025

Milky Way globular clusters (GCs) display chemical enrichment in a phenomenon called multiple stellar populations (MSPs). While the enrichment mechanism is not fully understood, there is a correlation between a cluster's mass and the fraction of enriched stars found therein. However, present-day GC masses are often smaller than their masses at the time of formation due to dynamical mass loss. In this work, we explore the relationship between mass and MSPs using the stellar stream 300S. We present the chemical abundances of eight red giant branch member stars in 300S with high-resolution spectroscopy from Magellan/MIKE. We identify one enriched star characteristic of MSPs and no detectable metallicity dispersion, confirming that the progenitor of 300S was a globular cluster. The fraction of enriched stars (12.5\%) observed in our 300S stars is less than the 50\% of stars found enriched in Milky Way GCs of comparable present-day mass ($\sim10^{4.5}$\msun). We calculate the mass of 300S's progenitor and compare it to the initial masses of intact GCs, finding that 300S aligns well with the trend between the system mass at formation and enrichment. 300S's progenitor may straddle the critical mass threshold for the formation of MSPs and can therefore serve as a benchmark for the stellar enrichment process. Additionally, we identify a CH star, with high abundances of \textit{s}-process elements, probably accreted from a binary companion. The rarity of such binaries in intact GCs may imply stellar streams permit the survival of binaries that would otherwise be disrupted., Comment: 15 pages, 7 tables, 8 figures. Accepted for publication in MNRAS

The AO327 drift survey for radio pulsars and transients used the Arecibo telescope from 2010 until its collapse in 2020. AO327 collected ~3100 hours of data at 327 MHz with a time resolution of 82 us and frequency resolution of 24 kHz. While the main motivation for such surveys is the discovery of new pulsars and new, even unforeseen, types of radio transients, they also serendipitously collect a wealth of data on known pulsars. We present an electronic catalog of data and data products on 206 pulsars whose periodic emission was detected by AO327 and are listed in the ATNF catalog of all published pulsars. The AO327 data products include dedispersed time series at full time resolution, average ("folded") pulse profiles, Gaussian pulse profile templates, and an absolute phase reference that allows phase-aligning the AO327 pulse profiles in a physically meaningful manner with profiles from data taken with other instruments. We also provide machine-readable tables with uncalibrated flux measurements at 327 MHz and pulse widths at 50% and 10% of the pulse peak determined from the fitted Gaussian profile templates. The AO327 catalog data set can be used in applications like population analysis of radio pulsars, pulse profile evolution studies in time and frequency, cone and core emission of the pulsar beam, scintillation, pulse intensity distributions, and others. It also constitutes a ready-made resource for teaching signal processing and pulsar astronomy techniques., Comment: Accepted for publication in the Astrophysical Journal

We present AstroLink, an efficient and versatile clustering algorithm designed to hierarchically classify astrophysically-relevant structures from both synthetic and observational data sets. We build upon CluSTAR-ND, a hierarchical galaxy/(sub)halo finder, so that AstroLink now generates a two-dimensional representation of the implicit clustering structure as well as ensuring that clusters are statistically distinct from the noisy density fluctuations implicit within the $n$-dimensional input data. This redesign replaces the three cluster extraction parameters from CluSTAR-ND with a single parameter, $S$ -- the lower statistical significance threshold of clusters, which can be automatically and reliably estimated via a dynamical model-fitting process. We demonstrate the robustness of this approach compared to AstroLink's predecessors by applying each algorithm to a suite of simulated galaxies defined over various feature spaces. We find that AstroLink delivers a more powerful clustering performance while being $\sim27\%$ faster and using less memory than CluSTAR-ND. With these improvements, AstroLink is ideally suited to extracting a meaningful set of hierarchical and arbitrarily-shaped astrophysical clusters from both synthetic and observational data sets -- lending itself as a great tool for morphological decomposition within the context of hierarchical structure formation., Comment: 12 pages, 10 figures. Accepted to MNRAS

The Green Bank North Celestial Cap survey is one of the largest and most sensitive searches for pulsars and transient radio objects. Observations for the survey have finished; priorities have shifted toward long-term monitoring of its discoveries. In this study, we have developed a pipeline to handle large datasets of archival observations and connect them to recent, high-cadence observations taken using the Canadian Hydrogen Intensity Mapping Experiment (CHIME) telescope. This pipeline handles data for 128 pulsars and has produced measurements of spin, positional, and orbital parameters that connect data over observation gaps as large as 2000 days. We have also measured glitches in the timing residuals for five of the pulsars included and proper motion for 19 sources (13 new). We include updates to orbital parameters for 19 pulsars, including 9 previously unpublished binaries. For two of these binaries, we provide updated measurements of post-Keplerian binary parameters, which result in much more precise estimates of the total masses of both systems. For PSR J0509+3801, the much improved measurement of the Einstein delay yields much improved mass measurements for the pulsar and its companion, 1.399(6)\Msun and 1.412(6)\Msun, respectively. For this system, we have also obtained a measurement of the orbital decay due to the emission of gravitational waves: $\dot{P}_{\rm B} = -1.37(7)\times10^{-12}$, which is in agreement with the rate predicted by general relativity for these masses., Comment: accepted for publication in The Astrophysical Journal

The oldest stars in the Milky Way (born in the first few billion years) are expected to have a high density in the inner few kpc, spatially overlapping with the Galactic bulge. We use spectroscopic data from the Pristine Inner Galaxy Survey (PIGS) to study the dynamical properties of ancient, metal-poor inner Galaxy stars. We compute distances using StarHorse, and orbital properties in a barred Galactic potential. With this paper, we release the spectroscopic AAT/PIGS catalogue (13 235 stars). We find that most PIGS stars have orbits typical for a pressure-supported population. The fraction of stars confined to the inner Galaxy decreases with decreasing metallicity, but many very metal-poor stars (VMP, [Fe/H] < -2.0) stay confined (~ 60% stay within 5 kpc). The azimuthal velocity v$_\phi$ also decreases between [Fe/H] = -1.0 and -2.0, but is constant for VMP stars (at ~ 40 km/s). The carbon-enhanced metal-poor (CEMP) stars in PIGS appear to have similar orbital properties compared to normal VMP stars. Our results suggest a possible transition between two spheroidal components - a more metal-rich, more concentrated, faster rotating component, and a more metal-poor, more extended and slower/non-rotating component. We propose that the former may be connected to pre-disc in-situ stars (or those born in large building blocks), whereas the latter may be dominated by contributions from smaller galaxies. This is an exciting era where large metal-poor samples, such as in this work (as well as upcoming surveys, e.g., 4MOST), shed light on the earliest evolution of our Galaxy., Comment: Accepted to MNRAS, new: Figures 4, 8 and 9 and Table 2 - Figure 9 shows two inner halo components

We present a Bayesian analysis of the Quaia sample of 1.3 million quasars as a test of the cosmological principle. This principle postulates that the universe is homogeneous and isotropic on sufficiently large scales, forming the basis of prevailing cosmological models. However, recent analyses of quasar samples have found a matter dipole inconsistent with the inferred kinematic dipole of the Cosmic Microwave Background (CMB), representing a tension with the expectations of the cosmological principle. Here, we explore various hypotheses for the distribution of quasars in Quaia, finding that the sample is influenced by selection effects with significant contamination near the galactic plane. After excising these regions, we find significant evidence that the Quaia quasar dipole is consistent with the CMB dipole, both in terms of the expected amplitude and direction. This result is in conflict with recent analyses, lending support to the cosmological principle and the interpretation that the observed dipole is due to our local departure from the Hubble flow., Comment: 14 pages, 7 figures, accepted for publication in MNRAS; updated to include published correction (2 pages, 2 figures)

Tidal streams and stellar shells are naturally formed in galaxy interactions and mergers. The Giant Stellar Stream (GSS), the North-East (NE), and Western (W) stellar shelves observed in Andromeda galaxy (M31) are examples of these structures and were formed through the merger of M31 and a satellite galaxy. Recent observational papers have provided strong evidence that the shells and GSS originate from a single progenitor. In this paper, we investigate the formation of these two stellar shelves and the detailed nature of their relationship to the GSS. We present numerical simulations of tidal disruption of a satellite galaxy assuming that it is a progenitor of the GSS and the shell system. We represent the progenitor as a dwarf spheroidal galaxy with the stellar mass of $10^{9} M_{\odot}$ and evolve its merger with M31 for 3 Gyrs to reproduce the chemodynamical properties of the NE and W shelves. We find that an initial metallicity of the progenitor with a negative radial gradient of $\Delta$ FeH = -0.3 $\pm$ 0.2, successfully reproduces observed metallicities of the NE, W shelves, and the GSS, showing that all these structures can originate from the same merger event., Comment: 9 pages, 10 figures, MNRAS, accepted

Whilst X-rays and Sunyaev-Zel'dovich observations allow to study the properties of the intra-cluster medium (ICM) of galaxy clusters, their gravitational potential may be constrained using strong gravitational lensing. Although being physically related, these two components are often described with different physical models. Here, we present a unified technique to derive the ICM properties from strong lensing for clusters in hydrostatic equilibrium. In order to derive this model, we present a new universal and self-similar polytropic temperature profile, which we fit using the X-COP sample of clusters. We subsequently derive an analytical model for the electron density, which we apply to strong lensing clusters MACS J0242.5-2132 and MACS J0949.8+1708. We confront the inferred ICM reconstructions to XMM-Newton and ACT observations. We contrast our analytical electron density reconstructions with the best canonical $\beta$-model. The ICM reconstructions obtained prove to be compatible with observations. However they appear to be very sensitive to various dark matter halo parameters constrained through strong lensing (such as the core radius), and to the halo scale radius (fixed in the lensing optimisations). With respect to the important baryonic effects, we make the sensitivity on the scale radius of the reconstruction an asset, and use the inferred potential to constrain the dark matter density profile using ICM observations. The technique here developed should allow to take a new, and more holistic path to constrain the content of galaxy clusters., Comment: 27 pages, 21 figures, published in MNRAS

Transient radio sources, such as fast radio bursts, intermittent pulsars, and rotating radio transients, can offer a wealth of information regarding extreme emission physics as well as the intervening interstellar and/or intergalactic medium. Vital steps towards understanding these objects include characterizing their source populations and estimating their event rates across observing frequencies. However, previous efforts have been undertaken mostly by individual survey teams at disparate observing frequencies and telescopes, and with non-uniform algorithms for searching and characterization. The Petabyte Project (TPP) aims to address these issues by uniformly reprocessing data from several petabytes of radio transient surveys covering two decades of observing frequency (300 MHz-20 GHz). The TPP will provide robust event rate analyses, in-depth assessment of survey and pipeline completeness, as well as revealing discoveries from archival and ongoing radio surveys. We present an overview of TPP's processing pipeline, scope, and our potential to make new discoveries., Comment: 5 pages, 1 figure. Accepted for publication in Proceedings for IAU Symposium 369: The Dawn of Cosmology & Multi-Messenger Studies with Fast Radio Bursts

The exploration of the redshift drift, a direct measurement of cosmological expansion, is expected to take several decades of observation with stable, sensitive instruments. We introduced a new method to probe cosmology which bypasses the long-period observation by observing the redshift difference, an accumulation of the redshift drift, in multiple-image gravitational lens systems. With this, the photons observed in each image will have traversed through different paths between the source and the observer, and so the lensed images will show different redshifts when observed at the same instance. Here, we consider the impact of the underlying cosmology on the observed redshift difference in gravitational lens systems, generating synthetic data for realistic lens models and exploring the accuracy of determined cosmological parameters. We show that, whilst the redshift difference is sensitive to the densities of matter and dark energy within a universe, it is independent of the Hubble constant. Finally, we determine the observational considerations for using the redshift difference as a cosmological probe, finding that one thousand lensed sources are enough to make robust determinations of the underlying cosmological parameters. Upcoming cluster lens surveys, such as the Euclid, are expected to detect a sufficient number of such systems., Comment: 10 pages, 12 figures, 1 table

As our nearest large neighbor, the Andromeda Galaxy provides a unique laboratory for investigating galaxy formation and the distribution and substructure properties of dark matter in a Milky Way-like galaxy. Here, we propose an initial 2-epoch ($\Delta t\approx 5$yr), 2-band Roman survey of the entire halo of Andromeda, covering 500 square degrees, which will detect nearly every red giant star in the halo (10$\sigma$ detection in F146, F062 of 26.5, 26.1AB mag respectively) and yield proper motions to $\sim$25 microarcsec/year (i.e., $\sim$90 km/s) for all stars brighter than F146 $\approx 23.6$ AB mag (i.e., reaching the red clump stars in the Andromeda halo). This survey will yield (through averaging) high-fidelity proper motions for all satellites and compact substructures in the Andromeda halo and will enable statistical searches for clusters in chemo-dynamical space. Adding a third epoch during the extended mission will improve these proper motions by $\sim t^{-1.5}$, to $\approx 11$ km/s, but this requires obtaining the first epoch in Year 1 of Roman operations. In combination with ongoing and imminent spectroscopic campaigns with ground-based telescopes, this Roman survey has the potential to yield full 3-d space motions of $>$100,000 stars in the Andromeda halo, including (by combining individual measurements) robust space motions of its entire globular cluster and most of its dwarf galaxy satellite populations. It will also identify high-velocity stars in Andromeda, providing unique information on the processes that create this population. These data offer a unique opportunity to study the immigration history, halo formation, and underlying dark matter scaffolding of a galaxy other than our own., Comment: Submitted in response to the call for Roman Space Telescope Core Community Survey white papers

We present the discovery and timing solutions of four millisecond pulsars (MSPs) discovered in the Arecibo 327 MHz Drift-Scan Pulsar Survey. Three of these pulsars are in binary systems, consisting of a redback (PSR J2055+1545), a black widow (PSR J1630+3550), and a neutron star-white dwarf binary (PSR J2116+1345). The fourth MSP, PSR J2212+2450, is isolated. We present the multiyear timing solutions as well as polarization properties across a range of radio frequencies for each pulsar. We perform a multiwavelength search for emission from these systems and find an optical counterpart for PSR J2055+1545 in Gaia DR3, as well as a gamma-ray counterpart for PSR J2116+1345 with the Fermi-LAT telescope. Despite the close colocation of PSR J2055+1545 with a Fermi source, we are unable to detect gamma-ray pulsations, likely due to the large orbital variability of the system. This work presents the first two binaries found by this survey with orbital periods shorter than a day; we expect to find more in the 40% of the survey data that have yet to be searched., Comment: 22 pages, 11 figures, 4 tables. Published in ApJ

A fundamental prediction of relativistic cosmologies is that, due to the expansion of space, observations of the distant cosmos should be time dilated and appear to run slower than events in the local universe. Whilst observations of cosmological supernovae unambiguously display the expected redshift-dependent time dilation, this has not been the case for other distant sources. Here we present the identification of cosmic time dilation in a sample of 190 quasars monitored for over two decades in multiple wavebands by assessing various hypotheses through Bayesian analysis. This detection counters previous claims that observed quasar variability lacked the expected redshift-dependent time dilation. Hence, as well as demonstrating the claim that the lack of the redshift dependence of quasar variability represents a significant challenge to the standard cosmological model, this analysis further indicates that the properties of quasars are consistent with them being truly cosmologically distant sources., Comment: 18 pages, 3 figures. Author version of a manuscript under review at Nature Astronomy

We present timing solutions for 21 pulsars discovered in 350 MHz surveys using the Green Bank Telescope (GBT). All were discovered in the Green Bank North Celestial Cap pulsar survey, with the exception of PSR J0957-0619, which was found in the GBT 350 MHz Drift-scan pulsar survey. The majority of our timing observations were made with the GBT at 820 MHz. With a spin period of 37 ms and a 528-day orbit, PSR J0032+6946 joins a small group of five other mildly recycled wide binary pulsars, for which the duration of recycling through accretion is limited by the length of the companion's giant phase. PSRs J0141+6303 and J1327+3423 are new disrupted recycled pulsars. We incorporate Arecibo observations from the NANOGrav pulsar timing array into our analysis of the latter. We also observed PSR J1327+3423 with the Long Wavelength Array, and our data suggest a frequency-dependent dispersion measure. PSR J0957-0619 was discovered as a rotating radio transient, but is a nulling pulsar at 820 MHz. PSR J1239+3239 is a new millisecond pulsar (MSP) in a 4-day orbit with a low-mass companion. Four of our pulsars already have published timing solutions, which we update in this work: the recycled wide binary PSR J0214+5222, the non-eclipsing black widow PSR J0636+5128, the disrupted recycled pulsar J1434+7257, and the eclipsing binary MSP J1816+4510, which is in an 8.7 hr orbit with a redback-mass companion., Comment: 32 pages, 17 figures, 9 tables. Submitted to ApJ

We present the lifetime star formation histories (SFHs) for six ultra-faint dwarf (UFD; $M_V>-7.0$, $ 4.9<\log_{10}({M_*(z=0)}/{M_{\odot}})<5.5$) satellite galaxies of M31 based on deep color-magnitude diagrams constructed from \textit{Hubble Space Telescope} imaging. These are the first SFHs obtained from the oldest main sequence turn-off of UFDs outside the halo of the Milky Way (MW). We find that five UFDs formed at least 50\% of their stellar mass by $z=5$ (12.6~Gyr ago), similar to known UFDs around the MW, but that 10-40\% of their stellar mass formed at later times. We uncover one remarkable UFD, \A{XIII}, which formed only 10\% of its stellar mass by $z=5$, and 75\% in a rapid burst at $z\sim2-3$, a result that is robust to choices of underlying stellar model and is consistent with its predominantly red horizontal branch. This "young" UFD is the first of its kind and indicates that not all UFDs are necessarily quenched by reionization, which is consistent with predictions from several cosmological simulations of faint dwarf galaxies. SFHs of the combined MW and M31 samples suggest reionization did not homogeneously quench UFDs. We find that the least massive MW UFDs ($M_*(z=5) \lesssim 5\times10^4 M_{\odot}$) are likely quenched by reionization, whereas more massive M31 UFDs ($M_*(z=5) \gtrsim 10^5 M_{\odot}$) may only have their star formation suppressed by reionization and quench at a later time. We discuss these findings in the context of the evolution and quenching of UFDs., Comment: 18 pages, 14 figures, 5 appendices, accepted for publication in ApJ

We present a novel test of the cosmological principle: the idea that, on sufficiently large scales, the universe should appear homogeneous and isotropic to observers comoving with the Hubble flow. This is a fundamental assumption in modern cosmology, underpinning the use of the Friedmann-Lema\^itre-Robertson-Walker metric as part of the concordance $\Lambda$CDM paradigm. However, the observed dipole imprinted on the Cosmic Microwave Background (CMB) is interpreted as our departure from the Hubble flow, and such a proper motion will induce a directionally-dependent time dilation over the sky. We illustrate the feasibility of detection of this 'time dilation dipole' and sketch the practical steps involved in its extraction from a catalogue of sources with intrinsic time-scales. In essence, whilst the scale of this dilation is small, being of order of 0.1%, it will in principle be detectable in large scale surveys of variable cosmological sources, such as quasars and supernovae. The degree of alignment of the time dilation dipole with the kinematic dipole derived from the CMB will provide a new assessment of the cosmological principle, and address the tension in dipole measures from other observations., Comment: 9 pages, 4 figures, accepted for publication in MNRAS

We use more than 872,000 mid-to-high resolution (R $\sim$ 20,000) spectra of stars from the GALAH survey to discern the spectra of diffuse interstellar bands (DIBs). We use four windows with the wavelength range from 4718 to 4903, 5649 to 5873, 6481 to 6739, and 7590 to 7890 \AA, giving a total coverage of 967 \AA. We produce $\sim$400,000 spectra of interstellar medium (ISM) absorption features and correct them for radial velocities of the DIB clouds. Ultimately, we combine the 33,115 best ISM spectra into six reddening bins with a range of $0.1 \,\mathrm{mag} < E\mathrm{(B-V)} < 0.7\, \mathrm{mag}$. A total of 183 absorption features in these spectra qualify as DIBs, their fitted model parameters are summarized in a detailed catalogue. From these, 64 are not reported in the literature, among these 17 are certain, 14 are probable and 33 are possible. We find that the broad DIBs can be fitted with a multitude of narrower DIBs. Finally, we create a synthetic DIB spectrum at unit reddening which should allow us to narrow down the possible carriers of DIBs and explore the composition of the ISM and ultimately better model dust and star formation as well as to correct Galactic and extragalactic observations. The majority of certain DIBs show a significant excess of equivalent width when compared to reddening. We explain this with observed lines of sight penetrating more uniform DIB clouds compared to clumpy dust clouds., Comment: 28 pages, 15 figures, 11 tables, accepted for publication in MNRAS

If the origin of life is rare and sensitive to the local conditions at the site of its emergence, then, using the principle of mediocrity within a multiverse framework, we may expect to find ourselves in a universe that is better than usual at creating these necessary conditions. We use this reasoning to investigate several origin of life scenarios to determine whether they are compatible with the multiverse, including the prebiotic soup scenario, hydrothermal vents, delivery of prebiotic material from impacts, and panspermia. We find that most of these scenarios induce a preference toward weaker-gravity universes, and that panspermia and scenarios involving solar radiation or large impacts as a disequilibrium source are disfavored. Additionally, we show that several hypothesized habitability criteria which are disfavored when the origin of life is not taken into account become compatible with the multiverse, and that the emergence of life and emergence of intelligence cannot both be sensitive to disequilibrium production conditions., Comment: 27 pages, 4 figures

Stellar activity and planetary atmospheric properties have the potential to strongly influence habitability. To date, neither have been adequately studied in the multiverse context, so there has been no assessment of how these effects impact the probabilities of observing our fundamental constants. Here, we consider the effects of solar wind, mass loss, and extreme ultra-violet (XUV) flux on planetary atmospheres, how these effects scale with fundamental constants, and how this affects the likelihood of our observations. We determine the minimum atmospheric mass that can withstand erosion, maintain liquid surface water, and buffer diurnal temperature changes. We consider two plausible sources of Earth's atmosphere, as well as the notion that only initially slowly rotating stars are habitable, and find that all are equally compatible with the multiverse. We consider whether planetary magnetic fields are necessary for habitability, and find five boundaries in parameter space where magnetic fields are precluded. We find that if an Earth-like carbon-to-oxygen ratio is required for life, atmospheric effects do not have much of an impact on multiverse calculations. If significantly different carbon-to-oxygen ratios are compatible with life, magnetic fields must not be essential for life, and planet atmosphere must not scale with stellar nitrogen abundance, or else the multiverse would be ruled out to a high degree of confidence., Comment: 23 pages, 2 figures

We build a statistical framework to infer the global properties of the satellite system of the Andromeda galaxy (M31) from the properties of individual dwarf galaxies located in the Pan-Andromeda Archaelogical Survey (PAndAS) and the previously determined completeness of the survey. Using forward modeling, we infer the slope of the luminosity function of the satellite system, the slope of its spatial density distribution, and the size-luminosity relation followed by the dwarf galaxies. We find that the slope of the luminosity function is $\beta=-1.5\pm0.1$. Combined with the spatial density profile, it implies that, when accounting for survey incompleteness, M31 hosts $92_{-26}^{+19}$ dwarf galaxies with $M_\textrm{V}<-5.5$ and a sky-projected distance from M31 between 30 and 300kpc. We conclude that many faint or distant dwarf galaxies remain to be discovered around Andromeda, especially outside the PAndAS footprint. Finally, we use our model to test if the higher number of satellites situated in the hemisphere facing the Milky Way could be explained simply by the detection limits of dwarf galaxy searches. We rule this out at $>99.9\%$ confidence and conclude that this anisotropy is an intrinsic feature of the M31 satellite system. The statistical framework we present here is a powerful tool to robustly constrain the properties of a satellite system and compare those across hosts, especially considering the upcoming start of the Euclid or Rubin large photometric surveys that are expected to uncover a large number of dwarf galaxies in the Local Volume., Comment: Submitted to ApJ - 12 pages, 6 figures, 2 tables

Recent detections of potentially habitable exoplanets around sunlike stars demand increased exploration of the physical conditions that can sustain life, by whatever methods available. Insight into these conditions can be gained by considering the multiverse hypothesis; in a multiverse setting, the probability of living in our universe depends on assumptions made about the factors affecting habitability. Various proposed habitability criteria can be systematically considered to rate each on the basis of their compatibility with the multiverse, generating predictions which can both guide expectations for life's occurrence and test the multiverse hypothesis. Here, we evaluate several aspects of planetary habitability, and show that the multiverse does indeed induce strong preferences among them. We find that the notion that a large moon is necessary for habitability is untenable in the multiverse scenario, as in the majority of parameter space, moons are not necessary to maintain stable obliquity. Further, we consider various proposed mechanisms for water delivery to the early Earth, including delivery from asteroids, both during giant planet formation and a grand tack, delivery from comets, and oxidation of a primary atmosphere by a magma ocean. We find that, depending on assumptions for how habitability depends on water content, some of these proposed mechanisms are disfavored in the multiverse scenario by Bayes factors of up to several hundred., Comment: 26 pages, 3 figures

We investigate the dependence of elemental abundances on physical constants, and the implications this has for the distribution of complex life for various proposed habitability criteria. We consider three main sources of abundance variation: differing supernova rates, alpha burning in massive stars, and isotopic stability, and how each affects the metal-to-rock ratio and the abundances of carbon, oxygen, nitrogen, phosphorus, sulfur, silicon, magnesium, and iron. Our analysis leads to several predictions for which habitability criteria are correct by determining which ones make our observations of the physical constants, as well as a few other observed features of our universe, most likely. Our results indicate that carbon-rich or carbon-poor planets are uninhabitable, slightly magnesium-rich planets are habitable, and life does not depend on nitrogen abundance too sensitively. We also find suggestive but inconclusive evidence that metal-rich planets and phosphorus-poor planets are habitable. These predictions can then be checked by probing regions of our universe that closely resemble normal environments in other universes. If any of these predictions are found to be wrong, the multiverse scenario would predict that the majority of observers are born in universes differing substantially from ours, and so can be ruled out, to varying degrees of statistical significance., Comment: 27 pages, 5 figures

We investigate the sensitivity of a universe's nuclear entropy after Big Bang nucleosynthesis (BBN) to variations in both the baryon-to-photon ratio and the temporal evolution of cosmological expansion. Specifically, we construct counterfactual cosmologies to quantify the degree by which these two parameters must vary from those in our Universe before we observe a substantial change in the degree of fusion, and thus nuclear entropy, during BBN. We find that, while the post-BBN nuclear entropy is indeed linked to baryogenesis and the Universe's expansion history, the requirement of leftover light elements does not place strong constraints on the properties of these two cosmological processes., Comment: 15 pages, 5 figures. It has been accepted for publishing but has not yet appeared. The journal is General Relativity and Gravitation

The Cosmological Principle, that the Universe is homogeneous and isotropic on sufficiently large scales, underpins the standard model of cosmology. However, a recent analysis of 1.36 million infrared-selected quasars has identified a significant tension in the amplitude of the number-count dipole compared to that derived from the CMB, thus challenging the Cosmological Principle. Here we present a Bayesian analysis of the same quasar sample, testing various hypotheses using the Bayesian evidence. We find unambiguous evidence for the presence of a dipole in the distribution of quasars with a direction that is consistent with the dipole identified in the CMB. However, the amplitude of the dipole is found to be 2.7 times larger than that expected from the conventional kinematic explanation of the CMB dipole, with a statistical significance of $5.7\sigma$. To compare these results with theoretical expectations, we sharpen the $\Lambda$CDM predictions for the probability distribution of the amplitude, taking into account a number of observational and theoretical systematics. In particular, we show that the presence of the Galactic plane mask causes a considerable loss of dipole signal due to a leakage of power into higher multipoles, exacerbating the discrepancy in the amplitude. By contrast, we show using probabilistic arguments that the source evolution of quasars improves the discrepancy, but only mildly so. These results support the original findings of an anomalously large quasar dipole, independent of the statistical methodology used., Comment: 15 pages, 7 figures; accepted for publication in MNRAS

We present timing solutions for 12 pulsars discovered in the Green Bank North Celestial Cap (GBNCC) 350 MHz pulsar survey, including six millisecond pulsars (MSPs), a double neutron star (DNS) system, and a pulsar orbiting a massive white dwarf companion. Timing solutions presented here include 350 and 820 MHz Green Bank Telescope data from initial confirmation and follow-up as well as a dedicated timing campaign spanning one year. PSR J1122$-$3546 is an isolated MSP, PSRs J1221$-$0633 and J1317$-$0157 are MSPs in black widow systems and regularly exhibit eclipses, and PSRs J2022+2534 and J2039$-$3616 are MSPs that can be timed with high precision and have been included in pulsar timing array experiments seeking to detect low-frequency gravitational waves. PSRs J1221$-$0633 and J2039$-$3616 have Fermi Large Area Telescope $\gamma$-ray counterparts and also exhibit significant $\gamma$-ray pulsations. We measure proper motion for three of the MSPs in this sample and estimate their space velocities, which are typical compared to those of other MSPs. We have detected the advance of periastron for PSR J1018$-$1523 and therefore measure the total mass of the double neutron star system, $m_{\rm tot}=2.3\pm0.3$ M$_{\odot}$. Long-term pulsar timing with data spanning more than one year is critical for classifying recycled pulsars, carrying out detailed astrometry studies, and shedding light on the wealth of information in these systems post-discovery., Comment: 21 pages, 5 figures, 7 tables

We present the results of a spectroscopic survey of the outskirts of 4 globular clusters -- NGC 1261, NGC 4590, NGC 1904, and NGC 1851 -- covering targets within 1 degree from the cluster centres, with 2dF/AAOmega on the Anglo-Australian Telescope (AAT) and FLAMES on the Very Large Telescope (VLT). We extracted chemo-dynamical information for individual stars, from which we estimated the velocity dispersion profile and the rotation of each cluster. The observations are compared to direct $N$-body simulations and appropriate {\sc limepy}/{\sc spes} models for each cluster to interpret the results. In NGC 1851, the detected internal rotation agrees with existing literature, and NGC 1261 shows some rotation signal beyond the truncation radius, likely coming from the escaped stars. We find that the dispersion profiles for both the observations and the simulations for NGC 1261, NGC 1851, and NGC 1904 do not decrease as the {\sc limepy}/{\sc spes} models predict beyond the truncation radius, where the $N$-body simulations show that escaped stars dominate; the dispersion profile of NGC 4590 follows the predictions of the {\sc limepy}/{\sc spes} models, though the data do not effectively extend beyond the truncation radius. The increasing/flat dispersion profiles in the outskirts of NGC 1261, NGC 1851 and NGC 1904, are reproduced by the simulations. Hence, the increasing/flat dispersion profiles of the clusters in question can be explained by the tidal interaction with the Galaxy without introducing dark matter., Comment: 17 pages, 10 figures. Accepted by MNRAS. Comments are welcome

Based upon a metallicity selection, we identify a significant sub-population of the inner halo globular clusters in the Andromeda Galaxy which we name the Dulais Structure. It is distinguished as a co-rotating group of 10-20 globular clusters which appear to be kinematically distinct from, and on average more metal-poor than, the majority of the inner halo population. Intriguingly, the orbital axis of this Dulais Structure is closely aligned with that of the younger accretion event recently identified using a sub-population of globular clusters in the outer halo of Andromeda, and this is strongly suggestive of a causal relationship between the two. If this connection is confirmed, a natural explanation for the kinematics of the globular clusters in the Dulais Structure is that they trace the accretion of a substantial progenitor (~10^11 Msun) into the halo of Andromeda during the last few billion years, that may have occurred as part of a larger group infall., Comment: 10 pages, 7 figures - accepted for publication in MNRAS

We present a 6-D map of the Orphan-Chenab (OC) stream by combining the data from Southern Stellar Stream Spectroscopic Survey ($S^5$) and {\it Gaia}. We reconstruct the proper motion, radial velocity, distance, on-sky track and stellar density along the stream with spline models. The stream has a total luminosity of $M_V=-8.2$ and metallicity of $\mathrm{[Fe/H]}=-1.9$, similar to classical Milky Way (MW) satellites like Draco. The stream shows drastic changes in its physical width varying from 200 pc to 1 kpc, but a constant line of sight velocity dispersion of 5 km/ss. Despite the large apparent variation in the stellar number density along the stream, the flow rate of stars along the stream is remarkably constant. We model the 6-D stream track by a Lagrange-point stripping method with a flexible MW potential in the presence of a moving extended Large Magellanic Cloud (LMC). This allows us to constrain the mass profile of the MW within the distance range 15.6 < r < 55.5 kpc, with the best measured enclosed mass of $(2.85\pm 0.1)\times 10^{11}\,M_\odot$ within 32.4 kpc. Our stream measurements are highly sensitive to the LMC mass profile with the most precise measurement of its enclosed mass made at 32.8 kpc, $(7.02\pm 0.9)\times10^{10}\, {\rm M}_\odot$. We also detect that the LMC dark matter halo extends to at least 53 kpc. The fitting of the OC stream allows us to constrain the past LMC trajectory and the degree of dynamical friction it experienced. We demonstrate that the stars in the OC stream show large energy and angular momentum spreads caused by LMC perturbation., Comment: accepted to MNRAS; data released with the paper is available on Zenodo https://doi.org/10.5281/zenodo.7222654

Reverberation mapping measurements have been used to constrain the relationship between the size of the broad-line region and luminosity of active galactic nuclei (AGN). This $R-L$ relation is used to estimate single-epoch virial black hole masses, and has been proposed for use to standardise AGN to determine cosmological distances. We present reverberation measurements made with H$\beta$ from the six-year Australian Dark Energy Survey (OzDES) Reverberation Mapping Program. We successfully recover reverberation lags for eight AGN at $0.12

Due to the expansion of our Universe, the redshift of distant objects changes with time. Although the amplitude of this redshift drift is small, it will be measurable with a decade-long campaigns on the next generation of telescopes. Here we present an alternative view of the redshift drift which captures the expansion of the universe in single epoch observations of the multiple images of gravitationally lensed sources. Considering a sufficiently massive lens, with an associated time delay of order decades, simultaneous photons arriving at a detector would have been emitted decades earlier in one image compared to another, leading to an instantaneous redshift difference between the images. We also investigate the effect of peculiar velocities on the redshift difference in the observed images. Whilst still requiring the observational power of the next generation of telescopes and instruments, the advantage of such a single epoch detection over other redshift drift measurements is that it will be less susceptible to systematic effects that result from requiring instrument stability over decade-long campaigns., Comment: 6 pages, 5 figures

The correlation between the broad line region radius and continuum luminosity ($R-L$ relation) of active galactic nuclei (AGN) is critical for single-epoch mass estimates of supermassive black holes (SMBHs). At $z \sim 1-2$, where AGN activity peaks, the $R-L$ relation is constrained by the reverberation mapping (RM) lags of the Mg II line. We present 25 Mg II lags from the Australian Dark Energy Survey (OzDES) RM project based on six years of monitoring. We define quantitative criteria to select good lag measurements and verify their reliability with simulations based on both the damped random walk stochastic model and the re-scaled, re-sampled versions of the observed lightcurves of local, well-measured AGN. Our sample significantly increases the number of Mg II lags and extends the $R-L$ relation to higher redshifts and luminosities. The relative iron line strength $\mathcal{R}_{\rm Fe}$ has little impact on the $R-L$ relation. The best-fit Mg II $R-L$ relation has a slope $\alpha = 0.39 \pm 0.08$ with an intrinsic scatter $\sigma_{\rm rl} = 0.15^{+0.03}_{-0.02}$. The slope is consistent with previous measurements and shallower than the H$\beta$ $R-L$ relation. The intrinsic scatter of the new $R-L$ relation is substantially smaller than previous studies and comparable to the intrinsic scatter of the H$\beta$ $R-L$ relation. Our new $R-L$ relation will enable more precise single-epoch mass estimates and SMBH demographic studies at cosmic noon., Comment: 16 pages, 8 figures; MNRAS, Volume 522, pp.4132

We present the first detailed study comparing the populations of stellar streams in cosmological simulations to observed Milky Way dwarf galaxy streams. In particular, we compare streams identified around Milky Way analogs in the FIRE-2 simulations to stellar streams observed by the Southern Stellar Stream Spectroscopic Survey (S5). For an accurate comparison between the stream populations, we produce mock Dark Energy Survey (DES) observations of the FIRE streams and estimate the detectability of their tidal tails and progenitors. The number and stellar mass distributions of detectable stellar streams is consistent between observations and simulations. However, there are discrepancies in the distributions of pericenters and apocenters, with the detectable FIRE streams, on average, forming at larger pericenters (out to > 110 kpc) and surviving only at larger apocenters (> 40 kpc) than those observed in the Milky Way. We find that the population of high-stellar mass dwarf galaxy streams in the Milky Way is incomplete. Interestingly, a large fraction of the FIRE streams would only be detected as satellites in DES-like observations, since their tidal tails are too low-surface brightness to be detectable. We thus predict a population of yet-undetected tidal tails around Milky Way satellites, as well as a population of fully undetected low surface brightness stellar streams, and estimate their detectability with the Rubin Observatory. Finally, we discuss the causes and implications of the discrepancies between the stream populations in FIRE and the Milky Way, and explore future avenues for tests of satellite disruption in cosmological simulations., Comment: 24 pages, 13 figures, 3 tables, submitted to ApJ

Recent analyses have found intriguing correlations between the colour ($c$) of type Ia supernovae (SNe Ia) and the size of their 'mass-step', the relationship between SN Ia host galaxy stellar mass ($M_\mathrm{stellar}$) and SN Ia Hubble residual, and suggest that the cause of this relationship is dust. Using 675 photometrically-classified SNe Ia from the Dark Energy Survey 5-year sample, we study the differences in Hubble residual for a variety of global host galaxy and local environmental properties for SN Ia subsamples split by their colour. We find a $3\sigma$ difference in the mass-step when comparing blue ($c<0$) and red ($c>0$) SNe. We observe the lowest r.m.s. scatter ($\sim0.14$ mag) in the Hubble residual for blue SNe in low mass/blue environments, suggesting that this is the most homogeneous sample for cosmological analyses. By fitting for $c$-dependent relationships between Hubble residuals and $M_\mathrm{stellar}$, approximating existing dust models, we remove the mass-step from the data and find tentative $\sim 2\sigma$ residual steps in rest-frame galaxy $U-R$ colour. This indicates that dust modelling based on $M_\mathrm{stellar}$ may not fully explain the remaining dispersion in SN Ia luminosity. Instead, accounting for a $c$-dependent relationship between Hubble residuals and global $U-R$, results in $\leq1\sigma$ residual steps in $M_\mathrm{stellar}$ and local $U-R$, suggesting that $U-R$ provides different information about the environment of SNe Ia compared to $M_\mathrm{stellar}$, and motivating the inclusion of galaxy $U-R$ colour in SN Ia distance bias correction., Comment: 19 pages, 8 figures. Published in MNRAS