The evolution of the structural and kinematic properties of early-type galaxies (ETGs), their scaling relations, as well as their stellar metallicity and age contain precious information on the assembly history of these systems. We present results on the evolution of the stellar mass-velocity dispersion relation of ETGs, focusing in particular on the effects of some selection criteria used to define ETGs. We also try to shed light on the role that in-situ and ex-situ stellar populations have in massive ETGs, providing a possible explanation of the observed metallicity distributions.

A large body of work based on collisionless cosmological N-body simulations going back over two decades has advanced the idea that collapsed dark matter haloes have simple and approximately universal forms for their mass density and pseudo-phase space density (PPSD) distributions. However, a general consensus on the physical origin of these results has not yet been reached. In the present study, we explore to what extent the apparent universality of these forms holds when we vary the initial conditions (i.e., the primordial power spectrum of density fluctuations) away from the standard CMB-normalised case, but still within the context of LCDM with a fixed expansion history. Using simulations that vary the initial amplitude and shape, we show that the structure of dark matter haloes retains a clear memory of the initial conditions. Specifically, increasing (lowering) the amplitude of fluctuations increases (decreases) the concentration of haloes and, if pushed far enough, the density profiles deviate strongly from the NFW form that is a good approximation for the CMB-normalised case. Although, an Einasto form works well. Rather than being universal, the slope of the PPSD (or pseudo-entropy) profile steepens (flattens) with increasing (decreasing) power spectrum amplitude and can exhibit a strong halo mass dependence. Our results therefore indicate that the previously identified universality of the structure of dark matter haloes is mostly a consequence of adopting a narrow range of (CMB-normalised) initial conditions for the simulations. Our new suite provides a useful test-bench against which physical models for the origin of halo structure can be validated.

We investigate the influence of environment on the cold-gas properties of galaxies at z=0 within the TNG100 cosmological, magnetohydrodynamic simulation, part of the IllustrisTNG suite. We extend previous post-processing methods for breaking gas cells into their atomic and molecular phases, and build detailed mocks to comprehensively compare to the latest surveys of atomic hydrogen (HI) in nearby galaxies, namely ALFALFA and xGASS. We use TNG100 to explore the HI content, star formation activity, and angular momentum of satellite galaxies, each as a function of environment, and find that satellites are typically a factor of ~3 poorer in HI than centrals of the same stellar mass, with the exact offset depending sensitively on parent halo mass. Due to the large physical scales on which HI measurements are made (~45--245 kpc), contributions from gas not bound to the galaxy of interest but in the same line of sight crucially lead to larger HI mass measurements in the mocks in many cases, ultimately aligning with observations. This effect is mass-dependent and naturally greater for satellites than centrals, as satellites are never isolated by definition. We also show that HI stripping in TNG100 satellites is closely accompanied by quenching, in tension with observational data that instead favour that HI is preferentially stripped before star formation is reduced.

The splashback radius ($R_{\rm sp}$) of dark matter halos has recently been detected using weak gravitational lensing and cross-correlations with galaxies. However, different methods have been used to measure $R_{\rm sp}$ 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 $R_{\rm sp}$ by comparing it with the value measured directly from particle dynamics. We show that estimates of $R_{\rm sp}$ from density and subhalo profiles correspond to different percentiles of the underlying $R_{\rm sp}$ 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 splashback-accretion rate relation directly.

We analytically derive the observed size-mass relation of galaxies' atomic hydrogen (HI), including limits on its scatter, based on simple assumptions about the structure of HI discs. We trial three generic profiles for HI surface density as a function of radius. Firstly, we assert that HI surface densities saturate at a variable threshold, and otherwise fall off exponentially with radius or, secondly, radius squared. Our third model assumes the total gas surface density is exponential, with the HI fraction at each radius depending on local pressure. These are tested against a compilation of 110 galaxies from the THINGS, LITTLE THINGS, LVHIS, and Bluedisk surveys, whose HI surface density profiles are well resolved. All models fit the observations well and predict consistent size-mass relations. Using an analytical argument, we explain why processes that cause gas disc truncation - such as ram-pressure stripping - scarcely affect the HI size-mass relation. This is tested with the IllustrisTNG(100) cosmological, hydrodynamic simulation and the Dark Sage semi-analytic model of galaxy formation, both of which capture radially resolved disc structure. For galaxies with m_*>10^9 M_solar and m_HI>10^8 M_solar, both simulations predict HI size-mass relations that align with observations, show no difference between central and satellite galaxies, and show only a minor, second-order dependence on host halo mass for satellites. Ultimately, the universally tight HI size-mass relation is mathematically inevitable and robust. Only by completely disrupting the structure of HI discs, e.g. through overly powerful feedback, could a simulation predict the relation poorly.

