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.

$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-\phi$ and the $v_{\mathrm{R}}-v_{\phi}$ 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.

[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-[$\alpha$/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 $\alpha$-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-[$\alpha$/Fe] populations, confirming their different chemo-dynamical histories. We then exploit the almost coeval $\alpha$-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 <$\Delta$M>= 0.10 $\pm$ 0.02 Msun and that the occurrence of massive (M $\gtrsim$ 1.1 Msun) $\alpha$-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.

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 [$\alpha$/Fe]-[Fe/H] space, with the GES population bending to lower [$\alpha$/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 [$\alpha$/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 [$\alpha$/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 [$\alpha$/Fe]. Thus the different positions of the downturns in the [$\alpha$/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.

We develop a Bayesian Machine Learning framework called BINGO (Bayesian INference for Galactic archaeOlogy) centred around a Bayesian neural network. After being trained on the APOGEE and \emph{Kepler} asteroseismic age data, BINGO is used to obtain precise relative stellar age estimates with uncertainties for the APOGEE stars. We carefully construct a training set to minimise bias and apply BINGO to a stellar population that is similar to our training set. We then select the 17,305 stars with ages from BINGO and reliable kinematic properties obtained from \textit{Gaia} DR2. By combining the age and chemo-kinematical information, we dissect the Galactic disc stars into three components, namely, the thick disc (old, high-[$\alpha$/Fe], [$\alpha$/Fe] $\gtrsim$ 0.12), the thin disc (young, low-[$\alpha$/Fe]) and the Bridge, which is a region between the thick and thin discs. Our results indicate that the thick disc formed at an early epoch only in the inner region, and the inner disc smoothly transforms to the thin disc. We found that the outer disc follows a different chemical evolution pathway from the inner disc. The outer metal-poor stars only start forming after the compact thick disc phase has completed and the star-forming gas disc extended outwardly with metal-poor gas accretion. We found that in the Bridge region the range of [Fe/H] becomes wider with decreasing age, which suggests that the Bridge region corresponds to the transition phase from the smaller chemically well-mixed thick to a larger thin disc with a metallicity gradient.

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.

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.

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\sigma$ discrepancy -- demonstrating that the tool presented in this paper is useful for meta-studies of astrophysical measurements from a large number of publications.

Fermi-LAT unidentified sources (unIDs) have proven to be compelling targets for performing indirect dark matter (DM) searches. In a previous work, we found that among the 1235 unIDs in Fermi-LAT catalogs (3FGL, 2FHL and 3FHL) only 44 of those are DM subhalos candidates. We now implement a spectral analysis to test whether these remaining sources are compatible or not with DM origin. This analysis is executed using almost 10 years of Pass 8 Fermi-LAT data. None of the unIDs are found to significantly prefer DM-induced emission compared to other, more conventional, astrophysical sources. In order to discriminate between pulsar and DM sources, we developed a new method which is based on the source's spectral curvature, peak energy, and its detection significance. We also look for spatial extension, which may be a hint for a DM origin according to our N-body simulation studies of the subhalo population. In addition, we used Gaia DR2 data to search for a potential stellar counterpart to our best DM subhalo candidates and, although no firm associations could be found, one of them coincides with the Sagittarius stream. Finally, previous constraints on the DM annihilation cross section are updated with the new number of remaining DM subhalo candidates among unIDs. Our limits now rule out canonical thermal WIMPs up to masses of 10 GeV for $b\bar{b}$ and 20 GeV for $\tau^+\tau^-$ annihilation channels, in this way being as sensitive and complementary to those obtained from other targets and probes.

In the last decade, the Kepler and CoRoT space-photometry missions have demonstrated the potential of asteroseismology as a novel, versatile and powerful tool to perform exquisite tests of stellar physics, and to enable precise and accurate characterisations of stellar properties, with impact on both exoplanetary and Galactic astrophysics. Based on our improved understanding of the strengths and limitations of such a tool, we argue for a new small/medium space mission dedicated to gathering high-precision, high-cadence, long photometric series in dense stellar fields. Such a mission will lead to breakthroughs in stellar astrophysics, especially in the metal poor regime, will elucidate the evolution and formation of open and globular clusters, and aid our understanding of the assembly history and chemodynamics of the Milky Way's bulge and few nearby dwarf galaxies.

A new all-sky visible and Near-InfraRed (NIR) space astrometry mission with a wavelength cutoff in the K-band is not just focused on a single or small number of key science cases. Instead, it is extremely broad, answering key science questions in nearly every branch of astronomy while also providing a dense and accurate visible-NIR reference frame needed for future astronomy facilities. For almost 2 billion common stars the combination of Gaia and a new all-sky NIR astrometry mission would provide much improved proper motions, answering key science questions -- from the solar system and stellar systems, including exoplanet systems, to compact galaxies, quasars, neutron stars, binaries and dark matter substructures. The addition of NIR will result in up to 8 billion newly measured stars in some of the most obscured parts of our Galaxy, and crucially reveal the very heart of the Galactic bulge region. In this white paper we argue that rather than improving on the accuracy, a greater overall science return can be achieved by going deeper than Gaia and by expanding the wavelength range to the NIR.

Gaia is a revolutionary space mission developed by ESA and is delivering 5 parameter astrometry, photometry and radial velocities over the whole sky with astrometric accuracies down to a few tens of micro-arcseconds. A weakness of Gaia is that it only operates at optical wavelengths. However, much of the Galactic centre and the spiral arm regions, important for certain studies, are obscured by interstellar extinction and this makes it difficult for Gaia to deeply probe. This problem can be overcome by switching to the Near Infra-Red (NIR) but this is not possible with silicon CCDs. Additionally, to scan the entire sky and make global absolute parallax measurements the spacecraft must have a constant rotation and this requires the detectors operate in Time Delayed Integration (TDI) mode or similar.

Using the cross-matched data of Gaia DR2 and 2MASS Point Source Catalog, we investigated the surface density distribution of stars aged ~1 Gyr in the thin disk in the range of 90{\deg} <= l <= 270{\deg}. We selected 4,654 stars above the turnoff corresponding to the age ~1 Gyr, that fall within a small box region in the color-magnitude diagram, (J-Ks)0 versus M(Ks), for which the distance and reddening are corrected. The selected sample shows an arm-like overdensity at 90{\deg} <= l <= 190{\deg}.This overdensity is located close to the Local arm traced by high-mass star forming regions (HMSFRs), but its pitch angle is slightly larger than that of the HMSFR-defined arm. Although the significance of the overdensity we report is marginal, its structure poses questions concerning both of the competing scenarios of spiral arms, the density-wave theory and the dynamic spiral arm model. The offset between the arms traced by stars and HMSFRs, i.e., gas, is difficult to be explained by the dynamic arm scenario. On the other hand, the pitch angle of the stellar Local arm, if confirmed, larger than that of the Perseus arm is difficult to be explained with the classical density-wave scenario. The dynamic arm scenario can explain it if the Local arm is in a growing up phase, while the Perseus arm is in a disrupting phase. Our result provides a new and complex picture of the Galactic spiral arms, and encourages further studies.

Gaia DR2 has provided an unprecedented wealth of information about the kinematics of stars in the Solar neighbourhood, and has highlighted the degree of features in the Galactic disc. We confront the data with a range of bar and spiral models in both action-angle space, and the $R_{\mathrm{G}}-v_{\phi}$ plane. We find that the phase mixing induced by transient spiral structure creates ridges and arches in the local kinematics which are consistent with the Gaia data. We are able to produce a qualitatively good match to the data when combined with a bar with a variety of pattern speeds, and show that it is non trivial to decouple the effects of the bar and the spiral structure.

We recommend configuring the LSST coverage of the inner Galactic plane to allow the production of the definitive age/metallicity map of the Galactic bulge from LSST data, matched to external surveys where appropriate. This will allow the formation history of the Galactic bulge to be reconstructed, as well as furnishing a huge legacy dataset to support one of the key LSST science goals (Mapping the Milky Way). We recommend precise multi-color $grizy$~photometry as deeply as spatial crowding will allow (ideally completing early in the project), complemented by a single-filter survey spread over the entire 10-year time baseline. Both strands should cover as broad an area within the Bulge as possible, with image quality sufficient to reach at least the bulge main sequence turn-off in seeing-limited observations. We specify metrics and a figure of merit by which candidate observing strategies could be evaluated with respect to `static' bulge science (proper motions and photometry)

The intracluster medium of galaxy clusters is permeated by {\mu}G magnetic fields. Observations with current and future facilities have the potential to illuminate the role of these magnetic fields play in the astrophysical processes of galaxy clusters. To obtain a greater understanding of how the initial seed fields evolve to the magnetic fields in the intracluster medium requires magnetohydrodynamic simulations. We critically assess the current Smoothed Particle Magneto-Hydrodynamics (SPMHD) schemes, especially highlighting the impact of a hyperbolic divergence cleaning scheme and artificial resistivity switch on the magnetic field evolution in cosmological simulations of the formation of a galaxy cluster using the N-body/SPMHD code gcmhd++. The impact and performance of the cleaning scheme and two different schemes for the artificial resistivity switch is demonstrated via idealized test cases and cosmological simulations. We demonstrate that the hyperbolic divergence cleaning scheme is effective at suppressing the growth of the numerical divergence error of the magnetic field and should be applied to any SPMHD simulation. Although the artificial resistivity is important in the strong field regime, it can suppress the growth of the magnetic field in the weak field regime, such as galaxy clusters. With sufficient resolution, simulations with divergence cleaning can reproduce observed magnetic fields. We conclude that the cleaning scheme alone is sufficient for galaxy cluster simulations, but our results indicate that the SPMHD scheme must be carefully chosen depending on the regime of the magnetic field.

We apply a simple axisymmetric disc model to 218 Galactic Cepheids whose accurate measurements of the distance and velocities are obtained by cross-matching an existing Cepheids catalogue with the Gaia DR2 data. Our model fit determines the "local centrifugal speed", $V_\mathrm{c}$ $-$ defined as the rotation speed required to balance the local radial gravitational force $-$ at the Sun's location to be $V_{c}(R_0)=236\pm 3$ km s$^{-1}$ and the Sun's azimuthal and radial peculiar motions to be $V_{\odot}=12.4\pm0.7$ km s$^{-1}$ and $U_{\odot}=7.7\pm0.9$ km s$^{-1}$, respectively. These results are obtained with strong priors on the solar radius, $R_0=8.2\pm0.1$ kpc, and Sun's angular rotation velocity, $\Omega_{\odot}=30.24\pm0.12$ km s$^{-1}$ kpc$^{-1}$. We also applied the axisymmetric model to mock data from an N-body/hydrodynamic simulation of a Milky Way-like galaxy with a bar and spiral arms. We find that our axisymmetric model fit to the young stars recovers the local centrifugal speed reasonably well, even in the face of significant non-axisymmetry. However, the local centrifugal speed determined from our Cepheid sample could suffer from systematic uncertainty as large as 6 km s$^{-1}$.

The second data release from ESA's Gaia mission has revealed many ridge-like structures in the velocity distribution of the Milky Way. We show that these can arise naturally from winding transient spiral structure that is commonly seen in N-body simulations of disk galaxies. We construct test particle models of the winding spiral structure, and compare the resulting distribution of orbits with the observed two-dimensional velocity distribution in the extended solar neighbourhood and with the distribution of rotational velocities over 8 kpc along the Sun--Galactic-centre--Galactic anti-centre line. We show that the ridges in these observations are well reproduced by the winding spiral model. Additionally, we demonstrate that the transient winding spiral potential can create a Hercules-like feature in the kinematics of the solar neighbourhood, either alone, or in combination with a long-slow bar potential.

We make a first attempt to find dwarf galaxies in eight \Fermi-LAT extended, unassociated, source fields using \Gaia\ DR2. We probe previously unexplored heliocentric distances of $d<20$~kpc with an extreme-deconvolution (XD) technique. We find no signature of a dwarf galaxy in any of these fields despite \Gaia's excellent astrometric accuracy. We estimate our detection limits by applying the XD method to mock data, obtaining a conservative limit on the stellar mass of $M_* < 10^4$~M$_{\sun}$ for $d < 20$\, kpc. Such a low stellar mass implies either a low-mass subhalo, or a massive stripped-down subhalo. We use an analytic model for stripped subhalos to argue that, given the sizes and fluxes of the \Fermi-LAT sources, we can reject the hypothesis that they owe to dark matter annihilation.

By taking advantage of the superb measurements of position and velocity for an unprecedented large number of stars provided in Gaia DR2, we have generated the first maps of the rotation velocity, $V_{\rm rot}$, and vertical velocity, $V_{\rm z}$, distributions as a function of the Galactocentric radius, $R_{\rm gal}$, across a radial range of $5<R_{\rm gal}<12$~kpc. In the $R-V_{\rm rot}$ map, we have identified many diagonal ridge features, which are compared with the location of the spiral arms and the expected outer Lindblad resonance of the Galactic bar. We have detected also radial wave-like oscillations of the peak of the vertical velocity distribution.

We have discovered a clear sign of the disruption phase of the Perseus arm in the Milky Way using Cepheid variables, taking advantage of the accurately measured distances of Cepheids and the proper motions from Gaia Data Release 1. Both the Galactocentric radial and rotation velocities of 77 Cepheids within 1.5 kpc of the Perseus arm are correlated with their distances from the locus of the Perseus arm, as the trailing side is rotating faster and moving inward compared to the leading side. We also found a negative vertex deviation for the Cepheids on the trailing side, $-27.6\pm2.4$ deg, in contrast to the positive vertex deviation in the solar neighborhood. This is, to our knowledge, the first direct evidence that the vertex deviation around the Perseus arm is affected by the spiral arm. We compared these observational trends with our $N$-body/hydrodynamics simulations based on a static density-wave spiral scenario and those based on a transient dynamic spiral scenario. Although our comparisons are limited to qualitative trends, they strongly favor a conclusion that the Perseus arm is in the disruption phase of a transient arm.

Investigations of the origin and evolution of the Milky Way disk have long relied on chemical and kinematic identification of its components to reconstruct our Galactic past. Difficulties in determining precise stellar ages have restricted most studies to small samples, normally confined to the solar neighbourhood. Here we break this impasse with the help of asteroseismic inference and perform a chronology of the evolution of the disk throughout the age of the Galaxy. We chemically dissect the Milky Way disk population using a sample of red giant stars spanning out to 2~kpc in the solar annulus observed by the {\it Kepler} satellite, with the added dimension of asteroseismic ages. Our results reveal a clear difference in age between the low- and high-$\alpha$ populations, which also show distinct velocity dispersions in the $V$ and $W$ components. We find no tight correlation between age and metallicity nor [$\alpha$/Fe] for the high-$\alpha$ disk stars. Our results indicate that this component formed over a period of more than 2~Gyr with a wide range of [M/H] and [$\alpha$/Fe] independent of time. Our findings show that the kinematic properties of young $\alpha$-rich stars are consistent with the rest of the high-$\alpha$ population and different from the low-$\alpha$ stars of similar age, rendering support to their origin being old stars that went through a mass transfer or stellar merger event, making them appear younger, instead of migration of truly young stars formed close to the Galactic bar.

We investigate the vertical metallicity gradients of five mono-age stellar populations between 0 and 11 Gyr for a sample of 18 435 dwarf stars selected from the cross-matched Tycho-Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiment (RAVE) Data Release 5. We find a correlation between the vertical metallicity gradients and age, with no vertical metallicity gradient in the youngest population and an increasingly steeper negative vertical metallicity gradient for the older stellar populations. The metallicity at disc plane remains almost constant between 2 and 8 Gyr, and it becomes significantly lower for the $8 < \tau \leqslant 11$ Gyr population. The current analysis also reveals that the intrinsic dispersion in metallicity increases steadily with age. We discuss that our results are consistent with a scenario that (thin) disc stars formed from a flaring (thin) star-forming disc.

Among dynamical modeling techniques, the made-to-measure (M2M) method for modeling steady-state systems is among the most flexible, allowing non-parametric distribution functions in complex gravitational potentials to be modeled efficiently using N-body particles. Here we propose and test various improvements to the standard M2M method for modeling observed data, illustrated using the simple setup of a one-dimensional harmonic oscillator. We demonstrate that nuisance parameters describing the modeled system's orientation with respect to the observer---e.g., an external galaxy's inclination or the Sun's position in the Milky Way---as well as the parameters of an external gravitational field can be optimized simultaneously with the particle weights. We develop a method for sampling from the high-dimensional uncertainty distribution of the particle weights. We combine this in a Gibbs sampler with samplers for the nuisance and potential parameters to explore the uncertainty distribution of the full set of parameters. We illustrate our M2M improvements by modeling the vertical density and kinematics of F-type stars in Gaia DR1. The novel M2M method proposed here allows full probabilistic modeling of steady-state dynamical systems, allowing uncertainties on the non-parametric distribution function and on nuisance parameters to be taken into account when constraining the dark and baryonic masses of stellar systems.

The stellar disk of the Milky Way shows complex spatial and abundance structure that is central to understanding the key physical mechanisms responsible for shaping our Galaxy. In this study, we use six very high resolution cosmological zoom simulations of Milky Way-sized haloes to study the prevalence and formation of chemically distinct disc components. We find that our simulations develop a clearly bimodal distribution in the $[\rm \alpha/Fe]$ -- $[\rm Fe/H]$ plane. We find two main pathways to creating this dichotomy which operate in different regions of the galaxies: a) an early ($z>1$) and intense high-$\rm[\alpha/Fe]$ star formation phase in the inner region ($R\lesssim 5$ kpc) induced by gas-rich mergers, followed by more quiescent low-$\rm[\alpha/Fe]$ star formation; and b) an early phase of high-$\rm[\alpha/Fe]$ star formation in the outer disc followed by a shrinking of the gas disc owing to a temporarily lowered gas accretion rate, after which disc growth resumes. In process b), a double-peaked star formation history around the time and radius of disc shrinking accentuates the dichotomy. If the early star formation phase is prolonged (rather than short and intense), chemical evolution proceeds as per process a) in the inner region, but the dichotomy is less clear. In the outer region, the dichotomy is only evident if the first intense phase of star formation covers a large enough radial range before disc shrinking occurs; otherwise, the outer disc consists of only low-$\rm[\alpha/Fe]$ sequence stars. We discuss the implication that both processes occurred in the Milky Way.

We apply the GC3 stream-finding method to RR Lyrae stars (RRLS) in the Catalina survey. We find two RRLS stream candidates at $>4\sigma$ confidence and another 12 at $>3.5\sigma$ confidence over the Galactocentric distance range $4 < D/{\rm kpc} < 26$. Of these, only two are associated with known globular clusters (NGC 1261 and Arp2). The remainder are candidate `orphan' streams, consistent with the idea that globular cluster streams are most visible close to dissolution. Our detections are likely a lower bound on the total number of dissolving globulars in the inner galaxy, since many globulars have few RRLS while only the brightest streams are visible over the Galactic RRLS background, particularly given the current lack of kinematical information. We make all of our candidate streams publicly available and provide a new GALSTREAMS Python library for the footprints of all known streams and overdensities in the Milky Way.

We have developed a novel Markov Chain Mote Carlo (MCMC) chemical "painting" technique to explore possible radial and vertical metallicity gradients for the thick disc progenitor. In our analysis we match an N-body simulation to the data from the Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. We assume that the thick disc has a constant scale-height and has completed its formation at an early epoch, after which time radial mixing of its stars has taken place. Under these assumptions, we find that the initial radial metallicity gradient of the thick disc progenitor should not be negative, but either flat or even positive, to explain the current negative vertical metallicity gradient of the thick disc. Our study suggests that the thick disc was built-up in an inside-out and upside-down fashion, and older, smaller and thicker populations are more metal poor. In this case, star forming discs at different epochs of the thick disc formation are allowed to have different radial metallicity gradients, including a negative one, which helps to explain a variety of slopes observed in high redshift disc galaxies. This scenario helps to explain the positive slope of the metallicity-rotation velocity relation observed for the Galactic thick disc. On the other hand, radial mixing flattens the slope of an existing gradient.

Using N-body simulations of the Galactic disks, we qualitatively study how the metallicity distributions of the thick and thin disk stars are modified by radial mixing induced by the bar and spiral arms. We show that radial mixing drives a positive vertical metallicity gradient in the mono-age disk population whose initial scale-height is constant and initial radial metallicity gradient is tight and negative. On the other hand, if the initial disk is flaring, with scale-height increasing with galactocentric radius, radial mixing leads to a negative vertical metallicity gradient, which is consistent with the current observed trend. We also discuss impacts of radial mixing on the metallicity distribution of the thick disk stars. By matching the metallicity distribution of N-body models to the SDSS/APOGEE data, we argue that the progenitor of the Milky Way's thick disk should not have a steep negative metallicity gradient.

Deciphering the assembly history of the Milky Way is a formidable task, which becomes possible only if one can produce high-resolution chrono-chemo-kinematical maps of the Galaxy. Data from large-scale astrometric and spectroscopic surveys will soon provide us with a well-defined view of the current chemo-kinematical structure of the Milky Way, but will only enable a blurred view on the temporal sequence that led to the present-day Galaxy. As demonstrated by the (ongoing) exploitation of data from the pioneering photometric missions CoRoT, Kepler, and K2, asteroseismology provides the way forward: solar-like oscillating giants are excellent evolutionary clocks thanks to the availability of seismic constraints on their mass and to the tight age-initial-mass relation they adhere to. In this paper we identify five key outstanding questions relating to the formation and evolution of the Milky Way that will need precise and accurate ages for large samples of stars to be addressed, and we identify the requirements in terms of number of targets and the precision on the stellar properties that are needed to tackle such questions. By quantifying the asteroseismic yields expected from PLATO for red-giant stars, we demonstrate that these requirements are within the capabilities of the current instrument design, provided that observations are sufficiently long to identify the evolutionary state and allow robust and precise determination of acoustic-mode frequencies. This will allow us to harvest data of sufficient quality to reach a 10% precision in age. This is a fundamental pre-requisite to then reach the more ambitious goal of a similar level of accuracy, which will only be possible if we have to hand a careful appraisal of systematic uncertainties on age deriving from our limited understanding of stellar physics, a goal which conveniently falls within the main aims of PLATO's core science.

We present a numerical study of the impact of AGN accretion and feedback on the star formation history of barred disc galaxies. Our goal is to determine whether the effect of feedback is positive (enhanced star formation) or negative (quenched star formation), and to what extent. We performed a series of 12 hydrodynamical simulations of disc galaxies, 10 barred and 2 unbarred, with various initial gas fractions and AGN feedback prescriptions. In barred galaxies, gas is driven toward the centre of the galaxy and causes a starburst, followed by a slow decay, while in unbarred galaxies the SFR increases slowly and steadily. AGN feedback suppresses star formation near the central black hole. Gas is pushed away from the black hole, and collides head-on with inflowing gas, forming a dense ring at a finite radius where star formation is enhanced. We conclude that both negative and positive feedback are present, and these effects mostly cancel out. There is no net quenching or enhancement in star formation, but rather a displacement of the star formation sites to larger radii. In unbarred galaxies, where the density of the central gas is lower, quenching of star formation near the black hole is more efficient, and enhancement of star formation at larger radii is less efficient. As a result, negative feedback dominates. Lowering the gas fraction reduces the star formation rate at all radii, whether or not there is a bar or an AGN.

Here we investigate the evolution of a Milky Way (MW) -like galaxy with the aim of predicting the properties of its progenitors all the way from $z \sim 20$ to $z = 0$. We apply GAMESH (Graziani et al. 2015) to a high resolution N-Body simulation following the formation of a MW-type halo and we investigate its properties at $z \sim 0$ and its progenitors in $0 < z < 4$. Our model predicts the observed galaxy main sequence, the mass-metallicity and the fundamental plane of metallicity relations in $0 < z < 4$. It also reproduces the stellar mass evolution of candidate MW progenitors in $0 \lesssim z \lesssim 2.5$, although the star formation rate and gas fraction of the simulated galaxies follow a shallower redshift dependence. We find that while the MW star formation and chemical enrichment are dominated by the contribution of galaxies hosted in Lyman $\alpha$-cooling halos, at z > 6 the contribution of star forming mini-halos is comparable to the star formation rate along the MW merger tree. These systems might then provide an important contribution in the early phases of reionization. A large number of mini-halos with old stellar populations, possibly Population~III stars, are dragged into the MW or survive in the Local Group. At low redshift dynamical effects, such as halo mergers, tidal stripping and halo disruption redistribute the baryonic properties among halo families. These results are critically discussed in light of future improvements including a more sophisticated treatment of radiative feedback and inhomogeneous metal enrichment.

We report on the detection of a small overdensity of stars in velocity space with systematically higher Galactocentric rotation velocity than the Sun by about 20 km s$^{-1}$ in the $Gaia$ Data Release 1 Tycho-Gaia astrometric solution (TGAS) data. We find these fast Galactic rotators more clearly outside of the Solar radius, compared to inside of the Solar radius. In addition, the velocity of the fast Galactic rotators is independent of the Galactocentric distance up to $R-R_{\odot}\sim0.6$ kpc. Comparing with numerical models, we qualitatively discuss that a possible cause of this feature is the co-rotation resonance of the Perseus spiral arm, where the stars in peri-centre phase in the trailing side of the Perseus spiral arm experience an extended period of acceleration owing to the torque from the Perseus arm.

The length and pattern speed of the Milky Way bar are still controversial. Photometric and spectroscopic surveys of the inner Galaxy, as well as gas kinematics, favour a long and slowly rotating bar, with corotation around a Galactocentric radius of 6 kpc. On the other hand, the existence of the Hercules stream in local velocity space favours a short and fast bar with corotation around 4 kpc. This follows from the fact that the Hercules stream looks like a typical signature of the outer Lindblad resonance of the bar. As we showed recently, reconciling this local stream with a slow bar would need to find a yet unknown alternative explanation, based for instance on the effect of spiral arms. Here, by combining the TGAS catalogue of the Gaia DR1 with LAMOST radial velocities, we show that the position of Hercules in velocity space as a function of radius in the outer Galaxy indeed varies exactly as predicted by fast bar models with a pattern speed no less than 1.8 times the circular frequency at the Sun's position.

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 [$\alpha$/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.

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.

Understanding the nature of spiral structure in disk galaxies is one of the main, and still unsolved questions in galactic astronomy. However, theoretical works are proposing new testable predictions whose detection is becoming feasible with recent development in instrumentation. In particular, streaming motions along spiral arms are expected to induce azimuthal variations in the chemical composition of a galaxy at a given galactic radius. In this letter we analyse the gas content in NGC 6754 with VLT/MUSE data to characterise its 2D chemical composition and H$\alpha$ line-of-sight velocity distribution. We find that the trailing (leading) edge of the NGC 6754 spiral arms show signatures of tangentially-slower, radially-outward (tangentially-faster, radially-inward) streaming motions of metal-rich (poor) gas over a large range of radii. These results show direct evidence of gas radial migration for the first time. We compare our results with the gas behaviour in a $N$-body disk simulation showing spiral morphological features rotating with a similar speed as the gas at every radius, in good agreement with the observed trend. This indicates that the spiral arm features in NGC 6754 may be transient and rotate similarly as the gas does at a large range of radii.

Using idealised N-body simulations of a Milky Way-sized disc galaxy, we qualitatively study how the metallicity distributions of the thin disc star particles are modified by the formation of the bar and spiral arm structures. The thin disc in our numerical experiments initially has a tight negative radial metallicity gradient and a constant vertical scale-height. We show that the radial mixing of stars drives a positive vertical metallicity gradient in the thin disc. On the other hand, if the initial thin disc is flared, with vertical scale-height increasing with galactocentric radius, the metal poor stars originally in the outer disc become dominant in regions above the disc plane at every radii. This process can drive a negative vertical metallicity gradient, which is consistent with the current observed trend. This model mimics a scenario where the star-forming thin disc was flared in the outer region at earlier epochs. Our numerical experiment with an initial flared disc predicts that the negative vertical metallicity gradient of the mono-age relatively young thin disc population should be steeper in the inner disc, and the radial metallicity gradient of the mono-age population should be shallower at greater heights above the disc plane. We also predict that the metallicity distribution function of mono-age young thin disc populations above the disc plane would be more positively skewed in the inner disc compared to the outer disc.

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.

This article is based on our discussion session on Milky Way models at the 592 WE-Heraeus Seminar, Reconstructing the Milky Way's History: Spectroscopic Surveys, Asteroseismology and Chemodynamical models. The discussion focused on the following question: "Are there distinct thick and thin disks?". The answer to this question depends on the definition one adopts for thin and thick disks. The participants of this discussion converged to the idea that there are at least two different types of disks in the Milky Way. However, there are still important open questions on how to best define these two types of disks (chemically, kinematically, geometrically or by age?). The question of what is the origin of the distinct disks remains open. The future Galactic surveys which are highlighted in this conference should help us answering these questions. The almost one-hour debate involving researchers in the field representing different modelling approaches (Galactic models such as TRILEGAL, Besancon and Galaxia, chemical evolution models, extended distribution functions method, chemodynamics in the cosmological context, and self-consistent cosmological simulations) illustrated how important is to have all these parallel approaches. All approaches have their advantages and shortcomings (also discussed), and different approaches are useful to address specific points that might help us answering the more general question above.

We present a new theoretical population synthesis model (the Galaxy Model) to examine and deal with large amounts of data from surveys of the Milky Way and to decipher the present and past structure and history of our own Galaxy. We assume the Galaxy to consist of a superposition of many composite stellar populations belonging to the thin and thick disks, the stellar halo and the bulge, and to be surrounded by a single dark matter halo component. A global model for the Milky Way's gravitational potential is built up self-consistently with the density profiles from the Poisson equation. In turn, these density profiles are used to generate synthetic probability distribution functions (PDFs) for the distribution of stars in colour-magnitude diagrams (CMDs). Finally, the gravitational potential is used to constrain the stellar kinematics by means of the moment method on a (perturbed)-distribution function. Spiral arms perturb the axisymmetric disk distribution functions in the linear response framework of density-wave theory where we present an analytical formula of the so-called `reduction factor' using Hypergeometric functions. Finally, we consider an analytical non-axisymmetric model of extinction and an algorithm based on the concept of probability distribution function to handle colour magnitude diagrams with a large number of stars. A genetic algorithm is presented to investigate both the photometric and kinematic parameter space. This galaxy model represents the natural framework to reconstruct the structure of the Milky Way from the heterogeneous data set of surveys such as Gaia-ESO, SEGUE, APOGEE2, RAVE and the Gaia mission.

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.

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.

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.

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.

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.

We have developed our original made-to-measure (M2M) algorithm, PRIMAL, with the aim of modelling the Galactic disc from upcoming Gaia data. From a Milky Way like N-body disc galaxy simulation, we have created mock Gaia data using M0III stars as tracers, taking into account extinction and the expected Gaia errors. In PRIMAL, observables calculated from the N-body model are compared with the target stars, at the position of the target stars. Using PRIMAL, the masses of the N-body model particles are changed to reproduce the target mock data, and the gravitational potential is automatically adjusted by the changing mass of the model particles. We have also adopted a new resampling scheme for the model particles to keep the mass resolution of the N-body model relatively constant. We have applied PRIMAL to this mock Gaia data and we show that PRIMAL can recover the structure and kinematics of a Milky Way like barred spiral disc, along with the apparent bar structure and pattern speed of the bar despite the galactic extinction and the observational errors.

Here we introduce GAMESH, a novel pipeline which implements self-consistent radiative and chemical feedback in a computational model of galaxy formation. By combining the cosmological chemical-evolution model GAMETE with the radiative transfer code CRASH, GAMESH can post process realistic outputs of a N-body simulation describing the redshift evolution of the forming galaxy. After introducing the GAMESH implementation and its features, we apply the code to a low-resolution N-body simulation of the Milky Way formation and we investigate the combined effects of self-consistent radiative and chemical feedback. Many physical properties, which can be directly compared with observations in the Galaxy and its surrounding satellites, are predicted by the code along the merger-tree assembly. The resulting redshift evolution of the Local Group star formation rates, reionisation and metal enrichment along with the predicted Metallicity Distribution Function of halo stars are critically compared with observations. We discuss the merits and limitations of the first release of GAMESH, also opening new directions to a full implementation of feedback processes in galaxy formation models by combining semi-analytic and numerical methods.

We compare the stellar motion around a spiral arm created in two different scenarios, transient/co-rotating spiral arms and density-wave-like spiral arms. We generate Gaia mock data from snapshots of the simulations following these two scenarios using our stellar population code, SNAPDRAGONS, which takes into account dust extinction and the expected Gaia errors. We compare the observed rotation velocity around a spiral arm similar in position to the Perseus arm, and find that there is a clear difference in the velocity features around the spiral arm between the co-rotating spiral arm and the density-wave-like spiral arm. Our result demonstrates that the volume and accuracy of the Gaia data are sufficient to clearly distinguish these two scenarios of the spiral arms.

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.

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.

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.

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.

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.

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]

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.

Chronos is our response to ESA's call for white papers to define the science for the future L2, L3 missions. Chronos targets the formation and evolution of galaxies, by collecting the deepest NIR spectroscopic data, from the formation of the first galaxies at z~10 to the peak of formation activity at z~1-3. The strong emission from the atmospheric background makes this type of survey impossible from a ground-based observatory. The spectra of galaxies represent the equivalent of a DNA fingerprint, containing information about the past history of star formation and chemical enrichment. The proposed survey will allow us to dissect the formation process of galaxies including the timescales of quenching triggered by star formation or AGN activity, the effect of environment, the role of infall/outflow processes, or the connection between the galaxies and their underlying dark matter haloes. To provide these data, the mission requires a 2.5m space telescope optimised for a campaign of very deep NIR spectroscopy. A combination of a high multiplex and very long integration times will result in the deepest, largest, high-quality spectroscopic dataset of galaxies from z=1 to 12, spanning the history of the Universe, from 400 million to 6 billion years after the big bang, i.e. covering the most active half of cosmic history.

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.

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.

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.

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.

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.

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.

Aims. We present here a new theoretical approach to population synthesis. The aim is to predict colour magnitude diagrams (CMDs) for huge numbers of stars. With this method we generate synthetic CMDs for N-body simulations of galaxies. Sophisticated hydrodynamic N-body models of galaxies require equal quality simulations of the photometric properties of their stellar content. The only prerequisite for the method to work is very little information on the star formation and chemical enrichment histories, i.e. the age and metallicity of all star-particles as a function of time. The method takes into account the gap between the mass of real stars and that of the star-particles in N-body simulations, which best correspond to the mass of star clusters with different age and metallicity, i.e. a manifold of single stellar sopulations (SSP). Methods. The theory extends the concept of SSP to include the phase-space (position and velocity) of each star. Furthermore, it accelerates the building up of simulated CMD by using a database of theoretical SSPs that extends to all ages and metallicities of interest. Finally, it uses the concept of distribution functions to build up the CMD. The technique is independent of the mass resolution and the way the N-body simulation has been calculated. This allows us to generate CMDs for simulated stellar systems of any kind: from open clusters to globular clusters, dwarf galaxies, or spiral and elliptical galaxies. Results. The new theory is applied to an N-body simulation of a disc galaxy to test its performance and highlight its flexibility.

Within a cosmological hydrodynamical simulation, we form a disc galaxy with sub- components which can be assigned to a thin stellar disc, thick disk, and a low mass stellar halo via a chemical decomposition. The thin and thick disc populations so selected are distinct in their ages, kinematics, and metallicities. Thin disc stars are young (<6.6 Gyr), possess low velocity dispersion ({\sigma}U,V,W = 41, 31, 25 km/s), high [Fe/H], and low [O/Fe]. The thick disc stars are old (6.6<age<9.8 Gyrs), lag the thin disc by \sim21 km/s, possess higher velocity dispersion ({\sigma}U,V,W = 49, 44, 35 km/s), relatively low [Fe/H] and high [O/Fe]. The halo component comprises less than 4% of stars in the "solar annulus" of the simulation, has low metallicity, a velocity ellipsoid defined by ({\sigma}U,V,W = 62, 46, 45 km/s) and is formed primarily in-situ during an early merger epoch. Gas-rich mergers during this epoch play a major role in fuelling the formation of the old disc stars (the thick disc). This is consistent with studies which show that cold accretion is the main source of a disc galaxy's baryons. Our simulation initially forms a relatively short (scalelength \sim1.7 kpc at z=1) and kinematically hot disc, primarily from gas accreted during the galaxy's merger epoch. Far from being a competing formation scenario, migration is crucial for reconciling the short, hot, discs which form at high redshift in {\Lambda}CDM, with the properties of the thick disc at z=0. The thick disc, as defined by its abundances maintains its relatively short scale-length at z = 0 (2.31 kpc) compared with the total disc scale-length of 2.73 kpc. The inside-out nature of disc growth is imprinted the evolution of abundances such that the metal poor {\alpha}-young population has a larger scale-length (4.07 kpc) than the more chemically evolved metal rich {\alpha}-young population (2.74 kpc).

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, [{\alpha}/Fe]-[Fe/H] and Type Ia/II supernovae rates.

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.

We examine radial and vertical metallicity gradients using a suite of disk galaxy simulations, supplemented with two classic chemical evolution approaches. We determine the rate of change of gradient and reconcile differences between extant models and observations within the `inside-out' disk growth paradigm. A sample of 25 disks is used, consisting of 19 from our RaDES (Ramses Disk Environment Study) sample, realised with the adaptive mesh refinement code RAMSES. Four disks are selected from the MUGS (McMaster Unbiased Galaxy Simulations) sample, generated with the smoothed particle hydrodynamics (SPH) code GASOLINE, alongside disks from Rahimi et al. (GCD+) and Kobayashi & Nakasato (GRAPE-SPH). Two chemical evolution models of inside-out disk growth were employed to contrast the temporal evolution of their radial gradients with those of the simulations. We find that systematic differences exist between the predicted evolution of radial abundance gradients in the RaDES and chemical evolution models, compared with the MUGS sample; specifically, the MUGS simulations are systematically steeper at high-redshift, and present much more rapid evolution in their gradients. We find that the majority of the models predict radial gradients today which are consistent with those observed in late-type disks, but they evolve to this self-similarity in different fashions, despite each adhering to classical `inside-out' growth. We find that radial dependence of the efficiency with which stars form as a function of time drives the differences seen in the gradients; systematic differences in the sub-grid physics between the various codes are responsible for setting these gradients. Recent, albeit limited, data at redshift z=1.5 are consistent with the steeper gradients seen in our SPH sample, suggesting a modest revision of the classical chemical evolution models may be required.

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.

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.

Radio observations of galaxy clusters show that the intra cluster medium is permeated by \mu 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.

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.

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.

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.

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.

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.

We study the age and metallicity distribution function (MDF) of metal-poor stars in the Milky Way halo as a function of galactocentric radius by combining N-body simulations and semi-analytical methods. We find that the oldest stars populate the innermost region, while extremely metal-poor stars are more concentrated within r < 60 kpc. The MDF of [Fe/H] < -2 stars varies only very weakly within the central 50 kpc, while the relative contribution of [Fe/H] < -2 stars strongly increases with r, varying from 16% within 7 kpc < r < 20 kpc up to > 40% for r > 20 kpc. This is due to the faster descent of the spatial distribution (as seen from Earth) of the more enriched population. This implies that the outer halo < 40 kpc is the best region to search for very metal-poor stars. Beyond ~ 60 kpc the density of [Fe/H] < -2 stars is maximum within dwarf galaxies. All these features are imprinted by a combination of (i) the virialization epoch of the star-forming haloes, and (ii) the metal enrichment history of the Milky Way environment.

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.

We analyse the kinematics and chemistry of the bulge stars of two simulated disc galaxies using our chemodynamical galaxy evolution code GCD+. First we compare stars that are born inside the galaxy with those that are born outside the galaxy and are accreted into the centre of the galaxy. Stars that originate outside of the bulge are accreted into it early in its formation within 3 Gyrs so that these stars have high [alpha/Fe] as well as having a high total energy reflecting their accretion to the centre of the galaxy. Therefore, higher total energy is a good indicator for finding accreted stars. The bulges of the simulated galaxies formed through multiple mergers separated by about a Gyr. Since [alpha/Fe] is sensitive to the first few Gyrs of star formation history, stars that formed during mergers at different epochs show different [alpha/Fe]. We show that the [Mg/Fe] against star formation time relation can be very useful to identify a multiple merger bulge formation scenario, provided there is sufficiently good age information available. Our simulations also show that stars formed during one of the merger events retain a systematically prograde rotation at the final time. This demonstrates that the orbit of the ancient merger that helped to form the bulge could still remain in the kinematics of bulge stars.

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.

Using high-resolution N-body simulations, we examine whether a major dry merger mitigates the difference in the radial density distributions between red and blue globular clusters (GCs). To this end, we study the relation between the density slope of the GCs in merger progenitors and that in a merger remnant, when the density distribution is described by $n_{\rm GC}\propto r^{-\alpha}$. We also study how our results depend on the merger orbit and the size of the core radius of the initial GC density distribution. We find that a major dry merger makes the GC profile flatter, and the steeper initial GC profile leads to more significant flattening, especially if the initial slope is steeper than $\alpha\sim3.5$. Our result suggests that if there is a major dry merger of elliptical galaxies whose red GCs have a steeper radial profile than the blue GCs, as currently observed, and their slopes are steeper than $\alpha\sim3.5$, the difference in the slopes between two populations becomes smaller after dry mergers. Therefore, the observed slopes of red and blue GCs can be a diagnostic of the importance of dry merger. The current observational data show that the red and blue GCs have more comparable and shallower slopes in some luminous galaxies, which may indicate that they have experienced dry mergers.

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.

We summarize the high-resolution science that has been done on high redshift galaxies with Adaptive Optics (AO) on the world's largest ground-based facilities and with the Hubble Space Telescope (HST). These facilities complement each other. Ground-based AO provides better light gathering power and in principle better resolution than HST, giving it the edge in high spatial resolution imaging and high resolution spectroscopy. HST produces higher quality, more stable PSF's over larger field-of-view's in a much darker sky-background than ground-based AO, and yields deeper wide-field images and low-resolution spectra than the ground. Faint galaxies have steadily decreasing sizes at fainter fluxes and higher redshifts, reflecting the hierarchical formation of galaxies over cosmic time. HST has imaged this process in great structural detail to z<~6, and ground-based AO and spectroscopy has provided measurements of their masses and other physical properties with cosmic time. Last, we review how the 6.5 meter James Webb Space Telescope (JWST) will measure First Light, reionization, and galaxy assembly in the near--mid-IR after 2013.

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)

We use a cosmological chemodynamical simulation to study how the group environment impacts the star formation properties of disk galaxies. The simulated group has a total mass of M~8x10^12 Msun and a total X-ray luminosity of L_X~10^41 erg s^-1. Our simulation suggests that ram pressure is not sufficient in this group to remove the cold disk gas from a V_rot~150 km s^-1 galaxy. However, the majority of the hot gas in the galaxy is stripped over a timescale of approximately 1 Gyr. Since the cooling of the hot gas component provides a source for new cold gas, the stripping of the hot component effectively cuts off the supply of cold gas. This in turn leads to a quenching of star formation. The galaxy maintains the disk component after the cold gas is consumed, which may lead to a galaxy similar to an S0. Our self-consistent simulation suggests that this strangulation mechanism works even in low mass groups, providing an explanation for the lower star formation rates in group galaxies relative to galaxies in the field.

A significant fraction of nearby elliptical galaxies are known to have high density gas disks in their circumnuclear (CN) region (0.1 to a few kpc). Yet, ellipticals, especially luminous ones, show little signs of recent star formation (SF). To investigate the possible cause of the dearth of SF in these systems, we study the gravitational stability of CN gas disks embedded within the potentials of both the stellar bulge and the central massive black hole (BH) in ellipticals. We find that CN disks in higher mass galaxies are generally more stable than those in lower mass galaxies, because higher mass galaxies tend to have more massive BHs and more centrally concentrated stellar density profiles. We also consider the case in which the central stellar density profile has a core, which is often observed for ellipticals whose total stellar mass is higher than about 10^11 Msun. Such a cored stellar density profile leads to more unstable CN disks than the power-law density profile characteristic of less massive galaxies. However, the more massive BHs in high-mass galaxies act to stabilize the CN disk. Our results demonstrate that the gravitational potentials of both the central BH and the stellar component should be taken into account when studying the properties of CN disks, as their stability is sensitive to both the BH mass and the stellar density profile. Our results could explain the observed trend that less luminous ellipticals have a greater tendency to exhibit ongoing SF than giant ellipticals.

The difference in density profiles of the contributions from different density peaks to dark matter halos results in certain expectations about the Milky Way's stellar halo. We cut our simulated halo stars into two populations: those forming before/during the last major merger, and those accreted after the last major merger. The former population are more centrally located at z=0, while stars forming in low mass late forming proto-galaxies are spread through the halo. A difference in observed binding energy distinguishes these two populations. We look at possible chemical abundance signatures of the two populations. We also show that galaxies forming in isolated low sigma peaks will form from primordial material. Thus, even though the oldest stars are centrally concentrated as they originated in the early collapsing, densest regions, primordial stars would be found distributed throughout the halo. Thus, the lack of observed metal free stars can be taken as directly constraining the Population III IMF, and the lowest metallicity observed stars can be interpreted as holding clues to the chemical yields of Pop III stars.

We carry out cosmological chemodynamical simulations with different strengths of supernova (SN) feedback and study how galactic winds from star-forming galaxies affect the features of hydrogen (HI) and metal (CIV and OVI) absorption systems in the intergalactic medium at high redshift. We find that the outflows tend to escape to low density regions, and hardly affect the dense filaments visible in HI absorption. As a result, the strength of HI absorption near galaxies is not reduced by galactic winds, but even slightly increases. We also find that a lack of HI absorption for lines of sight (LOS) close to galaxies, as found by Adelberger et al., can be created by hot gas around the galaxies induced by accretion shock heating. In contrast to HI, metal absorption systems are sensitive to the presence of winds. The models without feedback can produce the strong CIV and OVI absorption lines in LOS within 50 kpc from galaxies, while strong SN feedback is capable of creating strong CIV and OVI lines out to about twice that distance. We also analyze the mean transmissivity of HI, CIV, and OVI within 1 h$^{-1}$ Mpc from star-forming galaxies. The probability distribution of the transmissivity of HI is independent of the strength of SN feedback, but strong feedback produces LOS with lower transmissivity of metal lines. Additionally, strong feedback can produce strong OVI lines even in cases where HI absorption is weak. We conclude that OVI is probably the best tracer for galactic winds at high redshift.

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.

We study the spatial distribution of Galactic metal-free stars by combining an extremely high-resolution (7.8 X 10^5 solar masses per particle) Cold Dark Matter N-body simulation of the Milky-Way with a semi-analytic model of metal enrichment. This approach allows us to resolve halos with virial temperatures down to the 10^4K atomic cooling limit, and it is sufficiently flexible to make a number of robust conclusions, despite the extremely uncertain properties of the first stars. Galactic metal-free stars are formed over a large redshift range, which peaks at z~10, but continues down to z~5, contributing stars at wide range of Galactocentric radii. Stars containing only metals generated by primordial stars are similarly widespread. Neither changing the efficiency of metal dispersal by two orders of magnitude, nor drastically changing the approximations in our semi-analytical model can affect these result. Thus, if they have sufficiently long lifetimes, a significant number of stars formed in initially primordial star clusters should be found in the nearby Galactic halo regardless of the specifics of metal-free star formation. Observations of metal abundances in Galactic halo stars should be taken as directly constraining the properties of primordial stars, and the lack of metal-free halo stars today should be taken as strongly suggesting a 0.8 solar mass lower limit on the primordial initial mass function.

We employ N-body, smoothed particle hydrodynamical simulations, including detailed treatment of chemical enrichment, to follow a gas-rich merger which results in a galaxy with disk morphology. We trace the kinematic, structural and chemical properties of stars formed before, during, and after the merger. We show that such a merger produces two exponential disk components, with the older, hotter component having a scale-length 20% larger than the later-forming, cold disk. Rapid star formation during the merger quickly enriches the protogalactic gas reservoir, resulting in high metallicities of the forming stars. These stars form from gas largely polluted by Type II supernovae, which form rapidly in the merger-induced starburst. After the merger, a thin disk forms from gas which has had time to be polluted by Type Ia supernovae. Abundance trends are plotted, and we examine the proposal that increased star formation during gas-rich mergers may explain the high alpha-to-iron abundance ratios which exist in the relatively high-metallicity thick disk component of the Milky Way.

We review recent developments in the field of chemodynamical simulations of elliptical galaxies, highlighting (in an admittedly biased fashion) the work conducted with our cosmological N-body/SPH code GCD+. We have demonstrated previously the recovery of several primary integrated early-type system scaling relations (e.g. colour-magnitude relation, L_X-T_X-[Fe/H]_X) when employing a phenomenological AGN heating scheme in conjunction with a self-consistent treatment of star formation, supernovae feedback, radiative cooling, chemical enrichment, and stellar/X-ray population synthesis. Here we emphasise characteristics derived from the full spatial information contained within the simulated dataset, including stellar and coronal morphologies, metallicity distribution functions, and abundance gradients.

We report that neutral hydrogen (HI) gas clouds, resembling High Velocity Clouds (HVCs) observed in the Milky Way (MW), appear in MW-sized disk galaxies formed in high-resolution Lambda Cold Dark Matter (LCDM) cosmological simulations which include gas-dynamics, radiative cooling, star formation, supernova feedback, and metal enrichment. Two such disk galaxies are analyzed, and HI column density and velocity distributions in all-sky Aitoff projections are constructed. The simulations demonstrate that LCDM is able to create galaxies with sufficient numbers of anomalous velocity gas clouds consistent with the HVCs observed within the MW, and that they are found within a galactocentric radius of 150 kpc. We also find that one of the galaxies has a polar gas ring, with radius 30 kpc, which appears as a large structure of HVCs in the Aitoff projection. Such large structures may share an origin similar to extended HVCs observed in the MW, such as Complex C.

A suite of high-resolution N-body simulations of the Magellanic Clouds -- Milky Way system are presented and compared directly with newly available data from the HI Parkes All-Sky Survey (HIPASS). We show that the interaction between Small and Large Magellanic Clouds results in both a spatial and kinematical bifurcation of both the Stream and the Leading Arm. The spatial bifurcation of the Stream is readily apparent in the HIPASS data, and the kinematical bifurcation is also tentatively identified. This bifurcation provides strong support for the tidal disruption origin for the Magellanic Stream. A fiducial model for the Magellanic Clouds is presented upon completion of an extensive parameter survey of the potential orbital configurations of the Magellanic Clouds and the viable initial boundary conditions for the disc of the Small Magellanic Cloud. The impact of the choice of these critical parameters upon the final configurations of the Stream and Leading Arm is detailed.

We use a cosmological numerical simulation to study the tidal features produced by a minor merger with an elliptical galaxy. We find that the simulated tidal features are quantitatively similar to the red tidal features, i.e., dry tidal features, recently found in deep images of elliptical galaxies at intermediate redshifts. The minor merger in our simulation does not trigger star formation due to active galactic nuclei heating. Therefore, both the tidal features and the host galaxy are red, i.e. a dry minor merger. The stellar mass of the infalling satellite galaxy is about 10^10 Msun, and the tidal debris reach the surface brightness of mu_R~27 mag arcsec^-2. Thus, we conclude that tidal debris from minor mergers can explain the observed dry tidal features in ellipticals at intermediate redshifts, although other mechanisms (such as major dry mergers) may also be important.

We study the chemical and kinematic properties of the first galaxies which formed at a high redshift, using high resolution cosmological numerical simulations, and compared them with the recent observational results for the Sculptor dwarf spheroidal galaxy by Tolstoy et al., who found two distinct stellar populations: the lower metallicity stars are more spatially extended and possess a higher velocity dispersion than the higher metallicity stars. Our calculations reproduce these observations as the result of a steep metallicity gradient, within a single populations, induced by dissipative collapse of the gas component. We also predict strong [N/O] enhancements in the lowest metallicity stars in dwarf spheroidals, due to the preferential retention of ejected gas from intermediate mass stars, compared to Type II supernovae.

Increasingly large populations of disk galaxies are now being observed at increasingly high redshifts, providing new constraints on our knowledge of how such galaxies evolve. Are these observations consistent with a cosmology in which structures form hierarchically? To probe this question, we employ SPH/N-body galaxy scale simulations of late-type galaxies. We examine the evolution of these simulated disk galaxies from redshift 1 to 0, looking at the mass-size and luminosity-size relations, and the thickness parameter, defined as the ratio of scale-height to scale-length. The structural parameters of our simulated disks settle down quickly, and after redshift z=1 the galaxies evolve to become only slightly flatter. Our present day simulated galaxies are larger, more massive, less bright, and redder than at z=1. The inside-out nature of the growth of our simulated galaxies reduces, and perhaps eliminates, expectations of evolution in the size-mass relation.

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.

Observations suggest that the structural parameters of disk galaxies have not changed greatly since redshift 1. We examine whether these observations are consistent with a cosmology in which structures form hierarchically. We use SPH/N-body galaxy-scale simulations to simulate the formation and evolution of Milky-Way-like disk galaxies by fragmentation, followed by hierarchical merging. The simulated galaxies have a thick disk, that forms in a period of chaotic merging at high redshift, during which a large amount of alpha-elements are produced, and a thin disk, that forms later and has a higher metallicity. Our simulated disks settle down quickly and do not evolve much since redshift z~1, mostly because no major mergers take place between z=1 and z=0. During this period, the disk radius increases (inside-out growth) while its thickness remains constant. These results are consistent with observations of disk galaxies at low and high redshift.

The stellar halos of spiral galaxies bear important chemo-dynamical signatures of galaxy formation. We present here the analysis of 89 semi-cosmological spiral galaxy simulations, spanning ~ 4 magnitudes in total galactic luminosity. These simulations sample a wide variety of merging histories and show significant dispersion in halo metallicity at a given total luminosity - more than a factor of ten in metallicity. Our preliminary analysis suggests that galaxies with a more extended merging history possess halos which have younger and more metal rich stellar populations than the stellar halos associated with galaxies with a more abbreviated assembly. A correlation between halo metallicity and its surface brightness has also been found, reflecting the correlation between halo metallicity and its stellar mass. Our simulations are compared with recent Hubble Space Telescope observations of resolved stellar halos in nearby spirals.

Seven cosmological hydrodynamic simulations of disk galaxy formation are analyzed to determine the alignment of the disk within the dark matter halo and the internal structure of the halo. We find that the orientation of the outer halo, beyond ~0.1 r_vir, is unaffected by the presence of the disk. In contrast, the inner halo is aligned such that the halo minor axis aligns with the disk axis. The relative orientation of these two regions of the halo are uncorrelated. The alignment of the disk and inner halo appears to take place simultaneously through their joint evolution. The disconnect between these two regions of the halo should be taken into account when modelling tidal streams in the halos of disk galaxies and when calculating intrinsic alignments of disk galaxies based on the properties of dark matter halos.

Despite having entered the era of "precision cosmology," the formation of galaxies within the favored CDM cosmological paradigm remains problematic. By relating our N-body/smooth particle hydro-dynamical simulation to an extensive range of Galactic and extragalactic observations, we shed light on the formation and evolution of the Milky Way and other late type galaxies. In light of recent observations of the stellar populations of extragalactic thick disks, we examine the proposal that the thick disk formed in a high redshift period characterized by gas rich merging. We show such a scenario to be consistent with color observations. We then follow the evolution of structural parameters of the subsequently formed disk galaxy, from redshift ~1 to the present. Consistent with observation, little evolution in the ratio of scale-height to scale-length is found in our simulated galaxy, despite its somewhat chaotic origins. The simulated galaxy in this report forms part of a larger study of a suite of galaxies, with which these issues are to be studied in detail.

The recently emerging conviction that thick disks are prevalent in disk galaxies, and their seemingly ubiquitous old ages, means that the formation of the thick disk, perhaps more than any other component, holds the key to unravelling the evolution of the Milky Way, and indeed all disk galaxies. In Paper I, we proposed that the thick disk was formed in an epoch of gas rich mergers, at high redshift. This hypothesis was based on comparing N-body/SPH simulations to a variety of Galactic and extragalactic observations, including stellar kinematics, ages and chemical properties.Here examine our thick disk formation scenario in light of the most recent observations of extragalactic thick disks. In agreement, our simulted thick disks are old and relatively metal rich, with V-I colors that do not vary significantly with distance from the plane. Further, we show that our proposal results in an enhancement of alpha-elements in thick disk stars as compared with thin disk stars, consistent with observations of the relevant populations of the Milky Way. We also find that our scenario naturally leads to the formation of an old metal weak stellar halo population with high alpha-element abundances.

We study the effect of active galactic nuclei (AGN) heating on the chemodynamical evolution of elliptical galaxies and their X-ray and optical properties using high-resolution Lambda-dominated cold dark matter cosmological simulations. Our model considers an AGN as being ``active'' when a convergent gas inflow condition exists within a galaxy's nucleus; otherwise, the AGN is assumed to remain dormant. This induces a self-regulated activity for the AGN, the result of which leads to a stable hot corona and the suppression of significant late-time star formation - characteristics not encountered in traditional chemodynamical models of ellipticals. These properties of our AGN heating model leads to a system consistent with both the X-ray and optical properties of comparable elliptical galaxies.

We present a new chemical evolution model meant to be a first step in the self-consistent study of both optical and X-ray properties of elliptical galaxies. Detailed cooling and heating processes in the interstellar medium are taken into account using a mono-phase one-zone treatment which allows a more reliable modelling of the galactic wind regime with respect to previous work. The model successfully reproduces simultaneously the mass-metallicity, colour-magnitude, the L_X - L_B and the L_X - T relations, as well as the observed trend of the [Mg/Fe] ratio as a function of sigma, by adopting the prescriptions of Pipino & Matteucci (2004) for the gas infall and star formation timescales. We found that a late secondary accretion of gas from the environment plays a fundamental role in driving the L_X - L_B and L_X - T relations and can explain their large observational scatter. The iron discrepancy, namely the too high predicted iron abundance in X-ray haloes of ellipticals compared to observations, still persists. On the other hand, we predict [O/Fe] in the ISM which is in good agreement with the most recent observations. We suggest possible mechanisms acting on a galactic scale which may solve the iron discrepancy. In particular, mixing of gas driven by AGNs may preserve the gas mass (and thus the X-ray luminosity) while diluting the iron abundance. New predictions for the amounts of iron, oxygen and energy ejected into the intracluster medium (ICM) are presented and we conclude that Type Ia supernovae (SNe Ia) play a fundamental role in the ICM enrichment. SNe Ia activity, in fact, may power a galactic wind lasting for a considerable amount of the galactic lifetime, even in the case for which the efficiency of energy transfer into the ISM per SN Ia event is less than unity.

We present a detailed study of the real and integrals-of-motion space distributions of a satellite obtained from a self-consistent high-resolution simulation of a galaxy cluster and re-simulated using various analytical halo potentials. We found that the disrupted satellite appears as a coherent structure in integrals-of-motion space in all models (``live'' and analytical potential) although the distribution is significantly smeared for the live host halo. Further the primary mechanism for this smearing is the mass growth of the host, which changes both the energy and angular momentum of the satellite debris. Hence, this must be considered when searching for (stellar) streams with future observational experiments such as RAVE and GAIA.

The chemical evolution history of a galaxy hides clues about how it formed and has been changing through time. We have studied the chemical evolution history of the Milky Way (MW) and Andromeda (M31) to find which are common features in the chemical evolution of disc galaxies as well as which are galaxy-dependent. We use a semi-analytic multi-zone chemical evolution model. Such models have succeeded in explaining the mean trends of the observed chemical properties in these two Local Group spiral galaxies with similar mass and morphology. Our results suggest that while the evolution of the MW and M31 shares general similarities, differences in the formation history are required to explain the observations in detail. In particular, we found that the observed higher metallicity in the M31 halo can be explained by either a) a higher halo star formation efficiency or b) a larger reservoir of infalling halo gas with a longer halo formation phase. These two different pictures would lead to a) a higher [O/Fe] at low metallicities or b) younger stellar populations in the M31 halo, respectively. Both pictures result in a more massive stellar halo in M31, which suggests a possible correlation between the halo metallicity and its stellar mass.

Using self-consistent cosmological simulations of disc galaxy formation, we analyse the 1.4 GHz radio flux from high-redshift progenitors of present-day normal spirals within the context of present-day and planned next-generation observational facilities. We demonstrate that while current radio facilities such as the Very Large Array (VLA) are unlikely to trace these progenitors beyond redshifts z<0.2, future facilities such as the Square Kilometer Array (SKA) will readily probe their characteristics to redshifts z<2, and are likely to provide detections beyond z~3. We also demonstrate that the progenitors of present-day cD galaxies can emit in excess of 10 uJy of flux at redshifts z>1, and may be a non-negligible contributor to the micro-Jansky source counts derived from current deep VLA cm-wave surveys.

The Fifth Workshop on Galactic Chemodynamics (GCDV) is the pre-eminent workshop devoted exclusively to observational and theoretical research into the formation and evolution of galaxies. GCDV brought together 46 researchers from 13 countries for an intensive workshop at Swinburne University in July 2003. Each of the 28 papers included in the workshop proceedings have been formally refereed in their entirety, published as a special volume of the Publications of the Astronomical Society of Australia, and available through NASA's Astronomical Data System (ADS).

The rotation curves of low surface brightness (LSB) galaxies suggest that they possess significantly higher mass-to-light (M/L) ratios than their high surface brightness counterparts, indicating that LSB galaxies may be dark matter dominated. This interpretation is hampered by the difficulty of disentangling the disc and dark halo contributions from the disc dynamics of LSB galaxies. Recently, Fuchs (2002) has attempted such a disentanglement using spiral arm density wave and swing amplification theory, allowing an independent measurement of the disc mass; this work suggests that LSB discs are significantly more massive than previously believed. This would considerably reduce the amount of matter required in the dark halos in fitting the rotation curves. Interestingly, the high mass-to-light ratios derived for the discs appear inconsistent with standard stellar population synthesis models. In this paper, we investigate whether the high M/L ratios for the Fuchs LSB discs might be understood by adopting a very ``bottom heavy'' initial mass function (IMF). We find that an IMF with a power law exponent of around alpha=3.85 (compared to the standard Salpeter IMF, alpha=2.35) is sufficient to explain the unusually high M/L ratios of the Fuchs sample. Within the context of the models, the blue colours ((B-R)_0 < 1.0) of the sample galaxies result from being metal-poor ([Fe/H] = -1.5 ~ -1.0) and having undergone recent (~1-3 Gyr ago) star formation.

The disk galaxy simulated using our chemo-dynamical galaxy formation code, GCD+, is shown to have a thick disk component. This is evidenced by the velocity dispersion versus age relation for solar neighbourhood stars, which clearly shows an abrupt increase in velocity dispersion at lookback time of approximately 8 Gyrs, and is in excellent agreement with observation. These thick disk stars are formed from gas which is accreted to the galaxy during a chaotic period of hierarchical clustering at high redshift. This formation scenario is shown to be consistent with observations of both the Milky Way and extragalactic thick disks.

Hierarchical clustering represents the favoured paradigm for galaxy formation throughout the Universe; due to its proximity, the Magellanic system offers one of the few opportunities for astrophysicists to decompose the full six-dimensional phase-space history of a satellite in the midst of being cannibalised by its host galaxy. The availability of improved observational data for the Magellanic Stream and parallel advances in computational power has led us to revisit the canonical tidal model describing the disruption of the Small Magellanic Cloud and the consequent formation of the Stream. We suggest improvements to the tidal model in light of these recent advances.

Using stellar population synthesis techniques, we explore the photometric signatures of white dwarf progenitor dominated galactic halos, in order to constrain the fraction of halo mass that may be locked-up in white dwarf stellar remnants. We first construct a 10^9 M_sun stellar halo using the canonical Salpeter initial stellar mass distribution, and then allow for an additional component of low- and intermediate-mass stars, which ultimately give rise to white dwarf remnants. Microlensing observations towards the Large Magellanic Cloud, coupled with several ground-based proper motion surveys, have led to claims that in excess of 20% of the dynamical mass of the halo (10^12 M_sun) might be found in white dwarfs. Our results indicate that (1) even if only 1% of the dynamical mass of the dark halo today could be attributed to white dwarfs, their main sequence progenitors at high redshift (z ~ 3) would have resulted in halos more than 100 times more luminous than those expected from conventional initial mass functions alone, and (2) any putative halo white dwarf progenitor dominated initial mass function component, regardless of its dynamical importance, would be virtually impossible to detect at the present-day, due to its extremely faint surface brightness.

We study the chemodynamical evolution of elliptical galaxies and their X-ray and optical properties 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 abundances of heavy elements in both the stellar and gas components of galaxies. A giant elliptical galaxy formed in one of the simulations is quantitatively compared with the observational data in both the X-ray and optical regime. We implement a treatment of both pre-heating and AGN heating in this simulation, and examine the effect of these processes on elliptical galaxy formation. We find that the adopted pre-heating (T=10^7 K at z=4) is not strong enough to explain the observed X-ray or optical properties. On the other hand, our AGN heating model in which the gas inflow induces the AGN heating is consistent with both the X-ray and optical properties.

The primary present-day observables upon which theories of galaxy evolution are based are a system's morphology, dynamics, colour, and chemistry. Individually, each provides an important constraint to any given model; in concert, the four represent a fundamental (intractable) boundary condition for chemodynamical simulations. We review the current state-of-the-art semi-analytical and chemodynamical models for the Milky Way, emphasising the strengths and weaknesses of both approaches.

How do late type spiral galaxies form within the context of a CDM cosmology? We contrast N-body, smoothed particle hydrodynamical simulations of galaxy formation which employ two different supernova feedback mechanisms. Observed mass and metallicity distributions of the stellar halos of the Milky Way and M31 provide constraints on these models. A strong feedback model, incorporating an adiabatic phase in star burst regions, better reproduces these observational properties than our comparative thermal feedback model. This demonstrates the importance of energy feedback in regulating star formation in small systems, which collapse at early epochs, in the evolution of late type disk galaxies.

We explore a scenario where metal-poor globular clusters (GCs) are enriched by the first supernovae in the Universe. If the first stars in a 10^7 Msun dark halo were very massive (>180 Msun), then a pair instability supernova from a single massive star can produce sufficient iron to enrich 10^6 Msun of pristine, primordial gas to [Fe/H] ~ -2. In such a scenario, where a single massive star acts as a seed for halo GCs, the accurate abundance analysis of GC stars would allow a direct measurement of the Population III initial mass. Using the latest theoretical yields for zero metallicity stars in the mass range 140-260 Msun, we find that the metals expelled from a ~230 Msun star are consistent with [Si/Fe] and [Ca/Fe] observed in GC stars. However, no single star in this mass range can simultaneously explain all halo GC heavy-element abundance ratios, such as [V/Fe], [Ti/Fe] and [Ni/Fe]. These require a combination masses for the Population III stellar progenitors. The various observational consequences of this scenario are discussed.

We study the chemodynamical evolution of elliptical galaxies and their X-ray and optical properties 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 abundances of heavy element in both the stellar and gas components of galaxies. A giant elliptical galaxy formed in the simulations is quantitatively compared with the observational data in both the X-ray and optical regime. 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. Optical properties are derived by the population synthesis for the stellar component. We find that radiative cooling is important to interpret the observed X-ray luminosity, temperature, and metallicity of the hot gas 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. Time variation and radial dependence of X-ray properties and abundance ratios, such as [O/Fe] and [Si/Fe], of the X-ray emitting hot gas are also discussed.

We test the hypothesis that High-Velocity gas cloud Complex C is actually a high-latitude spiral arm extension in the direction of the Galactic warp, as opposed to the standard interpretation - that of a once extragalactic, but now infalling, gas cloud. A parallel Tree N-body code was employed to simulate the tidal interaction of a satellite perturber with the Milky Way. We find that a model incorporating a perturber of the mass of the Large Magellanic Cloud on a south-to-north polar orbit, crossing the disk at ~15 kpc, does yield a high- velocity, high-latitude, extension consistent with the spatial, kinematical and column density properties of Complex C. Unless this massive satellite remains undiscovered because of either a fortuitous alignment with the Galactic bulge (feasible within the framework of the model), or the lack of any associated baryonic component, we conclude that this alternative interpretation appears unlikely.

Observational evidence has suggested the possibility of a Galactic halo which is dominated by white dwarfs (WDs). While debate continues concerning the interpretation of this evidence, it is clear that an initial mass function (IMF) biased heavily toward WD precursors (1 < m/Msol < 8), at least in the early Universe, would be necessary in generating such a halo. Within the framework of homogeneous, closed-box models of Galaxy formation, such biased IMFs lead to an unavoidable overproduction of carbon and nitrogen relative to oxygen (as measured against the abundance patterns in the oldest stars of the Milky Way). Using a three-dimensional Tree N-body smoothed particle hydrodynamics code, we study the dynamics and chemical evolution of a galaxy with different IMFs. Both invariant and metallicity-dependent IMFs are considered. Our variable IMF model invokes a WD-precursor-dominated IMF for metallicities less than 5% solar (primarily the Galactic halo), and the canonical Salpeter IMF otherwise (primarily the disk). Halo WD density distributions and C,N/O abundance patterns are presented. While Galactic haloes comprised of ~5% (by mass) of WDs are not supported by our simulations, mass fractions of ~1-2% cannot be ruled out. This conclusion is consistent with the present-day observational constraints.

The present day chemical and dynamical properties of the Milky Way bear the imprint of the Galaxy's formation and evolutionary history. One of the most enduring and critical debates surrounding Galactic evolution is that regarding the competition between ``satellite accretion'' and ``monolithic collapse''; the apparent strong correlation between orbital eccentricity and metallicity of halo stars was originally used as supporting evidence for the latter. While modern-day unbiased samples no longer support the claims for a significant correlation, recent evidence has been presented by Chiba & Beers (2000,AJ,119,2843) for the existence of a minor population of high-eccentricity metal-deficient halo stars. It has been suggested that these stars represent the signature of a rapid (if minor) collapse phase in the Galaxy's history. Employing velocity- and integrals of motion-phase space projections of these stars, coupled with a series of N-body/Smoothed Particle Hydrodynamic (SPH) chemodynamical simulations, we suggest an alternative mechanism for creating such stars may be the recent accretion of a polar orbit dwarf galaxy.

We have developed a new galactic chemo-dynamical evolution code, called GCD+, for studies of galaxy formation and evolution. This code is based on our original three-dimensional tree N-body/smoothed particle hydrodynamics code which includes self-gravity, hydrodynamics, radiative cooling, star formation, supernova feedback, and metal enrichment. GCD+ includes a new Type II (SNe II) and Ia (SNe Ia) supernovae model taking into account the lifetime of progenitor stars, and chemical enrichment from intermediate mass stars. We apply GCD+ to simulations of elliptical galaxy formation, and examine the colour-magnitude relation (CMR), the Kormendy relation, and the [Mg/Fe]--magnitude relation of simulation end-products. GCD+ is a useful and unique tool which enables us to compare simulation results with the observational data directly and quantitatively. Our simulation confirm the results of Kawata (2001) who uses a simpler chemo-dynamical evolution code. We newly find that radiative cooling becomes more efficient and thus the gas infall rate increases, with decreasing mass of galaxies, which contributes to the slope of the CMR. In addition, the sophisticated treatments of both SNe II and SNe Ia in GCD+ show that feedback from SNe Ia plays a crucial role in the evolution of elliptical galaxies. We conclude that the feedback effect of SNe Ia should not be ignored in studying the evolution of 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.

The present-day chemical and dynamical properties of the Milky Way are signatures of the Galaxy's formation and evolution. Using a self consistent chemodynamical evolution code we examine these properties within the currently favoured paradigm for galaxy formation - hierarchical clustering within a CDM cosmology. Our Tree N-body/Smoothed Particle Hydrodynamics code includes a self-consistent treatment of gravity, hydrodynamics, radiative cooling, star formation, supernova feedback and chemical enrichment. Two models are described which explore the role of small-scale density perturbations in driving the evolution of structure within the Milky Way. The relationship between metallicity and kinematics of halo stars are quantified and the implications for galaxy formation discussed. While high-eccentricity halo stars have previously been considered a signature of ``rapid collapse'', we suggest that many such stars may have come from recently accreted satellites.

We are in a privileged location in the Universe which allows us to observe galactic interactions from close range -- the merger of our two nearest dwarf satellite galaxies, the LMC and SMC. It is important to understand the local merger process before we can have confidence in understanding mergers at high redshift. We present high resolution Nbody+SPH simulations of the disruption of the LMC and SMC and the formation of the Magellanic Stream, and discuss the implications for galaxy formation and evolution.

Halo initial mass functions (IMFs) heavily biased toward white dwarf (WD) precursors ~1-8 Msun have been suggested as a suitable mechanism for explaining MACHO events. However, by simple chemical evolution argument, Gibson & Mould (1997; GM97) pointed out that such WD-heavy IMFs (wdIMF) cause the overproduction of carbon and nitrogen. We re-examine this problem using numerical simulations.

We investigate the chemo-dynamical evolution of elliptical galaxies, to understand the origin of the mass-dependence of photometric properties such as the colour-magnitude relation (CMR). Our three-dimensional TREE N-body/SPH numerical simulation takes into account both Type II and Type Ia supernovae and follows the evolution of the abundances of several chemical elements. We derive the photometric properties of the simulation end-products and compare them with the observed CMR.

Motivated by the apparent order-of-magnitude discrepancy between the observed number of Local Group satellite galaxies, and that predicted by Lambda-CDM hierarchical clustering cosmologies, we explore an alternate suggestion - perhaps the missing satellites are not actually ``missing'', but are instead ``in disguise''. The disguise we consider here is that of the classical HI High-Velocity Clouds. Is it possible that what have been thought of traditionally as a ``Galactic'' phenomenon, are actually the building blocks of the Local Group? We discuss the strengths and weaknesses of this hypothesis, and highlight avenues of future research which may provide an unequivocal resolution to this contentious issue.

We numerically investigate the dynamical and chemical processes of the formation of elliptical galaxies in a cold dark matter (CDM) universe, in order to understand the origin of the mass-dependence of the photometric properties of elliptical galaxies. Our three-dimensional TREE N-body/SPH numerical simulations of elliptical galaxy formation take into account both Type II (SNe II) and Type Ia (SNe Ia) supernovae (SNe) and follow the time evolution of the abundances of several chemical elements (C, O, Ne, Mg, Si, and Fe). Moreover we compare different strengths of SNe feedback.In combination with stellar population synthesis, we derive the photometric properties of simulation end-products, including the magnitude, color, half-light radius, and abundance ratios, and compare them with the observed scaling relations directly and quantitatively. We find that the extremely strong influence of SNe is required to reproduce the observed color-magnitude relation (CMR), where we assume each SN yields energy of 4x10^51 ergs and that 90% of this energy is ejected as kinetic feedback. The feedback affects the evolution of lower mass systems more strongly and induces the galactic wind by which a larger fraction of gas is blown out in a lower mass system. Finally higher mass systems become more metal rich and have redder colors than lower mass systems. We emphasize based on our simulation results that the galactic wind is triggered mainly by SNe Ia rather than SNe II. In addition we examined the Kormendy relation, which prescribes the size of elliptical galaxies, and the [Mg/Fe]--magnitude relation, which provides a strong constraint on the star formation history.

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/\sigma). 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.

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.