Merger History Research Articles

The Imprint of Dark Matter on the Galactic Acceleration Field

Measurements of the accelerations of stars enabled by time-series extreme-precision spectroscopic observations, pulsar timing, and eclipsing binary stars in the solar neighborhood offer insights into the mass distribution of the Milky Way that do not rely on traditional equilibrium modeling. Given the measured accelerations, we can determine a total mass density and infer the amount of dark matter (DM) by accounting for the mass in stars, gas, and dust. Leveraging FIRE-2 simulations of Milky Way–mass galaxies we compare vertical acceleration profiles between cold DM (CDM) and self-interacting DM (SIDM) with a constant cross section of 1 cm2 g−1 across three halos with diverse assembly histories. Notably, significant asymmetries in vertical acceleration profiles near the midplane at fixed radii are observed in both CDM and SIDM, particularly in halos recently affected by mergers with satellites of Sagittarius/SMC-like masses or greater. These asymmetries offer a unique window into exploring the merger history of a galaxy. We show that SIDM halos manifest a more oblate shape and consistently exhibit higher local stellar and DM densities and steeper vertical acceleration gradients, up to 10%–30% steeper near the solar neighborhood. However, similar magnitude changes can arise from azimuthal variations in the baryonic components at a fixed radius and external influences like mergers, making it difficult to distinguish between CDM and SIDM using acceleration measurements in a single galaxy.

Action-based dynamical models of M31-like galaxies

ABSTRACT In this work, we present an action-based dynamical equilibrium model to constrain the phase-space distribution of stars in the stellar halo, present-day dark matter distribution, and the total mass distribution in M31-like galaxies. The model comprises a three-component gravitational potential (stellar bulge, stellar disc, and a dark matter halo), and a double power-law distribution function (DF), $f(\mathbf {J})$, which is a function of actions. A Bayesian model-fitting algorithm was implemented that enabled both parameters of the potential and DF to be explored. After testing the model-fitting algorithm on mock data drawn from the model itself, it was applied to a set of three M31-like haloes from the Auriga simulations (Auriga 21, Auriga 23, and Auriga 24). Furthermore, we tested the equilibrium assumption and the ability of a double power-law DF to represent the stellar halo stars. The model incurs an error in the total enclosed mass of around 10 per cent out to 100 kpc, thus justifying the equilibrium assumption. Furthermore, the double power-law DF used proves to be an appropriate description of the investigated M31-like haloes. The anisotropy profiles of the haloes were also investigated and discussed from a merger history point of view, however, our approach underscores to capture the non-equilibrium anisotropy information of stellar haloes in Auriga models, which is expected.

Galaxy archaeology for wet mergers: Globular cluster age distributions in the Milky Way and nearby galaxies