We are just starting to understand the physical processes driving the dramatic change in cosmic star-formation rate between $z\sim 2$ and the present day. A quantity directly linked to star formation is the molecular gas density, which should be measured through independent methods to explore variations due to cosmic variance and systematic uncertainties. We use intervening CO absorption lines in the spectra of mm-bright background sources to provide a census of the molecular gas mass density of the Universe. The data used in this work are taken from ALMACAL, a wide and deep survey utilizing the ALMA calibrator archive. While we report multiple Galactic absorption lines and one intrinsic absorber, no extragalactic intervening molecular absorbers are detected. However, thanks to the large redshift path surveyed ($\Delta z=182$), we provide constraints on the molecular column density distribution function beyond $z\sim 0$. In addition, we probe column densities of N(H$_2$) > 10$^{16}$ atoms~cm$^{-2}$, five orders of magnitude lower than in previous studies. We use the cosmological hydrodynamical simulation IllustrisTNG to show that our upper limits of $\rho ({\rm H}_2)\lesssim 10^{8.3} \text{M}_{\odot} \text{Mpc}^{-3}$ at $0 < z \leq 1.7$ already provide new constraints on current theoretical predictions of the cold molecular phase of the gas. These results are in agreement with recent CO emission-line surveys and are complementary to those studies. The combined constraints indicate that the present decrease of the cosmic star-formation rate history is consistent with an increasing depletion of molecular gas in galaxies compared to $z\sim 2$.

Using the IllustrisTNG simulations, we investigate the connection between galaxy morphology and star formation in central galaxies with stellar masses in the range $10^9-10^{11.5}~\mathrm{M}_{\odot}$. We quantify galaxy morphology by a kinematical decomposition of the stellar component into a spheroidal and a disc component (spheroid-to-total ratio, S/T) and by the concentration of the stellar mass density profile ($C_{82}$). S/T is correlated with stellar mass and star-formation activity, while $C_{82}$ correlates only with stellar mass. Overall, we find good agreement with observational estimates for both S/T and $C_{82}$. Low and high mass galaxies are dominated by random stellar motion, while only intermediate-mass galaxies ($M_{\star}\approx10^{10}-10^{10.5}~\mathrm{M}_{\odot}$) are dominated by ordered rotation. Whereas higher-mass galaxies are typical spheroids with high concentrations, lower-mass galaxies have low concentration, pointing to different formation channels. Although we find a correlation between S/T and star-formation activity, in the TNG model galaxies do not necessarily change their morphology when they transition through the green valley or when they cease their star formation, this depending on galaxy stellar mass and morphological estimator. Instead, the morphology (S/T and $C_{82}$) is generally set during the star-forming phase of galaxies. The apparent correlation between S/T and star formation arises because earlier-forming galaxies had, on average, a higher S/T at a given stellar mass. Furthermore, we show that mergers drive in-situ bulge formation in intermediate-mass galaxies and are responsible for the recent spheroidal mass assembly in the massive galaxies with $M_{\star}>10^{11}~\mathrm{M}_{\odot}$. In particular, these massive galaxies assemble about half of the spheroidal mass while star-forming and the other half through mergers while quiescent.

We have recently developed a post-processing framework to estimate the abundance of atomic and molecular hydrogen (HI and H2, respectively) in galaxies in large-volume cosmological simulations. Here we compare the HI and H2 content of IllustrisTNG galaxies to observations. We mostly restrict this comparison to $z \approx 0$ and consider six observational metrics: the overall abundance of HI and H2, their mass functions, gas fractions as a function of stellar mass, the correlation between H2 and star formation rate, the spatial distribution of gas, and the correlation between gas content and morphology. We find generally good agreement between simulations and observations, particularly for the gas fractions and the HI mass-size relation. The H2 mass correlates with star formation rate as expected, revealing an almost constant depletion time that evolves up to z = 2 as observed. However, we also discover a number of tensions with varying degrees of significance, including an overestimate of the total neutral gas abundance at z = 0 by about a factor of two and a possible excess of satellites with no or very little neutral gas. These conclusions are robust to the modelling of the HI/H2 transition. In terms of their neutral gas properties, the IllustrisTNG simulations represent an enormous improvement over the original Illustris run. All data used in this paper are publicly available as part of the IllustrisTNG data release.

We present the full public release of all data from the TNG100 and TNG300 simulations of the IllustrisTNG project. IllustrisTNG is a suite of large volume, cosmological, gravo-magnetohydrodynamical simulations run with the moving-mesh code Arepo. TNG includes a comprehensive model for galaxy formation physics, and each TNG simulation self-consistently solves for the coupled evolution of dark matter, cosmic gas, luminous stars, and supermassive blackholes from early time to the present day, z=0. Each of the flagship runs -- TNG50, TNG100, and TNG300 -- are accompanied by lower-resolution and dark-matter only counterparts, and we discuss scientific and numerical cautions and caveats relevant when using TNG. Full volume snapshots are available at 100 redshifts; halo and subhalo catalogs at each snapshot and merger trees are also released. The data volume now directly accessible online is ~750 TB, including 1200 full volume snapshots and ~80,000 high time-resolution subbox snapshots. This will increase to ~1.1 PB with the future release of TNG50. Data access and analysis examples are available in IDL, Python, and Matlab. We describe improvements and new functionality in the web-based API, including on-demand visualization and analysis of galaxies and halos, exploratory plotting of scaling relations and other relationships between galactic and halo properties, and a new JupyterLab interface. This provides an online, browser-based, near-native data analysis platform which supports user computation with fully local access to TNG data, alleviating the need to download large simulated datasets.

