Source author record

Suk-Jin Yoon

Suk-Jin Yoon appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

30works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

30 published item(s)

preprint2022arXiv

How the Galaxy Stellar Spins Acquire a Peculiar Tidal Connection?

We explore how the galaxy stellar spins acquire a peculiar tendency of being aligned with the major principal axes of the local tidal fields, in contrast to their DM counterparts which tend to be perpendicular to them, regardless of their masses. Analyzing the halo and subhalo catalogs from the IllustrisTNG 300 hydrodynamic simulations at $z\le 1$, we determine the cosines of the alignment angles, $\cosα$, between the galaxy stellar and DM spins. Creating four $\cosα$-selected samples of the galaxies and then controlling them to share the same density and mass distributions, we determine the average strengths of the alignments between the galaxy stellar spins and the tidal tensor major axes over each sample. It is clearly shown that at $z\le 0.5$ the more severely the galaxy stellar spin directions deviate from the DM counterparts, the stronger the peculiar tidal alignments become. Taking the ensemble averages of such galaxy properties as central blackhole to stellar mass ratio, specific star formation rate, formation epoch, stellar-to-total mass ratio, velocity dispersions, average metallicity, and degree of the cosmic web anisotropy over each sample, we also find that all of these properties exhibit either strong correlations or anti-correlations with $\cosα$. Our results imply that the peculiar tidal alignments of the galaxy stellar spins may be caused by anisotropic occurrence of some baryonic process responsible for discharging stellar materials from the galaxies along the tidal major directions at $z<1$.

preprint2022arXiv

Nonlinear Color-Metallicity Relations of Globular Clusters. XI. Nonlinearity Effect Revealed by NGC 5128 (Centaurus A) and NGC 4594 (Sombrero) Galaxies

Metallicity distributions (MDs) of globular clusters (GCs) provide crucial clues for the assembly and star formation history of their host galaxies. GC colors, when GCs are old, have been used as a proxy of GC metallicities. Bimodal GC color distributions (CDs) observed in most early-type galaxies have been interpreted as bimodal MDs for decades, suggesting the presence of merely two GC subpopulations within single galaxies. However, the conventional view has been challenged by a new theory that nonlinear metallicity-to-color conversion can cause bimodal CDs from unimodal MDs. The unimodal MDs seem natural given that MDs involved many thousand protogalaxies. The new theory has been tested and corroborated by various observational and theoretical studies. Here we examine the nonlinear nature of GC color-metallicity relations (CMRs) using photometric and spectroscopic GC data of NGC 5128 (Centaurus A) and NGC 4594 (Sombrero), in comparison with stellar population simulations. We find that, with a slight offset in color, the overall shapes of observed and modeled CMRs agree well for all available colors. Diverse color-depending morphologies of GC CDs of the two galaxies are well reproduced based on their observed spectroscopic MDs via our CMR models. The results corroborate the nonlinear CMR interpretation of the GC color bimodality, shedding further light on theories of galaxy formation.

preprint2022arXiv

Warped Disk Galaxies. I. Linking U type Warps in Groups/Clusters to Jellyfish Galaxies

arped disk galaxies are classified into two morphologies: S- and U-types. Conventional theories routinely attribute both types to galactic tidal interaction and/or gas accretion, but reproducing of U-types in simulations is extremely challenging. Here we investigate whether both types are governed by the same mechanisms using the most extensive sample of $\sim$8000 nearby (0.02\,$<$\,z\,$<$\,0.06) massive ($M_{*}/M_{\odot}$\,$>$\,$10^9$) edge-on disks from SDSS. We find that U-types show on average bluer optical colors and higher specific star formation rate (sSFR) than S-types, with more strongly warped U-types having higher sSFR. We also find that while the S-type warp properties correlate with the tidal force by the nearest neighbor regardless of the environment, there is no such correlation for U-types in groups/clusters, suggesting a non-tidal environmental could be at play for U-types, such as ram pressure stripping (RPS). Indeed, U-types are more common in groups/clusters than in fields and they have stellar mass, gas fraction, sSFR enhancement and phase-space distribution closely analogous to RPS-induced jellyfish galaxies in clusters. We furthermore show that the stellar disks of most RPS galaxies in the IllustirsTNG simulation are warped in U-shape and bent in opposite direction of stripped gas tails, satisfying theoretical expectations for stellar warps embeded in jellyfishes. We therefore suggest that despite the majority of U-types that live in fields being still less explained, RPS can be an alternative origin for those in groups/clusters.

preprint2021arXiv

Introducing a new multi-particle collision method for the evolution of dense stellar systems II. Core collapse

