Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “haloes”

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 433 records · Page 24

Measuring M31 globular cluster ages and metallicities using both photometry and spectroscopy

ABSTRACT The ages and metallicities of globular clusters play an important role not just in testing models for their formation and evolution but also in understanding the assembly history of their host galaxies. Here, we use a combination of imaging and spectroscopy to measure the ages and metallicities of globular clusters in M31, the closest massive galaxy to our own. We use the strength of the near-infrared calcium triplet spectral feature to provide a relatively age-insensitive prior on the metallicity when fitting stellar population models to the observed photometry. While the age–extinction degeneracy is an issue for globular clusters projected on to the disc of M31, we find generally old ages for globular clusters in the halo of M31 and in its satellite galaxy NGC 205 in line with previous studies. We measure ages for a number of outer halo globular clusters for the first time, finding that globular clusters associated with halo substructure extend to younger ages and higher metallicities than those associated with the smooth halo. This is in line with the expectation that the smooth halo was accreted earlier than the substructured halo.

Usher, Christopher (ORCID:0000000273837106)↗

Cosmic Frontier Project at SMU

The SMU group led by Professor Robert Kehoe continued to work to constrain dark energy measurements of the Dark Energy Spectroscopic Instrument (DESI). These constraints are obtained by measuring the expected peak in characteristic distance scale between galaxies, the ‘BAO peak’, as a function of redshift for four different types of galaxy. DESI will cover the redshift range to z<3.5 with approximately 40 million galaxy spectra. SMU’s primary emphasis during the funding period was on the early DESI data, including the development of software that would allow the measurement of the impact of incomplete targeting of galaxies in the survey. We provided software to calculate the weights that are critical to permit us to correctly compare data to the theoretical calculations of dark matter halos. SMU also began development of software that can be integrated into the analysis pipeline that uses alternate targeting to generate these weights. Baryonic physics has a critical impact on galaxy development and placement within dark matter halos. SMU also increased its efforts in the area of modeling this galaxy-halo connection in two areas. The primary emphasis was to work on a sub-halo abundance matching scheme using dark matter simulation with very fine mass resolution (‘Uchuu’). Our effort here is to produce models for use by DESI that more correctly account for the subhalo structure of dark matter distribution, and its manifestation in satellite galaxies. Our initial effort had been for emission line galaxies. In the funding period, we were able to produce first halo-occupation distribution models. This work is slated to contribute to DESI Key Project papers. We also furthered our ongoing effort to study galaxy placement in halos using semi-analytic methods. Kehoe’s group also continued to maintain its strong role in DESI Operations, with SMU postdoc serving as a Lead Observing Scientist.

79 ASTRONOMY AND ASTROPHYSICS↗

New limits from microlensing on Galactic black holes in the mass range 10 M ⊙ < M < 1000 M ⊙

We searched for long-duration microlensing events originating from intermediate-mass black holes (BH) in the halo of the Milky Way, using archival data from the EROS-2 and MACHO photometric surveys towards the Large Magellanic Cloud (LMC). We combined data from these two surveys to create a common database of light curves for 14.1 million objects in the LMC, covering a total duration of 10.6 years, with flux series measured in four wide passbands. We carried out a microlensing search on these light curves, complemented by the light curves of 22.7 million objects, observed only by EROS-2 or only by MACHO, over about 7 years, with flux series measured in only two passbands. A likelihood analysis, taking into account the LMC self-lensing and Milky Way disk contributions, allows us to conclude that compact objects with masses in the range 10 – 100 M ⊙ cannot make up more than ~15% of a standard halo total mass (at a 95% confidence level). Our analysis sensitivity weakens for heavier objects, although we still rule out the possibility of ~50% of the halo being made of ~1000 M ⊙ BHs. Combined with previous EROS results, an upper limit of ~15% of the total halo mass can be obtained for the contribution of compact halo objects in the mass range 10 –6 – 10 2 M ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗

Dissipative dark matter on FIRE – I. Structural and kinematic properties of dwarf galaxies

