Showing 1–50 of 170 results for author: Yi, S K

The SAMI Galaxy Survey: impact of star formation and AGN feedback processes on the ionized gas velocity dispersion

Authors: Sree Oh, Matthew Colless, Stefania Barsanti, Henry R. M. Zovaro, Scott M. Croom, Sukyoung K. Yi, Andrei Ristea, Jesse van de Sande, Francesco D'Eugenio, Joss Bland-Hawthorn, Julia J. Bryant, Sarah Casura, Hyunjin Jeong, Sarah M. Sweet, Tayyaba Zafar

Abstract: We investigate the influence of star formation and instantaneous AGN feedback processes on the ionized gas velocity dispersion in a sample of 1285 emission-line galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9$ from the integral-field spectroscopy SAMI Galaxy Survey. We fit both narrow and broad emission line components using aperture spectra integrated within one effective radius, while… ▽ More We investigate the influence of star formation and instantaneous AGN feedback processes on the ionized gas velocity dispersion in a sample of 1285 emission-line galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9$ from the integral-field spectroscopy SAMI Galaxy Survey. We fit both narrow and broad emission line components using aperture spectra integrated within one effective radius, while ensuring the elimination of velocity differences between the spectra of individual spaxels. Our analysis reveals that 386 (30%) galaxies can be adequately described using a single emission component while 356 (28%) galaxies require two (broad and narrow) components. Galaxies characterized by high mass, elevated star formation rate surface density, or type-2 AGN-like emissions tend to feature an additional broad emission-line component, leading to their classification as double-component galaxies. We explore the correlations between $M_*$ and gas velocity dispersions, highlighting that the prominence of the broad component significantly contributes to elevating the gas velocity dispersion. Galaxies displaying AGN-like emission based on optical definitions show enhanced gas velocity dispersions. In star-forming galaxies, both stellar mass and star-formation rate surface density substantially contribute to the velocity dispersion of the narrow component. Increased star-forming activity appears to elevate the velocity dispersion of the narrow component. The broad component exhibits a weaker dependence on stellar mass and is primarily driven by galactic outflows. We suggest that strong star forming activity leads to the formation of a broad emission-line component, but the impact on inflating gas velocity dispersion is moderate. On the other hand, AGN-driven outflows appear to be a more important contributor to the elevated velocity dispersion of the ionized gas. △ Less

Submitted 31 May, 2024; originally announced May 2024.

Comments: 17 pages, 11 figures

The Impact of Positive AGN Feedback on the Properties of Galaxies in a Semi-Analytic Model of Galaxy Formation

Authors: Emanuele Contini, Sukyoung K. Yi, Seyoung Jeon, Jinsu Rhee

Abstract: We introduce the state-of-the-art semi-analytic model FEGA (Formation and Evolution of GAlaxies), which incorporates updated prescriptions for key physical processes in galaxy formation. Notably, FEGA features an unprecedented semi-analytic modeling of positive Active Galactic Nuclei (AGN) feedback. The model combines the latest prescriptions for gas infall and cooling, a revised star formation re… ▽ More We introduce the state-of-the-art semi-analytic model FEGA (Formation and Evolution of GAlaxies), which incorporates updated prescriptions for key physical processes in galaxy formation. Notably, FEGA features an unprecedented semi-analytic modeling of positive Active Galactic Nuclei (AGN) feedback. The model combines the latest prescriptions for gas infall and cooling, a revised star formation recipe that incorporates the extended Kennicutt-Schmidt relation, disk instability, updated supernovae feedback, reincorporation of ejected gas, hot gas stripping from satellite galaxies, and the formation of diffuse light. A novel description of AGN feedback is introduced, describing the positive mode as a burst of star formation from a cooling gas fraction. FEGA is rigorously calibrated using an MCMC procedure to match the evolution of the stellar mass and K-band luminosity functions from high redshift to the present, and the local black hole-bulge mass relation. Subsequently, the model is tested against several observed and predicted scaling relations, including the star formation rate-mass, black hole-bulge and stellar mass, stellar-to-halo mass, and red fraction-mass relations. Additionally, we test FEGA against other galaxy properties such as the distribution of specific star formation rates, stellar metallicity and morphology. Our results demonstrate that the inclusion of positive AGN feedback can co-exist with its negative counterpart without drastic alterations to other prescriptions. Importantly, this inclusion improves the ability of the model to describe the primary scaling relations observed in galaxies. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: 30 pages, 12 figures and 1 table. Submitted to ApJS, comments very welcome

Formation pathways of the compact stellar systems

Authors: J. K. Jang, Sukyoung K. Yi, Soo-Chang Rey, Jinsu Rhee, Yohan Dubois, Taysun Kimm, Christophe Pichon, Katarina Kraljic, Suk Kim

Abstract: The formation pathways of compact stellar systems (CSSs) are still under debate. We utilize the \NH\ simulation to investigate the origins of such objects in the field environment. We identified 55 CSS candidates in the simulation whose properties are similar to those of the observed ultra-compact dwarfs and compact ellipticals. All but two most massive objects (compact elliptical candidates) are… ▽ More The formation pathways of compact stellar systems (CSSs) are still under debate. We utilize the \NH\ simulation to investigate the origins of such objects in the field environment. We identified 55 CSS candidates in the simulation whose properties are similar to those of the observed ultra-compact dwarfs and compact ellipticals. All but two most massive objects (compact elliptical candidates) are a result of a short starburst. Sixteen are formed by tidal stripping, while the other 39 are intrinsically compact from their birth. The stripped objects originate from dwarf-like galaxies with a dark halo, but most of their dark matter is stripped through their orbital motion around a more massive neighbor galaxy. The 39 intrinsically compact systems are further divided into ``associated'' or ``isolated'' groups, depending on whether they were born near a massive dark halo or not. The isolated intrinsic compact objects (7) are born in a dark halo and their stellar properties are older and metal-poor compared to the associated counterparts (32). The stripped compact objects occupy a distinct region in the age-metallicity plane from the intrinsic compact objects. The associated intrinsic compact objects in our sample have never had a dark halo; they are the surviving star clumps of a massive galaxy. △ Less

Submitted 16 May, 2024; originally announced May 2024.

Comments: 12 pages, 7 figures, 1 table. submitted to ApJ

Feeding Hidden Monsters: a Super-Eddington accreting Black Hole ~1.5 Gyr after the Big Bang

Authors: Hyewon Suh, Julia Scharwächter, Emanuele Paolo Farina, Federica Loiacono, Giorgio Lanzuisi, Günther Hasinger, Stefano Marchesi, Mar Mezcua, Roberto Decarli, Brian C. Lemaux, Marta Volonteri, Francesca Civano, Sukyoung K. Yi, San Han, Mark Rawlings, Denise Hung

Abstract: Recent James Webb Space Telescope (JWST) observations have revealed a surprisingly abundant population of faint, dusty active galactic nuclei (AGNs) at z~4-7. Together with the presence of supermassive black holes (SMBHs) at z>6, this raises questions about the formation and growth histories of early black holes. Current theories for the formation of seed black holes from the death of the first st… ▽ More Recent James Webb Space Telescope (JWST) observations have revealed a surprisingly abundant population of faint, dusty active galactic nuclei (AGNs) at z~4-7. Together with the presence of supermassive black holes (SMBHs) at z>6, this raises questions about the formation and growth histories of early black holes. Current theories for the formation of seed black holes from the death of the first stars (i.e. light seeds) and/or the direct collapse of primordial gas clouds (i.e. heavy seeds) still lack observational confirmation. Here, we present LID-568, a low-mass (7.2e6Msun) black hole hosting powerful outflows that is observed in an extreme phase of rapid growth at z~4. This object is similar to other JWST-discovered faint AGN populations, but is bright in X-ray emission and accreting at more than 4000% of the limit at which radiation pressure exceeds the force of gravitational attraction of the black hole (i.e. super-Eddington accretion). Analysis of JWST NIRSpec/IFU data reveals spatially extended Ha emission with velocities of ~ -600 - -500 km/s relative to the central black hole, indicative of robust nuclear-driven outflows. LID-568 represents an elusive low-mass black hole experiencing super-Eddington accretion as invoked by models of early black hole formation. This discovery showcases a previously undiscovered key parameter space and offers crucial insights into rapid black hole growth mechanisms in the early universe. △ Less

Submitted 8 May, 2024; originally announced May 2024.

Comments: Submitted

Brightest Cluster Galaxies and the Intracluster Light

Authors: Emanuele Contini, Sukyoung K. Yi, Seyoung Jeon

Abstract: In this chapter, we delve into the formation and primary characteristics of two significant components within galaxy clusters: the brightest cluster galaxies and the intracluster light. Drawing upon recent and pertinent studies in the field, we explore the mechanisms driving their growth from high redshift to the present day. We also examine how these formation mechanisms are intertwined with the… ▽ More In this chapter, we delve into the formation and primary characteristics of two significant components within galaxy clusters: the brightest cluster galaxies and the intracluster light. Drawing upon recent and pertinent studies in the field, we explore the mechanisms driving their growth from high redshift to the present day. We also examine how these formation mechanisms are intertwined with the dynamical state of their host clusters, suggesting their potential utility as luminous tracers of dark matter. △ Less

Submitted 1 April, 2024; originally announced April 2024.

Comments: 15 pages, 5 figures. Chapter of the book "Encyclopedia of Astrophysics", 1st Edition (Elsevier). Comments or any feedback welcome

Distribution of merging and post-merging galaxies in nearby galaxy clusters

Authors: Duho Kim, Yun-Kyeong Sheen, Yara L. Jaffé, Kshitija Kelkar, Adarsh Ranjan, Franco Piraino-Cerda, Jacob P. Crossett, Ana Carolina Costa Lourenço, Garreth Martin, Julie B. Nantais, Ricardo Demarco, Ezequiel Treister, Sukyoung K. Yi

Abstract: We study the incidence and spatial distribution of galaxies that are currently undergoing gravitational merging (M) or that have signs of a post merger (PM) in six galaxy clusters (A754, A2399, A2670, A3558, A3562, and A3716) within the redshift range, 0.05$\lesssim$$z$$\lesssim$0.08. To this aim, we obtained Dark Energy Camera (DECam) mosaics in $u^{\prime}$, $g^{\prime}$, and $r^{\prime}$-bands… ▽ More We study the incidence and spatial distribution of galaxies that are currently undergoing gravitational merging (M) or that have signs of a post merger (PM) in six galaxy clusters (A754, A2399, A2670, A3558, A3562, and A3716) within the redshift range, 0.05$\lesssim$$z$$\lesssim$0.08. To this aim, we obtained Dark Energy Camera (DECam) mosaics in $u^{\prime}$, $g^{\prime}$, and $r^{\prime}$-bands covering up to $3\times R_{200}$ of the clusters, reaching 28 mag/arcsec$^2$ surface brightness limits. We visually inspect $u^{\prime}$$g^{\prime}$$r^{\prime}$ color-composite images of volume-limited ($M_r < -20$) cluster-member galaxies to identify whether galaxies are of M or PM types. We find 4% M-type and 7% PM-type galaxies in the galaxy clusters studied. By adding spectroscopic data and studying the projected phase space diagram (PPSD) of the projected clustocentric radius and the line-of-sight velocity, we find that PM-type galaxies are more virialized than M-type galaxies, having 1--5% point higher fraction within the escape-velocity region, while the fraction of M-type was $\sim$10% point higher than PM-type in the intermediate environment. Similarly, on a substructure analysis, M types were found in the outskirt groups, while PM types populated groups in ubiquitous regions of the PPSD. Adopting literature-derived dynamical state indicator values, we observed a higher abundance of M types in dynamically relaxed clusters. This finding suggests that galaxies displaying post-merging features within clusters likely merged in low-velocity environments, including cluster outskirts and dynamically relaxed clusters. △ Less

Submitted 3 May, 2024; v1 submitted 11 March, 2024; originally announced March 2024.

Comments: 23 pages, 15 figures, 4 tables, Published in ApJ. For photometric catalogs and associated information, see https://data.kasi.re.kr/vo/DECam_catalogs/

Journal ref: ApJ 966 124 (2024)

Galaxies with grains: unraveling dust evolution and extinction curves with hydrodynamical simulations

Authors: Yohan Dubois, Francisco Rodríguez Montero, Corentin Guerra, Maxime Trebitsch, San Han, Ricarda Beckmann, Sukyoung K. Yi, Joseph Lewis, J. K. Jang