In a previous paper we introduced a new method for simulating collisional gravitational $N$-body systems with linear time scaling on $N$, based on the Multi-Particle Collision (MPC) approach. This allows us to simulate globular clusters with a realistic number of stellar particles in a matter of hours on a typical workstation. We evolve star clusters containing up to $10^6$ stars to core collapse and beyond. We quantify several aspects of core collapse over multiple realizations and different parameters, while always resolving the cluster core with a realistic number of particles. We run a large set of N-body simulations with our new code. The cluster mass function is a power-law with no stellar evolution, allowing us to clearly measure the effects of the mass spectrum. Leading up to core collapse, we find a power-law relation between the size of the core and the time left to core collapse. Our simulations thus confirm the theoretical self-similar contraction picture but with a dependence on the slope of the mass function. The time of core collapse has a non-monotonic dependence on the slope, which is well fit by a parabola. This holds also for the depth of core collapse and for the dynamical friction timescale of heavy particles. Cluster density profiles at core collapse show a broken power law structure, suggesting that central cusps are a genuine feature of collapsed cores. The core bounces back after collapse, and the inner density slope evolves to an asymptotic value. The presence of an intermediate-mass black hole inhibits core collapse. We confirm and expand on several predictions of star cluster evolution before, during, and after core collapse. Such predictions were based on theoretical calculations or small-size direct $N$-body simulations. Here we put them to the test on MPC simulations with a much larger number of particles, allowing us to resolve the collapsing core.

preprint2021arXiv

Living with Neighbors. III. The Origin of the Spin$-$Orbit Alignment of Galaxy Pairs: A Neighbor versus the Large-scale Structure

Recent observations revealed a coherence between the spin vector of a galaxy and the orbital motion of its neighbors. We refer to the phenomenon as "the spin$-$orbit alignment (SOA)" and explore its physical origin via the IllustrisTNG simulation. This is the first study to utilize a cosmological hydrodynamic simulation to investigate the SOA of galaxy pairs. In particular, we identify paired galaxies at $z = 0$ having the nearest neighbor with mass ratios from 1/10 to 10 and calculate the spin$-$orbit angle for each pair. Our results are as follows. (a) There exists a clear preference for prograde orientations (i.e., SOA) for galaxy pairs, qualitatively consistent with observations. (b) The SOA is significant for both baryonic and dark matter spins, being the strongest for gas and the weakest for dark matter. (c) The SOA is stronger for less massive targets and for targets having closer neighbors. (d) The SOA strengthens for galaxies in low-density regions, and the signal is dominated by central$-$satellite pairs in low-mass halos. (e) There is an explicit dependence of the SOA on the duration of interaction with its current neighbor. Taken together, we propose that the SOA witnessed at $z = 0$ has been developed mainly by interactions with a neighbor for an extended period of time, rather than tidal torque from the ambient large-scale structure.

preprint2020arXiv

Living with Neighbors. III. Scrutinizing the Spin$-$Orbit Alignment of Interacting Dark Matter Halo Pairs

We present that the spin$-$orbit alignment (SOA; i.e., the angular alignment between the spin vector of a halo and the orbital angular momentum vector of its neighbor) provides an important clue to how galactic angular momenta develop. In particular, we identify virial-radius-wise contact halo pairs with mass ratios from 1/3 to 3 in a set of cosmological $N$-body simulations, and divide them into merger and flyby subsamples according to their total (kinetic+potential) energy. In the spin$-$orbit angle distribution, we find a significant SOA in that $75.0\pm0.6$ % of merging neighbors and $58.7\pm0.6$ % of flybying neighbors are on the prograde orbit. The overall SOA of our sample is mainly driven by fast-rotating halos, corroborating that a well-aligned interaction spins a halo faster. More interestingly, we find for the first time a strong number excess of nearly perpendicular but still prograde interactions ($\sim75^{\circ}$) in the spin$-$orbit angle distribution for both the merger and flyby cases. Such prograde-polar interactions predominate for slow-rotating halos, testifying that misaligned interactions reduce the halos' spin. The frequency of the prograde-polar interactions correlates with the halo mass, yet anticorrelates with the large-scale density. This instantly invokes the spin-flip phenomenon that is conditional on the mass and environment. The prograde-polar interaction will soon flip the spin of a slow-rotator to align with its neighbor's orbital angular momentum. Finally, we propose a scenario that connects the SOA to the ambient large-scale structure based on the spin-flip argument.

preprint2020arXiv

MCG+07-20-052: Interacting dwarf pair in a group environment

We present an observational study of the interacting pair of dwarf galaxies, MCG+07-20-052, in the vicinity of Milky Way mass spiral galaxy NGC 2998. MCG+07-20-052 is located at a sky-projected distance of 105 kpc from NGC 2998 and the two have a relative line-of-sight velocity of 60 kms. We observed tidal tail-like extensions on both members (D1 and D2) of the interacting pair MCG+07-20-052. The interacting dwarf galaxies, D1 and D2, have B-band absolute magnitudes of $-$17.17 and $-$17.14 mag, respectively, and D2 is significantly bluer than D1. We obtained HI 21 cm line data of the NGC 2998 system using the Giant Metrewave Radio Telescope (GMRT) to get a more detailed view of the neutral hydrogen (HI) emission in the interacting dwarf galaxies and in the galaxy members of the NGC 2998 group. Evidence of a merger between the dwarf galaxies in the MCG+07-20-052 is also present in the HI kinematics and morphology where we find that the HI is mostly concentrated around D2, which also shows a higher level of star-forming activity and bluer $g-r$ color index compared to D1. In addition, we detect extended tenuous HI emission around another member galaxy, NGC 3006, located close to the MCG+07-20-052-pair at a sky-projected distance of 41 kpc. We compare here our results from the MCG+07-20-052 pair-NGC 2998 system with other known LMC-SMC-Milky Way type systems and discuss the possible origin of the dwarf-dwarf interaction.