The ALMA Spectroscopic Survey in the Hubble Ultra Deep Field (ASPECS) provides new constraints for galaxy formation models on the molecular gas properties of galaxies. We compare results from ASPECS to predictions from two cosmological galaxy formation models: the IllustrisTNG hydrodynamical simulations and the Santa Cruz semi-analytic model (SC SAM). We explore several recipes to model the H$_2$ content of galaxies, finding them to be consistent with one another, and take into account the sensitivity limits and survey area of ASPECS. For a canonical CO-to-H$_2$ conversion factor of $\alpha_{\rm CO} = 3.6\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$ the results of our work include: (1) the H$_2$ mass of $z>1$ galaxies predicted by the models as a function of their stellar mass is a factor of 2-3 lower than observed; (2) the models do not reproduce the number of H$_2$-rich ($M_{\rm H2} > 3\times 10^{10}\,\rm{M}_\odot$) galaxies observed by ASPECS; (3) the H$_2$ cosmic density evolution predicted by IllustrisTNG (the SC SAM) is in tension (only just agrees) with the observed cosmic density, even after accounting for the ASPECS selection function and field-to-field variance effects. The tension between models and observations at $z>1$ can be alleviated by adopting a CO-to-H$_2$ conversion factor in the range $\alpha_{\rm CO} = 2.0 - 0.8\,\rm{M}_\odot/(\rm{K}\,\rm{km/s}\,\rm{pc}^{2})$. Additional work on constraining the CO-to-H$_2$ conversion factor and CO excitation conditions of galaxies through observations and theory will be necessary to more robustly test the success of galaxy formation models.

We write in response to the call from the 2020 Decadal Survey to submit white papers illustrating the most pressing scientific questions in astrophysics for the coming decade. We propose exploration as the central question for the Decadal Committee's discussions.The history of astronomy shows that paradigm changing discoveries are not driven by well formulated scientific questions, based on the knowledge of the time. They were instead the result of the increase in discovery space fostered by new telescopes and instruments. An additional tool for increasing the discovery space is provided by the analysis and mining of the increasingly larger amount of archival data available to astronomers. Revolutionary observing facilities, and the state of the art astronomy archives needed to support these facilities, will open up the universe to new discovery. Here we focus on exploration for compact objects and multi messenger science. This white paper includes science examples of the power of the discovery approach, encompassing all the areas of astrophysics covered by the 2020 Decadal Survey.

The relation between halo mass, M, and concentration, c, is a critical component in our understanding of the structure of dark matter halos. While numerous models for this relation have been proposed, almost none of them attempt to derive the evolution of the relation analytically. We build on previous efforts to model the c-M relation as a function of physical parameters such as the peak height, $\nu$, and the effective power spectrum slope, $n_{\rm eff}$, which capture the dependence of $c$ on halo mass, redshift, and cosmology. We present three major improvements over previous models. First, we derive an analytical expression for the c-M relation that is valid under the assumption of pseudo-evolution, i.e., assuming that the density profiles of halos are static in physical coordinates while the definition of their boundary evolves. We find that this ansatz is highly successful in describing the evolution of the low-mass end of the c-M relation. Second, we employ a new physical variable, the effective exponent of linear growth, $\alpha_{\rm eff}$, to parameterize deviations from an Einstein-de Sitter expansion history. Third, we combine an updated definition of $n_{\rm eff}$ with the additional dependence on $\alpha_{\rm eff}$ and propose a phenomenological extension of our analytical framework to include all halo masses. This semianalytical model matches simulated concentrations in both scale-free models and LambdaCDM to 5% accuracy with very few exceptions and differs significantly from all previously proposed models. We present a publicly available code to compute the predictions of our model in the python toolkit Colossus, including updated parameters for the model of Diemer and Kravtsov.

This paper introduces Colossus, a public, open-source python package for calculations related to cosmology, the large-scale structure (LSS) of matter in the universe, and the properties of dark matter halos. The code is designed to be fast and easy to use, with a coherent, well-documented user interface. The cosmology module implements Friedman-Lemaitre-Robertson-Walker cosmologies including curvature, relativistic species, and different dark energy equations of state, and provides fast computations of the linear matter power spectrum, variance, and correlation function. The LSS module is concerned with the properties of peaks in Gaussian random fields and halos in a statistical sense, including their peak height, peak curvature, halo bias, and mass function. The halo module deals with spherical overdensity radii and masses, density profiles, concentration, and the splashback radius. To facilitate the rapid exploration of these quantities, Colossus implements more than 40 different fitting functions from the literature. I discuss the core routines in detail, with particular emphasis on their accuracy. Colossus is available at bitbucket.org/bdiemer/colossus.

