Stellar Component Research Articles

The distribution of globular clusters in kinematic spaces does not trace the accretion history of the host galaxy

Context. Reconstructing how all the stellar components of the Galaxy formed and assembled over time by studying the properties of the stars that form it is the aim of Galactic archaeology. Thanks to the launch of the ESA Gaia astrometric mission and the development of many spectroscopic surveys in recent years, we are for the first time in the position to delve into the layers of the past of the Galaxy. Globular clusters play a fundamental role in this research field since they are among the oldest stellar systems in the MW and thus bear witness to its entire past. Aims. As a natural result of galaxy formation, globular clusters did not necessarily all form in the Galaxy itself. Indeed, a fraction of them could have been formed in satellite galaxies accreted by the Milky Way over time. In recent years, there have been several attempts to constrain the nature of clusters (accreted or formed in the Milky Way itself) through the analysis of kinematic spaces, such as the E − Lz, Lperp − Lz, eccentricity − Lz, and the action space, as well as attempts to reconstruct the properties of the accretion events experienced by the Milky Way through time from this kind of analysis. This work aims to test a widely used assumption about the clustering of the accreted populations of globular clusters in the integrals of motions space. Methods. In this paper we analyse a set of dissipationless N-body simulations that reproduce the accretion of one or two satellites with their globular cluster population on a Milky Way-type galaxy. Results. Our results demonstrate that a significant overlap between accreted and ‘kinematically heated’ in situ globular clusters is expected in kinematic spaces for mergers with mass ratios of 1:10. In contrast with the standard assumptions made in the literature so far, we find that accreted globular clusters do not show dynamical coherence, that is, they do not cluster in kinematic spaces. In addition, we show that globular clusters can also be found in regions dominated by stars that have a different origin (i.e. a different progenitor). This casts doubt on the association between clusters and field stars that is generally made in the literature and is used to assign them to a common origin. By means of Gaussian mixture models, we demonstrate that the overlap of clusters is not only a projection effect on specific planes but is also found when the whole set of kinematic properties (i.e. E, Lz, Lperp, eccentricity, radial, and vertical actions) is taken into account. Overall, our findings severely question the recovered accretion history of the Milky Way based on the phase-space clustering of the globular cluster population.

The C/N ratio from FUV spectroscopy as a constraint on evolution of the dwarf nova HS 0218 + 3229

ABSTRACT White dwarfs that accrete from non-degenerate companions show anomalous carbon and nitrogen abundances in the photospheres of their stellar components have been postulated to be descendants of supersoft X-ray binaries. Measuring the carbon-to-nitrogen abundance ratio may provide constraints on their past evolution. We fit far-ultraviolet spectroscopy of the cataclysmic variable HS 0218 + 3229 taken with the Cosmic Origins Spectrograph using Markov chain Monte Carlo methods, and found the carbon-to-nitrogen ratio is about one tenth of the Solar value $(\rm{\log \mathrm{[C/N]}}=-0.56\pm 0.15)$. We also provide estimates of the silicon and aluminium abundances, and upper limits for iron and oxygen. Using the parameters we derived for HS 0218 + 3229 we reconstruct its past. We calculated a grid of mesa models and implemented Gaussian process fits in order to determine its most likely initial binary configuration. We found that an initial mass of the donor of $M_{\rm donor;i}=0.90-0.98,\rm{\mathrm{M}_{\odot }}$ and an initial orbital period of Porb; i = 2.88 d (Porb; i = 3.12–3.16 d) for an assumed initial white dwarf mass of $\rm{M_{\mathrm{WD}}}_\mathrm{;i}=0.83\, \rm{\mathrm{M}_{\odot }}$$(\rm{M_{\mathrm{WD}}}_{\rm ;i}=0.60\, \rm{\mathrm{M}_{\odot }})$ can replicate the measured parameters. The low mass ratio, $M_{\rm donor;i} / \rm{M_{\mathrm{WD}}}_{\rm ;i} =1.08-1.18\, (1.5-1.63)$, suggests that the system did not go through a phase of hydrogen-burning on the white dwarf’s surface. However, we can not exclude a phase of thermal time-scale mass transfer in the past. We predict that HS 0218 + 3229 will evolve below the ≃ 76.2 ± 1 min period minimum for normal cataclysmic variables.

