Source author record

Kenji Bekki

Kenji Bekki 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

85works
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

85 published item(s)

preprint2022arXiv

A Search for Annihilating Dark Matter in 47 Tucanae and Omega Centauri

A plausible formation scenario for the Galactic globular clusters 47 Tuc and Omega Cen is that they are tidally stripped remnants of dwarf galaxies, in which case they are likely to have retained a fraction of their dark matter cores. In this study, we have used the ultra-wide band receiver on the Parkes telescope (Murriyang) to place upper limits on the annihilation rate of exotic Light Dark Matter particles (X) via the XX/$e^+e^-$ channel using measurements of the recombination rate of positronium (Ps). This is an extension of a technique previously used to search for Ps in the Galactic Centre. However, by stacking of spectral data at multiple line frequencies, we have been able to improve sensitivity. Our measurements have resulted in 3-sigma flux density (recombination rate) upper limits of 1.7 mJy (1.4x10$^{43}$ s$^{-1}$) and 0.8 mJy (1.1x10$^{43}$ s$^{-1}$) for 47 Tuc and Omega Cen, respectively. Within the Parkes beam at the cluster distances, which varies from 10 to 23 pc depending on the frequency of the recombination line, we calculate upper limits to the dark matter mass and rms dark matter density of <1.2-1.3x10$^5$ fn$^{-0.5}$ (mX/MeV c$^{-2}$) solar masses and <48-54 fn$^{-0.5}$ (mX/MeV c$^{-2}$) solar masses pc$^{-3}$ for the clusters, where fn=Rn/Rp is the ratio of Ps recombination transitions to annihilations, estimated to be 0.01. The radio limits for Omega Cen suggest that, for a fiducial dark/luminous mass ratio of 0.05, any contribution from Light Dark Matter is small unless the cross section <7.9x10$^{-28}$ (mX/MeV c$^{-2}$)$^2$ cm$^3$ s$^{-1}$. Owing to the compactness and proximity of the clusters, archival 511-keV measurements suggest even tighter limits than permitted by CMB anisotropies, <8.6x10$^{-31}$ (mX/MeV c$^{-2}$)$^2$ cm$^3$ s$^{-1}$. Our recombination rate limits substantially improve on previous radio limits for the Milky Way.

preprint2022arXiv

The age gradients of galaxies in EAGLE: outside-in quenching as the origin of young bulges in cluster galaxies

Many disc galaxies in clusters have been found with bulges of similar age or younger than their surrounding discs, at odds with field galaxies of similar morphology and their expected inside-out formation. We use the EAGLE simulations to test potential origins for this difference in field and cluster galaxies. We find, in agreement with observations, that on average disc-dominated field galaxies in the simulations have older inner regions, while similar galaxies in groups and clusters have similarly aged or younger inner regions. This environmental difference is a result of outside-in quenching of the cluster galaxies. Prior to group/cluster infall, galaxies of a given present-day mass and morphology exhibit a similar evolution in their specific star formation rate (sSFR) profiles. Post-infall, the outer sSFRs of group and cluster galaxies significantly decrease due to interstellar medium stripping, while the central sSFR remains similar to field galaxies. Field disc galaxies instead generally retain radially increasing sSFR profiles. Thus, field galaxies continue to develop negative age gradients (younger discs), while cluster galaxies instead develop positive age gradients (younger bulges).

preprint2022arXiv

The present-day globular cluster kinematics of lenticular galaxies from the E-MOSAICS simulations and their relation to the galaxy assembly histories

We study the present-day rotational velocity ($V_{rot}$) and velocity dispersion ($σ$) profiles of the globular cluster (GC) systems in a sample of 50 lenticular (S0) galaxies from the E-MOSAICS galaxy formation simulations. We find that 82% of the galaxies have GCs that are rotating along the photometric major axis of the galaxy ($aligned$), while the remaining 18% of the galaxies do not ($misaligned$). This is generally consistent with the observations from the SLUGGS survey. For the $aligned$ galaxies, classified as $peaked$ $and$ $outwardly$ $decreasing$ (49%), $flat$ (24%) and $increasing$ (27%) based on the $V_{rot}/σ$ profiles out to large radii, we do not find any clear correlation between these present-day $V_{rot}/σ$ profiles of the GCs and the past merger histories of the S0 galaxies, unlike in previous simulations of galaxy stars. For just over half of the $misaligned$ galaxies, we find that the GC misalignment is the result of a major merger within the last 10 Gyr so that the $ex$-$situ$ GCs are misaligned by an angle between 0° (co-rotation) to 180° (counter-rotation) with respect to the $in$-$situ$ GCs, depending on the orbital configuration of the merging galaxies. For the remaining $misaligned$ galaxies, we suggest that the $in$-$situ$ metal-poor GCs, formed at early times, have undergone more frequent kinematic perturbations than the $in$-$situ$ metal-rich GCs. We also find that the GCs accreted early and the $in$-$situ$ GCs are predominantly located within 0.2 virial radii ($R_{200}$) from the centre of galaxies in 3D phase-space diagrams.

preprint2022arXiv

The VMC survey -- XLV. Proper motion of the outer LMC and the impact of the SMC

