Beecroft Institute for Particle Astrophysics and Cosmology

Abstract:

Radio continuum emission provides a unique opportunity to study star formation unbiased by dust obscuration. However, if radio observations are to be used to accurately trace star formation to high redshifts, it is crucial that the physical processes that affect the radio emission from star-forming galaxies are well understood. While inverse Compton (IC) losses from the cosmic microwave background (CMB) are negligible in the local universe, the rapid increase in the strength of the CMB energy density with redshift [] means that this effect becomes increasingly important at . Using a sample of high-redshift () Lyman-break galaxies selected in the rest-frame ultraviolet (UV), we have stacked radio observations from the MIGHTEE survey to estimate their 1.4-GHz flux densities. We find that for a given rest-frame UV magnitude, the 1.4-GHz flux density and luminosity decrease with redshift. We compare these results to the theoretical predicted effect of energy losses due to IC scattering off the CMB, and find that the observed decrease is consistent with this explanation. We discuss other possible causes for the observed decrease in radio flux density with redshift at a given UV magnitude, such as a top-heavy initial mass function at high redshift or an evolution of the dust properties, but suggest that IC scattering is the most compelling explanation.

Abstract:

With a new probabilistic technique for sampling interstellar object (ISO) orbits with high efficiency, we assess the observability of ISOs under a realistic cadence for the upcoming Vera C. Rubin Observatory’s Legacy Survey of Space and Time (LSST). Using the Ōtautahi–Oxford population model, we show that there will be complex on-sky structure in the pattern of direction and velocity revealed by the detected ISO population, with the expected enhanced northern flux complicating efforts to derive population parameters from the LSST’s predominately southern footprint. For reasonable luminosity functions with slopes of 2.5 ≤ qs ≤ 4.0, the most discoverable ISOs have Hr ≃ 14.6−20.7. The slope of the luminosity function of ISOs will be relatively quickly constrained by the characteristics of the LSST detected population, such as the distributions of perihelia, velocity at infinity, and discovery circumstances. Discoveries are evenly split around their perihelion passage and are biased to lower velocities. After their discovery by LSST, it will be rare for ISOs to be visible for less than a month; most will have mr ≤ 23 for months, and the window for spectroscopic characterization could be as long as 2 yr. While these probabilistic assessments are robust against model or spatial density refinements that change the absolute numbers of ISO discoveries, our simulations predict a yield of 6–51 asteroidal ISOs, which is similar to previous works and demonstrates the validity of our new methods.

Abstract:

The discovery of the third interstellar object (ISO), 3I/ATLAS (“3I”), provides a rare chance to directly observe a small body from another solar system. Studying its chemistry and dynamics will add to our understanding of how the processes of planetesimal formation and evolution happen across the Milky Way’s disk, and how such objects respond to the Milky Way’s potential. In this Letter, we present a first assessment of 3I in the context of the Ōtautahi–Oxford model, which uses data from Gaia in conjunction with models of protoplanetary disk chemistry and Galactic dynamics to predict the properties of the ISO population. The model shows that both the velocity and radiant of 3I are within the expected range. Its velocity predicts an age of over 7.6 Gyr and a high water mass fraction, which may become observable shortly. We also conclude that it is very unlikely that 3I shares an origin with either of the previous two ISO detections.

Abstract:

We present a detailed kinematic study of a sample of 32 massive (9.5 ⩽ log(M */M⊙) ⩽ 10.9) main sequence star-forming galaxies (MS SFGs) at 4 < z < 6 from the ALMA-CRISTAL programme. The data consist of deep (up to 15 hr observing time per target), high-resolution (∼1 kpc) ALMA observations of [C II]158 μm line emission. This dataset allowed us to carry out the first systematic, kiloparsec-scale (kpc-scale) characterisation of the kinematics nature of typical massive SFGs at these epochs. We find that ∼50% of the sample are disk-like, with a number of galaxies located in systems of multiple components. Kinematic modelling reveals these main sequence disks exhibit high-velocity dispersions (σ 0), with a median disk velocity dispersion of ∼70 km s⁻¹ and V rot/σ 0 ∼ 2, which is consistent with dominant gravity driving. The elevated disk dispersions are in line with the predicted evolution based on Toomre theory and the extrapolated trends from z ∼ 0–2.5 MS star-forming disks. The inferred dark matter (DM) mass fraction within the effective radius f DM(< R e) for the disk systems decreases with the central baryonic mass surface density. This is consistent with the trend reported by kinematic studies at z ≲ 3; roughly half the disks display f DM(< R e)≲ 30%. The CRISTAL sample of massive MS SFGs provides a reference of the kinematics of a representative population and extends the view onto typical galaxies beyond previous kpc-scale studies at z ≲ 3.

Abstract:

Context. Baryonic physics has a considerable impact on the distribution of matter in our Universe on scales probed by current and future cosmological surveys, acting as a key systematic in such analyses. Aims. We seek simple symbolic parametrisations for the impact of baryonic physics on the matter power spectrum for a range of physically motivated models, as a function of wavenumber, redshift, cosmology, and parameters controlling the baryonic feedback. Methods. We used symbolic regression to construct analytic approximations for the ratio of the matter power spectrum in the presence of baryons to that without such effects. We obtained separate functions of each of four distinct sub-grid prescriptions of baryonic physics from the CAMELS suite of hydrodynamical simulations (Astrid, IllustrisTNG, SIMBA, and Swift-EAGLE) as well as for a baryonification algorithm. We also provide functions that describe the uncertainty on these predictions, due to both the stochastic nature of baryonic physics and the errors on our fits. Results. The error on our approximations to the hydrodynamical simulations is comparable to the sample variance estimated through varying initial conditions, and our baryonification expression has a root mean squared error of better than one percent, although this increases on small scales. These errors are comparable to those of previous numerical emulators for these models. Our expressions are enforced to have the physically correct behaviour on large scales and at high redshift. Due to their analytic form, we are able to directly interpret the impact of varying cosmology and feedback parameters, and we can identify parameters that have little to no effect. Conlcusions. Each function is based on a different implementation of baryonic physics, and can therefore be used to discriminate between these models when applied to real data. We provide a publicly available code for all symbolic approximations found.