The progenitor galaxies of stellar haloes as ‘failed’ Milky Ways

ABSTRACT The stellar halo of the Milky Way records the history of its interactions with dwarf galaxies, whose subsequent destruction results in the formation of an extended stellar component. Recent works have suggested that galaxies with masses comparable to the Large Magellanic Cloud (LMC, $M_\star \sim 10^9$ M☉) may be the primary building blocks of the stellar halo of our Galaxy. We use cosmological simulations of the Lambda cold dark matter model to investigate LMC-mass galaxies at $z$ = 1–2 using a semi-analytical model of galaxy formation. We find that LMC analogues at $z$ = 2 evolve until the present day along three distinct pathways: (i) those that are destroyed in Milky Way-mass hosts; (ii) those that are themselves the main progenitors of Milky Way-mass galaxies; and (iii) those that survive until $z$ = 0, with stellar mass ∼1.0 dex lower than typical Milky Ways. We predict that the properties of these galaxies at $z$ = 2 (stellar metallicities, sizes, gas content, etc.) are largely indistinguishable, irrespective of which of these pathways is eventually taken; a survey targeting such galaxies in this redshift range would struggle to tell apart a ‘destroyed’ stellar halo progenitor from a ‘surviving’ LMC analogue. The only factor that determines the eventual fate of these galaxies is their proximity to a neighbouring Milky Way main progenitor at $z$ = 2: while the mean separation to a ‘surviving’ galaxy is around 7 Mpc, it is only 670 kpc to a ‘destroyed’ galaxy. This suggests that old stellar populations in the Milky Way may share intrinsic (i.e. non-dynamical) properties that are essentially indistinguishable from progenitors of its stellar halo.

Catching Tidal Dwarf Galaxies at a Later Evolutionary Stage with ALFALFA

We present deep optical imaging and photometry of four objects classified as “Almost-Dark” galaxies in the Arecibo Legacy Fast Arecibo L-band Feed Array (ALFALFA) survey because of their gas-rich nature and extremely faint or missing optical emission in existing catalogs. They have H i masses of 107–109 M ⊙ and distances of ∼9–100 Mpc. Observations with the WIYN 3.5 m telescope and One Degree Imager reveal faint stellar components with central surface brightnesses of ∼24–25 in the g band. We also present the results of H i synthesis observations with the Westerbork Synthesis Radio Telescope. These Almost-Dark galaxies have been identified as possible tidal dwarf galaxies (TDGs) based on their proximity to one or more massive galaxies. We demonstrate that AGC 229398 and AGC 333576 likely have the low dark matter content and large effective radii representative of TDGs. They are located much farther from their progenitors than previously studied TDGs, suggesting they are older and more evolved. AGC 219369 is likely dark matter dominated, while AGC 123216 has a dark matter content that is unusually high for a TDG, but low for a normal dwarf galaxy. We consider possible mechanisms for the formation of the TDG candidates such as a traditional major merger scenario and gas ejection from a high-velocity flyby. Blind H i surveys like ALFALFA enable the detection of gas-rich, optically faint TDGs that can be overlooked in other surveys, thereby providing a more complete census of the low-mass galaxy population and an opportunity to study TDGs at a more advanced stage of their life cycle.

Galaxy pairs inThe Three Hundredsimulations II: studying bound ones and identifying them via machine learning

ABSTRACTUsing the data set of The Three Hundred project, i.e. 324 hydrodynamical resimulations of cluster-sized haloes and the regions of radius 15 ${{h^{-1}\, {\rm Mpc}}}$ around them, we study galaxy pairs in high-density environments. By projecting the galaxies’ 3D coordinates onto a 2D plane, we apply observational techniques to find galaxy pairs. Based on a previous theoretical study on galaxy groups in the same simulations, we are able to classify the observed pairs into ‘true’ or ‘false’, depending on whether they are gravitationally bound or not. We find that the fraction of true pairs (purity) crucially depends on the specific thresholds used to find the pairs, ranging from around 30 to more than 80 per cent in the most restrictive case. Nevertheless, in these very restrictive cases, we see that the completeness of the sample is low, failing to find a significant number of true pairs. Therefore, we train a machine learning algorithm to help us identify these true pairs based on the properties of the galaxies that constitute them. With the aid of the machine learning model trained with a set of properties of all the objects, we show that purity and completeness can be boosted significantly using the default observational thresholds. Furthermore, this machine learning model also reveals the properties that are most important when distinguishing true pairs, mainly the size and mass of the galaxies, their spin parameter, gas content, and shape of their stellar components.