preprint2020arXiv

MOCCA SURVEY Database I: Binary Black Hole Mergers from Globular Clusters with Intermediate Mass Black Holes

The dynamical formation of black hole binaries in globular clusters that merge due to gravitational waves occurs more frequently in higher stellar density. Meanwhile, the probability to form intermediate mass black holes (IMBHs) also increases with the density. To explore the impact of the formation and growth of IMBHs on the population of stellar mass black hole binaries from globular clusters, we analyze the existing large survey of Monte-Carlo globular cluster simulation data (MOCCA SURVEY Database I). We show that the number of binary black hole mergers agrees with the prediction based on clusters' initial properties when the IMBH mass is not massive enough or the IMBH seed forms at a later time. However, binary black hole formation and subsequent merger events are significantly reduced compared to the prediction when the present-day IMBH mass is more massive than $\sim10^4 \rm M_{\odot}$ or the present-day IMBH mass exceeds about 1 per cent of cluster's initial total mass. By examining the maximum black hole mass in the system at the moment of black hole binary escaping, we find that $\sim$ 90 per cent of the merging binary black holes escape before the formation and growth of the IMBH. Furthermore, large fraction of stellar mass black holes are merged into the IMBH or escape as single black holes from globular clusters in cases of massive IMBHs, which can lead to the significant under-population of binary black holes merging with gravitational waves by a factor of 2 depending on the clusters' initial distributions.

preprint2020arXiv

Narrowband Ca Photometry for Dwarf Spheroidal Galaxies. I. Chemostructural Study on Draco, Sextans, and Canes Venatici I

A few dozen dwarf satellite galaxies of the Milky Way have been discovered, which are often viewed as the remaining building blocks of our Galaxy. The follow-up spectroscopy showed that dwarf galaxies have a sizeable spread in their metallicities. Several scenarios were suggested to explain the metallicity spread, which can be tested by the structural patterns of stellar subpopulations with distinct metallicities. However, such chemical plus structural examination, to which we refer to as "chemostructural study" is hindered by the lack of stars with spectroscopic metallicity. Here we propose the Ca$-$$by$ photometry as an alternative way to secure metallicities for a 2$-$3 orders of magnitude larger stellar sample than the spectroscopic sample and thus enable us to perform a chemostructural study on dwarf galaxies. In particular, we use the $hk$ index [$\equiv($Ca$-b)-(b-y)$], whose validity as a photometric metallicity indicator (and crass insensitivity to age) for red-giant-branch stars was upheld via Galactic globular clusters, and observe three dwarf spheroidal galaxies$-$Draco, Sextans, and Canes Venatici I (CVnI)$-$with Subaru/Suprime-Cam. We find that in all the galaxies the metal-rich stellar populations are more centrally concentrated than the metal-poor counterparts, suggesting that the central regions of the galaxies underwent extended star formation. Such a negative radial metallicity gradient for Sextans and CVnI opposes to the traditional spectroscopic results. We also find that their metallicity distribution functions (MDFs) can be characterized by a unimodal, skewed Gaussian shape with a metal-rich peak and a metal-poor tail. We discuss their features in the chemostructure and MDF in terms of dwarf galaxy formation theories.

preprint2020arXiv

Nonlinear Color-Metallicity Relations of Globular Clusters. IX. Different Radial Number Density Profiles between Blue and Red Clusters

The optical colors of globular clusters (GCs) in most large early-type galaxies are bimodal. Blue and red GCs show a sharp difference in the radial profile of their surface number density in the sense that red GCs are more centrally concentrated than blue GCs. An instant interpretation is that there exist two distinct GC subsystems having different radial distributions. This view, however, was challenged by a scenario in which, due to the nonlinear nature of the GC metallicity-to-color transformation for old ($\gtrsim$10 Gyr) GCs, a broad unimodal metallicity spread can exhibit a bimodal color distribution. Here we show, by simulating the radial trends in the GC color distributions of the four nearby giant elliptical galaxies (M87, M49, M60, and NGC 1399), that the difference in the radial profile between blue and red GCs stems naturally from the metallicity-to-color nonlinearity plus the well-known radial metallicity gradient of GC systems. The model suggests no or little radial variation in GC age even out to $\sim$20${R}_{\rm eff}$. Our results provide a simpler solution to the distinct radial profiles of blue and red GCs that does not necessarily invoke the presence of two GC subsystems and further fortify the nonlinearity scenario for the GC color bimodality phenomenon.

