ABSTRACT We present astrolink, an efficient and versatile clustering algorithm designed to hierarchically classify astrophysically relevant structures from both synthetic and observational data sets. We build upon clustar-nd, a hierarchical galaxy/(sub)halo finder, so that astrolink now generates a 2D representation of the implicit clustering structure as well as ensuring that clusters are statistically distinct from the noisy density fluctuations implicit within the n-dimensional input data. This redesign replaces the three cluster extraction parameters from clustar-nd with a single parameter, S – the lower statistical significance threshold of clusters, which can be automatically and reliably estimated via a dynamical model-fitting process. We demonstrate the robustness of this approach compared to astrolink’s predecessors by applying each algorithm to a suite of simulated galaxies defined over various feature spaces. We find that astrolink delivers a more powerful clustering performance while being $\sim 27~{{\ \rm per \, cent}}$ faster and using less memory than clustar-nd. With these improvements, astrolink is ideally suited to extracting a meaningful set of hierarchical and arbitrarily shaped astrophysical clusters from both synthetic and observational data sets – lending itself as a great tool for morphological decomposition within the context of hierarchical structure formation.

It has been shown in previous publications that the TNG100 simulation quantitatively reproduces the observed reduction in each of the total atomic and total molecular hydrogen gas for galaxies within massive halos, i.e., dense environments. In this Letter, we study how well TNG50 reproduces the resolved effects of a Virgo-like cluster environment on the gas surface densities of satellite galaxies with m * > 10 9 M ⊙ and star formation rate > 0.05 M ⊙ yr −1 . We select galaxies in the simulation that are analogous to those in the HERACLES and VERTICO surveys and mock-observe them to the common specifications of the data. Although TNG50 does not quantitatively match the observed gas surface densities in the centers of galaxies, the simulation does qualitatively reproduce the trends of gas truncation and central density suppression seen in VERTICO in both H i and H 2 . This result promises that modern cosmological hydrodynamic simulations can be used to reliably model the post-infall histories of cluster satellite galaxies.

We present CluSTAR-ND, a fast hierarchical galaxy/(sub)halo finder that produces Clustering Structure via Transformative Aggregation and Rejection in N-Dimensions. It is designed to improve upon Halo-OPTICS – an algorithm that automatically detects and extracts significant astrophysical clusters from the 3D spatial positions of simulation particles – by decreasing run-times, possessing the capability for metric adaptivity, and being readily applicable to data with any number of features. We directly compare these algorithms and find that not only does CluSTAR-ND produce a similarly robust clustering structure, it does so in a run-time that is at least 3 orders of magnitude faster. In optimising CluSTAR-ND's clustering performance, we have also carefully calibrated 4 of the 7 CluSTAR-ND parameters which – unless specified by the user – will be automatically and optimally chosen based on the input data. We conclude that CluSTAR-ND is a robust astrophysical clustering algorithm that can be leveraged to find stellar satellite groups on large synthetic or observational data sets.

Cosmological simulations are useful tools for studying the evolution of galaxies, and it is critical to accurately identify galaxies and their halos from raw simulation data. The friends-of-friends (FoF) algorithm has been widely adopted for this purpose because of its simplicity and expandability to higher dimensions. However, it is cost-inefficient when applied to high-resolution simulations because standard FoF implementation leads to too many distance calculations in dense regions. We confirm this through our exercise of applying the six-dimensional (6D) FoF galaxy finder code, VELOCI raptor , on the N ew H orizon simulation. The high particle resolution of N ew H orizon ( M star ∼ 10 4 M ⊙ ) allows a large central number density (10 6 kpc −3 ) for typical galaxies, resulting in a few days to weeks of galaxy searches for just one snapshot. Even worse, we observed a significant decrease in the FoF performance in the high-dimensional 6D searches: “the curse of dimensionality” problem. To overcome these issues, we have developed several implementations that can be readily applied to any tree-based FoF code. They include limiting visits to tree nodes, reordering the list of particles for searching neighbor particles, and altering the tree structure. Compared to the run with the original code, the new run with these implementations results in the identical galaxy detection with the ideal performance, O ( N log N ) , N being the number of particles in a galaxy—with a speed gain of a factor of 2700 in 3D or 12 in a 6D FoF search.

