Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Galaxy dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 307 records · Page 17

Memoirs of Mass Accretion: Probing the Edges of Intracluster Light in Simulated Galaxy Clusters

The diffuse starlight extending throughout massive galaxy clusters, known as intracluster light (ICL), has the potential to be read as a memoir of mass accretion: informative, individual, and yet imperfect. Here, we combine dark-matter-only zoom-in simulations from the Symphony suite with the Nimbus “star-tagging” model of the stellar halo to assess how much information about the mass assembly of an individual galaxy cluster can be gleaned from idealized measurements of ICL outskirts. We show that the edges of a cluster’s stellar profile—the primary (R sp⋆,1 ) and secondary (R sp⋆,2 ) stellar “splashback” radii—are sensitive to both continuous mass accretion histories (MAHs) and discrete merger events, making them potentially powerful probes of a cluster’s past. We find that R sp⋆,1 strongly correlates with the cluster’s mass ∼1 dynamical time ago, while R sp⋆,2 traces more recent MAH to a slightly lesser degree. In combination, these features can further distinguish between clusters that have and have not undergone a major merger within the past dynamical time. We use both to predict realistic cluster MAHs with the MultiCAM framework. These outer ICL features are significantly more sensitive to mass accretion and merger histories than the stellar mass gap and halo concentration, and perform comparably to the commonly used X-ray-based tracer of relaxedness, x off . While our analysis is idealized, the relevant ICL features are potentially detectable in next-generation deep imaging of nearby clusters. This work highlights the promise of ICL measurements and lays the groundwork for more detailed forecasts of their power.

79 ASTRONOMY AND ASTROPHYSICS↗

Gas phase synthesis of the C40 nano bowl C 40 H 10

Nanobowls represent vital molecular building blocks of end-capped nanotubes and fullerenes detected in combustion systems and in deep space such as toward the planetary nebula TC-1, but their fundamental formation mechanisms have remained elusive. By merging molecular beam experiments with electronic structure calculations, we reveal a complex chain of reactions initiated through the gas-phase preparation of benzocorannulene (C 24 H 12 ) via ring annulation of the corannulenyl radical (C 20 H 9 • ) by vinylacetylene (C 4 H 4 ) as identified isomer-selectively in situ via photoionization efficiency curves and photoion mass-selected threshold photoelectron spectra. In silico studies provided compelling evidence that the benzannulation mechanism can be expanded to pentabenzocorannulene (C 40 H 20 ) followed by successive cyclodehydrogenation to the C40 nanobowl (C 40 H 10 ) – a fundamental building block of buckminsterfullerene (C 60 ). This high-temperature pathway opens up isomer-selective routes to nanobowls via resonantly stabilized free-radical intermediates and ring annulation in circumstellar envelopes of carbon stars and planetary nebulae as their descendants eventually altering our insights of the complex chemistry of carbon in our Galaxy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational general relativistic force-free electrodynamics: I. Multi-coordinate implementation and testing

General relativistic force-free electrodynamics is one possible plasma-limit employed to analyze energetic outflows in which strong magnetic fields are dominant over all inertial phenomena. The amazing images of black hole (BH) shadows from the Galactic Center and the M87 galaxy provide a first direct glimpse into the physics of accretion flows in the most extreme environments of the universe. The efficient extraction of energy in the form of collimated outflows or jets from a rotating BH is directly linked to the topology of the surrounding magnetic field. We aim at providing a tool to numerically model the dynamics of such fields in magnetospheres around compact objects, such as BHs and neutron stars. To do so, we probe their role in the formation of high energy phenomena such as magnetar flares and the highly variable teraelectronvolt emission of some active galactic nuclei. In this work, we present numerical strategies capable of modeling fully dynamical force-free magnetospheres of compact astrophysical objects. Here, we provide implementation details and extensive testing of our implementation of general relativistic force-free electrodynamics in Cartesian and spherical coordinates using the infrastructure of the E INSTEIN T OOLKIT . The employed hyperbolic/parabolic cleaning of numerical errors with full general relativistic compatibility allows for fast advection of numerical errors in dynamical spacetimes. Such fast advection of divergence errors significantly improves the stability of the general relativistic force-free electrodynamics modeling of BH magnetospheres.

79 ASTRONOMY AND ASTROPHYSICS↗

Signatures of Velocity-dependent Dark Matter Self-interactions in Milky Way-mass Halos

