Augmenting astrophysical scaling relations with machine learning: Application to reducing the Sunyaev-Zeldovich flux-mass scatter.

Complex astrophysical systems often exhibit low-scatter relations between observable properties (e.g., luminosity, velocity dispersion, oscillation period). These scaling relations illuminate the underlying physics, and can provide observational tools for estimating masses and distances. Machine learning can provide a fast and systematic way to search for new scaling relations (or for simple extensions to existing relations) in abstract high-dimensional parameter spaces. We use a machine learning tool called symbolic regression (SR), which models patterns in a dataset in the form of analytic equations. We focus on the Sunyaev-Zeldovich flux-cluster mass relation (YSZ - M), the scatter in which affects inference of cosmological parameters from cluster abundance data. Using SR on the data from the IllustrisTNG hydrodynamical simulation, we find a new proxy for cluster mass which combines YSZ and concentration of ionized gas (cgas): M ∝ Yconc3/5 ≡ YSZ3/5(1 - A cgas). Yconc reduces the scatter in the predicted M by ∼20 - 30% for large clusters (M ≳ 1014 h-1 M⊙), as compared to using just YSZ. We show that the dependence on cgas is linked to cores of clusters exhibiting larger scatter than their outskirts. Finally, we test Yconc on clusters from CAMELS simulations and show that Yconc is robust against variations in cosmology, subgrid physics, and cosmic variance. Our results and methodology can be useful for accurate multiwavelength cluster mass estimation from upcoming CMB and X-ray surveys like ACT, SO, eROSITA and CMB-S4.

Anisotropic satellite galaxy quenching modulated by black hole activity.

The evolution of satellite galaxies is shaped by their constant interaction with the circumgalactic medium surrounding central galaxies, which in turn may be affected by gas and energy ejected from the central supermassive black hole1-6. The nature of such a coupling between black holes and galaxies is, however, much debated7-9 and observational evidence remains scarce10,11. Here we report an analysis of archival data on 124,163 satellite galaxies in the potential wells of 29,631 dark matter halos with masses between 1012 and 1014 solar masses. We find that quenched satellite galaxies are relatively less frequent along the minor axis of their central galaxies. This observation might appear counterintuitive given that black hole activity is expected to eject mass and energy preferentially in the direction of the minor axis of the host galaxy. We show, however, that the observed anisotropic signal results precisely from the ejective nature of black hole feedback in massive halos, as outflows powered by active galactic nuclei clear out the circumgalactic medium, reducing the ram pressure and thus preserving star formation in satellite galaxies. This interpretation is supported by the IllustrisTNG suite of cosmological numerical simulations, even though the model's sub-grid implementation of black hole feedback is effectively isotropic12.

First Star-Forming Structures in Fuzzy Cosmic Filaments.

In hierarchical models of structure formation, the first galaxies form in low-mass dark matter potential wells, probing the behavior of dark matter on kiloparsec scales. Even though these objects are below the detection threshold of current telescopes, future missions will open an observational window into this emergent world. In this Letter, we investigate how the first galaxies are assembled in a "fuzzy" dark matter (FDM) cosmology where dark matter is an ultralight ∼10^{-22} eV boson and the primordial stars are expected to form along dense dark matter filaments. Using a first-of-its-kind cosmological hydrodynamical simulation, we explore the interplay between baryonic physics and unique wavelike features inherent to FDM. In our simulation, the dark matter filaments show coherent interference patterns on the boson de Broglie scale and develop cylindrical solitonlike cores, which are unstable under gravity and collapse into kiloparsec-scale spherical solitons. Features of the dark matter distribution are largely unaffected by the baryonic feedback. On the contrary, the distributions of gas and stars, which do form along the entire filament, exhibit central cores imprinted by dark matter-a smoking gun signature of FDM.

Galaxy formation with BECDM - I. Turbulence and relaxation of idealized haloes.