We present the first set of cosmological baryonic zoom-in simulations of galaxies including dissipative self-interacting dark matter (dSIDM). These simulations utilize the Feedback In Realistic Environments galaxy formation physics, but allow the dark matter to have dissipative self-interactions analogous to standard model forces, parametrized by the self-interaction cross-section per unit mass, (σ/m), and the dimensionless degree of dissipation, 0 < f diss < 1. We survey this parameter space, including constant and velocity-dependent cross-sections, and focus on structural and kinematic properties of dwarf galaxies with $M_{\rm halo} \sim 10^{10-11}{\, \rm M_\odot }$ and $M_{\ast } \sim 10^{5-8}{\, \rm M_\odot }$. Central density profiles (parametrized as ρ ∝ rα) of simulated dwarfs become cuspy when $(\sigma /m)_{\rm eff} \gtrsim 0.1\, {\rm cm^{2}\, g^{-1}}$ (and f diss = 0.5 as fiducial). The power-law slopes asymptote to α ≈ -1.5 in low-mass dwarfs independent of cross-section, which arises from a dark matter ‘cooling flow’. Through comparisons with dark matter only simulations, we find the profile in this regime is insensitive to the inclusion of baryons. However, when $(\sigma /m)_{\rm eff} \ll 0.1\, {\rm cm^{2}\, g^{-1}}$, baryonic effects can produce cored density profiles comparable to non-dissipative cold dark matter (CDM) runs but at smaller radii. Simulated galaxies with $(\sigma /m) \gtrsim 10\, {\rm cm^{2}\, g^{-1}}$ and the fiducial fdiss develop significant coherent rotation of dark matter, accompanied by halo deformation, but this is unlike the well-defined thin ‘dark discs’ often attributed to baryon-like dSIDM. The density profiles in this high cross-section model exhibit lower normalizations given the onset of halo deformation. For our surveyed dSIDM parameters, halo masses and galaxy stellar masses do not show appreciable difference from CDM, but dark matter kinematics and halo concentrations/shapes can differ.

79 ASTRONOMY AND ASTROPHYSICS↗

The mass profiles of dwarf galaxies from Dark Energy Survey lensing

We present a novel approach to extracting dwarf galaxies from photometric data to measure their average halo mass profile with weak lensing. We characterize their stellar mass and redshift distributions with a spectroscopic calibration sample. By combining the ${\sim} 5000\,\mathrm{deg}^2$ multiband photometry from the Dark Energy Survey and redshifts from the Satellites Around Galactic Analogs Survey with an unsupervised machine learning method, we select a low-mass galaxy sample spanning redshifts $z\lt 0.3$ and divide it into three mass bins. From low to high median mass, the bins contain [146 420, 330 146, 275 028] galaxies and have median stellar masses of $\log _{10}(M_*/\text{M}_\odot)=\left[8.52\substack{+0.57 -0.76},\, 9.02\substack{+0.50 -0.64},\, 9.49\substack{+0.50 -0.58}\right]$ . We measure the stacked excess surface mass density profiles, $\Delta \Sigma (R)$, of these galaxies using galaxy–galaxy lensing with a signal-to-noise ratio of [14, 23, 28]. Through a simulation-based forward-modelling approach, we fit the measurements to constrain the stellar-to-halo mass relation and find the median halo mass of these samples to be $\log _{10}(M_{\rm halo}/\text{M}_\odot)$ = [$10.67\substack{+0.2 -0.4}$, $11.01\substack{+0.14 -0.27}$, $11.40\substack{+0.08 -0.15}$]. The cold dark matter profiles are consistent with NFW (Navarro, Frenk, and White) profiles over scales ${\lesssim} 0.15 \, {h}^{-1}$ Mpc. We find that ${\sim} 20$ per cent of the dwarf galaxy sample are satellites. This is the first measurement of the halo profiles and masses of such a comprehensive, low-mass galaxy sample. The techniques presented here pave the way for extracting and analysing even lower mass dwarf galaxies and for more finely splitting galaxies by their properties with future photometric and spectroscopic survey data.

