BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Daisuke Kawata

Daisuke Kawata appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.COastro-phastro-ph.SRastro-ph.EPastro-ph.IMComputation and Languagephysics.comp-ph

Catalog footprint

What is connected

52works
8topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

52 published item(s)

preprint2026arXiv

Towards Galactic Archaeology with Inferred Ages of Giant Stars From Gaia Spectra

In the era of Gaia, the accurate determination of stellar ages is transforming Galactic archaeology. We demonstrate the feasibility of inferring stellar ages from Gaia's RVS spectra and the BP/RP (XP) spectrophotometric data, specifically for red giant branch and high-mass red clump stars. We successfully train two machine learning models, dubbed SIDRA: Stellar age Inference Derived from Gaia spectRA to predict the age. The SIDRA-RVS model uses the RVS spectra and SIDRA-XP the stellar parameters obtained from the XP spectra. Both models use BINGO, an APOGEE-derived stellar age as the training data. SIDRA-RVS estimates ages of stars whose age is around $τ_\mathrm{BINGO}=10$~Gyr with a standard deviation of residuals of $\sim$ 0.12 dex in the unseen test dataset, while SIDRA-XP achieves higher precision with residuals $\sim$ 0.064 dex for stars around $τ_\mathrm{BINGO}=10$~ Gyr. Since SIDRA-XP outperforms SIDRA-RVS, we apply SIDRA-XP to analyse the ages for 2,218,154 stars. This allowed us to map the chronological and chemical properties of Galactic disc stars, reproducing the known distinct features such as the Gaia-Sausage-Enceladus merger and a potential gas-rich interaction event linked to the first infall of the Sagittarius dwarf galaxy. This study demonstrates that machine learning techniques applied to Gaia's spectra can provide valuable individual age information, particularly for giant stars, thereby enhancing our understanding of the Milky Way's formation and evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2025arXiv

Modelling the Milky Way's exoplanet population based on cosmological galaxy simulations