Context. Identifying past wet merger activity in galaxies has been a longstanding issue in extragalactic formation history studies. Gaia’s 6D kinematic measurements in our Milky Way (MW) have vastly extended the possibilities for Galactic archaeology, leading to the discovery of a multitude of early mergers in the MW’s past. As recent work has established a link between younger globular clusters (GCs; less than about 10–11 Gyr old) and wet galaxy merger events, the MW provides an ideal laboratory for testing which GC properties can be used to trace extragalactic galaxy formation histories. Aims. To test the hypothesis that GCs trace wet mergers, we relate the measured GC age distributions of the MW and three nearby galaxies, M 31, NGC 1407, and NGC 3115, to their merger histories and interpret the connection with wet mergers through an empirical model for GC formation. Methods. The GC ages of observed galaxies are taken from a variety of studies to analyze their age distributions side-by-side with the model. For the MW, we additionally cross-match the GCs with their associated progenitor host galaxies to disentangle the connection to the GC age distribution. For the modeled GCs, we take galaxies with similar GC age distributions as observed to compare their accretion histories with those inferred through observations. Results. We find that the MW GC age distribution is bimodal, mainly caused by younger GCs (10–11 Gyr old associated with Gaia-Sausage/Enceladus (GSE) and in part by unassociated high-energy GCs. The GSE GC age distribution also appears to be bimodal. We propose that the older GSE GCs (12–13 Gyr old) were accreted together with GSE, while the younger ones formed as a result of the merger. For the nearby galaxies, we find that clear peaks in the GC age distributions coincide with active early gas-rich merger phases. Even small signatures in the GC age distributions agree well with the expected wet formation histories of the galaxies inferred through other observed tracers. From the models, we predict that the involved cold gas mass can be estimated from the number of GCs found in the formation burst. Conclusions. Multimodal GC age distributions can trace massive wet mergers as a result of GCs being formed through them. From the laboratory of our own MW and nearby galaxies we conclude that the ages of younger GC populations of galaxies can be used to infer the wet merger history of a galaxy.

Formation of Superthin Galaxies in IllustrisTNG

Superthin galaxies are observed to have stellar disks with extremely small minor-to-major axis ratios. In this work, we investigate the formation of superthin galaxies in the TNG100 simulation. We trace the merger history and investigate the evolution of galaxy properties of a selected sample of superthin galaxies and a control sample of galaxies that share the same joint probability distribution in the stellar-mass and color diagram. Through making comparisons between the two galaxy samples, we find that present-day superthin galaxies had similar morphologies as the control sample counterparts at higher redshifts, but have developed extended flat “superthin” morphologies since z ∼ 1. During this latter evolution stage, superthin galaxies undergo an overwhelmingly higher frequency of prograde mergers (with orbit-spin angle θ orb ≤ 40°). Accordingly the spins of their dark matter halos have grown significantly and become noticeably higher than those of their normal disk counterparts. This further results in the buildup of their stellar disks at larger distances much beyond the regimes of normal disk galaxies. We also discuss the formation scenario of those superthin galaxies that live in larger dark matter halos as satellite galaxies therein.

Effect of Wave Dark Matter on Equal Mass Black Hole Mergers

For dark matter to be detectable with gravitational waves from binary black holes, it must reach higher than average densities in their vicinity. In the case of light (wavelike) dark matter, the density of dark matter between the binary can be significantly enhanced by accretion from the surrounding environment. Here we show that the resulting dephasing effect on the last ten orbits of an equal mass binary is maximized when the Compton wavelength of the scalar particle is comparable to the orbital separation, 2π/μ∼d. The phenomenology of the effect is different from the channels that are usually discussed, where dynamical friction (along the orbital path) and radiation of energy and angular momentum drive the dephasing, and is rather dominated by the radial force (the spacetime curvature in the radial direction) towards the overdensity between the black holes. While our numerical studies limit us to scales of the same order, this effect may persist at larger separations and/or particle masses, playing a significant role in the merger history of binaries. Published by the American Physical Society 2024

Detecting Wandering Intermediate-Mass Black Holes with AXIS in the Milky Way and Local Massive Galaxies

This white paper explores the detectability of intermediate-mass black holes (IMBHs) wandering in the Milky Way (MW) and massive local galaxies, with a particular emphasis on the role of AXIS. IMBHs, ranging within 103−6M⊙, are commonly found at the centers of dwarf galaxies and may exist, yet undiscovered, in the MW. By using model spectra for advection-dominated accretion flows (ADAFs), we calculated the expected fluxes emitted by a population of wandering IMBHs with masses of 105M⊙ in various MW environments and extrapolated our results to massive local galaxies. Around 40% of the potential population of wandering IMBHs in the MW can be detected in an AXIS deep field. We proposed criteria to aid with selecting IMBH candidates using already available optical surveys. We also showed that IMBHs wandering in >200 galaxies within 10 Mpc can be easily detected with AXIS when passing within dense galactic environments (e.g., molecular clouds and cold neutral medium). In summary, we highlighted the potential X-ray detectability of wandering IMBHs in local galaxies and provided insights for guiding future surveys. Detecting wandering IMBHs is crucial for understanding their demographics and evolution and the merging history of galaxies. This white paper is part of a series commissioned for the AXIS Probe Concept Mission; additional AXIS white papers can be found at the AXIS website.

Spatially resolved spectroscopic observations of gas emission in dwarf galaxies hosting accreting black hole candidate

ABSTRACT Recent studies on dwarf galaxies reveal that some of them harbour a massive black hole (BH), which is believed to have a similar mass of the supermassive BH ‘seeds’ at early times. The origin and growth of the primitive BHs are still open questions, since these BH seeds are hardly observed at high redshifts. Therefore, MBH of dwarf galaxies can be the perfect candidates to untangle BH ‘seeds’ properties and their influence on their host galaxy evolution, since MBH may preserve their initial conditions due to its quiet merger and accretion histories. We use optical integral field unit observations, obtained with the Gemini GMOS-IFU, to study the gas emission and kinematics in four dwarf galaxies, candidates to host MBH, based on the analysis of their [Fe x] luminosities measured from SDSS spectra. The [Fe x] emission line is not detected in our GMOS in any of the galaxies, prompting speculation that its absence in our recent data may stem from a past tidal disruption event coinciding with the observation period of the SDSS data. All galaxies exhibit extended gas emissions, and the spatially resolved emission-line ratio diagnostic diagrams present values that suggest active galactic nuclei (AGN) photoionization from the [S ii]–BPT diagram. The gas velocity fields of all galaxies are indicative of disturbed rotation patterns, with no detection of gas outflows in any of the sources. Although the [S ii]–BPT diagrams indicate AGN photoionization, further confirmation through multiwavelength observations is required to validate this scenario.

Characterizing ultra-high-redshift dark matter halo demographics and assembly histories with the gureft simulations

ABSTRACT Dark matter halo demographics and assembly histories are a manifestation of cosmological structure formation and have profound implications for the formation and evolution of galaxies. In particular, merger trees provide fundamental input for several modelling techniques, such as semi-analytic models (SAMs), sub-halo abundance matching (SHAM), and decorated halo occupation distribution models. Motivated by the new ultra-high-redshift (z ≳ 10) frontier enabled by JWST, we present a new suite of Gadget at Ultrahigh Redshift with Extra-Fine Timesteps (gureft) dark matter-only cosmological simulations that are carefully designed to capture halo merger histories and structural properties in the ultra-z universe. The simulation suite consists of four 10243-particle simulations with box sizes of 5, 15, 35, and 90 Mpc h−1, each with 170 snapshots stored between 40 ≥ z ≥ 6. With the unprecedented number of available snapshots and strategically chosen dynamic range covered by these boxes, gureft uncovers the emerging dark matter halo populations and their assembly histories in the earliest epochs of cosmic history. In this work, we present the halo mass functions (HMF) between z ∼ 20 and 6 down to log (Mvir/M⊙) ∼ 5, and show that at high redshift, these robust HMFs can differ substantially from commonly used analytic approximations or older fitting functions in the literature. We also present key physical properties of the ultra-high z halo population, such as concentration and spin, as well as their mass growth and merger rates, and again provide updated fitting functions.

Mutual information between galaxy properties and the initial predisposition

The immense diversity of the galaxy population in the universe is believed to stem from their disparate merging and star formation histories, and multi-scaleinfluences of diverse environments. No single causal factor of the initial state is known to explain how the galaxies formed and evolved to end uppossessing such various traits as they have at the present epoch. However, several observational studies have revealed that the key physical properties of the observedgalaxies in the local universe appeared to have a much simpler, lower-dimensional correlation structure than expected, the origin of which remains unexplained.Speculating that the emergence of such a simple correlation structure of the galaxy properties must be triggered by nature rather than by nurture,we explore if the present galaxy properties may be correlated with the initial precondition for protogalaxy angular momentum, τ, and test it against the datafrom the IllustrisTNG300-1 hydrodynamic simulation. Employing Shannon's information theory, we discover that τ shares a significantly large amount of mutual informationwith each of the four basic traits of the TNG galaxies at z = 0: the spin parameters, formation epochs, stellar-to-total mass ratios, and fraction of kinetic energy in ordered rotation.These basic traits except for the stellar-to-total mass ratios are found to contain even a larger amount of MI about τ than about the total masses and environmentsfor the case of giant galaxies with 11.5 ≤ log[M t/(h -1 M ⊙)] < 13. Our results imply that the initial condition of the universe must be more impactful on the galaxyevolution than conventionally thought.

Constraints on peculiar velocity distribution of binary black holes using gravitational waves with GWTC-3

Peculiar velocity encodes rich information about the formation, dynamics, evolution, and merging history of binary black holes. In this work, we employ a hierarchical Bayesian model to infer the peculiar velocity distribution of binary black holes. We use the data from GWTC-3 and assume a Maxwell-Boltzmann distribution for the peculiar velocities, but do not consider the dependence of peculiar velocity on the masses of black hole binaries. The constraint on the peculiar velocity distribution parameter, v 0, is weak and uninformative. However, the determination of peculiar velocity distribution can be significantly improved with next-generation ground-based gravitational wave detectors. For the Einstein Telescope, the relative uncertainty of v 0 will reduce to ∼ 10% using 103 golden binary black hole events. Our statistical approach thus provides a robust and prospective inference for determining the peculiar velocity distribution.

Fast forward modelling of galaxy spatial and statistical distributions

A forward modelling approach provides simple, fast and realistic simulations of galaxy surveys, without a complex underlying model. For this purpose, galaxy clustering needs to be simulated accurately, both for the usage of clustering as its own probe and to control systematics. We present a forward model to simulate galaxy surveys, where we extend the Ultra-Fast Image Generator to include galaxy clustering. We use the distribution functions of the galaxy properties, derived from a forward model adjusted to observations. This population model jointly describes the luminosity functions, sizes, ellipticities, SEDs and apparent magnitudes. To simulate the positions of galaxies, we then use a two-parameter relation between galaxies and halos with Subhalo Abundance Matching (SHAM). We simulate the halos and subhalos using the fast PINOCCHIO code, and a method to extract the surviving subhalos from the merger history. Our simulations contain a red and a blue galaxy population, for which we build a SHAM model based on star formation quenching. For central galaxies, mass quenching is controlled with the parameter Mlimit, with blue galaxies residing in smaller halos. For satellite galaxies, environmental quenching is implemented with the parameter tquench, where blue galaxies occupy only recently merged subhalos. We build and test our model by comparing to imaging data from the Dark Energy Survey Year 1. To ensure completeness in our simulations, we consider the brightest galaxies with i < 20. We find statistical agreement between our simulations and the data for two-point correlation functions on medium to large scales. Our model provides constraints on the two SHAM parameters Mlimit and tquench and offers great prospects for the quick generation of galaxy mock catalogues, optimized to agree with observations.

The rise and fall of bars in disc galaxies from z = 1 to z = 0

Context. Stellar bars are non-axisymmetric structures found in over 30 per cent of massive disc galaxies in the local Universe. The environment could play a significant role in determining whether or not a spiral galaxy is likely to develop a bar. Aims. We investigate the influence of the environment on the evolution of barred and unbarred disc galaxies with a mass of larger than 1010 M⊙ from z = 1 down to z = 0, employing the TNG50 magnetic-hydrodynamical simulation. Methods. We determined the fraction of barred galaxies that conserve their bar and the fraction of those that lost it by z = 0. We also estimate the fraction of unbarred galaxies at z = 1 that develop a bar at later times. We study the merger histories and the distance of close companions for each category to understand the role of the environment in the evolution of these galaxies. Results. We find that 49 per cent of z = 1 disc galaxies undergo a morphological transformation, transitioning into either a lenticular or spheroidal galaxy, while the other 51 per cent retain the large disc shape. The morphological alteration is mostly influenced by the environment. Lenticular and spheroidal galaxies tend to exist in denser environments and have more frequent mergers compared to disc galaxies. We find that bars are stable after they have formed, as over half of the barred galaxies (60.2 per cent) retain the bar structure and have experienced fewer mergers compared to those galaxies that lose their bars (5.6 per cent). These latter galaxies start with weaker and shorter bars at z = 1, are influenced by tidal interactions, and are frequently observed in more populated areas. Additionally, our study reveals that less than 20 per cent of unbarred galaxies will never develop a bar and exhibit the quietest merger history. Unbarred galaxies that undergo bar formation after z = 1 more frequently experience a merger event. Furthermore, tidal interactions with a close companion may account for bar formation in at least one-third of these instances. Conclusions. Our findings highlight that stable bars are prevalent in disc galaxies. Bar evolution may nonetheless be affected by the environment. Interactions with nearby companions or tidal forces caused by mergers have the capacity to disrupt the disc. This perturbance may materialise as the dissolution of the bar, the formation of a bar, or, in its most severe form, the complete destruction of the disc, resulting in morphological transformation. Bars that are weak and short at z = 1 and undergo major or minor mergers may eventually dissolve, whereas unbarred galaxies that enter crowded environments or experience a merger may develop a bar.

The high energy X-ray probe (HEX-P): constraining supermassive black hole growth with population spin measurements

Constraining the primary growth channel of supermassive black holes (SMBHs) remains one the most actively debated questions in the context of cosmological structure formation. Owing to the expected connection between SMBH spin parameter evolution and the accretion and merger history of individual black holes, population spin measurements offer a rare observational window into the cosmic growth of SMBHs. As of today, the most common method for estimating SMBH spin relies on modeling the relativistically broaden atomic profiles in the reflection spectrum observed in X-rays. In this paper, we study the observational requirements needed to confidently distinguish between the primary SMBH growth channels based on their distinct spin-mass distributions predicted by the Horizon-AGN cosmological simulation. Indoing so, we characterize outstanding limitations associated with the existing measurements and discuss the landscape of future observational campaigns which could be planned and executed with future X-ray observatories. We focus our attention on the High-Energy X-ray Probe (HEX-P), a proposed probe-class mission designed to serve the high-energy community in the 2030s.

Differences in the properties of disrupted and surviving satellites of Milky-Way-mass galaxies in relation to their host accretion histories

ABSTRACT From the chemodynamical properties of tidal debris in the Milky Way, it has been inferred that the dwarf satellites that have been disrupted had different chemical abundances from their present-day counterparts of similar mass that survive today, specifically, they had lower [Fe/H] and higher [Mg/Fe]. Here we use the Artemis simulations to study the relation between the chemical abundances of disrupted progenitors of MW-mass galaxies and their stellar mass, and the evolution of the stellar mass–metallicity relations (MZR) of this population with redshift. We find that these relations have significant scatter, which correlates with the accretion redshifts (zacc) of satellites, and with their cold gas fractions. We investigate the MZRs of dwarf populations accreted at different redshifts and find that they have similar slopes, and also similar with the slope of the MZR of the surviving population (≈0.32). However, the entire population of disrupted dwarfs displays a steeper MZR, with a slope of ≈0.48, which can be explained by the changes in the mass spectrum of accreted dwarf galaxies with redshift. We find strong relations between the (mass-weighted) 〈zacc〉 of the disrupted populations and their global chemical abundances (〈[Fe/H]〉 and 〈[Mg/Fe]〉), which suggests that chemical diagnostics of disrupted dwarfs can be used to infer the types of merger histories of their hosts. For the case of the MW, our simulations predict that the bulk of the disrupted population was accreted at 〈zacc〉 ≈ 2, in agreement with other findings. We also find that disrupted satellites form and evolve in denser environments, closer to their hosts, than their present-day counterparts.

The AGORA High-resolution Galaxy Simulations Comparison Project. V. Satellite Galaxy Populations in a Cosmological Zoom-in Simulation of a Milky Way–Mass Halo

We analyze and compare the satellite halo populations at z ∼ 2 in the high-resolution cosmological zoom-in simulations of a 1012 M ⊙ target halo (z = 0 mass) carried out on eight widely used astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Arepo-t, and Gizmo) for the AGORA High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near z = 2 for each code (hereafter “z ∼ 2”) at which the eight simulations are in the same stage in the target halo’s merger history. After identifying the matched pairs of halos between the CosmoRun simulations and the DMO simulations, we discover that each CosmoRun halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at z ∼ 2, the number of satellite galaxies is significantly fewer than that of dark matter halos in all participating AGORA simulations and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called “missing satellite problem” is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in AGORA and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution (≲100 proper pc at z = 2). We also compare other properties such as the stellar mass–halo mass relation and the mass–metallicity relation. Our work highlights the value of comparison studies such as AGORA, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms.