Large-scale cosmological simulations of galaxy formation currently do not resolve the densities at which molecular hydrogen forms, implying that the atomic-to-molecular transition must be modeled either on the fly or in postprocessing. We present an improved postprocessing framework to estimate the abundance of atomic and molecular hydrogen and apply it to the IllustrisTNG simulations. We compare five different models for the atomic-to-molecular transition, including empirical, simulation-based, and theoretical prescriptions. Most of these models rely on the surface density of neutral hydrogen and the ultraviolet (UV) flux in the Lyman-Werner band as input parameters. Computing these quantities on the kiloparsec scales resolved by the simulations emerges as the main challenge. We show that the commonly used Jeans length approximation to the column density of a system can be biased and exhibits large cell-to-cell scatter. Instead, we propose to compute all surface quantities in face-on projections and perform the modeling in two dimensions. In general, the two methods agree on average, but their predictions diverge for individual galaxies and for models based on the observed midplane pressure of galaxies. We model the UV radiation from young stars by assuming a constant escape fraction and optically thin propagation throughout the galaxy. With these improvements, we find that the five models for the atomic-to-molecular transition roughly agree on average but that the details of the modeling matter for individual galaxies and the spatial distribution of molecular hydrogen. We emphasize that the estimated molecular fractions are approximate due to the significant systematic uncertainties.

[Abridged] We study the abundance and clustering properties of HI at redshifts $z\leqslant5$ using TNG100, a large state-of-the-art magneto-hydrodynamic simulation of a 75 Mpc/h box size. We show that most of the HI lies within dark matter halos and quantify the average HI mass hosted by halos of mass M at redshift z. We find that only halos with circular velocities larger than $\simeq$ 30 km/s contain HI. While the density profiles of HI exhibit a large halo-to-halo scatter, the mean profiles are universal across mass and redshift. The HI in low-mass halos is mostly located in the central galaxy, while in massive halos is concentrated in the satellites. We show that the HI and matter density probability distribution functions differ significantly. Our results point out that for small halos the HI bulk velocity goes in the same direction and has the same magnitude as the halo peculiar velocity, while in large halos differences show up. We find that halo HI velocity dispersion follows a power-law with halo mass. We find a complicated HI bias, with HI becoming non-linear already at $k=0.3$ h/Mpc at $z\gtrsim3$. Our simulation reproduces the DLAs bias value from observations. We find that the clustering of HI can be accurately reproduced by perturbative methods. We identify a new secondary bias, by showing that the clustering of halos depends not only on mass but also on HI content. We compute the amplitude of the HI shot-noise and find that it is small at all redshifts. We study the clustering of HI in redshift-space, and show that linear theory can explain the ratio between the monopoles in redshift- and real-space down to small scales at high redshift. We find that the amplitude of the Fingers-of-God effect is larger for HI than for matter. We point out that accurate 21 cm maps can be created from N-body or approximate simulations rather than full hydrodynamic simulations.

We present the first stable release of Halotools (v0.2), a community-driven Python package designed to build and test models of the galaxy-halo connection. Halotools provides a modular platform for creating mock universes of galaxies starting from a catalog of dark matter halos obtained from a cosmological simulation. The package supports many of the common forms used to describe galaxy-halo models: the halo occupation distribution (HOD), the conditional luminosity function (CLF), abundance matching, and alternatives to these models that include effects such as environmental quenching or variable galaxy assembly bias. Satellite galaxies can be modeled to live in subhalos, or to follow custom number density profiles within their halos, including spatial and/or velocity bias with respect to the dark matter profile. The package has an optimized toolkit to make mock observations on a synthetic galaxy population, including galaxy clustering, galaxy-galaxy lensing, galaxy group identification, RSD multipoles, void statistics, pairwise velocities and others, allowing direct comparison to observations. Halotools is object-oriented, enabling complex models to be built from a set of simple, interchangeable components, including those of your own creation. Halotools has an automated testing suite and is exhaustively documented on http://halotools.readthedocs.io, which includes quickstart guides, source code notes and a large collection of tutorials. The documentation is effectively an online textbook on how to build and study empirical models of galaxy formation with Python.

Dark matter halo clustering depends not only on halo mass, but also on other properties such as concentration and shape. This phenomenon is known broadly as assembly bias. We explore the dependence of assembly bias on halo definition, parametrized by spherical overdensity parameter, $\Delta$. We summarize the strength of concentration-, shape-, and spin-dependent halo clustering as a function of halo mass and halo definition. Concentration-dependent clustering depends strongly on mass at all $\Delta$. For conventional halo definitions ($\Delta \sim 200\mathrm{m}-600\mathrm{m}$), concentration-dependent clustering at low mass is driven by a population of haloes that is altered through interactions with neighbouring haloes. Concentration-dependent clustering can be greatly reduced through a mass-dependent halo definition with $\Delta \sim 20\mathrm{m}-40\mathrm{m}$ for haloes with $M_{200\mathrm{m}} \lesssim 10^{12}\, h^{-1}\mathrm{M}_{\odot}$. Smaller $\Delta$ implies larger radii and mitigates assembly bias at low mass by subsuming altered, so-called backsplash haloes into now larger host haloes. At higher masses ($M_{200\mathrm{m}} \gtrsim 10^{13}\, h^{-1}\mathrm{M}_{\odot}$) larger overdensities, $\Delta \gtrsim 600\mathrm{m}$, are necessary. Shape- and spin-dependent clustering are significant for all halo definitions that we explore and exhibit a relatively weaker mass dependence. Generally, both the strength and the sense of assembly bias depend on halo definition, varying significantly even among common definitions. We identify no halo definition that mitigates all manifestations of assembly bias. A halo definition that mitigates assembly bias based on one halo property (e.g., concentration) must be mass dependent. The halo definitions that best mitigate concentration-dependent halo clustering do not coincide with the expected average splashback radii at fixed halo mass.

