Search

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\,$$\,$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., 28 pages, 15 figures, accpeted for publication in ApJ

We report the discovery of a rare early-type dwarf galaxy (dE), SDSS J125651.47+163024.2 (hereafter dE1256), possessing a tidal feature that was likely built up by accretion of an even smaller dwarf galaxy. dE1256 is located in a nearly isolated environment, at the outskirt of the Virgo cluster. A detailed morphological examination reveals that the accreted stellar population is mainly deposited in the outer part of dE1256, where the tidal tail is most prominent. The inner part of dE1256 is perfectly modeled with a simple S\'ersic function of index n = 0.63 and half-light radius R$_{h}$ = 0.6 kpc, but in contrast, the entire galaxy has a size of R$_{h}$ = 1.2 kpc. The mass ratio between the host and the putative accreted dwarf galaxy is calculated to be 5:1, assuming that the observed two components, inner S\'ersic, and outer tidal tail residual, represent the host's and accreted galaxy's stellar populations, respectively. We suggest that while the accretion contributes only 20% of the overall stellar population, the size of dE1256 grew by a factor of two via the accretion event. Our results provide, for the first time, strong observational evidence that a dE is undergoing a two-phase growth, a common phenomenon for massive galaxies., Comment: Accepted for publication in MNRAS

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\alpha$, between the galaxy stellar and DM spins. Creating four $\cos\alpha$-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\alpha$. 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, Comment: accepted for publication in ApJ, 11 figures, 1 tables, minor revision

We obtained spectra of some 140 globular clusters (GCs) associated with the Virgo central cD galaxy M87 with the Subaru/FOCAS MOS mode. The fundamental properties of GCs such as age, metallicity and $\alpha$-element abundance are investigated by using simple stellar population models. It is confirmed that the majority of M87 GCs are as old as, more metal-rich than, and more enhanced in $\alpha$-elements than the Milky Way GCs. Our high-quality, homogeneous dataset enables us to test the theoretical prediction of inflected color$-$metallicity relations (CMRs). The nonlinear-CMR hypothesis entails an alternative explanation for the widely observed GC color bimodality, in which even a unimodal metallicity spread yields a bimodal color distribution by virtue of nonlinear metallicity-to-color conversion. The newly derived CMRs of old, high-signal-to-noise-ratio GCs in M87 (the $V-I$ CMR of 83 GCs and the $M-T2$ CMR of 78 GCs) corroborate the presence of the significant inflection. Furthermore, from a combined catalog with the previous study on M87 GC spectroscopy, we find that a total of 185 old GCs exhibit a broad, unimodal metallicity distribution. The results corroborate the nonlinear-CMR interpretation of the GC color bimodality, shedding further light on theories of galaxy formation., Comment: 28 pages, 15 figures, 8 tables

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., Comment: 25 pages, 12 figures

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., Comment: 16 pages, 12 figures, accepted for publication in ApJ

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., Comment: 8 pages, 9 figs., 1 tab. Version accepted for publication in A&A

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., 18 pages, 11 figures, 3 tables, accepted for publication in ApJ

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., 9 pages, 8 figures, Accepted for publication in MNRAS

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., Accepted for publication in ApJS, data will be available at http://vo.kasi.re.kr/Phot_dSphs/

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$\alpha$ 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}., Comment: Accepted for publication in APJL

A numerical detection of the mass-dependent spin transition of the galaxies is presented. Analyzing a sample of the galaxies with stellar masses in the range of $10^{9}< (M_{\star}/M_{\odot})\le 10^{11}$ from the IllustrisTNG300-1 simulations, we explore the alignment tendency between the galaxy baryon spins and the three eigenvectors of the linearly reconstructed tidal field as a function of $M_{\star}$ and its evolution in the redshift range of $0\le z \le 1.5$. Detecting a significant signal of the occurrence of the mass-dependent transition of the galaxy spins, we show that the centrals differ from the satellites in their spin transition type. As $M_{\star}$ increases beyond a certain threshold mass, the preferred directions of the central galaxy spins transit from the minor to the intermediate tidal eigenvectors (type two) at $z=0.5$ and $1$, while those of the satellites transit from the minor to the major tidal eigenvectors (type one) at $z=1$ and $1.5$. It is also shown that the mass range and type of the spin transition depend on the galaxy morphology, degree of the alignments between the baryon and total spin vectors as well as on the environmental density. Meanwhile, the stellar spins of the galaxies are found to yield a weak signal of the T1 transitions at $z=0$, whose strength and trend depend on the degree of the alignments between the stellar and baryon spins. The possible mechanisms responsible for the T1 and T2 spin transitions are discussed., Accepted for publication in ApJ, 16 figures, 5 tables, improved analysis, focus switched to the baryon spins