Dark matter subhalo disturbance in the tidal field: Towards a new simulations accuracy

Dark matter (DM) identification is a pressing challenge for modern particle physics and cosmology. The success of this endeavor, however, is contingent on a number of parameters, some of which stem from the difficulty of predicting the DM sub-halo distribution in Milky Way-like galaxies and others from the effect of DM sub-halo signals on the gamma-ray instruments being utilized. This research compared the results of detailed simulations of individual DM sub-halos emerging in the external tidal fields with and without disc components and a baryonic bulge. The average DM particle is assumed to have a mass three times smaller than that assumed in cosmological zoom-in simulations of galaxy formation. The Via Lactea II simulation, which provides the basis for simulating the sub-haloes' tidal field in DM, was used to establish the baseline conditions for this investigation. A model of the Milky Way based on observations was used to simulate the tidal field; this model included a DM halo, a stellar bulge, and a stellar disc. When small (rp 24 kpc) percenters were included, the number of sub-haloes dropped to 39%, and the remaining sub-haloes were found to have a less dense orbital distribution. Since the total number of outer sub-haloes was only around 74% of the value anticipated by CDM, the addition of a stellar component to the tidal field had little effect on the sub-haloes with greater percenters. Previous research on the impacts of baryonic matter as mediated via the spurious disruption had reported ratios that were lower than the ones obtained in this study, suggesting that the latter's accuracy was to blame.

Lyman continuum leaker candidates among highly ionised, low-redshift dwarf galaxies selected from He II

Context. Contemporary research suggests that the reionisation of the intergalactic medium (IGM) in the early Universe was predominantly realised by star-forming (proto-)galaxies (SFGs). Due to observational constraints, our knowledge on the origins of sufficient amounts of ionising Lyman continuum (LyC) photons and the mechanisms facilitating their transport into the IGM remains sparse. Recent efforts have thus focussed on the study of local analogues to these high-redshift objects. Aims. We aim to acquire a set of very low-redshift SFGs that exhibit signs of a hard radiation field being present. A subsequent analysis of their emission line properties is intended to shed light on how the conditions prevalent in these objects compare to those predicted to be present in early SFGs that are thought to be LyC emitters (LCEs). Methods. We used archival spectroscopic SDSS DR12 data to select a sample of low-redshift He II 4686 emitters and restricted it to a set of SFGs with an emission line diagnostic sensitive to the presence of an active galactic nucleus, which serves as our only selection criterion. We performed a population spectral synthesis with FADO to reconstruct these galaxies’ star-formation histories (SFHs). Utilising the spectroscopic information at hand, we constrained the predominant ionisation mechanisms in these galaxies and inferred information on interstellar medium (ISM) conditions relevant for the escape of LyC radiation. Results. Our final sample consists of eighteen ionised, metal-poor galaxies (IMPs). These low-mass (6.2 ≤ log(M⋆/M⊙) ≤ 8.8), low-metallicity (7.54 ≤ log(O/H) + 12 ≤ 8.13) dwarf galaxies appear to be predominantly ionised by stellar sources. We find large [OIII] 5007/[OII] 3727 ratios and [SII] 6717,6731/Hα deficiencies, which provide strong indications for these galaxies to be LCEs. At least 40% of these objects are candidates for featuring cosmologically significant LyC escape fractions ≳10%. The IMPs’ SFHs exhibit strong similarities and almost all galaxies appear to contain an old (> 1 Gyr) stellar component, while also harbouring a young, two-stage (∼10 Myr and < 1 Myr) starburst, which we speculate might be related to LyC escape. Conclusions. The properties of the compact emission line galaxies presented here align well with those observed in many local LCEs. In fact, our sample may prove as an extension to the rather small catalogue of local LCEs, as the extreme ISM conditions we find are assumed to facilitate LyC leakage. Notably, all of our eighteen candidates are significantly closer (z < 0.1) than most established LCEs. If the inferred LyC photon loss is genuine, this demonstrates that selecting SFGs from He II 4686 is a powerful selection criterion in the search for LCEs.