Motivated by the recent proposal of the splashback radius as a physical boundary of dark matter halos, we present a parallel computer code for Subhalo and PARticle Trajectory Analysis (SPARTA). The code analyzes the orbits of all simulation particles in all host halos, billions of orbits in the case of typical cosmological N-body simulations. Within this general framework, we develop an algorithm that accurately extracts the location of the first apocenter of particles after infall into a halo, or splashback. We define the splashback radius of a halo as the smoothed average of the apocenter radii of individual particles. This definition allows us to reliably measure the splashback radii of 95% of host halos above a resolution limit of 1000 particles. We show that, on average, the splashback radius and mass are converged to better than 5% accuracy with respect to mass resolution, snapshot spacing, and all free parameters of the method.

The splashback radius $R_{\rm sp}$, the apocentric radius of particles on their first orbit after falling into a dark matter halo, has recently been suggested as a physically motivated halo boundary that separates accreting from orbiting material. Using the SPARTA code presented in Paper I, we analyze the orbits of billions of particles in cosmological simulations of structure formation and measure $R_{\rm sp}$ for a large sample of halos that span a mass range from dwarf galaxy to massive cluster halos, reach redshift 8, and include WMAP, Planck, and self-similar cosmologies. We analyze the dependence of $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ on the mass accretion rate $\Gamma$, halo mass, redshift, and cosmology. The scatter in these relations varies between 0.02 and 0.1 dex. While we confirm the known trend that $R_{\rm sp}/R_{\rm 200m}$ decreases with $\Gamma$, the relationships turn out to be more complex than previously thought, demonstrating that $R_{\rm sp}$ is an independent definition of the halo boundary that cannot trivially be reconstructed from spherical overdensity definitions. We present fitting functions for $R_{\rm sp}/R_{\rm 200m}$ and $M_{\rm sp}/M_{\rm 200m}$ as a function of accretion rate, peak height, and redshift, achieving an accuracy of 5% or better everywhere in the parameter space explored. We discuss the physical meaning of the distribution of particle apocenters and show that the previously proposed definition of $R_{\rm sp}$ as the radius of the steepest logarithmic density slope encloses roughly three-quarters of the apocenters. Finally, we conclude that no analytical model presented thus far can fully explain our results.

The density field in the outskirts of dark matter halos is discontinuous due to a caustic formed by matter at its first apocenter after infall. In this paper, we present an algorithm to identify the "splashback shell" formed by these apocenters in individual simulated halos using only a single snapshot of the density field. We implement this algorithm in the code SHELLFISH (SHELL Finding In Spheroidal Halos) and demonstrate that the code identifies splashback shells correctly and measures their properties with an accuracy of $<5\%$ for halos with more than 50,000 particles and mass accretion rates of $\Gamma_\textrm{DK14}>0.5$. Using SHELLFISH, we present the first estimates for several basic properties of individual splashback shells, such as radius, $R_\textrm{sp}$, mass, and overdensity, and provide fits to the distribution of these quantities as functions of $\Gamma_\textrm{DK14}$, $\nu_\textrm{200m}$, and $z.$ We confirm previous findings that $R_\textrm{sp}$ decreases with increasing $\Gamma_\textrm{DK14}$, but we show that independent of accretion rate, it also decreases with increasing $\nu_\textrm{200m}$. We also study the 3D structures of these shells and find that they generally have non-ellipsoidal oval shapes. We find that splashback radii estimated by SHELLFISH are $20\%-30\%$ larger than those estimated in previous studies from stacked density profiles at high accretion rates. We demonstrate that the latter are biased low due to the contribution of high-mass subhalos to these profiles and show that using the median instead of mean density in each radial bin mitigates the effect of substructure on density profiles and removes the bias.

We introduce The Fabric of the Universe, an art and science collaboration focused on exploring the cosmic web of dark matter with unconventional techniques and materials. We discuss two of our projects in detail. First, we describe a pipeline for translating three-dimensional (3D) density structures from N-body simulations into solid surfaces suitable for 3D printing, and present prints of a cosmological volume and of the infall region around a massive cluster halo. In these models, we discover wall-like features that are invisible in two-dimensional projections. Going beyond the sheer visualization of simulation data, we undertake an exploration of the cosmic web as a three-dimensional woven textile. To this end, we develop experimental 3D weaving techniques to create sphere-like and filamentary shapes and radically simplify a region of the cosmic web into a set of filaments and halos. We translate the resulting tree structure into a series of commands that can be executed by a digital weaving machine, and present a large-scale textile installation.