In this work, we explore the impact of elastic, anisotropic, velocity-dependent dark matter (DM) self-interactions on the host halo and subhalos of Milky Way (MW)--mass systems. We consider a generic self-interacting dark matter (SIDM) model parameterized by the masses of a light mediator and the DM particle. The ratio of these masses, w, sets the velocity scale above which momentum transfer due to DM self-interactions becomes inefficient. We perform high-resolution zoom-in simulations of an MW-mass halo for values of w that span scenarios in which self-interactions either between the host and its subhalos or only within subhalos efficiently transfer momentum, and we study the effects of self-interactions on the host halo and on the abundance, radial distribution, orbital dynamics, and density profiles of subhalos in each case. The abundance and properties of surviving subhalos are consistent with being determined primarily by subhalo--host halo interactions. In particular, subhalos on radial orbits in models with larger values of the cross section at the host halo velocity scale are more susceptible to tidal disruption owing to mass loss from ram pressure stripping caused by self-interactions with the host. This mechanism suppresses the abundance of surviving subhalos relative to collisionless DM simulations, with stronger suppression for larger values of w. Thus, probes of subhalo abundance around MW-mass hosts can be used to place upper limits on the self-interaction cross section at velocity scales of ~200 km s-1, and combining these measurements with the orbital properties and internal dynamics of subhalos may break degeneracies among velocity-dependent SIDM models.

79 ASTRONOMY AND ASTROPHYSICS↗

Hydrodynamical Simulations of Black Hole Binary Formation in AGN Disks

Abstract We study close encounters between two single black holes (BHs) embedded in an AGN disk using a series of global 2D hydrodynamics simulations. We find that when the disk density is sufficiently high, bound BH binaries can be formed by the collision of their circum-single disks. Our analysis demonstrates that, after a BH pair passes the pericenter of their relative trajectory, post-collision gas drag may slow down the BHs, possibly forcing the two BHs to stay tightly bound. A binary formed by a close encounter can have a compact semimajor axis, large eccentricity, and retrograde orbital angular momentum. We provide a fitting formula that can accurately predict whether a close encounter can form a binary based on the gas mass and the incoming energy of the encounter. This fitting formula can be easily implemented in other long-term simulations that study the dynamical evolution of BHs in active galactic nucleus disks.

79 ASTRONOMY AND ASTROPHYSICS↗

Star-forming S0 galaxies in the SDSS-IV MaNGA survey

ABSTRACT To investigate star-forming activities in early-type galaxies, we select a sample of 52 star-forming S0 galaxies (SFS0s) from the SDSS-IV MaNGA survey. We find that SFS0s have smaller stellar masses compared to normal S0s in MaNGA. After matching the stellar mass to select the control sample, we find that the mean Sérsic index of SFS0 bulges (1.76 ± 0.21) is significantly smaller than that of the control sample (2.57 ± 0.20), suggesting the existence of pseudo-bulges in SFS0s. After introducing the environmental information, SFS0s show smaller spin parameters in the field than in groups, while the control sample has no obvious difference in different environments, which may suggest different dynamical processes in SFS0s. Furthermore, with derived N/O and O/H abundance ratios, SFS0s in the field show nitrogen enrichment, providing evidence for the accretion of metal-poor gas in the field environment. To study the star-formation relation, we show that the slope of the spatially resolved star-formation main sequence is nearly 1.0 with MaNGA integral field unit data, confirming the self-regulation of star-formation activity on kpc scales.

79 ASTRONOMY AND ASTROPHYSICS↗

Core formation in high-z massive haloes: heating by post-compaction satellites and response to AGN outflows