Representing low-mass black hole seeds in cosmological simulations: A new sub-grid stochastic seed model

ABSTRACT The physical origin of the seeds of supermassive black holes (SMBHs), with postulated initial masses ranging from ∼105 M⊙ to as low as ∼102 M⊙, is currently unknown. Most existing cosmological hydrodynamic simulations adopt very simple, ad hoc prescriptions for BH seeding and seed at unphysically high masses ∼105–106 M⊙. In this work, we introduce a novel sub-grid BH seeding model for cosmological simulations that is directly calibrated to high-resolution zoom simulations that explicitly resolve ∼103 M⊙ seeds forming within haloes with pristine, dense gas. We trace the BH growth along galaxy merger trees until their descendants reach masses of ∼104 or 105 M⊙. The results are used to build a new stochastic seeding model that directly seeds these descendants in lower resolution versions of our zoom region. Remarkably, we find that by seeding the descendants simply based on total galaxy mass, redshift and an environmental richness parameter, we can reproduce the results of the detailed gas-based seeding model. The baryonic properties of the host galaxies are well reproduced by the mass-based seeding criterion. The redshift-dependence of the mass-based criterion captures the combined influence of halo growth, dense gas formation, and metal enrichment on the formation of ∼103 M⊙ seeds. The environment-based seeding criterion seeds the descendants in rich environments with higher numbers of neighbouring galaxies. This accounts for the impact of unresolved merger dominated growth of BHs, which produces faster growth of descendants in richer environments with more extensive BH merger history. Our new seed model will be useful for representing a variety of low-mass seeding channels within next-generation larger volume uniform cosmological simulations.