Using high-resolution Hubble Space Telescope imaging data, we perform a visual morphological classification of ∼ 36, 000 galaxies at z < 1 in the DEVILS/COSMOS region. As the main goal of this study, we derive the stellar mass function (SMF) and stellar mass density (SMD) sub-divided by morphological types. We find that visual morphological classification using optical imaging is increasingly difficult at z > 1 as the fraction of irregular galaxies and merger systems (when observed at rest-frame UV/blue wavelengths) dramatically increases. We determine that roughly two-thirds of the total stellar mass of the Universe today was in place by z ∼ 1. Double-component galaxies dominate the SMD at all epochs and increase in their contribution to the stellar mass budget to the present day. Elliptical galaxies are the second most dominant morphological type and increase their SMD by ∼ 2.5 times, while by contrast, the pure-disk population significantly decreases by ∼ 85%. According to the evolution of both highand low-mass ends of the SMF, we find that mergers and in-situ evolution in disks are both present at z < 1, and conclude that double-component galaxies are predominantly being built by the in-situ evolution in disks (apparent as the growth of the low-mass end with time), while mergers are likely responsible for the growth of ellipticals (apparent as the increase of intermediate/high-mass end).

Predicting the merger timescale ($\tau_{\rm merge}$) of merging dark matter halos, based on their orbital parameters and the structural properties of their hosts, is a fundamental problem in gravitational dynamics that has important consequences for our understanding of cosmological structure formation and galaxy formation. Previous models predicting $\tau_{\rm merge}$ have shown varying degrees of success when compared to the results of cosmological $N$-body simulations. We build on this previous work and propose a new model for $\tau_{\rm merge}$ that draws on insights derived from these simulations. We find that published predictions can provide reasonable estimates for $\tau_{\rm merge}$ based on orbital properties at infall, but tend to underpredict $\tau_{\rm merge}$ inside the host virial radius ($R_{200}$) because tidal stripping is neglected, and overpredict it outside $R_{200}$ because the host mass is underestimated. Furthermore, we find that models that account for orbital angular momentum via the circular radius $R_{\rm circ}$ underpredict (overpredict) $\tau_{\rm merge}$ for bound (unbound) systems. By fitting for the dependence of $\tau_{\rm merge}$ on various orbital and host halo properties,we derive an improved model for $\tau_{\rm merge}$ that can be applied to a merging halo at any point in its orbit. Finally, we discuss briefly the implications of our new model for $\tau_{\rm merge}$ for semi-analytical galaxy formation modelling.

We build upon Ordering Points To Identify the Clustering Structure (optics ), a hierarchical clustering algorithm well known to be a robust data miner, in order to produce halo-optics , an algorithm designed for the automatic detection and extraction of all meaningful clusters between any two arbitrary sizes. We then apply halo-optics to the 3D spatial positions of halo particles within four separate synthetic Milky Way-type galaxies, classifying the stellar and dark matter structural hierarchies. Through visualization of the halo-optics output, we compare its structure identification to the state-of-the-art galaxy/(sub)halo finder VELOCIraptor , finding excellent agreement even though halo-optics does not consider kinematic information in this current implementation. We conclude that halo-optics is a robust hierarchical halo finder, although its determination of lower spatial-density features such as the tails of streams could be improved with the inclusion of extra localized information such as particle kinematics and stellar metallicity into its distance metric.