dark matter↗

Modest set of simulations of local-type primordial non-Gaussianity in the DESI era

A measurement of a primordial non-Gaussianity (PNG) signal through late- or early-Universe probes has the potential to transform our understanding of the physics of the primordial Universe. While large-scale structure observables in principle contain vital information, interpreting these measurements is challenging due to poorly understood astrophysical effects. Luckily, N-body simulations, such as the abacuspng set presented in this study, consisting of 9 boxes, each of size L box = 2 Gpc/h and particle mass of 1.01 × 10 10 M ⊙ /h, provide a viable path forward. As validation, we find good agreement between the simulations and our expectations from one-loop perturbation theory (PT) and the "separate universe" method for the matter bispectrum, matter power spectrum and the halo bias parameter associated with PNG, b φ . As a science application, here we investigate the link between halo assembly bias and b φ for halo properties known to play a vital role in accurately predicting galaxy clustering: concentration, shear (environment), and accretion rate. We find a strong response for all three parameters, suggesting that the connection between b φ and the assembly history of halos needs to be taken into account by future PNG analyses. We further perform the first study of the b φ parameter from fits to early DESI data of the luminous red galaxy (LRG) and quasistellar object (QSO) samples and comment on the effect on f NL constraints for the allowed galaxy-halo models (note that σ[f NL ]∝ $\frac{σ[b_φ]}{b_φ}$). We find that the error on f NL is 15, 8, 7 for the LRGs at z = 0.5 and z = 0.8 and QSOs at z = 1.4, respectively, suggesting that a thorough understanding of galaxy assembly bias is warranted so as to perform robust high-precision analysis of local-type PNG with future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Gravothermal collapse and the diversity of galactic rotation curves

The rotation curves of spiral galaxies exhibit a great diversity that challenges our understanding of galaxy formation and the nature of dark matter. Previous studies showed that in self-interacting dark matter (SIDM) models with a cross section per unit mass of σ/m ≈ $\mathscr{O}$(1) cm 2 /g, the predicted dark matter central densities are a good match to the observed densities in galaxies. Here, in this work, we explore a regime with a larger cross section of σ/m ≈ 20−40 cm 2 /g in dwarf galactic halos. We will show that such strong dark matter self-interactions can further amplify the diversity of halo densities inherited from their assembly history. High concentration halos can enter the gravothermal collapse phase within 10 Gyr, resulting in a high density, while low concentration ones remain in the expansion phase and have a low density. We fit the rotation curves of 14 representative low surface brightness galaxies and demonstrate how the large range of observed central densities is naturally accommodated in the strong SIDM regime of σ/m ≈ 20 – 40 cm 2 /g. Galaxies that are outliers in the previous studies, due to their high halo central densities, are no longer outliers in this SIDM regime as their halos would be in the collapse phase. For galaxies with a low density, the SIDM fits are robust to the variation of the cross section. Our findings open up a new window for testing gravothermal collapse, the unique signature of strong dark matter self-interactions, and exploring a broader SIDM model space. As an example, we illustrate how the larger cross sections favored by our fits, together with upper limits from strong lensing observations in clusters, pick out the preferred SIDM model space for a dark matter particle coupled to a light gauge boson in the Born regime.

dark matter↗

Distant Echoes of the Milky Way’s Last Major Merger