A trail of the invisible: blue globular clusters trace the radial density distribution of the dark matter – case study of NGC 4278

ABSTRACTWe present new, deep optical observations of the early-type galaxy NGC 4278, which is located in a small loose group. We find that the galaxy lacks fine substructure, that is, it appears relaxed, out to a radius of ∼70 kpc. Our g- and i-band surface brightness profiles are uniform down to our deepest levels of ∼28 mag arcsec−2. This spans an extremely large radial range of more than 14 half-mass radii. Combined with archival globular cluster (GC) number density maps and a new analysis of the total mass distribution obtained from archival Chandra X-ray data, we find that the red GC subpopulation traces well the stellar mass density profile from 2.4 out to even 14 half-mass radii, while the blue GC subpopulation traces the total mass density profile of the galaxy over a large radial range. Our results reinforce the scenario that red GCs form mostly in situ along with the stellar component of the galaxy, while the blue GCs are more closely aligned with the total mass distribution in the halo and were accreted along with halo matter. We conclude that for galaxies where the X-ray emission from the hot halo is too faint to be properly observable and as such is not available to measure the dark matter profile, the blue GC population can be used to trace this dark matter component out to large radii.

Conditions for clump survival in High-zdisc galaxies

ABSTRACTWe study the survival of giant clumps in high-redshift disc galaxies, short-lived (S) versus long-lived (L), and two L subtypes, via analytic modelling and simulations. We develop a criterion for clump survival, with/without gas, based on a survivability parameter S. It compares the energy sources by supernova feedback and gravitational contraction to the clump binding energy and losses by outflows and turbulence dissipation. The clump properties are derived from Toomre instability, approaching virial/Jeans equilibrium, and the supernova energy deposit uses an up-to-date bubble analysis. For moderate feedback, we find L clumps with circular velocities ${\sim}50\, {\rm km}\, {\rm s}^{-1}$ and masses ≥108 M⊙. They favour galaxies with circular velocities ${\ge}200\, {\rm km\,s}^{-1}$, consistent at z ∼ 2 with the typical disc stellar mass, ≥109.3 M⊙. L clumps favour disc gas fractions ≥0.3, low-mass bulges, and z ∼ 2. They disfavour more effective feedback due to, e.g. supernova clustering, very strong radiative feedback, top-heavy stellar mass function, or particularly high star-formation-rate (SFR) efficiency. A subtype of L clumps (LS), which lose their gas in several free-fall times but retain bound stellar components, may be explained by less contraction and stronger gravitational effects, where clump mergers increase the SFR efficiency. These may give rise to globular clusters. The more massive L clumps (LL) retain most of their baryons for tens of free-fall times with a roughly constant star-formation rate.

Age-divided mean stellar populations from full spectrum fitting as the simplified star formation and chemical evolution history of a galaxy: methodology and reliability

ABSTRACT We introduce a practical methodology for investigating the star formation and chemical evolution history of a galaxy: age-divided mean stellar populations (ADPs) from full spectrum fitting. In this method, the mass-weighted mean stellar populations and mass fractions (fmass) of young and old stellar components in a galaxy are separately estimated, which are divided with an age cut (selected to be 109.5 yr ≈3.2 Gyr in this paper). To examine the statistical reliability of ADPs, we generate 10 000 artificial galaxy spectra, each of which consists of five random simple stellar population components. Using the Penalized PiXel-Fitting (ppxf) package, we conduct full spectrum fitting to the artificial spectra with noise as a function of wavelength, imitating the real noise of Sydney-Australian Astronomical Observatory Multi-object Integral field spectrograph (SAMI) galaxies. As a result, the Δ (= output − input) of age and metallicity appears to significantly depend on not only signal-to-noise ratio (S/N), but also luminosity fractions (flum) of young and old components. At given S/N and flum, Δ of young components tends to be larger than Δ of old components; e.g. σ(Δ[M/H]) ∼ 0.40 versus 0.23 at S/N = 30 and flum = 50 per cent. The age-metallicity degeneracy appears to be insignificant, but Δlog(age/yr) shows an obvious correlation with Δfmass for young stellar components ($\mathcal {R}\sim 0.6$). The impact of dust attenuation and emission lines appears to be mostly insignificant. We discuss how this methodology can be applied to spectroscopic studies of the formation histories of galaxies, with a few examples of SAMI galaxies.