Characterising the intra-cluster light in The Three Hundred simulations

We characterise the intra-cluster light (ICL) in ensembles of full-physics cluster simulations from THE THREE HUNDRED project, a suite of 324 hydrodynamical resimulations of cluster-sized halos. We identify the ICL as those stellar particles bound to the potential of the cluster itself, but not to any of its substructures, and separate the brightest cluster galaxy (BCG) by means of a fixed 50 kpc aperture. We find the total BCG+ICL mass to be in agreement with state-of-the-art observations of galaxy clusters. The ICL mass fraction of our clusters is between 30 and 50% of the total stellar mass within R500, while the BCG represents around 10%. We further find no trend of the ICL fraction with cluster halo mass, at least not in the range [0.2, 3] 1015 h−1 M⊙ considered here. For the dynamical state, characterised both by theoretical estimators and by the recent merging history of the cluster, there is a clear correlation, such that more relaxed clusters and those that have undergone fewer recent mergers have a higher ICL fraction. Finally, we investigate the possibility of using the ICL to explore the dark matter (DM) component of galaxy clusters. We compute the volumetric density profile for the DM and ICL components and show that, up to R500, the ratio between the two can be described by a power law. Working with the velocity dispersion profiles instead, we show that the ratio can be fit by a straight line. Providing the parameters of these fits, we show how the ICL can be used to infer DM properties.