The majority of the Milky Way’s stellar halo consists of debris from our galaxy’s last major merger, the Gaia-Sausage-Enceladus (GSE). In the past few years, stars from the GSE have been kinematically and chemically studied in the inner 30 kpc of our galaxy. However, simulations predict that accreted debris could lie at greater distances, forming substructures in the outer halo. Here we derive metallicities and distances using Gaia DR3 XP spectra for an all-sky sample of luminous red giant stars, and map the outer halo with kinematics and metallicities out to 100 kpc. We obtain follow-up spectra of stars in two strong overdensities—including the previously identified outer Virgo Overdensity—and find them to be relatively metal rich and on predominantly retrograde orbits, matching predictions from simulations of the GSE merger. We argue that these are apocentric shells of GSE debris, forming 60–90 kpc counterparts to the 15–20 kpc shells that are known to dominate the inner stellar halo. Extending our search across the sky with literature radial velocities, we find evidence for a coherent stream of retrograde stars encircling the Milky Way from 50 to 100 kpc, in the same plane as the Sagittarius Stream but moving in the opposite direction. These are the first discoveries of distant and structured imprints from the GSE merger, cementing the picture of an inclined and retrograde collision that built up our galaxy’s stellar halo.

79 ASTRONOMY AND ASTROPHYSICS↗

The DESI One-percent Survey: Evidence for Assembly Bias from Low-redshift Counts-in-cylinders Measurements

Abstract We explore the galaxy-halo connection information that is available in low-redshift samples from the early data release of the Dark Energy Spectroscopic Instrument (DESI). We model the halo occupation distribution (HOD) from z = 0.1 to 0.3 using Survey Validation 3 (SV3; a.k.a., the One-Percent Survey) data of the DESI Bright Galaxy Survey. In addition to more commonly used metrics, we incorporate counts-in-cylinders (CiC) measurements, which drastically tighten HOD constraints. Our analysis is aided by the Python package, galtab , which enables the rapid, precise prediction of CiC for any HOD model available in halotools . This methodology allows our Markov chains to converge with much fewer trial points, and enables even more drastic speedups due to its GPU portability. Our HOD fits constrain characteristic halo masses tightly and provide statistical evidence for assembly bias, especially at lower luminosity thresholds: the HOD of central galaxies in z ∼ 0.15 samples with limiting absolute magnitude M r < −20.0 and M r < −20.5 samples is positively correlated with halo concentration with a significance of 99.9% and 99.5%, respectively. Our models also favor positive central assembly bias for the brighter M r < −21.0 sample at z ∼ 0.25 (94.8% significance), but there is no significant evidence for assembly bias with the same luminosity threshold at z ∼ 0.15. We provide our constraints for each threshold sample’s characteristic halo masses, assembly bias, and other HOD parameters. These constraints are expected to be significantly tightened with future DESI data, which will span an area 100 times larger than that of SV3.

79 ASTRONOMY AND ASTROPHYSICS↗

Milky Way-est: Cosmological Zoom-in Simulations with Large Magellanic Cloud and Gaia–Sausage–Enceladus Analogs

Abstract We present Milky Way-est, a suite of 20 cosmological cold-dark-matter-only zoom-in simulations of Milky Way (MW)-like host halos. Milky Way-est hosts are selected such that they (i) are consistent with the MW’s measured halo mass and concentration, (ii) accrete a Large Magellanic Cloud (LMC)-like (≈10 11 M ⊙ ) subhalo within the last 2 Gyr on a realistic orbit, placing them near 50 kpc from the host center at z ≈ 0, and (iii) undergo a >1:5 sub-to-host halo mass ratio merger with a Gaia–Sausage–Enceladus (GSE)-like system at early times (0.67 < z < 3). Hosts satisfying these LMC and GSE constraints constitute <1% of all halos in the MW’s mass range, and their total masses grow rapidly at late times due to LMC analog accretion. Compared to hosts of a similar final halo mass that are not selected to include LMC and GSE analogs, Milky Way-est hosts contain 22% more subhalos with present-day virial masses above 10 8 M ⊙ throughout the virial radius, on average. This enhancement reaches ≈80% in the inner 100 kpc and is largely, if not entirely, due to LMC-associated subhalos. These systems also induce spatial anisotropy in Milky Way-est subhalo populations, with ≈60% of the total subhalo population within 100 kpc found in the current direction of the LMC. Meanwhile, we find that GSE-associated subhalos do not significantly contribute to present-day Milky Way-est subhalo populations. These results provide context for our Galaxy’s dark matter structure and subhalo population and will help interpret a range of measurements that are currently only possible in the MW.

