Explore the vast range of titles available.

The cold dark matter (CDM) cosmological model has been remarkably successful in explaining cosmic structure over an enormous span of redshift, but it has faced persistent challenges from observations that probe the innermost regions of dark matter halos and the properties of the Milky Way’s dwarf galaxy satellites. We review the current observational and theoretical status of these “small-scale controversies.” Cosmological simulations that incorporate only gravity and collisionless CDM predict halos with abundant substructure and central densities that are too high to match constraints from galaxy dynamics. The solution could lie in baryonic physics: Recent numerical simulations and analytical models suggest that gravitational potential fluctuations tied to efficient supernova feedback can flatten the central cusps of halos in massive galaxies, and a combination of feedback and low star formation efficiency could explain why most of the dark matter subhalos orbiting the Milky Way do not host visible galaxies. However, it is not clear that this solution can work in the lowest mass galaxies, where discrepancies are observed. Alternatively, the small-scale conflicts could be evidence of more complex physics in the dark sector itself. For example, elastic scattering from strong dark matter self-interactions can alter predicted halo mass profiles, leading to good agreement with observations across a wide range of galaxy mass. Gravitational lensing and dynamical perturbations of tidal streams in the stellar halo provide evidence for an abundant population of low-mass subhalos in accord with CDM predictions. These observational approaches will get more powerful over the next few years.

The scale-free gravothermal fluid formalism has long proved effective in describing the evolution of self-interacting dark matter halos with a constant dark matter particle cross section. However, whether the gravothermal fluid solutions match numerical simulations for velocity-dependent cross-section scenarios remains untested. In this work, we provide a fast mapping method that relates the constant-cross-section gravothermal solution to models with arbitrary velocity dependence in the cross section. We show that the gravothermal solutions after mapping are in good agreement with Arepo N-body simulation results. We illustrate the power of this approach by applying this fast mapping method to a halo hosting a low-surface-brightness galaxy, UGC 128. We show that this fast mapping method can be used to constrain free parameters in a physically motivated cross-section model and illustrate the parameter space favored by the rotation curve measurement.

Testing the standard cosmological model (ΛCDM) at small scales is challenging. Galaxies that inhabit low-mass dark matter halos provide an ideal test bed for dark matter models by linking observational properties of galaxies at small scales (low mass, low velocity) to low-mass dark matter halos. However, the observed kinematics of these galaxies do not align with the kinematics of the dark matter halos predicted to host them, obscuring our understanding of the low-mass end of the galaxy–halo connection. We use deep H i observations of low-mass galaxies at high spectral resolution in combination with cosmological simulations of dwarf galaxies to better understand the connection between dwarf galaxy kinematics and low-mass halos. Specifically, we use H i line widths to directly compare to the maximum velocities in a dark matter halo and find that each deeper measurement approaches the expected one-to-one relationship between the observed kinematics and the predicted kinematics in ΛCDM. We also measure baryonic masses and place these on the baryonic Tully–Fisher relation (BTFR). Again, our deepest measurements approach the theoretical predictions for the low-mass end of this relation, a significant improvement on similar measurements based on line widths measured at 50% and 20% of the peak. Our data also hint at the rollover in the BTFR predicted by hydrodynamical simulations of ΛCDM for low-mass galaxies.

The Sun is a bright gamma-ray source due to hadronic cosmic-ray interactions with solar gas. While it is known that incoming cosmic rays must generally first be reflected by solar magnetic fields to produce outgoing gamma rays, theoretical models have yet to reproduce the observed spectra. We introduce a simplified model of the solar magnetic fields that captures the main elements relevant to gamma-ray production. These are a flux tube, representing the network elements, and a flux sheet, representing the intergranular sheets. Both the tube and sheet have a horizontal size of order 100 km and serve as sites where cosmic rays are reflected and gamma rays are produced. While our simplified double-structure model does not capture all the complexities of the solar-surface magnetic fields, such as Alfvén turbulence from wave interactions or magnetic fluctuations from convection motions, it improves on previous models by reasonably producing both the hard spectrum seen by Fermi Large Area Telescope at 1–200 GeV and the considerably softer spectrum seen by the High Altitude Water Cherenkov Observatory (HAWC) at near 103 GeV. We show that lower-energy (≲10 GeV) gamma rays are primarily produced in the network elements and higher-energy (≳few × 10 GeV) gamma rays in the intergranular sheets. Notably, the spectrum softening observed by HAWC results from the limited effectiveness of capturing and reflecting ∼104 GeV cosmic rays by the finite-sized intergranular sheets. Our study is important for understanding cosmic-ray transport in the solar atmosphere and will lead to insights into small-scale magnetic fields at the photosphere.