Abstract: We introduce a model for dust evolution in the RAMSES code for simulations of galaxies with a resolved multiphase interstellar medium. Dust is modelled as a fluid transported with the gas component, and is decomposed into two sizes, 5 nm and 0.1 $μ\rm m$, and two chemical compositions for carbonaceous and silicate grains. Using a suite of isolated disc simulations with different masses and metalli… ▽ More We introduce a model for dust evolution in the RAMSES code for simulations of galaxies with a resolved multiphase interstellar medium. Dust is modelled as a fluid transported with the gas component, and is decomposed into two sizes, 5 nm and 0.1 $μ\rm m$, and two chemical compositions for carbonaceous and silicate grains. Using a suite of isolated disc simulations with different masses and metallicities, the simulations can explore the role of these processes in shaping the key properties of dust in galaxies. The simulated Milky Way analogue reproduces the dust-to-metal mass ratio (DTM), depletion factors, size distribution and extinction curves of the Milky Way. Galaxies with lower metallicities reproduce the observed decrease in the DTM with metallicity at around a few 0.1 $\rm Z_\odot$. This break in the DTM corresponds to a galactic gas metallicity threshold that marks the transition from an ejecta-dominated to an accretion-dominated grain growth, and that is different for silicate and carbonaceous grains, with $\simeq$ 0.1 $\rm Z_\odot$ and $\simeq$ 0.5 $\rm Z_\odot$ respectively. This leads to more Magellanic Cloud-like extinction curves, i.e. with steeper slopes in the ultraviolet and a weaker bump feature at 217.5 nm, in galaxies with lower masses and lower metallicities. Steeper slopes in these galaxies are caused by the combination of the higher efficiency of gas accretion by silicate relative to carbonaceous grains and by the low rates of coagulation that preserves the amount of small silicate grains. Weak bumps are due to the overall inefficient accretion growth of carbonaceous dust at low metallicity, whose growth is mostly supported by the release of large grains in SN ejecta. We also show that the formation of CO molecules is a key component to limit the ability of carbonaceous dust to grow, in particular in low-metallicity gas-rich galaxies. △ Less

Submitted 4 June, 2024; v1 submitted 28 February, 2024; originally announced February 2024.

Comments: A&A accepted

The SAMI Galaxy Survey: galaxy spin is more strongly correlated with stellar population age than mass or environment

Authors: S. M. Croom, J. van de Sande, S. P. Vaughan, T. H. Rutherford, C. P. Lagos, S. Barsanti, J. Bland-Hawthorn, S. Brough, J. J. Bryant, M. Colless, L. Cortese, F. D'Eugenio, A. Fraser-McKelvie, M. Goodwin, N. P. F. Lorente, S. N. Richards, A. Ristea, S. M. Sweet, S. K. Yi, T. Zafar

Abstract: We use the SAMI Galaxy Survey to examine the drivers of galaxy spin, $λ_{R_e}$, in a multi-dimensional parameter space including stellar mass, stellar population age (or specific star formation rate) and various environmental metrics (local density, halo mass, satellite vs. central). Using a partial correlation analysis we consistently find that age or specific star formation rate is the primary p… ▽ More We use the SAMI Galaxy Survey to examine the drivers of galaxy spin, $λ_{R_e}$, in a multi-dimensional parameter space including stellar mass, stellar population age (or specific star formation rate) and various environmental metrics (local density, halo mass, satellite vs. central). Using a partial correlation analysis we consistently find that age or specific star formation rate is the primary parameter correlating with spin. Light-weighted age and specific star formation rate are more strongly correlated with spin than mass-weighted age. In fact, across our sample, once the relation between light-weighted age and spin is accounted for, there is no significant residual correlation between spin and mass, or spin and environment. This result is strongly suggestive that present-day environment only indirectly influences spin, via the removal of gas and star formation quenching. That is, environment affects age, then age affects spin. Older galaxies then have lower spin, either due to stars being born dynamically hotter at high redshift, or due to secular heating. Our results appear to rule out environmentally dependent dynamical heating (e.g. galaxy-galaxy interactions) being important, at least within $1R_e$ where our kinematic measurements are made. The picture is more complex when we only consider high-mass galaxies ($M_*\gtrsim 10^{11}$M$_{\odot}$). While the age-spin relation is still strong for these high-mass galaxies, there is a residual environmental trend with central galaxies preferentially having lower spin, compared to satellites of the same age and mass. We argue that this trend is likely due to central galaxies being a preferred location for mergers. △ Less

Submitted 9 February, 2024; originally announced February 2024.

Comments: 24 pages, 9 figures. Accepted for publication in MNRAS

Journal ref: MNRAS, Volume 529, Issue 4, April 2024, Pages 3446-3468

Diffuse Light in Milky-Way like Haloes

Authors: E. Contini, S. Han, S. Jeon, J. Rhee, S. K. Yi

Abstract: We investigate the diffuse light (DL) content of dark matter haloes in the mass range $11.5\leq \log M_{halo}\leq13$, a range that includes also the dark matter halo of the Milky-Way, taking advantage of a state-of-the-art semi-analytic model run on the merger trees extracted from a set of high-resolution cosmological simulations. The fraction of DL in such relatively small haloes is found to prog… ▽ More We investigate the diffuse light (DL) content of dark matter haloes in the mass range $11.5\leq \log M_{halo}\leq13$, a range that includes also the dark matter halo of the Milky-Way, taking advantage of a state-of-the-art semi-analytic model run on the merger trees extracted from a set of high-resolution cosmological simulations. The fraction of DL in such relatively small haloes is found to progressively decrease from the high to the low mass end, in good agreement with analytic (\citealt{purcell2007}) and numerical results from simulations (\citealt{proctor2023,ahvazi2023}), in good agreement also with the fraction of the DL observed in the Milky-Way (\citealt{deason2019}) and M31 (\citealt{harmsen2017}). Haloes with different masses have a different efficiency in producing DL: $\log M_{halo} \simeq 13$ is found to be the characteristic halo mass where the production of DL is the most efficient, while the overall efficiency decreases at both larger (\citealt{contini2024}) and smaller scales (this work). The DL content in this range of halo mass is the result of stellar stripping due to tidal interaction between satellites and its host (95\%) and mergers between satellites and the central galaxy (5\%), with pre-processed material, sub-channel of mergers and stripping and so already included in the 100\%, that contributes no more than 8\% on average. The halo concentration is the main driver of the DL formation: more concentrated haloes have higher DL fractions that come from stripping of more massive satellites in the high halo mass end, while dwarfs contribute mostly in the low halo mass end. △ Less

Submitted 28 February, 2024; v1 submitted 26 January, 2024; originally announced January 2024.

Comments: 9 pages, 6 figures. Small corrections for typos

Cosmic evolution of black hole-spin and galaxy orientations: clues from the NewHorizon and Galactica simulations

Authors: Sebastien Peirani, Yasushi Suto, Ricarda S. Beckmann, Marta Volonteri, Yen-Ting Lin, Yohan Dubois, Sukyoung K. Yi, Christophe Pichon, Katarina Kraljic, Minjung Park, Julien Devriendt, San Han, Wei-Huai Chen

Abstract: (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the a… ▽ More (Reduced) Using the recent cosmological high-resolution zoom-in simulations, NewHorizon and Galactica, in which the evolution of black hole spin is followed on the fly, we have tracked the cosmic history of a hundred of black holes (BHs) with a mass greater than 2x10^4 Ms. For each of them, we have studied the variations of the three dimensional angle (Psi) subtended between the BH spins and the angular momentum vectors of their host galaxies. The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low mass BHs (<3x10^4 Ms) are rapidly spun up by gas accretion and their spin tends to be aligned with their host galaxy spin. Then follows a second phase in which the accretion of gas onto low mass BHs (<10^5 Ms) is quite chaotic and inefficient, reflecting the complex and disturbed morphologies of forming proto-galaxies at high redshifts. The variations of Psi are rather erratic during this phase and are mainly driven by the rapid changes of the direction of the galaxy angular momentum. Then, in a third and long phase, BHs are generally well settled in the center of galaxies around which the gas accretion becomes much more coherent (>10^5 Ms). In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We have also derived the distributions of cos(Psi) at different redshifts and found that BHs and galaxy spins are generally aligned. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda. In particular, the distribution of lambda traces well the alignment tendency in the 3-d analysis. Such predictions provide an interesting background for future observational analyses. △ Less

Submitted 25 March, 2024; v1 submitted 8 January, 2024; originally announced January 2024.

Comments: Accepted for publication in A&A (12/03/2024)

Journal ref: A&A 686, A233 (2024)

The Connection between the Intracluster Light and its Host Halo: Formation Time and Contribution from Different Channels

Authors: Emanuele Contini, Jinsu Rhee, San Han, Seyoung Jeon, Sukyoung K. Yi

Abstract: We extend the analysis presented in \cite{contini2023a} to higher redshifts, up to $z=2$, by focusing on the relation between the intracluster light (ICL) fraction and the halo mass, its dependence with redshift, role played by the halo concentration and formation time, in a large sample of simulated galaxy groups/clusters with $13\lesssim \log M_{halo} \lesssim 15$. Moreover, a key focus is to is… ▽ More We extend the analysis presented in \cite{contini2023a} to higher redshifts, up to $z=2$, by focusing on the relation between the intracluster light (ICL) fraction and the halo mass, its dependence with redshift, role played by the halo concentration and formation time, in a large sample of simulated galaxy groups/clusters with $13\lesssim \log M_{halo} \lesssim 15$. Moreover, a key focus is to isolate the relative contributions provided by the main channels for the ICL formation to the total amount. The ICL fraction at higher redshift is weakly dependent on halo mass, and comparable with that at the present time, in agreement with recent observations. Stellar stripping, mergers and pre-processing are the major responsible channels of the ICL formation, with stellar stripping that accounts for $\sim 90\%$ of the total ICL, regardless of halo mass and redshift. Pre-processing is an important process for clusters to accrete already formed ICL. The diffuse component forms very early, $z\sim 0.6$, and its formation depends on both concentration and formation time of the halo, with more concentrated and earlier formed haloes that assemble their ICL earlier than later formed ones. The efficiency of this process is independent of halo mass, but increases with decreasing redshift, which implies that stellar stripping becomes more important with time as the concentration increases. This highlights the link between the ICL and the dynamical state of a halo: groups/clusters that have a higher fraction of diffuse light are more concentrated, relaxed and in an advanced stage of growth. △ Less

Submitted 14 November, 2023; v1 submitted 30 October, 2023; originally announced October 2023.

Comments: Corrected for several bugs and typos. Clean now

On the Origin of the Variety of Velocity Dispersion Profiles of Galaxies

Authors: San Han, Sukyoung K. Yi, Sree Oh, Mina Pak, Scott M. Croom, Julien Devriendt, Yohan Dubois, Taysun Kimm, Katarina Kraljic, Christophe Pichon, Marta Volonteri

Abstract: Observed and simulated galaxies exhibit a significant variation in their velocity dispersion profiles. We examine the inner and outer slopes of stellar velocity dispersion profiles using integral field spectroscopy data from two surveys, SAMI (for $z < 0.115$) and CALIFA (for $z < 0.03$), comparing them with results from two cosmological hydrodynamic simulations: Horizon-AGN (for $z = 0.017$) and… ▽ More Observed and simulated galaxies exhibit a significant variation in their velocity dispersion profiles. We examine the inner and outer slopes of stellar velocity dispersion profiles using integral field spectroscopy data from two surveys, SAMI (for $z < 0.115$) and CALIFA (for $z < 0.03$), comparing them with results from two cosmological hydrodynamic simulations: Horizon-AGN (for $z = 0.017$) and NewHorizon (for $z\lesssim1$). The simulated galaxies closely reproduce the variety of velocity dispersion slopes and stellar mass dependence of both inner and outer radii ($0.5\,r_{50}$ and $3\,r_{50}$) as observed, where $r_{50}$ stands for half-light radius. The inner slopes are mainly influenced by the relative radial distribution of the young and old stars formed in-situ: a younger center shows a flatter inner profile. The presence of accreted (ex-situ) stars has two effects on the velocity dispersion profiles. First, because they are more dispersed in spatial and velocity distributions compared to in-situ formed stars, it increases the outer slope of the velocity dispersion profile. It also causes the velocity anisotropy to be more radial. More massive galaxies have a higher fraction of stars formed ex-situ and hence show a higher slope in outer velocity dispersion profile and a higher degree of radial anisotropy. The diversity in the outer velocity dispersion profiles reflects the diverse assembly histories among galaxies. △ Less

Submitted 31 May, 2024; v1 submitted 24 October, 2023; originally announced October 2023.

Comments: 22 pages, 16 figures, accepted for publication in ApJ

The formation of cores in galaxies across cosmic time -- the existence of cores is not in tension with the LCDM paradigm

Authors: R. A. Jackson, S. Kaviraj, S. K. Yi, S. Peirani, Y. Dubois, G. Martin, J. E. G. Devriendt, A. Slyz, C. Pichon, M. Volonteri, T. Kimm, K. Kraljic

Abstract: The `core-cusp' problem is considered a key challenge to the LCDM paradigm. Halos in dark matter only simulations exhibit `cuspy' profiles, where density continuously increases towards the centre. However, the dark matter profiles of many observed galaxies (particularly in the dwarf regime) deviate strongly from this prediction, with much flatter central regions (`cores'). We use NewHorizon (NH),… ▽ More The `core-cusp' problem is considered a key challenge to the LCDM paradigm. Halos in dark matter only simulations exhibit `cuspy' profiles, where density continuously increases towards the centre. However, the dark matter profiles of many observed galaxies (particularly in the dwarf regime) deviate strongly from this prediction, with much flatter central regions (`cores'). We use NewHorizon (NH), a hydrodynamical cosmological simulation, to investigate core formation, using a statistically significant number of galaxies in a cosmological volume. Halos containing galaxies in the upper (M* > 10^10.2 MSun) and lower (M* < 10^8 MSun) ends of the stellar mass distribution contain cusps. However, halos containing galaxies with intermediate (10^8 MSun < M* < 10^10.2 MSun) stellar masses are generally cored, with typical halo masses between 10^10.2 MSun and 10^11.5 MSun. Cores form through supernova-driven gas removal from halo centres, which alters the central gravitational potential, inducing dark matter to migrate to larger radii. While all massive (M* > 10^9.5 MSun) galaxies undergo a cored-phase, in some cases cores can be removed and cusps reformed. This happens if a galaxy undergoes sustained star formation at high redshift, which results in stars (which, unlike the gas, cannot be removed by baryonic feedback) dominating the central gravitational potential. After cosmic star formation peaks, the number of cores, and the mass of the halos they are formed in, remain constant, indicating that cores are being routinely formed over cosmic time after a threshold halo mass is reached. The existence of cores is, therefore, not in tension with the standard paradigm. △ Less