Swarming in stellar streams: Unveiling the structure of the Jhelum stream with ant colony-inspired computation

The halo of the Milky Way galaxy hosts multiple dynamically coherent substructures known as stellar streams that are remnants of tidally disrupted orbiting systems such as globular clusters (GCs) and dwarf galaxies (DGs). A particular case is that of the Jhelum stream, which is known for its unusual and complex morphology. Using the available data from the Gaia DR3 catalog, we extracted a region on the sky that contains Jhelum, and fine-tuned this selection by enforcing limits on the magnitude and proper motion of the selected stars. We then applied the novel Locally Aligned Ant Technique (LAAT) on the position and proper motion space of stars belonging to the selected region to highlight the stars that are closely aligned with a local manifold in the data and the stars belonging to regions of high local density. We find that the overdensity representing the stream in proper motion space is composed of two components, and show the correspondence of these two signals to the previously reported narrow and broad spatial components of Jhelum. We then made use of the radial velocity measurements provided by the S5 survey and confirm, for the first time, a separation between the stars belonging to the two components in radial velocity. We show that the narrow and broad components have velocity dispersions of 4.84+1.23-0.79 km s-1 and 19.49+2.19-1.84 km s−1, and metallicity dispersions of 0.15+0.18-0.10 and 0.34+0.13-0.09, respectively. These measurements, as well as the given difference in component widths, could be explained with a probable scenario where Jhelum is the remnant of a GC embedded within a DG and where both were accreted onto the Milky Way during their infall. Although the properties of Jhelum could be explained with this merger scenario, other progenitors of the narrow component remain possible such as a nuclear star cluster or a DG. To rule these possibilities out, we would need more observational data of member stars of the stream. Our analysis shows that the internal structure of streams holds great information on their past formation history, and therefore provides further insight into the merger history of the Milky Way.