The Large Magellanic Cloud (LMC) is the most luminous satellite galaxy of the Milky Way and owing to its companion, the Small Magellanic Cloud (SMC), represents an excellent laboratory to study the interaction of dwarf galaxies. The aim of this study is to investigate the kinematics of the outer regions of the LMC by using stellar proper motions to understand the impact of interactions, e.g. with the SMC about 250 Myr ago. {We calculate proper motions using multi-epoch $K_\mathrm{s}$-band images from the VISTA survey of the Magellanic Clouds system (VMC). Observations span a time baseline of 2$-$5 yr. We combine the VMC data with data from the Gaia early Data Release 3 and introduce a new method to distinguish between Magellanic and Milky Way stars based on a machine learning algorithm. This new technique enables a larger and cleaner sample selection of fainter sources as it reaches below the red clump of the LMC. We investigate the impact of the SMC on the rotational field of the LMC and find hints of stripped SMC debris. The south east region of the LMC shows a slow rotational speed compared to the overall rotation. $N$-body simulations suggest that this could be caused by a fraction of stripped SMC stars, located in that particular region, that move opposite to the expected rotation.

preprint2021arXiv

Simulations of globular clusters within their parent galaxies: multiple stellar populations and internal kinematics

Using three-dimensional smoothed particle hydrodynamics simulations, we investigate the formation of multiple stellar populations (MSPs) in globular clusters (GCs) within the context of their parent galaxies. In our scenario, the second generation (2G) of stars originate from both asymptotic giant branch (AGB) polluters and pristine gas accreted from the host galaxy. Previous theoretical and numerical studies have demonstrated that this 'AGB with dilution' model has the potential to alleviate several problems faced by the classical AGB scenario. However, the accretion of pristine gas on to the GC has yet to be investigated within the context of the parent galaxy. This paper presents the preliminary results from our original simulation code which models GC formation from giant molecular clouds in a host galaxy, and subsequent gas accretion on to the GC. By simulating the genesis of the 2G over a 370 Myr time frame, we demonstrate that the fraction of 2G stars are inextricably linked to the GC's environment. Our simulations rationalize the wide variety of abundance patterns, kinematics, and 2G concentrations by altering the initial conditions of both the GC progenitor and the host galaxy itself. Most notably, we reproduce a positive correlation between the fraction of 2G stars and the initial mass of the cluster. We discuss the physical implications of our scenario and compare our simulations with observations of the Galactic GC 47 Tucanae (47 Tuc). Finally, we present scaling relations that encompass the wider GC population and serve as a reference for future observations.

preprint2020arXiv

A new way to constrain the densities of intra-group medium in groups of galaxies with convolutional neural networks

Ram pressure (RP) can influence the evolution of cold gas content and star formation rates of galaxies. One of the key parameters for the strength of RP is the density of intra-group medium ($ρ_{\rm igm}$), which is difficult to estimate if the X-ray emission from it is too weak to be observed. We propose a new way to constrain $ρ_{\rm igm}$ through an application of convolutional neural networks (CNNs) to simulated gas density and kinematic maps galaxies under strong RP. We train CNNs using $9\times{}10^4$ 2D images of galaxies under various RP conditions, then validate performance with $10^4$ new test images. This new method can be applied to real observational data from ongoing WALLABY and SKA surveys to quickly obtain estimates of $ρ_{\rm igm}$. Simulated galaxy images have $1.0$ kpc resolution, which is consistent with that expected from the future WALLABY survey. The trained CNN models predict the normalised IGM density, $\hatρ_{\rm igm}$ where $0.0 \le \hatρ_{\rm igm, n} < 10.0$, accurately with root mean squared error values ($\rm RMSE$) of $0.72$, $0.83$ and $0.74$ for the density, kinematic and joined 2D maps, respectively. Trained models are unable to predict the relative velocity of galaxies with respect to the IGM ($v_{\rm rel}$) precisely, and struggle to generalise for different RP conditions. We apply our CNNs to the observed HI column density map of NGC 1566 in the Dorado group to estimate its IGM density.

preprint2020arXiv

Bars formed in galaxy merging and their classification with deep learning

Stellar bars are a common morphological feature of spiral galaxies. While it is known that they can form in isolation, or be induced tidally, few studies have explored the production of stellar bars in galaxy merging. We look to investigate bar formation in galaxy merging using methods from deep learning to analyse our N-body simulations. The primary aim is to determine the constraints on the mass ratio and orientations of merging galaxies that are most conducive to bar formation. We further aim to explore whether it is possible to classify simulated barred spiral galaxies based on the mechanism of their formation. We test the feasibility of this new classification schema with simulated galaxies. Using a set of 29,400 images obtained from our simulations, we first trained a convolutional neural network to distinguish between barred and non-barred galaxies. We then tested the network on simulations with different mass ratios and spin angles. We adapted the core neural network architecture for use with our additional aims. We find that a strong inverse relationship between mass ratio and the number of bars produced. We also identify two distinct phases in the bar formation process; (1) the initial, tidally induced formation pre-merger, and (2) the destruction and/or regeneration of the during and after the merger. Mergers with low mass ratios and closely-aligned orientations are considerably more conducive to bar formation compared to equal-mass mergers. We demonstrate the flexibility of our deep learning approach by showing it is feasible to classify bars based on their formation mechanism.

preprint2020arXiv

Formation of massive black holes in ultra-compact dwarf galaxies: migration of primordial intermediate-mass black holes in N-body simulation

Recent observational studies of ultra-compact dwarf galaxies (UCDs) have discovered massive black holes (MBHs) with masses of more than ${10^6~\rm M_\odot}$, in their central regions. We here consider that these MBHs can be formed through merging of intermediate-mass black holes (IMBH), with masses of ${[10^3-10^5]~{\rm M}_{\odot}}$, within the stellar nuclei of dwarf galaxies, which are progenitors of UCDs. We numerically investigate this formation process for a wide range of model parameters using N-body simulations. This means that IMBH growth and feedback is neglected in this study. We find that only massive IMBHs of $10^5~\rm M_\odot$ sink into the central regions of their host dwarf ($\approx 10^{10}~\rm M_\odot$) to be gravitationally trapped by its stellar nucleus within less than 1 Gyr in most dwarf models. We also find that lighter IMBHs with $[1 - 30] \times 10^3~\rm M_\odot$ sink into the centre in low-mass dwarfs ($\approx 10^{9}~\rm M_\odot$) due to more efficient dynamical friction (DF). Additionally, we show that the IMBHs can form binaries in the centre and, rarely, before they reach the centre, which may lead to the IMBHs merging and thus emitting gravitational waves that could be detected by LISA. Finally, we discuss the required number of IMBHs for the MBH formation in UCDs and the physical roles of stellar nuclei in IMBH binaries and mergers.

preprint2020arXiv

Formation of massive globular clusters with dark matter and its implication on dark matter annihilation

Recent observational studies of $γ$-ray emission from massive globular clusters (GCs) have revealed possible evidence of dark matter (DM) annihilation within GCs. It is, however, still controversial whether the emission comes from DM or from milli-second pulsars. We here present the new results of numerical simulations, which demonstrate that GCs with DM can originate from nucleated dwarfs orbiting the ancient MW. The simulated stripped nuclei (i.e., GCs) have the central DM densities ranging from 0.1 to several ${\rm M_\odot pc^{-3}}$, depending on the orbits and the masses of the host dwarf galaxies. However, GCs born outside the central regions of their hosts can have no/little DM after their hosts are destroyed and the GCs become the Galactic halo GCs. These results suggest that only GCs originating from stellar nuclei of dwarfs can possibly have DM. We further calculate the expected $γ$-ray emission from these simulated GCs and compare them to observations of $ω$ Cen. Given the large range of DM densities in the simulated GCs, we suggest that the recent possible detection of DM annihilation from GCs should be more carefully interpreted.

preprint2020arXiv

Photoelectric heating effects on the evolution of luminous disk galaxies

Photoelectric heating (PEH) influences the temperature and density of the interstellar medium (ISM), and potentially also affecting star formation. PEH is expected to have a stronger effect on massive galaxies, as they host larger dust reservoirs compared to dwarf systems. Accordingly, in this paper, we study PEH effects in Milky Way-like galaxies using smoothed particle hydrodynamics (SPH) code which self-consistently implements the evolution of the gas, dust, and interstellar radiation field (ISRF). Dust evolution includes dust formation by stars, destruction by SNe, and growth in dense media. We find that PEH suppresses star formation due to the excess heating that reduces the ISM density. This suppression is seen across the entire range of gas fractions, star formation recipes, dust models, and PEH efficiencies investigated by our code. The suppression ranges from negligible values to approximately a factor of five depending on the specific implementation. Galaxy models having higher gas fraction experience higher star formation suppression. The adopted dust model also alters the extent of star formation suppression. Moreover, when PEH is switched on, galaxy models show higher gas outflow rates and have higher loading factors indicative of enhanced SNe feedback. In gas-rich models (i.e. a gas fraction of 0.5), we also find that PEH suppresses the formation of disk clumps via violent disk instabilities, and thus suppresses bulge formation via clumps migration to the central regions.

preprint2020arXiv

The assembly history of the nearest S0 galaxy NGC 3115 from its kinematics out to six half-light radii

Using new and archival data, we study the kinematic properties of the nearest field S0 galaxy, NGC 3115, out to $\sim6.5$ half-light radii ($R_\mathrm{e}$) from its stars (integrated starlight), globular clusters (GCs) and planetary nebulae (PNe). We find evidence of three kinematic regions with an inner transition at $\sim0.2\ R_\mathrm{e}$ from a dispersion-dominated bulge ($V_\mathrm{rot}/σ<1$) to a fast-rotating disk ($V_\mathrm{rot}/σ>1$), and then an additional transition from the disk to a slowly rotating spheroid at $\sim2-2.5\, R_\mathrm{e}$, as traced by the red GCs and PNe (and possibly by the blue GCs beyond $\sim5\, R_\mathrm{e}$). From comparison with simulations, we propose an assembly history in which the original progenitor spiral galaxy undergoes a gas-rich minor merger that results in the embedded kinematically cold disk that we see today in NGC 3115. At a later stage, dwarf galaxies, in mini mergers (mass-ratio $<$ 1:10), were accreted building-up the outer slowly rotating spheroid, with the central disk kinematics largely unaltered. Additionally, we report new spectroscopic observations of a sample of ultra-compact dwarfs (UCDs) around NGC 3115 with the Keck/KCWI instrument. We find that five UCDs are inconsistent with the general rotation field of the GCs, suggesting an \textit{ex-situ} origin for these objects, i.e. perhaps the remnants of tidally stripped dwarfs. A further seven UCDs follow the GC rotation pattern, suggesting an \textit{in-situ} origin and, possibly a GC-like nature.

preprint2020arXiv

The role of dust destruction and dust growth in the evolution of the interstellar medium

We use Milky Way-like chemodynamical simulations with a new treatment for dust destruction and growth to investigate how these two processes affect the properties of the interstellar medium in galaxies. We focus on the role of two specific parameters: f_des (a new parameter that determines the fraction of dust destroyed in a single gas particle surrounding supernova) and C_s (the probability that a metal atom or ion sticks to the dust grain after colliding, i.e., the sticking coefficient) in regulating the amount and distribution of dust, cold gas and metals in galaxies. We find that simulated galaxies with low f_des and/or high C_s values produce not only more dust, but they also have a shallower correlation between dust surface density and total gas surface density, and a steeper correlation between dust-to-gas ratio and metallicity. Only for values of f_des between 0.01 and 0.02, and of C_s between 0.5 and 1 our simulations produce an average slope of the dust-to-gas ratio versus metallicity relation consistent with observations. f_des values correspond to a range of a total fraction of dust destroyed by a single supernova between 0.42 and 0.44. Lastly, we compare predictions of several simulations (with different star formation recipes, gas fractions, central metallicities, and metallicity gradients) to the spatially resolved M101 galaxy, and conclude that metallicity is the primary driver of the spatial distribution of dust, while dust-to-gas ratio controls the cold gas distribution, as it regulates the atomic-to-molecular hydrogen conversion rate.

preprint2020arXiv

The SAMI Galaxy Survey: A Range in S0 Properties Indicating Multiple Formation Pathways

It has been proposed that S0 galaxies are either fading spirals or the result of galaxy mergers. The relative contribution of each pathway, and the environments in which they occur remains unknown. Here we investigate stellar and gas kinematics of 219 S0s in the SAMI Survey to look for signs of multiple formation pathways occurring across the full range of environments. We identify a large range of rotational support in their stellar kinematics, which correspond to ranges in their physical structure. We find that pressure-supported S0s with $v/σ$ below 0.5 tend to be more compact and feature misaligned stellar and gas components, suggesting an external origin for their gas. We postulate that these S0s are consistent with being formed through a merger process. Meanwhile, comparisons of ellipticity, stellar mass and Sérsic index distributions with spiral galaxies shows that the rotationally supported S0s with $v/σ$ above 0.5 are more consistent with a faded spiral origin. In addition, a simulated merger pathway involving a compact elliptical and gas-rich satellite results in an S0 that lies within the pressure-supported group. We conclude that two S0 formation pathways are active, with mergers dominating in isolated galaxies and small groups, and the faded spiral pathway being most prominent in large groups ($10^{13} < M_{halo} < 10^{14}$).

preprint2020arXiv

The Sami Galaxy Survey: stellar populations of passive spiral galaxies in different environment

We investigate the stellar populations of passive spiral galaxies as a function of mass and environment, using integral field spectroscopy data from the Sydney-AAO Multi-object Integral field spectrograph Galaxy Survey. Our sample consists of $52$ cluster passive spirals and $18$ group/field passive spirals, as well as a set of S0s used as a control sample. The age and [Z/H] estimated by measuring Lick absorption line strength indices both at the center and within $1R_{\rm e}$ do not show a significant difference between the cluster and the field/group passive spirals. However, the field/group passive spirals with log(M$_\star$/M$_\odot)\gtrsim10.5$ show decreasing [$α$/Fe] along with stellar mass, which is $\sim0.1$ dex smaller than that of the cluster passive spirals. We also compare the stellar populations of passive spirals with S0s. In the clusters, we find that passive spirals show slightly younger age and lower [$α$/Fe] than the S0s over the whole mass range. In the field/group, stellar populations show a similar trend between passive spirals and S0s. In particular, [$α$/Fe] of the field/group S0s tend to be flattening with increasing mass above log(M$_\star$/M$_\odot)\gtrsim10.5$, similar to the field/group passive spirals. We relate the age and [$α$/Fe] of passive spirals to their mean infall time in phase-space; we find a positive correlation, in agreement with the prediction of numerical simulations. We discuss the environmental processes that can explain the observed trends. The results lead us to conclude that the formation of the passive spirals and their transformation into S0s may significantly depend on their environments.

preprint2020arXiv

The VMC survey -- XXXIX: Mapping metallicity trends in the Small Magellanic Cloud using near-infrared passbands

We have derived high spatial resolution metallicity maps covering $\sim$42 deg$^2$ across the Small Magellanic Cloud (SMC) in an attempt to understand its metallicity distribution and gradients up to a radius of $\sim$ 4$^{\circ}$. Using the near-infrared VISTA Survey of the Magellanic Clouds, our data cover a thrice larger area compared with previous studies. We identify red giant branch (RGB) stars in spatially distinct $Y, (Y-K_{\rm s})$ colour--magnitude diagrams. In any of our selected subregions, the RGB slope is used as an indicator of the average metallicity, based on calibration to metallicity using spectroscopic data. The metallicity distribution across the SMC is unimodal and can be fitted by a Gaussian distribution with a peak at [Fe/H] = $-$0.97 dex ($σ$[Fe/H] = 0.05 dex). We find evidence of a shallow gradient in metallicity ($-0.031 \pm 0.005$ dex deg$^{-1}$) from the galactic centre to radii of 2$^{\circ}$--2.5$^{\circ}$, followed by a flat metallicity trend from $\sim$ 3.5$^{\circ}$ to 4$^{\circ}$. We find that the SMC's metallicity gradient is radially asymmetric. It is flatter towards the East than to the West, hinting at mixing and/or distortion of the spatial metallicity distribution (within the inner 3$^{\circ}$), presumably caused by tidal interactions between the Magellanic Clouds.

preprint2020arXiv

The VMC survey -- XXXVIII. Proper motion of the Magellanic Bridge

The Magellanic Clouds are a nearby pair of interacting dwarf galaxies and satellites of the Milky Way. Studying their kinematic properties is essential to understanding their origin and dynamical evolution. They have prominent tidal features and the kinematics of these features can give hints about the formation of tidal dwarfs, galaxy merging and the stripping of gas. In addition they are an example of dwarf galaxies that are in the process of merging with a massive galaxy. The goal of this study is to investigate the kinematics of the Magellanic Bridge, a tidal feature connecting the Magellanic Clouds, using stellar proper motions to understand their most recent interaction. We calculated proper motions based on multi-epoch $K_{s}$-band aperture photometry, which were obtained with the Visible and Infrared Survey Telescope for Astronomy (VISTA), spanning a time of 1-3 yr, and we compared them with $Gaia$ Data Release 2 (DR2) proper motions. We tested two methods for removing Milky Way foreground stars using $Gaia$~DR2 parallaxes in combination with VISTA photometry or using distances based on Bayesian inference. We obtained proper motions for a total of 576,411 unique sources over an area of $23$ deg$^{2}$ covering the Magellanic Bridge including mainly Milky Way foreground stars, background galaxies, and a small population of possible Magellanic Bridge stars ($<$15,000). The first proper motion measurement of the Magellanic Bridge centre is $1.80\pm0.25$ mas yr$^{-1}$ in right ascension and $-0.72\pm0.13$ mas yr$^{-1}$ in declination. The proper motion measurements confirm a flow motion from the Small to the Large Magellanic Cloud. This flow can now be measured all across the entire length of the Magellanic Bridge. Our measurements indicate that the Magellanic Bridge is stretching.

preprint2019arXiv

A new formation model for omega Centauri: the crossroad of astrophysical processes

We investigate the formation processes of the Galactic globular cluster (GC) omega Cen with multiple stellar populations based on our original hydrodynamical simulations with chemical enrichment by Type II supernovae (SNe II), asymptotic giant branch (AGB) stars, and neutron star mergers (NSMs). The principal results are as follows. Multiple stellar populations with a wide range of [Fe/H] can be formed from rather massive and compact molecular cloud with a mass of 2 * 10^7 M_sun in the central region of its dwarf galaxy within less than a few hundred Myr. Gas ejected from SNe II and AGB stars can mix well to form new stars with higher He abundances (Y) and higher [Fe/H]. The He-rich stars are strongly concentrated in the GC's central region so that the GC can show a steep negative gradient of Y. Relative ratios of light elements to Fe show bimodal distributions for a given [Fe/H] owing to star formation from original gas and AGB ejecta. [La/Fe] and [Ba/Fe] can rapidly increase until [Fe/H]~-1.5 and then decrease owing to Fe ejection from SNe II. Although AGB ejecta can be almost fully retained in intra-cluster medium, NSM ejecta can be retained only partially. This difference in the retention capability is responsible for the observed unique [Eu/Fe]-[Fe/H] and [La/Eu]-[Fe/H] relations in omega Cen. Some observational results such as the [O/Na]$-$[Fe/H] relation and radial [Fe/H] gradient are yet to be well reproduced in the present model.

preprint2016arXiv

Formation of anomalous globular clusters with metallicity spreads: A unified picture

Recent observations have revealed that at least 8 globular clusters (GCs) in the Galaxy show internal abundance spreads in [Fe/H]. We investigate the origin of these `anomalous' GCs using numerical simulations of GCs in the dwarfs orbiting around the Galaxy and chemical evolution model of dwarfs hosting the GCs. The principal results are as follows. GCs formed in a host dwarf galaxy with a total mass of ~ 10^10 M_sun can merge to form a single nuclear GC before the host is completely destroyed by the Galaxy, if they are massive (> 3*10^5 M_sun) and if they are formed in the inner region (R<400 pc). The GC merger remnants can capture field stars during its spiral-in to nuclear regions. If two GCs are formed from star formation events separated by ~300 Myr in their host dwarf, then the new GC formed from GC merging can have [Fe/H] spread of 0.2 dex and [Ba/Fe] spread of 0.3 dex. GCs formed from GC merging can show variety of internal abundance spreads depending on the details of their hosts' chemical evolution. We suggest that anomalous GCs were formed from GC merging that occurred before the destruction of GC host dwarfs yet after self-enrichment processes responsible for the observed anti-correlations between chemical abundances of light elements. We also suggest that the observed no/little dependence of [Eu/Fe] on [Fe/H] in the Galactic GC M22 is evidence of massive dwarf galaxies hosting these anomalous GCs.

preprint2016arXiv

Radio Astrometry towards the Nearby Universe with the SKA

This chapter summarizes radio astrometry in relation to current very long baseline interferometry (VLBI) projects and describes its perspectives with the SKA. The scientific goals of the astrometry with the SKA have been discussed in the international and Japanese communities of researchers, whose major issues are shown here. We have demonstrated some of the issues, such as censuses of possible targets and the technical feasibility of astrometry in the SKA frequency bands. The preliminary results of our case studies on SKA astrometry are also presented. In addition, possible synergy and commensality of the SKA astrometric projects with those in the optical and infrared astrometric missions, especially JASMINE (Japan Astrometry Satellite Mission for INfrared Exploration) are discussed.

preprint2016arXiv

Significant enhancement of ${\rm H_2}$ formation in disk galaxies under strong ram pressure

We show, for the first time, that ${\rm H_2}$ formation on dust grains can be enhanced in disk galaxies under strong ram-pressure (RP). We numerically investigate how the time evolution, of ${\rm H}$ {\sc i} and ${\rm H_2}$ components in disk galaxies orbiting a group/cluster of galaxies, can be influenced by hydrodynamical interaction between the gaseous components of the galaxies and the hot intra-cluster medium (ICM). We find that compression of ${\rm H}$ {\sc i} caused by RP increases ${\rm H_2}$ formation in disk galaxies, before RP rapidly strips ${\rm H}$ {\sc i}, cutting off the fuel supply and causing a drop in ${\rm H_2}$ density. We also find that the level of this ${\rm H_2}$ formation enhancement in a disk galaxy under RP depends on the mass of its host cluster dark matter (DM) halo, initial positions and velocities of the disk galaxy, and disk inclination angle with respect to the orbital plane. We demonstrate that dust growth is a key factor in the evolution of the ${\rm H}$ {\sc i} and ${\rm H_2}$ mass in disk galaxies under strong RP. We discuss how the correlation between ${\rm H_2}$ fractions and surface gas densities of disk galaxies evolves with time in the galaxies under RP. We also discuss whether or not galaxy-wide star formation rates (SFRs) in cluster disk galaxies can be enhanced by RP if the SFRs depend on ${\rm H_2}$ densities.

preprint2016arXiv

The SAMI galaxy survey: Galaxy Interactions and Kinematic Anomalies in Abell 119

Galaxy mergers are important events that can determine the fate of a galaxy by changing its morphology, star-formation activity and mass growth. Merger systems have commonly been identified from their disturbed morphologies, and we now can employ Integral Field Spectroscopy to detect and analyze the impact of mergers on stellar kinematics as well. We visually classified galaxy morphology using deep images ($μ_{\rm r} = 28\,\rm mag\,\, arcsec^{-2}$) taken by the Blanco 4-m telescope at the Cerro Tololo Inter-American Observatory. In this paper we investigate 63 bright ($M_{\rm r}<-19.3$) spectroscopically-selected galaxies in Abell 119; of which 53 are early type and 20 galaxies show a disturbed morphology by visual inspection. A misalignment between the major axes in the photometric image and the kinematic map is conspicuous in morphologically-disturbed galaxies. Our sample is dominated by early-type galaxies, yet it shows a surprisingly tight Tully-Fisher relation except for the morphologically-disturbed galaxies which show large deviations. Three out of the eight slow rotators in our sample are morphology disturbed. The visually-selected morphologically-disturbed galaxies are generally more asymmetric, visually as well as kinematically. Our findings suggest that galaxy interactions, including mergers and perhaps fly-bys, play an important role in determining the orientation and magnitude of galaxy's angular momentum.

preprint2016arXiv

The VMC Survey - XXI. New star cluster candidates discovered from infrared photometry in the Small Magellanic Cloud

We report the first search for new star clusters performed using the VISTA near-infrared YJKs Magellanic Clouds survey (VMC) data sets. We chose a pilot field of ~ 0.4 deg^2 located in the South-West of the Small Magelllanic Cloud (SMC) bar, where the star field is among the densest and highest reddened region in the galaxy. In order to devise an appropriate automatic procedure we made use of dimensions and stellar densities observed in the VMC data sets of the known clusters in this area. We executed different kernel density estimations over a sample of more than 358000 stars with magnitudes measured in the three YJKs filters. We analysed the new cluster candidates whose colour-magnitude diagrams (CMDs), cleaned from field star contamination, were used to assess the clusters' reality and estimate reddenings and ages of the genuine systems. As a result 38 objects (~ a 55% increase in the known star clusters located in the surveyed field) of 0.15 - 0.40 arcmin (2.6 - 7.0 pc) in radius resulted to have near-infrared CMD features which resemble those of star clusters of young to moderate intermediate age (log(t /yr ) ~ 7.5-9.0). Most of the new star cluster candidates are hardly recognizable in optical images without the help of a sound star field decontaminated CMD analysis. For highly reddened star cluster candidates (E(B - V ) >= 0.6 mag) the VMC data sets were necessary in order to recognize them.

preprint2015arXiv

Dust-regulated galaxy formation and evolution:A new chemodynamical model with live dust particles

Interstellar dust plays decisive roles in the conversion of neutral to molecular hydrogen (H_2), the thermodynamical evolution of interstellar medium (ISM), and the modification of spectral energy distributions (SEDs) of galaxies. These important roles of dust have not been self-consistently included in previous numerical simulations of galaxy formation and evolution. We have therefore developed a new model by which one can investigate whether and how galaxy formation and evolution can be influenced by dust-related physical processes such as photo-electric heating, H_2 formation on dust, and stellar radiation pressure on dust in detail. A novel point of the model is that different dust species in a galaxy are represented by `live dust' particles (i.e., not test particles). Therefore, dust particles in a galaxy not only interact gravitationally with all four components of the galaxy (i.e., dark matter, stars, gas, and dust) but also are grown and destroyed through physical processes of ISM. First we describe a way to include dust-related physical processes in Nbody+hydrodynamical simulations of galaxy evolution in detail. Then we show some preliminary results of dust-regulated galaxy evolution. The preliminary results suggest that the evolution of dust distributions driven by radiation pressure of stars is very important for the evolution of star formation rates, chemical abundances, H_2 fractions, and gas distributions in galaxies.

preprint2015arXiv

Formation of emission line dots and extremely metal-deficient dwarfs from almost dark galaxies

Recent observations have discovered a number of extremely gas-rich very faint dwarf galaxies possibly embedded in low-mass dark matter halos. We investigate star formation histories of these gas-rich dwarf ("almost dark") galaxies both for isolated and interacting/merging cases. We find that although star formation rates (SFRs) are very low (<10^-5 M_sun/yr) in the simulated dwarfs in isolation for the total halo masses (M_h) of 10^8-10^9 M_sun, they can be dramatically increased to be ~ 10^{-4} M_sun/yr when they interact or merge with other dwarfs. These interacting faint dwarfs with central compact HII regions can be identified as isolated emission line dots ("ELdots") owing to their very low surface brightness envelopes of old stars. The remnant of these interacting and merging dwarfs can finally develop central compact stellar systems with very low metallicities (Z/Z_sun<0.1), which can be identified as extremely metal-deficient ("XMD") dwarfs. These results imply that although there would exist many faint dwarfs that can be hardly detected in the current optical observations, they can be detected as isolated ELdots or XMD dwarfs, when they interact with other galaxies and their host environments. We predict that nucleated ultra-faint dwarfs formed from merging of almost dark dwarfs can be identified as globular clusters owing to the very low surface brightness stellar envelopes, if they are distant objects.

preprint2015arXiv

Formation of ultra-compact blue dwarf galaxies and their evolution into nucleated dwarfs

We propose that there is an evolutionary link between ultra-compact blue dwarf galaxies (UCBDs) with active star formation and nucleated dwarfs based on the results of numerical simulations of dwarf-dwarf merging. We consider the observational fact that low-mass dwarfs can be very gas-rich, and thereby investigate the dynamical and chemical evolution of very gas-rich, dissipative dwarf-dwarf mergers. We find that the remnants of dwarf-dwarf mergers can be dominated by new stellar populations formed from the triggered starbursts and consequently can have blue colors and higher metallicities (Z~[0.2-1]Z_sun). We also find that the remnants of these mergers can have rather high mass-densities (10^4 M_sun pc^-3) within the central 10 pc and small half-light radii (40-100 pc). The radial stellar structures of some merger remnants are similar to those of nucleated dwarfs. Star formation can continue in nuclear gas disks (R<100 pc) surrounding stellar galactic nuclei (SGNs) so that the SGNs can finally have multiple stellar populations with different ages and metallicities. These very compact blue remnants can be identified as UCBDs soon after merging and as nucleated dwarfs after fading of young stars. We discuss these results in the context of the origins of metal-rich ultra-compact dwarfs (UCDs) and SGNs.

preprint2015arXiv

Kinematics of dwarf galaxies in gas-rich groups, and the survival and detectability of tidal dwarf galaxies

We present DEIMOS multi-object spectroscopy (MOS) of 22 star-forming dwarf galaxies located in four gas-rich groups, including six newly-discovered dwarfs. Two of the galaxies are strong tidal dwarf galaxy (TDG) candidates based on our luminosity-metallicity relation definition. We model the rotation curves of these galaxies. Our sample shows low mass-to-light ratios (M/L=0.73$\pm0.39M_\odot/L_\odot$) as expected for young, star-forming dwarfs. One of the galaxies in our sample has an apparently strongly-falling rotation curve, reaching zero rotational velocity outside the turnover radius of $r_{turn}=1.2r_e$. This may be 1) a polar ring galaxy, with a tilted bar within a face-on disk; 2) a kinematic warp. These scenarios are indistinguishable with our current data due to limitations of slit alignment inherent to MOS-mode observations. We consider whether TDGs can be detected based on their tidal radius, beyond which tidal stripping removes kinematic tracers such as H$α$ emission. When the tidal radius is less than about twice the turnover radius, the expected falling rotation curve cannot be reliably measured. This is problematic for as much as half of our sample, and indeed more generally, galaxies in groups like these. Further to this, the H$α$ light that remains must be sufficiently bright to be detected; this is only the case for three (14%) galaxies in our sample. We conclude that the falling rotation curves expected of tidal dwarf galaxies are intrinsically difficult to detect.