Resolving faint galaxies in large volumes is critical for accurate cosmic reionisation simulations. While less demanding than hydrodynamical simulations, semi-analytic reionisation models still require very large N-body simulations in order to resolve the atomic cooling limit across the whole reionisation history within box sizes $\gtrsim 100 \, h^{-1} {\rm Mpc}$. To facilitate this, we extend the mass resolution of N-body simulations using a Monte Carlo algorithm. We also propose a method to evolve positions of Monte Carlo halos, which can be an input for semi-analytic reionisation models. To illustrate, we present an extended halo catalogue that reaches a mass resolution of $M_\text{halo} = 3.2 \times 10^7 \, h^{-1} \text{M}_\odot$ in a $105 \, h^{-1} {\rm Mpc}$ box, equivalent to an N-body simulation with $\sim 6800^3$ particles. The resulting halo mass function agrees with smaller volume N-body simulations with higher resolution. Our results also produce consistent two-point correlation functions with analytic halo bias predictions. The extended halo catalogues are applied to the \textsc{meraxes} semi-analytic reionisation model, which improves the predictions on stellar mass functions, star formation rate densities and volume-weighted neutral fractions. Comparison of high resolution large volume simulations with both small volume or low resolution simulations confirms that both low resolution and small volume simulations lead to reionisation ending too rapidly. Lingering discrepancies between the star formation rate functions predicted with and without our extensions can be traced to the uncertain contribution of satellite galaxies.

We present the Virgo Environment Traced in CO (VERTICO) survey, a new effort to map 12CO (2–1), 13CO (2–1), and C18O (2–1) in 51 Virgo Cluster galaxies with the Atacama Compact Array, part of the Atacama Large Millimeter/submillimeter Array. The primary motivation of VERTICO is to understand the physical mechanisms that perturb molecular gas disks, and therefore star formation and galaxy evolution, in dense environments. This first paper contains an overview of VERTICO's design and sample selection, 12CO (2–1) observations, and data reduction procedures. We characterize global 12CO (2–1) fluxes and molecular gas masses for the 49 detected VERTICO galaxies, provide upper limits for the two nondetections, and produce resolved 12CO (2–1) data products (median resolution = 8″ ≈ 640 pc). Azimuthally averaged 12CO (2–1) radial intensity profiles are presented along with derived molecular gas radii. We demonstrate the scientific power of VERTICO by comparing the molecular gas size–mass scaling relation for our galaxies with a control sample of field galaxies, highlighting the strong effect that radius definition has on this correlation. We discuss the drivers of the form and scatter in the size–mass relation and highlight areas for future work. VERTICO is an ideal resource for studying the fate of molecular gas in cluster galaxies and the physics of environment-driven processes that perturb the star formation cycle. Upon public release, the survey will provide a homogeneous legacy data set for studying galaxy evolution in our closest cluster.

For astronomers to make a significant contribution to the reduction of climate change-inducing greenhouse gas emissions, we first must quantify the sources of our emissions and review the most effective approaches for reducing them. Here we estimate that Australian astronomers’ total greenhouse gas emissions from their regular work activities are ≳25 ktCO2e yr–1 (equivalent kilotonnes of carbon dioxide per year). This can be broken into ~15 ktCO2e yr–1 from supercomputer usage, ~4.2 ktCO2e yr–1 from flights (where individuals’ flight emissions correlate with seniority), >3.3 ktCO2e yr–1 from the operation of observatories, and 2.6 ± 0.4 ktCO2e yr–1 from powering office buildings. Split across faculty scientists, postdoctoral researchers and PhD students, this averages to ≳37 tCO2e yr–1 per astronomer, more than 40% greater than the average Australian non-dependant’s emissions in total, and equivalent to around five times the global average. To combat these environmentally unsustainable practices, we suggest that astronomers should strongly preference the use of supercomputers, observatories and office spaces that are predominantly powered by renewable energy sources. Where current facilities do not meet this requirement, their funders should be lobbied to invest in renewables, such as solar or wind farms. Air travel should also be reduced wherever possible, replaced primarily by video conferencing, which should also promote inclusivity. Australian astronomers generate more greenhouse emissions than the average Australian citizen, thereby exacerbating the climate crisis. By quantifying contributions from different activities such as supercomputing and air travel, as presented here, astronomers can focus on reducing emissions by changing their practices in the most critical areas.