preprint2020arXiv

Starbursting Nuclei in Old Dwarf Galaxies

Nuclei of early-type dwarf galaxies (dEs) are usually younger than the galaxy main body, and such discrepancy in age has been a puzzle. To explore the origin of young nuclei in dEs, we study a sample of dEs having compact star-forming blobs that are visually similar to dEs' nuclei but by far bluer. We find that (1) the compact star-forming blobs have a typical stellar mass of one percent of the host galaxy stellar mass; (2) some of the blobs are positioned slightly off from the center of the galaxies; (3) the H$α$ equivalent width measured from the publicly available Sloan Digital Sky Survey fiber spectroscopy shows their formation ages being an order of few Mega-year; and (4) their emission line metallicities, 12\,+\,log(O/H), are as high as the solar value, while the underlying galaxies have the typical stellar populations of dEs, i.e., log(Z/Z$_{\sun}$)\,$\sim$\,$-0.8$. Based on the results, we argue that the central star-forming blobs can provide a caught-in-the-act view of nuclei formation in dEs, and discussing possible formation mechanisms of young nuclei in old dEs. We particularly propose that these off-centered compact star-forming regions may act as seeds of nuclei as proposed in the `wet migration' scenario of \cite{Guillard16}.

preprint2016arXiv

A new catalog of homogenised absorption line indices for Milky Way globular clusters from high-resolution integrated spectroscopy

We perform integrated spectroscopy of 24 Galactic globular clusters. Spectra are observed from one core radius for each cluster with a high wavelength resolution of ~2.0 A FWHM. In combination with two existing data sets from Puzia et al. (2002) and Schiavon et al. (2005), we construct a large database of Lick spectral indices for a total of 53 Galactic globular clusters with a wide range of metallicities, -2.4 < [Fe/H] < 0.1, and various horizontal-branch morphologies. The empirical index-to-metallicity conversion relationships are provided for the 20 Lick indices for the use of deriving metallicities for remote, unresolved stellar systems.

preprint2016arXiv

Nonlinear Color--Metallicity Relations of Globular Clusters. VI. On Calcium II Triplet Based Metallicities of Globular Clusters in Early-type Galaxies

The metallicity distribution function of globular clusters (GCs) in galaxies is a key to understanding galactic formation and evolution. The calcium II triplet (CaT) index has recently become a popular metal abundance indicator thanks to its sensitivity to GC metallicity. Here we revisit and assess the reliability of CaT as a metallicity indicator using our new stellar population synthesis simulations based on empirical, high-resolution fluxes. The model shows that the CaT strength of old ($>$ 10 Gyr) GCs is proportional to ${\rm [Fe/H]}$ below $-0.5$. In the modest metal-rich regime, however, CaT does not increase anymore with ${\rm [Fe/H]}$ due to the little contribution from coolest red giant stars to the CaT absorption. The nonlinear nature of the color--$CaT$ relation is confirmed by the observations of GCs in nearby early-type galaxies. This indicates that the CaT should be used carefully when deriving metallicities of metal-rich stellar populations. Our results offer an explanation for the observed sharp difference between the color and $CaT$ distributions of GCs in the same galaxies. We take this as an analogy to the view that metallicity--color and metallicity--Lick index nonlinearity of GCs is primarily responsible for their observed "bimodal" distributions of colors and absorption indices.

preprint2016arXiv

Odyssey: A Public GPU-Based Code for General-Relativistic Radiative Transfer in Kerr Spacetime

General-relativistic radiative transfer (GRRT) calculations coupled with the calculation of geodesics in the Kerr spacetime are an essential tool for determining the images, spectra and light curves from matter in the vicinity of black holes. Such studies are especially important for ongoing and upcoming millimeter/submillimeter (mm/sub-mm) Very Long Baseline Interferometry (VLBI) observations of the supermassive black holes at the centres of Sgr A^{*} and M87. To this end we introduce Odyssey, a Graphics Processing Unit(GPU)-based code for ray tracing and radiative transfer in the Kerr spacetime. On a single GPU, the performance of Odyssey can exceed 1 nanosecond per photon, per Runge-Kutta integration step. Odyssey is publicly available, fast, accurate, and flexible enough to be modified to suit the specific needs of new users. Along with a Graphical User Interface (GUI) powered by a video-accelerated display architecture, we also present an educational software tool, Odyssey_Edu, for showing in real time how null geodesics around a Kerr black hole vary as a function of black hole spin and angle of incidence onto the black hole.

preprint2015arXiv

Power spectra of outflow-driven turbulence