Evolution of binary stars in the early evolutionary phases of ultra-faint dwarf galaxies

ABSTRACT The dynamics of binary stars provides a unique avenue to gather insight into the study of the structure and dynamics of star clusters and galaxies. In this paper, we present the results of a set of N-body simulations aimed at exploring the evolution of binary stars during the early evolutionary phases of ultra-faint dwarf galaxies (UFD). In our simulations, we assume that the stellar component of the UFD is initially dynamically cold and evolves towards its final equilibrium after undergoing the violent relaxation phase. We show that the early evolutionary phases of the UFD significantly enhance the disruption of wide binaries and leave their dynamical fingerprints on the semimajor axis distribution of the surviving binaries as compared to models initially in equilibrium. An initially thermal eccentricity distribution is preserved except for the widest binaries for which it evolves towards a superthermal distribution; for a binary population with an initially uniform eccentricity distribution, memory of this initial distribution is rapidly lost for most binaries as wider binaries evolve to approach a thermal/superthermal distribution. The evolution of binaries is driven both by tidal effects due to the potential of the UFD dark matter halo and collisional effects associated to binary–binary/single star encounters. Collisional effects are particularly important within the clumpy substructure characterizing the system during its early evolution; in addition to enhancing binary ionization and evolution of the binary orbital parameters, encounters may lead to exchanges of either of the primordial binary components with one of the interacting stars.

Rubin LSST Observing Strategies to Maximize Volume and Uniformity Coverage of Star-forming Regions in the Galactic Plane

A complete map of the youngest stellar populations of the Milky Way in the era of all-sky surveys is one of the most challenging goals in modern astrophysics. The characterization of the youngest stellar components is crucial not only for a global overview of the Milky Way’s structure, of the Galactic thin disk, and its spiral arms, but also for local studies. In fact, the identification of star-forming regions (SFRs) and the comparison with the environment in which they form are also fundamental to put SFRs in the context of the surrounding giant molecular clouds and to understand still unknown physical mechanisms related to star and planet formation processes. In 10 yr of observations, the Vera C. Rubin Legacy Survey of Space and Time (Rubin LSST) will achieve an exquisite photometric depth that will allow us to significantly extend the volume within which we will be able to discover new SFRs and to enlarge the region of our own Galaxy we have detailed knowledge about. We describe here a metric that estimates the total number of young stars with ages t < 10 Myr and masses >0.3 M ⊙ that will be detected with the Rubin LSST observations in the gri bands at a 5σ magnitude significance. We examine the results of our metric adopting the most recent simulated Rubin LSST survey strategies in order to evaluate the impact that different observing strategies might have on our science case.

Assessing indirect methods to determine black hole masses using NGC 4151

ABSTRACT Accurately determining the black hole mass (MBH) in active galactic nuclei (AGNs) is crucial to constraining their properties and to studying their evolution. While direct methods yield reliable measurements of MBH in unobscured type 1 AGNs, where the dynamics of stellar or gas components can be directly observed, only indirect methods can be applied to the vast majority of heavily absorbed type 2 AGNs, which represent most of the AGN population. Since it is difficult to evaluate the accuracy and precision of these indirect methods, we utilize the nearby X-ray bright Seyfert galaxy NGC 4151, whose MBH has been tightly constrained with several independent direct methods, as a laboratory to assess the reliability of three indirect methods that have been applied to obscured AGNs. All three, the X-ray scaling method, the Fundamental Plane of BH activity, and the M–σ correlation, yield MBH values consistent with those inferred from direct methods and can therefore be considered accurate. However, only the X-ray scaling method and the M–σ correlation are precise because the substantial scatter in the Fundamental Plane of BH activity allows only for crude estimates. Of the four M–σ correlations we used, only the one from Kormendy and Ho yields a value consistent with the dynamical estimates. This study suggests that the best approach to estimating the BH mass in systems where direct dynamical methods cannot be applied is to utilize a combination of indirect methods, taking into account their different ranges of applicability.