preprint2015arXiv

New quasars behind the Magellanic Clouds. Spectroscopic confirmation of near-infrared selected candidates

Quasi--stellar objects (quasars) located behind nearby galaxies provide an excellent absolute reference system for astrometric studies, but they are difficult to identify because of fore- and background contamination. Deep wide--field, high angular resolution surveys spanning the entire area of nearby galaxies are needed to obtain a complete census of such quasars. We embarked on a program to expand the quasar reference system behind the Large and the Small Magellanic Clouds, the Magellanic Bridge, and the Magellanic Stream, connecting the Clouds with the Milky Way. Hundreds of quasar candidates were selected based on their near--infrared colors and variability properties from the ongoing public ESO VISTA Magellanic Clouds survey. A subset of 49 objects was followed up with optical spectroscopy. We confirmed the quasar nature of 37 objects (34 new identifications), four are low redshift objects, three are probably stars, and the remaining three lack prominent spectral features for a secure classification; bona fide quasars, judging from their broad absorption lines are located, as follows: 10 behind the LMC, 13 behind the SMC, and 14 behind the Bridge. The quasars span a redshift range from z~0.5 to z~4.1. Upon completion the VMC survey is expected to yield a total of ~1500 quasars with Y<19.32 mag, J<19.09 mag, and Ks<18.04 mag.

preprint2015arXiv

The Next Generation Virgo Cluster Survey. VI. The Kinematics of Ultra-compact Dwarfs and Globular Clusters in M87

The origin of ultra-compact dwarfs (UCDs)--objects larger and more massive than typical globular clusters (GCs), but more compact than typical dwarf galaxies--has been hotly debated in the 15 years since their discovery. Even whether UCDs should be considered galactic in origin, or simply the most extreme GCs, is not yet settled. We present the dynamical properties of 97 spectroscopically confirmed UCDs (rh >~10 pc) and 911 GCs associated with central cD galaxy of the Virgo cluster, M87. Our UCDs, of which 89% have M_star > ~2X10^6 M_sun and 92% are as blue as the classic blue GCs, nearly triple the sample of previous confirmed Virgo UCDs, providing by far the best opportunity for studying the global dynamics of a UCD system. We found that (1) UCDs have a surface number density profile that is shallower than that of the blue GCs in the inner ~ 70 kpc and as steep as that of the red GCs at larger radii; (2) UCDs exhibit a significantly stronger rotation than the GCs, and the blue GCs seem to have a velocity field that is more consistent with that of the surrounding dwarf ellipticals than with that of UCDs; (3) UCDs have a radially increasing orbital anisotropy profile, and are tangentially-biased at radii < ~ 40 kpc and radially-biased further out. In contrast, the blue GCs become more tangentially-biased at larger radii beyond ~ 40 kpc; (4) GCs with M_star > 2X10^6 M_sun have rotational properties indistinguishable from the less massive ones, suggesting that it is the size, instead of mass, that differentiates UCDs from GCs as kinematically distinct populations. We conclude that most UCDs in M87 are not consistent with being merely the most luminous and extended examples of otherwise normal GCs. The radially-biased orbital structure of UCDs at large radii is in general agreement with the "tidally threshed dwarf galaxy" scenario.

preprint2015arXiv

The origin of dust extinction curves with or without the 2175 A bump in galaxies: The case of the Magellanic Clouds

The Large and Small Magellanic Clouds (LMC and SMC, respectively) are observed to have characteristic dust extinction curves that are quite different from those of the Galaxy (e.g., strength of the 2175 A bump). Although the dust composition and size distribution of the Magellanic Clouds (MCs) that can self-consistently explain their observed extinction curves have been already proposed, it remain unclear whether and how the required dust properties can be achieved in the formation histories of the MCs. We therefore investigate the time evolution of the dust properties of the MCs and thereby derive their extinction curves using one-zone chemical evolution models with formation and evolution of small and large silicate and carbonaceous dust grains and dusty winds associated with starburst events. We find that the observed SMC extinction curve without a conspicuous 2175 A bump can be reproduced well by our SMC model, if the small carbon grains can be selectively lost through the dust wind during the latest starburst about 0.2 Gyr ago. We also find that the LMC extinction curve with a weak 2175 A bump can be reproduced by our LMC model with less efficient removal of dust through dust wind. We discuss possible physical reasons for different dust wind efficiencies between silicate and graphite and among galaxies.

preprint2015arXiv

The transformation and quenching of simulated gas-rich dwarf satellites within a group environment

The underlying mechanisms driving the quenching of dwarf-mass satellite galaxies remain poorly constrained, but recent studies suggest they are particularly inefficient for those satellites with stellar mass 10$^{\rm 9}$ M$_{\odot}$. We investigate the characteristic evolution of these systems with chemodynamical simulations and idealised models of their tidal/hydrodynamic interactions within the 10$^{\rm 13-13.5}$ M$_{\odot}$ group-mass hosts in which they are preferentially quenched. Our fiducial simulations highlight the role played by secular star formation and stellar bars, and demonstrate a transition from a gas-rich to passive, HI-deficient state (i.e. $Δ$SFR$\le$-1, def$_{\rm HI}$$\ge$0.5) within 6 Gyr of first infall. Furthermore, in the 8-10 Gyr in which these systems have typically been resident within group hosts, the bulge-to-total ratio of an initially bulgeless disc can increase to 0.3$<$B/T$<$0.4, its specific angular momentum $λ_{\rm R}$ reduce to $\sim$0.5, and strong bisymmetries formed. Ultimately, this scenario yields satellites resembling dwarf S0s, a result that holds for a variety of infall inclinations/harassments albeit with broad scatter. The key assumptions here lie in the rapid removal of the satellite's gaseous halo upon virial infall, and the satellite's local intra-group medium density being defined by the host's spherically-averaged profile. We demonstrate how quenching can be greatly enhanced if the satellite lies in an overdensity, consistent with recent cosmological-scale simulations but contrasting with observationally-inferred quenching mechanisms/timescales; an appraisal of these results with respect to the apparent preferential formation of dS0s/S0s in groups is also given.

preprint2015arXiv

The VMC Survey - XIV. First results on the look-back time star-formation rate tomography of the Small Magellanic Cloud