We investigate the power spectra of outflow-driven turbulence through high-resolution three-dimensional isothermal numerical simulations where the turbulence is driven locally in real-space by a simple spherical outflow model. The resulting turbulent flow saturates at an average Mach number of ~2.5 and is analysed through density and velocity power spectra, including an investigation of the evolution of the solenoidal and compressional components. We obtain a shallow density power spectrum with a slope of ~-1.2 attributed to the presence of a network of localised dense filamentary structures formed by strong shock interactions. The total velocity power spectrum slope is found to be ~-2.0, representative of Burgers shock dominated turbulence model. The density weighted velocity power spectrum slope is measured as ~-1.6, slightly less than the expected Kolmogorov scaling value (slope of -5/3) found in previous works. The discrepancy may be caused by the nature of our real space driving model and we suggest there is no universal scaling law for supersonic compressible turbulence. We find that on average, solenoidal modes slightly dominate in our turbulence model as the interaction between strong curved compressible shocks generates solenoidal modes, and compressible modes decay faster.

preprint2014arXiv

Formation of Warped Disks by Galactic Fly-by Encounters. I. Stellar Disks

Warped disks are almost ubiquitous among spiral galaxies. Here we revisit and test the `fly-by scenario' of warp formation, in which impulsive encounters between galaxies are responsible for warped disks. Based on N-body simulations, we investigate the morphological and kinematical evolution of the stellar component of disks when galaxies undergo fly-by interactions with adjacent dark matter halos. We find that the so-called `S'-shaped warps can be excited by fly-bys and sustained for even up to a few billion years, and that this scenario provides a cohesive explanation for several key observations. We show that disk warp properties are governed primarily by the following three parameters; (1) the impact parameter, i.e., the minimum distance between two halos, (2) the mass ratio between two halos, and (3) the incident angle of the fly-by perturber. The warp angle is tied up with all three parameters, yet the warp lifetime is particularly sensitive to the incident angle of the perturber. Interestingly, the modeled S-shaped warps are often non-symmetric depending on the incident angle. We speculate that the puzzling U- and L-shaped warps are geometrically superimposed S-types produced by successive fly-bys with different incident angles, including multiple interactions with a satellite on a highly elongated orbit.

preprint2013arXiv

Density distributions of outflow driven turbulence

Protostellar jets and outflows are signatures of star formation and promising mechanisms for driving supersonic turbulence in molecular clouds. We quantify outflow-driven turbulence through three-dimensional numerical simulations using an isothermal version of the robust total variation diminishing code. We drive turbulence in real-space using a simplified spherical outflow model, analyse the data through density probability distribution functions (PDF), and investigate the Core Formation Rate per free-fall time (CFR_ff). The real-space turbulence driving method produces a negatively skewed density PDF possessing an enhanced tail on the low-density side. It deviates from the log-normal distributions typically obtained from Fourier-space turbulence driving at low densities, but can provide a good fit at high-densities, particularly in terms of mass weighted rather than volume weighted density PDF. Due to this fact, we suggest that the CFR_ff determined from a Fourier-driven turbulence model could be comparable to that of our particular real-space driving model, which has a ratio of solenoidal to compressional components from the resulting turbulence velocity fields of ~0.6.

preprint2013arXiv

Nonlinear Color-Metallicity Relations of Globular Clusters. IV. Testing the Nonlinearity Scenario for Color Bimodality via HST/WFC3 u-band Photometry of M84 (NGC 4374)

Color distributions of globular clusters (GCs) in most massive galaxies are bimodal. Assuming linear color-to-metallicity conversions, bimodality is viewed as the presence of merely two GC subsystems with distinct metallicities, which serves as a critical backbone of various galaxy formation theories. Recent studies, however, revealed that the color-metallicity relations (CMRs) often used to derive GC metallicities (e.g., CMRs of g-z, V-I and C-T1) are in fact inflected. Such inflection can create bimodal color distributions if the underlying GC metallicity spread is simply broad as expected from the hierarchical merging paradigm of galaxy formation. In order to test the nonlinear-CMR scenario for GC color bimodality, the u-band photometry is proposed because the u-related CMRs (e.g., CMRs of u-g and u-z) are theoretically predicted to be least inflected and most distinctive among commonly used optical CMRs. Here, we present Hubble Space Telescope (HST)/WFC3 F336W (u-band) photometry of the GC system in M84, a giant elliptical in the Virgo galaxy cluster. Combining the u data with the existing HST ACS/WFC g and z data, we find that the u-z and u-g color distributions are different from the g-z distribution in a very systematic manner and remarkably consistent with our model predictions based on the nonlinear-CMR hypothesis. The results lend further confidence to validity of the nonlinear-CMR scenario as an explanation for GC color bimodality. There are some GC systems showing bimodal spectroscopic metallicity, and in such systems the inflected CMRs often create stronger bimodality in the color domain.

preprint2013arXiv

Nonlinear Color-Metallicity Relations of Globular Clusters. V. Nonlinear Absorption-line Index versus Metallicity Relations and Bimodal Index Distributions of M31 Globular Clusters