The solar disk is a continuous source of GeV–TeV gamma rays. The emission is thought to originate from hadronic Galactic cosmic rays (GCRs) interacting with the gas in the photosphere and uppermost convection zone after being reflected by solar magnetic fields. Despite this general understanding, existing theoretical models have yet to match observational data. At the photosphere and the uppermost convection zone, granular convection drives a multiscale magnetic field, forming a larger-scale filamentary structure while also generating turbulence-scale Alfvén wave turbulence. Here, we demonstrate that the larger-scale filamentary field shapes the overall gamma-ray emission spectrum, and the Alfvén wave turbulence is critical for further suppressing the gamma-ray emission spectrum below ∼100 GeV. For a standard Alfvén wave turbulence level, our model’s predicted spectrum slope from 1 GeV to 1 TeV is in excellent agreement with observations from the Fermi Large Area Telescope and HAWC, an important achievement. The predicted absolute flux is a factor of 2–5 lower than the observed data; we outline future directions to resolve this discrepancy. The key contribution of our work is providing a new theoretical framework for using solar disk gamma-ray observations to probe hadronic GCR transport in the lower solar atmosphere.

The observed rotation curves of dwarf galaxies exhibit significant diversity at fixed halo mass, challenging galaxy formation within the cold dark matter (CDM) model. Previous cosmological galaxy formation simulations with baryonic physics fail to reproduce the full diversity of rotation curves, suggesting that there is a flaw in baryonic feedback models, observational bias, or that an alternative to CDM must be invoked. In this work, we use the Marvelous Massive Dwarf zoom-in simulations, a suite of high-resolution dwarf simulations with M200 ∼ 1010–1011 M⊙and M* ∼ 107–109 M⊙, designed to target the mass range where the galaxy rotation curve diversity is maximized, i.e., between Vmax∼70 and 100 km s−1. We add to this a set of low-mass galaxies from the Marvel Dwarf Zoom Volumes to extend the galaxy mass range to lower values. Our fiducial star formation and feedback models produce simulated dwarfs with a broader range of rotation curve shapes, similar to observations. These simulations both create dark matter cores via baryonic feedback, reproducing the slower-rising rotation curves, while also allowing for compact galaxies and steeply rising rotation curves. Our simulated dwarfs also reproduce the observed size–M* relation, including scatter, producing both extended and compact dwarfs for the first time in simulated field dwarfs. However, the slowly rising and high baryon mass fraction, as well as the steeply rising and low baryon mass fraction, remain missing. We explore star formation and feedback models and conclude that previous simulations may have had feedback that was too strong to produce compact dwarfs.

We present a spectroscopic analysis of a broad-line active galactic nucleus (AGN) selected as a little red dot (LRD) at z = 5.3 behind the Bullet cluster (Bz5.3) based on JWST/NIRCam and NIRSpec data. The detection of strong Fe ii, O i, and Ca ii triplet emission lines, along with the evidence of broad Balmer lines, provides strong evidence of a broad-line region (BLR) and an accreting supermassive black hole. Notably, we report the first detection of the λ1304 bump (i.e., a blend of O i λ1304 and Si ii) at high redshift, a feature commonly seen in local AGNs but not yet reported in the early Universe. The O i λ1304/λ8446 photon ratio provides an independent measurement of dust attenuation in galaxies. In Bz5.3, this ratio is highly suppressed (0.1–0.3), implying significant internal dust extinction, with estimated dust attenuation AV ∼ 0.4–1.0. We identify Lyβ fluorescence as the dominant excitation mechanism of the low-ionization lines, with additional contributions from collisional excitation. High O i λ8446 equivalent width and weak O i λ7774 support this interpretation. The detection of iron emission, whether from broad permitted or narrow forbidden lines, supports the presence of a stratified BLR, as also recently proposed in local LRDs. Photoionization modeling of O i λ8446 and Ca ii further suggests the coexistence of multiple gas phases with distinct densities and ionization states, highlighting the complexity of the BLR. Bz5.3 thus offers a rare window into early AGN activity and BLR physics at early times.