Mass Models of the Milky Way and Estimation of Its Mass from the Gaia DR3 Data Set

We use data from the Gaia DR3 data set to estimate the mass of the Milky Way (MW) by analyzing the rotation curve in the range of distances 5 to 28 kpc. We consider three mass models: The first model adds a spherical dark matter (DM) halo, following the Navarro–Frenk–White (NFW) profile, to the known stellar components. The second model assumes that DM is confined to the Galactic disk, following the idea that the observed density of gas in the Galaxy is related to the presence of a more massive DM disk (DMD), similar to the observed correlation between DM and gas in other galaxies. The third model only uses the known stellar-mass components and is based on the Modified Newton Dynamics (MOND) theory. Our results indicate that the DMD model is comparable in accuracy to the NFW and MOND models and fits the data better at large radii where the rotation curve declines but has the largest errors. For the NFW model, we obtain a virial mass M vir = (6.5 ± 0.3) × 1011 M ⊙ with concentration parameter c = 14.5, which is lower than what is typically reported. In the DMD case, we find that the MW mass is M d = (1.6 ± 0.5) × 1011 M ⊙ with a disk’s characteristic radius of R d = 17 kpc.

ISM metallicity variations across spiral arms in disk galaxies

Chemical abundance variations in the interstellar medium provide important information about the galactic evolution, star-formation, and enrichment histories. Recent observations of disk galaxies suggest that if large-scale azimuthal metallicity variations appear in the ISM, they are linked to the spiral arms. In this work, using a set of chemodynamical simulations of the Milky Way-like spiral galaxies, we quantify the impact of gas radial motions (migration) in the presence of a pre-existing radial metallicity gradient and the local ISM enrichment on both global and local variations of the mean ISM metallicity in the vicinity of the spiral arms. In all the models, we find the scatter of the gas metallicity of ≈0.04 − 0.06 dex at a given galactocentric distance. On large scales, we observe the presence of spiral-like metallicity patterns in the ISM which are more prominent in models with the radial metallicity gradient. However, in our simulations, the morphology of the large-scale ISM metallicity distributions significantly differs from the spiral arm structure in stellar and gas components resulting in both positive and negative residual (after subtraction of the radial gradient) metallicity trends along spiral arms. We discuss the correlations of the residual ISM metallicity values with the star formation rate, gas kinematics and offset to the spiral arms, concluding that the presence of a radial metallicity gradient is essential for the azimuthal variations of metallicity. At the same time, the local enrichment alone is unlikely to drive systematic variations of the metallicity across the spirals.

Dynamics of Molecular Gas in the Central Region of the Quasar I Zwicky 1

We present a study of the molecular gas distribution and kinematics in the cicumnuclear region (radii ≲2 kpc) of the z ≈ 0.061 quasar I Zwicky 1 using a collection of available Atacama Large Millimeter/submillimeter Array observations of the carbon monoxide (CO) emission. With an angular resolution of ∼0.″36 (corresponding to ∼400 pc), the host-galaxy substructures including the nuclear molecular gas disk, spiral arms, and a compact bar-like component are resolved. We analyzed the gas kinematics based on the CO image cube and obtained the rotation curve and radial distribution of velocity dispersion. The velocity dispersion is about 30 km s−1 in the outer CO disk region and rises up to ≳100 km s−1 at radius ≲1 kpc, suggesting that the central region of the disk is dynamically hot. We constrain the CO-to-H2 conversion factor, α CO, by modeling the cold gas disk dynamics. We find that, with prior knowledge about the stellar and dark matter components, the α CO value in the circumnuclear region of this quasar host galaxy is , which is between the value reported in ultraluminous infrared galaxies and in the Milky Way. The central 1 kpc region of this quasar host galaxy has significant star formation activity, which can be identified as a nuclear starburst. We further investigate the high-velocity dispersion in the central region. We find that the interstellar medium (ISM) turbulent pressure derived from the gas velocity dispersion is in equilibrium with the weight of the ISM. This argues against extra power from active galactic nuclei feedback that significantly affects the kinematics of the cold molecular gas.

The impact of gravity on the evolution of cavity in the cluster of stars