In this paper we aim to simulate realistic exoplanet populations across different regions of the MW by combining state-of-the-art cosmological simulations of our Galaxy with exoplanet formation models and observations. We model the exoplanet populations around single stars, using planet occurrence rates and multiplicity depending on stellar mass, metallicity, and planet type, and assign them physical parameters such as mass and orbital period. Focussing first on the solar vicinity, we find mostly metallicity-driven differences in the distributions of non-hosting and planet-hosting single stars. In our simulated solar neighbourhood, 52.5% of all planets are Earth-like (23% of them located in the Habitable Zone), 44% are super-Earths/Neptunes, and 3.5% are giant planets. A comparison with the census of Kepler exoplanets and candidates shows that, when taking into account the most relevant selection effects, we obtain a similar distribution of exoplanets compared to the observed population. However, we also detect significant differences in the exoplanet and host star distributions (e.g. more planets around F-type and red-giant stars compared to observations) that we attribute mostly to a too strong recent star formation and a too large disc scale height in the simulation, as well as to some caveats in our exoplanet population synthesis that will be addressed in future work. Extending our analysis to other regions of the simulated MW and to other simulated galaxies, we find that the relative percentages of planet types remain largely consistent as long as the simulated galaxy matches the morphology and mass of the MW. We have created a fast and flexible framework to produce exoplanet populations based on MW-like simulations that can easily be adapted to produce predictions for the yields of future exoplanet detection missions. (abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Age distribution of stars in boxy/peanut/X-shaped bulges formed without bar buckling

Some barred galaxies, including the Milky Way, host a boxy/peanut/X-shaped bulge (BPX-shaped bulge). Previous studies suggested that the BPX-shaped bulge can either be developed by bar buckling or by vertical inner Lindblad resonance (vILR) heating without buckling. In this paper, we study the observable consequence of an BPX-shaped bulge built up quickly after bar formation via vILR heating without buckling, using an N-body/hydrodynamics simulation of an isolated Milky Way-like galaxy. We found that the BPX-shaped bulge is dominated by stars born prior to bar formation. This is because the bar suppresses star formation, except for the nuclear stellar disc (NSD) region and its tips. The stars formed near the bar ends have higher Jacobi energy, and when these stars lose their angular momentum, their non-circular energy increases to conserve Jacobi energy. This prevents them from reaching the vILR to be heated to the BPX region. By contrast, the NSD forms after the bar formation. From this simulation and general considerations, we expect that the age distributions of the NSD and BPX-shaped bulge formed without bar buckling do not overlap each other. Then, the transition age between these components betrays the formation time of the bar, and is testable in future observations of the Milky Way and extra-galactic barred galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Constraining Ultra Light Dark Matter with the Galactic Nuclear Star Cluster

We use the Milky Way's nuclear star cluster (NSC) to test the existence of a dark matter 'soliton core', as predicted in ultra-light dark matter (ULDM) models. Since the soliton core size is proportional to mDM^{-1}, while the core density grows as mDM^{2}, the NSC (dominant stellar component within about 3 pc) is sensitive to a specific window in the dark matter particle mass, mDM. We apply a spherical isotropic Jeans model to fit the NSC line-of-sight velocity dispersion data, assuming priors on the precisely measured Milky Way's supermassive black hole (SMBH) mass and the well-measured NSC density profile. We find that the current observational data reject the existence of a soliton core for a single ULDM particle with mass in the range 10^{-20.4} < mDM < 10^{-18.5} eV, assuming that the soliton core structure is not affected by the Milky Way's SMBH. We test our methodology on mock data, confirming that we are sensitive to the same range in ULDM mass as for the real data. Dynamical modelling of a larger region of the Galactic centre, including the nuclear stellar disc, promises tighter constraints over a broader range of mDM. We will consider this in future work.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Mira variables in the Milky Way's nuclear stellar disc: discovery and classification

The properties of the Milky Way's nuclear stellar disc give crucial information on the epoch of bar formation. Mira variables are promising bright candidates to study the nuclear stellar disc, and through their period-age relation dissect its star formation history. We report on a sample of $1782$ Mira variable candidates across the central $3\times3\,\mathrm{deg}^2$ of the Galaxy using the multi-epoch infrared VISTA Variables in Via Lactea (VVV) survey. We describe the algorithms employed to select candidate variable stars and then model their light curves using periodogram and Gaussian process methods. By combining with WISE, 2MASS and other archival photometry, we model the multi-band light curves to refine the periods and inspect the amplitude variation between different photometric bands. The infrared brightness of the Mira variables means many are too bright and missed by VVV. However, our sample follows a well-defined selection function as expected from artificial star tests. The multi-band photometry is modelled using stellar models with circumstellar dust that characterise the mass loss rates. We demonstrate how $\gtrsim90$ per cent of our sample is consistent with O-rich chemistry. Comparison to period-luminosity relations demonstrates that the bulk of the short period stars are situated at the Galactic Centre distance. Many of the longer period variables are very dusty, falling significantly under the O-rich Magellanic Cloud and solar neighbourhood period-luminosity relations and exhibit high mass-loss rates of $\sim2.5\times10^{-5}M_\odot\,\mathrm{yr}^{-1}$. The period distribution appears consistent with the nuclear stellar disc forming $\gtrsim8\,\mathrm{Gyr}$ ago although it is not possible to disentangle the relative contributions of the nuclear stellar disc and the contaminating bulge.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Constraining stellar population parameters from narrow band photometric surveys using convolutional neural networks

Upcoming large-area narrow band photometric surveys, such as J-PAS, will enable us to observe a large number of galaxies simultaneously and efficiently. However, it will be challenging to analyse the spatially-resolved stellar populations of galaxies from such big data to investigate galaxy formation and evolutionary history. We have applied a convolutional neural network (CNN) technique, which is known to be computationally inexpensive once it is trained, to retrieve the metallicity and age from J-PAS-like narrow band images. The CNN was trained using mock J-PAS data created from the CALIFA IFU survey and the age and metallicity at each data point, which are derived using full spectral fitting to the CALIFA spectra. We demonstrate that our CNN model can consistently recover age and metallicity from each J-PAS-like spectral energy distribution. The radial gradients of the age and metallicity for galaxies are also recovered accurately, irrespective of their morphology. However, it is demonstrated that the diversity of the dataset used to train the neural networks has a dramatic effect on the recovery of galactic stellar population parameters. Hence, future applications of CNNs to constrain stellar populations will rely on the availability of quality spectroscopic data from samples covering a wide range of population parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Age Dating the Galactic Bar with the Nuclear Stellar Disc

From the decades of the theoretical studies, it is well known that the formation of the bar triggers the gas funnelling into the central sub-kpc region and leads to the formation of a kinematically cold nuclear stellar disc (NSD). We demonstrate that this mechanism can be used to identify the formation epoch of the Galactic bar, using an N-body/hydrodynamics simulation of an isolated Milky Way-like galaxy. As shown in many previous literature, our simulation shows that the bar formation triggers an intense star formation for ~1 Gyr in the central region, and forms a NSD. As a result, the oldest age limit of the NSD is relatively sharp, and the oldest population becomes similar to the age of the bar. Therefore, the age distribution of the NSD tells us the formation epoch of the bar. We discuss that a major challenge in measuring the age distribution of the NSD in the Milky Way is contamination from other non-negligible stellar components in the central region, such as a classical bulge component. We demonstrate that because the NSD is kinematically colder than the other stellar populations in the Galactic central region, the NSD population can be kinematically distinguished from the other stellar populations, if the 3D velocity of tracer stars are accurately measured. Hence, in addition to the line-of-sight velocities from spectroscopic surveys, the accurate measurements of the transverse velocities of stars are necessary, and hence the near-infrared space astrometry mission, JASMINE, would play a crutial role to identify the formation epoch of the Galactic bar. We also discuss that the accuracy of stellar age estimation is also crucial to measure the oldest limit of the NSD stellar population.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Age dissection of the Milky Way discs: red giants in the Kepler field

[Abridged] Ensemble studies of red-giant stars with exquisite asteroseismic, spectroscopic, and astrometric constraints offer a novel opportunity to recast and address long-standing questions concerning the evolution of stars and of the Galaxy. Here, we infer masses and ages for nearly 5400 giants with available Kepler light curves and APOGEE spectra, and discuss some of the systematics that may affect the accuracy of the inferred stellar properties. First, we look at age-chemical-abundances relations. We find a dearth of young, metal-rich stars, and the existence of a significant population of old (8-9 Gyr), low-[$α$/Fe], super-solar metallicity stars, reminiscent of the age and metallicity of the well-studied open cluster NGC6791. The age-chemo-kinematic properties of these stars indicate that efficient radial migration happens in the thin disk. We find that ages and masses of the nearly 400 $α$-element-rich red-giant-branch (RGB) stars in our sample are compatible with those of an old (~11 Gyr), nearly coeval, chemical-thick disk population. Using a statistical model, we show that 95% of the population was born within ~1.5 Gyr. Moreover, we find a difference in the vertical velocity dispersion between low- and high-[$α$/Fe] populations, confirming their different chemo-dynamical histories. We then exploit the almost coeval $α$-rich population to gain insight into processes that may have altered the mass of a star along its evolution, which are key to improve the mapping of the observed stellar mass to age. We find evidence for a mean integrated RGB mass loss <$Δ$M>= 0.10 $\pm$ 0.02 Msun and that the occurrence of massive (M $\gtrsim$ 1.1 Msun) $α$-rich stars is of the order of 5% on the RGB, and significantly higher in the RC, supporting the scenario in which most of these stars had undergone interaction with a companion.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Explaining the chemical trajectories of accreted and in-situ halo stars of the Milky Way

The Milky Way underwent its last significant merger ten billion years ago, when the Gaia-Enceladus-Sausage (GES) was accreted. Accreted GES stars and progenitor stars born prior to the merger make up the bulk of the inner halo. Even though these two main populations of halo stars have similar $durations$ of star formation prior to their merger, they differ in [$α$/Fe]-[Fe/H] space, with the GES population bending to lower [$α$/Fe] at a relatively low value of [Fe/H]. We use cosmological simulations of a 'Milky Way' to argue that the different tracks of the halo stars through the [$α$/Fe]-[Fe/H] plane are due to a difference in their star formation history and efficiency, with the lower mass GES having its low and constant star formation regulated by feedback whilst the higher mass main progenitor has a higher star formation rate prior to the merger. The lower star formation efficiency of GES leads to lower gas pollution levels, pushing [$α$/Fe]-[Fe/H] tracks to the left. In addition, the increasing star formation rate maintains a higher relative contribution of Type~II SNe to Type~Ia SNe for the main progenitor population that formed during the same time period, thus maintaining a relatively high [$α$/Fe]. Thus the different positions of the downturns in the [$α$/Fe]-[Fe/H] plane for the GES stars are not reflective of different star formation durations, but instead reflect different star formation efficiencies. We argue that cosmological simulations match a wide range of independent observations, breaking degeneracies that exist in simpler models.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Sausage & Mash: The dual origin of the Galactic thick disc and halo from the gas-rich Gaia-Enceladus-Sausage merger

We analyse a set of cosmological magneto-hydrodynamic simulations of the formation of Milky Way-mass galaxies identified to have a prominent radially anisotropic stellar halo component similar to the so-called "Gaia Sausage" found in the Gaia data. We examine the effects of the progenitor of the Sausage (the Gaia-Enceladus-Sausage, GES) on the formation of major galactic components analogous to the Galactic thick disc and inner stellar halo. We find that the GES merger is likely to have been gas-rich and contribute 10-50$\%$ of gas to a merger-induced centrally concentrated starburst that results in the rapid formation of a compact, rotationally supported thick disc that occupies the typical chemical thick disc region of chemical abundance space. We find evidence that gas-rich mergers heated the proto-disc of the Galaxy, scattering stars onto less-circular orbits such that their rotation velocity and metallicity positively correlate, thus contributing an additional component that connects the Galactic thick disc to the inner stellar halo. We demonstrate that the level of kinematic heating of the proto-galaxy correlates with the kinematic state of the population before the merger, the progenitor mass and orbital eccentricity of the merger. Furthermore, we show that the mass and time of the merger can be accurately inferred from local stars on counter-rotating orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The power of co-ordinate transformations in dynamical interpretations of Galactic structure

$Gaia$ DR2 has provided an unprecedented wealth of information about the positions and motions of stars in our Galaxy, and has highlighted the degree of disequilibria in the disc. As we collect data over a wider area of the disc it becomes increasingly appealing to start analysing stellar actions and angles, which specifically label orbit space, instead of their current phase space location. Conceptually, while $\bar{x}$ and $\bar{v}$ tell us about the potential and local interactions, grouping in action puts together stars that have similar frequencies and hence similar responses to dynamical effects occurring over several orbits. Grouping in actions and angles refines this further to isolate stars which are travelling together through space and hence have shared histories. Mixing these coordinate systems can confuse the interpretation. For example, it has been suggested that by moving stars to their guiding radius, the Milky Way spiral structure is visible as ridge-like overdensities in the $Gaia$ data \citep{Khoperskov+19b}. However, in this work, we show that these features are in fact the known kinematic moving groups, both in the $L_z-ϕ$ and the $v_{\mathrm{R}}-v_ϕ$ planes. Using simulations we show how this distinction will become even more important as we move to a global view of the Milky Way. As an example, we show that the radial velocity wave seen in the Galactic disc in $Gaia$ and APOGEE should become stronger in the action-angle frame, and that it can be reproduced by transient spiral structure.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Towards Machine-assisted Meta-Studies: The Hubble Constant

We present an approach for automatic extraction of measured values from the astrophysical literature, using the Hubble constant for our pilot study. Our rules-based model -- a classical technique in natural language processing -- has successfully extracted 298 measurements of the Hubble constant, with uncertainties, from the 208,541 available arXiv astrophysics papers. We have also created an artificial neural network classifier to identify papers in arXiv which report novel measurements. From the analysis of our results we find that reporting measurements with uncertainties and the correct units is critical information when distinguishing novel measurements in free text. Our results correctly highlight the current tension for measurements of the Hubble constant and recover the $3.5σ$ discrepancy -- demonstrating that the tool presented in this paper is useful for meta-studies of astrophysical measurements from a large number of publications.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Asymmetric mass models of disk galaxies - I. Messier 99

Mass models of galactic disks traditionally rely on axisymmetric density and rotation curves, paradoxically acting as if their most remarkable asymmetric features, such as lopsidedness or spiral arms, were not important. In this article, we relax the axisymmetry approximation and introduce a methodology that derives 3D gravitational potentials of disk-like objects and robustly estimates the impacts of asymmetries on circular velocities in the disk midplane. Mass distribution models can then be directly fitted to asymmetric line-of-sight velocity fields. Applied to the grand-design spiral M99, the new strategy shows that circular velocities are highly nonuniform, particularly in the inner disk of the galaxy, as a natural response to the perturbed gravitational potential of luminous matter. A cuspy inner density profile of dark matter is found in M99, in the usual case where luminous and dark matter share the same center. The impact of the velocity nonuniformity is to make the inner profile less steep, although the density remains cuspy. On another hand, a model where the halo is core dominated and shifted by 2.2-2.5 kpc from the luminous mass center is more appropriate to explain most of the kinematical lopsidedness evidenced in the velocity field of M99. However, the gravitational potential of luminous baryons is not asymmetric enough to explain the kinematical lopsidedness of the innermost regions, irrespective of the density shape of dark matter. This discrepancy points out the necessity of an additional dynamical process in these regions: possibly a lopsided distribution of dark matter.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Metal diffusion in smoothed particle hydrodynamics simulations of dwarf galaxies

We perform a series of smoothed particle hydrodynamics simulations of isolated dwarf galaxies to compare different metal mixing models. In particular, we examine the role of diffusion in the production of enriched outflows, and in determining the metallicity distributions of gas and stars. We investigate different diffusion strengths, by changing the pre-factor of the diffusion coefficient, by varying how the diffusion coefficient is calculated from the local velocity distribution, and by varying whether the speed of sound is included as a velocity term. Stronger diffusion produces a tighter [O/Fe]-[Fe/H] distribution in the gas, and cuts off the gas metallicity distribution function at lower metallicities. Diffusion suppresses the formation of low-metallicity stars, even with weak diffusion, and also strips metals from enriched outflows. This produces a remarkably tight correlation between "metal mass-loading" (mean metal outflow rate divided by mean metal production rate) and the strength of diffusion, even when the diffusion coefficient is calculated in different ways. The effectiveness of outflows at removing metals from dwarf galaxies and the metal distribution of the gas is thus dependent on the strength of diffusion. By contrast, we show that the metallicities of stars are not strongly dependent on the strength of diffusion, provided that some diffusion is present.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spiral-induced velocity and metallicity patterns in a cosmological zoom simulation of a Milky Way-sized galaxy

We use a high resolution cosmological zoom simulation of a Milky Way-sized halo to study the observable features in velocity and metallicity space associated with the dynamical influence of spiral arms. For the first time, we demonstrate that spiral arms, that form in a disc in a fully cosmological environment with realistic galaxy formation physics, drive large-scale systematic streaming motions. In particular, on the trailing edge of the spiral arms the peculiar galacto-centric radial and azimuthal velocity field is directed radially outward and azimuthally backward, whereas it is radially inward and azimuthally forward on the leading edge. Owing to the negative radial metallicity gradient, this systematic motion drives, at a given radius, an azimuthal variation in the residual metallicity that is characterised by a metal rich trailing edge and a metal poor leading edge. We show that these signatures are theoretically observable in external galaxies with Integral Field Unit instruments such as VLT/MUSE, and if detected, would provide evidence for large-scale systematic radial migration driven by spiral arms.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The effects of bar-spiral coupling on stellar kinematics in the Galaxy

We investigate models of the Milky Way disc taking into account simultaneously the bar and a two-armed quasi-static spiral pattern. Away from major resonance overlaps, the mean stellar radial motions in the plane are essentially a linear superposition of the isolated effects of the bar and spirals. Thus, provided the bar is strong enough, even in the presence of spiral arms, these mean radial motions are predominantly affected by the Galactic bar for large scale velocity fluctuations. This is evident when comparing the peculiar line-of-sight velocity power spectrum of our coupled models with bar-only models. However, we show how forthcoming spectroscopic surveys could disentangle bar-only non-axisymmetric models of the Galaxy from models in which spiral arms have a significant amplitude. We also point out that overlaps of low-order resonances are sufficient to enhance stellar churning within the disc, even when the spirals amplitude is kept constant. Nevertheless, for churning to be truly non-local, stronger or (more likely) transient amplitudes would be needed: otherwise the disc is actually mostly unaffected by churning in the present models. Finally, regarding vertical breathing modes, the combined effect of the bar and spirals on vertical motions is a clear non-linear superposition of the isolated effects of both components, significantly superseding the linear superposition of modes produced by each perturber separately, thereby providing an additional effect to consider when analysing the observed breathing mode of the Galactic disc in the extended Solar neighbourhood.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Mass Dependence of Star Formation Histories in Barred Spiral Galaxies

We performed a series of 29 gasdynamical simulations of disc galaxies, barred and unbarred, with various stellar masses, to study the impact of the bar on star formation history. Unbarred galaxies evolve very smoothly, with a star formation rate (SFR) that varies by at most a factor of three over a period of 2 Gyr. The evolution of barred galaxies is much more irregular, especially at high stellar masses. In these galaxies, the bar drives a substantial amount of gas toward the centre, resulting in a high SFR, and producing a starburst in the most massive galaxies. Most of the gas is converted into stars, and gas exhaustion leads to a rapid drop of star formation after the starburst. In massive barred galaxies (stellar mass M* > 2x10^10 Msun) the large amount of gas funnelled toward the centre is completely consumed by the starburst, while in lower-mass barred galaxies it is only partially consumed. Gas concentration is thus higher in lower-mass barred galaxies than it is in higher-mass ones. Even though unbarred galaxies funnelled less gas toward their centre, the lower SFR allows this gas to accumulate. At late times, the star formation efficiency is higher in barred galaxies than unbarred ones, enabling these galaxies to maintain a higher SFR with a smaller gas supply. Several properties, such as the global SFR, central SFR, or central gas concentration, vary monotonically with time for unbarred galaxies, but not for barred galaxies. Therefore one must be careful when comparing barred and unbarred galaxies that share one observational property, since these galaxies might be at very different stages of their respective evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Rotation-Metallicity Relation for the Galactic Disk as Measured in the Gaia DR1 TGAS and APOGEE Data

Previous studies have found that the Galactic rotation velocity-metallicity (V-[Fe/H]) relations for the thin and thick disk populations show negative and positive slopes, respectively. The first Gaia Data Release includes the Tycho-Gaia Astrometric Solution (TGAS) information, which we use to analyze the V-[Fe/H] relation for a strictly selected sample with high enough astrometric accuracy. We aim to arrive at an explanation for the slopes of the V-[Fe/H] relationship. We measure the V-[Fe/H] relation for thin and thick disk stars classified on the basis of their [$α$/Fe] and [Fe/H] abundances. We find dV/d[Fe/H]= -18 +/- 2 km/s/dex for stars in the thin disk and dV/d[Fe/H]= +23 +/- 10 km/s/dex for thick disk stars, so we confirm the different signs for the slopes. The negative value of dV/d[Fe/H] for thick disk stars is consistent with previous studies, but the combination of TGAS and APOGEE data provide higher precision, even though systematic errors could exceed +/-5 km/s/dex . Our average measurement of dV/d[Fe/H] for local thick disk stars shows a somewhat flatter slope than the previous studies, but we confirm a significant spread and a dependence of the slope on the [alpha/Fe] ratio of the stars. Using a simple N-body model, we demonstrate that the observed trend for the thick and thin disk can be explained by the observed radial metallicity gradients and the correlation between orbital eccentricity and metallicity in the thick disk.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Impact of radial migration on stellar and gas radial metallicity distribution

Radial migration is defined as the change in guiding centre radius of stars and gas caused by gains or losses of angular momentum that result from gravitational interaction with non-axisymmetric structure. This has been shown to have significant impact on the metallicity distribution in galactic discs, and therefore affects the interpretation of Galactic archeology. We use a simulation of a Milky Way-sized galaxy to examine the effect of radial migration on the star and gas radial metallicity distribution. We find that both the star and gas component show significant radial migration. The stellar radial metallicity gradient remains almost unchanged but the radial metallicity distribution of the stars is broadened to produce a greater dispersion at all radii. However, the metallicity dispersion of the gas remains narrow. We find that the main drivers of the gas metallicity distribution evolution are metal enrichment and mixing: more efficient metal enrichment in the inner region maintains a negative slope in the radial metallicity distribution, and the metal mixing ensures the tight relationship of the gas metallicity with the radius. The metallicity distribution function reproduces the trend in the age-metallicity relation found from observations for stars younger than 1.0 Gyr in the Milky Way.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radial migration in numerical simulations of Milky Way-sized galaxies

We show that in N-body simulations of isolated spiral discs, spiral arms appear to transient, recurring features that co-rotate with the stellar disc stars at all radii. As a consequence, stars around the spiral arm continually feel a tangential force from the spiral and gain/lose angular momentum at all radii where spiral structure exists, without gaining significant amounts of random energy. We demonstrate that the ubiquitous radial migration in these simulations can be seen as outward (inward) systematic streaming motions along the trailing (leading) side of the spiral arms. We characterise these spiral induced peculiar motions and compare with those of the Milky Way obtained from APOGEE red clump data. We find that transient, co-rotating spiral arms are consistent with the data, in contrast with density wave-like spirals which are qualitatively inconsistent. In addition, we show that, in our simulations, radial migration does not change the radial metallicity gradient significantly, and broadens the metallicity distribution function at all radii.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Spiral and bar driven peculiar velocities in Milky Way sized galaxy simulations

We investigate the kinematic signatures induced by spiral and bar structure in a set of simulations of Milky Way-sized spiral disc galaxies. The set includes test particle simulations that follow a quasi-stationary density wave-like scenario with rigidly rotating spiral arms, and $N$-body simulations that host a bar and transient, co-rotating spiral arms. From a location similar to that of the Sun, we calculate the radial, tangential and line-of-sight peculiar velocity fields of a patch of the disc and quantify the fluctuations by computing the power spectrum from a two-dimensional Fourier transform. We find that the peculiar velocity power spectrum of the simulation with a bar and transient, co-rotating spiral arms fits very well to that of APOGEE red clump star data, while the quasi-stationary density wave spiral model without a bar does not. We determine that the power spectrum is sensitive to the number of spiral arms, spiral arm pitch angle and position with respect to the spiral arm. However, it is necessary to go beyond the line of sight velocity field in order to distinguish fully between the various spiral models with this method. We compute the power spectrum for different regions of the spiral discs, and discuss the application of this analysis technique to external galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The stellar kinematics of co-rotating spiral arms in Gaia mock observations

We have observed an N-body/Smoothed Particle Hydrodynamics simulation of a Milky Way like barred spiral galaxy. We present a simple method that samples N-body model particles into mock Gaia stellar observations and takes into account stellar populations, dust extinction and Gaia's science performance estimates. We examine the kinematics around a nearby spiral arm at a similar position to the Perseus arm at three lines of sight in the disc plane; (l,b)=(90,0), (120,0) and (150,0) degrees. We find that the structure of the peculiar kinematics around the co-rotating spiral arm, which is found in Kawata et al. (2014b), is still visible in the observational data expected to be produced by Gaia despite the dust extinction and expected observational errors of Gaia. These observable kinematic signatures will enable testing whether the Perseus arm of the Milky Way is similar to the co-rotating spiral arms commonly seen in N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gas and Stellar Motions and Observational Signatures of Co-Rotating Spiral Arms

We have observed a snapshot of our N-body/Smoothed Particle Hydrodynamics simulation of a Milky Way-sized barred spiral galaxy in a similar way to how we can observe the Milky Way. The simulated galaxy shows a co-rotating spiral arm, i.e. the spiral arm rotates with the same speed as the circular speed. We observed the rotation and radial velocities of the gas and stars as a function of the distance from our assumed location of the observer at the three lines of sight on the disc plane, (l, b) = (90, 0), (120, 0) and (150,0) deg. We find that the stars tend to rotate slower (faster) behind (at the front of) the spiral arm and move outward (inward), because of the radial migration. However, because of their epicycle motion, we see a variation of rotation and radial velocities around the spiral arm. On the other hand, the cold gas component shows a clearer trend of rotating slower (faster) and moving outward (inward) behind (at the front of) the spiral arm, because of the radial migration. We have compared the results with the velocity of the maser sources from Reid et al. (2014), and find that the observational data show a similar trend in the rotation velocity around the expected position of the spiral arm at l = 120 deg. We also compared the distribution of the radial velocity from the local standard of the rest, V_LSR, with the APOGEE data at l = 90 deg as an example.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

M2M modelling of the Galactic disc via PRIMAL: Fitting to Gaia Error Added Data

We have adapted our made-to-measure (M2M) algorithm PRIMAL to use mock Milky Way like data constructed from an N-body barred galaxy with a boxy bulge in a known dark matter potential. We use M0 giant stars as tracers, with the expected error of the ESA space astrometry mission Gaia. We demonstrate the process of constructing mock Gaia data from an N-body model, including the conversion of a galactocentric Cartesian coordinate N-body model into equatorial coordinates and how to add error to it for a single stellar type. We then describe the modifications made to PRIMAL to work with observational error. This paper demonstrates that PRIMAL can recover the radial profiles of the surface density, radial velocity dispersion, vertical velocity dispersion and mean rotational velocity of the target disc, along with the pattern speed of the bar, to a reasonable degree of accuracy despite the lack of accurate target data. We also construct mock data which take into account dust extinction and show that PRIMAL recovers the structure and kinematics of the disc reasonably well. In other words, the expected accuracy of the Gaia data is good enough for PRIMAL to recover these global properties of the disc, at least in a simplified condition, as used in this paper.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Numerical Simulation of a possible origin of the positive radial metallicity gradient of the thick disk

We analyze the radial and vertical metallicity and [alpha/Fe] gradients of the disk stars of a disk galaxy simulated in a fully cosmological setting with the chemodynamical galaxy evolution code, GCD+. We study how the radial abundance gradients vary as a function of height above the plane and find that the metallicity ([alpha/Fe]) gradient becomes more positive (negative) with increasing height, changing sign around 1.5 kpc above the plane. At the largest vertical height (2 < |z| < 3 kpc), our simulated galaxy shows a positive radial metallicity gradient. We find that the positive metallicity gradient is caused by the age-metallicity and age-velocity dispersion relation, where the younger stars have higher metallicity and lower velocity dispersion. Due to the age-velocity dispersion relation, a greater fraction of younger stars reach |z| > 2 kpc at the outer region, because of the lower gravitational restoring force of the disk, i.e. flaring. As a result, the fraction of younger stars with higher metallicity due to the age-metallicity relation becomes higher at the outer radii, which makes the median metallicity higher at the outer radii. Combining this result with the recently observed age-metallicity and age-velocity dispersion relation for the Milky Way thick disk stars suggested by Haywood et al. (2013), we argue that the observed (small) positive radial metallicity gradient at large heights of the Milky Way disk stars can be explained by the flaring of the younger thick and/or thin disk stars.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Orbits of radial migrators and non-migrators around a spiral arm in N-body simulations

Recent numerical N-body simulations of spiral galaxies have shown that spiral arms in N-body simulations seem to rotate at a similar speed to the local rotation speed of the stellar disc material. This in turn yields winding, transient and recurrent spiral structure, whose co-rotating nature gives rise to changes in the angular momentum (radial migration) of star particles close to the spiral arm at many radii. From high resolution N-body simulations, we highlight the evolution of strongly migrating star particles (migrators) and star particles that do not migrate (non-migrators) around a spiral arm. We investigate the individual orbit histories of migrators and non-migrators and find that there are several types of migrator and non-migrator, each with unique radial evolution. We contrast each type of orbit to establish the reasons for the differences between them. We find that the positive (negative) migrators sustain a position behind (in front of) the spiral arm, and feel continuous tangential force as long as the spiral arm persists. On the other hand, non-migrators stay close to the spiral arm, and pass or are passed by the spiral arm one or two times. Although they gain or lose the angular momentum when they are behind or in front of the spiral arm, their net angular momentum change becomes close to zero. We discuss also the long term effects of radial migration on the radial metallicity distribution and radial angular momentum and mass profiles. [Abridged]

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Disc galaxy modelling with a particle-by-particle M2M method

We have developed the initial version of a new particle-by-particle adaptation of the made-to-measure (M2M) method, aiming to model the Galactic disc from upcoming Galactic stellar survey data. In our new particle-by-particle M2M, the observables of the target system are compared with those of the model galaxy at the position of the target stars (i.e. particles). The weights of the model particles are changed to reproduce the observables of the target system, and the gravitational potential is automatically adjusted by the changing weights of the particles. This paper demonstrates, as the initial work, that the particle-by-particle M2M can recreate a target disc system created by an N-body simulation in a known dark matter potential, with no error in the observables. The radial profiles of the surface density, velocity dispersion in the radial and perpendicular directions, and the rotational velocity of the target disc are all well reproduced from the initial disc model, whose scale length is different from that of the target disc. We also demonstrate that our M2M can be applied to an incomplete data set and recreate the target disc reasonably well when the observables are restricted to a part of the disc. We discuss our calibration of the model parameters and the importance of regularization.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Investigating Bar Structure of Disc Galaxies via PRIMAL: A PaRtIcle-by-particle M2M ALgorithm

We have modified our particle-by-particle adaptation of the made-to-measure (M2M) method, with the aim of modelling the Galactic disc from upcoming Galactic stellar survey data. In our new particle-by-particle M2M algorithm, PRIMAL, the observables of the target system are compared with those of the model galaxy at the position of the target stars, i.e. particles. The mass of the model particles are adjusted to reproduce the observables of the target system, and the gravitational potential is automatically adjusted by the changing mass of the particles. This paper builds upon our previous work, introducing likelihood-based velocity constraints in PRIMAL. In this paper we apply PRIMAL to barred disc galaxies created by a N-body simulation in a known dark matter potential, with no error in the observables. This paper demonstrates that PRIMAL can recover the radial profiles of the surface density, velocity dispersion in the radial and perpendicular directions, and the rotational velocity of the target discs, along with the apparent bar structure and pattern speed of the bar, especially when the reference frame is adjusted so that the bar angle of the target galaxy is aligned to that of the model galaxy at every timestep.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Numerical Simulations of Bubble-Induced Star Formation in Dwarf Irregular Galaxies with a Novel Stellar Feedback Scheme

To study the star formation and feedback mechanism, we simulate the evolution of an isolated dwarf irregular galaxy (dIrr) in a fixed dark matter halo, similar in size to WLM, using a new stellar feedback scheme. We use the new version of our original N-body/smoothed particle chemodynamics code, GCD+, which adopts improved hydrodynamics, metal diffusion between the gas particles and new modelling of star formation and stellar wind and supernovae (SNe) feedback. Comparing the simulations with and without stellar feedback effects, we demonstrate that the collisions of bubbles produced by strong feedback can induce star formation in a more widely spread area. We also demonstrate that the metallicity in star forming regions is kept low due to the mixing of the metal-rich bubbles and the metal-poor inter-stellar medium. Our simulations also suggest that the bubble-induced star formation leads to many counter-rotating stars. The bubble-induced star formation could be a dominant mechanism to maintain star formation in dIrrs, which is different from larger spiral galaxies where the non-axisymmetric structures, such as spiral arms, are a main driver of star formation.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Spiral arm pitch angle and galactic shear rate in N-body simulations of disc galaxies

Spiral galaxies are observed to exhibit a range of morphologies, in particular in the shape of spiral arms. A key diagnostic parameter is the pitch angle, which describes how tightly wound the spiral arms are. Observationally and analytically, a correlation between pitch angle and galactic shear rate has been detected. For the first time, we examine whether this effect is detected in N-body simulations by calculating and comparing pitch angles of both individual density waves and overall spiral structure in a suite of N-body simulations. We find that higher galactic shear rates produce more tightly wound spiral arms, both in individual mode patterns (density waves) and in the overall density enhancement. Although the mode pattern pitch angles by construction remain constant with time, the overall logarithmic spiral arm winds over time, which could help to explain the scatter in the relation between pitch angle versus shear seen from observations. The correlation between spiral arm pitch angle and galactic shear rate that we find in N-body simulations may also explain why late Hubble type of spiral galaxies tend to have more open arms.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Connection between Star Formation and Metallicity Evolution in Barred Spiral Galaxies

We have performed a series of chemodynamical simulations of barred disc galaxies. Our goal is to determine the physical processes responsible for the increase in the central gas-phase metallicity and of the central SFR observed in SDSS. We define a 2kpc diameter central aperture to approximate the integrated spectroscopic fibre measurements from the SDSS. The chemical evolution observed within this central region depends critically upon the relative size of the bar and the aperture, which evolves strongly with time. At t~0.5Gyr, a strong bar forms, whose length is considerably longer than the 2kpc aperture. The stars and gas lose angular momentum and follow elongated orbits that cause an intense mixing of the gas between the central region and its surroundings. During the next 1.5Gyr, the orbits of the gas contract significantly until the entire gas bar is contained in the 2kpc aperture, resulting in a net flux of gas into the central region. During this period, the metallicity in the central region increases steadily, and this enrichment is dominated by metal-rich gas that is flowing into the central region. The main result of this work is therefore that the observed enrichment in the centres of barred galaxies is not dominated by in-situ enrichment by stars formed in the centre. Rather, star formation occurs along the full length of the bar, much of which occurs initially outside the 2kpc aperture. About 50% of the metals that end up in the central region originate from this extended bar-long star formation, but flow into the central region due to loss of angular momentum. We conclude that there is no direct connection between central SFR and central metallicity. The central metallicity does not originate exclusively from central stars. Instead, the global SFR (especially along the bar) and the large-scale flow of enriched gas play a major role.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Calibrating an updated SPH scheme within GCD+

We adapt a modern scheme of smoothed particle hydrodynamics (SPH) to our tree N-body/SPH galactic chemodynamics code GCD+. The applied scheme includes imple- mentations of the artificial viscosity switch and artificial thermal conductivity pro- posed by Morris & Monaghan (1997), Rosswog & Price (2007) and Price (2008), to model discontinuities and Kelvin-Helmholtz instabilities more accurately. We first present hydrodynamics test simulations and contrast the results to runs undertaken without artificial viscosity switch or thermal conduction. In addition, we also explore the different levels of smoothing by adopting larger or smaller smoothing lengths, i.e. a larger or smaller number of neighbour particles, Nnb. We demonstrate that the new version of GCD+ is capable of modelling Kelvin-Helmholtz instabilities to a simi- lar level as the mesh code, Athena. From the Gresho vortex, point-like explosion and self-similar collapse tests, we conclude that setting the smoothing length to keep the number of neighbour particles as high as Nnb~58 is preferable to adopting smaller smoothing lengths. We present our optimised parameter sets from the hydrodynamics tests.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Dynamics of stars around spiral arms in an N-body/SPH simulated barred-spiral galaxy

We run N-body smoothed particle hydrodynamics (SPH) simulations of a Milky Way sized galaxy. The code takes into account hydrodynamics, self-gravity, star formation, supernova and stellar wind feedback, radiative cooling and metal enrichment. The simulated galaxy is a barred-spiral galaxy consisting of a stellar and gas disc, enveloped in a static dark matter halo. Similar to what is found in our pure N-body simulation of a non-barred galaxy in Grand et. al. (2012), we find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is similar to the rotation of star particles. Compared to the non-barred case, we find that the spiral arm pattern speed is slightly faster than the rotation speed of star particles: the bar appears to boost the pattern speed ahead of the rotational velocity. We trace particle motion around the spiral arms at different radii, and demonstrate that there are star particles that are drawn towards and join the arm from behind (in front of) the arm and migrate toward the outer (inner) regions of the disc until the arm disappears as a result of their transient nature. We see this migration over the entire radial range analysed, which is a consequence of the spiral arm rotating at similar speeds to star particles at all radii, which is inconsistent with the prediction of classical density wave theory. The bar does not prevent this systematic radial migration, which is shown to largely preserve circular orbits. We also demonstrate that there is no significant offset of different star forming tracers across the spiral arm, which is also inconsistent with the prediction of classical density wave theory.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

RAMSES-CH: A New Chemodynamical Code for Cosmological Simulations

We present a new chemodynamical code - Ramses-CH - for use in simulating the self-consistent evolution of chemical and hydrodynamical properties of galaxies within a fully cosmological framework. We build upon the adaptive mesh refinement code Ramses, which includes a treatment of self-gravity, hydrodynamics, star formation, radiative cooling, and supernovae feedback, to trace the dominant isotopes of C, N, O, Ne, Mg, Si, and Fe. We include the contribution of Type Ia and II supernovae, in addition to low- and intermediate-mass asymptotic giant branch stars, relaxing the instantaneous recycling approximation. The new chemical evolution modules are highly flexible and portable, lending themselves to ready exploration of variations in the underpining stellar and nuclear physics. We apply Ramses-CH to the cosmological simulation of a typical L* galaxy, demonstrating the successful recovery of the basic empirical constraints regarding, [α/Fe]-[Fe/H] and Type Ia/II supernovae rates.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Swift/UVOT catalogue of NGC4321 star forming sources: A case against density wave theory

We study the star forming regions in the spiral galaxy NGC4321, taking advantage of the spatial resolution (2.5 arcsec FWHM) of the Swift/UVOT camera and the availability of three UV passbands in the region 1600-3000 A, in combination with optical and IR imaging from SDSS, KPNO/Ha and Spitzer/IRAC, to obtain a catalogue of 787 star forming regions out to three disc scale lengths. We determine the properties of the young stellar component and its relationship with the spiral arms. The Ha luminosities of the sources have a strong decreasing radial trend, suggesting more massive star forming regions in the central part of the galaxy. When segregated with respect to NUV-optical colour, blue sources have a significant excess of flux in the IR at 8 micron, revealing the contribution from PAHs, although the overall reddening of these sources stays below E(B-V)=0.2 mag. The distribution of distances to the spiral arms is compared for subsamples selected according to Ha luminosity, NUV-optical colour, or ages derived from a population synthesis model. An offset is expected between these subsamples as a function of radius if the pattern speed of the spiral arm were constant - as predicted by classic density wave theory. No significant offsets are found, favouring instead a mechanism where the pattern speed has a radial dependence.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Comparison of Galaxy Group Luminosity Functions from Semi-Analytic Models

Semi-analytic models (SAMs) are currently one of the primary tools with which to model statistically significant ensembles of galaxies. The underlying physical prescriptions inherent to each SAM are, in many cases, different from one another. Several SAMs have been applied to the dark matter merger trees extracted from the Millennium Run, including those associated with the well-known Munich and Durham lineages. We compare the predicted luminosity distributions of galaxy groups using four publicly available SAMs (De Lucia et al. 2006; Bower et al. 2006; Bertone et al. 2007; Font et al. 2008), in order to explore a galactic environment in which the models have not been explored to the same degree as they have in the field or in rich clusters. We identify a characteristic "wiggle" in the group galaxy luminosity function generated using the De Lucia et al. (2006) SAM, that is not present in the Durham-based models, consistent to some degree with observations. However, a comparison between conditional luminosity functions of groups between the models and observations of Yang et al. (2007) suggest that neither model is a particularly good match. The luminosity function wiggle is interpreted as the result of the two-mode AGN feedback implementation used in the Munich models, which itself results in flattened magnitude gap distribution. An associated analysis of the magnitude gap distribution between first- and second-ranked group galaxies shows that while the Durham models yield distributions with approximately equal luminosity first- and second-ranked galaxies, in agreement with observations, the De Lucia et al. models favours the scenario in which the second-ranked galaxy is approximately one magnitude fainter than the primary,especially when the dynamic range of the mock data is limited to 3 magnitudes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Chemical signature of gas-rich disc-disc mergers at high redshift

We performed numerical simulations of mergers between gas-rich disc galaxies, which result in the formation of late-type galaxies. Stars formed during the merger end up in a thick disc that is partially supported by velocity dispersion and has high [alpha/Fe] ratios at all metallicities. Stars formed later end up in a thin, rotationally supported disc which has lower [alpha/Fe] ratios. While the structural and kinematical properties of the merger remnants depend strongly upon the orbital parameters of the mergers, we find a clear chemical signature of gas-rich mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Cosmological Simulations using GCMHD+

Radio observations of galaxy clusters show that the intra cluster medium is permeated by μG magnetic fields. The origin and evolution of these cosmological magnetic fields is currently not well understood and so their impact on the dynamics of structure formation is not known. Numerical simulations are required to gain a greater understanding and produce predictions for the next generation of radio telescopes. We present the galactic chemodynamics smoothed particle magnetohydrodynamic (SPMHD) code (GCMHD+), which is an MHD implementation for the cosmological smoothed particle hydrodynamic code GCD+. The results of 1, 2 and 3 dimensional tests are presented and the performance of the code is shown relative to the ATHENA grid code. GCMHD+ shows good agreement with the reference solutions produced by ATHENA. The code is then used to simulate the formation of a galaxy cluster with a simple primordial magnetic field embedded in the gas. A homogeneous seed field of 10^-11 G is amplified by a factor of 10^3 during the formation of the cluster. The results show good agreement with the profiles found in other magnetic cluster simulations of similar resolution.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Disc Heating: Comparing the Milky Way with Cosmological Simulations

We present the analysis of a suite of simulations run with different particle-and grid-based cosmological hydrodynamical codes and compare them with observational data of the Milky Way. This is the first study to make comparisons of properties of galaxies simulated with particle and grid-based codes. Our analysis indicates that there is broad agreement between these different modelling techniques. We study the velocity dispersion - age relation for disc stars at z=0 and find that four of the simulations are more consistent with observations by Holmberg et al. (2008) in which the stellar disc appears to undergo continual/secular heating. Two other simulations are in better agreement with the Quillen & Garnett (2001) observations that suggest a "saturation" in the heating profile for young stars in the disc. None of the simulations have thin discs as old as that of the Milky Way. We also analyse the kinematics of disc stars at the time of their birth for different epochs in the galaxies' evolution and find that in some simulations old stars are born cold within the disc and are subsequently heated, while other simulations possess old stellar populations which are born relatively hot. The models which are in better agreement with observations of the Milky Way's stellar disc undergo significantly lower minor-merger/assembly activity after the last major merger - i.e. once the disc has formed. All of the simulations are significantly "hotter" than the Milky Way disc; on top of the effects of mergers, we find a "floor" in the dispersion that is related to the underlying treatment of the heating and cooling of the interstellar medium, and the low density threshold which such codes use for star formation. This finding has important implications for all studies of disc heating that use hydrodynamical codes.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Metallicity gradients of disc stars for a cosmologically simulated galaxy

We analyse for the first time the radial abundance gradients of the disc stars of a disc galaxy simulated with our three dimensional, fully cosmological chemodynamical galaxy evolution code GCD+. We study how [Fe/H], [N/O], [O/Fe], [Mg/Fe] and [Si/Fe] vary with galactocentric radius. For the young stars of the disc, we found a negative slope for [Fe/H] and [N/O] but a positive [O/Fe], [Mg/Fe] and [Si/Fe] slope with radius. By analysing the star formation rate (SFR) at different radii, we found that the simulated disc contains a greater fraction of young stars in the outer regions, while the old stars tend to be concentrated in the inner parts of the disc. This can explain the positive [alpha/Fe] gradient as well as the negative [N/O] gradient with radius. This radial trend is a natural outcome of an inside-out formation of the disc, regardless of its size and can thus explain the recently observed positive [alpha/Fe] gradients in the Milky Way disc open clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Stellar dynamics around transient co-rotating spiral arms

Spiral density wave theory attempts to describe the spiral pattern in spiral galaxies in terms of a long-lived wave structure with a constant pattern speed in order to avoid the winding dilemma. The pattern is consequently a rigidly rotating, long-lived feature. We run an N-body/SPH simulation of a Milky Way-sized barred disk, and find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is almost equal to the rotation curve of the galaxy. We trace particle motion around the spiral arms. We show that particles from behind and in front of the spiral arm are drawn towards and join the arm. Particles move along the arm in the radial direction and we find a clear trend that they migrate toward the outer (inner) radii on the trailing (leading) side of the arm. Our simulations demonstrate that at all radii where there is a co-rotating spiral arm the particles continue to be accelerated (decelerated) by the spiral arm for long periods, which leads to strong migration.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The dynamics of stars around spiral arms

Spiral density wave theory attempts to describe the spiral pattern in spiral galaxies in terms of a long-lived wave structure with a constant pattern speed in order to avoid the winding dilemma. The pattern is consequently a rigidly rotating, long-lived feature. We run N-body simulations of a giant disc galaxy consisting of a pure stellar disc and a static dark matter halo, and find that the spiral arms are transient features whose pattern speeds decrease with radius, in such a way that the pattern speed is almost equal to the rotation curve of the galaxy. We trace particle motion around the spiral arms. We show that particles from behind and in front of the spiral arm are drawn towards and join the arm. Particles move along the arm in the radial direction and we find a clear trend that they migrate toward the outer (inner) radii on the trailing (leading) side of the arm. Our simulations demonstrate that because the spiral arm feature is co-rotating, the particles continue to be accelerated (decelerated) by the spiral arm for long periods, which leads to strong and efficient migration, at all radii in the disc.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Towards a self-consistent numerical model of late-type galaxies: Calibrating the effects of sub-grid physics on galactic models

We carry out several isolated galaxy evolution simulations in a fixed dark matter halo gravitational potential using the new version of our N-body/Smoothed Particle Hydrodynamics (SPH) code GCD+. The new code allows us to more accurately model and follow the evolution of the gas and stellar components of the system including powerful supernovae feedback and its effects on the inter-stellar medium (ISM). Here we present the results of six simulations of an M33-sized late-type disc galaxy each with varying values for our model parameters which include the star formation efficiency (C*), the energy released per supernovae explosion (E_SN) and the energy released per unit time from stellar winds (E_SW). We carry out both a pixel-by-pixel and radial ring analysis method for each of our galaxies comparing our results to the observed Schmidt-Kennicutt Law and vertical gas velocity dispersion versus radius relation amongst others. We find that our models with higher feedback more closely resemble the observations and that feedback plays a pivotal role in obtaining both the observed Schmidt-Kennicutt and gas velocity dispersion relations.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Are Dry Mergers Dry, Moist, or Wet?

We present a spectral analysis of a sample of red-sequence galaxies identified by van Dokkum (2005) as dry merger remnants and ongoing dry mergers. Kinematics, stellar population absorption features, and ionisation from emission lines, are derived. We find that approximately half of the sample showing strong tidal features have younger stellar populations than a control sample at a given velocity dispersion. Conversely, galaxies with weak tidal tails and/or ongoing mergers -- with the exception of one galaxy -- do not show this young component. This seems to indicate that the young stellar populations observed in a significant fraction of ellipticals is the consequence of star formation triggered by mergers. This young component is consistent with a light "frosting" of young stars (<2% by mass) superimposed upon a dominant, old (~11 Gyr), stellar population. In terms of stellar populations, these mergers are, in fact, fairly dry. We found, however, that merger remnants with young stellar populations are supported by rotation, contrary to the expectations of a major dry merger. This suggests that the small amount of gas involved has been sufficient to produce a dynamically cold stellar component. Half of the galaxies with strong tidal distortion, however, are slow rotating and have stellar populations compatible with the control sample at a given velocity dispersion. Remarkably, none of the galaxies with velocity dispersions in excess of 250 km/s have a young stellar component, independent of the merger stage.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Structure, Kinematics, and Chemical Enrichment Patterns after Major Gas-Rich Disc-Disc Mergers

We used an N-body smoothed particle hydrodynamics algorithm, with a detailed treatment of star formation, supernovae feedback, and chemical enrichment, to perform eight simulations of mergers between gas-rich disc galaxies. We vary the mass ratio of the progenitors, their rotation axes, and their orbital parameters and analyze the kinematic, structural, and chemical properties of the remnants. Six of these simulations result in the formation of a merger remnant with a disc morphology as a result of the large gas-fraction of the remnants. We show that stars formed during the merger (a sudden starburst occur in our simulation and last for 0.2-0.3 Gyr) and those formed after the merger have different kinematical and chemical properties. The first ones are located in thick disc or the halo. They are partially supported by velocity dispersion and have high [alpha/Fe] ratios even at metallicities as high as [Fe/H]=-0.5. The former ones -- the young component -- are located in a thin disc rotationally supported and have lower [alpha/Fe] ratios. The difference in the rotational support of both components results in the rotation of the thick disc lagging that of the thin disc by as much as a factor of two, as recently observed.We find that, while the kinematic and structural properties of the merger remnant depends strongly upon the orbital parameters of the mergers, there is a remarkable uniformity in the chemical properties of the mergers. This suggests that general conclusions about the chemical signature of gas-rich mergers can be drawn.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Hydrodynamical Adaptive Mesh Refinement Simulations of Disk Galaxies

To date, fully cosmological hydrodynamic disk simulations to redshift zero have only been undertaken with particle-based codes, such as GADGET, Gasoline, or GCD+. In light of the (supposed) limitations of traditional implementations of smoothed particle hydrodynamics (SPH), or at the very least, their respective idiosyncrasies, it is important to explore complementary approaches to the SPH paradigm to galaxy formation. We present the first high-resolution cosmological disk simulations to redshift zero using an adaptive mesh refinement (AMR)-based hydrodynamical code, in this case, RAMSES. We analyse the temporal and spatial evolution of the simulated stellar disks' vertical heating, velocity ellipsoids, stellar populations, vertical and radial abundance gradients (gas and stars), assembly/infall histories, warps/lopsideness, disk edges/truncations (gas and stars), ISM physics implementations, and compare and contrast these properties with our sample of cosmological SPH disks, generated with GCD+. These preliminary results are the first in our long-term Galactic Archaeology Simulation program.

0 citations0 reviewsNo code linkedOpen access
preprint2008arXiv

Simulating the Mass-Metallicity Relation from z=1

We use 112 N-body/hydrodynamical simulations in the standard Cold Dark Matter universe, to follow the formation of galaxy-sized halos and investigate the chemical enrichment of both the stellar component and the interstellar medium of galaxies, with stellar masses larger than 1e9 Msun. The resulting chemical properties of the simulated galaxies are broadly consistent with the observations. The predicted relationship between the mean metallicity and the galaxy stellar mass for both the stellar and the gaseous components at z=0 are in agreement with the relationships observed locally. The predicted scatter about these relationships, which is traced to the differing merging histories amongst the simulated galaxies with similar final masses, is similar to that observed. The predicted correlations between the total mass and the stellar mass of galaxies in our simulated sample from the present epoch up to z=1 agree with observed ones. The stellar mass versus metallicity relation and its associated scatter are reproduced by the simulations as consequences of the increasing efficiency of the conversion of gas into stars with stellar mass, and the differing merging histories amongst the galaxies with similar masses. The old ages of simulated low mass galaxies at z=0, and the weak level of chemical evolution for massive galaxies suggest however that our modeling of the supernova feedback may be incomplete, or that other feedback processes have been neglected. (Abridged)

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

The Spatial Distribution of the Galactic First Stars II: SPH Approach

We use cosmological, chemo-dynamical, smoothed particle hydrodynamical simulations of Milky-Way-analogue galaxies to find the expected present-day distributions of both metal-free stars that formed from primordial gas and the oldest star populations. We find that metal-free stars continue to form until z~4 in halos that are chemically isolated and located far away from the biggest progenitor of the final system. As a result, if the Population III initial mass function allows stars with low enough mass to survive until z=0 (< 0.8 Msol), they would be distributed throughout the Galactic halo. On the other hand, the oldest stars form in halos that collapsed close to the highest density peak of the final system, and at z=0 they are located preferentially in the central region of the Galaxy, i.e., in the bulge. According to our models, these trends are not sensitive to the merger histories of the disk galaxies or the implementation of supernova feedback. Furthermore, these full hydrodynamics results are consistent with our N-body results in Paper I, and lend further weight to the conclusion that surveys of low-metallicity stars in the Galactic halo can be used to directly constrain the properties of primordial stars. In particular, they suggest that the current lack of detections of metal-free stars implies that their lifetimes were shorter than a Hubble time, placing constraints on the metal-free initial mass function.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Gas Rich Mergers in Disk Formation

In order to explain disk galaxy formation within the hierarchical structure formation, it seems that gas rich mergers must play an important role. We review here our previous studies which have shown the importance of mergers at high redshift being gas rich, in the formation of both the stellar halo and thick disk components of disk galaxies. Regulation of star formation in the building blocks of our galaxy is required to form a low mass low metallicity stellar halo. This regulation results in high redshift, gas rich mergers during which the thick disk forms. In these proceedings, we categorise stars from our simulated disk galaxy into thin and thick disk components by using the Toomre diagram.

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Simulating the Hot X-ray Emitting Gas in Elliptical Galaxies

We study the chemo-dynamical evolution of elliptical galaxies and their hot X-ray emitting gas using high-resolution cosmological simulations. Our Tree N-body/SPH code includes a self-consistent treatment of radiative cooling, star formation, supernovae feedback, and chemical enrichment. We present a series of LCDM cosmological simulations which trace the spatial and temporal evolution of heavy element abundance patterns in both the stellar and gas components of galaxies. X-ray spectra of the hot gas are constructed via the use of the vmekal plasma model, and analysed using XSPEC with the XMM EPN response function. Simulation end-products are quantitatively compared with the observational data in both the X-ray and optical regime. We find that radiative cooling is important to interpret the observed X-ray luminosity, temperature, and metallicity of the interstellar medium of elliptical galaxies. However, this cooled gas also leads to excessive star formation at low redshift, and therefore results in underlying galactic stellar populations which are too blue with respect to observations.

0 citations0 reviewsNo code linkedOpen access
preprint2000arXiv

The Role of Clustering of Sub-Clumps in Bright Elliptical Galaxy Formation from a Low-Spin Seed Galaxy

We reveal the role of clustering of sub-clumps, which is expected in the cold dark matter (CDM) universe, in forming a bright elliptical galaxy (BEG) from a low-spin seed galaxy. This can be done by comparing the evolution of a low-spin seed galaxy including small-scale density fluctuations expected in the CDM universe (Model 1) with that of a completely uniform one (Model 2), using numerical experiments. We show that Model 2 cannot reproduce the properties of BEGs and forms a disk which is too compact and too bright due to the conservation of the initial-small angular momentum. In Model 1 clustering of the sub-clumps caused by initial small-scale density fluctuations leads to angular momentum transfer from the baryon component to the dark matter and consequently a nearly spherical system supported by random motions is formed. Moreover the collisionless property of the stars formed in the sub-clumps prevents the dissipative contraction of the system, leading to a large measured half-light radius. As a result, the end-product is quite well reproduces the observed properties of BEGs, such as the de Vaucouleurs light-profile, typical color and metallicity gradients, the large half-light radius, the small ratio of the rotational velocity to the velocity dispersion (V/σ). We conclude that the clustering of sub-clumps, i.e., the hierarchical clustering, plays a crucial role in the formation of BEGs from a low-spin seed galaxy.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

Galaxy Formation from a Low-Spin Density Perturbation in a CDM Universe

In order to understand the formation process of elliptical galaxies which are not rotationally supported, we have carried out numerical simulations of the galaxy formation from the density perturbation with a rotation corresponding to a small spin parameter. The three-dimensional TREE N-Body/SPH simulation code used in this paper includes the dark matter and gas dynamics, radiative cooling, star formation, supernova feedback, and metal enrichment. The initial condition is a slowly rotating, top-hat over-dense sphere on which the perturbations expected in a cold dark matter (CDM) universe are superposed. By means of the stellar population synthesis, we calculated the surface brightness profile, the metallicity distribution, and the photometric properties of the end-product, and found that these properties quantitatively agree with the observed properties of bright elliptical galaxies. Thus, we conclude that, in a CDM universe, the proto-galaxy which has a spin-parameter as small as 0.02 evolves into an elliptical galaxy.

0 citations0 reviewsNo code linkedOpen access