Recent spectroscopy on the globular cluster (GC) system of M31 with unprecedented precision witnessed a clear bimodality in absorption-line index distributions of old GCs. Such division of extragalactic GCs, so far asserted mainly by photometric color bimodality, has been viewed as the presence of merely two distinct metallicity subgroups within individual galaxies and forms a critical backbone of various galaxy formation theories. Given that spectroscopy is a more detailed probe into stellar population than photometry, the discovery of index bimodality may point to the very existence of dual GC populations. However, here we show that the observed spectroscopic dichotomy of M31 GCs emerges due to the nonlinear nature of metallicity-to-index conversion and thus one does not necessarily have to invoke two separate GC subsystems. We take this as a close analogy to the recent view that metallicity-color nonlinearity is primarily responsible for observed GC color bimodality. We also demonstrate that the metallicity-sensitive magnesium line displays non-negligible metallicity-index nonlinearity and Balmer lines show rather strong nonlinearity. This gives rise to bimodal index distributions, which are routinely interpreted as bimodal metallicity distributions, not considering metallicity-index nonlinearity. Our findings give a new insight into the constitution of M31's GC system, which could change much of the current thought on the formation of GC systems and their host galaxies.

preprint2013arXiv

The effect of second generation populations on the integrated colors of metal-rich globular clusters in early-type galaxies

The mean color of globular clusters (GCs) in early-type galaxies is in general bluer than the integrated color of halo field stars in host galaxies. Metal-rich GCs often appear more associated with field stars than metal-poor GCs, yet show bluer colors than their host galaxy light. Motivated by the discovery of multiple stellar populations in Milky Way GCs, we present a new scenario in which the presence of second-generation (SG) stars in GCs is responsible for the color discrepancy between metal-rich GCs and field stars. The model assumes that the SG populations have enhanced helium abundance as evidenced by observations, and it gives a good explanation of the bluer optical colors of metal-rich GCs than field stars as well as strong Balmer lines and blue UV colors of metal-rich GCs. Ours may be complementary to the recent scenario suggesting the difference in stellar mass functions (MFs) as an origin for the GC-to-star color offset. A quantitative comparison is given between the SG and MF models.

preprint2013arXiv

Two distinct red giant branch populations in the globular cluster NGC 2419 as tracers of a merger event in the Milky Way

Recent spectroscopic observations of the outer halo globular cluster (GC) NGC 2419 show that it is unique among GCs, in terms of chemical abundance patterns, and some suggest that it was originated in the nucleus of a dwarf galaxy. Here we show, from the Subaru narrow-band photometry employing a calcium filter, that the red giant-branch (RGB) of this GC is split into two distinct subpopulations. Comparison with spectroscopy has confirmed that the redder RGB stars in the $hk$[=(Ca$-b)-(b-y)$] index are enhanced in [Ca/H] by $\sim$0.2 dex compared to the bluer RGB stars. Our population model further indicates that the calcium-rich second generation stars are also enhanced in helium abundance by a large amount ($Δ$Y = 0.19). Our photometry, together with the results for other massive GCs (e.g., $ω$ Cen, M22, and NGC 1851), suggests that the discrete distribution of RGB stars in the $hk$ index might be a universal characteristic of this growing group of peculiar GCs. The planned narrow-band calcium photometry for the Local Group dwarf galaxies would help to establish an empirical connection between these GCs and the primordial building blocks in the hierarchical merging paradigm of galaxy formation.

preprint2013arXiv

Wide-Field Multiband Photometry of Globular Cluster Systems in the Fornax Galaxy Cluster

We present wide-field multiband photometry of globular cluster (GC) systems in NGC 1399, NGC 1404, and NGC 1387 located at the central region of the Fornax galaxy cluster. Observation was carried out through U, B, V, and I bands, which marks one of the widest and deepest U-band studies on extragalactic GC systems. The present U-band photometry enables us to significantly reduce the contamination by a factor of two for faint sources (Vo~23.5). The main results based on some 2000 GC candidates around NGC 1399, NGC 1404, and NGC 1387 are as follows: (1) the GC system in each galaxy exhibits bimodal color distributions in all colors examined, but the shape of color histograms varies systematically depending on colors; (2) NGC 1399 shows that the mean colors of both blue and red GCs become bluer with increasing galactocentric radius; (3) NGC 1399 shows overabundance of GCs in the directions of NGC 1404 and NGC 1387, indicating their ongoing interactions; and (4) NGC 1399 also exhibits a ~0.5' offset between the center of the inner GC distribution and the galaxy's optical center, suggesting that NGC 1399 is not yet dynamically relaxed and may be undergoing merger events.

preprint2012arXiv

A Keck/LRIS Spatially-Resolved Spectroscopic Study of a LINER Galaxy SDSS J091628.05+420818.7