We present a theoretical analysis of some unexplored aspects of relaxed Bose-Einstein condensate dark matter (BECDM) haloes. This type of ultralight bosonic scalar field dark matter is a viable alternative to the standard cold dark matter (CDM) paradigm, as it makes the same large-scale predictions as CDM and potentially overcomes CDM's small-scale problems via a galaxy-scale de Broglie wavelength. We simulate BECDM halo formation through mergers, evolved under the Schrödinger-Poisson equations. The formed haloes consist of a soliton core supported against gravitational collapse by the quantum pressure tensor and an asymptotic r-3 NFW-like profile. We find a fundamental relation of the core-to-halo mass with the dimensionless invariant Ξ ≡ |E|/M3/(Gm/h)2 or Mc/M ≃ 2.6Ξ1/3, linking the soliton to global halo properties. For r ≥ 3.5 rc core radii, we find equipartition between potential, classical kinetic and quantum gradient energies. The haloes also exhibit a conspicuous turbulent behaviour driven by the continuous reconnection of vortex lines due to wave interference. We analyse the turbulence 1D velocity power spectrum and find a k-1.1 power law. This suggests that the vorticity in BECDM haloes is homogeneous, similar to thermally-driven counterflow BEC systems from condensed matter physics, in contrast to a k-5/3 Kolmogorov power law seen in mechanically-driven quantum systems. The mode where the power spectrum peaks is approximately the soliton width, implying that the soliton-sized granules carry most of the turbulent energy in BECDM haloes.

The formation of the first stars and galaxies.

Observations made using large ground-based and space-borne telescopes have probed cosmic history from the present day to a time when the Universe was less than one-tenth of its present age. Earlier still lies the remaining frontier, where the first stars, galaxies and massive black holes formed. They fundamentally transformed the early Universe by endowing it with the first sources of light and chemical elements beyond the primordial hydrogen and helium produced in the Big Bang. The interplay of theory and upcoming observations promises to answer the key open questions in this emerging field.

Protostar formation in the early universe.

The nature of the first generation of stars in the universe remains largely unknown. Observations imply the existence of massive primordial stars early in the history of the universe, and the standard theory for the growth of cosmic structure predicts that structures grow hierarchically through gravitational instability. We have developed an ab initio computer simulation of the formation of primordial stars that follows the relevant atomic and molecular processes in a primordial gas in an expanding universe. The results show that primeval density fluctuations left over from the Big Bang can drive the formation of a tiny protostar with a mass 1% that of the Sun. The protostar is a seed for the subsequent formation of a massive primordial star.

Numerical simulations unravel the cosmic web.

The universe is permeated by a network of filaments, sheets, and knots collectively forming a "cosmic web." The discovery of the cosmic web, especially through its signature of absorption of light from distant sources by neutral hydrogen in the intervening intergalactic medium, exemplifies the interplay between theory and experiment that drives science and is one of the great examples in which numerical simulations have played a key and decisive role. We recount the milestones in our understanding of cosmic structure; summarize its impact on astronomy, cosmology, and physics; and look ahead by outlining the challenges faced as we prepare to probe the cosmic web at new wavelengths.

Energy input from quasars regulates the growth and activity of black holes and their host galaxies.

In the early Universe, while galaxies were still forming, black holes as massive as a billion solar masses powered quasars. Supermassive black holes are found at the centres of most galaxies today, where their masses are related to the velocity dispersions of stars in their host galaxies and hence to the mass of the central bulge of the galaxy. This suggests a link between the growth of the black holes and their host galaxies, which has indeed been assumed for a number of years. But the origin of the observed relation between black hole mass and stellar velocity dispersion, and its connection with the evolution of galaxies, have remained unclear. Here we report simulations that simultaneously follow star formation and the growth of black holes during galaxy-galaxy collisions. We find that, in addition to generating a burst of star formation, a merger leads to strong inflows that feed gas to the supermassive black hole and thereby power the quasar. The energy released by the quasar expels enough gas to quench both star formation and further black hole growth. This determines the lifetime of the quasar phase (approaching 100 million years) and explains the relationship between the black hole mass and the stellar velocity dispersion.

Brownian motion in gravitationally interacting systems.

We derive a model that describes the dynamics of a Brownian particle, such as a massive black hole, in a stellar system dominated by gravitational forces, and examine whether it achieves a state of equipartition of kinetic energy with the stars. This problem has been considered before only for stellar systems with an isothermal Maxwellian distribution of velocities; here we study other examples and confirm our calculations with N-body simulations. We show that in certain cases the black hole's steady state kinetic energy can be very far from equipartition.