79 ASTRONOMY AND ASTROPHYSICS↗

Subhalos in Galaxy Clusters: Coherent Accretion and Internal Orbits

Subhalo dynamics in galaxy cluster host halos govern the observed distribution and properties of cluster member galaxies. We use the IllustrisTNG simulation to investigate the accretion and orbits of subhalos found in cluster-size halos. We find that the median change in the major axis direction of cluster-size host halos is approximately 80° between a ∼ 0.1 and the present day. We identify coherent regions in the angular distribution of subhalo accretion, and ∼68% of accreted subhalos enter their host halo through ∼38% of the surface area at the virial radius. The majority of galaxy clusters in the sample have ∼2 such coherent regions. We further measure angular orbits of subhalos with respect to the host major axis and use a clustering algorithm to identify distinct orbit modes with varying oscillation timescales. The orbit modes correlate with subhalo accretion conditions. Subhalos in orbit modes with shorter oscillations tend to have lower peak masses and accretion directions somewhat more aligned with the major axis. One orbit mode, exhibiting the least oscillatory behavior, largely consists of subhalos that accrete near the plane perpendicular to the host halo major axis. Our findings are consistent with expectations from inflow from major filament structures and internal dynamical friction: most subhalos accrete through coherent regions, and more massive subhalos experience fewer orbits after accretion. Our work offers a unique quantification of subhalo dynamics that can be connected to how the intracluster medium strips and quenches cluster galaxies.

galaxy clusters↗

COZMIC. III. Cosmological Zoom-in Simulations of Self-interacting Dark Matter with Suppressed Initial Conditions

We present eight cosmological dark matter (DM)-only zoom-in simulations of a Milky Way–like system that include suppression of the linear matter power spectrum P(k), and/or velocity-dependent DM self-interactions, as the third installment of the COZMIC suite. We consider a model featuring a massive dark photon that mediates DM self-interactions and decays into massless dark fermions. The dark photon and dark fermions suppress linear matter perturbations, resulting in dark acoustic oscillations in P(k), which ultimately affect dwarf galaxy scales. The model also features a velocity-dependent elastic self-interaction between DM particles (SIDM), with a cross section that can alleviate small-scale structure anomalies. For the first time, our simulations test the impact of P(k) suppression on gravothermal evolution in an SIDM scenario that leads to core collapse in (sub)halos with present-day virial masses below ≈10 9 M ⊙ . In simulations with P(k) suppression and self-interactions, the lack of low-mass (sub)halos and the delayed growth of structure reduce the fraction of core-collapsed systems relative to SIDM simulations without P(k) suppression. In particular, P(k) suppression that saturates current warm DM constraints almost entirely erases core collapse in isolated halos. Models with less extreme P(k) suppression produce core collapse in ≈20% of subhalos and ≈5% of isolated halos above 10 8 M ⊙ , and also increase the abundance of extremely low-concentration isolated low-mass halos relative to SIDM. These results reveal a complex interplay between early and late-Universe DM physics, revealing new discovery scenarios in the context of upcoming small-scale structure measurements.

dark matter↗

Flux-driven algebraic damping of m = 2 diocotron mode

Experiments with pure electron plasmas in a Malmberg–Penning trap have observed linear in time, algebraic damping of m = 2 diocotron modes. Transport due to small field asymmetries produces a low-density halo of electrons moving radially outward from the plasma core, and the mode damping begins when the halo reaches the resonant radius of the mode. The damping rate is proportional to the flux of halo particles through the resonant layer. The damping is related to, but distinct from spatial Landau damping in which a linear wave–particle resonance produces exponential damping. This paper reports an analytic theory that captures the main signatures reported for this novel damping, namely, that the damping begins when the halo particles reach the resonant radius and that the damping is algebraic in time with nearly constant damping rate. The model also predicts a nonlinear frequency shift. The model provides two ways to think about the damping. It results from a transfer of canonical angular momentum from the mode to halo particles being swept by the mode field through the nonlinear cat's eye orbits of the resonant region. More mechanistically, the electric field produced by the perturbed charge density of the resonant particles acts back on the plasma core causing E × B drift that gives rise to the damping and nonlinear frequency shift.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