ABSTRACT Observed rotation curves in star-forming galaxies indicate a puzzling dearth of dark matter in extended flat cores within haloes of mass ≥1012M⊙ at z ∼ 2. This is not reproduced by current cosmological simulations, and supernova-driven outflows are not effective in such massive haloes. We address a hybrid scenario where post-compaction merging satellites heat up the dark-matter cusps by dynamical friction, allowing active galactic nucleus (AGN)-driven outflows to generate cores. Using analytic and semi-analytic models (SatGen), we estimate the dynamical friction heating as a function of satellite compactness for a cosmological sequence of mergers. Cosmological simulations (VELA) demonstrate that satellites of initial virial masses >1011.3M⊙, which undergo wet compactions, become sufficiently compact for significant heating. Constituting a major fraction of the accretion on to haloes ≥1012M⊙, these satellites heat up the cusps in half a virial time at z ∼ 2. Using a model for outflow-driven core formation (CuspCore), we demonstrate that the heated dark-matter cusps develop extended cores in response to removal of half the gas mass, while the more compact stellar systems remain intact. The mergers keep the dark matter hot, while the gas supply, fresh and recycled, is sufficient for the AGN outflows. AGNs indeed become effective in haloes ≥1012M⊙, where the black hole growth is no longer suppressed by supernovae and its compaction-driven rapid growth is maintained by a hot circumgalactic medium. For simulations to reproduce the dynamical friction effects, they should resolve the compaction of the massive satellites and avoid artificial tidal disruption. AGN feedback could be boosted by clumpy black hole accretion and clumpy response to AGN.

Dekel, Avishai↗

Jets with a Twist: The Emergence of FR0 Jets in a 3D GRMHD Simulation of Zero-angular-momentum Black Hole Accretion

Abstract Spinning supermassive black holes (BHs) in active galactic nuclei magnetically launch relativistic collimated outflows, or jets. Without angular momentum supply, such jets are thought to perish within 3 orders of magnitude in distance from the BH, well before reaching kiloparsec scales. We study the survival of such jets at the largest scale separation to date, via 3D general relativistic magnetohydrodynamic simulations of rapidly spinning BHs immersed into uniform zero-angular-momentum gas threaded by a weak vertical magnetic field. We place the gas outside the BH sphere of influence, or the Bondi radius, chosen to be much larger than the BH gravitational radius, R B = 10 3 R g . The BH develops dynamically important large-scale magnetic fields, forms a magnetically arrested disk (MAD), and launches relativistic jets that propagate well outside R B and suppress BH accretion to 1.5% of the Bondi rate, M ̇ B . Thus, low-angular-momentum accretion in the MAD state can form large-scale jets in Fanaroff–Riley (FR) type I and II galaxies. Subsequently, the disk shrinks and exits the MAD state: barely a disk (BAD), it rapidly precesses, whips the jets around, globally destroys them, and lets 5%–10% of M ̇ B reach the BH. Thereafter, the disk starts rocking back and forth by angles 90°–180°: the rocking accretion disk (RAD) launches weak intermittent jets that spread their energy over a large area and suppress BH accretion to ≲2% M ̇ B . Because the BAD and RAD states tangle up the jets and destroy them well inside R B , they are promising candidates for the more abundant, but less luminous, class of FR0 galaxies.

Astronomy & Astrophysics↗

Science Validation for Dark Energy Research with Optical Imaging Surveys

The Universe has been expanding at an accelerating rate over the past several billion years, as though an unknown form of dark energy permeates all of space. The ultimate scientific goal of the proposed research is to distinguish between different physical mechanisms that could account for this observed accelerated expansion, for example, the zero-point energy of the vacuum, a dynamical form of energy that varies in time and/or space, or a modification to our theory of gravity. Wide- field optical imaging surveys of the night sky can test these competing models by measuring both the cosmic expansion history and the growth of large-scale structure. To this end, the Dark Energy Survey (DES) has cataloged several hundred million galaxies and thousands of supernovae. The Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST), will enlarge the census to billions of galaxies and hundreds of thousands of supernovae. A critical question for these dark energy experiments is whether systematic uncertainties can continue to be controlled at a level to keep pace with the statistical precision offered by such enormous datasets. The immediate research objectives of this project were (1) to prepare and validate input datasets that are the foundation of cosmological analyses with DES, and (2) to prepare for value-added characterization of Rubin Observatory commissioning data to inform early operations and accelerate the realization of dark energy science from LSST data products. For DES, we assembled and curated cosmology-ready data releases that include value-added components such as enhanced photometric and astrometric calibrations, alternative source extraction algorithms, maps of the survey coverage and survey conditions, object classifications, object quality selections, galaxy shapes, and photometric redshifts. We used the galaxy clustering technique to validate the photometric redshift distributions of various galaxy samples to be used as lenses in combined studies of galaxy clustering and weak gravitational lensing. The galaxy clustering redshift analysis was enhanced by use of a larger sample of reference galaxies from the eBOSS spectroscopic survey that extends to higher redshifts. For LSST, we prepared for science validation studies of commissioning data aimed at dark energy science capability that extend beyond the normative system-level tests to be done by the Rubin Observatory Construction Project. We identified a set of proposed survey strategies and candidate target fields that could be observed during the commissioning period to enhance science validation activities related to studies of dark energy.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cosmology at high redshift — a probe of fundamental physics