Submitted 8 January, 2024; v1 submitted 19 October, 2023; originally announced October 2023.

Comments: 13 pages, 11 figures, Accepted to MNRAS

The Intracluster Light and its Link with the Dynamical State of the Host Group/Cluster: the Role of the Halo Concentration

Authors: E. Contini, S. Jeon, J. Rhee, S. Han, S. K. Yi

Abstract: We investigate on the role of the halo concentration in the formation of the intra-cluster light (ICL) in galaxy groups and clusters, as predicted by a state-of-art semi-analytic model of galaxy formation, coupled with a set of high-resolution dark matter only simulations. The analysis focuses on how the fraction of ICL correlates with halo mass, concentration and fraction of early-type galaxies (… ▽ More We investigate on the role of the halo concentration in the formation of the intra-cluster light (ICL) in galaxy groups and clusters, as predicted by a state-of-art semi-analytic model of galaxy formation, coupled with a set of high-resolution dark matter only simulations. The analysis focuses on how the fraction of ICL correlates with halo mass, concentration and fraction of early-type galaxies (ETGs) in a large sample of groups and clusters with $13.0\leq \log M_{halo} \leq 15.0$. The fraction of ICL follows a normal distribution, a consequence of the stochastic nature of the physical processes responsible for the formation of the diffuse light. The fractional budget of ICL depends on both halo mass (very weakly) until group scales, and concentration (remarkably). More interestingly, the ICL fraction is higher in more concentrated objects, a result of the stronger tidal forces acting in the innermost regions of the haloes where the concentration is the quantity playing the most relevant role. Our model predictions do not show any dependence between the ICL and ETGs fractions and so, we instead suggest the concentration rather than the mass, as recently claimed, to be the main driver of the ICL formation. The diffuse light starts to form in groups via stellar stripping and mergers and later assembled in more massive objects. However, the formation and assembly keep going on group/cluster scales at lower redshift through the same processes, mainly via stellar stripping in the vicinity of the central regions where tidal forces are stronger. △ Less

Submitted 12 October, 2023; v1 submitted 4 October, 2023; originally announced October 2023.

Comments: 16 pages, 6 figures, 1 table. Accepted for publication in ApJ, small corrections

Emergence and cosmic evolution of the Kennicutt-Schmidt relation driven by interstellar turbulence

Authors: Katarina Kraljic, Florent Renaud, Yohan Dubois, Christophe Pichon, Oscar Agertz, Eric Andersson, Julien Devriendt, Jonathan Freundlich, Sugata Kaviraj, Taysun Kimm, Garreth Martin, Sébastien Peirani, Álvaro Segovia Otero, Marta Volonteri, Sukyoung K. Yi

Abstract: The scaling relations between the gas content and star formation rate of galaxies provide useful insights into processes governing their formation and evolution. We investigate the emergence and the physical drivers of the global Kennicutt-Schmidt (KS) relation at $0.25 \leq z \leq 4$ in the cosmological hydrodynamic simulation NewHorizon capturing the evolution of a few hundred galaxies with a re… ▽ More The scaling relations between the gas content and star formation rate of galaxies provide useful insights into processes governing their formation and evolution. We investigate the emergence and the physical drivers of the global Kennicutt-Schmidt (KS) relation at $0.25 \leq z \leq 4$ in the cosmological hydrodynamic simulation NewHorizon capturing the evolution of a few hundred galaxies with a resolution of $\sim$ 40 pc. The details of this relation vary strongly with the stellar mass of galaxies and the redshift. A power-law relation $Σ_{\rm SFR} \propto Σ_{\rm gas}^{a}$ with $a \approx 1.4$, like that found empirically, emerges at $z \approx 2 - 3$ for the most massive half of the galaxy population. However, no such convergence is found in the lower-mass galaxies, for which the relation gets shallower with decreasing redshift. At the galactic scale, the star formation activity correlates with the level of turbulence of the interstellar medium, quantified by the Mach number, rather than with the gas fraction (neutral or molecular), confirming previous works. With decreasing redshift, the number of outliers with short depletion times diminishes, reducing the scatter of the KS relation, while the overall population of galaxies shifts toward low densities. Using pc-scale star formation models calibrated with local Universe physics, our results demonstrate that the cosmological evolution of the environmental and intrinsic conditions conspire to converge towards a significant and detectable imprint in galactic-scale observables, in their scaling relations, and in their reduced scatter. △ Less

Submitted 12 September, 2023; originally announced September 2023.

Comments: 26 pages, 22 figures

On the Significance of the Thick Disks of Disk Galaxies

Authors: Sukyoung K. Yi, J. K. Jang, Julien Devriendt, Yohan Dubois, San Han, Taysun Kimm, Katarina Kraljic, Minjung Park, Sebastien Peirani, Christophe Pichon, Jinsu Rhee

Abstract: 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 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 thick disk mass fraction. We also use… ▽ More 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 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 thick disk mass fraction. We also use the NewHorizon2 simulation that was run on the same initial conditions but additionally traced nine chemical elements. Based on a sample of 27 massive disk galaxies with M* > 10^10 M_{\odot} in NewHorizon, the contribution of the thick disk was found to be 34 \pm 15% in r-band luminosity or 48 \pm 13% 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 severely contaminated by stars of different kinematic properties while missing a 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 mid-plane with high circularity and get heated with time via fluctuation in the force field. Depending on the star formation and merger histories, galaxies may naturally develop thick disks with significantly different properties. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 28 pages, 24 figures. Submitted to ApJS

The MAGPI Survey: Impact of environment on the total internal mass distribution of galaxies in the last 5 Gyr

Authors: Caro Derkenne, Richard M. McDermid, Adriano Poci, J. Trevor Mendel, Francesco D'Eugenio, Seyoung Jeon, Rhea-Silvia Remus, Sabine Bellstedt, Andrew J. Battisti, Joss Bland-Hawthorn, Anna Ferre-Mateu, Caroline Foster, K. E. Harborne, Claudia D. P. Lagos, Yingjie Peng, Piyush Sharda, Gauri Sharma, Sarah Sweet, Kim-Vy H. Tran, Lucas M. Valenzuela, Sam Vaughan, Emily Wisnioski, Sukyoung K. Yi

Abstract: We investigate the impact of environment on the internal mass distribution of galaxies using the Middle Ages Galaxy Properties with Integral field spectroscopy (MAGPI) survey. We use 2D resolved stellar kinematics to construct Jeans dynamical models for galaxies at mean redshift $z \sim 0.3$, corresponding to a lookback time of $3-4$ Gyr. The internal mass distribution for each galaxy is parameter… ▽ More We investigate the impact of environment on the internal mass distribution of galaxies using the Middle Ages Galaxy Properties with Integral field spectroscopy (MAGPI) survey. We use 2D resolved stellar kinematics to construct Jeans dynamical models for galaxies at mean redshift $z \sim 0.3$, corresponding to a lookback time of $3-4$ Gyr. The internal mass distribution for each galaxy is parameterised by the combined mass density slope $γ$ (baryons $+$ dark matter), which is the logarithmic change of density with radius. We use a MAGPI sample of 28 galaxies from low-to-mid density environments and compare to density slopes derived from galaxies in the high density Frontier Fields clusters in the redshift range $0.29 <z < 0.55$, corresponding to a lookback time of $\sim 5$ Gyr. We find a median density slope of $γ= -2.22 \pm 0.05$ for the MAGPI sample, which is significantly steeper than the Frontier Fields median slope ($γ= -2.01 \pm 0.04$), implying the cluster galaxies are less centrally concentrated in their mass distribution than MAGPI galaxies. We also compare to the distribution of density slopes from galaxies in Atlas3D at $z \sim 0$, because the sample probes a similar environmental range as MAGPI. The Atlas3D median total slope is $γ= -2.25 \pm 0.02$, consistent with the MAGPI median. Our results indicate environment plays a role in the internal mass distribution of galaxies, with no evolution of the slope in the last 3-4 Gyr. These results are in agreement with the predictions of cosmological simulations. △ Less

Submitted 16 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 522, Issue 3, July 2023, Pages 3602 - 3626

The SAMI Galaxy Survey: Environmental analysis of the orbital structures of passive galaxies

Authors: Giulia Santucci, Sarah Brough, Jesse van de Sande, Richard McDermid, Stefania Barsanti, Joss Bland-Hawthorn, Julia J. Bryant, Scott M. Croom, Claudia Lagos, Jon S. Lawrence, Matt S. Owers, Glenn van de Ven, Sam P. Vaughan, Sukyoung K. Yi

Abstract: Most dynamical models of galaxies to date assume axisymmetry, which is not representative of a significant fraction of massive galaxies. We have built triaxial orbit-superposition Schwarzschild models of galaxies observed by the SAMI Galaxy Survey, in order to reconstruct their inner orbital structure and mass distribution. The sample consists of 153 passive galaxies with total stellar masses in t… ▽ More Most dynamical models of galaxies to date assume axisymmetry, which is not representative of a significant fraction of massive galaxies. We have built triaxial orbit-superposition Schwarzschild models of galaxies observed by the SAMI Galaxy Survey, in order to reconstruct their inner orbital structure and mass distribution. The sample consists of 153 passive galaxies with total stellar masses in the range $10^{9.5}$ to $10^{12} M_{\odot}$. We present an analysis of the internal structures and intrinsic properties of these galaxies as a function of their environment. We measure their environment using three proxies: central or satellite designation, halo mass and local $5^{th}$ nearest neighbour galaxy density. We find that although these intrinsic properties correlate most strongly with stellar mass, environment does play a secondary role: at fixed stellar mass, galaxies in the densest regions are more radially anisotropic. In addition, central galaxies, and galaxies in high local densities show lower values of edge-on spin parameter proxy \lam. We also find suggestions of a possible trend of the fractions of orbits with environment for lower-mass galaxies (between $10^{9.5}$ and $10^{11} M_{\odot}$) such that, at fixed stellar mass, galaxies in higher local densities and halo mass have higher fractions of hot orbits and lower fractions of warm orbits. Our results demonstrate that after stellar mass, environment does play a role in shaping present-day passive galaxies. △ Less

Submitted 7 March, 2023; originally announced March 2023.

Comments: 21 pages. Accepted for publication in MNRAS

The SAMI Survey: Evidence for dynamical coupling of ionised gas and young stellar populations

Authors: Caroline Foster, Sam Vaughan, Amelia Fraser-McKelvie, Sarah Brough, Julia J. Bryant, Scott M. Croom, Francesco D'Eugenio, Brent Groves, Iraklis S. Konstantopoulos, Ángel R. López-Sánchez, Sree Oh, Matt S. Owers, Sarah M. Sweet, Jesse van de Sande, Emily Wisnioski, Sukyoung K. Yi, Henry R. M. Zovaro

Abstract: We explore local and global dynamical differences between the kinematics of ionised gas and stars in a sample of galaxies from Data Release 3 of the SAMI Galaxy Survey. We find better agreement between local (i.e., comparing on a spaxel-to-spaxel basis) velocities and dispersion of gas and stars in younger systems as with previous work on the asymmetric drift in galaxies, suggesting that the dynam… ▽ More We explore local and global dynamical differences between the kinematics of ionised gas and stars in a sample of galaxies from Data Release 3 of the SAMI Galaxy Survey. We find better agreement between local (i.e., comparing on a spaxel-to-spaxel basis) velocities and dispersion of gas and stars in younger systems as with previous work on the asymmetric drift in galaxies, suggesting that the dynamics of stars and ionised gas are initially coupled. The intrinsic scatter around the velocity and dispersion relations increases with increasing stellar age and mass, suggesting that subsequent mechanisms such as internal processes, divergent star formation and assembly histories also play a role in setting and altering the dynamics of galaxies. The global (flux-weighted) dynamical support of older galaxies is hotter than in younger systems. We find that the ionised gas in galaxies is almost always dynamically colder than the stars with a steeper velocity gradient. In absolute terms, the local difference in velocity dispersion is more pronounced than the local difference in velocity, possibly reflecting inherent differences in the impact of turbulence, inflow and/or feedback on gas compared to stars. We suggest how these findings may be taken into account when comparing high and low redshift galaxy samples to infer dynamical evolution. △ Less

Submitted 13 February, 2023; originally announced February 2023.

Comments: 11 pages, 9 figures, 1 table, accepted for publication in MNRAS

Translators of galaxy morphology indicators between observation and simulation

Authors: J. K. Jang, Sukyoug K. Yi, Yohan Dubois, Jinsu Rhee, Christophe Pichon, Taysun Kimm, Julien Devriendt, Marta Volonteri, Sugata Kaviraj, Sebastien Peirani, Sree Oh, Scott Croom

Abstract: Based on the recent advancements in the numerical simulations of galaxy formation, we anticipate the achievement of realistic models of galaxies in the near future. Morphology is the most basic and fundamental property of galaxies, yet observations and simulations still use different methods to determine galaxy morphology, making it difficult to compare them. We hereby perform a test on the recent… ▽ More Based on the recent advancements in the numerical simulations of galaxy formation, we anticipate the achievement of realistic models of galaxies in the near future. Morphology is the most basic and fundamental property of galaxies, yet observations and simulations still use different methods to determine galaxy morphology, making it difficult to compare them. We hereby perform a test on the recent NewHorizon simulation which has spatial and mass resolutions that are remarkably high for a large-volume simulation, to resolve the situation. We generate mock images for the simulated galaxies using SKIRT that calculates complex radiative transfer processes in each galaxy. We measure morphological indicators using photometric and spectroscopic methods following observer's techniques. We also measure the kinematic disk-to-total ratios using the Gaussian mixture model and assume that they represent the true structural composition of galaxies. We found that spectroscopic indicators such as $V/σ$ and $λ_{R}$ closely trace the kinematic disk-to-total ratios. In contrast, photometric disk-to-total ratios based on the radial profile fitting method often fail to recover the true kinematic structure of galaxies, especially for small galaxies. We provide translating equations between various morphological indicators. △ Less

Submitted 8 May, 2023; v1 submitted 2 November, 2022; originally announced November 2022.

Comments: 16 pages, 15 figures (including appendix), accepted by ApJ. April 13, 2023

On the Viability of Determining Galaxy Properties from Observations I: Star Formation Rates and Kinematics

Authors: Kearn Grisdale, Laurence Hogan, Dimitra Rigopoulou, Niranjan Thatte, Miguel Pereira-Santaella, Julien Devriendt, Adrianne Slyz, Ismael García-Bernete, Yohan Dubois, Sukyoung K. Yi, Katarina Kraljic

Abstract: We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes… ▽ More We explore how observations relate to the physical properties of the emitting galaxies by post-processing a pair of merging $z\sim2$ galaxies from the cosmological, hydrodynamical simulation NewHorizon using LCARS (Light from Cloudy Added to RAMSES) to encode the physical properties of the simulated galaxy into H$α$ emission line. By carrying out mock observations and analysis on these data cubes we ascertain which physical properties of the galaxy will be recoverable with the HARMONI spectrograph on the European Extremely Large Telescope (ELT). We are able to estimate the galaxy's star formation rate and dynamical mass to a reasonable degree of accuracy, with values within a factor of $1.81$ and $1.38$ of the true value. The kinematic structure of the galaxy is also recovered in mock observations. Furthermore, we are able to recover radial profiles of the velocity dispersion and are therefore able to calculate how the dynamical ratio varies as a function of distance from the galaxy centre. Finally, we show that when calculated on galaxy scales the dynamical ratio does not always provide a reliable measure of a galaxy's stability against gravity or act as an indicator of a minor merger. △ Less

Submitted 25 April, 2022; originally announced April 2022.

Comments: 20 pages, 17 figures, Accepted for Publication in MNRAS, Comments and questions welcome

Massive central galaxies of galaxy groups in the Romulus simulations: an overview of galaxy properties at z=0

Authors: Seoyoung Lyla Jung, Douglas Rennehan, Vida Saeedzadeh, Arif Babul, Michael Tremmel, Thomas R. Quinn, S. Ilani Loubser, E. O'Sullivan, Sukyoung K. Yi

Abstract: Contrary to many stereotypes about massive galaxies, observed brightest group galaxies (BGGs) are diverse in their star formation rates, kinematic properties, and morphologies. Studying how they evolve into and express such diverse characteristics is an important piece of the galaxy formation puzzle. We use a high-resolution cosmological suite of simulations Romulus and compare simulated central g… ▽ More Contrary to many stereotypes about massive galaxies, observed brightest group galaxies (BGGs) are diverse in their star formation rates, kinematic properties, and morphologies. Studying how they evolve into and express such diverse characteristics is an important piece of the galaxy formation puzzle. We use a high-resolution cosmological suite of simulations Romulus and compare simulated central galaxies in group-scale halos at $z=0$ to observed BGGs. The comparison encompasses the stellar mass-halo mass relation, various kinematic properties and scaling relations, morphologies, and the star formation rates. Generally, we find that Romulus reproduces the full spectrum of diversity in the properties of the BGGs very well, albeit with a tendency toward lower than the observed fraction of quenched BGGs. We find both early-type S0 and elliptical galaxies as well as late-type disk galaxies; we find Romulus galaxies that are fast-rotators as well as slow-rotators; and we observe galaxies transforming from late-type to early-type following strong dynamical interactions with satellites. We also carry out case studies of selected Romulus galaxies to explore the link between their properties, and the recent evolution of the stellar system as well as the surrounding intragroup/circumgalactic medium. In general, mergers/strong interactions quench star-forming activity and disrupt the stellar disk structure. Sometimes, however, such interactions can also trigger star-formation and galaxy rejuvenation. Black hole feedback can also lead to a decline of the star formation rate but by itself, it does not typically lead to complete quenching of the star formation activity in the BGGs. △ Less

Submitted 10 June, 2022; v1 submitted 28 February, 2022; originally announced March 2022.

Comments: 27 pages, 18 figures, Accepted to MNRAS

The SAMI Galaxy Survey: the difference between ionised gas and stellar velocity dispersions

Authors: Sree Oh, Matthew Colless, Francesco D'Eugenio, Scott M. Croom, Luca Cortese, Brent Groves, Lisa J. Kewley, Jesse van de Sande, Henry Zovaro, Mathew R. Varidel, Stefania Barsanti, Sarah Brough, Julia J. Bryant, Sarah Casura, Jon S. Lawrence, Nuria P. F. Lorente, Anne M. Medling, Matt S. Owers, Sukyoung K. Yi

Abstract: We investigate the mean locally-measured velocity dispersions of ionised gas ($σ_{\rm gas}$) and stars ($σ_*$) for 1090 galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9.5$ from the SAMI Galaxy Survey. For star-forming galaxies, $σ_*$ tends to be larger than $σ_{\rm gas}$, suggesting that stars are in general dynamically hotter than the ionised gas (asymmetric drift). The difference betwe… ▽ More We investigate the mean locally-measured velocity dispersions of ionised gas ($σ_{\rm gas}$) and stars ($σ_*$) for 1090 galaxies with stellar masses $\log\,(M_*/M_{\odot}) \geq 9.5$ from the SAMI Galaxy Survey. For star-forming galaxies, $σ_*$ tends to be larger than $σ_{\rm gas}$, suggesting that stars are in general dynamically hotter than the ionised gas (asymmetric drift). The difference between $σ_{\rm gas}$ and $σ_*$ ($Δσ$) correlates with various galaxy properties. We establish that the strongest correlation of $Δσ$ is with beam smearing, which inflates $σ_{\rm gas}$ more than $σ_*$, introducing a dependence of $Δσ$ on both the effective radius relative to the point spread function and velocity gradients. The second-strongest correlation is with the contribution of active galactic nuclei (AGN) (or evolved stars) to the ionised gas emission, implying the gas velocity dispersion is strongly affected by the power source. In contrast, using the velocity dispersion measured from integrated spectra ($σ_{\rm aper}$) results in less correlation between the aperture-based $Δσ$ ($Δσ_{\rm aper}$) and the power source. This suggests that the AGN (or old stars) dynamically heat the gas without causing significant deviations from dynamical equilibrium. Although the variation of $Δσ_{\rm aper}$ is much smaller than that of $Δσ$, a correlation between $Δσ_{\rm aper}$ and gas velocity gradient is still detected, implying there is a small bias in dynamical masses derived from stellar and ionised gas velocity dispersions. △ Less

Submitted 21 February, 2022; originally announced February 2022.

Comments: 18 pages, 13 figures

Performance Enhancement of Tree-based Friends-of-friend Galaxy-finder for High-resolution Simulations of Galaxy Formation

Authors: Jinsu Rhee, Pascal Elahi, Sukyoung K. Yi

Abstract: 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-friend (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 beca… ▽ More 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-friend (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 6-dimensional (6D) FoF galaxy finder code, VELOCIraptor (Elahi et al. 2019), on the NewHorizon simulation (Dubois et al. 2021). The high particle resolution of NewHorizon ($M_{\rm star} \sim 10^4 M_{\odot}$) allows a large central number density ($10^{6}\,{\rm 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, re-ordering 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 6D FoF search. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: 24 pages, 16 figures, and 2 tables. Accepted for publication in ApJ

The NewHorizon Simulation -- To Bar Or Not To Bar

Authors: J. Reddish, K. Kraljic, M. S. Petersen, K. Tep, Y. Dubois, C. Pichon, S. Peirani, F. Bournaud, H. Choi, J. Devriendt, R. Jackson, G. Martin, M. J. Park, M. Volonteri, S. K. Yi

Abstract: We use the NewHorizon simulation to study the redshift evolution of bar properties and fractions within galaxies in the stellar masses range $M_{\star} = 10^{7.25} - 10^{11.4} \ \rm{M}_{\odot}$ over the redshift range $z = 0.25 - 1.3$. We select disc galaxies using stellar kinematics as a proxy for galaxy morphology. We employ two different automated bar detection methods, coupled with visual insp… ▽ More We use the NewHorizon simulation to study the redshift evolution of bar properties and fractions within galaxies in the stellar masses range $M_{\star} = 10^{7.25} - 10^{11.4} \ \rm{M}_{\odot}$ over the redshift range $z = 0.25 - 1.3$. We select disc galaxies using stellar kinematics as a proxy for galaxy morphology. We employ two different automated bar detection methods, coupled with visual inspection, resulting in observable bar fractions of $f_{\rm bar} = 0.070_{-0.012}^{+0.018}$ at $z\sim$ 1.3, decreasing to $f_{\rm bar} = 0.011_{-0.003}^{+0.014}$ at $z\sim$ 0.25. Only one galaxy is visually confirmed as strongly barred in our sample. This bar is hosted by the most massive disk and only survives from $z=1.3$ down to $z=0.7$. Such a low bar fraction, in particular amongst Milky Way-like progenitors, highlights a missing bars problem, shared by literally all cosmological simulations with spatial resolution $<$100 pc to date. The analysis of linear growth rates, rotation curves and derived summary statistics of the stellar, gas and dark matter components suggest that galaxies with stellar masses below $10^{9.5}-10^{10} \ \rm{M}_{\odot}$ in NewHorizon appear to be too dominated by dark matter to form a bar, while more massive galaxies typically have formed large bulges that prevent bar persistence at low redshift. This investigation confirms that the evolution of the bar fraction puts stringent constraints on the assembly history of baryons and dark matter onto galaxies. △ Less

Submitted 24 February, 2022; v1 submitted 4 June, 2021; originally announced June 2021.

Comments: 25 pages, 13 figures, 5 tables, accepted by Monthly Notices of Royal Astronomical Society (MNRAS) on 16/02/22

Predicting the Observability of Population III Stars with ELT-HARMONI via the Helium $1640{\rmÅ}$ emission line

Authors: Kearn Grisdale, Niranjan Thatte, Julien Devriendt, Miguel Pereira-Santaella, Adrianne Slyz, Taysun Kimm, Yohan Dubois, Sukyoung K. Yi

Abstract: Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-process… ▽ More Population III (Pop. III) stars, as of yet, have not been detected, however as we move into the era of extremely large telescopes this is likely to change. One likely tracer for Pop. III stars is the HeII$\lambda1640$ emission line, which will be detectable by the HARMONI spectrograph on the European Extremely Large Telescope (ELT) over a broad range of redshifts ($2\leq z\leq14$). By post-processing galaxies from the cosmological, AMR-hydrodynamical simulation NewHorizon with theoretical spectral energy distributions (SED) for Pop. III stars and radiative transfer (i.e. the Yggdrasil Models and CLOUDY look-up tables respectively) we are able to compute the flux of HeII$\lambda1640$ for individual galaxies. From mock 10 hour observations of these galaxies we show that HARMONI will be able to detect Pop. III stars in galaxies up to $z\sim10$ provided Pop. III stars have a top heavy Initial Mass Function (IMF). Furthermore, we find that should Pop. III stars instead have an IMF similar to those of the Pop. I stars, the HeII$\lambda1640$ line would only be observable for galaxies with Pop. III stellar masses in excess of $10^{7}{\,M}_\odot\,{\rm yr}^{-1}$, average stellar age $<1{\rm Myr}$ at $z=4$. Finally, we are able to determine the minimal intrinsic flux required for HARMONI to detect Pop. III stars in a galaxy up to $z=10$. △ Less

Submitted 5 January, 2021; originally announced January 2021.

Comments: 23 Pages, 17 Figures, Accepted for publication in Monthly Notices of the Royal Astronomical Society

Star-Gas Misalignment in Galaxies: II. Origins Found from the Horizon-AGN Simulation

Authors: Donghyeon J. Khim, Sukyoung K. Yi, Christophe Pichon, Yohan Dubois, Julien Devriendt, Hoseung Choi, Julia J. Bryant, Scott M. Croom

Abstract: There have been many studies aiming to reveal the origins of the star-gas misalignment found in galaxies, but there still is a lack of understanding of the contribution from each formation channel candidate. We aim to answer the question by investigating the misaligned galaxies in Horizon-AGN, a cosmological large-volume simulation of galaxy formation. There are 27,903 galaxies of stellar mass… ▽ More There have been many studies aiming to reveal the origins of the star-gas misalignment found in galaxies, but there still is a lack of understanding of the contribution from each formation channel candidate. We aim to answer the question by investigating the misaligned galaxies in Horizon-AGN, a cosmological large-volume simulation of galaxy formation. There are 27,903 galaxies of stellar mass $M_* > 10^{10} M_\odot$ in our sample, of which 5,984 are in a group of the halo mass of $M_{200} > 10^{12} M_\odot$. We have identified four main formation channels of misalignment and quantified their level of contribution: mergers (35%), interaction with nearby galaxies (23%), interaction with dense environments or their central galaxies (21%), and secular evolution including smooth accretion from neighboring filaments (21%). We found in the simulation that the gas, rather than stars, is typically more vulnerable to dynamical disturbances; hence, misalignment formation is mainly due to the change in the rotational axis of the gas rather than stars, regardless of the origin. We have also inspected the lifetime (duration) of the misalignment. The decay timescale of the misalignment shows a strong anti-correlation with the kinematic morphology ($V/σ$) and the cold gas fraction of the galaxy. The misalignment has a longer lifetime in denser regions, which is linked with the environmental impact on the host galaxy. There is a substantial difference in the length of the misalignment lifetime depending on the origin, and it can be explained by the magnitude of the initial position angle offset and the physical properties of the galaxies. △ Less

Submitted 8 December, 2020; originally announced December 2020.

Comments: 27 pages, 20 figures, 1 table. Submitted to ApJ Supplement. Comments are welcome

The MAGPI Survey -- science goals, design, observing strategy, early results and theoretical framework

Authors: C. Foster, J. T. Mendel, C. D. P. Lagos, E. Wisnioski, T. Yuan, F. D'Eugenio, T. M. Barone, K. E. Harborne, S. P. Vaughan, F. Schulze, R. -S. Remus, A. Gupta, F. Collacchioni, D. J. Khim, P. Taylor, R. Bassett, S. M. Croom, R. M. McDermid, A. Poci, A. J. Battisti, J. Bland-Hawthorn, S. Bellstedt, M. Colless, L. J. M. Davies, C. Derkenne , et al. (18 additional authors not shown)

Abstract: We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive… ▽ More We present an overview of the Middle Ages Galaxy Properties with Integral Field Spectroscopy (MAGPI) survey, a Large Program on ESO/VLT. MAGPI is designed to study the physical drivers of galaxy transformation at a lookback time of 3-4 Gyr, during which the dynamical, morphological, and chemical properties of galaxies are predicted to evolve significantly. The survey uses new medium-deep adaptive optics aided MUSE observations of fields selected from the GAMA survey, providing a wealth of publicly available ancillary multi-wavelength data. With these data, MAGPI will map the kinematic and chemical properties of stars and ionised gas for a sample of 60 massive (> 7 x 10^10 M_Sun) central galaxies at 0.25 < z <0.35 in a representative range of environments (isolated, groups and clusters). The spatial resolution delivered by MUSE with Ground Layer Adaptive Optics (GLAO, 0.6-0.8 arcsec FWHM) will facilitate a direct comparison with Integral Field Spectroscopy surveys of the nearby Universe, such as SAMI and MaNGA, and at higher redshifts using adaptive optics, e.g. SINS. In addition to the primary (central) galaxy sample, MAGPI will deliver resolved and unresolved spectra for as many as 150 satellite galaxies at 0.25 < z <0.35, as well as hundreds of emission-line sources at z < 6. This paper outlines the science goals, survey design, and observing strategy of MAGPI. We also present a first look at the MAGPI data, and the theoretical framework to which MAGPI data will be compared using the current generation of cosmological hydrodynamical simulations including EAGLE, Magneticum, HORIZON-AGN, and Illustris-TNG. Our results show that cosmological hydrodynamical simulations make discrepant predictions in the spatially resolved properties of galaxies at z ~ 0.3. MAGPI observations will place new constraints and allow for tangible improvements in galaxy formation theory. △ Less

Submitted 14 June, 2021; v1 submitted 27 November, 2020; originally announced November 2020.

Comments: 31 pages, 13 figures, PASA accepted

YZiCS: On the Mass Segregation of Galaxies in Clusters

Authors: Seonwoo Kim, Emanuele Contini, Hoseung Choi, San Han, Jaehyun Lee, Sree Oh, Xi Kang, Sukyoung K. Yi

Abstract: Mass segregation, a tendency of more massive galaxies being distributed closer to the cluster center, is naturally expected from dynamical friction, but its presence is still controversial. Using deep optical observations of 14 Abell clusters (KYDISC) and a set of hydrodynamic simulations (YZiCS), we find in some cases a hint of mass segregation inside the virial radius. Segregation is visible mor… ▽ More Mass segregation, a tendency of more massive galaxies being distributed closer to the cluster center, is naturally expected from dynamical friction, but its presence is still controversial. Using deep optical observations of 14 Abell clusters (KYDISC) and a set of hydrodynamic simulations (YZiCS), we find in some cases a hint of mass segregation inside the virial radius. Segregation is visible more clearly when the massive galaxy fraction is used instead of mean stellar mass. The trend is more significant in the simulations than in the observations. To find out the mechanisms working on mass segregation, we look into the evolution of individual clusters simulated. We find that the degree of mass segregation is different for different clusters: the trend is visible only for low-mass clusters. We compare the masses of galaxies and their dark haloes at the time of infall and at the present epoch to quantify the amount of tidal stripping. We then conclude that satellites that get accreted at earlier epochs or galaxies in more massive clusters go through more tidal stripping. These effects in combination result in a correlation between the host halo mass and the degree of stellar mass segregation. △ Less

Submitted 11 October, 2020; originally announced October 2020.

Comments: 10 pages, 8 figures, 2 tables. Accepted for publication in ApJ

Dark-matter-deficient dwarf galaxies form via tidal stripping of dark matter in interactions with massive companions

Authors: R. A. Jackson, S. Kaviraj, G. Martin, J. E. G. Devriendt, A. Slyz, J. Silk, Y. Dubois, S. K. Yi, C. Pichon, M. Volonteri, H. Choi, T. Kimm, K. Kraljic, S. Peirani

Abstract: In the standard Lambda-CDM paradigm, dwarf galaxies are expected to be dark-matter-rich, as baryonic feedback is thought to quickly drive gas out of their shallow potential wells and quench star formation at early epochs. Recent observations of local dwarfs with extremely low dark matter content appear to contradict this picture, potentially bringing the validity of the standard model into questio… ▽ More In the standard Lambda-CDM paradigm, dwarf galaxies are expected to be dark-matter-rich, as baryonic feedback is thought to quickly drive gas out of their shallow potential wells and quench star formation at early epochs. Recent observations of local dwarfs with extremely low dark matter content appear to contradict this picture, potentially bringing the validity of the standard model into question. We use NewHorizon, a high-resolution cosmological simulation, to demonstrate that sustained stripping of dark matter, in tidal interactions between a massive galaxy and a dwarf satellite, naturally produces dwarfs that are dark-matter-deficient, even though their initial dark matter fractions are normal. The process of dark matter stripping is responsible for the large scatter in the halo-to-stellar mass relation in the dwarf regime. The degree of stripping is driven by the closeness of the orbit of the dwarf around its massive companion and, in extreme cases, produces dwarfs with halo-to-stellar mass ratios as low as unity, consistent with the findings of recent observational studies. ~30 per cent of dwarfs show some deviation from normal dark matter fractions due to dark matter stripping, with 10 per cent showing high levels of dark matter deficiency (Mhalo/M*<10). Given their close orbits, a significant fraction of dark-matter-deficient dwarfs merge with their massive companions (e.g. ~70 per cent merge over timescales of ~3.5 Gyrs), with the dark-matter-deficient population being constantly replenished by new interactions between dwarfs and massive companions. The creation of these galaxies is, therefore, a natural by-product of galaxy evolution and their existence is not in tension with the standard paradigm. △ Less

Submitted 7 January, 2021; v1 submitted 5 October, 2020; originally announced October 2020.

Comments: Accepted for publication in MNRAS, 13 pages, 11 figures

Exploring the origin of thick disks using the NewHorizon and Galactica simulations

Authors: Minjung J. Park, Sukyoung K. Yi, Sebastien Peirani, Christophe Pichon, Yohan Dubois, Hoseung Choi, Julien Devriendt, Sugata Kaviraj, Taysun Kimm, Katarina Kraljic, Marta Volonteri

Abstract: Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than $10^{10}\,\rm M_\odot$ in recent cosmological high-resolution zoom-in simulations: Galactica and NewHorizon. The thin and thick disks are reasonably reproduced by… ▽ More Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than $10^{10}\,\rm M_\odot$ in recent cosmological high-resolution zoom-in simulations: Galactica and NewHorizon. The thin and thick disks are reasonably reproduced by the simulations with scale heights and luminosity ratios as observed. We then spatially classify the thin and thick disks and find that the thick disk stars are older, metal-poorer, kinematically-hotter, and higher in accreted star fraction, while both disks are dominated by the stars formed in situ. Half of the in-situ stars in the thick disks are formed before the galaxies develop their disks, and the rest are formed in spatially and kinematically thinner disks and then thickened with time by heating. However, the 19 galaxies have various properties and evolutionary routes, highlighting the need for statistically-large samples to draw general conclusions. We conclude from our simulations that the thin and thick disk components are not entirely distinct in terms of formation processes, but rather markers of the evolution of galactic disks. Moreover, as the combined result of the thickening of the existing disk stars and the continued formation of young thin-disk stars, the vertical distribution of stars does not change much after the disks settle, pointing to the modulation of both orbital diffusion and star formation by the same confounding factor: the proximity of galaxies to marginal stability. △ Less

Submitted 25 February, 2021; v1 submitted 25 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJS

Introducing the NewHorizon simulation: Galaxy properties with resolved internal dynamics across cosmic time

Authors: Yohan Dubois, Ricarda Beckmann, Frédéric Bournaud, Hoseung Choi, Julien Devriendt, Ryan Jackson, Sugata Kaviraj, Taysun Kimm, Katarina Kraljic, Clotilde Laigle, Garreth Martin, Min-Jung Park, Sébastien Peirani, Christophe Pichon, Marta Volonteri, Sukyoung K. Yi

Abstract: Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NewHorizon project in which we simulate at high resolution a zoom-in region of… ▽ More Hydrodynamical cosmological simulations are increasing their level of realism by considering more physical processes and having greater resolution or larger statistics. However, usually either the statistical power of such simulations or the resolution reached within galaxies are sacrificed. Here, we introduce the NewHorizon project in which we simulate at high resolution a zoom-in region of $\sim(16\,\rm Mpc)^3$ that is larger than a standard zoom-in region around a single halo and is embedded in a larger box. A resolution of up to 34 pc is reached within galaxies; this allows the simulation to capture the multi-phase nature of the interstellar medium and the clumpy nature of the star formation process in galaxies. In this introductory paper, we present several key fundamental properties of galaxies and their black holes, including the galaxy mass function, cosmic star formation rate, galactic metallicities, the Kennicutt-Schmidt relation, the stellar-to-halo mass relation, galaxy sizes, stellar kinematics and morphology, gas content within galaxies and its kinematics, and the black hole mass and spin properties over time. The various scaling relations are broadly reproduced by NewHorizon with some differences with the standard observables. Owing to its exquisite spatial resolution, NewHorizon captures the inefficient process of star formation in galaxies, which evolve over time from being more turbulent, gas rich, and star bursting at high redshift. These high-redshift galaxies are also more compact, and they are more elliptical and clumpier until the level of internal gas turbulence decays enough to allow for the formation of discs. The NewHorizon simulation gives access to a broad range of galaxy formation and evolution physics at low-to-intermediate stellar masses, which is a regime that will become accessible in the near future through surveys such as the LSST. △ Less

Submitted 28 June, 2021; v1 submitted 22 September, 2020; originally announced September 2020.

Comments: 29 pages, 28 figures, A&A accepted

The role of mergers and interactions in driving the evolution of dwarf galaxies over cosmic time

Authors: G. Martin, R. A. Jackson, S. Kaviraj, H. Choi, J. E. G. Devriendt, Y. Dubois, T. Kimm, K. Kraljic, S. Peirani, C. Pichon, M. Volonteri, S. K. Yi

Abstract: Dwarf galaxies (M*<10^9 Msun) are key drivers of mass assembly in high mass galaxies, but relatively little is understood about the assembly of dwarf galaxies themselves. Using the \textsc{NewHorizon} cosmological simulation (40 pc spatial resolution), we investigate how mergers and fly-bys drive the mass assembly and structural evolution of around 1000 field and group dwarfs up to z=0.5. We find… ▽ More Dwarf galaxies (M*<10^9 Msun) are key drivers of mass assembly in high mass galaxies, but relatively little is understood about the assembly of dwarf galaxies themselves. Using the \textsc{NewHorizon} cosmological simulation (40 pc spatial resolution), we investigate how mergers and fly-bys drive the mass assembly and structural evolution of around 1000 field and group dwarfs up to z=0.5. We find that, while dwarf galaxies often exhibit disturbed morphologies (5 and 20 per cent are disturbed at z=1 and z=3 respectively), only a small proportion of the morphological disturbances seen in dwarf galaxies are driven by mergers at any redshift (for 10^9 Msun, mergers drive only 20 per cent morphological disturbances). They are instead primarily the result of interactions that do not end in a merger (e.g. fly-bys). Given the large fraction of apparently morphologically disturbed dwarf galaxies which are not, in fact, merging, this finding is particularly important to future studies identifying dwarf mergers and post-mergers morphologically at intermediate and high redshifts. Dwarfs typically undergo one major and one minor merger between z=5 and z=0.5, accounting for 10 per cent of their total stellar mass. Mergers can also drive moderate star formation enhancements at lower redshifts (3 or 4 times at z=1), but this accounts for only a few per cent of stellar mass in the dwarf regime given their infrequency. Non-merger interactions drive significantly smaller star formation enhancements (around two times), but their preponderance relative to mergers means they account for around 10 per cent of stellar mass formed in the dwarf regime. △ Less

Submitted 1 November, 2020; v1 submitted 15 July, 2020; originally announced July 2020.

Comments: 22 pages, 19 figures, version accepted for publication in MNRAS. Some minor changes and additional appendices

The origin of low-surface-brightness galaxies in the dwarf regime

Authors: R. A. Jackson, G. Martin, S. Kaviraj, M. Ramsøy, J. E. G. Devriendt, T. Sedgwick, C. Laigle, H. Choi, R. S. Beckmann, M. Volonteri, Y. Dubois, C. Pichon, S. K. Yi, A. Slyz, K. Kraljic, T. Kimm, S. Peirani, I. Baldry

Abstract: Low-surface-brightness galaxies (LSBGs) -- defined as systems that are fainter than the surface-brightness limits of past wide-area surveys -- form the overwhelming majority of galaxies in the dwarf regime (M* < 10^9 MSun). Using NewHorizon, a high-resolution cosmological simulation, we study the origin of LSBGs and explain why LSBGs at similar stellar mass show the large observed spread in surfac… ▽ More Low-surface-brightness galaxies (LSBGs) -- defined as systems that are fainter than the surface-brightness limits of past wide-area surveys -- form the overwhelming majority of galaxies in the dwarf regime (M* < 10^9 MSun). Using NewHorizon, a high-resolution cosmological simulation, we study the origin of LSBGs and explain why LSBGs at similar stellar mass show the large observed spread in surface brightness. New Horizon galaxies populate a well-defined locus in the surface brightness -- stellar mass plane, with a spread of ~3 mag arcsec^-2, in agreement with deep SDSS Stripe data. Galaxies with fainter surface brightnesses today are born in regions of higher dark-matter density. This results in faster gas accretion and more intense star formation at early epochs. The stronger resultant supernova feedback flattens gas profiles at a faster rate which, in turn, creates shallower stellar profiles (i.e. more diffuse systems) more rapidly. As star formation declines towards late epochs (z<1), the larger tidal perturbations and ram pressure experienced by these systems (due to their denser local environments) accelerate the divergence in surface brightness, by increasing their effective radii and reducing star formation respectively. A small minority of dwarfs depart from the main locus towards high surface brightnesses, making them detectable in past wide surveys. These systems have anomalously high star-formation rates, triggered by recent, fly-by or merger-driven starbursts. We note that objects considered extreme/anomalous at the depth of current datasets, e.g. `ultra-diffuse galaxies', actually dominate the predicted dwarf population and will be routinely visible in future surveys like LSST. △ Less

Submitted 6 January, 2021; v1 submitted 13 July, 2020; originally announced July 2020.

Comments: 16 pages, 15 figures, Accepted for publication in MNRAS

Black hole mergers from dwarf to massive galaxies with the NewHorizon and Horizon-AGN simulations

Authors: Marta Volonteri, Hugo Pfister, Ricarda S. Beckman, Yohan Dubois, Monica Colpi, Christopher J. Conselice, Massimo Dotti, Garreth Martin, Ryan Jackson, Katarina Kraljic, Christophe Pichon, Maxime Trebitsch, Sukyoung. K. Yi, Julien Devriendt, Sebastien Peirani

Abstract: Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large… ▽ More Massive black hole (MBH) coalescences are powerful sources of low-frequency gravitational waves. To study these events in the cosmological context we need to trace the large-scale structure and cosmic evolution of a statistical population of galaxies, from dim dwarfs to bright galaxies. To cover such a large range of galaxy masses, we analyse two complementary simulations: Horizon-AGN with a large volume and low resolution which tracks the high-mass (> 1e7 Msun) MBH population, and NewHorizon with a smaller volume but higher resolution that traces the low-mass (< 1e7 Msun) MBH population. While Horizon-AGN can be used to estimate the rate of inspirals for Pulsar Timing Arrays, NewHorizon can investigate MBH mergers in a statistical sample of dwarf galaxies for LISA, which is sensitive to low-mass MBHs. We use the same method to analyse the two simulations, post-processing MBH dynamics to account for time delays mostly determined by dynamical friction and stellar hardening. In both simulations, MBHs typically merge long after the galaxies do, so that the galaxy morphology at the time of the MBH merger is no longer determined by the galaxy merger from which the MBH merger originated. These time delays cause a loss of high-z MBH coalescences, shifting the peak of the MBH merger rate to z~1-2. This study shows how tracking MBH mergers in low-mass galaxies is crucial to probing the MBH merger rate for LISA and investigate the properties of the host galaxies. △ Less

Submitted 15 September, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: MNRAS in press

YZiCS: Unveiling Quenching History of Cluster Galaxies Using Phase-space Analysis

Authors: Jinsu Rhee, Rory Smith, Hoseung Choi, Emanuele Contini, S. Lyla Jung, San Han, Sukyoung K. Yi

Abstract: We used the time since infall (TSI) of galaxies, obtained from the Yonsei Zoom-in Cluster Simulation, and the star formation rate (SFR) from the Sloan Digital Sky Survey (SDSS) Data Release 10 to study how quickly star formation of disk galaxies is quenched in cluster environments. We first confirm that both simulated and observed galaxies are consistently distributed in phase space. We then hypot… ▽ More We used the time since infall (TSI) of galaxies, obtained from the Yonsei Zoom-in Cluster Simulation, and the star formation rate (SFR) from the Sloan Digital Sky Survey (SDSS) Data Release 10 to study how quickly star formation of disk galaxies is quenched in cluster environments. We first confirm that both simulated and observed galaxies are consistently distributed in phase space. We then hypothesize that the TSI and SFR are causally connected; thus, both the TSI and SFR of galaxies at each position of phase space can be associated through abundance matching. Using a flexible model, we derive the star formation history (SFH) of cluster galaxies that best reproduces the relationship between the TSI and SFR at $z\sim 0.08$. According to this SFH, we find that the galaxies with $M_{*} > 10^{9.5} M_{\odot}$ generally follow the so-called "delayed-then-rapid" quenching pattern. Our main results are as following: (i) Part of the quenching takes place outside clusters through mass quenching and pre-processing. The e-folding timescale of this "$ex\text{-}situ$ quenching phase" is roughly 3 Gyr with a strong inverse mass dependence. (ii) The pace of quenching is maintained roughly for 2 Gyr ("delay time") during the first crossing time into the cluster. During the delay time, quenching remains gentle probably because gas loss happens primarily on hot and neutral gases. (iii) Quenching becomes more dramatic (e-folding timescale of roughly 1 Gyr) after delay time, probably because ram pressure stripping is strongest near the cluster center. Counter-intuitively, more massive galaxies show shorter quenching timescales mainly because they enter their clusters with lower gas fractions due to $ex\text{-}situ$ quenching. △ Less

Submitted 11 February, 2020; originally announced February 2020.

Comments: 24 pages, 11 figures, 1 table, accepted to ApJS

The Roles of Mass and Environment in the Quenching of Galaxies. II

Authors: E. Contini, Q. Gu, X. Ge, J. Rhee, S. K. Yi, X. Kang

Abstract: We take advantage of an analytic model of galaxy formation coupled to the merger tree of an N-body simulation to study the roles of environment and stellar mass in the quenching of galaxies. The model has been originally set in order to provide the observed evolution of the stellar mass function as well as reasonable predictions of the star formation rate-stellar mass relation, from high redshift… ▽ More We take advantage of an analytic model of galaxy formation coupled to the merger tree of an N-body simulation to study the roles of environment and stellar mass in the quenching of galaxies. The model has been originally set in order to provide the observed evolution of the stellar mass function as well as reasonable predictions of the star formation rate-stellar mass relation, from high redshift to the present time. We analyse the stellar mass and environmental quenching efficiencies and their dependence on stellar mass, halo mass (taken as a proxy for the environment) and redshift. Our analysis shows that the two quenching efficiencies are redshift, stellar and halo mass dependent, and that the halo mass is also a good proxy for the environment. The environmental quenching increases with decreasing redshift and is inefficient below $\log M_* \sim 9.5$, reaches the maximum value at $\log M_* \sim 10.5$, and decreases again, becoming poorly efficient at very high stellar mass ($\log M_* \gtrsim 11.5$). Central and satellites galaxies are mass quenched differently: for the former, the quenching efficiency depends very weakly on redshift, but strongly on stellar mass; for the latter, it strongly depends on both stellar mass and redshift in the range $10\lesssim \log M_* \lesssim 11$. According to the most recent observational results, we find that the two quenching efficiencies are not separable: intermediate mass galaxies are environmental quenched faster, as well as intermediate/massive galaxies in more massive haloes. At stellar masses lower than $\log M_* \lesssim 9.5$ both quenching mechanisms become inefficient, independently of the redshift. △ Less

Submitted 22 January, 2020; v1 submitted 5 January, 2020; originally announced January 2020.

Comments: 13 pages, 8 figures, accepted for publication in ApJ, minor revision after Proof correction

Star-gas misalignment in galaxies: I. The properties of galaxies from the Horizon-AGN simulation and comparisons to SAMI

Authors: Donghyeon J. Khim, Sukyoung K. Yi, Yohan Dubois, Julia J. Bryant, Christophe Pichon, Scott M. Croom, Joss Bland-Hawthorn, Sarah Brough, Hoseung Choi, Julien Devriendt, Brent Groves, Matt S. Owers, Samuel N. Richards, Jesse van de Sande, Sarah M. Sweet

Abstract: Recent integral field spectroscopy observations have found that about 11% of galaxies show star-gas misalignment. The misalignment possibly results from external effects such as gas accretion, interaction with other objects, and other environmental effects, hence providing clues to these effects. We explore the properties of misaligned galaxies using Horizon-AGN, a large-volume cosmological simula… ▽ More Recent integral field spectroscopy observations have found that about 11% of galaxies show star-gas misalignment. The misalignment possibly results from external effects such as gas accretion, interaction with other objects, and other environmental effects, hence providing clues to these effects. We explore the properties of misaligned galaxies using Horizon-AGN, a large-volume cosmological simulation, and compare the result with the result of the Sydney-AAO Multi-object integral field spectrograph (SAMI) Galaxy Survey. Horizon-AGN can match the overall misalignment fraction and reproduces the distribution of misalignment angles found by observations surprisingly closely. The misalignment fraction is found to be highly correlated with galaxy morphology both in observations and in the simulation: early-type galaxies are substantially more frequently misaligned than late-type galaxies. The gas fraction is another important factor associated with misalignment in the sense that misalignment increases with decreasing gas fraction. However, there is a significant discrepancy between the SAMI and Horizon-AGN data in the misalignment fraction for the galaxies in dense (cluster) environments. We discuss possible origins of misalignment and disagreement. △ Less

Submitted 8 December, 2020; v1 submitted 15 November, 2019; originally announced November 2019.

Comments: 23 pages with 15 figures. Accepted for publication in ApJ

The SAMI Galaxy Survey: rules of behaviour for spin-ellipticity radial tracks in galaxies

Authors: Alexander Rawlings, Caroline Foster, Jesse van de Sande, Dan S. Taranu, Scott M. Croom, Joss Bland-Hawthorn, Sarah Brough, Julia J. Bryant, Matthew Colless, Claudia del P. Lagos, Iraklis S. Konstantopoulos, Jon S. Lawrence, Ángel R. López-Sánchez, Nuria P. F. Lorente, Anne M. Medling, Sree Oh, Matt S. Owers, Samuel N. Richards, Nicholas Scott, Sarah M. Sweet, Sukyoung K. Yi

Abstract: We study the behaviour of the spin-ellipticity radial tracks for 507 galaxies from the Sydney AAO Multi-object Integral Field (SAMI) Galaxy Survey with stellar kinematics out to $\geq1.5R_\text{e}$. We advocate for a morpho-dynamical classification of galaxies, relying on spatially-resolved photometric and kinematic data. We find the use of spin-ellipticity radial tracks is valuable in identifying… ▽ More We study the behaviour of the spin-ellipticity radial tracks for 507 galaxies from the Sydney AAO Multi-object Integral Field (SAMI) Galaxy Survey with stellar kinematics out to $\geq1.5R_\text{e}$. We advocate for a morpho-dynamical classification of galaxies, relying on spatially-resolved photometric and kinematic data. We find the use of spin-ellipticity radial tracks is valuable in identifying substructures within a galaxy, including embedded and counter-rotating discs, that are easily missed in unilateral studies of the photometry alone. Conversely, bars are rarely apparent in the stellar kinematics but are readily identified on images. Consequently, we distinguish the spin-ellipticity radial tracks of seven morpho-dynamical types: elliptical, lenticular, early spiral, late spiral, barred spiral, embedded disc, and 2-sigma galaxies. The importance of probing beyond the inner radii of galaxies is highlighted by the characteristics of galactic features in the spin-ellipticity radial tracks present at larger radii. The density of information presented through spin-ellipticity radial tracks emphasises a clear advantage to representing galaxies as a track, rather than a single point, in spin-ellipticity parameter space. △ Less

Submitted 4 October, 2019; v1 submitted 3 October, 2019; originally announced October 2019.

Comments: 23 pages, 15 figures

The Roles of Mass and Environment in the Quenching of Galaxies

Authors: E. Contini, Q. Gu, X. Kang, J. Rhee, S. K. Yi

Abstract: We study the roles of stellar mass and environment in quenching the star formation activity of a large set of simulated galaxies by taking advantage of an analytic model coupled to the merger tree extracted from an N-body simulation. The analytic model has been set to match the evolution of the global stellar mass function since redshift $z\sim 2.3$ and give reasonable predictions of the star form… ▽ More We study the roles of stellar mass and environment in quenching the star formation activity of a large set of simulated galaxies by taking advantage of an analytic model coupled to the merger tree extracted from an N-body simulation. The analytic model has been set to match the evolution of the global stellar mass function since redshift $z\sim 2.3$ and give reasonable predictions of the star formation history of galaxies at the same time. We find that stellar mass and environment play different roles: the star formation rate/specific star formation rate-$M_*$ relations are independent of the environment (defined as the halo mass) at any redshift probed, $0<z<1.5$, for both star forming and quiescent galaxies, while the star formation rate-$M_{halo}$ relation strongly depends on stellar mass in the same redshift range, for both star forming and quiescent galaxies. Moreover, the star formation rate and the specific star formation rate are strongly dependent on stellar mass even when the distance from the cluster core is used as a proxy for the environment, rather than the halo mass. We then conclude that stellar mass is the main driver of galaxy quenching at any redshift probed in this study, not just at $z>1$ as generally claimed, while the environment has a minimal role. All the physical processes linked to the environment must act on very short timescales, such that they do not influence the star formation of active galaxies, but increase the probability of a given galaxy to become quiescent. △ Less

Submitted 21 August, 2019; v1 submitted 6 August, 2019; originally announced August 2019.

Comments: 13 pages, 8 figures, accepted for publication in ApJ

Impact of Galaxy Mergers on the Colours of Cluster Galaxies

Authors: Sree Oh, Keunho Kim, Joon Hyeop Lee, Minjin Kim, Yun-Kyeong Sheen, Jinsu Rhee, Chang H. Ree, Hyunjin Jeong, Luis C. Ho, Jaemann Kyeong, Byeong-Gon Park, Sukyoung K. Yi

Abstract: We examine the ultraviolet and optical colours of 906 cluster galaxies from the KASI-Yonsei Deep Imaging Survey of Clusters (KYDISC). Galaxies have been divided into two categories, morphologically-disturbed and undisturbed galaxies, based on the visual signatures related to recent mergers. We find that galaxies with signatures of recent mergers show significantly bluer colours than undisturbed ga… ▽ More We examine the ultraviolet and optical colours of 906 cluster galaxies from the KASI-Yonsei Deep Imaging Survey of Clusters (KYDISC). Galaxies have been divided into two categories, morphologically-disturbed and undisturbed galaxies, based on the visual signatures related to recent mergers. We find that galaxies with signatures of recent mergers show significantly bluer colours than undisturbed galaxies. Disturbed galaxies populate more on the cluster outskirts, suggesting recent accretion into the cluster environment, which implies that disturbed galaxies can be less influenced by the environmental quenching process and remain blue. However, we still detect bluer colours of disturbed galaxies in all locations (cluster core and outskirts) for the fixed morphology, which is difficult to understand just considering the difference in time since infall into a cluster. Moreover, blue disturbed galaxies show features seemingly related to recent star formation. Therefore, we suspect that mergers make disturbed galaxies keep their blue colour longer than undisturbed galaxies under the effect of the environmental quenching through either merger-induced star formation or central gas concentration which is less vulnerable for gas stripping. △ Less

Submitted 21 July, 2019; v1 submitted 10 July, 2019; originally announced July 2019.

Comments: 13 pages, 15 figures; Accepted for publication in MNRAS

The SAMI Galaxy Survey: mass-kinematics scaling relations

Authors: Dilyar Barat, Francesco D'Eugenio, Matthew Colless, Sarah Brough, Barbara Catinella, Luca Cortese, Scott M. Croom, Anne M. Medling, Sree Oh, Jesse van de Sande, Sarah M. Sweet, Sukyoung K. Yi, Joss Bland-Hawthorn, Julia Bryant, Michael Goodwin, Brent Groves, Jon Lawrence, Matt S. Owers, Samuel N. Richards, Nicholas Scott

Abstract: We use data from the Sydney-AAO Multi-object Integral-field spectroscopy (SAMI) Galaxy Survey to study the dynamical scaling relation between galaxy stellar mass $M_*$ and the general kinematic parameter $S_K = \sqrt{K V_{rot}^2 + σ^2}$ that combines rotation velocity $V_{rot}$ and velocity dispersion $σ$. We show that the $\log M_* - \log S_K$ relation: (1)~is linear above limits set by propertie… ▽ More We use data from the Sydney-AAO Multi-object Integral-field spectroscopy (SAMI) Galaxy Survey to study the dynamical scaling relation between galaxy stellar mass $M_*$ and the general kinematic parameter $S_K = \sqrt{K V_{rot}^2 + σ^2}$ that combines rotation velocity $V_{rot}$ and velocity dispersion $σ$. We show that the $\log M_* - \log S_K$ relation: (1)~is linear above limits set by properties of the samples and observations; (2)~has slightly different slope when derived from stellar or gas kinematic measurements; (3)~applies to both early-type and late-type galaxies and has smaller scatter than either the Tully-Fisher relation ($\log M_* - \log V_{rot}$) for late types or the Faber-Jackson relation ($\log M_* - \logσ$) for early types; and (4)~has scatter that is only weakly sensitive to the value of $K$, with minimum scatter for $K$ in the range 0.4 and 0.7. We compare $S_K$ to the aperture second moment (the `aperture velocity dispersion') measured from the integrated spectrum within a 3-arcsecond radius aperture ($σ_{3^{\prime\prime}}$). We find that while $S_{K}$ and $σ_{3^{\prime\prime}}$ are in general tightly correlated, the $\log M_* - \log S_K$ relation has less scatter than the $\log M_* - \log σ_{3^{\prime\prime}}$ relation. △ Less

Submitted 29 May, 2019; originally announced May 2019.

Comments: 14 pages, 8 figures, Accepted 2019 May 22. Received 2019 May 18; in original form 2019 January 6

New Horizon: On the origin of the stellar disk and spheroid of field galaxies at $z=0.7$

Authors: Min-Jung Park, Sukyoung K. Yi, Yohan Dubois, Christophe Pichon, Taysun Kimm, Julien Devriendt, Hoseung Choi, Marta Volonteri, Sugata Kaviraj, Sebastien Peirani

Abstract: The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$ at $z=0.7$. The origins of stellar particles are classified according to their birthplac… ▽ More The origin of the disk and spheroid of galaxies has been a key open question in understanding their morphology. Using the high-resolution cosmological simulation, New Horizon, we explore kinematically decomposed disk and spheroidal components of 144 field galaxies with masses greater than $\rm 10^9\,M_{\odot}$ at $z=0.7$. The origins of stellar particles are classified according to their birthplace (in situ or ex situ) and their orbits at birth. Before disk settling, stars form mainly through chaotic mergers between proto-galaxies and become part of the spheroidal component. When disk settling starts, we find that more massive galaxies begin to form disk stars from earlier epochs; massive galaxies commence to develop their disks at $z\sim1-2$, while low-mass galaxies do after $z\sim1$. The formation of disks is affected by accretion as well, as mergers can trigger gas turbulence or induce misaligned gas infall that prevents galaxies from forming co-rotating disk stars. The importance of accreted stars is greater in more massive galaxies, especially in developing massive spheroids. A significant fraction of the spheroids comes from the disk stars that are perturbed, which becomes more important at lower redshifts. Some ($\sim12.5\%$) of our massive galaxies develop counter-rotating disks from the gas infall misaligned with the existing disk plane, which can last for more than a Gyr until they become the dominant component, and flip the angular momentum of the galaxy in the opposite direction. The final disk-to-total ratio of a galaxy needs to be understood in relation to its stellar mass and accretion history. We quantify the significance of the stars with different origins and provide them as guiding values. △ Less

Submitted 13 August, 2019; v1 submitted 6 May, 2019; originally announced May 2019.

Comments: Accepted for publication in ApJ. 23 pages with 13 figures

The SAMI Galaxy Survey: Quenching of star formation in clusters I. Transition galaxies

Authors: Matt S. Owers, Michael J. Hudson, Kyle A. Oman, Joss Bland-Hawthorn, S. Brough, Julia J. Bryant, Luca Cortese, Warrick J. Couch, Scott M. Croom, Jesse van de Sande, Christoph Federrath, Brent Groves, A. M. Hopkins, J. S. Lawrence, Nuria P. F. Lorente, Richard M. McDermid, Anne M. Medling, Samuel N. Richards, Nicholas Scott, Dan S. Taranu, Charlotte Welker, Sukyoung K. Yi

Abstract: We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}δ$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with ste… ▽ More We use integral field spectroscopy from the SAMI Galaxy Survey to identify galaxies that show evidence for recent quenching of star formation. The galaxies exhibit strong Balmer absorption in the absence of ongoing star formation in more than 10% of their spectra within the SAMI field of view. These $\rm{H}δ$-strong galaxies (HDSGs) are rare, making up only $\sim 2$% (25/1220) of galaxies with stellar mass ${\rm log(}M_*/M_{\odot})>10$. The HDSGs make up a significant fraction of non-passive cluster galaxies (15%; 17/115) and a smaller fraction (2.0%; 8/387) of the non-passive population in low-density environments. The majority (9/17) of cluster HDSGs show evidence for star formation at their centers, with the HDS regions found in the outer parts of the galaxy. Conversely, the $\rm{H}δ$-strong signal is more evenly spread across the galaxy for the majority (6/8) of HDSGs in low-density environments, and is often associated with emission lines that are not due to star formation. We investigate the location of the HDSGs in the clusters, finding that they are exclusively within 0.6$R_{200}$ of the cluster centre, and have a significantly higher velocity dispersion relative to the cluster population. Comparing their distribution in projected-phase-space to those derived from cosmological simulations indicates that the cluster HDSGs are consistent with an infalling population that have entered the central 0.5$r_{200, 3D}$ cluster region within the last $\sim 1\,$Gyr. In the 8/9 cluster HDSGs with central star formation, the extent of star formation is consistent with that expected of outside-in quenching by ram-pressure stripping. Our results indicate that the cluster HDSGs are currently being quenched by ram-pressure stripping on their first passage through the cluster. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: 44 pages (including 8 pages of appendices), 18 figures, 3 tables. Accepted for publicatioin in ApJ

Physical properties of SDSS satellite galaxies in projected phase-space

Authors: A. Pasquali, R. Smith, A. Gallazzi, G. De Lucia, S. Zibetti, M. Hirschmann, S. K. Yi

Abstract: We investigate how environment affects satellite galaxies using their location within the projected phase-space of their host haloes from the Wang et al.'s group catalogue. Using the Yonsei Zoom in Cluster Simulations, we derive zones of constant mean infall time T_inf in projected phase-space, and catalogue in which zone each observed galaxy falls. Within each zone we compute the mean observed ga… ▽ More We investigate how environment affects satellite galaxies using their location within the projected phase-space of their host haloes from the Wang et al.'s group catalogue. Using the Yonsei Zoom in Cluster Simulations, we derive zones of constant mean infall time T_inf in projected phase-space, and catalogue in which zone each observed galaxy falls. Within each zone we compute the mean observed galaxy properties including specific star formation rate, luminosity-weighted age, stellar metallicity and [alpha/Fe] abundance ratio. By comparing galaxies in different zones, we inspect how shifting the mean infall time from recent infallers (mean T_inf < 3 Gyr) to ancient infallers (mean T_inf > 5 Gyr) impacts galaxy properties at fixed stellar and halo mass. Ancient infallers are more quenched, and the impact of environmental quenching is visible down to low host masses (< group masses). Meanwhile, the quenching of recent infallers is weakly dependent on host mass, indicating they have yet to respond strongly to their current environment. [alpha/Fe] and especially metallicity are less dependent on host mass, but show a dependence on mean T_inf. We discuss these results in the context of longer exposure times for ancient infallers to environmental effects, which grow more efficient in hosts with a deeper potential well and a denser intracluster medium. We also compare our satellites with a control field sample, and find that even the most recent infallers (mean T_inf < 2 Gyr) are more quenched than field galaxies, in particular for cluster mass hosts. This supports the role of pre-processing and/or faster quenching in satellites. △ Less

Submitted 14 January, 2019; originally announced January 2019.

Comments: 24 pages, 14 figures, accepted for publication in MNRAS

The SAMI Galaxy Survey: Observing the environmental quenching of star formation in GAMA groups

Authors: A. L. Schaefer, S. M. Croom, N. Scott, S. Brough, J. T. Allen, K. Bekki, J. Bland-Hawthorn, J. V. Bloom, J. J. Bryant, L. Cortese, L. J. M. Davies, C. Federrath, L. M. R. Fogarty, A. W. Green, B. Groves, A. M. Hopkins, I. S. Konstantopoulos, A. R. López-Sánchez, J. S. Lawrence, R. E. McElroy, A. M. Medling, M. S. Owers, M. B. Pracy, S. N. Richards, A. S. G. Robotham , et al. (3 additional authors not shown)

Abstract: We explore the radial distribution of star formation in galaxies in the SAMI Galaxy Survey as a function of their local group environment. Using a sample of galaxies in groups (with halo masses less than $ \simeq 10^{14} \, \mathrm{M_{\odot}}$) from the Galaxy And Mass Assembly Survey, we find signatures of environmental quenching in high-mass groups ($M_{G} > 10^{12.5} \, \mathrm{M_{\odot}}$). Th… ▽ More We explore the radial distribution of star formation in galaxies in the SAMI Galaxy Survey as a function of their local group environment. Using a sample of galaxies in groups (with halo masses less than $ \simeq 10^{14} \, \mathrm{M_{\odot}}$) from the Galaxy And Mass Assembly Survey, we find signatures of environmental quenching in high-mass groups ($M_{G} > 10^{12.5} \, \mathrm{M_{\odot}}$). The mean integrated specific star formation rate of star-forming galaxies in high-mass groups is lower than for galaxies in low-mass groups or that are ungrouped, with $Δ\log(sSFR/\mathrm{yr^{-1}}) = 0.45 \pm 0.07$. This difference is seen at all galaxy stellar masses. In high-mass groups, star-forming galaxies more massive than $M_{*} \sim 10^{10} \, \mathrm{M_{\odot}}$ have centrally-concentrated star formation. These galaxies also lie below the star-formation main sequence, suggesting they may be undergoing outside-in quenching. Lower mass galaxies in high-mass groups do not show evidence of concentrated star formation. In groups less massive than $M_{G} = 10^{12.5} \, \mathrm{M_{\odot}}$ we do not observe these trends. In this regime we find a modest correlation between centrally-concentrated star formation and an enhancement in total star formation rate, consistent with triggered star formation in these galaxies. △ Less

Submitted 28 November, 2018; originally announced November 2018.

Comments: 22 Pages, 11 Figures, accepted for publication in MNRAS

Theoretical Predictions of Colors and Metallicity of the Intra-Cluster Light

Authors: E. Contini, S. K. Yi, X. Kang

Abstract: We study colors and metallicities of the Brightest Cluster Galaxies (BCGs) and Intra-Cluster Light (ICL) in galaxy groups and clusters, as predicted by a semi-analytic model of galaxy formation, coupled with a set of high-resolution N-body simulations. The model assumes stellar stripping and violent relaxation processes during galaxy mergers to be the main channels for the formation of the ICL. We… ▽ More We study colors and metallicities of the Brightest Cluster Galaxies (BCGs) and Intra-Cluster Light (ICL) in galaxy groups and clusters, as predicted by a semi-analytic model of galaxy formation, coupled with a set of high-resolution N-body simulations. The model assumes stellar stripping and violent relaxation processes during galaxy mergers to be the main channels for the formation of the ICL. We find that BCGs are more metal-rich and redder than the ICL, at all redshifts since the ICL starts to form ($z\sim 1$). In good agreement with several observed data, our model predicts negative radial metallicity and color gradients in the BCG+ICL system. By comparing the typical colors of the ICL with those of satellite galaxies, we find that mass and metals in the ICL come from galaxies of different mass, depending on the redshift. Stripping of low mass galaxies, $9<\log M_* <10$, is the most important contributor in the early stage of the ICL formation, but the bulk of the mass/metals contents are given by intermediate/massive galaxies, $10<\log M_* <11$, at lower redshift. Our analysis supports the idea that stellar stripping is more important than galaxy mergers in building-up the ICL, and highlights the importance of colors/metallicity measurements for understanding the formation and evolution of the ICL. △ Less

Submitted 25 November, 2018; v1 submitted 7 November, 2018; originally announced November 2018.

Comments: 12 pages, 6 figures, accepted for publication in ApJ after minor revision

YZiCS: Preprocessing of dark halos in the hydrodynamic zoom-in simulation of clusters

Authors: San Han, Rory Smith, Hoseung Choi, Luca Cortese, Barbara Catinella, Emanuele Contini, Sukyoung K. Yi

Abstract: To understand the galaxy population in clusters today, we should also consider the impact of previous environments prior to cluster infall, namely preprocessing. We use YZiCS, a hydrodynamic high-resolution zoom-in simulation of 15 clusters, and focus on the tidal stripping suffered by the dark matter halos of cluster members due to preprocessing. We find ~48% of today's cluster members were once… ▽ More To understand the galaxy population in clusters today, we should also consider the impact of previous environments prior to cluster infall, namely preprocessing. We use YZiCS, a hydrodynamic high-resolution zoom-in simulation of 15 clusters, and focus on the tidal stripping suffered by the dark matter halos of cluster members due to preprocessing. We find ~48% of today's cluster members were once satellites of other hosts. This is slightly higher than previous estimates, in part, because we consider not just group-mass hosts, but hosts of all masses also. Thus, we find the preprocessed fraction is poorly correlated with cluster mass and is instead related to each cluster's recent mass growth rate. Hosts less massive than groups are significant contributors, providing more than one-third of the total preprocessed fraction. We find that halo mass loss is a clear function of the time spent in hosts. However, two factors can increase the mass loss rate considerably; the mass ratio of a satellite to its host, and the cosmological epoch when the satellite was hosted. The latter means we may have previously underestimated the role of high redshift groups. From a sample of heavily tidally stripped members in clusters today, nearly three quarters were previously in a host. Thus, visibly disturbed cluster members are more likely to have experienced preprocessing. Being hosted before cluster infall enables cluster members to experience tidal stripping for extended durations compared to direct cluster infall and at earlier epochs when hosts were more destructive. △ Less

Submitted 8 September, 2018; originally announced September 2018.

Comments: 13 pages, 12 figures, Accepted for publication in ApJ

On the origin of gas-poor galaxies in galaxy clusters using cosmological hydrodynamic simulations

Authors: Seoyoung L. Jung, Hoseung Choi, O. Ivy Wong, Taysun Kimm, Aeree Chung, Sukyoung K. Yi

Abstract: The environmental effect is commonly used to explain the excess of gas-poor galaxies in galaxy clusters. Meanwhile, the presence of gas-poor galaxies at cluster outskirts, where galaxies have not spent enough time to feel the cluster environmental effect, hints for the presence of pre-processing. Using cosmological hydrodynamic simulations on 16 clusters, we investigate the mechanisms of gas deple… ▽ More The environmental effect is commonly used to explain the excess of gas-poor galaxies in galaxy clusters. Meanwhile, the presence of gas-poor galaxies at cluster outskirts, where galaxies have not spent enough time to feel the cluster environmental effect, hints for the presence of pre-processing. Using cosmological hydrodynamic simulations on 16 clusters, we investigate the mechanisms of gas depletion of galaxies found inside clusters. The gas depletion mechanisms can be categorized into three channels based on where and when they took place. First, 34$\%$ of our galaxies are gas poor before entering clusters (`pre-processing'). They are mainly satellites that have undergone the environmental effect inside group halos. Second, 43$\%$ of the sample became quickly gas deficient in clusters before the first pericentric pass (`fast cluster processing'). Some of them were group satellites that are low in gas at the time of cluster entry compared to the galaxies directly coming from the field. Even the galaxies with large gas fractions take this channel if they fall into massive clusters ($> 10^{14.5}\, \rm M_{\odot}$) or approach cluster centers through radial orbits. Third, 24$\%$ of our sample retain gas even after their first pericentric pass (`slow cluster processing') as they fall into the less massive clusters and/or have circular orbits. The relative importance of each channel varies with a cluster's mass, while the exact degree of significance is subject to large uncertainties. Group pre-processing accounts for a third of the total gas depletion; but it also determines the gas fraction of galaxies at their cluster entry which in turn determines whether a galaxy should take the fast or the slow cluster processing. △ Less

Submitted 5 September, 2018; originally announced September 2018.

Comments: Accepted to ApJ, 18 pages, 17 figures

Wobbling galaxy spin axes in dense environments

Authors: Jaehyun Lee, Suk Kim, Hyunjin Jeong, Rory Smith, Hoseung Choi, Ho Seong Hwang, Seok-Joo Joo, Hak-Sub Kim, Yongdae Lee, Sukyoung K. Yi

Abstract: The orientation of galaxy spin vectors within the large scale structure has been considered an important test of our understanding of structure formation. We investigate the angular changes of galaxy spin vectors in clusters - denser environments than are normally focused upon, using hydrodynamic zoomed simulations of 17 clusters YZiCS and a set of complementary controlled simulations. The magnitu… ▽ More The orientation of galaxy spin vectors within the large scale structure has been considered an important test of our understanding of structure formation. We investigate the angular changes of galaxy spin vectors in clusters - denser environments than are normally focused upon, using hydrodynamic zoomed simulations of 17 clusters YZiCS and a set of complementary controlled simulations. The magnitude by which galaxies change their spin vector is found to be a function of their rotational support with larger cumulative angular changes of spin vectors when they have initially lower $V_θ/σ$. We find that both mergers and tidal perturbations can significantly swing spin vectors, with larger changes in spin vector for smaller pericentre distances. Strong tidal perturbations are also correlated with the changes in stellar mass and specific angular momentum of satellite galaxies. However, changes in spin vector can often result in a canceling out of previous changes. As a result, the integrated angular change is always much larger than the angular change measured at any instant. Also, overall the majority of satellite galaxies do not undergo mergers or sufficiently strong tidal perturbation after infall into clusters, and thus they end up suffering little change to their spin vectors. Taken as a whole, these results suggest that any signatures of spin alignment from the large scale structure will be preserved in the cluster environment for many gigayears. △ Less

Submitted 23 July, 2018; originally announced July 2018.

Comments: 17 pages, 10 figures, Accepted for publication in ApJ