Spin$-$orbit alignment (SOA; i.e., the vector alignment between the halo spin and the orbital angular momentum of neighboring halos) provides an important clue to how galactic angular momenta develop. For this study, we extract virial-radius-wise contact halo pairs with mass ratios between 1/10 and 10 from a set of cosmological $N$-body simulations. In the spin--orbit angle distribution, we find a significant SOA in that 52.7%$\pm$0.2% of neighbors are on the prograde orbit. The SOA of our sample is mainly driven by low-mass target halos ($, 18 pages, 12 figures, accepted for publication in ApJ. arXiv admin note: text overlap with arXiv:2005.06479

Galaxies in pairs show enhanced star formation (SF) compared to their counterparts in isolation, which is often explained by the tidal effect of neighboring galaxies. Recent observations, however, reported that galaxies paired with early-type neighbors do not undergo the SF enhancement. Here we revisit the influence of neighbors using a large sample of paired galaxies from the Sloan Digital Sky Survey and a carefully constructed control sample of isolated counterparts. We find that star-forming neighbors enhance SF, and even more so for more star-forming (and closer) neighbors, which can be attributed to collisions of interstellar medium (ISM) leading to SF. We further find that, contrary to the anticipated tidal effect, quiescent neighbors quench SF, and even more so for more quiescent (and closer) neighbors. This seems to be due to removal of gas reservoirs via ram pressure stripping and gas accretion cutoff by hot gas halos of quiescent neighbors, on top of their paucity of ISM to collide to form stars. Our findings, especially the intimate connection of SF to the status and strength of neighbors' SF, imply that the hydrodynamic mechanisms, along with the tidal effect, play a crucial role during the early phase of galactic interactions., 15 pages, 9 figures, accepted for publication in ApJ

We report a possible discovery of a relic star cluster in the Sextans dwarf spheroidal galaxy. Using the \textit{hk} index ($\equiv$($Ca-b$)$-$($b-y$)) as a photometric metallicity indicator, we discriminate the metal-poor and metal-rich stars in the galaxy and find unexpected number density excess of metal-poor stars located 7'.7 ($\sim$190 pc in projected distance) away from the known galactic center. The $V-I$ color$-$magnitude diagram (CMD) for stars around the density excess reveals that both the main sequence and the giant branch are considerably narrower and redder than the bulk of field stars in Sextans. Our stellar population models show (a) that the narrow CMD is best reproduced by a simple stellar population with an age of $\sim$13 Gyr and [Fe/H] of $\sim$$-$2.3 dex, and (b) that the redder $V-I$ color of the $hk$-weak population is explained $only$ if it is $\sim$2 Gyr older than the field stars. The results lead us to conclude that the off-centered density peak is likely associated with an old, metal-poor globular cluster. The larger spatial extent ($>$ 80 pc in radius) and the smaller number of stars ($\sim$1000) than typical globular clusters point to a star cluster that is in the process of dissolution. The finding serves as the first detection of a surviving star cluster in Sextans, supporting previous suggestions of the presence of star cluster remnants in the galaxy. If confirmed, the survival of a relic star cluster until now implies a $cored$ dark matter halo profile for this dwarf galaxy., 10 pages, 5 figures, published in ApJL

We present an evolutionary population synthesis model of the surface brightness fluctuation (SBF) for normal and He-enriched simple stellar populations (SSPs). While our SBF model for the normal-He population agrees with other existing models, the He-rich population, containing hotter horizontal-branch stars and brighter red-clump stars than the normal-He population, entails a substantial change in the SBF of SSPs. We show that the SBF magnitudes are affected by He-rich populations at least $\sim$0.3~mag even in $I$- and near-IR bands at given colors, from which the SBF-based distances are often derived. Due to uncertainties both in observations and models, however, the SBFs of Galactic globular clusters and early-type galaxies do not allow verifying the He-enriched model. We propose that when combined with independent metallicity and age indicators such as ${\rm Mg}_2$ and ${\rm H}\beta$, the UV and optical SBFs can readily detect underlying He-rich populations in unresolved stellar systems at a distance out to $\gtrsim 20$\,Mpc. A full set of the spectro-photometric and SBF data for SSPs from the Yonsei Evolutionary Population Synthesis (YEPS) model is available for download at http://cosmic.yonsei.ac.kr/YEPS.htm., Comment: 24 pages, 14 figures, 4 tables, Accepted for publication in ApJS

The color distributions of globular clusters (GCs) in individual early-type galaxies show great diversity in their morphology. Based on the conventional "linear" relationship between the colors and metallicities of GCs, the GC metallicity distributions inferred from colors and in turn their formation histories, should be as diverse as they appear. In contrast, here we show that an alternative scenario rooted in the "nonlinear" nature of the metallicity-to-color transformation points to a simpler and more coherent picture. Our simulations of the color distributions for $\sim$80 GC systems in early-type galaxies from the ACS Virgo and Fornax Cluster Surveys suggest that the majority ($\sim$70$\%$) of early-type galaxies have old ($\sim$13 Gyr) and coeval GCs. Their variety in the color distribution morphology stems mainly from one parameter, the mean metallicity of a GC system. Furthermore, the color distributions of the remaining ($\sim$30$\%$) GC systems are also explained by the nonlinearity scenario, assuming additional young or intermediate-age GCs with a number fraction of $\sim$20$\%$ of underlying old GCs. Our results reinforce the nonlinearity explanation for the GC color bimodality and provide a new perspective on early-type galaxy formation in the cluster environment, such as the Virgo and Fornax galaxy clusters., 33 pages, 7 figures, 3 tables, accepted for publication in ApJS

We present a statistical analysis of flybys of dark matter halos compared to mergers using cosmological $N$-body simulations. We mainly focus on gravitationally interacting target halos with mass of $10^{10.8}-10^{13.0}h^{-1}M_{\odot}$, and their neighbors are counted only when the mass ratio is 1:3$-$3:1 and the distance is less than the sum of the virial radii of target and neighbor. The neighbors are divided into the flyby or merger samples if the pair's total energy is greater or smaller, respectively, than the capture criterion with consideration of dynamical friction. The main results are as follows: (a) The flyby fraction increases by up to a factor of 50 with decreasing halo mass and by up to a factor of 400 with increasing large-scale density, while the merger fraction does not show any significant dependencies on these two parameters; (b) The redshift evolution of the flyby fraction is twofold, increasing with redshift at $01$, while the merger fraction increases monotonically with redshift at $z=0\sim4$; (c) The multiple interactions with two or more neighbors are on average flyby-dominated, and their fraction has a mass and environment dependence similar to that for the flyby fraction; (d) Given that flybys substantially outnumber mergers toward $z=0$ (by a factor of five) and the multiple interactions are flyby-dominated, the flyby's contribution to galactic evolution is stronger than ever at the present epoch, especially for less massive halos and in the higher density environment. We propose a scenario that connects the evolution of the flyby and merger fractions to the hierarchical structure formation process., Comment: 21 pages, 12 figures, and 1 table, accepted for publication in ApJ

We present the largest publicly available catalog of interacting dwarf galaxies. It includes 177 nearby merging dwarf galaxies of stellar mass M$_{*}$ $, Comment: 26 Pages, Accepted for publication in ApJS

To investigate AGN outflows as a tracer of AGN feedback on star-formation, we perform integral-field spectroscopy of 20 type 2 AGNs at z10$^{41.5}$ erg/s, and exhibit strong outflow signatures in the [O III] kinematics. By decomposing the emission-line profile, we obtain the maps of the narrow and broad components of [O III] and H$\alpha$ lines, respectively. The broad components in both [O III] and H$\alpha$ represent the non-gravitational kinematics, i.e., gas outflows, while the narrow components, especially in H$\alpha$, represent the gravitational kinematics, i.e., rotational disk. By using the integrated spectra within the flux-weighted size of the narrow-line region, we estimate the energetics of the gas outflows. The ionized gas mass is 1.0-38.5$\times 10^5$ $M_{\odot}$, and the mean mass outflow rate is 4.6$\pm$4.3 $M_{\odot}$/yr, which is a factor of ~260 higher than the mean mass accretion rate 0.02$\pm$0.01 $M_{\odot}$/yr. The mean energy injection rate of the sample is 0.8$\pm$0.6% of the AGN bolometric luminosity, while the momentum flux is (5.4$\pm$3.6)$\times$ $L_{bol}$/c on average, except for two most kinematically energetic AGNs with low $L_{bol}$, which are possibly due to the dynamical timescale of the outflows. The estimated outflow energetics are consistent with the theoretical expectations for energy-conserving outflows from AGNs, yet we find no supporting evidence of instantaneous quenching of star formation due to the outflows., Comment: 25 pages, 27 figures, 5 tables; Accepted for publication in ApJ

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., 21 pages, 10 figures, accepted for publication in ApJS

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., 5 pages, 3 figures, accepted for publication in MNRASL

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., Comment: 12 pages, 12 figures. Websites are ready for downloading the source code (Odyssey) and the educational software (Odyssey_Edu). Accepted by ApJ

We study a gas rich merging dwarf system KUG 0200-096. Deep optical image reveals an optically faint tail with a length of 20 kpc, giving a visual impression of tidal antenna similar to NGC 4038/39. The interacting dwarf galaxies have B-band absolute magnitudes -18.06 and -16.63 mag. We identify a young stellar clump of stellar mass of 2$\times$10$^{7}$ M$_{\sun}$ at the tip of the antenna, possibly a Tidal Dwarf Galaxy (TDG). The putative TDG candidate is quite blue with $g-r$ color index of -0.07 mag, whereas the interacting dwarf galaxies have $g-r$ color indices 0.29 and 0.19 mag. The TDG is currently forming star at the rate of 0.02 M$_{\sun}$/yr. We obtained HI 21 cm line data of KUG 0200-096 using the GMRT to get a more detailed view of neutral hydrogen (HI) emission in interacting dwarf galaxies and its TDG. Evidence of merger between the dwarf galaxy pair is also presence in HI kinematics and morphology where we find the HI contents of interacting pair is disturbed, forming a extended tail toward the TDG. The HI velocity field shows strong gradient along the HI tidal tail extension. We present a comparative study between the Antennae galaxy, NGC 4038/39, and KUG 0200-096 in both optical and HI gas properties and discuss possible origin of KUG 0200-096 TDG., Accepted for publication in AJ

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., Comment: 13 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters

Strong Balmer absorption lines and the lack of Ha and [OII] emission lines signify that E+As are post-starburst systems. Recent studies suggest that E+As may undergo the transition from the `blue cloud' to the `red sequence' and eventually migrate to red sequence ETGs. An observational validation of this scenario is to identify the intervening galaxy population between E+As and the red-sequence. Motivated by recent findings with GALEX that a large fraction of ETGs exhibit UV-excess as a sign of RSF, we investigate the possible connection of the UV-excess galaxies to E+As. In particular, we examine the FP scaling relations of the largest sample of ~1,000 E+As selected from the SDSS and ~20,000 morphologically-selected SDSS ETGs with GALEX UV data. The FP parameters, combined with stellar population indicators, reveal a certain group of UV-excess ETGs that bridges between E+As and quiescent red galaxies. The newly identified galaxies are the post-starburst systems characterized by UV-excess but no Ha emission. This is a conceptual generalisation of "E+A", in that the Balmer absorption line in the "E+A" definition is replaced with UV-optical colours that are far more sensitive to RSF than the Balmer lines. We refer to these UV-excess galaxies as "E+a" galaxies, which stands for elliptical ("E") galaxies with a minority of A-type ("a") young stars. The species are either (1) galaxies that experienced starbursts weaker than those observed in E+As (1~10% of E+As, "mild E+As") or (2) the products of passively evolved E+As after quenching star formation quite a while ago (~1 Gyr, "old E+As"). We suggest that the latter type of E+a galaxies represents the most recent arrival to the red sequence in the final phase of the "E+A" to "red early-type" transition. (Abridged), Comment: 15 pages, 15 figures, Accepted for publication in MNRAS

Despite the efforts of the past decade, the origin of the bimodal horizontal-branch (HB) found in some globular clusters (GCs) remains a conundrum. Inspired by the discovery of multiple stellar populations in the {\it most massive} Galactic GC, $\omega$ Centauri, we investigate the possibility that two distinct populations may coexist and are responsible for the bimodal HBs in the {\it third} and {\it fifth} brightest GCs, NGC 6388 and NGC 6441. Using the population synthesis technique, we examine two different chemical ``self-enrichment'' hypotheses in which a primordial GC was sufficiently massive to contain two or more distinct populations as suggested by the populations found in $\omega$ Cen: (1) the age-metallicity relation scenario in which two populations with different metallicity and age coexist, following an internal age-metallicity relation, and (2) the super-helium-rich scenario in which GCs contain a certain fraction of helium-enhanced stars, for instance, the second generation stars formed from the helium-enriched ejecta of the first. The comparative study indicates that the detailed color-magnitude diagram morphologies and the properties of the RR Lyrae variables in NGC 6388 and NGC 6441 support the latter scenario; i.e., the model which assumes a minor fraction ($\sim$ 15 %) of helium-excess (Y $\simeq$ 0.3) stars. The results suggest that helium content is the main driver behind the HB bimodality found most often in massive GCs. If confirmed, the GC-to-GC variation of helium abundance should be considered a {\it local} effect, further supporting the argument that age is the {\it global} second parameter of HB morphology., Comment: Accepted in ApJ, 33 pages, 5 figures, 4 tables

We present ultraviolet (UV) photometry of M31 globular clusters (GCs) found in 23 Galaxy Evolution Explorer (GALEX) images covering the entirety of M31. We detect 485 and 273 GCs (and GC candidates) in the near-ultraviolet (NUV; 2267 A) and far-ultraviolet (FUV; 1516 A), respectively. Comparing M31 data with those of Galactic GCs in the UV with the aid of population models, we find that the age ranges of old GCs in M31 and the Galactic halo are similar. Three metal-rich ([Fe/H]>-1) GCs in M31 produce significant FUV flux making their FUV-V colors unusually blue for their metallicities. These are thought to be analogs of the two peculiar Galactic GCs NGC 6388 and NGC 6441 with extended blue HB stars. Based on the models incorporating helium enriched subpopulations in addition to the majority of the population that have a normal helium abundance, we suggest that even small fraction of super-helium-rich subpopulations in GCs can reproduce the observed UV bright metal-rich GCs. Young clusters in M31 show distinct UV and optical properties from GCs in Milky Way. Population models indicate that their typical age is less than ~ 2 Gyrs. A large fraction of young GCs have the kinematics of the thin, rapidly rotating disk component. However, a subset of the old GCs also shares the thin-disk kinematics of the younger clusters. The existence of young GCs on the outskirts of M31 disk suggests the occurrence of a significant recent star formation in the thin-disk of M31. Old thin-disk GCs may set constraints on the epoch of early formation of the M31 thin-disk. We detect 12 (10) intermediate-age GC candidates in NUV (FUV). We suggest that some of spectroscopically identified intermediate-age GCs may not be truly intermediate in age, but rather older GCs that possess developed HB., 43 pages, 14 figures, accepted for the GALEX special issue of ApJS

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., 10 pages, 10 figures, accepted for publication in MNRAS

Recent observations for the color-magnitude diagrams (CMDs) of the massive globular cluster Omega Centauri have shown that it has a striking double main sequence (MS), with a minority population of bluer and fainter MS well separated from a majority population of MS stars. Here we confirm, with the most up-to-date Y2 isochrones, that this special feature can only be reproduced by assuming a large variation (Delta Y = 0.15) of primordial helium abundance among several distinct populations in this cluster. We further show that the same helium enhancement required for this special feature on the MS can by itself reproduce the extreme horizontal-branch (HB) stars observed in Omega Cen, which are hotter than normal HB stars. Similarly, the complex features on the HBs of other globular clusters, such as NGC 2808, are explained by large internal variations of helium abundance. Supporting evidence for the helium-rich population is also provided by the far-UV (FUV) observations of extreme HB stars in these clusters, where the enhancement of helium can naturally explain the observed fainter FUV luminosity for these stars. The presence of super helium-rich populations in some globular clusters suggests that the third parameter, other than metallicity and age, also influences CMD morphology of these clusters., 4 pages, 4 figures, accepted for publication in the Astrophysical Journal Letters

The UV upturn is the rising flux with decreasing wavelength between the Lyman limit and 2500��found virtually in all bright spheroidal galaxies. It has been a mystery ever since it was first detected by the OAO-2 space telescope (Code & Welch 1979) because such old metal-rich populations were not expected to contain any substantial number of hot stars. It was confirmed by following space missions, ANS (de Boer 1982), IUE (Bertola et al. 1982) and HUT (Brown et al. 1997). The positive correlation between the UV-to-optical colour (i.e., the strength of the UV upturn) and the Mg2 line strength found by Burstein et al. (1987) through IUE observations has urged theorists to construct novel scenarios in which metal-rich ($\gtrsim Z_{\odot}$) old ($\gtrsim$ a few Gyr) stars become UV bright (Greggio & Renzini 1990; Horch et al. 1992). Also interesting was to find using HUT that, regardless of the UV strength, the UV spectral slopes at 1000--2000��in the six UV bright galaxies were nearly identical suggesting a very small range of temperatures of the UV sources in these galaxies (Brown et al. 1997), which corresponds to $T_{\rm eff} \approx 20,000 \pm 3,000$ K. This, together with other evidence, effectively ruled out young stars as the main driver of the UV upturn. A good review on the observational side of the story is given in the next article by Tom Brown, as well as in the recent articles of Greggio & Renzini (1999) and O'Connell (1999)., 6 figures; belated paper from Keele Conference

The discovery of multiple stellar populations in Milky Way globular clusters (GCs) has stimulated various follow-up studies on helium-enhanced stellar populations. Here we present the evolutionary population synthesis models for the spectro-photometric evolution of simple stellar populations (SSPs) with varying initial helium abundance ($Y_{\rm ini}$). We show that $Y_{\rm ini}$ brings about {dramatic} changes in spectro-photometric properties of SSPs. Like the normal-helium SSPs, the integrated spectro-photometric evolution of helium-enhanced SSPs is also dependent on metallicity and age for a given $Y_{\rm ini}$. {We discuss the implications and prospects for the helium-enhanced populations in relation to the second-generation populations found in the Milky Way GCs.} All of the models are available at \url{http://web.yonsei.ac.kr/cosmic/data/YEPS.htm}., Accepted for publication in ApJ, 21 pages, 12 figures, 3 tables

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., Comment: 16 pages, 13 figures, 3 tables. Accepted for publication in ApJ

The strength of the H$\beta$ index has been computed for the integrated spectra of model globular clusters from the evolutionary population synthesis. These models take into account, for the first time, the detailed systematic variation of horizontal-branch (HB) morphology with age and metallicity. Our models show that the H$\beta$ index is significantly affected by the presence of blue HB stars. Because of the contribution from blue HB stars, the H$\beta$ does not monotonically decrease as metallicity increases at a given age. Instead, it reaches a maximum strength when the distribution of HB stars is centered around 9500 K, the temperature where the H$\beta$ index becomes strongest. Our models indicate that the strength of the H$\beta$ index increases as much as 0.75 {\AA} due to the presence of blue HB stars. The comparison of the recent Keck observations of the globular cluster system in the Milky Way Galaxy with those in giant elliptical galaxies, NGC 1399 and M87, shows a systematic shift in the H$\beta$ against metallicity plane. Our models suggest that this systematic difference is understood if the globular cluster systems in giant elliptical galaxies are several billion years older, in the mean, than the Galactic counterpart. Further observations of globular cluster systems in the external galaxies from the large ground-based telescopes and space UV facilities will enable to clarify whether this difference is indeed due to the age difference or other explanations are also possible., Comment: 21 pages, 8 figures, Accepted for publication in the Astronomical Journal, August 2000 issue, Postscript files are available at http://csaweb.yonsei.ac.kr/~hclee/Hbeta

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., 12 pages, 3 figures, 1 table, accepted for publication in the ApJL

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., Comment: Accepted for publication in ApJ, 35 pages, 11 figures, 8 tables

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 ($\Delta$Y = 0.19). Our photometry, together with the results for other massive GCs (e.g., $\omega$ 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., Comment: 4 pages, 4 figures, 1 table, accepted for the publication in ApJL

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., 9 pages, 8 figures, 2 tables

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., 27 pages, 10 figures, 2 tables, accepted for publication in ApJ

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], 34 pages, 6 figures, 6 tables, Accepted for publication in ApJ

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., 55 pages, 7 figures, 5 tables; Accepted for publication in ApJ

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., 12 pages, 4 figures, 1 table, Accepted for publication in ApJL

We explore the age distribution of the globular cluster (GC) system of the nearby elliptical galaxy NGC 5128 using ultraviolet (UV) photometry from Galaxy Evolution Explorer (GALEX) observations, with UV - optical colors used as the age indicator. Most GCs in NGC 5128 follow the general trends of GCs in M31 and Milky Way in UV - optical color-color diagram, which indicates that the majority of GCs in NGC 5128 are old similar to the age range of old GCs in M31 and Milky Way. A large fraction of spectroscopically identified intermediate-age GC (IAGC) candidates with ~ 3-8 Gyr are not detected in the FUV passband. Considering the nature of intermediate-age populations being faint in the far-UV (FUV) passband, we suggest that many of the spectroscopically identified IAGCs may be truly intermediate in age. This is in contrast to the case of M31 where a large fraction of spectroscopically suggested IAGCs are detected in FUV and therefore may not be genuine IAGCs but rather older GCs with developed blue horizontal branch stars. Our UV photometry strengthens the results previously suggesting the presence of GC and stellar subpopulation with intermediate age in NGC 5128. The existence of IAGCs strongly indicates the occurrence of at least one more major star formation episode after a starburst at high redshift., 8 pages, 3 figures, accepted for ApJ Letter

In recent years, increasing evidence for chemical complexity and multiple stellar populations in massive globular clusters (GCs) has emerged, including extreme horizontal branches (EHBs) and UV excess. Our goal is to improve our understanding of UV excess in the regime of both massive GCs and ultra-compact dwarf galaxies (UCDs). To this end, we use deep archival GALEX data of the central Fornax cluster to measure NUV and FUV magnitudes of UCDs and massive GCs. We obtain NUV photometry for a sample of 35 compact objects with -13.5, Comment: Research Note, 7 pages, 3 figures, accepted for publication in A&A

We have studied ~2100 early-type galaxies in the SDSS DR3 which have been detected by the GALEX Medium Imaging Survey (MIS), in the redshift range 0 < z < 0.11. Combining GALEX UV photometry with corollary optical data from the SDSS, we find that, at a 95 percent confidence level, at least ~30 percent of galaxies in this sample have UV to optical colours consistent with some recent star formation within the last Gyr. In particular, galaxies with a NUV - r colour less than 5.5 are very likely to have experienced such recent star formation, taking into account the possibility of a contribution to NUV flux from the UV upturn phenomenon. We find quantitative agreement between the observations and the predictions of a semi-analytical LCDM hierarchical merger model and deduce that early-type galaxies in the redshift range 0 < z < 0.11 have ~1 to 3 percent of their stellar mass in stars less than 1 Gyr old. The average age of this recently formed population is ~300 to 500 Myrs. We also find that monolithically evolving galaxies, where recent star formation can be driven solely by recycled gas from stellar mass loss, cannot exhibit the blue colours (NUV - r < 5.5) seen in a significant fraction (~30 percent) of our observed sample., To appear in GALEX dedicated ApJ issue in December 2007 (accepted June 2006). A version with high-resolution figures can be downloaded at http://www-astro.physics.ox.ac.uk/~skaviraj/PAPERS/uv_cmr.pdf

We study the UV and optical properties of 38 massive GCs in M87, imaged using the HST/STIS and WFPC2 instruments. The majority of these GCs appear extremely bright in the FUV - roughly a magnitude brighter than their Galactic counterparts with similar metallicities. The observed FUV flux is several times larger than predictions of canonical old stellar population models. These canonical models, which assume a fiducial helium enrichment parameter, dY/dZ=2, are able to reproduce the observed FUV fluxes only if ages are 3--5 Gyr larger than the WMAP age of the Universe, although the same models fit the UV photometry of Galactic and M31 GCs for ages less than the WMAP age. A similar discrepancy is found between the mass-weighted and UV-luminosity weighted ages of the massive Galactic GC omega Cen, whose CMD (including peculiar features on its well-populated horizontal branch) can be accurately reproduced by invoking a small super-He-rich (dY/dZ ~ 90) stellar component. By comparison to omega Cen, we propose that the majority of M87 GCs in our sample contain strong signatures of similarly minor super-He-rich sub-components. Although we cannot prove or disprove the extreme helium scenario at the moment, we show that the same phenomenon that causes the extended horizontal branch of omega Cen explains the UV brightness of our sample. If this is indeed due to the extreme helium, this study would be the first to find its signatures in extragalactic objects., Comment: 19 pages, 10 figures, accepted for publication in MNRAS

The colors of globular clusters in most of large elliptical galaxies are bimodal. This is generally taken as evidence for the presence of two cluster subpopulations that have different geneses. However, here we find that, because of the non-linear nature of the metallicity-to-color transformation, a coeval group of old clusters with a unimodal metallicity spread can exhibit color bimodality. The models of cluster colors indicate that the horizontal-branch stars are the main drivers behind the empirical non-linearity. We show that the scenario gives simple and cohesive explanations for all the key observations, and could simplify theories of elliptical galaxy formation., Science, 311, 1129; Minor changes to text to match the published version (9 pages, 3 figures)

We present the results of our analysis of the photometric data of globular clusters in the elliptical galaxy NGC 5128 (Cen A). We show that the integrated colour U-B can be an effective metallicity indicator for simple stellar populations. This is because it is sensitive to metallicity via the opacity effect but relatively insensitive to the effective main sequence turn-off temperature of the population (and thus to age) when Teff ~ 7000-12000K, that is, when the Balmer Jump is strong. This flat U-B vs Teff relation is a result of the fact that the blueing effect of the optical continuum with increasing temperature is temporarily stopped by the Balmer Jump which becomes greater with temperature in this range. In this study we use U-B and B-V as metallicity and age indicators, respectively. We first show that the use of the U-B vs B-V two-colour diagram roughly yields the metallicities and ages of the Milky Way globular clusters independently determined, and then apply the technique to the clusters in NGC 5128. There is a large range in U-B, which corresponds to [Fe/H] of -2.0 through over +0.3. The large uncertainties from the data and the models currently prevent us from pinning down their ages and metallicities. Although a constant age for all these clusters cannot be ruled out, there is a hint of the metal-rich clusters being younger. Significance of these results and caveats of the analysis are discussed., 14 pages, 8 figures, Accepted MNRAS, Table 1 available on email request

We deduce the globular cluster formation history of the nearby elliptical NGC 5128 by using a phenomenological chemical enrichment model to reproduce its observed metallicity distribution function (MDF). We derive the observed MDF using recently obtained U and B photometry of the NGC 5128 GC system, with (U-B) used as the metallicity indicator. We test whether the observed MDF can be reproduced using a single starburst at high redshift followed by passive evolution (monolithic scenario). Our results strongly indicate that this scenario is incapable of reproducing the observed MDF of the NGC 5128 GC system. At least one more major starburst is required. A persistent feature of our best-fit models is a strong GC formation episode in the first 1-2 Gyrs, followed by a second formation episode typically after 10-12 Gyrs. Roughly one-third of the stellar mass in the GC system is created in this late starburst, producing clusters that are potentially less than 1 Gyr old. Although we do not have adequate resolution in the model or in the observed MDF to make a reliable claim about the possibility of additional GC formation episodes, our results suggest that the GCs in this galaxy did not form solely in a monolithic collapse and seem to support observational evidence which indicates that NGC 5128 is a merger remnant., Comment: replaced with accepted version

One of the long-standing problems in modern astronomy is the curious division of Galactic globular clusters, the "Oosterhoff dichotomy," according to the properties of their RR Lyrae stars. Here, we find that most of the lowest-metallicity ([Fe/H], Comment: 10 pages, 4 figures, published in Science