Using spatially-resolved spectra obtained with the Low Resolution Imaging Spectrometer at the Keck I telescope, we investigate the nature of ionizing sources and kinematic properties of emission-line gas in a LINER galaxy SDSS J091628.05+420818.7, which is a nearby (z = 0.0241) and bright (M_r = -20.2) early-type galaxy. After subtracting stellar absorption features using a combination of simple stellar population models, we measure the flux, line-of-sight velocity, and velocity dispersion of four emission lines, i.e., Hα, Hβ, [O III] λ5007, and [N II] λ6584, to study radial change of emission-line fluxes and velocities. Compared to the point-spread-function of the observation, the emission-line region is slightly extended but comparable to the seeing size. The central concentration of emission-line gas suggests that ionization is triggered by a nuclear source, excluding old stellar population as ionizing sources. We find that emission-line gas is counter-rotating with respect to stellar component and that the [O III] λ5007 line is blueshifted compared to other emission lines, possibly due to an outflow.

preprint2012arXiv

Initial Size Distribution of the Galactic Globular Cluster System

Despite the importance of their size evolution in understanding the dynamical evolution of globular clusters (GCs) of the Milky Way, studies are rare that focus specifically on this issue. Based on the advanced, realistic Fokker-Planck (FP) approach, we predict theoretically the initial size distribution (SD) of the Galactic GCs along with their initial mass function and radial distribution. Over one thousand FP calculations in a wide parameter space have pinpointed the best-fit initial conditions for the SD, mass function, and radial distribution. Our best-fit model shows that the initial SD of the Galactic GCs is of larger dispersion than today's SD, and that typical projected half-light radius of the initial GCs is ~4.6 pc, which is 1.8 times larger than that of the present-day GCs (~2.5 pc). Their large size signifies greater susceptibility to the Galactic tides: the total mass of destroyed GCs reaches 3-5x10^8 M_sun$, several times larger than the previous estimates. Our result challenges a recent view that the Milky Way GCs were born compact on the sub-pc scale, and rather implies that (1) the initial GCs are generally larger than the typical size of the present-day GCs, (2) the initially large GCs mostly shrink and/or disrupt as a result of the galactic tides, and (3) the initially small GCs expand by two-body relaxation, and later shrink by the galactic tides.

preprint2012arXiv

UV Properties of Galactic Globular Clusters with GALEX I. The Color-Magnitude Diagrams

We present GALEX data for 44 Galactic globular clusters obtained during 3 GALEX observing cycles between 2004 and 2008. This is the largest homogeneous data set on the UV photometric properties of Galactic globular clusters ever collected. The sample selection and photometric analysis are discussed, and color-magnitude diagrams are presented. The blue and intermediate-blue horizontal branch is the dominant feature of the UV color-magnitude diagrams of old Galactic globular clusters. Our sample is large enough to display the remarkable variety of horizontal branch shapes found in old stellar populations. Other stellar types that are obviously detected are blue stragglers and post core-He burning stars. The main features of UV color-magnitude diagrams of Galactic globular clusters are briefly discussed. We establish the locus of post-core He burning stars in the UV color-magnitude diagram and present a catalog of candidate AGB-manqu é, post early-AGB, and post-AGB stars within our cluster sample.

preprint2012arXiv

Yonsei Evolutionary Population Synthesis (YEPS) Model. I. Spectroscopic Evolution of Simple Stellar Populations

We present a series of papers on the year-2012 version of Yonsei Evolutionary Population Synthesis (YEPS) model which is constructed on over 20 years of heritage. This first paper delineates the spectroscopic aspect of integrated light from stellar populations older than 1 Gyr. The standard YEPS is based on the most up-to-date Yonsei-Yale stellar evolutionary tracks and BaSel 3.1 flux libraries, and provides absorption line indices of the Lick/IDS system and high-order Balmer lines for simple stellar populations as functions of stellar parameters, such as metallicity, age and α-element mixture. Special care has been taken to incorporate systematic contribution from horizontal-branch stars which alters the temperature-sensitive Balmer lines significantly, resulting in up to 5 Gyr difference in age estimation of old, metal-poor stellar populations. We also find that the horizontal branches exert an appreciable effect not only on the Balmer lines but also on the metallicity-sensitive lines including the magnesium index. This is critical to explain the intriguing bimodality found in index distributions of globular clusters in massive galaxies and to derive spectroscopic metallicities accurately from various indices. A full set of the spectroscopic and photometric YEPS model data of the entire parameter space is currently downloadable at http://web.yonsei.ac.kr/cosmic/data/YEPS.htm

preprint2011arXiv

Nonlinear Color-Metallicity Relations of Globular Clusters. II. A Test on the Nonlinearity Scenario for Color Bimodality Using the u-band Colors: The Case of M87 (NGC 4486)