Galaxy cluster outskirts mark the transition region from the mildly non-linear cosmic web to the highly non-linear, virialised, cluster interior. It is in this transition region that the intra-cluster medium (ICM) begins to influence the properties of accreting galaxies and groups, as ram pressure impacts a galaxy’s cold gas content and subsequent star formation rate. Conversely, the thermodynamical properties of the ICM in this transition region should also feel the influence of accreting substructure (i.e. galaxies and groups), whose passage can drive shocks. In this paper, we use a suite of cosmological hydrodynamical zoom simulations of a single galaxy cluster, drawn from the nIFTy comparison project, to study how the dynamics of substructure accreted from the cosmic web influences the thermodynamical properties of the ICM in the cluster’s outskirts. We demonstrate how features evident in radial profiles of the ICM (e.g. gas density and temperature) can be linked to strong shocks, transient and short-lived in nature, driven by the passage of substructure. The range of astrophysical codes and galaxy formation models in our comparison are broadly consistent in their predictions (e.g. agreeing when and where shocks occur, but differing in how strong shocks will be); this is as we would expect of a process driven by large-scale gravitational dynamics and strong, inefficently radiating, shocks. This suggests that mapping such shock structures in the ICM in a cluster’s outskirts (via e.g. radio synchrotron emission) could provide a complementary measure of its recent merger and accretion history.

Upcoming measurements of the highly redshifted 21 cm line with next-generation radio telescopes such as the Hydrogen Epoch of Reionization Array (HERA) and Square Kilometer Array will provide the intriguing opportunity to probe dark matter (DM) physics during the Epoch of Reionization, Cosmic Dawn, and the Dark Ages. With HERA already under construction, there is a pressing need to thoroughly understand the impact of DM physics on the intergalactic medium (IGM) during these epochs. We present first results of a hydrodynamic simulation suite with 2 x 512(3) particles in a (100 h(-1) Mpc)(3) box with DM annihilation and baryonic cooling physics. We focus on redshift z similar to 11, just before reionization starts in our simulations, and discuss the imprint of DM annihilation on the IGM and on structure formation. We find that whereas structure formation is not affected by thermal weakly interacting massive particles heavier than m(chi) greater than or similar to 100 MeV, heating from O(GeV) DM particles may leave a significant imprint on the IGM that alters the 21 cm signal. Cold gas in low-density regions is particularly susceptible to the effects of DM heating. We note, however, that delayed energy deposition is not currently accounted for in our simulations.

We study the Intra-Halo Stellar Component (IHSC) of Milky Way-mass systems up to galaxy clusters in the Horizon-AGN cosmological hydrodynamical simulation. We identify the IHSC using an improved phase-space galaxy finder algorithm which provides an adaptive, physically motivated and shape-independent definition of this stellar component, that can be applied to halos of arbitrary masses. We explore the IHSC mass fraction-total halo's stellar mass, $f_{M*,IHSC}-M*$, relation and the physical drivers of its scatter. We find that on average the $f_{M*,IHSC}$ increases with $M_{*,tot}$, with the scatter decreasing strongly with mass from 2 dex at $M_{*,tot}\sim10^{11}M_\odot$ to 0.3 dex at group masses. At high masses, $M_{*,tot}>10^{11.5}M_\odot$, $f_{M*,IHSC}$ increases with the number of substructures, and with the mass ratio between the central galaxy and largest satellite, at fixed $M_{*,tot}$. From mid-size groups and systems below $M_{*,tot}<10^{12}M_\odot$, we find that the central galaxy's stellar rotation-to-dispersion velocity ratio, V/{\sigma}, displays the strongest (anti)-correlation with $f_{M*,IHSC}$ at fixed $M_{*,tot}$ of all the galaxy and halo properties explored, transitioning from $f_{M*,IHSC}$<0.1% for high V/{\sigma}, to $f_{M*,IHSC}\sim5$% for low V/{\sigma} galaxies. By studying the $f_{M*,IHSC}$ temporal evolution, we find that, in the former, mergers not always take place, but if they did, they happened early (z>1), while the high $f_{M*,IHSC}$ population displays a much more active merger history. In the case of massive groups and galaxy clusters, $M_{*,tot}>10^{12}M_\odot$, a fraction $f_{M*,IHSC}\sim$10-20% is reached at $z\sim1$ and then they evolve across lines of constant $f_{M*,IHSC}$ modulo some small perturbations. Because of the limited simulation's volume, the latter is only tentative and requires a larger sample of simulated galaxy clusters to confirm.