Gladders et al. have recently suggested that the star formation histories (SFHs) of individual galaxies are characterized by a log-normal function in time, implying a slow decline rather than rapid quenching. We test their conjecture on theoretical SFHs from the cosmological simulation Illustris and on observationally inferred SFHs. While the log-normal form necessarily ignores short-lived features such as starbursts, it fits the overall shape of the majority of SFHs very well. In particular, 85% of the cumulative SFHs are fitted to within a maximum error of 5% of the total stellar mass formed, and 99% to within 10%. The log-normal performs systematically better than the commonly used delayed-tau model, and is superseded only by functions with more than three free parameters. Poor fits are mostly found in galaxies that were rapidly quenched after becoming satellites. We explore the log-normal parameter space of normalization, peak time, and full width at half maximum, and find that the simulated and observed samples occupy similar regions, though Illustris predicts wider, later-forming SFHs on average. The ensemble of log-normal fits correctly reproduces complex metrics such as the evolution of Illustris galaxies across the star formation main sequence, but overpredicts their quenching timescales. SFHs in Illustris are a diverse population not determined by any one physical property of galaxies, but follow a tight relation, where $\mathrm{width}\propto\mathrm{(peak\ time)}^{3/2}$. We show that such a relation can be explained qualitatively (though not quantitatively) by a close connection between the growth of dark matter halos and their galaxies.

We carry out a systematic investigation of the total mass density profile of massive (Mstar>2e11 Msun) early-type galaxies and its dependence on galactic properties and host halo mass with the aid of a variety of lensing/dynamical data and large mock galaxy catalogs. The latter are produced via semi-empirical models that, by design, are based on just a few basic input assumptions. Galaxies, with measured stellar masses, effective radii and S\'{e}rsic indices, are assigned, via abundance matching relations, host dark matter halos characterized by a typical LCDM profile. Our main results are as follows: (i) In line with observational evidence, our semi-empirical models naturally predict that the total, mass-weighted density slope at the effective radius gamma' is not universal, steepening for more compact and/or massive galaxies, but flattening with increasing host halo mass. (ii) Models characterized by a Salpeter or variable initial mass function and uncontracted dark matter profiles are in good agreement with the data, while a Chabrier initial mass function and/or adiabatic contractions/expansions of the dark matter halos are highly disfavored. (iii) Currently available data on the mass density profiles of very massive galaxies (Mstar>1e12 Msun), with Mhalo>3e14 Msun, favor instead models with a stellar profile flatter than a S\'{e}rsic one in the very inner regions (r<3-5 kpc), and a cored NFW or Einasto dark matter profile with median halo concentration a factor of ~2 or <1.3, respectively, higher than those typically predicted by N-body numerical simulations.

The lensing signal around galaxy clusters can, in principle, be used to test detailed predictions for their average mass profile from numerical simulations. However, the intrinsic shape of the profiles can be smeared out when a sample that spans a wide range of cluster masses is averaged in physical length units. This effect especially conceals rapid changes in gradient such as the steep drop associated with the splashback radius, a sharp edge corresponding to the outermost caustic in accreting halos. We optimize the extraction of such local features by scaling individual halo profiles to a number of spherical overdensity radii, and apply this method to 16 X-ray-selected high-mass clusters targeted in the Cluster Lensing And Supernova survey with Hubble. By forward-modeling the weak- and strong-lensing data presented by Umetsu et al., we show that, regardless of the scaling overdensity, the projected ensemble density profile is remarkably well described by an NFW or Einasto profile out to $R \sim 2.5h^{-1}$Mpc, beyond which the profiles flatten. We constrain the NFW concentration to $c_{200c} = 3.66 \pm 0.11$ at $M_{200c} \simeq 1.0 \times 10^{15}h^{-1}M_\odot$, consistent with and improved from previous work that used conventionally stacked lensing profiles, and in excellent agreement with theoretical expectations. Assuming the profile form of Diemer & Kravtsov and generic priors calibrated from numerical simulations, we place a lower limit on the splashback radius of the cluster halos, if it exists, of $R_{sp}/r_{200m} > 0.89$ ($R_{sp} > 1.83h^{-1}$Mpc) at 68% confidence. The corresponding density feature is most pronounced when the cluster profiles are scaled by $r_{200m}$, and smeared out when scaled to higher overdensities.

The fact that the clustering of dark matter halos depends not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies.

