The 100 most-cited articles in Monthly Notices of the Royal Astronomical Society, 2024

ABSTRACT We present a new determination of the evolving galaxy ultraviolet (UV) luminosity function (LF) over the redshift range $8.5< z< 15.5$ using a combination of several major Cycle-1 JWST imaging programmes – Public Release IMaging for Extragalactic Research, JWST Advanced Deep Extragalactic Survey, and Next Generation Deep Extragalactic Exploratory Public Survey. This multifield approach yields a total of $\simeq 370$ arcmin2 of JWST/NIRCam imaging, reaching (5-$\sigma$) depths of $\simeq 30$ AB mag in the deepest regions. We select a sample of 2548 galaxies with a significant probability of lying at high redshift ($p(z> 8.5)> 0.05$) to undertake a statistical calculation of the UV LF. Our new measurements span $\simeq 4$ mag in UV luminosity at $z=9-12.5$, placing new constraints on both the shape and evolution of the LF at early times. Our measurements yield a new estimate of the early evolution of cosmic star-formation rate density ($\rho _{\rm {SFR}}$) confirming the gradual decline deduced from early JWST studies, at least out to $z \simeq 12$. Finally we show that the…

ABSTRACT Nearly a decade ago, we began to see indications that reionization-era galaxies power hard radiation fields rarely seen at lower redshift. Most striking were detections of nebular C iv emission in what appeared to be typical low-mass galaxies, requiring an ample supply of 48 eV photons to triply ionize carbon. We have obtained deep JWST/NIRSpec R = 1000 spectroscopy of the two z > 6 C iv-emitting galaxies known prior to JWST. Here, we present a rest-UV to optical spectrum of one of these two systems, the multiply-imaged z = 6.1 lensed galaxy RXCJ2248-ID. NIRCam imaging reveals two compact (<22 pc) clumps separated by 220 pc, with one comprising a dense concentration of massive stars (>10 400 M⊙ yr−1 kpc−2) formed in a recent burst. We stack spectra of 3 images of the galaxy (J = 24.8–25.9), yielding a very deep spectrum providing a high-S/N template of strong emission line sources at z > 6. The spectrum reveals narrow high-ionization lines (He ii, C iv, N iv]) with line ratios consistent with powering by massive stars. The rest-optical spectrum is dominated by very strong…

ABSTRACT We report ultra-deep, medium-resolution spectroscopic observations for four quiescent galaxies with log$_{10}(M_*/\mathrm{M_\odot })\gt 11$ at $3 \lt z \lt 5$. These data were obtained with JWST NIRSpec as part of the Early eXtragalactic Continuum and Emission Line Science (EXCELS) survey, which we introduce in this work. The first two galaxies are newly selected from PRIMER UDS imaging, both at $z=4.62$ and separated by 860 pkpc on the sky, within a larger structure for which we confirm several other members. Both formed at $z\simeq 8-10$. These systems could plausibly merge by the present day to produce a local massive elliptical galaxy. The other two ultra-massive quiescent galaxies are previously known at $z=3.99$ and 3.19, with the latter (ZF-UDS-7329) having been the subject of debate as potentially too old and too massive to be accommodated by the $\Lambda$-CDM halo-mass function. Both exhibit high stellar metallicities, and for ZF-UDS-7329 we are able to measure the $\alpha -$enhancement, obtaining [Mg/Fe] = $0.42^{+0.19}_{-0.17}$. We finally evaluate whether these…

ABSTRACT In this paper, we construct the circular velocity curve of the Milky Way out to ∼30 kpc, providing an updated model of the dark matter density profile. We derive precise parallaxes for 120 309 stars with a data-driven model, using APOGEE DR17 spectra combined with GaiaDR3, 2MASS, and WISE photometry. At outer galactic radii up to 30 kpc, we find a significantly faster decline in the circular velocity curve compared to the inner parts. This decline is better fit with a cored Einasto profile with a slope parameter $0.91^{+0.04}_{-0.05}$ than a generalized Navarro–Frenk–White (NFW) profile. The virial mass of the best-fitting dark matter halo profile is only $1.81^{+0.06}_{-0.05}\times 10^{11}$ M⊙, significantly lower than what a generalized NFW profile delivers. We present a study of the potential systematics, affecting mainly large radii. Such a low mass for the Galaxy is driven by the functional forms tested, given that it probes beyond our measurements. It is found to be in tension with mass measurements from globular clusters, dwarf satellites, and streams. Our…

ABSTRACT We identify a low-metallicity ($12+\log ({\rm O}/{\rm H})=7.59$) Ly $\alpha$-emitting galaxy at $z=5.943$ with evidence of a strong Balmer jump, arising from nebular continuum. While Balmer jumps are sometimes observed in low-redshift star-forming galaxies, this galaxy also exhibits a steep turnover in the UV continuum. Such turnovers are typically attributed to absorption by a damped Ly $\alpha$ system (DLA); however, the shape of the turnover and the high observed Ly $\alpha$ escape fraction ($f_{\rm esc,Ly\alpha }~\sim 27~{{\ \rm per\ cent}}$) is also consistent with strong nebular two-photon continuum emission. Modelling the UV turnover with a DLA requires extreme column densities ($N_{\rm HI}\,\,\gt\,\, 10^{23}$ cm$^{-2}$), and simultaneously explaining the high $f_{\rm esc,Ly\alpha }$ requires a fine-tuned geometry. In contrast, modelling the spectrum as primarily nebular provides a good fit to both the continuum and emission lines, motivating scenarios in which (a) we are observing only nebular emission or (b) the ionizing source is powering extreme nebular emission…

ABSTRACT Reionization is thought to be driven by faint star-forming galaxies, but characterizing this population has long remained very challenging. Here, we utilize deep nine-band JADES (JWST Advanced Deep Extragalactic Survey)/NIRCam (Near-Infrared Camera) imaging to study the star-forming and ionizing properties of 756 $z\sim 6-9$ galaxies, including hundreds of very ultraviolet (UV)-faint objects ($M_\mathrm{UV}\gt -18$). The faintest ($m\sim 30$) galaxies in our sample typically have stellar masses of $M_\ast \sim (1-3)\times 10^7\ \mathrm{ M}_\odot$ and young light-weighted ages ($\sim$50 Myr), though some show strong Balmer breaks implying much older ages ($\sim$500 Myr). We find no evidence for extremely massive galaxies ($\gt 3\times 10^{10}\ \mathrm{ M}_\odot$) in our sample. We infer a strong (factor $\gt $2) decline in the typical [O iii]$+$H $\beta$ equivalent widths (EWs) towards very faint $z\sim 6-9$ galaxies, yet a weak UV luminosity dependence on the H $\alpha$ EWs at $z\sim 6$. We demonstrate that these EW trends can be explained if fainter galaxies have…

ABSTRACT Combining the public JWST/NIRCam imaging programs CEERS, PRIMER, and JADES, spanning a total area of $\sim 500\, {\rm arcmin}^2$, we obtain a sample of $\gt $30 000 galaxies at $z_{\rm phot}\sim 4\!-\!9$ that allows us to perform a complete, rest-optical-selected census of the galaxy population at $z\gt 3$. Comparing the stellar mass $M_*$ and the UV-slope $\beta$ distributions between JWST- and HST-selected samples, we generally find very good agreement and no significant biases. Nevertheless, JWST enables us to probe a new population of UV-red galaxies that was missing from previous HST-based Lyman-break galaxy (LBG) samples. We measure galaxy stellar mass functions (SMFs) at $z\sim 4\!-\!9$ down to limiting masses of $10^{7.5}\!-\!10^{8.5}\, {\rm M_\odot }$, finding steep low-mass slopes over the entire redshift range, reaching values of $\alpha \approx -2$ at $z\gtrsim 6$. At the high-mass end, UV-red galaxies dominate at least out to $z\sim 6$. The implied redshift evolution of the SMF suggests a rapid build-up of massive dust-obscured or quiescent galaxies from $z\sim…

ABSTRACT Numerical simulations have become one of the key tools used by theorists in all the fields of astrophysics and cosmology. The development of modern tools that target the largest existing computing systems and exploit state-of-the-art numerical methods and algorithms is thus crucial. In this paper, we introduce the fully open-source highly-parallel, versatile, and modular coupled hydrodynamics, gravity, cosmology, and galaxy-formation code Swift. The software package exploits hybrid shared- and distributed-memory task-based parallelism, asynchronous communications, and domain-decomposition algorithms based on balancing the workload, rather than the data, to efficiently exploit modern high-performance computing cluster architectures. Gravity is solved for using a fast-multipole-method, optionally coupled to a particle mesh solver in Fourier space to handle periodic volumes. For gas evolution, multiple modern flavours of Smoothed Particle Hydrodynamics are implemented. Swift also evolves neutrinos using a state-of-the-art particle-based method. Two complementary networks of…

ABSTRACT Pulsar timing arrays search for nanohertz-frequency gravitational waves by regularly observing ensembles of millisecond pulsars over many years to look for correlated timing residuals. Recently the first evidence for a stochastic gravitational wave background has been presented by the major arrays, with varying levels of significance (${\sim} 2\sigma \!-\! 4\sigma$). In this paper, we present the results of background searches with the MeerKAT Pulsar Timing Array. Although of limited duration (4.5 yr), the ${\sim} 250\,000$ arrival times with a median error of just $3 \, \mu {\rm s}$ on 83 pulsars make it very sensitive to spatial correlations. Detection of a gravitational wave background requires careful modelling of noise processes to ensure that any correlations represent a fit to the underlying background and not other misspecified processes. Under different assumptions about noise processes, we can produce either what appear to be compelling Hellings–Downs correlations of high significance ($3\sigma \!-\! 3.4\sigma$) with a spectrum close to that which is predicted, or…

ABSTRACT The power-law slope of the rest-ultraviolet (UV) continuum (fλ ∝ λβ) is a key metric of early star-forming galaxies, providing one of our only windows into the stellar populations and physical conditions of z ≳ 10 galaxies. Expanding upon previous studies with limited sample sizes, we leverage deep imaging from the JWST Advanced Deep Extragalactic Survey (JADES) to investigate the UV slopes of 179 z ≳ 9 galaxies with apparent magnitudes of mF200W ≃ 26–31, which display a median UV slope of β = −2.4. We compare to a statistical sample of z ≃ 5–9 galaxies, finding a shift towards bluer rest-UV colours at all $M_{\rm UV}$. The most UV-luminous z ≳ 9 galaxies are significantly bluer than their lower redshift counterparts, representing a dearth of moderately red galaxies within the first 500 Myr. At yet earlier times, the z ≳ 11 galaxy population exhibits very blue UV slopes, implying very low impact from dust attenuation. We identify a robust sample of 44 galaxies with β ≲ −2.8, which have spectral energy distributions requiring models of density-bounded H ii regions and median…

ABSTRACT A surprising finding of recent studies is the large number of Active Galactic Nuclei (AGN) associated with moderately massive black holes ($\rm \log (M_\bullet /M_\odot)\sim 6-8$), in the first billion years after the big bang (z > 5). In this context, a relevant finding has been the large fraction of candidate dual AGN, both at large separations (several kpc) and in close pairs (less than a kpc), likely in the process of merging. Frequent black hole merging may be a route for black hole growth in the early universe; however, previous findings are still tentative and indirect. We present JWST/NIRSpec-IFU observations of a galaxy at z = 7.15 in which we find evidence for a $\rm \log (M_\bullet /M_\odot)\sim 7.7$ accreting black hole, as traced by a broad component of H β emission, associated with the Broad Line Region (BLR) around the black hole. This BLR is offset by 620 pc in projection from the centroid of strong rest-frame optical emission, with a velocity offset of ∼40 km s–1. The latter region is also characterized by (narrow) nebular emission features typical of AGN,…