On the statistical theory of self-gravitating collisionless dark matter flow: Scale and redshift variation of velocity and density distributions

The statistics of velocity and density fields are crucial for cosmic structure formation and evolution. Here, this paper extends our previous work on the two-point second-order statistics for the velocity field [Phys. Fluids 35, 077105 (2023)] to one-point probability distributions for both density and velocity fields. The scale and redshift variation of density and velocity distributions are studied by a halo-based non-projection approach. First, all particles are divided into halo and out-of-halo particles so that the redshift variation can be studied via generalized kurtosis of distributions for halo and out-of-halo particles, respectively. Second, without projecting particle fields onto a structured grid, the scale variation is analyzed by identifying all particle pairs on different scales $r$. We demonstrate that: (i) Delaunay tessellation can be used to reconstruct the density field. The density correlation, spectrum, and dispersion functions were obtained, modeled, and compared with the N-body simulation; (ii) the velocity distributions are symmetric on both small and large scales and are non-symmetric with a negative skewness on intermediate scales due to the inverse energy cascade on small scales with a constant rate $\varepsilon_u$; (iii) On small scales, the even order moments of pairwise velocity $\Delta u_L$ follow a two-thirds law $\propto{(-\varepsilon_ur)}^{2/3}$, while the odd order moments follow a linear scaling $\langle(\Delta u_L)^{2n+1}\rangle=(2n+1)\langle(\Delta u_L)^{2n}\rangle\langle\Delta u_L\rangle\propto{r}$; (iv) The scale variation of the velocity distributions was studied for longitudinal velocities $u_L$ or $u_L^{'}$, pairwise velocity (velocity difference) $\Delta u_L$=$u_L^{'}$-$u_L$ and velocity sum $\Sigma u_L$=$u^{'}_L$+$u_L$. Fully developed velocity fields are never Gaussian on any scale, despite that they can initially be Gaussian; (v) On small scales, $u_L$ and $\Sigma u_L$ can be modeled by a $X$ distribution to maximize the entropy of the system. The distribution of $\Delta u_L$ can be different; (vi) On large scales, $\Delta u_L$ and $\Sigma u_L$ can be modeled by a logistic or a $X$ distribution, while $u_L$ has a different distribution; (vii) the redshift variation of the velocity distributions follows the evolution of the $X$ distribution involving a shape parameter $\alpha(z)$ decreasing with time.

79 ASTRONOMY AND ASTROPHYSICS↗

Multi-tracer intensity mapping: cross-correlations, line noise & decorrelation

Line intensity mapping (LIM) is a rapidly emerging technique for constraining cosmology and galaxy formation using multi-frequency, low angular resolution maps. Many LIM applications crucially rely on cross-correlations of two line intensity maps, or of intensity maps with galaxy surveys or galaxy/CMB lensing. We present a consistent halo model to predict all these cross-correlations and enable joint analyses, in 3D redshift-space and for 2D projected maps. We extend the conditional luminosity function formalism to the multi-line case, to consistently account for correlated scatter between multiple galaxy line luminosities. This allows us to model the scale-dependent decorrelation between two line intensity maps, a key input for foreground rejection and for approaches that estimate auto-spectra from cross-spectra. This also enables LIM cross-correlations to reveal astrophysical properties of the interstellar medium inacessible with LIM auto-spectra. We expose the different sources of luminosity scatter or "line noise" in LIM, and clarify their effects on the 1-halo and galaxy shot noise terms. In particular, we show that the effective number density of halos can in some cases exceed that of galaxies, counterintuitively. Using observational and simulation input, we implement this halo model for the Hα, [Oiii], Lyman-α, CO and [Cii] lines. Here, we encourage observers and simulators to measure galaxy luminosity correlation coefficients for pairs of lines whenever possible. Our code is publicly available at https://github.com/EmmanuelSchaan/HaloGen/tree/LIM. In a companion paper, we use this halo model formalism and code to highlight the degeneracies between cosmology and astrophysics in LIM, and to compare the LIM observables to galaxy detection for a number of surveys.

79 ASTRONOMY AND ASTROPHYSICS↗

Velocity-dependent annihilation radiation from dark matter subhalos in cosmological simulations

Here, we use the suite of Milky Way-like galaxies in the Auriga simulations to determine the contribution to annihilation radiation from dark matter subhalos in three velocity-dependent dark matter annihilation models: Sommerfeld, p-wave, and d-wave models. We compare these to the corresponding distribution in the velocity-independent s-wave annihilation model. For both the hydrodynamical and dark-matter-only simulations, only in the case of the Sommerfeld-enhanced annihilation does the total annihilation flux from subhalos exceed the total annihilation flux from the smooth halo component within the virial radius of the halo. Progressing from Sommerfeld to the s, p, and d-wave models, the contribution from the smooth component of the halo becomes more dominant, implying that for the p-wave and d-wave models the smooth component is by far the dominant contribution to the radiation. Comparing to the Galactic center excess observed by Fermi-LAT, for all simulated halos the emission is dominated by the smooth halo contribution. However, it is possible that for Sommerfeld models, extrapolation down to mass scales below the current resolution limit of the simulation would imply a non-negligible contribution to the gamma-ray emission from the Galactic Center region.

79 ASTRONOMY AND ASTROPHYSICS↗

Non-parametric Lagrangian biasing from the insights of neural nets

We present a Lagrangian model of galaxy clustering bias in which we train a neural net using the local properties of the smoothed initial density field to predict the late-time mass-weighted halo field.By fitting the mass-weighted halo field in the AbacusSummit simulations at z = 0.5, we find that including three coarsely spaced smoothing scales gives the best recovery of the halo power spectrum. Adding more smoothing scales may lead to 2–5% underestimation of the large-scale power and can cause the neural net to overfit.We find that the fitted halo-to-mass ratio can be well described by two directions in the original high-dimension feature space.Projecting the original features into these two principal components and re-training the neural net either reproduces the original training result, or outperforms it with a better match of the halo power spectrum. The elements of the principal components are unlikely to be assigned physical meanings, partly owing to the features being highly correlated between different smoothing scales.Our work illustrates a potential need to include multiple smoothing scales when studying galaxy bias, and this can be done easily with machine-learning methods that can take in high dimensional input feature space.

79 ASTRONOMY AND ASTROPHYSICS↗

Joint velocity and density reconstruction of the Universe with nonlinear differentiable forward modeling

Reconstructing the initial conditions of the Universe from late-time observations has the potential to optimally extract cosmological information. Due to the high dimensionality of the parameter space, a differentiable forward model is needed for convergence, and recent advances have made it possible to perform reconstruction with nonlinear models based on galaxy (or halo) positions. In addition to positions, future surveys will provide measurements of galaxies' peculiar velocities through the kinematic Sunyaev-Zel'dovich effect (kSZ), type Ia supernovae, the fundamental plane relation, and the Tully-Fisher relation. Here we develop the formalism for including halo velocities, in addition to halo positions, to enhance the reconstruction of the initial conditions. We show that using velocity information can significantly improve the reconstruction accuracy compared to using only the halo density field. Here we study this improvement as a function of shot noise, velocity measurement noise, and angle to the line of sight. We also show how halo velocity data can be used to improve the reconstruction of the final nonlinear matter overdensity and velocity fields. We have built our pipeline into the differentiable Particle-Mesh FlowPM package, paving the way to perform field-level cosmological inference with joint velocity and density reconstruction. This is especially useful given the increased ability to measure peculiar velocities in the near future.

79 ASTRONOMY AND ASTROPHYSICS↗