An observational program focused on the high redshift (2<z<6) Universe has the opportunity to dramatically improve over upcoming LSS and CMB surveys on measurements of both the standard cosmological model and its extensions. Using a Fisher matrix formalism that builds upon recent advances in Lagrangian perturbation theory, we forecast constraints for future spectroscopic and 21-cm surveys on the standard cosmological model, curvature, neutrino mass, relativistic species, primordial features, primordial non-Gaussianity, dynamical dark energy, and gravitational slip. We compare these constraints with those achievable by current or near-future surveys such as DESI and Euclid, all under the same forecasting formalism, and compare our formalism with traditional linear methods. Here, our Python code FishLSS — used to calculate the Fisher information of the full shape power spectrum, CMB lensing, the cross-correlation of CMB lensing with galaxies, and combinations thereof — is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗

How accurately can we detect the splashback radius of dark matter haloes and its correlation with accretion rate?

ABSTRACT The splashback radius (Rsp) of dark matter haloes has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure Rsp and to assess the significance of its detection. In this paper, we use simulations to study the precision and accuracy to which we can detect the splashback radius with 3D density, 3D subhalo, and weak lensing profiles. We study how well various methods and tracers recover Rsp by comparing it with the value measured directly from particle dynamics. We show that estimates of Rsp from density and subhalo profiles correspond to different percentiles of the underlying Rsp distribution of particle orbits. At low accretion rates, a second caustic appears and can bias results. Finally, we show that upcoming lensing surveys may be able to constrain the Rsp–Γdyn relation directly.

79 ASTRONOMY AND ASTROPHYSICS↗

Diverse Chemo-Dynamical Properties of Nitrogen-Rich Stars Identified from Low-Resolution Spectra

The second generation of stars in the globular clusters (GCs) of the Milky Way (MW) exhibit unusually high N, Na, or Al, compared to typical Galactic halo stars at similar metallicities. The halo field stars enhanced with such elements are believed to have originated in disrupted GCs or escaped from existing GCs. We identify such stars in the metallicity range –3.0 < [Fe/H] < 0.0 from a sample of ~ 36,800 giant stars observed in the Sloan Digital Sky Survey and Large Sky Area Multi-Object Fiber Spectroscopic Telescope survey, and present their dynamical properties. The N-rich population (NRP) and N-normal population (NNP) among our giant sample do not exhibit similarities in either in their metallicity distribution function (MDF) or dynamical properties. We find that, even though the MDF of the NRP looks similar to that of the MW’s GCs in the range of [Fe/H] < –1.0, our analysis of the dynamical properties does not indicate similarities between them in the same metallicity range, implying that the escaped members from existing GCs may account for a small fraction of our N-rich stars, or the orbits of the present GCs have been altered by the dynamical friction of the MW. We also find a significant increase in the fraction of N-rich stars in the halo field in the very metal-poor (VMP; [Fe/H] < –2.0) regime, comprising up to ~ 20% of the fraction of the N-rich stars below [Fe/H] = –2.5, hinting that partially or fully destroyed VMP GCs may have in some degree contributed to the Galactic halo. A more detailed dynamical analysis of the NRP reveals that our sample of N-rich stars do not share a single common origin. Although a substantial fraction of the N-rich stars seem to originate from the GCs formed in situ, more than 60% of them are not associated with those of typical Galactic populations, but probably have extragalactic origins associated with Gaia Sausage/Enceladus, Sequoia, and Sagittarius dwarf galaxies, as well as with presently unrecognized progenitors.

79 ASTRONOMY AND ASTROPHYSICS↗

CODEX: Role of velocity substructure in the scaling relations of galaxy clusters