ABSTRACT The energy injection through Hawking evaporation has been used to put strong constraints on primordial black holes as a dark matter candidate at masses below 1017 g. However, Hawking’s semiclassical approximation breaks down at latest after half-decay. Beyond this point, the evaporation could be significantly suppressed, as was shown in recent work. In this study we review existing cosmological and astrophysical bounds on primordial black holes, taking this effect into account. We show that the constraints disappear completely for a reasonable range of parameters, which opens a new window below 1010 g for light primordial black holes as a dark matter candidate.

ABSTRACT JWST has discovered a large population of active galactic nuclei (AGN) at high redshift, which are weak in the X-rays. Here we present the NIRSpec spectrum of the most extreme of these objects, GN-28074, an AGN at $z=2.26$ with prominent hydrogen and He i broad lines, and with the highest limit on the bolometric to X-ray luminosity ratio among all spectroscopically confirmed AGN in GOODS. This source is also characterized by a mid-IR excess, likely associated with the AGN torus’ hot dust. The high bolometric luminosity and moderate redshift of this AGN allow us to explore its properties more in depth relative to other JWST-discovered AGN. The NIRSpec spectrum reveals prominent, slightly blueshifted absorption of H$\alpha$, H$\beta$, and He i $\lambda$10830. The Balmer absorption lines require gas with densities of $n_{\rm H}\gt 10^8~{\rm cm}^{-3}$, consistent with clouds in the broad-line region (BLR). This finding suggests that part of the X-ray weakness is due to high (Compton-thick) X-ray absorption by clouds in the BLR, or in its outer regions. GN-28074 is also…

ABSTRACT We use the Cosmic Archaeology Tool (CAT) semi-analytical model to explore the contribution of Population (Pop) III/II stars and active galactic nuclei (AGNs) to the galaxy ultraviolet (UV) luminosity function (LF) evolution at 4 ≤ z ≤ 20. We compare in particular with recent JWST data in order to explore the apparent tension between observations and theoretical models in the number density of bright galaxies at z ≳ 10. The model predicts a star formation history dominated by UV faint (MUV > −18) galaxies, with a Pop III contribution of $\lesssim 10~{{\ \rm per\ cent}}$ ($\lesssim 0.5~{{\ \rm per\ cent}}$) at z ≃ 20 (z ≃ 10). Stars are the primary sources of cosmic reionization, with $5~{{\ \rm per\ cent}}\!-\!10~{{\ \rm per\ cent}}$ of ionizing photons escaping into the intergalatic medium at 5 ≤ z ≤ 10, while the contribution of unobscured AGNs becomes dominant only at z ≲ 5. The predicted stellar and AGN UV LFs reproduce the observational data at 5 ≲ z ≲ 9–10. At higher redshift, CAT predicts a steeper evolution in the faint-end slope (MUV > –18), and a number density of…

ABSTRACT We present an analysis of the ultraviolet (UV) continuum slopes (β) for a sample of 172 galaxy candidates at 8 < zphot < 16 selected from a combination of JWST NIRCam imaging and COSMOS/UltraVISTA ground-based near-infrared imaging. Focusing primarily on a new sample of 121 galaxies at 〈z〉 ≃ 11 selected from ≃ 320 arcmin2 of public JWST imaging data across 15 independent data sets, we investigate the evolution of β in the galaxy population at z ≥ 9. We find a significant trend between β and redshift, with the inverse-variance weighted mean UV slope evolving from 〈β〉 = −2.17 ± 0.06 at z = 9.5 to 〈β〉 = −2.59 ± 0.06 at z = 11.5. Based on a comparison with stellar population models including nebular continuum emission, we find that at z > 10.5 the average UV continuum slope is consistent with the intrinsic blue limit of dust-free stellar populations (βint ≃ −2.6). These results suggest that the moderately dust-reddened galaxy population at z < 10 was essentially unattenuated at z ≃ 11. The extremely blue galaxies being uncovered at z > 10 place important constraints on dust…

ABSTRACT We use deep JWST/NIRSpec R ∼ 1000 slit spectra of 113 galaxies at $1.7 < z < 3.5$, selected from the mass-complete Blue Jay survey, to investigate the prevalence and typical properties of neutral gas outflows at cosmic noon. We detect excess Na id absorption (beyond the stellar contribution) in 46 per cent of massive galaxies (log M*/M⊙ > 10), with similar incidence rates in star-forming and quenching systems. Half of the absorption profiles are blueshifted by at least 100 km s−1, providing unambiguous evidence for neutral gas outflows. Galaxies with strong Na id absorption are distinguished by enhanced emission line ratios consistent with AGN ionization. We conservatively measure mass outflow rates of 3–100 M⊙ yr−1; comparable to or exceeding ionized gas outflow rates measured for galaxies at similar stellar mass and redshift. The outflows from the quenching systems (log(sSFR)[yr−1] ≲ −10) have mass loading factors of 4–360, and the energy and momentum outflow rates exceed the expected injection rates from supernova explosions, suggesting that these galaxies could possibly…

ABSTRACT We present the first comprehensive halo occupation distribution (HOD) analysis of the Dark Energy Spectroscopic Instrument (DESI) One-Percent Survey luminous red galaxy (LRG) and Quasi Stellar Object (QSO) samples. We constrain the HOD of each sample and test possible HOD extensions by fitting the redshift-space galaxy 2-point correlation functions in 0.15 < r < 32 h−1 Mpc in a set of fiducial redshift bins. We use AbacusSummit cubic boxes at Planck 2018 cosmology as model templates and forward model galaxy clustering with the AbacusHOD package. We achieve good fits with a standard HOD model with velocity bias, and we find no evidence for galaxy assembly bias or satellite profile modulation at the current level of statistical uncertainty. For LRGs in 0.4 < z < 0.6, we infer a satellite fraction of $f_\mathrm{sat} = 11\pm 1~{y{\ \mathrm{per\,cent}}}$, a mean halo mass of $\log _{10}\overline{M}_h/M_\odot =13.40^{+0.02}_{-0.02}$, and a linear bias of $b_\mathrm{lin} = 1.93_{-0.04}^{+0.06}$. For LRGs in 0.6 < z < 0.8, we find $f_\mathrm{sat}=14\pm 1~{{\ \mathrm{per\,cent}}}$,…

ABSTRACT We present the SARAO MeerKAT Galactic Plane Survey (SMGPS), a 1.3 GHz continuum survey of almost half of the Galactic Plane (251° ≤l ≤ 358° and 2° ≤l ≤ 61° at |b| ≤ $1{_{.}^{\circ}}5$). SMGPS is the largest, most sensitive, and highest angular resolution 1 GHz survey of the plane yet carried out, with an angular resolution of 8 arcsec and a broad-band root-mean-square sensitivity of ∼10–20 μJy beam−1. Here, we describe the first publicly available data release from SMGPS which comprises data cubes of frequency-resolved images over 908–1656 MHz, power-law fits to the images, and broad-band zeroth moment integrated intensity images. A thorough assessment of the data quality and guidance for future usage of the data products are given. Finally, we discuss the tremendous potential of SMGPS by showcasing highlights of the Galactic and extragalactic science that it permits. These highlights include the discovery of a new population of non-thermal radio filaments; identification of new candidate supernova remnants, pulsar wind nebulae and planetary nebulae; improved…

ABSTRACT The presence of a stellar bar in a disc galaxy indicates that the galaxy hosts in its main part a dynamically settled disc and that bar-driven processes are taking place in shaping its evolution. Studying the cosmic evolution of the bar fraction in disc galaxies is therefore essential to understand galaxy evolution in general. Previous studies have found, using the Hubble Space Telescope (HST), that the bar fraction significantly declines from the local Universe to redshifts near one. Using the first four pointings from the JWST Cosmic Evolution Early Release Science Survey and the initial public observations for the Public Release Imaging for Extragalactic Research, we extend the studies of the bar fraction in disc galaxies to redshifts 1 ≤ $z$ ≤ 3, that is, for the first time beyond redshift two. We only use galaxies that are also present in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey on the Extended Groth Strip and Ultra Deep Survey HST observations. An optimized sample of 368 close-to-face-on galaxies is visually classified to find the fraction of bars…

ABSTRACT We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical ‘zoom-in’ simulations of 40 Milky Way-mass haloes and 26 dwarf galaxy–mass haloes run with the moving-mesh code arepo. Auriga adopts the Lambda cold dark matter cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully–Fisher relation, the star-forming (SF) main sequence, and H i gas fraction/disc thickness. Finally, we…

ABSTRACT We present high-resolution ($\sim 0.14$ arcsec = 710 pc) Atacama Large Millimetre/submillimetre Array [C ii] 158 $\mu$m and dust continuum follow-up observations of REBELS-25, a [C ii]-luminous ($L_{\mathrm{[CII]}}=(1.7\pm 0.2)\times 10^9\, \mathrm{L_{\odot }}$) galaxy at redshift $z=7.3065\pm 0.0001$. These high-resolution, high signal-to-noise observations allow us to study the sub-kpc morphology and kinematics of this massive ($M_* = 8^{+4}_{-2} \times 10^9 \mathrm{{\rm M}_{\odot }}$) star-forming (SFR$_{\mathrm{UV+IR}} = 199^{+101}_{-63} \mathrm{{\rm M}_{\odot }} \mathrm{yr}^{-1}$) galaxy in the Epoch of Reionization. By modelling the kinematics with $^{\mathrm{3D}}$BAROLO, we find it has a low-velocity dispersion ($\bar{\sigma } = 33^{+9}_{-7}$ km s$^{-1}$) and a high ratio of ordered-to-random motion ($V_{\mathrm{rot, ~max}}/\bar{\sigma } = 11 ^{+6}_{-5}$), indicating that REBELS-25 is a dynamically cold disc. Additionally, we find that the [C ii] distribution is well fit by a near-exponential disc model, with a Sérsic index, n, of $1.3 \pm 0.2$, and we see tentative…

ABSTRACT We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest resolution run of the IllustrisTNG suite of lambda-cold dark matter magnetohydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogues, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogues, that is, galaxies with stellar discy morphology, with a stellar mass in the range of $M_* = 10^{10.5 - 11.2}~\rm {\rm M}_{\odot }$, and within an MW-like 500 kpc-scale environment at $z=0$. These are resolved with baryonic (dark matter) mass resolution of $8.5\times 10^4\rm {\rm M}_{\odot }$ ($4.5\times 10^5\rm {\rm M}_{\odot }$) and $\sim 150$ pc of average gas spatial resolution in the star-forming regions. The majority of TNG50 MW/M31 analogues at $z=0$ exhibit a bar, 60 per cent are star-forming, the sample includes three Local Group-like systems, and a number of galaxies host one or more satellites as massive as, for example, the Magellanic Clouds. Even within such a relatively…

ABSTRACT Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev–Zel’dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, ‘baryonification’. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model, and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters, with the joint WL + kSZ analysis constraining $S_8=0.823^{+0.019}_{-0.020}$. We measure the suppression of the…

ABSTRACT We present AstroCLIP, a single, versatile model that can embed both galaxy images and spectra into a shared, physically meaningful latent space. These embeddings can then be used – without any model fine-tuning – for a variety of downstream tasks including (1) accurate in-modality and cross-modality semantic similarity search, (2) photometric redshift estimation, (3) galaxy property estimation from both images and spectra, and (4) morphology classification. Our approach to implementing AstroCLIP consists of two parts. First, we embed galaxy images and spectra separately by pre-training separate transformer-based image and spectrum encoders in self-supervised settings. We then align the encoders using a contrastive loss. We apply our method to spectra from the Dark Energy Spectroscopic Instrument and images from its corresponding Legacy Imaging Survey. Overall, we find remarkable performance on all downstream tasks, even relative to supervised baselines. For example, for a task like photometric redshift prediction, we find similar performance to a specifically trained…