This paper investigates the evolution of cavities for the cluster of stars in the context of extended Gauss–Bonnet gravity. We assume a spherically symmetric geometry along with locally anisotropic fluid distribution. It is assumed that the proper radial distance among neighboring stellar components stays unchanged during purely areal evolution stage. We provide some analytical solutions by using general formulism in gravitation theory. The thick-shells cavities at one or both boundary surfaces are found to satisfy the Darmois conditions. Moreover, we also investigate the physical behavior of cavity models by considering the stellar , and . We conclude that the dark matter has a strong impact on the evolution of cavities in the cluster of stars.

Unravelling the mass spectrum of destroyed dwarf galaxies with the metallicity distribution function

ABSTRACTAccreted stellar populations are comprised of the remnants of destroyed galaxies, and often dominate the ‘stellar haloes’ of galaxies such as the Milky Way (MW). This ensemble of external contributors is a key indicator of the past assembly history of a galaxy. We introduce a novel statistical method that uses the unbinned metallicity distribution function (MDF) of a stellar population to estimate the mass spectrum of its progenitors. Our model makes use of the well-known mass–metallicity relation of galaxies and assumes Gaussian MDF distributions for individual progenitors: the overall MDF is thus a mixture of MDFs from smaller galaxies. We apply the method to the stellar halo of the MW, as well as the classical MW satellite galaxies. The stellar components of the satellite galaxies have relatively small sample sizes, but we do not find any evidence for accreted populations with L &amp;gt; Lhost/100. We find that the MW stellar halo has N ∼ 1−3 massive progenitors (L ≳ 108L⊙) within 10 kpc, and likely several hundred progenitors in total. We also test our method on simulations of MW-mass haloes, and find that our method is able to recover the true accreted population within a factor of 2. Future data sets will provide MDFs with orders of magnitude more stars, and this method could be a powerful technique to quantify the accreted populations down to the ultra-faint dwarf mass scale for both the MW and its satellites.

Baryonic Effects on Lagrangian Clustering and Angular Momentum Reconstruction

Recent studies illustrate the correlation between the angular momenta of cosmic structures and their Lagrangian properties. However, only baryons are observable and it is unclear whether they reliably trace the cosmic angular momenta. We study the Lagrangian mass distribution, spin correlation, and predictability of dark matter, gas, and stellar components of galaxy–halo systems using IllustrisTNG, and show that the primordial segregations between components are typically small. Their protoshapes are also similar in terms of the statistics of moment of inertia tensors. Under the common gravitational potential they are expected to exert the same tidal torque and the strong spin correlations are not destroyed by the nonlinear evolution and complicated baryonic effects, as confirmed by the high-resolution hydrodynamic simulations. We further show that their late-time angular momenta traced by total gas, stars, or the central galaxies, can be reliably reconstructed by the initial perturbations. These results suggest that baryonic angular momenta can potentially be used in reconstructing the parameters and models related to the initial perturbations.

The Formation of the Brightest Cluster Galaxy and Intracluster Light in Cosmological N-body Simulations with the Galaxy Replacement Technique

We investigate the formation channels of the intracluster light (ICL) and the brightest cluster galaxy (BCG) in clusters at z = 0. For this, we perform multi-resolution cosmological N-body simulations using the “galaxy replacement technique.” We study the formation channels of the ICL and BCG as a function of distance from the cluster center and the dynamical state of the clusters at z = 0. To do this, we trace back the stars of the ICL and BCG, and identify the stellar components in which they existed when they first fell into the clusters. We find that the progenitors of the ICL and BCG in the central region of the cluster fell earlier and with a higher total mass ratio of the progenitors to the cluster compared to the outer region. This causes a negative radial gradient in the infall time and total mass ratio of the progenitors. Although stellar mass of the progenitors does not show the same radial gradient in all clusters, massive galaxies (M gal > 1010 M ⊙ h−1) are the dominant formation channel of the ICL and BCG for all clusters, except for our most relaxed cluster. For clusters that are dynamically more unrelaxed, we find that the progenitors of the ICL and BCG fall into their clusters more recently, and with a higher mass and mass ratio. Furthermore, we find that the diffuse material of massive galaxies and group-mass halos that is formed by preprocessing contributes significantly to the ICL in the outer region of the unrelaxed clusters.