We analyse deep images from the VISTA survey of the Magellanic Clouds in the YJKs filters, covering 14 sqrdeg (10 tiles), split into 120 subregions, and comprising the main body and Wing of the Small Magellanic Cloud (SMC). We apply a colour--magnitude diagram reconstruction method that returns their best-fitting star formation rate SFR(t), age-metallicity relation (AMR), distance and mean reddening, together with 68% confidence intervals. The distance data can be approximated by a plane tilted in the East-West direction with a mean inclination of 39 deg, although deviations of up to 3 kpc suggest a distorted and warped disk. After assigning to every observed star a probability of belonging to a given age-metallicity interval, we build high-resolution population maps. These dramatically reveal the flocculent nature of the young star-forming regions and the nearly smooth features traced by older stellar generations. They document the formation of the SMC Wing at ages <0.2 Gyr and the peak of star formation in the SMC Bar at 40 Myr. We clearly detect periods of enhanced star formation at 1.5 Gyr and 5 Gyr. The former is possibly related to a new feature found in the AMR, which suggests ingestion of metal-poor gas at ages slightly larger than 1 Gyr. The latter constitutes a major period of stellar mass formation. We confirm that the SFR(t) was moderately low at even older ages.

preprint2015arXiv

The VMC survey - XVI. Spatial variation of the cluster-formation activity in the innermost regions of the Large Magellanic Cloud

We present results based on $YJK_{\rm s}$ photometry of star clusters in the Large Magellanic Cloud (LMC), distributed throughout the central part of the galaxy's bar and the 30 Doradus region. We analysed the field-star decontaminated colour--magnitude diagrams of 313 clusters to estimate their reddening values and ages. The clusters are affected by a mean reddening of $E(B-V) \in [0.2,0.3]$ mag, where the average internal LMC reddening amounts to $\sim$ 0.1--0.2 mag. The region covering 30 Doradus includes clusters with reddening values in excess of $E(B-V)$ = 0.4 mag. Our cluster sample spans the age range $7.0 \le \log(t$ yr$^{-1}) < 9.0$, represents an increase of 30 per cent in terms of the number of clusters with robust age estimates and comprises a statistically complete sample in the LMC regions of interest here. The resulting cluster frequencies suggest that the outermost regions of the LMC bar first experienced enhanced cluster formation -- $\log(t$ yr$^{-1}) \in [8.5,9.0]$ -- before the activity proceeded, although in a patchy manner, to the innermost regions, for $\log(t$ yr$^{-1}) < 7.7$. Cluster frequencies in the 30 Doradus region show that the area is dominated by very recent cluster formation. The derived star-formation frequencies suggest that the cluster and field-star populations do not seem to have fully evolved as fully coupled systems during the last $\sim$ 100 Myr.

preprint2015arXiv

The VMC survey - XVII. The proper motion of the Small Magellanic Cloud and of the Milky Way globular cluster 47 Tucanae

In this study we use multi-epoch near-infrared observations from the VISTA survey of the Magellanic Cloud system (VMC) to measure the proper motion of different stellar populations in a tile of 1.5 deg sq. in size in the direction of the Galactic globular cluster 47 Tuc. We obtain the proper motion of the cluster itself, of the Small Magellanic Cloud (SMC), and of the field Milky Way stars. Stars of the three main stellar components are selected from their spatial distribution and their distribution in colour-magnitude diagrams. Their average coordinate displacement is computed from the difference between multiple Ks-band observations for stars as faint as Ks=19 mag. Proper motions are derived from the slope of the best-fitting line among 10 VMC epochs over a time baseline of ~1 yr. Background galaxies are used to calibrate the absolute astrometric reference frame. The resulting absolute proper motion of 47 Tuc is (mu_alpha cos(delta), mu_delta)=(+7.26+/-0.03, -1.25+/-0.03) mas/yr. This measurement refers to about 35000 sources distributed between 10 and 60 arcmin from the cluster centre. For the SMC we obtain (mu_alpha cos(delta), mu_delta)=(+1.16+/-0.07, -0.81+/-0.07) mas/yr from about 5250 red clump and red giant branch stars. The absolute proper motion of the Milky Way population in the line-of-sight (l =305.9, b =-44.9) of this VISTA tile is (mu_alpha cos(delta), mu_delta)=(+10.22+/-0.14, -1.27+/-0.12) mas/yr and results from about 4000 sources. Systematic uncertainties associated to the astrometric reference system are 0.18 mas/yr. Thanks to the proper motion we detect 47 Tuc stars beyond its tidal radius.

preprint2014arXiv

Chemical evolution of galaxies with radiation-driven dust wind

We discuss how the removal of interstellar dust by radiation pressure of stars influences the chemical evolution of galaxies by using a new one-zone chemical evolution models with dust wind. The removal efficiency of an element (e.g., Fe, Mg, and Ca) through radiation-driven dust wind in a galaxy is assumed to depend both on the dust depletion level of the element in interstellar medium and the total luminosity of the galaxy in the new model. We particularly focus on the time evolution of [alpha/Fe] and its dependence on model parameters for dust wind in this study. The principal results are as follows. The time evolution of [Ca/Fe] is significantly different between models with and without dust wind in the sense that [Ca/Fe] can be systematically lower in the models with dust wind. The time evolution of [Mg/Fe], on the other hand, can not be so different between the models with and without dust wind owing to the lower level of dust depletion for Mg. As a result of this, [Mg/Ca] can be systematically higher in the models with dust wind. We compare these results with the observed elemental features of stars in the Large Magellanic Cloud (LMC), because a growing number of observational studies on [alpha/Fe] for the LMC have been recently accumulated for a detailed comparison. Based on the present new results, we also discuss the origins of [alpha/Fe] in the Fornax dwarf galaxy and elliptical galaxies in the context of radiation-driven dust wind.

preprint2014arXiv

Choirs HI Galaxy Groups: The metallicity of dwarf galaxies

We present a recalibration of the luminosity-metallicity relation for gas-rich, star-forming dwarfs to magnitudes as faint as M$_R$ ~ -13. We use the Dopita et al. (2013) metallicity calibrations to calibrate the relation for all of the data in this analysis. In metallicity-luminosity space we find two sub-populations within a sample of high-confidence SDSS DR8 star-forming galaxies; 52% are metal-rich giants and 48% are metal-medium galaxies. Metal-rich dwarfs classified as tidal dwarf galaxy (TDG) candidates in the literature are typically of metallicity 12 + log(O/H) = 8.70 $\pm$ 0.05, while SDSS dwarfs fainter than M$_R$ = -16 have a mean metallicity of 12 + log(O/H) = 8.28 $\pm$ 0.10, regardless of their luminosity, indicating that there is an approximate floor to the metallicity of low luminosity galaxies. Our hydrodynamical simulations predict that TDGs should have metallicities elevated above the normal luminosity-metallicity relation. Metallicity can therefore be a useful diagnostic for identifying TDG candidate populations in the absence of tidal tails. At magnitudes brighter than M$_R$ ~ -16 our sample of 53 star-forming galaxies in 9 HI gas-rich groups is consistent with the normal relation defined by the SDSS sample. At fainter magnitudes there is an increase in dispersion in metallicity of our sample, suggestive of a wide range of HI content and environment. In our sample we identify three (16% of dwarfs) strong TDG candidates (12 + log(O/H) > 8.6), and four (21%) very metal poor dwarfs (12 + log(O/H) < 8.0), which are likely gas-rich dwarfs with recently ignited star formation.

preprint2014arXiv

Cosmic evolution of dust in galaxies: Methods and preliminary results

We investigate the redshift (z) evolution of dust properties, its dependences on initial conditions of galaxy formation, and physical correlations between dust, gas, and stellar contents at different z based on our original chemodynamical simulations of galaxy formation with dust growth and destruction. In this preliminary investigation, we first determine the reasonable ranges of the most important two parameters for dust evolution, i.e., the timescales of dust growth and destruction, by comparing the observed and simulated dust properties and molecular hydrogen H2 content of the Galaxy. We then investigate the z-evolution of dust-to-gas-ratios (D) and, H2 gas fraction (f_H2), and gas-phase chemical abundances (e.g., A_O=12+log(O/H)) in the simulated disk and dwarf galaxies. The principal results are as follows. Both D and f_H2 can rapidly increase during the early dissipative formation of galactic disks (z ~ 2-3) and the z-evolution of these depends on initial mass densities, spin parameters, and masses of galaxies. The observed A_O-D relation can be qualitatively reproduced, but the simulated dispersion of D at a given A_O is smaller. The simulated galaxies with larger total dust masses show larger H2 and stellar masses and higher f_H2. Disk galaxies show negative radial gradients of D and the gradients are steeper for more massive galaxies. Both dust-to-metal ratios and gas-phase [S/Fe] can be significantly different within a single galaxy and between different galaxies at different z, which means that fixed dust-to-metal ratios should not be used in investigating H2 contents and spectral energy distributions of galaxies.

preprint2014arXiv

Evolution of dust and molecular hydrogen in the Magellanic System

We investigate the evolution of the interstellar medium (ISM) in self-consistent, chemodynamical simulations of the Magellanic Clouds (MCs) during their recent (z<0.3) past. An explicit modelling of dust and molecular hydrogen lifecycles enables the comparison of our models against the observed properties of the ISM, including elemental depletion from the gas-phase. Combining this model with a tidal-dominated paradigm for the formation for the Magellanic Stream and Bridge, we reproduce the age-metallicity relations, long gas depletion timescales, and presently observed dust and molecular hydrogen masses of the MCs to within their respective uncertainties. We find that these models' enrichment depends sensitively on the processing of dust within the ISM and the dynamical influence of external tides/stellar bars. The ratio of characteristic dust destruction timescales in our SMC and LMC models, a governing parameter of our models' evolution, is consistent with estimates based on observed supernova (SN) rates. Our reference MC models tend to exhibit the disputed universal dust-to-metal ratio, which we argue stems from the adoption of high SNe II condensation efficiencies. Our models are the first to reproduce the one-tenth solar metallicity of the Stream/Leading Arm following tidal stripping of the SMC; the hypothesis that the LMC contributes a metal-rich filament to the Stream, as implied by recent kinematic and abundance analyses, is also appraised in this study.

preprint2014arXiv

The VMC Survey. XI. Radial Stellar Population Gradients in the Galactic Globular Cluster 47 Tucanae

We present a deep near-infrared color--magnitude diagram of the Galactic globular cluster 47 Tucanae, obtained with the Visible and Infrared Survey Telescope for Astronomy (VISTA) as part of the VISTA near-infrared Y, J, K_s survey of the Magellanic System (VMC). The cluster stars comprising both the subgiant and red-giant branches exhibit apparent, continuous variations in color--magnitude space as a function of radius. Subgiant-branch stars at larger radii are systematically brighter than their counterparts closer to the cluster core; similarly, red-giant-branch stars in the cluster's periphery are bluer than their more centrally located cousins. The observations can very well be described by adopting an age spread of ~0.5 Gyr as well as radial gradients in both the cluster's helium abundance (Y) and metallicity (Z), which change gradually from Y = 0.28, Z = 0.005 in the cluster core to Y = 0.25, Z = 0.003 in its periphery. We conclude that the cluster's inner regions host a significant fraction of second-generation stars, which decreases with increasing radius; the stellar population in the 47 Tuc periphery is well approximated by a simple stellar population.

preprint2014arXiv

Tidal-induced lopsidedness in Magellanic-type galaxies

We investigate the tidally-induced conversion of barred late-type spirals to Magellanic-type discs with numerical simulations, to establish how the lifetime of lopsidedness (asymmetry) varies with numerical parametrizations. Using a reference model based on observed properties of the Large Magellanic Cloud (LMC), we show that its off-centre bar, one-arm spiral and one-sided star formation can be induced by a brief tidal interaction. We thereafter perform a detailed parameter study, and find that at the dynamical mass of LMC-type discs ($\sim$10$^{\rm 10}$ \md{}), stellar lopsidedness (as quantified by the m$=$1 Fourier mode) and bar off-centredness can vary widely in amplitude, but are generally short-lived ($\sim$Gyr). Tidal interactions induce more persistent lopsidedness in lower mass galaxies (several Gyr), in particular those with large halo-to-disc mass ratios as implied by recent halo occupation models. We suggest that the tidal interactions play a larger role in the observed ubiquity of lopsidedness than the presently favoured gas accretion mechanism for lower mass galaxies. Other characteristics of Magellanic-type galaxies, such as one-sided star formation, are quantified and tend to arise more prominently in discs with later-type spiral structure (more halo dominated inner disc, weaker bars) following retrograde orbital encounters.

preprint2013arXiv

Choirs, HI Galaxy Groups: Catalogue and Detection of Star-forming Dwarf Group Members

Hα observations centred on galaxies selected from the HI Parkes All Sky Survey (HIPASS, Barnes et al. 2001) typically show one and sometimes two star-forming galaxies within the approximately 15-arcminute beam of the Parkes 64-m HI detections. In our Survey of Ionization in Neutral Gas Galaxies (SINGG, Meurer et al. 2006) we found fifteen cases of HIPASS sources containing four or more emission line galaxies (ELGs). We name these fields Choir groups. In the most extreme case we found a field with at least nine ELGs. In this paper we present a catalogue of Choir group members in the context of the wider SINGG sample. The dwarf galaxies in the Choir groups would not be individually detectable in HIPASS at the observed distances if they were isolated, but are detected in SINGG narrow-band imaging due to their membership of groups with sufficiently large total HI mass. The ELGs in these groups are similar to the wider SINGG sample in terms of size, Hα equivalent width, and surface brightness. Eight of these groups have two large spiral galaxies with several dwarf galaxies and may be thought of as morphological analogues of the Local Group. However, on average our groups are not significantly HI-deficient, suggesting that they are at an early stage of assembly, and more like the M81 group. The Choir groups are very compact at typically only 190 kpc in projected distance between the two brightest members. They are very similar to SINGG fields in terms of star formation efficiency (the ratio of star formation rate to HI mass; SFE), showing an increasing trend in SFE with stellar mass.

preprint2013arXiv

Coevolution of dust, gas, and stars in galaxies - I. Spatial distributions and scaling-relations of dust and molecular hydrogen