We show that the projected number density profiles of SDSS photometric galaxies around galaxy clusters displays strong evidence for the splashback radius, a sharp halo edge corresponding to the location of the first orbital apocenter of satellite galaxies after their infall. We split the clusters into two subsamples with different mean projected radial distances of their members, $\langle R_{\rm mem}\rangle$, at fixed richness and redshift, and show that the sample with smaller $\langle R_{\rm mem}\rangle$ has a smaller ratio of the splashback radius to the traditional halo boundary $R_{\rm 200m}$, than the subsample with larger $\langle R_{\rm mem}\rangle$, indicative of different mass accretion rates for the two subsamples. The same cluster samples were recently used by Miyatake et al. to show that their large-scale clustering differs despite their similar weak lensing masses, demonstrating strong evidence for halo assembly bias. We expand on this result by presenting a 6.6-$\sigma$ detection of halo assembly bias using the cluster-photometric galaxy cross-correlations. Our measured splashback radii are smaller, while the strength of the assembly bias signal is stronger, than expectations from N-body simulations based on the $\Lambda$-dominated, cold dark matter structure formation model. Dynamical friction or cluster-finding systematics such as miscentering or projection effects are not likely to be the sole source of these discrepancies.

The boundaries of cold dark matter halos are commonly defined to enclose a density contrast $\Delta$ relative to a reference (mean or critical) density. We argue that a more physical boundary of halos is the radius at which accreted matter reaches its first orbital apocenter after turnaround. This splashback radius, $R_{sp}$, manifests itself as a sharp density drop in the halo outskirts, at a location that depends upon the mass accretion rate. We present calibrations of $R_{sp}$ and the enclosed mass, $M_{sp}$, as a function of the accretion rate and alternatively peak height. We find that $R_{sp}$ varies between $\approx0.8-1R_{200m}$ for rapidly accreting halos and $\approx1.5R_{200m}$ for slowly accreting halos. The extent of a halo and its associated environmental effects can thus extend well beyond the conventionally defined "virial" radius. We show that $M_{sp}$ and $R_{sp}$ evolve relatively strongly compared to other commonly used definitions. In particular, $M_{sp}$ evolves significantly even for the smallest dwarf-sized halos at $z=0$. We also contrast $M_{sp}$ with the mass enclosed within four scale radii of the halo density profile, $M_{<4rs}$, which characterizes the inner halo. During the early stages of halo assembly, $M_{sp}$ and $M_{<4rs}$ evolve similarly, but in the late stages $M_{<4rs}$ stops increasing while $M_{sp}$ continues to grow significantly. This illustrates that halos at low $z$ can have "quiet" interiors while continuing to accrete mass in their outskirts. We discuss potential observational estimates of the splashback radius and show that it may already have been detected in galaxy clusters.

We present a numerical study of dark matter halo concentrations in $\Lambda$CDM and self-similar cosmologies. We show that the relation between concentration, $c$, and peak height, $\nu$, exhibits the smallest deviations from universality if halo masses are defined with respect to the critical density of the universe. These deviations can be explained by the residual dependence of concentration on the local slope of the matter power spectrum, $n$, which affects both the normalization and shape of the $c$-$\nu$ relation. In particular, there is no well-defined floor in the concentration values. Instead, the minimum concentration depends on redshift: at fixed $\nu$, halos at higher $z$ experience steeper slopes $n$, and thus have lower minimum concentrations. We show that the concentrations in our simulations can be accurately described by a universal seven-parameter function of only $\nu$ and $n$. This model matches our $\Lambda$CDM results to $\lesssim 5\%$ accuracy up to $z = 6$, and matches scale-free $\Omega_{\rm m} = 1$ models to $\lesssim 15\%$. The model also reproduces the low concentration values of Earth--mass halos at $z \approx 30$, and thus correctly extrapolates over $16$ orders of magnitude in halo mass. The predictions of our model differ significantly from all models previously proposed in the literature at high masses and redshifts. Our model is in excellent agreement with recent lensing measurements of cluster concentrations.

We present a systematic study of the density profiles of LCDM halos, focusing on the outer regions, 0.1 < r/Rvir < 9. We show that the median and mean profiles of halo samples of a given peak height exhibit significant deviations from the universal analytic profiles discussed previously in the literature, such as the Navarro-Frenk-White and Einasto profiles, at radii r > 0.5 R200m. In particular, at these radii the logarithmic slope of the median density profiles of massive or rapidly accreting halos steepens more sharply than predicted. The steepest slope of the profiles occurs at r ~ R200m, and its absolute value increases with increasing peak height or mass accretion rate, reaching slopes of -4 and steeper. Importantly, we find that the outermost density profiles at r > R200m are remarkably self-similar when radii are rescaled by R200m. This self-similarity indicates that radii defined with respect to the mean density are preferred for describing the structure and evolution of the outer profiles. However, the inner density profiles are most self-similar when radii are rescaled by R200c. We propose a new fitting formula that describes the median and mean profiles of halo samples selected by their peak height or mass accretion rate with accuracy < 10% at all radii, redshifts and masses we studied, r < 9 Rvir, 0 < z < 6 and Mvir > 1.7E10 Msun/h. We discuss observational signatures of the profile features described above, and show that the steepening of the outer profile should be detectable in future weak-lensing analyses of massive clusters. Such observations could be used to estimate the mass accretion rate of cluster halos.