Galaxy–galaxy strong gravitational lenses can constrain dark matter models and the Lambda cold dark matter cosmological paradigm at subgalactic scales. Currently, there is a dearth of images of these rare systems with high signal-to-noise ratio (SNR) and angular resolution. The Nancy Grace Roman Space Telescope (hereafter Roman), scheduled for launch in late 2026, will play a transformative role in strong-lensing science with its planned wide-field surveys. With its remarkable 0.281 square degree field of view and diffraction-limited angular resolution of ~0 .″ 1, Roman is uniquely suited to characterizing dark matter substructure from a robust population of strong lenses. We present a yield simulation of detectable strong lenses in Roman’s planned High Latitude Wide Area Survey (HLWAS). We simulate a population of galaxy–galaxy strong lenses across cosmic time with cold dark matter subhalo populations, select those detectable in the HLWAS, and generate simulated images accounting for realistic Wide Field Instrument detector effects. For a fiducial case of single 146 s exposures, we predict around 160,000 detectable strong lenses in the HLWAS, of which about 500 will have sufficient SNR to be amenable to detailed substructure characterization. We investigate the effect of variation of the point-spread function across Roman’s field of view on detecting individual subhalos and the suppression of the subhalo mass function at low masses. Our simulation products are available to support strong-lens science with Roman, such as training neural networks and validating dark matter substructure analysis pipelines.

Stellar streams in the Milky Way are promising detectors of low-mass dark matter (DM) subhalos predicted by ΛCDM. Passing subhalos induce perturbations in streams that indicate the presence of the subhalos. Understanding how known DM-dominated satellites impact streams is a crucial step toward using stream perturbations to constrain the properties of dark perturbers. Here, we cross-match a Gaia Early Data Release 3 and SEGUE member catalog of the Cetus-Palca stream (CPS) with H3 for additional radial velocity measurements and fit the orbit of the CPS using this six-dimensional (6D) data. We demonstrate for the first time that the ultra-faint dwarf Segue 2 had a recent (77 ± 5 Myr ago) close flyby (within the stream's 2σ width) with the CPS. This interaction enables constraints on Segue 2’s mass and density profile at larger radii ( O(1) kpc) than are probed by its stars ( O(10) pc). While Segue 2 is not expected to strongly affect the portion of the stream covered by our 6D data, we predict that if Segue 2’s mass within ∼ 6 kpc is 5 × 109 M⊙, the CPS's velocity dispersion will be ∼ 40 km s−1 larger at ϕ1 > 20° than at ϕ1 < 0°. If no such heating is detected, Segue 2’s mass cannot exceed 109 M⊙ within ∼ 6 kpc. The proper motion distribution of the CPS near the impact site is mildly sensitive to the shape of Segue 2’s density profile. This study presents a critical test for frameworks designed to constrain properties of dark subhalos from stream perturbations.

Using medium-band imaging from the newly released Merian Survey, we conduct a nonparametric morphological analysis of Hα emission maps and stellar continua for a sample of galaxies with 8≲log(M⋆/M⊙)<10.3 at 0.064 < z < 0.1. We present a novel method for estimating the stellar continuum emission through the Merian Survey’s N708 medium-band filter, which we use to measure Hα emission and produce Hα maps for our sample of galaxies with seven-band Merian photometry and available spectroscopy. We measure nonparametric morphological statistics for the Hα and stellar continuum images, explore how the morphology of the Hα differs from the continuum, and investigate how the parameters evolve with the galaxies’ physical properties. In agreement with previous results for more massive galaxies, we find that the asymmetry of the stellar continuum increases with specific star formation rate (sSFR), and we extend the trend to lower masses, also showing that it holds for the asymmetry of the Hα emission. We find that the lowest-mass galaxies with the highest sSFR have Hα emission that is consistently heterogeneous and compact, while the less active galaxies in this mass range have Hα emission that appears diffuse. At higher masses, our data do not span a sufficient range in sSFR to evaluate whether similar trends apply. We conclude that high sSFRs in low-mass galaxies likely result from dynamical instabilities that compress a galaxy’s molecular gas to a dense region near the center.