The optical color distributions of globular clusters (GCs) in most large elliptical galaxies are bimodal. Based on the assumed linear relationship between GC colors and their metallicities, the bimodality has been taken as evidence of two GC subsystems with different metallicities in each galaxy and led to a number of theories in the context of galaxy formation. More recent observations and modeling of GCs, however, suggest that the color-metallicity relations (CMRs) are inflected, and thus colors likely trace metallicities in a nonlinear manner. The nonlinearity could produce bimodal color distributions from a broad underlying metallicity spread, even if it is unimodal. Despite the far-reaching implications, whether CMRs are nonlinear and whether the nonlinearity indeed causes the color bimodality are still open questions. Given that the spectroscopic refinement of CMRs is still very challenging, we here propose a new photometric technique to probe the possible nonlinear nature of CMRs. In essence, a color distribution of GCs is a "projected" distribution of their metallicities. Since the form of CMRs hinges on which color is used, the shape of color distributions varies depending significantly on the colors. Among other optical colors, the u-band related colors (e.g., u - g and u - z) are theoretically predicted to exhibit significantly less inflected CMRs than other preferred CMRs (e.g., for g - z). As a case study, we performed the HST/WFPC2 archival u-band photometry for the M87 (NGC 4486) GC system with confirmed color bimodality. We show that the u-band color distributions are significantly different from that of g - z, and consistent with our model predictions. With more u-band measurements, this method will support or rule out the nonlinear-CMR scenario for the origin of GC color bimodality with high confidence. [Abridged]

preprint2011arXiv

Nonlinear Color-Metallicity Relations of Globular Clusters. III. On the Discrepancy in Metallicity between Globular Cluster Systems and their Parent Elliptical Galaxies

One of the conundrums in extragalactic astronomy is the discrepancy in observed metallicity distribution functions (MDFs) between the two prime stellar components of early-type galaxies-globular clusters (GCs) and halo field stars. This is generally taken as evidence of highly decoupled evolutionary histories between GC systems and their parent galaxies. Here we show, however, that new developments in linking the observed GC colors to their intrinsic metallicities suggest nonlinear color-to-metallicity conversions, which translate observed color distributions into strongly-peaked, unimodal MDFs with broad metal-poor tails. Remarkably, the inferred GC MDFs are similar to the MDFs of resolved field stars in nearby elliptical galaxies and those produced by chemical evolution models of galaxies. The GC MDF shape, characterized by a sharp peak with a metal-poor tail, indicates a virtually continuous chemical enrichment with a relatively short timescale. The characteristic shape emerges across three orders of magnitude in the host galaxy mass, suggesting a universal process of chemical enrichment among various GC systems. Given that GCs are bluer than field stars within the same galaxy, it is plausible that the chemical enrichment processes of GCs ceased somewhat earlier than that of field stellar population, and if so, GCs preferentially trace the major, vigorous mode of star formation events in galactic formation. We further suggest a possible systematic age difference among GC systems, in that the GC systems in more luminous galaxies are older. This is consistent with the downsizing paradigm of galaxies and supports additionally the similar nature shared by GCs and field stars. Our findings suggest that GC systems and their parent galaxies have shared a more common origin than previously thought, and hence greatly simplify theories of galaxy formation.

preprint2011arXiv

The effect of helium-enhanced stellar populations on the ultraviolet-upturn phenomenon of early-type galaxies

Recent observations and modeling of globular clusters with multiple populations strongly indicate the presence of super helium-rich subpopulations in old stellar systems. Motivated by this, we have constructed new population synthesis models with and without helium-enhanced subpopulations to investigate their impact on the UV-upturn phenomenon of quiescent early-type galaxies. We find that our models with helium- enhanced subpopulations can naturally reproduce the strong UV-upturns observed in giant elliptical galaxies assuming an age similar to that of old globular clusters in the Milky Way. The major source of far-UV (FUV) flux, in this model, is relatively metal- poor and helium-enhanced hot horizontal branch stars and their progeny. The Burstein et al. (1988) relation of the F U V - V color with metallicity is also explained either by the variation of the fraction of helium-enhanced subpopulations or by the spread in mean age of stellar populations in early-type galaxies.

preprint2009arXiv

The Presence of Two Distinct Red Giant Branches in the Globular Cluster NGC 1851

There is a growing body of evidence for the presence of multiple stellar populations in some globular clusters, including NGC 1851. For most of these peculiar globular clusters, however, the evidence for the multiple red giant-branches (RGBs) having different heavy elemental abundances as observed in Omega Centauri is hitherto lacking, although spreads in some lighter elements are reported. It is therefore not clear whether they also share the suggested dwarf galaxy origin of Omega Cen or not. Here we show from the CTIO 4m UVI photometry of the globular cluster NGC 1851 that its RGB is clearly split into two in the U - I color. The two distinct RGB populations are also clearly separated in the abundance of heavy elements as traced by Calcium, suggesting that the type II supernovae enrichment is also responsible, in addition to the pollutions of lighter elements by intermediate mass asymptotic giant branch stars or fast-rotating massive stars. The RGB split, however, is not shown in the V - I color, as indicated by previous observations. Our stellar population models show that this and the presence of bimodal horizontal-branch distribution in NGC 1851 can be naturally reproduced if the metal-rich second generation stars are also enhanced in helium.