We present a systematic study of the pure deflagration model of Type Ia supernovae using three-dimensional, high-resolution, full-star hydrodynamical simulations, nucleosynthetic yields calculated using Lagrangian tracer particles, and light curves calculated using radiation transport. We evaluate the simulations by comparing their predicted light curves with many observed SNe Ia using the SALT2 data-driven model and find that the simulations may correspond to under-luminous SNe Iax. We explore the effects of the initial conditions on our results by varying the number of randomly selected ignition points from 63 to 3500, and the radius of the centered sphere they are confined in from 128 to 384 km. We find that the rate of nuclear burning depends on the number of ignition points at early times, the density of ignition points at intermediate times, and the radius of the confining sphere at late times. The results depend primarily on the number of ignition points, but we do not expect this to be the case in general. The simulations with few ignition points release more nuclear energy $E_{\mathrm{nuc}}$, have larger kinetic energies $E_{\mathrm{K}}$, and produce more $^{56}$Ni than those with many ignition points, and differ in the distribution of $^{56}$Ni, Si, and C/O in the ejecta. For these reasons, the simulations with few ignition points exhibit higher peak B-band absolute magnitudes $M_\mathrm{B}$ and light curves that rise and decline more quickly; their $M_\mathrm{B}$ and light curves resemble those of under-luminous SNe Iax, while those for simulations with many ignition points are not.

Understanding the evolution of scaling relations between the observable properties of clusters and their total mass is key to realizing their potential as cosmological probes. In this study, we investigate whether the evolution of cluster scaling relations is affected by the spurious evolution of mass caused by the evolving reference density with respect to which halo masses are defined (pseudo-evolution). We use the relation between mass, M, and velocity dispersion, sigma, as a test case, and show that the deviation from the M-sigma relation of cluster-sized halos caused by pseudo-evolution is smaller than 10% for a wide range of mass definitions. The reason for this small impact is a tight relation between the velocity dispersion and mass profiles, sigma(<r) = const * (GM(<r) / r)^1/2, which holds across a wide range of radii. We show that such a relation is generically expected for a variety of density profiles, as long as halos are in approximate Jeans equilibrium. Thus, as the outer "virial" radius used to define the halo mass, R, increases due to pseudo-evolution, halos approximately preserve their M-sigma relation. This result highlights the fact that tight scaling relations are the result of tight equilibrium relations between radial profiles of physical quantities. We find exceptions at very small and very large radii, where the profiles deviate from the relations they exhibit at intermediate radii. We discuss the implications of these results for other cluster scaling relations, and argue that pseudo-evolution should have a small effect on most scaling relations, except for those that involve the stellar masses of galaxies. In particular, we show that the relation between stellar-mass fraction and total mass is affected by pseudo-evolution and is largely shaped by it for halo masses of less than 10^14 solar masses.

We propose a robust, quantitative method to compare the synthetic light curves of a Type Ia Supernova (SNIa) explosion model with a large set of observed SNeIa, and derive a figure of merit for the explosion model's agreement with observations. The synthetic light curves are fit with the data-driven model SALT2 which returns values for stretch, color, and magnitude at peak brightness, as well as a goodness-of-fit parameter. Each fit is performed multiple times with different choices of filter bands and epoch range in order to quantify the systematic uncertainty on the fitted parameters. We use a parametric population model for the distribution of observed SNIa parameters from large surveys, and extend it to represent red, dim, and bright outliers found in a low-redshift SNIa data set. We discuss the potential uncertainties of this population model and find it to be reliable given the current uncertainties on cosmological parameters. Using our population model, we assign each set of fitted parameters a likelihood of being observed in nature, and a figure of merit based on this likelihood. We define a second figure of merit based on the quality of the light curve fit, and combine the two measures into an overall figure of merit for each explosion model. We compute figures of merit for a variety of 1D, 2D and 3D explosion models and show that our evaluation method allows meaningful inferences across a wide range of light curve quality and fitted parameters.

A dark matter halo is commonly defined as a spherical overdensity of matter with respect to a reference density, such as the critical density or the mean matter density of the Universe. Such definitions can lead to a spurious pseudo-evolution of halo mass simply due to redshift evolution of the reference density, even if its physical density profile remains constant over time. We estimate the amount of such pseudo-evolution of mass between z=1 to 0 for halos identified in a large N-body simulation, and show that it accounts for almost the entire mass evolution of the majority of halos with M200 of about 1E12 solar masses and can be a significant fraction of the apparent mass growth even for cluster-sized halos. We estimate the magnitude of the pseudo-evolution assuming that halo density profiles remain static in physical coordinates, and show that this simple model predicts the pseudo-evolution of halos identified in numerical simulations to good accuracy, albeit with significant scatter. We discuss the impact of pseudo-evolution on the evolution of the halo mass function and show that the non-evolution of the low-mass end of the halo mass function is the result of a fortuitous cancellation between pseudo-evolution and the absorption of small halos into larger hosts. We also show that the evolution of the low mass end of the concentration-mass relation observed in simulations is almost entirely due to the pseudo-evolution of mass. Finally, we discuss the implications of our results for the interpretation of the evolution of various scaling relations between the observable properties of galaxies and galaxy clusters and their halo masses.