We investigate the time evolution of dust properties, molecular hydrogen (H_2) contents, and star formation histories in galaxies by using our original chemodynamical simulations. The simulations include the formation of dust in the stellar winds of supernovae (SNe) and asymptotic giant branch (AGB) stars, the growth and destruction processes of dust in the interstellar medium (ISM), the formation of polycyclic aromatic hydrocarbon (PAH) dust in carbon-rich AGB stars, the H_2 formation on dust grains, and the H_2 photo-dissociation due to far ultra-violet (FUV) light in a self-consistent manner. We focus mainly on disk galaxies with different total masses in this preliminary study. The principle results are as follows: The star formation histories of disk galaxies can be regulated by the time evolution of interstellar dust, mainly because the formation rates of H_2 can be controlled by dust properties. The observed correlation between dust-to-gas-ratios (D) and gas-phase oxygen abundances (A_O = 12+log(O/H)) can be reproduced reasonably well in the present models. The disks show negative radial gradients (i.e., larger in inner regions) of H_2 fraction (f_H2), PAH-to-dust mass ratio (f_PAH), D, and A_O and these gradients evolve with time. The surface-mass densities of dust (Sigma_dust) are correlated more strongly with the total surface gas densities than with those of H_2. Local gaseous regions with higher D are more likely to have higher f_H2 in individual disks and total H_2 masses correlate well with total dust masses. More massive disk galaxies are more likely to have higher D, f_PAH, and f_H2 and smaller dust-to-stellar mass ratios. We also compare between galactic star formation histories in the metallicity-dependent and dust-dependent star formation models and find no major differences.

preprint2013arXiv

Evolution of the high-mass end of the stellar initial mass functions in starburst galaxies

We investigate the time evolution and spatial variation of the stellar initial mass function (IMF) in star-forming disk galaxies by using chemodynamical simulations with an IMF model depending both on local densities and metallicities ([Fe/H]) of the interstellar medium (ISM). We find that the slope (alpha) of a power-law IMF (N(m) ~ m^-alpha) for stellar masses larger than 1M_sun evolves from the canonical Salpeter IMF (alpha ~ 2.35) to be moderately top-heavy one (alpha ~ 1.9) in the simulated disk galaxies with starbursts triggered by galaxy interaction. We also find that alpha in star-forming regions correlates with star formation rate densities (Sigma_SFR in units of M_sun yr^{-1} kpc^{-2}). Feedback effects of Type Ia and II supernovae are found to prevent IMFs from being too top-heavy (alpha < 1.5). The simulation predicts alpha ~ 0.23 log Sigma_SFR + 1.7 for log Sigma_SFR > -2 (i.e., more top-heavy in higher Sigma_SFR), which is reasonably consistent well with corresponding recent observational results. The present study also predicts that inner regions of starburst disk galaxies have smaller alpha thus are more top-heavy (d alpha/d R ~ 0.07 kpc^{-1} for R < 5 kpc). The predicted radial alpha gradient can be tested against future observational studies of the alpha variation in star-forming galaxies.

preprint2013arXiv

Galactic star formation enhanced and quenched by ram pressure in groups and clusters