The use of galaxy clusters as cosmological probes relies on a detailed understanding of their properties. They define cluster selection and ranking linked to a cosmologically significant cluster mass function. Previous studies have employed small samples of clusters, concentrating on achieving the first calibrations of cluster properties with mass, while the diversity of cluster properties has been revealed via detailed studies. The large spectroscopic follow-up on the CODEX cluster sample with SDSS and NOT enables a detailed study of hundreds of clusters, lifting the limitations of previous samples. We aim to update the spectroscopic cluster identification of CODEX by running the spectroscopic group finder on the follow-up spectroscopy results and connecting the dynamical state of clusters to their scaling relations. We implemented a reproducible spectroscopic membership determination and cleaning procedures, based on the redMaPPer membership, running the spectroscopic group finder on the follow-up spectroscopy results and cleaning the membership for spectroscopic outliers. We applied the Anderson-Darling test for velocity substructure and analysed its influence on the scaling relations. We also tested the effect of the X-ray-to-optical centre offset on the scaling relations. We report on the scaling relations between richness, X-ray luminosity, and velocity dispersion for a complete sample of clusters with at least 15 members. Clusters with velocity substructure exhibit enhanced velocity dispersion for a given richness and are characterized by 2.5 times larger scatter. Clusters that have a strong offset in X-ray-to-optical centres have comparable scaling relations as clusters with substructure. We demonstrate that there is a consistency in the parameters of the scaling relations for the low- and high-richness galaxy clusters. Splitting the clusters by redshift, we note a decrease in scatter with redshift in all scaling relations. We localize the redshift range where a high scatter is observed to z < 0.15, which is in agreement with the literature results on the scatter. We note that the increase in scatter for both high- and low-luminosity clusters is z < 0.15, suggesting that both cooling and the resulting active galactic nucleus feedback are at the root of this scatter.

79 ASTRONOMY AND ASTROPHYSICS↗

Low-temperature formation of pyridine and (iso)quinoline via neutral–neutral reactions

Aromatic molecules represent fundamental building blocks in prebiotic chemistry and are contemplated as vital precursors to DNA and RNA nitrogen bases. However, despite the identification of some 300 molecules in extraterrestrial environments, the pathways to pyridine (C 5 H 5 N), pyridinyl (C 5 H 4 N•), and (iso)quinoline (C 9 H 7 N) – the simplest represen-tative of mono and bicyclic aromatic molecule carrying nitrogen – are elusive. Here, we afford compelling evidence on the gas-phase formation of methylene amidogen (H 2 CN•) and cyanomethyl (H 2 CCN•) radicals via molecular beam studies and electronic structure calcula-tions. The modeling of the chemistries of Taurus Molecular Cloud (TMC-1) and Titan’s atmosphere contemplates a complex chain of reactions synthesizing pyridine, pyridinyl, and (iso)quinoline from H 2 CN• and H 2 CCN• at levels of up to 75%. Furthermore, this study affords unique entry points to precursors of DNA and RNA nitrogen bases in hydrocarbon-rich extraterrestrial environments thus changing the way we think about the origin of prebiotic molecules in our Galaxy.

79 ASTRONOMY AND ASTROPHYSICS↗

Delving into the depths of NGC 3783 with XRISM V. Broad-band modelling of ionized outflows

The Seyfert 1 galaxy NGC 3783 hosts a multi-phase warm absorber (WA) that has been extensively studied in the X-ray band. High-resolution spectra from 2000−2001 revealed a complex outflow with multiple ionization and velocity components. Two decades later, new XMM-Newton and XRISM observations allow us to investigate the long-term evolution of these outflows. We performed joint spectral modelling of the XMM-Newton /RGS and XRISM/Resolve time-averaged spectra using the pion photoionization code within SPEX. We derived the ionization parameter, column density, turbulent velocity, and outflow velocity for each absorption component, and investigated their thermal stability and absorption measure distribution (AMD) to characterize the physical and dynamical properties of the WA in NGC 3783 in 2024. We compared these results with the 2000−2001 epoch to assess long-term variability, stability, and possible changes in the absorber population. We identify eight WA components spanning log ξ = 1.08 − 3.38 and outflow velocities of 480−1230 km s −1 . The ranges of column densities and turbulent velocities remain broadly consistent with the WAs from 2000−2001, but the earlier data contained more low-ionization, high-velocity components. The total column density of all the outflows in 2024 is 1.5 times larger than in 2000−2001, which means that it has been replenished by fresh material. The dominant unresolved transition array (UTA) absorber (component B3) has increased its column density by a factor of three while maintaining a similar ionization parameter. The WAs in NGC 3783 have undergone significant structural and dynamical evolution over the past 24 years.

79 ASTRONOMY AND ASTROPHYSICS↗