Hierarchical models of structure formation predict that dark matter halo assembly histories are characterised by episodic mergers and interactions with other haloes. An accurate description of this process will provide insights into the dynamical evolution of haloes and the galaxies that reside in them. Using large cosmological N-body simulations, we characterise halo orbits to study the interactions between substructure haloes and their hosts, and how different evolutionary histories map to different classes of orbits. We use two new software tools - WhereWolf, which uses halo group catalogues and merger trees to ensure that haloes are tracked accurately in dense environments, and OrbWeaver, which quantifies each halo's orbital parameters. We demonstrate how WhereWolf improves the accuracy of halo merger trees, and we use OrbWeaver to quantify orbits of haloes. We assess how well analytical prescriptions for the merger timescale from the literature compare to measured merger timescales from our simulations and find that existing prescriptions perform well, provided the ratio of substructure-to-host mass is not too small. In the limit of small substructure-to-host mass ratio, we find that the prescriptions can overestimate the merger timescales substantially, such that haloes are predicted to survive well beyond the end of the simulation. This work highlights the need for a revised analytical prescription for the merger timescale that more accurately accounts for processes such as catastrophic tidal disruption.

ABSTRACT The standard cosmological paradigm currently lacks a detailed account of physics in the dark sector, the dark matter and energy that dominate cosmic evolution. In this paper, we consider the distinguishing factors between three alternative models – warm dark matter, quintessence, and coupled dark matter–energy – and lambda cold dark matter (ΛCDM) through numerical simulations of cosmological structure formation. Key halo statistics – halo spin/velocity alignment between large-scale structure and neighbouring haloes, halo formation time, and migration – were compared across cosmologies within the redshift range 0 ≤ z ≤ 2.98. We found the alignment of halo motion and spin to large-scale structures and neighbouring haloes to be similar in all cosmologies for a range of redshifts. The search was extended to low-density regions, avoiding non-linear disturbances of halo spins, yet very similar alignment trends were found between cosmologies, which are difficult to characterize and use as a probe of cosmology. We found that haloes in quintessence cosmologies form earlier than their ΛCDM counterparts. Relating this to the fact that such haloes originate in high-density regions, such findings could hold clues to distinguishing factors for the quintessence cosmology from the standard model. However, in general, halo statistics are not an accurate probe of the dark sector physics.

ABSTRACT We model the large-scale linear galaxy bias bg(x, z) as a function of redshift z and observed absolute magnitude threshold x for broad-band continuum emission from the far-infrared to ultraviolet, as well as for prominent emission lines, such as the H α, H β, Ly a, and [O ii] lines. The modelling relies on the semi-analytic galaxy formation model galform, run on the state-of-the-art N-body simulation surfs with the Planck 2015 cosmology. We find that both the differential bias at observed absolute magnitude x and the cumulative bias for magnitudes brighter than x can be fitted with a five-parameter model: bg(x, z) = a + b(1 + z)e(1 + exp [(x − c)d]). We also find that the bias for the continuum bands follows a very similar form regardless of wavelength due to the mixing of star-forming and quiescent galaxies in a magnitude-limited survey. Differences in bias only become apparent when an additional colour separation is included, which suggest extensions to this work could look at different colours at fixed magnitude limits. We test our fitting formula against observations, finding reasonable agreement with some measurements within 1σ statistical uncertainties, and highlighting areas of improvement. We provide the fitting parameters for various continuum bands, emission lines, and intrinsic galaxy properties, enabling a quick estimation of the linear bias in any typical survey of large-scale structure.