ABSTRACT The MM ro-vibrational line list for methane (12CH4) is presented; MM covers wavelengths λ > 0.83 μm (wavenumbers up to 12 000 cm−1) and contains over 50 billion transitions between 9155 208 states with total angular momentum J ≤ 60. MM was generated through solution of the nuclear motion Schrödinger equation using variational program trove for an empirically derived potential energy surface (PES) and a new high-level ab initio dipole moment surface. The PES was constructed by fitting the ro-vibrational energies of CH4 to a set of highly accurate, experimentally derived energies. Molecular states are classified using the Td(M) symmetry group and are fully assigned with rotation and vibration quantum numbers. The MM line list is adapted to high-resolution applications by replacing the calculated ro-vibrational energies with the experimentally derived values where available, namely for 23 208 states with J ≤ 27 below 9986 cm−1. Doing so leads to over 1000 000 experimentally derived CH4 lines compared to approximately $330\, 000$ lines of 12CH4 in the HITRAN data base. The MM…

ABSTRACT Extreme emission line galaxies (EELGs) exhibit large equivalent widths (EW) in their rest-optical emission lines ([O iii]$\lambda 5007$ or H $\alpha$ rest-frame EW$\gt 750$ Å) which can be tied to a recent upturn in star formation rate (SFR), due to the sensitivity of the nebular line emission and the rest-optical continuum to young ($\lt 10$ Myr) and evolved stellar populations, respectively. By studying a sample of 85 star-forming galaxies (SFGs), spanning the redshift and magnitude interval $3 \lt z\lt 9.5$ and $-16\gt $M$\rm _{UV}\gt -21$, in the JWST Advanced Deep Extragalactic Survey (JADES) with NIRSpec/prism spectroscopy, we determine that SFGs initiate an EELG phase when entering a significant burst of star formation, with the highest EWs observed in EELGs with the youngest luminosity-weighted ages ($\lt 5$ Myr) and the highest burst intensity (those with the greatest excess between their current and long-term average SFR). We spectroscopically confirm that a greater proportion of SFGs are in an EELG phase at high redshift in our UV-selected sample ($61\pm 4~{{\…

ABSTRACT We use NIRCam imaging from the JWST Advanced Deep Extragalactic Survey (JADES) to study the ionizing properties of a sample of 14 652 galaxies at $3 \le z_{\rm {phot}} \le 9$, 90 per cent complete in stellar mass down to log(M$_{\star }$/[M$_{\odot }])\approx 7.5$. Out of the full sample, 1620 of the galaxies have spectroscopic redshift measurements from the literature. We use the spectral energy distribution fitting code Prospector to fit all available photometry and infer galaxy properties. We find a significantly milder evolution of the ionizing photon production efficiency ($\xi _{\rm {ion}}$) with redshift and UV magnitude than previously reported. Interestingly, we observe two distinct populations in $\xi _{\rm {ion}}$, distinguished by their burstiness (given by SFR$_{10}$/SFR$_{100}$). Both populations show the same evolution with z and M$_{\rm {UV}}$, but have a different $\xi _{\rm {ion}}$ normalization. We convolve the more representative $\log (\xi _{\rm {ion}} (z,\text{M}_{\rm {UV}}))$ relations (accounting for $\sim 97$ % of the sample), with luminosity…

ABSTRACT The origin of cosmic high-energy neutrinos remains largely unexplained. For high-energy neutrino alerts from IceCube, a coincidence with time-variable emission has been seen for three different types of accreting black holes: (1) a gamma-ray flare from a blazar (TXS 0506+056), (2) an optical transient following a stellar tidal disruption event (TDE; AT2019dsg), and (3) an optical outburst from an active galactic nucleus (AGN; AT2019fdr). For the latter two sources, infrared follow-up observations revealed a powerful reverberation signal due to dust heated by the flare. This discovery motivates a systematic study of neutrino emission from all supermassive black hole with similar dust echoes. Because dust reprocessing is agnostic to the origin of the outburst, our work unifies TDEs and high-amplitude flares from AGN into a population that we dub accretion flares. Besides the two known events, we uncover a third flare that is coincident with a PeV-scale neutrino (AT2019aalc). Based solely on the optical and infrared properties, we estimate a significance of 3.6σ for this…

ABSTRACT We report the chemical abundance pattern of GS_3073, a galaxy hosting an overmassive active black hole at $z=5.55$, by leveraging observations from JWST/NIRSpec and Very Large Telescope/VIsible Multi-Object Spectrograph. Based on the rest-frame ultraviolet (UV) emission lines, which trace high-density ($\sim 10^5~{\rm cm^{-3}}$) and highly ionized gas, we derive $\rm \log (N/O) = 0.42^{+0.13}_{-0.10}$. At an estimated metallicity of $0.2~Z_{\odot }$, this is the most extreme nitrogen-rich object found by JWST thus far. In comparison, the relative carbon abundance derived is $\rm \log (C/O) = -0.38^{+0.13}_{-0.11}$, which is not significantly higher than those in local galaxies and stars with similar metallicities. We also found potential detection of [Fe vii]$\lambda 6087$ and [Fe xiv]$\lambda 5303$, both blended with [Ca v]. We inferred a range of Fe abundances compatible with those in local stars and galaxies. Overall, the chemical abundance pattern of GS_3073 is compatible with enrichment by supermassive stars with $M_* \gtrsim 1000~M_\odot$, asymptotic giant branch…

ABSTRACT Fulfilling the rich promise of rapid advances in time-domain astronomy is only possible through confronting our observations with physical models and extracting the parameters that best describe what we see. Here, we introduce redback; a Bayesian inference software package for electromagnetic transients. redback provides an object-orientated python interface to over 12 different samplers and over 100 different models for kilonovae, supernovae, gamma-ray burst afterglows, tidal disruption events, engine-driven transients among other explosive transients. The models range in complexity from simple analytical and semi-analytical models to surrogates built upon numerical simulations accelerated via machine learning. redback also provides a simple interface for downloading and processing data from various catalogues such as Swift and FINK. The software can also serve as an engine to simulate transients for telescopes such as the Zwicky Transient Facility and Vera Rubin with realistic cadences, limiting magnitudes, and sky coverage or a hypothetical user-constructed survey or a…

ABSTRACT In recent years, cosmological hydrodynamical simulations have proven their utility as key interpretative tools in the study of galaxy formation and evolution. In this work, we present a comparative analysis of the baryon cycle in three publicly available, leading cosmological simulation suites: EAGLE, IllustrisTNG, and SIMBA. While these simulations broadly agree in terms of their predictions for the stellar mass content and star formation rates of galaxies at $z\approx 0$, they achieve this result for markedly different reasons. In EAGLE and SIMBA, we demonstrate that at low halo masses ($M_{\rm 200c}\lesssim 10^{11.5}\, \mathrm{M}_{\odot }$), stellar feedback (SF)-driven outflows can reach far beyond the scale of the halo, extending up to $2\!-\!3\times R_{\rm 200c}$. In contrast, in TNG, SF-driven outflows, while stronger at the scale of the interstellar medium, recycle within the circumgalactic medium (within $R_{\rm 200c}$). We find that active galactic nucleus (AGN)-driven outflows in SIMBA are notably potent, reaching several times $R_{\rm 200c}$ even at halo masses…

ABSTRACT The finding of active galactic nuclei (AGNs) in dwarf galaxies has important implications for galaxy evolution and supermassive black hole formation models. Yet, how AGNs in dwarf galaxies form is still debated, in part due to scant demographics. We use the MaNGA survey, comprising ∼10 000 galaxies at z < 0.15, to identify AGN dwarf galaxies using a spaxel-by-spaxel classification in three spatially-resolved emission-line diagnostic diagrams (the [N ii]-, [S ii]-, and [O i]-BPT) and the WHAN diagram. This yields a sample of 664 AGN dwarf galaxies, the largest to date, and an AGN fraction of $\sim 20~{{\ \rm per\ cent}}$ that is significantly larger than that of single-fibre-spectroscopy studies (i.e. $\sim 1~{{\ \rm per\ cent}}$). This can be explained by the lower bolometric luminosity (<1042 erg s−1) and accretion rate (sub-Eddington) of the MaNGA AGN dwarf galaxies. We additionally identify 1176 SF-AGNs (classified as star-forming in the [N ii]-BPT but as AGNs in the [S ii]- and [O i]-BPT), 122 Composite, and 173 LINER sources. The offset between the optical centre of…

ABSTRACT We introduce version 2.0 of the shark semi-analytic model of galaxy formation after many improvements to the physics included. The most significant being (i) a model describing the exchange of angular momentum (AM) between the interstellar medium and stars; (ii) a new active galactic nuclei feedback model which has two modes, a wind and a jet mode, with the jet mode tied to the jet energy production; (iii) a model tracking the development of black hole (BH) spins; (iv) more sophisticated modelling of environmental effects on satellite galaxies; and (v) automatic parameter exploration using Particle Swarm Optimization. We focus on two timely research topics: the structural properties of galaxies and the quenching of massive galaxies. For the former, sharkv2.0 is capable of producing a more realistic stellar size–mass relation with a plateau marking the transition from disc- to bulge-dominated galaxies, and scaling relations between specific AM and mass that agree well with observations. For the quenching of massive galaxies, sharkv2.0 produces massive galaxies that are more…

ABSTRACT The observation of the global 21 cm signal produced by neutral hydrogen gas in the intergalactic medium (IGM) during the Dark Ages, Cosmic Dawn, and Epoch of Reionization requires measurements with extremely well-calibrated wideband radiometers. We describe the design and characterization of the Mapper of the IGM Spin Temperature (MIST), which is a new ground-based, single-antenna, global 21 cm experiment. The design of MIST was guided by the objectives of avoiding systematics from an antenna ground plane and cables around the antenna, as well as maximizing the instrument’s on-sky efficiency and portability for operations at remote sites. We have built two MIST instruments, which observe in the range 25–105 MHz. For the 21 cm signal, this frequency range approximately corresponds to redshifts 55.5 > z > 12.5, encompassing the Dark Ages and Cosmic Dawn. The MIST antenna is a horizontal blade dipole of 2.42 m in length, 60 cm in width, and 52 cm in height above the ground. This antenna operates without a metal ground plane. The instruments run on 12 V batteries and have a…

ABSTRACT We describe JWST/NIRSpec prism measurements of Ly α emission in z ≳ 5 galaxies. We identify Ly α detections in 10 out of 69 galaxies with robust rest-optical emission-line redshift measurements at 5 ≤ z < 7 in the Cosmic Evolution Early Release Science (CEERS) and DDT-2750 observations of the Extended Groth Strip field. Galaxies at z ≃ 6 with faint continuum (F150W=27–29 mag) are found with extremely large rest-frame Ly α equivalent widths (EWs; ranging up to 286 Å). Likely Ly α detections are also seen in two new z > 7 galaxies (z = 7.49 and 7.17) from the second epoch of CEERS observations, both showing large Ly α EWs that likely indicate significant transmission through the intergalactic medium (IGM). We measure high Ly α escape fractions in the 12 Ly α emitters in our sample (median 0.28), two of which show $f_{\rm esc}^{ {\rm Ly}\alpha }$ near unity (>0.80). We find that $50_{-11}^{+11}$ per cent of z ≃ 6 galaxies with [O iii] + H β EW>1000 Å have $f_{\rm esc}^{ {\rm Ly}\alpha }$ >0.2, consistent with the fractions found in lower redshift samples with matched [O iii] +…

ABSTRACT Accretion discs in active galactic nuclei (AGNs) foster black hole (BH) formation, growth, and mergers. Stellar mass BHs migrate inwards under the influence of hydrodynamical torques unless they encounter a region where the torque flips sign. At these migration traps, BHs accumulate and merge via dynamical or gas-assisted interactions, producing high-frequency LIGO/Virgo/KAGRA (LVK) gravitational wave (GW) sources and potentially cutting off the supply of extreme mass ratio inspirals that would otherwise make low-frequency, LISA-band GWs. In this paper, we study the interplay between different types of migration torques, focusing especially on the ‘thermal torques’ generated by the thermal response of the AGN to embedded stellar-mass BHs that accrete through their own mini-discs. In contrast to previous work, we find that Type I torques cannot produce migration traps on their own, but thermal torques often do, particularly in low-mass AGN. The migration traps produced by thermal torques exist at much larger distances (∼103−5 gravitational radii) than do previously…

ABSTRACT The active galactic nucleus (AGN) disc has been proposed as a potential channel for the merger of binary black holes. The population of massive stars and black holes in AGN discs captured from the nuclei cluster plays a crucial role in determining the efficiency of binary formation and final merger rate within the AGN discs. In this paper, we investigate the capture process using analytical and numerical approaches. We discover a new constant integral of motion for one object’s capture process. Applying this result to the whole population of the nuclei cluster captured by the AGN disc, we find that the population of captured objects depends on the angular density and eccentricity distribution of the nuclei clusters and is effectively independent of the radial density profile of the nuclei cluster and disc models. An isotropic nuclei cluster with thermal eccentricity distribution predicts a captured profile dN/dr ∝ r−1/4. The captured objects are found to be dynamically crowded within the disc. Direct binary formation right after the capture would be promising, especially…

ABSTRACT Reliable indirect diagnostics of LyC photon escape from galaxies are required to understand which sources were the dominant contributors to reionization. While multiple LyC escape fraction (fesc) indicators have been proposed to trace favourable conditions for LyC leakage from the interstellar medium of low-redshift ‘analogue’ galaxies, it remains unclear whether these are applicable at high redshifts where LyC emission cannot be directly observed. Using a library of 14 120 mock spectra of star-forming galaxies with redshifts 4.64 ≤ z ≤ 10 from the SPHINX20 cosmological radiation hydrodynamics simulation, we develop a framework for the physics that leads to high fesc. We investigate LyC leakage from our galaxies based on the criteria that successful LyC escape diagnostics must (i) track a high-specific star formation rate, (ii) be sensitive to stellar population age in the range 3.5–10 Myr representing the times when supernova first explode to when LyC production significantly drops, and (iii) include a proxy for neutral gas content and gas density in the interstellar…

ABSTRACT This paper provides a comprehensive overview of how fitting of baryon acoustic oscillations (BAO) is carried out within the upcoming Dark Energy Spectroscopic Instrument’s (DESI) 2024 results using its DR1 data set, and the associated systematic error budget from theory and modelling of the BAO. We derive new results showing how non-linearities in the clustering of galaxies can cause potential biases in measurements of the isotropic ($\alpha _{\mathrm{iso}}$) and anisotropic ($\alpha _{\mathrm{ap}}$) BAO distance scales, and how these can be effectively removed with an appropriate choice of reconstruction algorithm. We then demonstrate how theory leads to a clear choice for how to model the BAO and develop, implement, and validate a new model for the remaining smooth-broad-band (i.e. without BAO) component of the galaxy clustering. Finally, we explore the impact of all remaining modelling choices on the BAO constraints from DESI using a suite of high-precision simulations, arriving at a set of best practices for DESI BAO fits, and an associated theory and modelling…

ABSTRACT We introduce Project GIBLE (Gas Is Better resoLved around galaxiEs), a suite of cosmological zoom-in simulations where gas in the circumgalactic medium (CGM) is preferentially simulated at ultra-high numerical resolution. Our initial sample consists of eight galaxies, all selected as Milky Way-like galaxies at z = 0 from the TNG50 simulation. Using the same galaxy formation model as IllustrisTNG, and the moving-mesh code arepo, we re-simulate each of these eight galaxies maintaining a resolution equivalent to TNG50-2 (mgas ∼ 8 × 105 M⊙). However, we use our super-Lagrangian refinement scheme to more finely resolve gas in the CGM around these galaxies. Our highest resolution runs achieve 512 times better mass resolution (∼103 M⊙). This corresponds to a median spatial resolution of ∼75 pc at 0.15 R200, c, which coarsens with increasing distance to ∼700 pc at the virial radius. We make predictions for the covering fractions of several observational tracers of multiphase CGM gas: H i, Mg ii, C iv, and O vii. We then study the impact of improved resolution on small scale…

ABSTRACT A key event in the history of the Milky Way is the formation of the bar. This event affects the subsequent structural and dynamical evolution of the entire Galaxy. When the bar formed, gas was likely rapidly funnelled to the centre of the Galaxy settling in a star-forming nuclear disc. The Milky Way bar formation can then be dated by considering the age distribution of the oldest stars in the formed nuclear stellar disc. In this highly obscured and crowded region, reliable age tracers are limited, but bright, high-amplitude Mira variables make useful age indicators as they follow a period–age relation. We fit dynamical models to the proper motions of a sample of Mira variables in the Milky Way’s nuclear stellar disc region. Weak evidence for inside-out growth and both radial and vertical dynamical heating with time of the nuclear stellar disc is presented, suggesting that the nuclear stellar disc is dynamically well-mixed. Furthermore, for Mira variables around a ∼350-d period, there is a clear transition from nuclear stellar disc-dominated kinematics to background…

ABSTRACT We present the census of H $\beta$ + [${\rm O\, {\small III}}$] $4960,5008\rm{\mathring{\rm\,A}}$ emitters at $6.8\lt z\lt 9.0$ from the JWST FRESCO survey over 124 arcmin$^2$ in the GOODS-North and GOODS-South fields. Our unbiased spectroscopic search results in 137 spectroscopically confirmed galaxies at $6.8\lt z\lt 9.0$ with observed [${\rm O\, {\small III}}$] fluxes $f_{\rm { [{\rm O\, {\small III}}]}}\gtrsim 1\times 10^{-18}\ \rm {ergs}\ \rm {s}^{-1} \ \rm {cm}^{-2}$. The rest-frame optical line ratios of the median stacked spectrum (median $M_{\rm {UV}}=-19.65^{+0.59}_{-1.05}$) indicate negligible dust attenuation, low metallicity ($12+\log (\rm {O/H})= 7.2-7.7$) and a high ionization parameter $\log _{10}U \simeq -2.5$. We find a factor $\times 1.3$ difference in the number density of $6.8\lt z\lt 9.0$ galaxies between GOODS-South and GOODS-North, which is caused by a single overdensity at $7.0\lt z\lt 7.2$ in GOODS-North. The bright end of the UV luminosity function of spectroscopically confirmed [${\rm O\, {\small III}}$] emitters is in good agreement with HST…

ABSTRACT The stellar initial mass function (IMF) in the early universe is essential to understand the formation of ancient galaxies. To this end, we conduct a series of long-term radiation hydrodynamic simulations following star cluster formation, varying the metallicity from Z/Z⊙ = 10−4 to 1. We particularly consider the effects of protostellar radiative feedback, which modify the exact shape of the IMF and determine the star formation efficiency (SFE), i.e. the ratio between the mass in stars and the initial gas mass in the parental cloud. Our results show that the IMF changes from a Salpeter-type to a top-heavy function as the metallicity decreases. When Z/Z⊙ ≲ 10−2, the IMF becomes log-flat and distinct from a Salpeter-like IMF. Stellar feedback is effective in shaping both the low- and high-mass ends of the IMF. Heating of dust grains by stellar radiation suppresses small-scale fragmentation and reduces the number of low-mass stars with M* ≲ 1 M⊙ at all metallicities. The ionizing radiation hinders the growth of massive stars, steepening the slope of the IMF at the high-mass…

ABSTRACT Measurements of neutron star mass and radius or tidal deformability deliver unique insight into the equation of state (EOS) of cold dense matter. EOS inference is very often done using generalized parametric or non-parametric models, which deliver no information on composition. In this paper, we consider a microscopic nuclear EOS model based on a field theoretical approach. We show that current measurements from NICER and gravitational wave observations constrain primarily the symmetric nuclear matter EOS. We then explore what could be delivered by measurements of mass and radius at the level anticipated for future large-area X-ray timing telescopes. These should be able to place very strong limits on the symmetric nuclear matter EOS, in addition to constraining the nuclear symmetry energy that determines the proton fraction inside the neutron star.

ABSTRACT The Dark Energy Spectroscopic Instrument (DESI) collaboration recently released its first year of data (DR1) on baryon acoustic oscillations (BAO) in galaxy, quasar, and Lyman-$\alpha$ forest tracers. When combined with cosmic microwave background (CMB) and Type Ia supernovae (SNIa) data, DESI BAO results suggest potential thawing behaviour in dark energy. Cosmological analyses utilize comoving distances along ($D_\mathrm{ H}$) and perpendicular to ($D_\mathrm{ M}$) the line of sight. Notably, there are $1\sim 2\sigma$ deviations in $D_\mathrm{ M}$ and $D_\mathrm{ H}$ from Planck cosmology values in the luminous red galaxies (LRG) bins LRG1 and LRG2.This study examines the role of LRG1 and LRG2 in diverging DESI 2024 BAO cosmology from Planck cosmology. We use angle-averaged distance $D_\mathrm{ V}$ and the ratio $F_{\rm AP}=D_\mathrm{ M}/D_\mathrm{ H}$, which are more directly related to the measured monopole and quadrupole components of the galaxy power spectrum or correlation function, instead of the officially adopted $D_\mathrm{ M}$ and $D_\mathrm{ H}$. This…

ABSTRACT JWST has revealed a large population of UV-bright galaxies at $z\gtrsim 10$ and possibly overly massive galaxies at $z\gtrsim 7$, challenging standard galaxy formation models in the ΛCDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $4\lesssim z \lesssim 10$ in both ΛCDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $z\simeq 12$ based on spectroscopically confirmed galaxies (eight galaxies at $z\simeq 11\!-\!13.5$) exhibit no tension with the benchmark model. Photometric constraints at $12 \lesssim z\lesssim 16$ can be fully explained within EDE via either moderately increased star-formation efficiencies ($\epsilon _{\ast}\sim 3\!-\!10\ \hbox{per cent}$ at $M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$) or enhanced UV variabilities ($\sigma _{\rm UV}\sim 0.8\!-\!1.3$ mag at $M_{\rm halo}\sim 10^{10.5}{\,…

ABSTRACT We employ Hubble data and Gaussian Processes in order to reconstruct the dynamical connection function in $f(Q)$ cosmology beyond the coincident gauge. In particular, there exist three branches of connections that satisfy the torsionless and curvatureless conditions, parametrized by a new dynamical function $\gamma$. We express the redshift dependence of $\gamma$ in terms of the $H(z)$ function and the $f(Q)$ form and parameters, and then we reconstruct it using 55 $H(z)$ observation data. First, we investigate the case where ordinary conservation law holds, and we reconstruct the $f(Q)$ function, which is very well described by a quadratic correction on top of symmetric teleparallel equivalent of general relativity. Proceeding to the general case, we consider two of the most studied $f(Q)$ models of the literature, namely the square-root and the exponential one. In both cases we reconstruct $\gamma (z)$, and we show that according to Akaike Information Criterion and Bayesian Information Criterion information criteria its inclusion is favoured compared to both $\Lambda$cold…

ABSTRACT Early JWST observations that targeted so-called double-break sources (attributed to Lyman and Balmer breaks at z > 7), reported a previously unknown population of very massive, evolved high-redshift galaxies. This surprising discovery led to a flurry of attempts to explain these objects’ unexpected existence including invoking alternatives to the standard Λ cold dark matter (ΛCDM) cosmological paradigm. To test these early results, we adopted the same double-break candidate galaxy selection criteria to search for such objects in the JWST images of the CAnadian Near Infrared Imager and Slitless Spectrograph (NIRISS) Unbiased Cluster Survey (CANUCS), and found a sample of 19 sources over five independent CANUCS fields that cover a total effective area of ∼60 arcmin2 at z ∼ 8. However, (1) our spectral energy distribution fits do not yield exceptionally high-stellar masses for our candidates, while (2) spectroscopy of five of the candidates shows that while all five are at high redshifts, their red colours are due to high-equivalent width emission lines in star-forming…

ABSTRACT Pulsar timing arrays (PTAs) are ensembles of regularly observed millisecond pulsars timed to high precision. Each pulsar in an array could be affected by a suite of noise processes, most of which are astrophysically motivated. Analysing them carefully can be used to understand these physical processes. However, the primary purpose of these experiments is to detect signals that are common to all pulsars, in particular signals associated with a stochastic gravitational wave background. To detect this, it is paramount to appropriately characterize other signals that may otherwise impact array sensitivity or cause a spurious detection. Here, we describe the second data release and first detailed noise analysis of the pulsars in the MeerKAT Pulsar Timing Array, comprising high-cadence and high-precision observations of 83 millisecond pulsars over 4.5 yr. We use this analysis to search for a common signal in the data, finding a process with an amplitude of $\log _{10}{A_{\mathrm{ CURN}}} = -14.25^{+0.21}_{-0.36}$ and spectral index $\gamma _\mathrm{CURN} = 3.60^{+1.31}_{-0.89}$.…

ABSTRACT We present rest-frame optical data of the $z\sim 4$ submillimetre galaxy GN20 obtained with the JWST Near Infrared Spectrograph (NIRSpec) in integral field spectroscopy mode. The H$\alpha$ emission is asymmetric and clumpy and extends over a projected distance of >15 kpc. To first order, the large-scale ionized gas kinematics are consistent with a turbulent ($\sigma \sim 90$ km s$^{-1}$), rotating disc ($v_{\rm rot}\sim 500$ km s$^{-1}$), congruent with previous studies of its molecular and ionized gas kinematics. However, we also find clear evidence for non-circular motions in the H$\alpha$ kinematics. We discuss their possible connection with various scenarios, such as external perturbations, accretion, or radial flows. In the centre of GN20, we find broad-line emission (full width at half-maximum $\sim 1000{-}2000$ km s$^{-1}$) in the H$\alpha$ + [N ii] complex, suggestive of fast, active galactic nucleus-driven winds or, alternatively, of the broad-line region of an active black hole. Elevated values of [N ii] $\lambda 6583$/H$\alpha \ \gt\ 0.4$ and of the Hα equivalent…

ABSTRACT Observations and high-resolution hydrodynamical simulations indicate that massive star clusters assemble hierarchically from subclusters with a universal power-law cluster mass function. We study the consequences of such assembly for the formation of intermediate-mass black holes (IMBHs) at low metallicities (Z = 0.01 Z⊙) with our updated N-body code bifrost based on the hierarchical fourth-order forward integrator. bifrost integrates few-body systems using secular and regularized techniques including post-Newtonian equations of motion up to order PN3.5 and gravitational-wave recoil kicks for BHs. Single stellar evolution is treated using the fast population synthesis code sevn. We evolve three cluster assembly regions with Ntot = 1.70–2.35 × 106 stars following a realistic IMF in ∼1000 subclusters for t = 50 Myr. IMBHs with masses up to m• ∼ 2200 M⊙ form rapidly mainly via the collapse of very massive stars (VMSs) assembled through repeated collisions of massive stars followed by growth through tidal disruption events and BH mergers. No IMBHs originate from the stars in…

ABSTRACT Stars move away from their birthplaces over time via a process known as radial migration, which blurs chemo–kinematic relations used for reconstructing the Milky Way (MW) formation history. To understand the true time evolution of the MW, one needs to take into account the effects of this process. We show that stellar birth radii can be derived directly from the data with minimum prior assumptions on the Galactic enrichment history. This is done by first recovering the time evolution of the stellar birth metallicity gradient, $\mathrm{ d}\mathrm{[Fe/H]}(R, \tau)/\mathrm{ d}R$, through its inverse relation to the metallicity range as a function of age today, allowing us to place any star with age and metallicity measurements back to its birthplace, R$_b$. Applying our method to a large high-precision data set of MW disc subgiant stars, we find a steepening of the birth metallicity gradient from 11 to 8 Gyr ago, which coincides with the time of the last massive merger, Gaia–Sausage–Enceladus (GSE). This transition appears to play a major role in shaping both the…

ABSTRACT Here, we present high spectral resolution ${\it JWST}$ NIRSpec observations of GN 42437, a low-mass (log(M$_\ast /{\rm M}_\odot)=7.9$), compact ($r_e \lt 500$pc), extreme starburst galaxy at $z=5.59$ with 13 emission-line detections. GN 42437 has a low metallicity (5–10 per cent Z$_\odot$) and its rest-frame H$\alpha$ equivalent width suggests nearly all of the observed stellar mass formed within the last 3 Myr. GN 42437 has an extraordinary 7$\sigma$ significant [Ne v] 3427 Å detection. The [Ne v] line has a rest-frame equivalent width of $11\pm 2$Å, [Ne v]/H$\alpha =0.04\pm 0.007$, [Ne v]/[Ne iii] 3870Å $= 0.26\pm 0.04$, and [Ne v]/He ii 4687Å $ = 1.2\pm 0.5$. Ionization from massive stars, shocks, or high-mass X-ray binaries cannot simultaneously produce these [Ne v] and other low-ionization line ratios. Reproducing the complete nebular structure requires both massive stars and accretion on to a black hole. We do not detect broad lines nor do the traditional diagnostics indicate that GN 42437 has an accreting black hole. Thus, the very high-ionization emission lines…

ABSTRACT The mass distribution in massive elliptical galaxies encodes their evolutionary history, thus providing an avenue to constrain the baryonic astrophysics in their evolution. The power-law assumption for the radial mass profile in ellipticals has been sufficient to describe several observables to the noise level, including strong lensing and stellar dynamics. In this paper, we quantitatively constrained any deviation, or the lack thereof, from the power-law mass profile in massive ellipticals through joint lensing–dynamics analysis of a large statistical sample with 77 galaxy–galaxy lens systems. We performed an improved and uniform lens modelling of these systems from archival Hubble Space Telescope imaging using the automated lens modelling pipeline dolphin. We combined the lens model posteriors with the stellar dynamics to constrain the deviation from the power law after accounting for the line-of-sight lensing effects, a first for analyses on galaxy–galaxy lenses. We find that the Sloan Lens ACS Survey lens galaxies with a mean redshift of 0.2 are consistent with the…

ABSTRACT We present a comprehensive study of the star formation histories of massive-quenched galaxies at $z=3$ in three semi-analytical models (Shark, gaea, Galform) and three cosmological hydrodynamical simulations (Eagle, IllustrisTNG, Simba). We study the predicted number density and stellar mass function of massive-quenched galaxies, their formation and quenching time-scales and star formation properties of their progenitors. Predictions are disparate in all these diagnostics, for instance: (i) some simulations reproduce the observed number density of very massive-quenched galaxies ($\gt 10^{11}\, \rm {\rm M}_{\odot }$) but underpredict the high density of intermediate-mass ones, while others fit well the lower masses but underpredict the higher ones; (ii) in most simulations, except for gaea and Eagle, most massive-quenched galaxies had starburst periods, with the most intense ones happening at $4\lt z\lt 5$; however, only in Shark and IllustrisTNG we do find a large number of progenitors with star formation rates $\gt 300\rm \, {\rm M}_{\odot }\, yr^{-1}$; (iii) quenching…

ABSTRACT We present simulation-based cosmological wcold dark matter (wCDM) inference using dark energy survey year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky pkdgrav3 dark matter simulations, with cosmological model parameters sampled with a mixed…

ABSTRACT Pulse profile modelling (PPM) is a comprehensive relativistic ray-tracing technique employed to determine the properties of neutron stars. In this study, we apply this technique to the Type I X-ray burster and accretion-powered millisecond pulsar XTE J1814−338, extracting its fundamental properties using PPM of its thermonuclear burst oscillations. Using data from its 2003 outburst, and a single uniform temperature hotspot model, we infer XTE J1814−338 to be located at a distance of $7.2^{+0.3}_{-0.4}$ kpc, with a mass of $1.21^{+0.05}_{-0.05}$ M$_\odot$ and an equatorial radius of $7.0^{+0.4}_{-0.4}$ km. Our results also offer insight into the time evolution of the hotspot but point to some potential shortcomings of the single uniform temperature hotspot model. We explore the implications of this result, including what we can learn about thermonuclear burst oscillation mechanisms and the importance of modelling the accretion contribution to the emission during the burst.

ABSTRACT We report results from 8 h of JWST/MIRI low resolution spectroscopic (LRS) monitoring directly followed by 7 h of JWST/NIRSpec prism spectroscopic monitoring of the benchmark binary brown dwarf WISE 1049AB, the closest, brightest brown dwarfs known. We find water, methane, and CO absorption features in both components, including the 3.3 μm methane absorption feature and a tentative detection of small grain ($\lt $ 1μm) silicate absorption at $\gt $8.5 μm in WISE 1049A. Both components vary significantly ($\gt 1~{{\rm per\ cent}}$), with WISE 1049B displaying larger variations than WISE 1049A. Using K-means clustering, we find three main transition points in wavelength for both components of the binary: (1) change in behaviour at $\sim$2.3 μm coincident with a CO absorption bandhead, (2) change in behaviour at 4.2 μm, close to the CO fundamental band at $\lambda \gt $ 4.4 µm, and (3) change in behaviour at 8.3–8.5 µm, potentially corresponding to silicate absorption. We interpret the light curves observed with both NIRSpec and MIRI as likely stemming from (1) a deep pressure…

ABSTRACT We investigate close encounters by stellar mass black holes (BHs) in the gaseous discs of active galactic nuclei (AGNs) as a potential formation channel of binary black holes (BBHs). We perform a series of 2D isothermal viscous hydrodynamical simulations within a shearing box prescription using the Eulerian grid code Athena++. We co-evolve the embedded BHs with the gas keeping track of the energetic dissipation and torquing of the BBH by gas gravitation and inertial forces. To probe the dependence of capture on the initial conditions, we discuss a suite of 345 simulations spanning BBH impact parameter (b) and local AGN disc density (ρ0). We identify a clear region in b − ρ0 space where gas assisted BBH capture is efficient. We find that the presence of gas leads to strong energetic dissipation during close encounters between unbound BHs, forming stably bound eccentric BBHs. We find that the gas dissipation during close encounters increases for systems with increased disc density and deeper periapsis passages rp, fitting a power law such that $\Delta E \propto \rho…

ABSTRACT 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 ($\sim 60~{{\ \rm per \, cent}}$ stay within 5 kpc). The azimuthal velocity vϕ also decreases between [Fe/H] = −1.0 and −2.0, but is constant for VMP stars (at ∼+40 km s−1). 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…

ABSTRACT The 1D power spectrum P1D of the Ly α forest provides important information about cosmological and astrophysical parameters, including constraints on warm dark matter models, the sum of the masses of the three neutrino species, and the thermal state of the intergalactic medium. We present the first measurement of P1D with the quadratic maximum likelihood estimator (QMLE) from the Dark Energy Spectroscopic Instrument (DESI) survey early data sample. This early sample of 54 600 quasars is already comparable in size to the largest previous studies, and we conduct a thorough investigation of numerous instrumental and analysis systematic errors to evaluate their impact on DESI data with QMLE. We demonstrate the excellent performance of the spectroscopic pipeline noise estimation and the impressive accuracy of the spectrograph resolution matrix with 2D image simulations of raw DESI images that we processed with the DESI spectroscopic pipeline. We also study metal line contamination and noise calibration systematics with quasar spectra on the red side of the Ly α emission line. In…

ABSTRACT Several studies have detected Lyman-alpha (Ly α) from bright ($M_{\small UV}\lesssim -21.5$) galaxies during the early stages of reionization despite the significantly neutral intergalactic medium. To explain these detections, it has been suggested that z > 7 Ly α emitters (LAEs) inhabit physical Mpc (pMpc)-scale ionized regions powered by overdensities of faint galaxies; however, systematic searches for these overdensities near LAEs have been challenging. Here, we use Cosmic Evolution Early Release Science JWST/Near Infrared Camera imaging to search for large-scale galaxy overdensities near two very ultraviolet (UV)-bright, z = 8.7 LAEs in the Extended Groth Strip (EGS) field. We colour select 27 z = 8.4–9.1 candidates, including the one LAE in the footprint (EGSY8p7). From spectral energy distribution models, we infer moderately faint UV luminosities ($-21.2\lesssim {M_{\small UV}}\lesssim -19.1$) and stellar masses of M* ≈ 107.5–8.8 M⊙. All are efficient ionizing agents ($\xi _{\text{ion}}^{*}\approx 10^{25.5-26.0}$ Hz erg−1) and are generally morphologically simple with…

ABSTRACT Given the sensitivity of the resonant Lyman $\alpha$ (Ly $\alpha $) transition to absorption by neutral hydrogen, observations of Ly $\alpha$ emitting galaxies (LAEs) have been widely used to probe the ionizing capabilities of reionization-era galaxies and their impact on the intergalactic medium (IGM). However, prior to JWST our understanding of the contribution of fainter sources and of ionized ‘bubbles’ at earlier stages of reionization remained uncertain. Here, we present the characterization of three exceptionally distant LAEs at $z \gt 8$, newly discovered by JWST/Near-Infrared Spectrograph in the JADES survey. These three similarly bright ($M_\text{UV} \approx -20 \, \mathrm{mag}$) LAEs exhibit small Ly $\alpha$ velocity offsets from the systemic redshift, $\Delta v_\rm{{Ly\alpha }} \lesssim 200 \, \mathrm{km \, s^{-1}}$, yet span a range of Ly $\alpha$ equivalent widths (15, 31, and $132 \, \mathring{\rm A}$). The former two show moderate Ly $\alpha$ escape fractions ($f_\rm{esc, {Ly\alpha }} \approx 10~{{\rm per\,cent}}$), whereas Ly $\alpha$ escapes remarkably…

ABSTRACT The atmospheres of small, close-in exoplanets are vulnerable to rapid mass loss during protoplanetary disc dispersal via a process referred to as ‘boil-off’, in which confining pressure from the local gas disc reduces, inducing atmospheric loss and contraction. We construct self-consistent models of planet evolution during gaseous core accretion and boil-off. As the surrounding disc gas dissipates, we find that planets lose mass via subsonic breeze outflows which allow causal contact to exist between disc and planet. Planets initially accrete of order $\sim 10~{{\ \rm per\ cent}}$ in atmospheric mass, however, boil-off can remove $\gtrsim 90~{{\ \rm per\ cent}}$ of this mass during disc dispersal. We show that a planet’s final atmospheric mass fraction is strongly dictated by the ratio of cooling time-scale to disc dispersal time-scale, as well as the planet’s core mass and equilibrium temperature. With contributions from core cooling and radioactivity, we show that core luminosity eventually leads to the transition from boil-off to core-powered mass loss. We find that…

ABSTRACT Lyman-alpha (Ly α) emission from galaxies is currently our most promising probe for constraining when and how reionization began, and thus when the first galaxies formed. At z > 7, the majority of galaxies detected with Ly α are in candidate overdensities. Here, we quantify the probability of these galaxies residing in large ionized bubbles. We create (1.6 Gpc)3 intergalactic medium (IGM) simulations: sufficient volume to robustly measure bubble size distributions around UV-bright galaxies and rare overdensities. We find ${M_{\small UV}}\lesssim -16$ galaxies and overdensities are ≳10–1000 × more likely to trace ionized bubbles compared to randomly selected positions. The brightest galaxies and strongest overdensities have bubble size distributions with highest characteristic size and least scatter. We compare two models: gradual reionization driven by numerous UV-faint galaxies versus rapid reionization by rarer brighter galaxies, producing larger bubbles at fixed neutral fraction. We demonstrate that recently observed z ∼ 7 overdensities are highly likely to trace large…

ABSTRACT We identify the progenitor star of SN 2023ixf in Messier 101 using Keck/NIRC2 adaptive optics imaging and pre-explosion Hubble Space Telescope (HST)/Advanced Camera for Surveys (ACS) images. The supernova, localized with diffraction spikes and high-precision astrometry, unambiguously coincides with a progenitor candidate of $m_\text{F814W}=24.87\pm 0.05$ (AB). Given its reported infrared excess and semiregular variability, we fit a time-dependent spectral energy distribution (SED) model of a dusty red supergiant (RSG) to a combined data set of HST optical, ground-based near-infrared, and Spitzer Infrared Array Camera (IRAC) [3.6], [4.5] photometry. The progenitor resembles an RSG of $T_\text{eff}=3488\pm 39$ K and $\log (L/\mathrm{L}_\odot)=5.15\pm 0.02$, with a $0.13\pm 0.01$ dex ($31.1\pm 1.7$ per cent) luminosity variation at a period of $P=1144.7\pm 4.8$ d, obscured by a dusty envelope of $\tau =2.92\pm 0.02$ at $1\, \mu \text{m}$ in optical depth (or $A_\text{V}=8.43\pm 0.11$ mag). The signatures match a post-main-sequence star of $18.2_{-0.6}^{+1.3}\,…

ABSTRACT For decades, the boundary of cosmic filaments has been a subject of debate. In this work, we determine the physically motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at $z=0$ in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from $z=4$ to 1, and then slightly increases. Moreover, we find that the radius of the filament depends on the length of the filament, the distance from the connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experienced rapid contraction before $z=1$, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have…

ABSTRACT Models of accretion discs surrounding active galactic nuclei (AGNs) find vast applications in high-energy astrophysics. The broad strategy is to parametrize some of the key disc properties such as gas density and temperature as a function of the radial coordinate from a given set of assumptions on the underlying physics. Two of the most popular approaches in this context were presented by Sirko & Goodman and Thompson et al. We present a critical reanalysis of these widely used models, detailing their assumptions and clarifying some steps in their derivation that were previously left unsaid. Our findings are implemented in the pAGN module for the Python programming language, which is the first public implementation of these accretion-disc models. We further apply pAGN to the evolution of stellar-mass black holes embedded in AGN discs, addressing the potential occurrence of migration traps.

ABSTRACT Recent observations with JWST have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the big bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox$^{\it HR}$, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments (FIRE) project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox$^{\it HR}$ re-simulates the cosmic volume ($L=22.1$ cMpc) of the original FIREbox run with eight times higher mass resolution ($m_{\rm b}\sim {}7800\, M_\odot$), but with identical physics, down to $z\sim {}6$. FIREbox$^{\it HR}$ predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at $z\sim {}6{\!-\!}14$, demonstrating that relatively high star formation…

ABSTRACT The distribution of mass in galaxy-scale strong gravitational lenses is often modelled as an elliptical power-law plus ‘external shear’, which notionally accounts for neighbouring galaxies and cosmic shear along our line of sight. A small amount of external shear could come from these sources, but we show that the vast majority does not. Except in a handful of rare systems, the best-fitting values do not correlate with independent measurements of line-of-sight shear: from weak lensing in 45 Hubble Space Telescope images, or in 50 mock images of lenses with complex distributions of mass. Instead, the best-fit external shear is aligned with the major or minor axis of 88 per cent of lens galaxies; and the amplitude of the external shear increases if that galaxy is discy. We conclude that ‘external shear’ attached to a power-law model is not physically meaningful, but a fudge to compensate for lack of model complexity. Since it biases other model parameters that are interpreted as physically meaningful in several science analyses (e.g. measuring galaxy evolution, dark matter…

ABSTRACT The MeerKAT International GHz Tiered Extragalactic Exploration Survey (MIGHTEE) is one of the large survey projects using the MeerKAT telescope, covering four fields that have a wealth of ancillary data available. We present Data Release 1 of the MIGHTEE continuum survey, releasing total intensity images and catalogues over $\sim$20 deg$^2$, across three fields at $\sim$1.2–1.3 GHz. This includes 4.2 deg$^2$ over the Cosmic Evolution Survey (COSMOS) field, 14.4 deg$^2$ over the XMM Large-Scale Structure (XMM-LSS) field and deeper imaging over 1.5 deg$^2$ of the Extended Chandra Deep Field South (CDFS). We release images at both a lower resolution (7–9 arcsec) and higher resolution ($\sim 5$ arcsec). These images have central rms sensitivities of $\sim 1.3-$2.7 $\mu$Jy beam$^{-1}$ ($\sim 1.2-$3.6 $\mu$Jy beam$^{-1}$) in the lower (higher) resolution images, respectively. We also release catalogues comprised of $\sim$144 000 ($\sim$114 000) sources using the lower (higher) resolution images. We compare the astrometry and flux-density calibration with the Early Science data in…

ABSTRACT We report observations of the highly active FRB 20240114A with MeerKAT using the ultra-high frequency (UHF; $544\!-\!1088$ MHz) and L band ($856\!-\!1712$ MHz) receivers. A total of 62 bursts were detected in coherent tied-array beams using the MeerTRAP real-time transient detection pipeline. We measure a structure-optimizing dispersion measure of $527.65\pm 0.01\, \text{pc}\, \text{cm}^{-3}$ using the brightest burst in the sample. We find the bursts of FRB 20240114A are generally detected in part of the broad-band of MeerKAT, $\sim 40~{{\ \rm per\ cent}}$ in the UHF and $\sim 30~{{\ \rm per\ cent}}$ in the L band, indicating the band limited nature. We analyse the fluence distribution of the 44 bursts detected at UHF, constraining the fluence completeness limit to $\sim 1$ Jy ms, above which the cumulative burst rate follows a power law $R (\gt F)\propto (F/1\, \text{Jy}\, \text{ms})^\gamma$ with $\gamma =-1.8\pm 0.2$. Using channelized telescope data captured in our transient buffer we localize FRB 20240114A in the image domain to RA = 21h27m39.86s, Dec. = +04d19m45.01s…

ABSTRACT Jetted astrophysical phenomena with black hole engines, including binary mergers, jetted tidal disruption events, and X-ray binaries, require a large-scale vertical magnetic field for efficient jet formation. However, a dynamo mechanism that could generate these crucial large-scale magnetic fields has not been identified and characterized. We have employed three-dimensional global general relativistic magnetohydrodynamical simulations of accretion discs to quantify, for the first time, a dynamo mechanism that generates large-scale magnetic fields. This dynamo mechanism primarily arises from the non-linear evolution of the magnetorotational instability (MRI). In this mechanism, large non-axisymmetric MRI-amplified shearing wave modes, mediated by the axisymmetric azimuthal magnetic field, generate and sustain the large-scale vertical magnetic field through their non-linear interactions. We identify the advection of magnetic loops as a crucial feature, transporting the large-scale vertical magnetic field from the outer regions to the inner regions of the accretion disc. This…

ABSTRACT The radius valley, i.e. a dearth of planets with radii between 1.5 and 2 Earth radii, provides insights into planetary formation and evolution. Using homogenously revised planetary parameters from Kepler 1-min short cadence light curves, we remodel transits of 72 small planets mostly orbiting low-mass stars, improving the precision and accuracy of planet parameters. By combining this sample with a similar sample of planets around higher mass stars, we determine the depth of the radius valley as a function of stellar mass. We find that the radius valley is shallower for low-mass stars compared to their higher mass counterparts. Upon comparison, we find that theoretical models of photoevaporation underpredict the number of planets observed inside the radius valley for low-mass stars: with decreasing stellar mass, the predicted fraction of planets inside the valley remains approximately constant whereas the observed fraction increases. We argue that this provides evidence for the presence of icy planets around low-mass stars. Alternatively, planets orbiting low-mass stars…

ABSTRACT The thermal continuum emission observed from accreting black holes across X-ray bands has the potential to be leveraged as a powerful probe of the mass and spin of the central black hole. The vast majority of existing ‘continuum fitting’ models neglect emission sourced at and within the innermost stable circular orbit (ISCO) of the black hole. Numerical simulations, however, find non-zero emission sourced from these regions. In this work, we extend existing techniques by including the emission sourced from within the plunging region, utilizing new analytical models that reproduce the properties of numerical accretion simulations. We show that in general the neglected intra-ISCO emission produces a hot-and-small quasi-blackbody component, but can also produce a weak power-law tail for more extreme parameter regions. A similar hot-and-small blackbody component has been added in by hand in an ad hoc manner to previous analyses of X-ray binary spectra. We show that the X-ray spectrum of MAXI J1820+070 in a soft-state outburst is extremely well described by a full Kerr black…

ABSTRACT We apply the hierarchical probabilistic spectral energy distribution (SED) model bayesn to analyse a sample of 475 type Ia supernovae (0.015 < z < 0.4) from Foundation, DES3YR and PS1MD to investigate the properties of dust in their host galaxies. We jointly infer the dust law RV population distributions at the SED level in high- and low-mass galaxies simultaneously with dust-independent, intrinsic differences. We find an intrinsic mass step of −0.049 ± 0.016 mag, at a significance of 3.1σ, when allowing for a constant intrinsic, achromatic magnitude offset. We additionally apply a model allowing for time- and wavelength-dependent intrinsic differences between SNe Ia in different mass bins, finding ∼2σ differences in magnitude and colour around peak and 4.5σ differences at later times. These intrinsic differences are inferred simultaneously with a difference in population mean RV of ∼2σ significance, demonstrating that both intrinsic and extrinsic differences may play a role in causing the host galaxy mass step. We also consider a model which allows the mean of the RV…

ABSTRACT With the advent of JWST, we can probe the rest-frame optical emission of galaxies at $z\gt 3$ with high sensitivity and spatial resolution, making it possible to accurately characterize red, optically faint galaxies and thus move towards a more complete census of the galaxy population at high redshifts. To this end, we present a sample of 148 massive, dusty galaxies from the JWST/Cosmic Evolution Early Release Science survey, colour-selected using solely JWST bands. With deep JWST/NIRCam data from 1.15 to 4.44 $\mu$m and ancillary HST/ACS and WFC3 data, we determine the physical properties of our sample using spectral energy distribution fitting with bagpipes. We demonstrate that our selection method efficiently identifies massive ($\mathrm{\langle \log {\it M}_\star /M_\odot \rangle \sim 10}$) and dusty ($\mathrm{\langle {\it A}_V\rangle \sim 2.7\ mag}$) sources, with a majority at $z\,\gt\, 3$ and predominantly lying on the galaxy main sequence. The main results of this work are the stellar mass functions (SMFs) of red, optically faint galaxies from redshifts between…

ABSTRACT We implemented Population III (Pop. III) star formation in mini-haloes within the Meraxes semi-analytic galaxy formation and reionization model, run on top of a N-body simulation with L = 10 h−1 cMpc with 20483 particles resolving all dark matter haloes down to the mini-haloes (∼105 M⊙). Our modelling includes the chemical evolution of the IGM, with metals released through supernova-driven bubbles that expand according to the Sedov–Taylor model. We found that SN-driven metal bubbles are generally small, with radii typically of 150 ckpc at z = 6. Hence, the majority of the first galaxies are likely enriched by their own star formation. However, as reionization progresses, the feedback effects from the UV background become more pronounced, leading to a halt in star formation in low-mass galaxies, after which external chemical enrichment becomes more relevant. We explore the sensitivity of the star formation rate density and stellar mass functions to the unknown values of free parameters. We also discuss the observability of Pop. III dominated systems with JWST, finding that…

ABSTRACT We present cosmological constraints derived from peak counts, minimum counts, and the angular power spectrum of the Subaru Hyper Suprime-Cam first-year (HSC Y1) weak lensing shear catalogue. Weak lensing peak and minimum counts contain non-Gaussian information and hence are complementary to the conventional two-point statistics in constraining cosmology. In this work, we forward-model the three summary statistics and their dependence on cosmology, using a suite of N-body simulations tailored to the HSC Y1 data. We investigate systematic and astrophysical effects including intrinsic alignments, baryon feedback, multiplicative bias, and photometric redshift uncertainties. We mitigate the impact of these systematics by applying cuts on angular scales, smoothing scales, signal-to-noise ratio bins, and tomographic redshift bins. By combining peaks, minima, and the power spectrum, assuming a flat-ΛCDM model, we obtain $S_{8} \equiv \sigma _8\sqrt{\Omega _m/0.3}= 0.810^{+0.022}_{-0.026}$, a 35 per cent tighter constraint than that obtained from the angular power spectrum alone.…

ABSTRACT We combine stellar orbits with the abundances of the heavy, r-process element europium and the light, $\alpha$-element, silicon to separate in situ and accreted populations in the Milky Way (MW) across all metallicities. At high orbital energy, the accretion-dominated halo shows elevated values of [Eu/Si], while at lower energies, where many of the stars were born in situ, the levels of [Eu/Si] are lower. These systematically different levels of [Eu/Si] in the MW and the accreted halo imply that the scatter in [Eu/$\alpha$] within a single galaxy is smaller than previously thought. At the lowest metallicities, we find that both accreted and in situ populations trend down in [Eu/Si], consistent with enrichment via neutron star mergers. Through compiling a large data set of abundances for 54 globular clusters (GCs), we show that differences in [Eu/Si] extend to populations of in situ/accreted GCs. We interpret this consistency as evidence that in r-process elements GCs trace the star formation history of their hosts, motivating their use as sub-Gyr timers of galactic…

ABSTRACT Spectra of the highest redshift galaxies taken with JWST are now allowing us to see into the heart of the reionization epoch. Many of these observed galaxies exhibit strong damping wing absorption redward of their Lyman-$\alpha$ emission. These observations have been used to measure the redshift evolution of the neutral fraction of the intergalactic medium and sizes of ionized bubbles. However, these estimates have been made using a simple analytic model for the intergalactic damping wing. We explore the recent observations with models of inhomogeneous reionization from the Sherwood–Relics simulation suite. We carry out a comparison between the damping wings calculated from the simulations and from the analytic model. We find that although the agreement is good on the red side of the Lyman-$\alpha$ emission, there is a discrepancy on the blue side due to residual neutral hydrogen present in the simulations, which saturates the intergalactic absorption. For this reason, we find that it is difficult to reproduce the claimed observations of large bubble sizes at $z \sim 7$,…

ABSTRACT Intensive broad-band reverberation mapping campaigns have shown that AGN variability is significantly more complex than expected from disc reverberation of the variable X-ray illumination. The UV/optical variability is highly correlated and lagged, with longer lags at longer wavelengths as predicted, but the observed time-scales are longer than expected. Worse, the UV/optical light curves are not well correlated with the X-rays, which should drive them. Instead, we consider an intrinsically variable accretion disc, where slow mass accretion rate fluctuations are generated in the optical-UV disc, propagating down to modulate intrinsically faster X-ray variability from the central regions. We match our model to Fairall 9, a well-studied AGN with L ∼ 0.1LEdd, where the spectrum is dominated by the UV/EUV. Our model produces light curves where the X-rays and UV have very different fast variability, yet are well correlated on longer time-scales, as observed. It predicts that the intrinsic variability has optical/UV leading the X-rays, but including reverberation of the variable…

ABSTRACT The population of hot Jupiters is extremely diverse, with large variations in their irradiation, period, gravity, and chemical composition. To understand the intrinsic planet diversity through the observed population level trends, we explore the a priori scatter in the population created by the different responses of atmospheric circulation to planetary parameters. We use the SPARC/MITgcm three-dimensional global circulation model to simulate 345 planets spanning a wide range of instellation, metallicity, gravity, and rotation periods typical for hot Jupiters, while differentiating between models with and without titanium oxide (TiO)/vanadium oxide (VO) in their atmosphere. We show that the combined effect of the planetary parameters leads to a large diversity in the ability of atmospheres to transport heat from day-side to night-side at a given equilibrium temperature. We further show that the hotspot offset is a non-monotonic function of planetary rotation period and explain our findings by a competition between the rotational and divergent parts of the circulation. As a…

ABSTRACT Observed accretion rates onto the Milky Way and other local spirals fall short of that required to sustain star formation for cosmological timescales. A potential avenue for this unseen accretion is a rotating inflow in the volume-filling hot phase ($\sim 10^6\, {\rm K}$) of the circumgalactic medium (CGM), as suggested by some cosmological simulations. Using hydrodynamic simulations and a new analytic solution valid in the slow-rotation limit, we show that a hot inflow spins up as it approaches the galaxy, while remaining hot, subsonic, and quasi-spherical. Within the radius of angular momentum support ($\sim 15\, {\rm kpc}$ for the Milky Way) the hot flow flattens into a disc geometry and then cools from $\sim 10^6$ to $\sim 10^4\, {\rm K}$ at the disc–halo interface. Cooling affects all hot gas, rather than just a subset of individual gas clouds, implying that accretion via hot inflows does not rely on local thermal instability in contrast with ‘precipitation’ models for galaxy accretion. Prior to cooling and accretion the inflow completes ≈tcool/tff radians of rotation,…

ABSTRACT The third data release of Gaia was the first to include orbital solutions assuming non-single stars. Here, we apply the astrometric triage technique of Shahaf et al. to identify binary star systems with companions that are not single main-sequence stars. Gaia’s synthetic photometry of these binaries is used to distinguish between systems likely to have white-dwarf companions and those that may be hierarchical triples. The study uncovered a population of nearly $3\, 200$ binaries, characterized by orbital separations on the order of an astronomical unit, in which the faint astrometric companion is probably a white dwarf. This sample increases the number of orbitally solved binary systems of this type by about two orders of magnitude. Remarkably, over 110 of these systems exhibit significant ultraviolet excess flux, confirming this classification and, in some cases, indicating their relatively young cooling ages. We show that the sample is not currently represented in synthetic binary populations, and is not easily reproduced by available binary population synthesis codes.…

ABSTRACT Compact binary mergers forming in star clusters may exhibit distinctive features that can be used to identify them among observed gravitational-wave sources. Such features likely depend on the host cluster structure and the physics of massive star evolution. Here, we dissect the population of compact binary mergers in the dragon-II simulation data base, a suite of 19 direct N-body models representing dense star clusters with up to 106 stars and $\lt 33~{{\ \rm per\ cent}}$ of stars in primordial binaries. We find a substantial population of black hole binary (BBH) mergers, some of them involving an intermediate-mass BH (IMBH), and a handful mergers involving a stellar BH and either a neutron star (NS) or a white dwarf (WD). Primordial binary mergers, $\sim 30~{{\ \rm per\ cent}}$ of the whole population, dominate ejected mergers. Dynamical mergers, instead, dominate the population of in-cluster mergers and are systematically heavier than primordial ones. Around 20 per cent of dragon-II mergers are eccentric in the Laser Interferometer Space Antenna (LISA) band and 5 per…

ABSTRACT We present the most comprehensive catalogue to date of Type I superluminous supernovae (SLSNe), a class of stripped-envelope supernovae (SNe) characterized by exceptionally high luminosities. We have compiled a sample of 262 SLSNe reported through 2022 December 31. We verified the spectroscopic classification of each SLSN and collated an exhaustive data set of ultraviolet, optical, and infrared photometry totalling over 30 000 photometric detections. Using these data, we derive observational parameters such as the peak absolute magnitudes, rise and decline time-scales, as well as bolometric luminosities, temperature, and photospheric radius evolution for all SLSNe. Additionally, we model all light curves using a hybrid model that includes contributions from both a magnetar central engine and the radioactive decay of $^{56}$Ni. We explore correlations among various physical and observational parameters, and recover the previously found relation between ejecta mass and magnetar spin, as well as the overall progenitor pre-explosion mass distribution with a peak at $\approx…

ABSTRACT The ability of Big Bang Nucleosynthesis theory to accurately predict the primordial abundances of helium and deuterium, as well as the baryon content of the Universe, is considered one of the most significant achievements in modern physics. In the present study, we consider two highly motivated hybrid $f(T)$ models and constrain them using the observations from the Big Bang Nucleosynthesis era. In addition, using late-time observations of Cosmic Chronometers and Gamma-Ray Bursts, the ranges of the model parameters are confined which are in good agreement with early time bounds. Subsequently, the common ranges obtained from the analysis for early and late time are summarized. Further, we verify the intermediating epochs by investigating the profiles of cosmographic parameters using the model parameter values from the common range. From this study, we find the considered teleparallel models are viable candidates to explain the primordial-intermediating-present epochs.

ABSTRACT Galaxy clusters are important probes for both cosmology and galaxy formation physics. We test the cosmological, hydrodynamical FLAMINGO (Full-hydro large-scale structure simulations with all-sky mapping for the interpretation of next generation observations) simulations by comparing to observations of the gaseous properties of clusters measured from X-ray observations. FLAMINGO contains unprecedented numbers of massive galaxy groups ($\gt 10^6$) and clusters ($\gt 10^5$) and includes variations in both cosmology and galaxy formation physics. We predict the evolution of cluster scaling relations as well as radial profiles of the temperature, density, pressure, entropy, and metallicity for different masses and redshifts. We show that the differences between volume-, and X-ray-weighting of particles in the simulations, and between cool- and non-cool-core samples, are similar in size as the differences between simulations for which the stellar and AGN (active galactic nucleus) feedback has been calibrated to produce significantly different gas fractions. Compared to thermally…

ABSTRACT Strongly lensed supernovae are rare and valuable probes of cosmology and astrophysics. Upcoming wide-field time-domain surveys, such as the Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST), are expected to discover an order-of-magnitude more lensed supernovae than have previously been observed. In this work, we investigate the cosmological prospects of lensed type Ia supernovae (SNIa) in LSST by quantifying the expected annual number of detections, the impact of stellar microlensing, follow-up feasibility, and how to best separate lensed and unlensed SNIa. We simulate SNIa lensed by galaxies, using the current LSST baseline v3.0 cadence, and find an expected number of 44 lensed SNIa detections per year. Microlensing effects by stars in the lensing galaxy are predicted to lower the lensed SNIa detections by ∼8 per cent. The lensed events can be separated from the unlensed ones by jointly considering their colours and peak magnitudes. We define a ‘gold sample’ of ∼10 lensed SNIa per year with time delay >10 d, >5 detections before light curve peak, and…

ABSTRACT JWST has recently sparked a new era of Lyα spectroscopy, delivering the first measurements of the Lyα escape fraction and velocity profile in typical galaxies at z ≃ 6−10. These observations offer new prospects for insight into the earliest stages of reionization. But to realize this potential, we need robust models of Lyα properties in galaxies at z ≃ 5−6 when the IGM is mostly ionized. Here, we use new JWST observations from the JADES and FRESCO surveys combined with VLT/MUSE and Keck/DEIMOS data to characterize statistical distributions of Lyα velocity offsets, escape fractions, and EWs in z ≃ 5−6 galaxies. We find that galaxies with large Lyα escape fractions (>0.2) are common at z ≃ 5−6, comprising 30 per cent of Lyman break selected samples. Comparing to literature studies, our census suggests that Lyα becomes more prevalent in the galaxy population towards higher redshift from z ∼ 3 to z ∼ 6, although we find that this evolution slows considerably between z ∼ 5 and z ∼ 6, consistent with modest attenuation from residual H i in the mostly ionized IGM at z ≃ 5−6. We…

ABSTRACT We present the results obtained from an X-ray timing study of the new black hole candidate Swift J1727.8−1613. The work is based on Hard X-ray Modulation Telescope (Insight-HXMT) observations carried out during the 2023 outburst. Prominent type-C low-frequency Quasi-periodic Oscillations (LFQPOs) are detected throughout the observations. With the substantial effective area of the Insight-HXMT at high energies, we examine the energy dependence of various parameters, including the centroid frequency, fractional rms, and phase lags of the type-C QPOs. Our findings align closely with those observed in high-inclination systems. During the initial stage of the outburst, a peaked noise component is also detected, the frequency of which is highly correlated with the LFQPO frequency, aligning with the Psaltis–Belloni–van der Klis relation. By assuming that the peaked noise originates from the precession of the accretion disc, the spin of this source can be constrained. Our results suggest that this source may possess a high spin.

ABSTRACT The ΛCDM paradigm successfully explains the large-scale structure of the Universe, but is less well constrained on subgalactic scales. Gravitational lens modelling has been used to measure the imprints of dark substructures on lensed arcs, testing the small-scale predictions of ΛCDM. However, the methods required for these tests are subject to degeneracies among the lens mass model and the source light profile. We present a case study of the unique compound gravitational lens SDSSJ0946+1006, wherein a dark, massive substructure has been detected, whose reported high concentration would be unlikely in a ΛCDM universe. For the first time, we model the first two background sources in both I- and U-band HST imaging, as well as VLT-MUSE emission line data for the most distant source. We recover a lensing perturber at a 5.9σ confidence level with mass $\log _{10}(M_\mathrm{sub}/{\rm M}_{\odot })=9.2^{+0.4}_{-0.1}$ and concentration $\log _{10}c=2.4^{+0.5}_{-0.3}$. The concentration is more consistent with CDM subhaloes than previously reported, and the mass is compatible with…

ABSTRACT Hot Jupiters present a unique opportunity for measuring how planet formation history shapes present-day atmospheric composition. However, due to the myriad pathways influencing composition, a well-constructed sample of planets is needed to determine whether formation history can be accurately traced back from atmospheric composition. To this end, the BOWIE-ALIGN survey (A spectral Light Investigation into hot gas Giant origiNs by the collaboration of Bristol, Oxford, Warwick, Imperial, Exeter, +) will compare the compositions of eight hot Jupiters around F stars, four with orbits aligned with the stellar rotation axis, and four misaligned. Using the alignment as an indicator for planets that underwent disc migration or high-eccentricity migration, one can determine whether migration history produces notable differences in composition between the two samples of planets. This paper describes the planet formation model that motivates our observing programme. Our model traces the accretion of chemical components from the gas and dust in the disc over a broad parameter space to…

ABSTRACT We present results from fully general relativistic (GR), three-dimensional (3D), neutrino-radiation magneto-hydrodynamic (MHD) simulations of stellar core collapse of a 20 M⊙ star with spectral neutrino transport. Our focus is to study the gravitational-wave (GW) signatures from the magnetorotationally (MR)-driven models. By parametrically changing the initial angular velocity and the strength of the magnetic fields in the core, we compute four models. Among our models, only those with cores having an initial magnetic field strength of 1012 G and rotation rates of 1 or 2 rad s−1 produce MHD jets. Seen from the direction perpendicular to the rotational axis, a characteristic waveform is obtained exhibiting a monotonic time increase in the wave amplitude. As previously identified, this stems from the propagating MHD outflows along the axis. We show that the GW amplitude from anisotropic neutrino emission becomes more than one order-of-magnitude bigger than that from the matter contribution, whereas seen from the rotational axis, both of the two components are in the same…

ABSTRACT Galaxy environment plays an important role in driving the transformation of galaxies from blue and star forming to red and quenched. Recent works have focused on the role of cosmic web filaments in galaxy evolution and have suggested that stellar mass segregation, quenching of star formation, and gas-stripping may occur within filaments. We study the relationship between distance to filament and the stellar mass, colour, and H i gas content of galaxies using data from the REsolved Spectroscopy of a Local VolumE survey and Environmental COntext (ECO) catalogue, two overlapping census-style, volume-complete surveys. We use the Discrete Persistence Structures Extractor to identify cosmic web filaments over the full ECO area. We find that galaxies close to filaments have higher stellar masses, in agreement with previous results. Controlling for stellar mass, we find that galaxies also have redder colours and are more gas poor closer to filaments. When accounting for group membership and halo mass, we find that these trends in colour and gas content are dominated by the…

ABSTRACT We use numerical simulations to model Gaia DR3 data with the aim of constraining the Milky Way (MW) bar and spiral structure parameters. We show that both the morphology and the velocity field in MW-like galactic disc models are strong functions of time, changing dramatically over a few tens of Myr. This suggests that by finding a good match to the observed radial velocity field, vR(x, y), we can constrain the bar-spiral orientation. Incorporating uncertainties into our models is necessary to match the data; most importantly, a heliocentric distance uncertainty above 10–15 per cent distorts the bar’s shape and vR quadrupole pattern morphology, and decreases its apparent angle with respect to the Sun-Galactocentric line. An excellent match to the Gaia DR3 vR(x, y) field is found for a simulation with a bar length Rb ≈ 3.6 kpc. We argue that the data are consistent with an MW bar as short as ∼3 kpc, for moderate strength inner disc spiral structure (A2/A0 ≈ 0.25) or, alternatively, with a bar length up to ∼5.2 kpc, provided that spiral arms are quite weak (A2/A0 ≈ 0.1), and…

ABSTRACT We report constraints on a variety of non-standard cosmological models using the full 5-yr photometrically classified type Ia supernova sample from the Dark Energy Survey (DES-SN5YR). Both Akaike Information Criterion (AIC) and Suspiciousness calculations find no strong evidence for or against any of the non-standard models we explore. When combined with external probes, the AIC and Suspiciousness agree that 11 of the 15 models are moderately preferred over Flat-$\Lambda$CDM suggesting additional flexibility in our cosmological models may be required beyond the cosmological constant. We also provide a detailed discussion of all cosmological assumptions that appear in the DES supernova cosmology analyses, evaluate their impact, and provide guidance on using the DES Hubble diagram to test non-standard models. An approximate cosmological model, used to perform bias corrections to the data holds the biggest potential for harbouring cosmological assumptions. We show that even if the approximate cosmological model is constructed with a matter density shifted by $\Delta \Omega…

ABSTRACT Over the past year, JWST has uncovered galaxies at record-breaking distances up to z ∼ 13. The JWST UNCOVER (ultra-deep NIRSpec and NIRcam observations before the epoch of reionization) program has obtained ultra-deep multiwavelength NIRCam imaging of the massive galaxy cluster A2744 over ∼45 arcmin2 down to ∼29.5 AB mag. Here, we present a robust ultraviolet (UV) luminosity function derived through lensing clusters at 9 < z < 12. Using comprehensive end-to-end simulations, we account for all lensing effects and systematic uncertainties in deriving both the amplification factors and the effective survey volume. Our results confirm the intriguing excess of UV-bright galaxies (MUV <−20 AB mag) previously reported at z > 9 in recent JWST studies. In particular, a double power-law (DPL) describes better the bright end of the luminosity function compared to the classical Schechter form. The number density of these bright galaxies is 10–100 times larger than theoretical predictions and previous findings based on Hubble Space Telescope (HST) observations. Additionally, we measure…

ABSTRACT Dark matter halo demographics and assembly histories are a manifestation of cosmological structure formation and have profound implications for the formation and evolution of galaxies. In particular, merger trees provide fundamental input for several modelling techniques, such as semi-analytic models (SAMs), sub-halo abundance matching (SHAM), and decorated halo occupation distribution models. Motivated by the new ultra-high-redshift (z ≳ 10) frontier enabled by JWST, we present a new suite of Gadget at Ultrahigh Redshift with Extra-Fine Timesteps (gureft) dark matter-only cosmological simulations that are carefully designed to capture halo merger histories and structural properties in the ultra-z universe. The simulation suite consists of four 10243-particle simulations with box sizes of 5, 15, 35, and 90 Mpc h−1, each with 170 snapshots stored between 40 ≥ z ≥ 6. With the unprecedented number of available snapshots and strategically chosen dynamic range covered by these boxes, gureft uncovers the emerging dark matter halo populations and their assembly histories in the…