We investigate how ram pressure of intragroup and intracluster medium can influence the spatial and temporal variations of star formation (SF) of disk galaxies with halo masses (M_h) ranging from 10^10 M_sun to 10^12 M_sun (i.e., from dwarf irregular to Milky Way-type) in groups and clusters with 10^13 < M_h/M_sun < 10^15 by using numerical simulations with a new model for time-varying ram pressure. The long-term evolution of SF rates and H_alpha morphologies corresponding to the distributions of star-forming regions are particularly investigated for different model parameters. The principal results are as follows. Whether ram pressure can enhance or reduce SF depends on M_h of disk galaxies and inclination angles of gas disks with respect to their orbital directions for a given orbit and a given environment. For example, SF can be moderately enhanced in disk galaxies with M_h=10^12 M_sun at the pericenter passages in a cluster with M_h=10^14 M_sun whereas it can be completely shut down (`quenching') for low-mass disks with M_h=10^10 M_sun. Ram pressure can reduce the H_alpha-to-optical-disk-size ratios of disks and the revel of the reduction depends on M_h and orbits of disk galaxies for a given environment. Disk galaxies under strong ram pressure show characteristic H_alpha morphologies such as ring-like, one-sided, and crescent-like distributions.

preprint2013arXiv

Origin of a bottom-heavy stellar initial mass function in elliptical galaxies

We investigate the origin of a bottom-heavy stellar initial mass function (IMF) recently observed in elliptical galaxies by using chemical evolution models with a non-universal IMF. We adopt the Kroupa IMF with the three slopes (alpha_1, alpha_2, and alpha_3) dependent on metallicities ([Fe/H]) and densities (rho_g) of star-forming gas clouds and thereby search for the best IMF model that can reproduce (i) the observed steep IMF slope (alpha_2 ~ 3, i.e., bottom-heavy) for low stellar masses (m < 1 M_sun) and (ii) the correlation of alpha_2 with chemical properties of elliptical galaxies in a self-consistent manner. We find that if the IMF slope alpha_2 depends both on [Fe/H] and rho_g, then elliptical galaxies with higher [Mg/Fe] can have steeper alpha_2 (~3) in our models. We also find that the observed positive correlation of stellar mass-to-light ratios (M/L) with [Mg/Fe] in elliptical galaxies can be quantitatively reproduced in our models with alpha_2 = beta[Fe/H]}+gamma log rho_g, where beta ~ 0.5 and gamma ~ 2. We discuss whether the IMF slopes for low-mass (alpha_2) and high-mass stars (alpha_3) need to vary independently from each other to explain a number of IMF-related observational results self-consistently. We also briefly discuss why alpha_2 depends differently on [Fe/H] in dwarf and giant elliptical galaxies.

preprint2013arXiv

Simulating galaxy evolution with a non-universal stellar initial mass function

We consider that the stellar initial mass function (IMF) depends on physical properties of star-forming molecular clouds in galaxies and thereby investigate how such a non-universal IMF (NUIMF) influences galaxy evolution. We incorporate a NUIMF model into galaxy-scale chemodynamical simulations in order to investigate the differences in chemical and dynamical evolution of disk galaxies between the NUIMF and universal IMF (UIMF) models. In the adopted NUIMF model, the three slopes of the Kroupa IMF depend independently on densities and metallicities ([Fe/H]) of molecular gas clouds, and production rates of metals and dust from massive and AGB stars can vary according to the time evolution of the three IMF slopes. The preliminary results of the simulations are as follows. Star formation rates in actively star-forming disk galaxies can be significantly lower in the NUIMF model than in the UIMF model, and the differences between the two models can be more remarkable in galaxies with higher SFRs. Chemical enrichment can proceed more rapidly in the NUIMF model and [Mg/Fe] for a given metallicity is higher in the NUIMF model. The evolution of H2 fraction (f_H2) and dust-to-gas ratio (D) is more rapid in the NUIMF model so that the final f_H2 and $D$ can be higher in the NUIMF model. Formation of massive stellar clumps in gas-rich disks is more strongly suppressed owing to the stronger SN feedback effect in the NUIMF model. The radial density profiles of new stars within the central 1kpc are shallower in the NUIMF model.

preprint2013arXiv

Stripping of nitrogen-rich AGB ejecta from interacting dwarf irregular galaxies

Dwarf irregular galaxies (dIrrs) including the Magellanic Clouds in the local Universe, in many cases, exhibit an unusually low N/O abundance ratio (log N/O ~ -1.5) in H II regions as compared with the solar value (~-0.9). This ratio is broadly equivalent to the average level of extremely metal-poor stars in the Galactic halo, suggesting that N released from asymptotic giant branch (AGB) stars is missing in the present-day interstellar matter of these dIrrs. We find evidence for past tidal interactions in the properties of individual dIrrs exhibiting low N/O ratios, while a clear signature of interactions is unseen for dIrrs with high N/O ratios. Accordingly, we propose that the ejecta of massive AGB stars that correspond to a major production site of N can be stripped from dIrrs that have undergone a strong interaction with a luminous galaxy. The physical process of its stripping is made up of two stages: (i) the ejecta of massive AGB stars in a dIrr are first merged with those of the bursting prompt SNe Ia and pushed up together to the galaxy halo, and (ii) subsequently through tidal interactions with a luminous galaxy, these ejecta are stripped from a dwarf galaxy's potential well. Our new chemical evolution models with stripping of AGB ejecta succeed in reproducing the observed low N/O ratio. Furthermore, we perform N-body + hydrodynamical simulations to trace the fate of AGB ejecta inside a dIrr orbiting the Milky Way, and confirm that a tidal interaction is responsible for the efficient stripping of AGB ejecta from dIrrs.

preprint2012arXiv

Chemical evolution of the Large Magellanic Cloud

We adopt a new chemical evolution model for the Large Magellanic Cloud (LMC) and thereby investigate its past star formation and chemical enrichment histories. The delay time distribution of type Ia supernovae recently revealed by type Ia supernova surveys is incorporated self-consistently into the new model. The principle results are summarized as follows. The present gas mass fraction and stellar metallicity as well as the higher [Ba/Fe] in metal-poor stars at [Fe/H]<-1.5 can be more self-consistently explained by models with steeper initial mass functions. The observed higher [Mg/Fe] (> 0.3) at [Fe/H] ~ -0.6 and higher [Ba/Fe] (>0.5) at [Fe/H] ~ -0.3 can be due to significantly enhanced star formation about 2 Gyr ago. The observed overall [Ca/Fe]-[Fe/H] relation and remarkably low [Ca/Fe] (<-0.2) at [Fe/H]>-0.6 are consistent with models with short-delay supernova Ia and with the more efficient loss of Ca possibly caused by an explosion mechanism of type II supernovae. Although the metallicity distribution functions do not show double peaks in the models with a starburst about 2 Gyr ago, they show characteristic double peaks in the models with double starbursts at ~200 Myr and ~2 Gyr ago. The observed apparent dip of [Fe/H] around ~1.5 Gyr ago in the age--metallicity relation can be reproduced by models in which a large amount (~10^9 M_{sun}) of metal-poor ([Fe/H]<-1) gas can be accreted onto the LMC.

preprint2012arXiv

Feedback effects of aspherical supernovae explosions on galaxies

We investigate how explosions of aspherical supernovae (A-SNe) can influence star formation histories and chemical evolution of dwarf galaxies by using a new chemodynamical model. We mainly present the numerical results of two comparative models so that the A-SN feedback effects on galaxies can be more clearly seen. SNe originating from stars with masses larger than 30M_sun are A-SNe in the "ASN" model whereas all SNe are spherical ones (S-SNe) in the "SSN" model. Each S-SN and A-SN are assumed to release feedback energy of 10^{51} erg and 10^{52} erg, respectively, and chemical yields and feedback energy of A-SN ejecta depend on angles between the axis of symmetry and the ejection directions. We find that star formation can become at least by a factor of ~3 lower in the ASN model in comparison with the SSN one owing to the more energetic feedback of A-SNe. As a result of this, chemical evolution can proceed very slowly in the ASN model. A-SN feedback effects can play a significant role in the formation of giant gaseous holes and energetic gaseous outflow and unique chemical abundances (e.g., high [Mg/Ca]). Based on these results, we provide a number of implications of the A-SN feedback effects on galaxy formation and evolution.

preprint2012arXiv

First Evidence of Globular Cluster Formation from the Ejecta of Prompt Type Ia Supernovae

Recent spectroscopic observations of globular clusters (GCs) in the Large Magellanic Cloud (LMC) have discovered that one of the intermediate-age GC, NGC 1718 with [Fe/H]=-0.7 has an extremely low [Mg/Fe] ratio of ~-0.9. We propose that NGC 1718 was formed from the ejecta of type Ia supernovae (SNe Ia) mixed with very metal-poor ([Fe/H] <-1.3) gas about ~ 2 Gyr ago. The proposed scenario is shown to be consistent with the observed abundances of Fe-group elements such as Cr, Mn, and Ni. In addition, compelling evidence for asymptotic giant branch stars playing a role in chemical enrichment during this GC formation is found. We suggest that the origin of the metal-poor gas is closely associated with the efficient gas-transfer from the outer gas disk of the Small Magellanic Cloud to the LMC disk. We anticipate that the outer part of the LMC disk contains field stars exhibiting significantly low [Mg/Fe] ratios, formed through the same process as NGC 1718.

preprint2012arXiv

The influences of the Magellanic Clouds on the Galaxy: Pole shift, warp, and star formation history

We investigate how the Large Magellanic Cloud (LMC) influences the evolution of the Galaxy after the LMC enters into the virial radius of the dark matter halo of the Galaxy for the first time. Both the Galaxy and the LMC are modeled as N-body particles in our models so that the dynamical influences of the LMC on the Galaxy can be investigated in a fully self-consistent manner. Furthermore, the orbital parameters for the LMC are carefully chosen such that the present location of the LMC in the Galaxy can be rather precisely reproduced in our simulations. We particularly investigate the influences of the LMC on the precession rate, the outer stellar and gaseous structures, and the star formation history of the Galaxy. Our principals results are summarized as follows. The LMC-Galaxy dynamical interaction can cause "pole shift" (or irregular precession/nutation) of the Galaxy and the typical rate of pole shift (dot θ_{d}) is ~2 degree Gyr^{-1} corresponding to ~ 7 muas yr^{-1}. The LMC-Galaxy interaction induces the formation of the outer warp structures of the Galaxy, which thus confirms the results of previous numerical simulations on the formation of the Galactic warp. The LMC-Galaxy interaction also induces the formation of outer (R>20 kpc) spiral arms and increases the vertical velocity dispersion of the outer disk significantly. The mean star formation rate of the Galaxy for the last several Gyrs can be hardly influenced by the LMC's tidal force. The age and metallicity distribution of stars in the solar-neighborhood (7 kpc < R < 10 kpc) for the last several Gyr can be only slightly changed by the past LMC-Galaxy interaction. Based on these results, we discuss how the possible ongoing Galactic pole shift with 10 muas yr^{-1} can be detected by future observational studies by GAIA.

preprint2012arXiv

The Tidal Origin of the Magellanic Stream and the Possibility of a Stellar Counterpart

We present an N-body model that reproduces the morphology and kinematics of the Magellanic Stream (MS), a vast neutral hydrogen (HI) structure that trails behind the Large and Small Magellanic Clouds (LMC and SMC, respectively) in their orbit about the Milky Way. After investigating 8 million possible orbits consistent with the latest proper motions, we adopt an orbital history in which the LMC and SMC have only recently become a strongly interacting binary pair. We find that their first close encounter 2 Gyr ago provides the necessary tidal forces to disrupt the disk of the SMC and thereby create the MS. The model also reproduces the on-sky bifurcation of the two filaments of the MS, and we suggest that a bound association with the Milky Way is required to reproduce the bifurcation. Additional HI structures are created during the tidal evolution of the SMC disk, including the Magellanic Bridge, the "Counter-Bridge", and two branches of leading material. Insights into the chemical evolution of the LMC are also provided, as a substantial fraction of the material stripped away from the SMC is engulfed by the LMC. Lastly we compare three different N-body realizations of the stellar component of the SMC, which we model as a pressure-supported spheroid motivated by recent kinematical observations. We find that an extended spheroid is better able to explain the stellar periphery of the SMC, and the tidal evolution of the spheroid may imply the existence of a stellar stream akin to the gaseous MS.

preprint2012arXiv

Two-Component Galactic Bulge Probed with Renewed Galactic Chemical Evolution Model

Results of recent observations of the Galactic bulge demand that we discard a simple picture of its formation, suggesting the presence of two stellar populations represented by two peaks of stellar metallicity distribution (MDF) in the bulge. To assess this issue, we construct Galactic chemical evolution models that have been updated in two respects: First, the delay time distribution (DTD) of type Ia supernovae (SNe Ia) recently revealed by extensive SN Ia surveys is incorporated into the models. Second, the nucleosynthesis clock, the s-processing in asymptotic giant branch (AGB) stars, is carefully considered in this study. This novel model first shows that the Galaxy feature tagged by the key elements, Mg, Fe, Ba for the bulge as well as thin and thick disks is compatible with a short-delay SN Ia. We present a successful modeling of a two-component bulge including the MDF and the evolutions of [Mg/Fe] and [Ba/Mg], and reveal its origin as follows. A metal-poor component (<[Fe/H]>~-0.5) is formed with a relatively short timescale of ~1 Gyr. These properties are identical to the thick disk's characteristics in the solar vicinity. Subsequently from its remaining gas mixed with a gas flow from the disk outside the bulge, a metal-rich component (<[Fe/H]>~+0.3) is formed with a longer timescale (~4 Gyr) together with a top-heavy initial mass function that might be identified with the thin disk component within the bulge.

preprint2011arXiv

A possible physical connection between helium-rich stellar populations of massive globular clusters and the UV upturn of galactic spheroids

We discuss a possible physical connection between helium-rich (Y > 0.35) stellar populations of massive globular clusters (GCs) and the ultraviolet (UV) upturn of galactic spheroids by using analytical and numerical models. In our model, all stars are initially formed as bound or unbound star clusters (SCs) formed from giant molecular clouds (GMCs) and the SCs can finally become GCs, open clusters, and field stars depending on physical properties of their host GMCs. An essential ingredient of the model is that helium-rich stars are formed almost purely from gas ejected from massive asymptotic giant branch (AGB) stars. The helium-rich star formation is assumed to occur within massive SCs if the masses of the progenitor GMCs are larger than a threshold mass (M_thres). These massive SCs can finally become either massive GCs or helium-rich field stars depending on whether they are disintegrated or not. Using this model, we show that if the initial mass functions (IMFs) in galactic spheroids are mildly top-heavy, then the mass fractions of helium-rich main-sequence stars (F_He) can be as large as ~ 0.1 for M_thres=10^7 M_sun. F_He is found to depend on IMFs and M_thres such that it can be larger for shallower IMFs and smaller M_thres. The inner regions of galactic spheroids show larger F_He in almost all models. Based on these results, we suggest that if the UV upturn of elliptical galaxies is due to the larger fractions of helium-rich stars, then the origin can be closely associated with top-heavy IMFs in the galaxies.

preprint2011arXiv

Constraining the orbital history of the Magellanic Clouds: A new bound scenario suggested by the tidal origin of the Magellanic Stream

Bound orbits have traditionally been assigned to the Large and Small Magellanic Clouds (LMC and SMC, respectively) in order to provide a formation scenario for the Magellanic Stream (MS) and its Leading Arm (LA), two prominent neutral hydrogen (HI) features connected to the LMC and SMC. However, Hubble Space Telescope (HST) measurements of the proper motions of the LMC and SMC have challenged the plausibility of bound orbits, causing the origin of the MS to re-emerge as a contested issue. We present a new tidal model in which structures resembling the bifurcated MS and elongated LA are able to form in a bound orbit consistent with the HST proper motions. The LMC and SMC have remained bound to each other only recently in our model despite being separately bound to the Milky Way for more than 5 Gyr. We find that the MS and LA are able to form as a consequence of LMC-dominated tidal stripping during the recent dynamical coupling of the LMC and SMC. Our orbital model depends on our assumption that the Milky Way has a constant circular velocity of V_cir=250 km/s up to 160 kpc, which implies a massive isothermal halo that is not completely rejected by observations.

preprint2011arXiv

Formation of massive globular clusters with heavy element abundance spread in the Galactic building blocks

A growing number of recent observations have revealed that the Galactic globular cluster (GC) omega Cen is not the only GC that shows abundance spread in heavy elements (e.g., Fe). In order to understand the origin of the Galactic GCs with heavy element abundance spread ("HEAS"), we investigate the formation processes of massive GCs (MGCs) with masses larger than 10^6 M_sun in gas-rich dwarf galaxies interacting and merging with the very young Galaxy. We find that massive and compact stellar clumps with masses larger than 10^6 M_sun, which can be regarded as progenitors of MGCs, can form from massive gas clumps that are developed through merging of gaseous regions initially at different regions and thus with different metallicities. Therefore it is inevitable that MGCs formed in dwarfs have HEAS. The abundance spread in each individual MGC depends on the radial metallicity gradient of the host dwarf such that it can be larger for the steeper metallicity gradient. For example, MGCs formed in a dwarf with a central metallicity of [Fe/H]=-1.1 and the radial gradient of -0.2 dex kpc^{-1} can have the abundance spread of Delta [Fe/H] ~.2. The simulated MGCs appear to be significantly flattened owing to their dissipative formation from gas disks of their host dwarfs. Based on these results, we discuss possibly diverse formation mechanisms for the Galactic GCs such as M22, M54, NGC 2419, omega Cen, and Terzan 5.

preprint2011arXiv

Formation of the Galactic bulge from a two-component stellar disk: Explaining cylindrical rotation and vertical metallicity gradient

Recent observational studies have revealed that the Galactic bulge has cylindrical rotation and a steeper vertical metallicity gradient. We adopt two representative models for the bulge formation and thereby investigate whether the two models can explain both the observed cylindrical rotation and vertical metallicity gradient in a self-consistent manner. One is the "pure disk scenario" (PDS) in which the bulge is formed from a pure thin stellar disk through spontaneous bar instability. The other is the "two-component disk scenario" (TCDS) in which the bulge is formed from a disk composed of thin and thick disks through bar instability. Our numerical simulations show that although PDS can reproduce the cylindrical rotation, it shows a rather flatter vertical metallicity gradient that is inconsistent with observations. The derived flatter metallicity gradient is due to the vertical mixing of stars with different initial metallicities by the stellar bar. This result implies that the bulge can not be simply formed from a pure thin stellar disk. On the other hand, the bulge formed from the two-component disk in TCDS can explain both the observed cylindrical rotation and vertical metallicity gradient of the Galactic bulge reasonably well. In TCDS, more metal-poor stars at higher |z| (vertical distance) which originate from the already dynamically hotter thick disk can not be strongly influenced by vertical mixing of the bar so that they can stay in situ for longer timescales and thus keep the lower metallicity at higher |z|. Consequently, the vertical metallicity gradient of the bulge composed of initially thin and thick disk stars can not be so flattened, even if the gradient of the thin disk can be flattened significantly by the bar in TCDS.

preprint2011arXiv

Global Properties of `Ordinary' Early-type Galaxies: photometry and spectroscopy of stars and globular clusters in NGC 4494

We present a comprehensive analysis of the spatial, kinematic, and chemical properties of stars and globular clusters (GCs) in the `ordinary' elliptical galaxy NGC 4494 using data from the Keck and Subaru telescopes. We derive galaxy surface brightness and colour profiles out to large galactocentric radii. We compare the latter to metallicities derived using the near-infrared Calcium Triplet. We obtain stellar kinematics out to ~3.5 effective radii. The latter appear flattened or elongated beyond ~1.8 effective radii in contrast to the relatively round photometric isophotes. In fact, NGC 4494 may be a flattened galaxy, possibly even an S0, seen at an inclination of ~45 degrees. We publish a catalogue of 431 GC candidates brighter than i=24 based on the photometry, of which 109 are confirmed spectroscopically and 54 have measured spectroscopic metallicities. We also report the discovery of 3 spectroscopically confirmed ultra-compact dwarfs around NGC 4494 with measured metallicities of -0.4<[Fe/H]<-0.3. Based on their properties, we conclude that they are simply bright GCs. The metal-poor globular clusters are found to be rotating with similar amplitude as the galaxy stars, while the metal-rich globular clusters show marginal rotation. We supplement our analysis with available literature data and results. Using model predictions of galaxy formation, and a suite of merger simulations, we find that many of the observational properties of NGC 4494 may be explained by formation in a relatively recent gas-rich major merger. Complete studies of individual galaxies incorporating a range of observational avenues and methods such as the one presented here will be an invaluable tool for constraining the fine details of galaxy formation models, especially at large galactocentric radii.

preprint2011arXiv

On interpretation of recent proper motion data for the Large Magellanic Cloud

Recent observational studies using the Hubble Space Telescope (HST) have derived the center-of-mass proper motion (CMPM) of the Large Magellanic Cloud (LMC). Although these studies carefully treated both rotation and perspective effects in deriving the proper motion for each of the sampled fields, they did not consider the effects of local random motion in the derivation. This means that the average PM of the fields (i.e., the observed CMPM) could significantly deviate from the true CMPM, because the effect of local random motion can not be close to zero in making the average PM for the small number of the fields (~10). We discuss how significantly the observationally derived CMPM can deviate from the true CMPM by applying the same method as used in the observations for a dynamical model of the LMC with a known true CMPM. We find that the deviation can be as large as ~ 50 km/s (~0.21 mas/yr), if the LMC has a thick disk and a maximum circular velocity of ~120$ km/s. We also find that the deviation depends both on the total number of the sampled fields and on structure and kinematics of the LMC. We therefore suggest that there is a possibility that the observed CMPM of the LMC deviates from the true one to some extent. We also show that a simple mean of PM for a large number of the LMC fields (~1000) can be much closer to the true CMPM.

preprint2011arXiv

On the origin of the stellar halo and multiple stellar populations in the globular cluster NGC 1851

We propose that the observed stellar halo around the globular cluster (GC) NGC 1851 is evidence for its formation in the central region of its defunct host dwarf galaxy. We numerically investigate the long-term dynamical evolution of a nucleated dwarf galaxy embedded in a massive dark matter halo under the strong tidal field of the Galaxy. The dwarf galaxy is assumed to have a stellar nucleus (or a nuclear star cluster) that could be the progenitor for NGC 1851. We find that although the dark matter halo and the stellar envelope of the host dwarf of NGC 1851 can be almost completely stripped during its orbital evolution around the Galaxy, a minor fraction of stars in the dwarf can remain trapped by the gravitational field of the nucleus. The stripped nucleus can be observed as NGC 1851 with no/little dark matter whereas stars around the nucleus can be observed as a diffuse stellar halo around NGC 1851. The simulated stellar halo has a symmetric distribution with a power-law density slope of ~ -2 and shows no tidal tails within ~200pc from NGC 1851. We show that two GCs can merge with each other to form a new nuclear GC embedded in field stars owing to the low stellar velocity dispersion of the host dwarf. This result makes no assumption on the ages and/or chemical abundances of the two merging GCs. Thus the observed stellar halo and characteristic multiple stellar populations in NGC 1851 suggest that NGC 1851 could have formed initially in the central region of an ancient dwarf galaxy.

preprint2011arXiv

Origin of chemical and dynamical properties of the Galactic thick disk

We adopt a scenario in which the Galactic thick disk was formed by minor merging between the first generation of the Galactic thin disk (FGTD) and a dwarf galaxy about 9 Gyr ago and thereby investigate chemical and dynamical properties of the Galactic thick disk. In this scenario, the dynamical properties of the thick disk have long been influenced both by the mass growth of the second generation of the Galactic thin disk (i.e., the present thin disk) and by its non-axisymmetric structures. On the other hand, the early star formation history and chemical evolution of the thin disk was influenced by the remaining gas of the thick disk. Based on N-body simulations and chemical evolution models, we investigate the radial metallicity gradient, structural and kinematical properties, and detailed chemical abundance patterns of the thick disk. Our numerical simulations show that the ancient minor merger event can significantly flatten the original radial metallicity gradient of the FGTD, in particular, in the outer part, and also can be responsible for migration of inner metal-rich stars into the outer part (R>10kpc). The simulations show that the central region of the thick disk can develop a bar due to dynamical effects of a separate bar in the thin disk. The simulated orbital eccentricity distributions in the thick disk for models with higher mass-ratios (~0.2) and lower orbital eccentricities (~ 0.5) of minor mergers are in good agreement with the corresponding observations. The simulated V_{phi}-|z| relation of the thick disk in models with low orbital inclination angles of mergers are also in good agreement with the latest observational results. Our Galactic chemical evolution models can explain both the observed metallicity distribution functions (MDFs) and correlations between [Mg/Fe] and [Fe/H] for the two disks in a self-consistent manner.

preprint2011arXiv

Stellar population gradients in the cores of nearby field E+A galaxies

We have selected a sample of local E+A galaxies from the Sloan Digital Sky Survey (SDSS) Data Release 7 for follow up integral field spectroscopy with the Wide Field Spectrograph (WiFeS) on the ANU 2.3-m telescope. The sample was selected using the Halpha line in place of the [OII]3727 line as the indicator of on-going star formation (or lack thereof). This allowed us to select a lower redshift sample of galaxies than available in the literature since the [OII]3727 falls off the blue end of the wavelength coverage in the SDSS for the very lowest redshift objects. This low redshift selection means that the galaxies have a large angular to physical scale which allows us to resolve the central ~1kpc region of the galaxies; the region where stellar population gradients are expected. Such observations have been difficult to make using other higher redshift samples because even at redshifts z~0.1 the angular to physical scale is similar to the resolution provided by ground based seeing. Our integral field spectroscopy has enabled us to make the first robust detections of Balmer line gradients in the centres of E+A galaxies. Six out of our sample of seven, and all the galaxies with regular morphologies, are observed to have compact and centrally-concentrated Balmer line absorption. This is evidence for compact young cores and stellar population gradients which are predicted from models of mergers and tidal interactions which funnel gas into the galaxy core. Given the generally isolated nature of our sample this argues for the galaxies being seen in the late stage of a merger where the progenitors have already coalesced.

preprint2011arXiv

The Effect of Drag from the Galactic Hot Halo on the Magellanic Stream and Leading Arm

We study the effect of drag induced by the Galactic hot halo on the two neutral hydrogen (HI) cloud complexes associated with the Large and Small Magellanic Clouds: the Magellanic Stream (MS) and the Leading Arm (LA). In particular, we adopt the numerical models of previous studies and re-simulate the tidal formation of the MS and LA with the inclusion of a drag term. We find that the drag has three effects which, although model-dependent, may bring the tidal formation scenario into better agreement with observations: correcting the LA kinematics, reproducing the MS column density gradient, and enhancing the formation of MS bifurcation. We furthermore propose a two-stage mechanism by which the bifurcation forms. In general, the inclusion of drag has a variety of both positive and negative effects on the global properties of the MS and LA, including their on-sky positions, kinematics, radial distances, and column densities. We also provide an argument which suggests that ram pressure stripping and tidal stripping are mutually exclusive candidates for the formation of the MS and LA.

preprint2011arXiv

Transformation from spirals into S0s with bulge growth in groups of galaxies

Recent observations have revealed that the time evolution of the S0 number fraction at intermediate and high redshifts (0.2<z<0.8) is more dramatic in groups of galaxies than in clusters. In order to understand the origin of S0s in groups, we investigate numerically the morphological transformation of spirals into S0s through group-related physical processes. Our chemodynamical simulations show that spirals in group environments can be strongly influenced by repetitive slow encounters with group member galaxies so that those with thin disks and prominent spiral arm structures can be transformed into S0s with thick disks and without prominent spiral arm structure. Such tidal interactions can also trigger repetitive starbursts within the bulges of spirals and consequently increase significantly the masses of their bulges. Owing to rapid consumption of gas initially in spirals during the bulge growth, the S0s can become gas-poor. The S0s transformed from spirals in this way have young and metal-rich stellar populations in the inner regions of their bulges. The simulated S0s have lower maximum rotational velocities and flatter radial line-of-sight velocity dispersion profiles in comparison to their progenitor spirals. The formation processes of S0s due to tidal interactions depend not only on the masses and orbits of the progenitor spirals, but also on group mass. A significant fraction (10-30%) of stars and gas can be stripped during this spiral to S0 morphological transformation so that intragroup stars and gas can be formed. Based on these results, we discuss structures, kinematics, chemical properties, and the Tully-Fisher relation of S0s in groups.

preprint2011arXiv

When was the Large Magellanic Cloud accreted onto the Galaxy ?

Using fully self-consistent N-body models for the dynamical evolution of the Large Magellanic Cloud (LMC) in the Galaxy, we show that if the LMC initially has an extended old stellar halo before its commencement of tidal interaction with the Galaxy, physical properties of the stars stripped from the LMC stellar halo can have fossil information as to when and where the LMC was accreted onto the Galaxy for the first time. If the epoch of the first LMC accretion onto the Galaxy from outside its viriral radius is more than ~4 Gyr ago (i.e., at least two pericenter passages), the stars stripped from the stellar halo of the LMC can form an irregular polar ring or a thick disk with a size of ~100 kpc and rotational kinematics. On the other hand, if the LMC was first accreted onto the Galaxy quite recently (~ 2 Gyr ago), the stripped stars form shorter leading and trailing stellar stream at R=50-120 kpc. Also distributions of the stripped stars in phase space between the two cases can be significantly different. The derived differences in structure and kinematics of the stripped stars therefore suggest that if we compare the observed three-dimensional (3D) distribution and kinematics of the outer Galactic stellar halo along the polar-axis, then we can give strong constraints on the past orbit of the LMC. We find that the orbital properties of the LMC in the successful formation models of the Magellanic Stream (MS) are consistent with those predicted from recent LambdaCDM cosmological simulations. We conclude that the LMC was accreted onto the Galaxy more than ~4 Gyr ago so that interaction between the LMC, the Small Magellanic Cloud (SMC), and the Galaxy could form the MS and its Leading Arms (LAs).

preprint2010arXiv

On the transition from nuclear-cluster to black-hole dominated galaxy cores

Giant elliptical galaxies, believed to be built from the merger of lesser galaxies, are known to house a massive black hole at their center rather than a compact star cluster. If low- and intermediate-mass galaxies do indeed partake in the hierarchical merger scenario, then one needs to explain why their dense nuclear star clusters are not preserved in merger events. A valuable clue may the recent revelation that nuclear star clusters and massive black holes frequently co-exist in intermediate mass bulges and elliptical galaxies. In an effort to understand the physical mechanism responsible for the disappearance of nuclear star clusters, we have numerically investigated the evolution of merging star clusters with seed black holes. Using black holes that are 1-5% of their host nuclear cluster mass, we reveal how their binary coalescence during a merger dynamically heats the newly wed star cluster, expanding it, significantly lowering its central stellar density, and thus making it susceptible to tidal destruction during galaxy merging. Moreover, this mechanism provides a pathway to explain the observed reduction in the nucleus-to-galaxy stellar mass ratio as one proceeds from dwarf to giant elliptical galaxies.

preprint2010arXiv

Origin of the unusually low nitrogen abundances in young populations of the Large Magellanic Cloud

It is a longstanding problem that HII regions and very young stellar populations in the Large Magellanic Cloud (LMC) have the nitrogen abundances ([N/H]) by a factor of ~7 lower than the solar value. We here discuss a new scenario in which the observed unusually low nitrogen abundances can be closely associated with recent collision and subsequent accretion of HI high velocity clouds (HVCs) that surround the Galaxy and have low nitrogen abundances. We show that if the observed low [N/H] is limited to very young stars with ages less than ~10^7 yr, then the collision/accretion rate of the HVCs onto the LMC needs to be ~ 0.2 M_sun/yr (corresponding to the total HVC mass of 10^6-10^7 M_sun) to dilute the original interstellar medium (ISM) before star formation. The required accretion rate means that even if the typical mass of HVCs accreted onto the LMC is ~ 10^7 M_sun, the Galaxy needs to have ~2500 massive HVCs within the LMC's orbital radius with respect to the Galactic center. The required rather large number of massive HVCs drives us to suggest that the HVCs are not likely to efficiently dilute the ISM of the LMC and consequently lower the [N/H]. We thus suggest the transfer of gas with low [N/H] from the Small Magellanic Cloud (SMC) to the LMC as a promising scenario that can explain the observed low [N/H].

preprint2010arXiv

Rotation and multiple stellar population in globular clusters

We investigate structure and kinematics of the second generation of stars (SG) formed from gaseous ejecta of the first generation of stars (FG) in forming globular clusters (GCs). We consider that SG can be formed from gaseous ejecta from AGB stars of FG with the initial total mass of 10^6-10^8 M_sun to explain the present masses of the Galactic GCs. Our 3D hydrodynamical simulations with star formation show that SG formed in the central regions of FG can have a significant amount of rotation (V/sigma ~0.8-2.5). The rotational amplitude of SG can depend strongly on the initial kinematics of FG. We thus propose that some GCs composed of FG and SG had a significant amount of rotation when they were formed. We also suggest that although later long-term (~10 Gyr) dynamical evolution of stars can smooth out the initial structural and kinematical differences between FG and SG to a large extent, initial flattened structures and rotational kinematics of SG can be imprinted on shapes and internal rotation of the present GCs. We discuss these results in terms of internal rotation observed in the Galactic GCs.

preprint2010arXiv

Secondary star formation within massive star clusters: Origin of multiple stellar populations in globular clusters

We numerically investigate whether and how gaseous ejecta from AGB stars can be converted into new stars within originally massive star clusters (MSCs) in order to understand the origin of multiple stellar populations in globular clusters (GCs). We adopt a scenario in which (i) MSCs with masses of M_s can be formed from high-mass, high-density giant molecular clouds (GMCs) in their host galactic building blocks embedded in dark matter halos at high redshifts and (ii) their evolution therefore can be significantly influenced by M_s, their initial locations, and physical properties of their hosts. Our 3D hydrodynamical simulations show that gaseous ejecta from AGB stars can be retained within MSCs and consequently converted into new stars very efficiently in the central regions of MSCs, only if M_s exceed a threshold mass (M_th) of ~10^6 M_sun. The new stars can correspond to the ``second generation (SG)'' of stars with higher Na and lower O abundances observed in GCs. Star formation efficiencies during the formation of SG stars within MSCs with M_s > M_th can be rather high (0.3-0.9) so that very compact new clusters within original MSCs can be formed. M_s should be as large as 10^6-10^7 M_sun to explain the observed large fraction of SG stars in the present ordinary Galactic GCs, because new stars can consist of only 1-4% among all stars for the standard IMF. Nuclear MSCs are found to retain much more effectively the AGB ejecta and convert more efficiently the gas into new stars owing to much deeper gravitational potential of their hosts. We suggest that both M_s and their locations within their hosts can determine whether abundance spread can be seen only in light elements or even in heavy ones.

preprint2010arXiv

Synchronized formation of starburst and poststarburst galaxies in merging clusters of galaxies

We propose that synchronized triggering of star formation in gas-rich galaxies is possible during major mergers of cluster of galaxies, based on new numerical simulations of the time evolution of the physical properties of the intracluster medium (ICM) during such a merger event. Our numerical simulations show that the external pressure of the ICM in which cluster member galaxies are embedded, can increase significantly during cluster merging. As such, efficient star formation can be triggered in gas-rich members as a result of the strong compression of their cold gas by the increased pressure. We also suggest that these star-forming galaxies can subsequently be transformed into poststarburst galaxies, with their spatial distribution within the cluster being different to the rest of its population. We discuss whether this possible merger-induced enhancement in the number of star-forming and post-star-forming cluster galaxies is consistent with the observed evolution of galaxies in merging clusters.

preprint2010arXiv

The Star Cluster Population of the Collisional Ring Galaxy NGC 922

We present a detailed study of the star cluster population detected in the galaxy NGC922, one of the closest collisional ring galaxies known to date, using HST/WFPC2 UBVI photometry, population synthesis models, and N-body/SPH simulations.We find that 69% of the clusters are younger than 7Myr, and that most of them are located in the ring or along the bar, consistent with the strong Halpha emission. The cluster luminosity function slope of 2.1-2.3 for NGC922 is in agreement with those of young clusters in nearby galaxies. Models of the cluster age distribution match the observations best when cluster disruption is considered. We also find clusters with ages (>50Myr) and masses (>10^5 Msun) that are excellent progenitors for faint fuzzy clusters. The images also show a tidal plume pointing toward the companion. Its stellar age from our analysis is consistent with pre-existing stars that were stripped off during the passage of the companion. Finally, a comparison of the star-forming complexes observed in NGC922 with those of a distant ring galaxy from the GOODS field indicates very similar masses and sizes, suggesting similar origins.

preprint2009arXiv

Dark impact and galactic star formation: Origin of the Gould belt

The Milky Way has a giant stellar structure in the solar neighborhood, which has a size of $\sim 1$ kpc, a mass of $\sim 10^6 {\rm M}_{\odot}$, and a ring-like distribution of young stars. Fundamental physical properties of this local enigmatic structure, known as the Gould belt (GB), have not been reproduced by previously proposed models. We first show that the local enigmatic structure can be formed about 30 Myr ago as a result of a high-speed, oblique collision between a gas cloud with a mass of $\sim 10^6 {\rm M}_{\odot}$ and a dark matter clump with a mass of $\sim 10^7 {\rm M}_{\odot}$ based on numerical simulations of the collision. We find that strong dynamical impact of the clump transforms the flattened cloud into a ring-like stellar structure after induced star formation within the cloud. Our simulations furthermore demonstrate that the stellar structure is moderately elongated and significantly inclined with respect to the disk of the Milky Way owing to the strong tidal torque by the colliding clump. We thus suggest that the GB is one of stellar substructures formed from collisions between gas clouds and dark matter clumps predicted in the hierarchical clustering scenario of galaxy formation. We also suggest that collisions of dark matter clumps with their host galaxies can significantly change star formation histories for some of their gas clouds thus influence galactic global star formation histories to some extent.

preprint2009arXiv

Origin of rotational kinematics in the globular cluster system of M31: A new clue to the bulge formation

We propose that the rotational kinematics of the globular cluster system (GCS) in M31 can result from a past major merger event that could have formed its bulge component. We numerically investigate kinematical properties of globular clusters (GCs) in remnants of galaxy mergers between two disks with GCs in both their disk and halo components. We find that the GCS formed during major merging can show strongly rotational kinematics with the maximum rotational velocities of 140 - 170 km/s for a certain range of orbital parameters of merging. We also find that a rotating stellar bar, which can be morphologically identified as a boxy bulge if seen edge-on, can be formed in models for which the GCSs show strongly rotational kinematics. We thus suggest that the observed rotational kinematics of GCs with different metallicities in M31 can be closely associated with the ancient major merger event. We discuss whether the formation of the rotating bulge/bar in M31 can be due to the ancient merger.

preprint2009arXiv

Ram pressure stripping of halo gas in disk galaxies: Implications on galactic star formation in different environments

We numerically investigate evolution of gaseous halos around disk galaxies in different environments ranging from small groups to rich clusters in order to understand galaxy evolution in these environments. Our simulations self-consistently incorporate effects of ram pressure of intergalactic medium (IGM) on disk and halo gas of galaxies and hydrodynamical interaction between disk and halo gas so that mass fractions of halos gas stripped by ram pressure of IGM (F_strip) can be better estimated. We mainly investigate how F_strip depends on total masses of their host environments (M_host}), galactic masses (M_gal), densities and temperature of IGM (T_IGM and rho_IGM, respectively), relative velocities between IGM and galaxies (V_r), and physical properties of disks (e.g., gas mass fraction). We find that typically 60-80% of halo gas can be efficiently stripped from Milky Way-type disk galaxies by ram pressure in clusters with M_host 10^{14} M_sun We also find that F_strip depends on M_host such that F_strip is higher for larger M_host. Furthermore it is found that F_strip can be higher in disk galaxies with smaller M_gal for a given environment. Our simulations demonstrate that the presence of disk gas can suppress ram pressure stripping of halo gas owing to hydrodynamical interaction between halo and disk gas. Ram pressure stripping of halo gas is found to be efficient (i.e., F_strip>0.5) even in small and/or compact groups, if rho_IGM ~ 10^5 M_sun/kpc^3 and V_r ~ 400 km/s. Based on the derived radial distributions of remaining halo gas after ram pressure stripping, we propose that truncation of star formation after halo gas stripping can occur outside-in in disk galaxies.

preprint2009arXiv

Young Stellar Populations in the Collisional Ring Galaxy NGC 922

We studied the star cluster population properties in the nearby collisional ring galaxy NGC 922 using HST/WFPC2 photometry and population synthesis modeling. We found that 69% of the detected clusters are younger than 7 Myr, and that most of them are located in the ring or along the bar, consistent with the strong Halpha emission. The images also show a tidal plume pointing toward the companion. Its stellar age is consistent with pre-existing stars that were probably stripped off during the passage of the companion. We compared the star-forming complexes observed in NGC 922 with those of a distant ring galaxy from the GOODS field. It indicates very similar masses and sizes, suggesting similar origins. Finally, we found clusters that are excellent progenitor candidates for faint fuzzy clusters.

preprint2007arXiv

Evolution of the Small Magellanic Cloud

Based on the results of N-body simulations on the last 2.5 Gyr evolution of the Large and Small Magellanic Clouds (LMC and SMC, respectively) interacting with the Galaxy, we firstly show when and where the leading arms (LAs) of the Magellanic stream (MS) can pass through the Galactic plane after the MS formation. We secondly show collisions between the outer Galactic HI disk and the LAs of the MS can create giant HI holes and chimney-like structures in the disk about 0.2 Gyr ago. We thirdly show that a large amount of metal-poor gas is stripped from the SMC and transfered to the LMC during the tidal interaction between the Clouds and the Galaxy about 0.2 and 1.3 Gyr ago. We thus propose that this metal-poor gas can closely be associated with the origin of LMC's young and intermediate-age stars and star clusters with distinctively low-metallicities with [Fe/H] < -0.6.

preprint2006arXiv

The U-shaped distribution of globular cluster specific frequencies in a biased globular cluster formation scenario

Using high-resolution numerical simulations, we investigate mass- and luminosity-normalized specific frequencies (T_N and S_N, respectively) of globular cluster systems (GCSs) in order to understand the origin of the observed U-shaped relation between S_N and V-band magnitude (M_V) of their host galaxies. We adopt a biased GC formation scenario in which GC formation is truncated in galaxy halos that are virialized at a later redshift, z_trun. T_N is derived for galaxies with GCs today and converted into S_N for reasonable galaxy mass-to-light-ratios (M/L). We find that T_N depends on halo mass (M_h) in the sense that T_N can be larger in more massive halos with M_h > 10^9 M_sun, if z_trun is as high as 15. We however find that the dependence is too weak to explain the observed S_N-M_V relation and the wide range of S_N in low-mass early-type galaxies with -20.5 < M_V < -16.0 mag for a reasonable constant M/L. The M_V-dependence of S_N for the low-mass galaxies can be well reproduced, if the mass-to-light-ratio M_h/L_V \propto M_h^α, where αis as steep as -1. Based on these results, we propose that the origin of the observed U-shaped S_N-M_V relation of GCSs can be understood in terms of the bimodality in the dependence of M_h/L_V on M_h of their host galaxies. We also suggest that the observed large dispersionin S_N in low-mass galaxies is due partly to the large dispersion in T_N.

preprint2006arXiv

Triggered star formation in the Magellanic Clouds

We discuss how tidal interaction between the Large Magellanic Cloud (LMC), the Small Magellanic Cloud (SMC), and the Galaxy triggers galaxy-wide star formation in the Clouds for the last ~0.2 Gyr based on our chemodynamical simulations on the Clouds. Our simulations demonstrate that the tidal interaction induces the formation of asymmetric spiral arms with high gas densities and consequently triggers star formation within the arms in the LMC. Star formation rate in the present LMC is significantly enhanced just above the eastern edge of the LMC's stellar bar owing to the tidal interaction. The location of the enhanced star formation is very similar to the observed location of 30 Doradus, which suggests that the formation of 30 Doradus is closely associated with the last Magellanic collision about 0.2 Gyr ago. The tidal interaction can dramatically compress gas initially within the outer part of the SMC so that new stars can be formed from the gas to become intergalactic young stars in the inter-Cloud region (e.g., the Magellanic Bridge). The metallicity distribution function of the newly formed stars in the Magellanic Bridge has a peak of [Fe/H] ~ -0.8, which is significantly lower than the stellar metallicity of the SMC.

preprint2002arXiv

Extended HI spiral structure and the figure rotation of triaxial dark halos

The HI disk of the blue compact dwarf (BCD) galaxy NGC 2915 extends to 22 optical scalelengths and shows spiral arms reaching far beyond the optical component. None of the previous theories for spiral structure provide likely explanations for these very extended spiral arms. Our numerical simulations first demonstrate that such large spiral arms can form in an extended gas disk embedded in a massive triaxial dark matter halo with slow figure rotation, through the strong gravitational torque of the rotating halo. We then show that the detailed morphological properties of the developed spirals and rings depend strongly on the pattern speed of the figure rotation, the shape of the triaxial halo, and the inclination of the disk with respect to the plane including the triaxial halo's long and middle axes. These results strongly suggest that the dark matter halo of NGC 2915 is triaxial and has figure rotation. Based on these results, we also suggest that dynamical effects of triaxial halos with figure rotation are important in various aspect of galaxy formation and evolution, such as formation of polar ring galaxies, excitation of non-axisymmetric structures in low surface-brightness galaxies, and gas fueling to the central starburst regions of BCDs.

preprint2002arXiv

Galactic Recycling: The HI Ring Around NGC 1533

We report the discovery of a new HI ring around the S0 galaxy NGC 1533. The ring orbits at a radius of 35 kpc, well outside the optical extent of the galaxy. We have conducted N-body/SPH numerical simulations to show this HI ring could be the merger remnant of a tidally destroyed galaxy. We find no optical component associated with the HI ring. However, observations hint at H-alpha emission associated with the SE part of the ring only. The H-alpha is in the form of a few very small isolated emission line regions. The large HI velocity dispersions (up to 30 km/s) and velocity gradients (up to 50 km/s/kpc) in this region indicate the H-alpha emission could be due to star formation triggered by clouds colliding within the ring.

preprint2001arXiv

Tidal truncation of gas replenishment and global suppression of galactic star formation in distant clusters

Recent spectroscopic observations of galaxies in distant clusters have revealed that the rate of star formation in star-forming galaxies is significantly suppressed with respect to their counterparts in the field at a similar redshift. It is, however, highly uncertain which physical processes are responsible for this suppression. We present the results of a numerical investigation of how the global tidal field of a cluster dynamically influences the reservoir of halo gas surrounding a disk galaxy as it falls into the cluster from the surrounding field. We find that the tidal field of the cluster efficiently removes the halo gas from the galaxy, thereby halting its accretion onto the disk, and thus the fueling of star formation within. This effectively truncates the galaxy's star formation. We also find that this tidal truncation does not depend very strongly on the orbit of the disk with respect to the center of the cluster. These results suggest that the global tidal field of clusters is capable of causing a widespread and uniform suppression of star formation in galaxies accreted by the cluster. In light of these results, we discuss the importance of this tidal truncation of the gas supply in the formation of S0 galaxies in clusters.

preprint2000arXiv

ALMA and high redshift dusty starburst mergers

By using a new numerical code for deriving the spectral energy distributions of galaxies, we have investigated the time evolution of morphological properties, the star formation rate, and the submillimeter flux at 850 $μ$m in high-redshift ($z$) dusty starburst mergers with mass ratio ($m_{2}$) of two disks ranging from 0.1 (minor merger) to 1.0 (major one). We found that the maximum star-formation rate, the degree of dust extinction, and the 850 $μ$m flux are larger for mergers with larger $m_{2}$. The 850 $μ$m flux from mergers at 1.5 $\le z \le$ 3.0 in the observer frame is found to be a few mJy for major merger cases, and at most $\sim$ 100 $μ$Jy for minor ones. This result suggests that only high-redshift major mergers are now detected by SCUBA with the current 850 $μm$ detection limit of a few mJy. These results imply that LMSA (ALMA) with the expected detection limit of the order of 10 $μ$Jy at 850 $μm$ can be used to study high-redshift mergers with variously different $m_2$, and thus provide an important clue to the formation of galaxies in a high-redshift universe.

preprint1999arXiv

Morphological and photometric evolution of ultra-luminous infrared galaxies: Nature of faint SCUBA sources

We investigate when and how a dusty starburst galaxy merger can be heavily obscured by dust and consequently becomes an ultra-luminous infrared galaxy (ULIRG), based on numerical simulations of chemodynamical and photometric evolution of dusty gas-rich major galaxy mergers. We found that a major galaxy merger is more likely to become an ULIRG preferentially in the merger late phase, when the two disks become very close and the very high-density dusty gas can obscure heavily the central secondary massive starburst. We furthermore show how the optical and near-infrared morphology of a simulated ULIRG at intermediate (z=0.4) and high redshift (z=1 - 2) can be observed by the Hubble Space Telescope (HST) in order to present a plausible explanation for the origin of some host galaxies of the faint SCUBA sources recently observed by Smail et al. (1998). The results of our numerical simulations imply that some SCUBA sources with apparently faint and compact HST morphology can be higher redshift dust-enshrouded starburst mergers with their outer low surface brightness tidal features hardly detectable in the present optics of the HST.

preprint1998arXiv

Group-cluster merging and the formation of starburst galaxies

A significant fraction of clusters of galaxies are observed to have substructure, which implies that merging between clusters and subclusters is a rather common physical process of cluster formation. It still remains unclear how cluster merging affects the evolution of cluster member galaxies. We report the results of numerical simulations, which show the dynamical evolution of a gas-rich late-type spiral in a merger between a small group of galaxies and a cluster. The simulations demonstrate that time-dependent tidal gravitational field of the merging excites non-axisymmetric structure of the galaxy, subsequently drives efficient transfer of gas to the central region, and finally triggers a secondary starburst. This result provides not only a new mechanism of starbursts but also a close physical relationship between the emergence of starburst galaxies and the formation of substructure in clusters. We accordingly interpret post-starburst galaxies located near substructure of the Coma cluster as one observational example indicating the global tidal effects of group-cluster merging. Our numerical results furthermore suggest a causal link between the observed excess of blue galaxies in distant clusters and cluster virialization process through hierarchical merging of subclusters.

preprint1998arXiv

Stellar populations in gas-rich galaxy mergers II. Feedback effects of Type Ia and II supernovae

We numerically investigate chemodynamical evolution of major disk-disk galaxy mergers in order to explore the origin of mass-dependent chemical, photometric, and spectroscopic properties observed in elliptical galaxies. We particularly investigate the dependence of the fundamental properties on merger progenitor disk mass (M_d). Main results obtained in this study are the following three. (1) More massive (luminous) ellipticals formed by galaxy mergers between more massive spirals have larger metallicity (Z) and thus show redder colors: The typical metallicity ranges from ~ 1.0 solar abundance (Z~ 0.02) for ellipticals formed by mergers with M_d = 10^10 M_solar to ~ 2.0 solar (Z ~ 0.04) for those with M_d = 10^12 M_solar. (2) The absolute magnitude of negative metallicity gradients developed in galaxy mergers is more likely to be larger for massive ellipticals. Absolute magnitude of metallicity gradient correlates with that of age gradient in ellipticals in the sense that an elliptical with steeper negative metallicity gradient is more likely to show steeper age gradient. (3) Both M/L_B and M/L_K, where M, L_B, and L_K are total stellar mass of galaxy mergers, B-band luminosity, and K-band one, respectively, depend on galactic mass in such a way that more massive ellipticals have larger M/L_B and smaller M/L_K.

preprint1998arXiv

Unequal-mass galaxy mergers and the creation of cluster S0 galaxies

It is a longstanding and remarkable problem when and how red S0 galaxies were formed in clusters of galaxies. We here propose that the major mechanism for the S0 creation is galaxy merging between two spirals with unequal mass. Our numerical simulations demonstrate that galaxy merging exhausts a large amount of interstellar medium of two gas-rich spirals owing to the moderately enhanced star formation, and subsequently transforms the two into one gas-poor S0 galaxy with structure and kinematics strikingly similar to the observed ones. This secondary S0 formation with enhanced star formation explains a smaller fraction of S0 population recently observed in some distant clusters of galaxies. Unequal-mass galaxy mergers thus provide an evolutionary link between a larger number of blue spirals observed in intermediate redshift clusters and red S0s prevalent in the present-day ones.