Researcher profile

Kenji Bekki

Kenji Bekki contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
17works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

17 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 &#39;AGB with dilution&#39; 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&#39;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&#39;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&#39;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.