Impact of Varying Dark Energy on Future Large-Scale Structure Studies

The cosmos withholds multiple mysteries such as dark forms of matter and energy that are yet beyond human comprehension. In the standard cosmological model, $\Lambda$CDM, the cosmological constant $\Lambda$ is thought to be responsible for the late accelerated expansion of the Universe. However, recent results from the Dark Energy Spectroscopic Instrument (DESI) suggest the possibility of evolving dark energy, which warrants further exploration. Future surveys such as the Rubin Observatory Legacy Survey of Space and Time (LSST) will map the large-scale structure (LSS) with unprecedented precision, giving us valuable statistical information about the cosmos. The goal of this research is to investigate the potential impact of an evolving dark energy scenario on cosmological parameters constrained by LSS probes, as will be mapped by the LSST. For this, we examine the power spectra of lens galaxies, source galaxies, and cross-power spectra between lens galaxies and source galaxies. The combination of these statistics is commonly referred to as '3 $\times$ 2 points'. For this investigation, we created a set of simulations resembling LSST data and used them to perform cosmological parameter inference in two scenarios: one in which the simulated data is based on the fiducial model and another on evolving dark energy. We then examined the degeneracy between the cosmological parameters and checked for potential shifts in the parametric space when the data contains dynamical dark energy but the modeling assumes $\Lambda$CDM. Our findings indicate that mismodeling the dark energy equation of state can significantly impact parameter inference, particularly affecting the total matter density, $\Omega_m$, and the growth of structures, as represented by the $S_8$ parameter. These results highlight the importance of further exploring extensions of the $\Lambda$CDM model in future LSS studies.

Yaman Acharya, A.↗

Reconstructing the Stripping History of the Sagittarius Stream with Neural Networks

The Sagittarius (Sgr) Stream is produced by the ongoing disruption of the Sgr dwarf spheroidal (dSph) galaxy and is thought to contain multiple wraps that were stripped during different pericentric passages. In this study, we introduce a neural-network–based method trained on N-body simulations to infer the stripping time of Sgr Stream stars directly from their phase-space coordinates. We combine spectroscopic data from SEGUE, APOGEE DR17, and LAMOST DR7 low-resolution spectroscopic (LRS) survey with Gaia EDR3 astrometry and distance estimates from the latest StarHorse catalog to identify high-quality Sgr Stream members. Applying our method to these stars, we measure a clear metallicity gradient with stripping time, well described by a linear relation with slope ∼0.3 dex Gyr −1 . We further predict the stripping times of globular clusters previously suggested to originate from the Sgr dSph. M 54, Terzan 7, Terzan 8, and Arp 2 exhibit stripping times consistent with being currently bound to the Sgr remnant. Pal 12, Whiting 1, and NGC 2419 are inferred to have been stripped 0.9 ± 0.1, 1.1 ± 0.2, and 2.1 ± 0.2 Gyr ago, respectively. For NGC 4147 and NGC 5634, whose membership in the Sgr system remains uncertain, our analysis suggests stripping times of 1.1 ± 0.4 and 1.1 ± 0.1 Gyr, respectively, if they are ultimately confirmed as genuine Sgr members. These results demonstrate that data-driven models of dynamical stripping histories offer a promising approach for reconstructing the formation and chemical evolution of the Sgr Stream.

79 ASTRONOMY AND ASTROPHYSICS↗

Halo formation from Yukawa forces in the very early Universe

If long-range attractive forces exist and are stronger than gravity then cosmic halo formation can begin in the radiation-dominated era. We study a simple realization of this effect in a system where dark matter fermions have Yukawa interactions mediated by scalar particles, analogous to the Higgs boson in the standard model. We develop a self-consistent description of the system including exact background dynamics of the scalar field, and precise modeling of the fermion density fluctuations. For the latter, we provide accurate approximations for the linear growth as well as quantitative modeling of the nonlinear evolution using N-body simulations. We find that halo formation occurs exponentially fast and on scales substantially larger than simple estimates predict. Here, the final fate of these halos remains uncertain, but could be annihilation, dark stars, primordial black holes, or even the existence of galaxy-sized halos at matter-radiation equality. More generally, our results demonstrate the importance of mapping scalar-mediated interactions onto structure formation outcomes and constraints for beyond the standard model theories.

79 ASTRONOMY AND ASTROPHYSICS↗