On the Significance of the Thick Disks of Disk Galaxies

Thick disks are a prevalent feature observed in numerous disk galaxies, including our own Milky Way. Their significance has been reported to vary widely, ranging from a few percent to 100% of the disk mass, depending on the galaxy and the measurement method. We use the NewHorizon simulation, which has high spatial and stellar mass resolutions, to investigate the issue of the thick-disk mass fraction. We also use the NewHorizon2 simulation, which was run on the same initial conditions, but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M * > 1010 M ⊙ in NewHorizon, the contribution of the thick disk was found to be 20% ± 11% in r-band luminosity or 35% ± 15% in mass to the overall galactic disk, which seems in agreement with observational data. The vertical profiles of 0, 22, and 5 galaxies are best fitted by 1, 2, or 3 sech2 components, respectively. The NewHorizon2 data show that the selection of thick-disk stars based on a single [α/Fe] cut is contaminated by stars of different kinematic properties, while missing the bulk of kinematically thick disk stars. Vertical luminosity profile fits recover the key properties of thick disks reasonably well. The majority of stars are born near the galactic midplane with high circularity and get heated with time via fluctuations in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties.

The SAMI Galaxy Survey: using tidal streams and shells to trace the dynamical evolution of massive galaxies

ABSTRACT Slow rotator galaxies are distinct amongst galaxy populations, with simulations suggesting that a mix of minor and major mergers are responsible for their formation. A promising path to resolve outstanding questions on the type of merger responsible, is by investigating deep imaging of massive galaxies for signs of potential merger remnants. We utilize deep imaging from the Subaru-Hyper Suprime Cam Wide data to search for tidal features in massive [log10(M*/M⊙) &amp;gt; 10] early-type galaxies (ETGs) in the SAMI Galaxy Survey. We perform a visual check for tidal features on images where the galaxy has been subtracted using a Multi-Gauss Expansion (MGE) model. We find that 31$^{+2}_{-2}$ per cent of our sample show tidal features. When comparing galaxies with and without features, we find that the distributions in stellar mass, light-weighted mean stellar population age, and H${\alpha}$ equivalent width are significantly different, whereas spin ($\lambda _{R_{\rm {e}}}$), ellipticity, and bulge-to-total ratio have similar distributions. When splitting our sample in age, we find that galaxies below the median age (10.8 Gyr) show a correlation between the presence of shells and lower $\lambda _{R_{\rm {e}}}$, as expected from simulations. We also find these younger galaxies which are classified as having ‘strong’ shells have lower $\lambda _{R_{\rm {e}}}$. However, simulations suggest that merger features become undetectable within ∼2–4 Gyr post-merger. This implies that the relationship between tidal features and merger history disappears for galaxies with older stellar ages, i.e. those that are more likely to have merged long ago.