ABSTRACT Linking the properties of galaxies to the assembly history of their dark matter haloes is a central aim of galaxy evolution theory. This paper introduces a dimensionless parameter s ∈ [0, 1], the ‘tree entropy’, to parametrize the geometry of a halo’s entire mass assembly hierarchy, building on a generalization of Shannon’s information entropy. By construction, the minimum entropy (s = 0) corresponds to smoothly assembled haloes without any mergers. In contrast, the highest entropy (s = 1) represents haloes grown purely by equal-mass binary mergers. Using simulated merger trees extracted from the cosmological N-body simulation SURFS, we compute the natural distribution of s, a skewed bell curve peaking near s = 0.4. This distribution exhibits weak dependences on halo mass M and redshift z, which can be reduced to a single dependence on the relative peak height δc/σ(M, z) in the matter perturbation field. By exploring the correlations between s and global galaxy properties generated by the SHARK semi-analytic model, we find that s contains a significant amount of information on the morphology of galaxies – in fact more information than the spin, concentration, and assembly time of the halo. Therefore, the tree entropy provides an information-rich link between galaxies and their dark matter haloes.

ABSTRACT We present the first detection of mass-dependent galactic spin alignments with local cosmic filaments with >2σ confidence using IFS kinematics. The 3D network of cosmic filaments is reconstructed on Mpc scales across GAlaxy and Mass Assembly fields using the cosmic web extractor DisPerSe. We assign field galaxies from the SAMI survey to their nearest filament segment in 3D and estimate the degree of alignment between SAMI galaxies’ kinematic spin axis and their nearest filament in projection. Low-mass galaxies align their spin with their nearest filament while higher mass counterparts are more likely to display an orthogonal orientation. The stellar transition mass from the first trend to the second is bracketed between $10^{10.4}$ and $10^{10.9}\, \mathrm{ M}_{\odot }$, with hints of an increase with filament scale. Consistent signals are found in the Horizon-AGN cosmological hydrodynamic simulation. This supports a scenario of early angular momentum build-up in vorticity rich quadrants around filaments at low stellar mass followed by progressive flip of spins orthogonal to the cosmic filaments through mergers at high stellar mass. Conversely, we show that dark matter only simulations post-processed with a semi-analytical model treatment of galaxy formation struggles to reproduce this alignment signal. This suggests that gas physics is key in enhancing the galaxy-filament alignment.

We present new cosmological hydrodynamic simulations that incorporate dark matter annihilation feedback (DMAF), whereby energy released from the annihilation of dark matter particles through decay channels such as photon or positron-electron pairs provides additional heating sources for local baryonic material. For annihilation rates comparable to weakly interacting massive particle (WIMP)-like particles, we find that the key influence of DMAF is to inhibit gas accretion on to haloes. Such diminished gas accretion early in the lifetimes of haloes results in reduced gas fractions in smaller haloes, and the delayed halo formation times of larger structures, suggesting that DMAF could impact the stellar age distribution in galaxies and morphology of dwarfs. For a dark matter particle mass of m(chi) similar to 10 MeV, there is a 'critical halo mass' of similar to 10(13) M-circle dot at z = 0, below which there are large differences when compared to Lambda cold dark matter, such as a reduction in the abundance of halo structures as large as 25 per cent, reduced gas content by 50 per cent and central gas densities reduced down to 10 per cent within haloes of mass similar to 10(12) M-circle dot but with increasing effects in smaller haloes. Higher dark matter particle mass models have a smaller 'critical halo mass'. For a m(chi) similar to 100 MeV model, we find differences start appearing below halo masses of similar to 10(12) M-circle dot and a m(chi) greater than or similar to 1 GeV model, this mass scale lies below the resolution of our simulations, though we still observe changes in the morphology of dwarf galaxies.