We present a set of self-consistent chemo-dynamical simulations of MW-type galaxies formation to study the origin of the bimodality of $\alpha$-elements in stellar populations. We explore how the bimodality is related to the geometrically and kinematically defined stellar discs, gas accretion and radial migration. We find that the two $\alpha$-sequences are formed in quite different physical environments. The high-$\alpha$ sequence is formed early from a burst of star formation (SF) in a turbulent, compact gaseous disc which forms a thick disc. The low-$\alpha$ stellar populations is the result of quiescent SF supported by the slow accretion of enriched gas onto a radially extended thin disc. Stellar feedback-driven outflows during the formation of the thick disc are responsible for the enrichment of the surrounding gaseous halo, which subsequently feeds the disc on a longer time-scale. During the thin disc phase, chemical evolution reaches an equilibrium metallicity and abundance, where the stars pile-up. This equilibrium metallicity decreases towards the outer disc, generating the ridgeline that forms the low-$\alpha$ sequence. We identify a second mechanism capable of creating a low-$\alpha$ sequence in one of our simulations. Rapid shutdown of the SF, provoked by the feedback at the end of the thick disc phase, suppresses the chemical enrichment of the halo gas, which, once accreted onto the star-forming disc, dilutes the ISM at the beginning of the thin disc formation. Both mechanisms can operate in a galaxy, but the former is expected to occur when SF efficiency ceases to be dominated by the formation of the thick disc, while the latter can occur in the inner regions. Being the result of the presence of low and high gas density environments, the bimodality is independent of any particular merger history, suggesting that it could be much more widespread than has been claimed.

Understanding radial migration is a crucial point to build relevant chemical and dynamical evolution models of the Milky Way disk. In this paper, we analyze a high-resolution N-body simulation of a Milky Way-type galaxy to study the role that the slowing down of a stellar bar has is generating migration from the inner to the outer disk. Stellar particles are trapped by the main resonances (corotation and Outer Lindblad resonance) which then propagate outwards across the disk due to the bar slowing down. Once the bar strength reaches its maximal amplitude, some of the stars, delivered to the outer disk, escape the resonances and some of them settle on nearly circular orbits. The number of the escaped stars gradually increases also due to the decrease of the bar strength when the boxy/peanut bulge forms. We show that this mechanism is not limited only to stars on nearly circular orbits: also stars initially on more eccentric orbits can be transferred outwards (out to the OLR location) and can end up on nearly circular orbits. Therefore, the propagation of the bar resonances outwards can induce the circularization of the orbits of some of the migrating stars. The mechanism investigated in this paper can explain the presence of metal-rich stars at the solar vicinity and more generally in the outer galactic disk. Our dynamical model predicts that up to 3% of stars in between of corotation and the OLR can be formed in the innermost region of the Milky Way. The epoch of the Milky Way bar formation can be potentially constrained by analyzing the age distribution of the most metal-rich stars at the solar vicinity.

We analyse a set of very metal-poor stars for which accurate chemical abundances have been obtained as part of the ESO Large Program "First stars" in the light of the Gaia DR2 data. The kinematics and orbital properties of the stars in the sample show they probably belong to the thick disc, partially heated to halo kinematics, and to the accreted Nissen & Schuster-Gaia Sausage-Enceladus satellite. The continuity of these properties with stars at both higher ($\rm [Fe/H]>-2$) and lower metallicities ($\rm [Fe/H]<-4.$) suggests that the Galaxy at $\rm [Fe/H] \lesssim -0.5$ and down to at least $\rm [Fe/H]\sim-6$ is dominated by these two populations. In particular, we show that the disc extends continuously from $\rm [Fe/H] \le -4$ (where stars with disc-like kinematics have been recently discovered) up to $\rm [Fe/H] \ge -2$, the metallicity regime of the Galactic thick disc. There exists indeed an "ultra-metal poor thick disc", which constitutes the extremely metal-poor tail of the canonical Galactic thick disc, and which extends this latter from $\rm [Fe/H] \sim -0.5$ up to the most metal-poor stars discovered in the Galaxy up to date. These results suggest that the disc may be the main, and possibly the only stellar population that has formed in the Galaxy at these metallicities. This requires that the dissipative collapse that led to the formation of the old Galactic disc must have been extremely fast. We also discuss these results in the light of recent simulation efforts made to reproduce the first stages of Milky Way-type galaxies.

Boxy, peanut- or X-shaped "bulges" are observed in a large fraction of barred galaxies viewed in, or close to, edge-on projection, as well as in the Milky Way. They are the product of dynamical instabilities occurring in stellar bars, which cause the latter to buckle and thicken vertically. Recent studies have found nearby galaxies that harbour two such features arising at different radial scales, in a nested configuration. In this paper we explore the formation of such double peanuts, using a collisionless N-body simulation of a pure disc evolving in isolation within a live dark matter halo, which we analyse in a completely analogous way to observations of real galaxies. In the simulation we find a stable double configuration consisting of two X/peanut structures associated to the same galactic bar - rotating with the same pattern speed - but with different morphology, formation time, and evolution. The inner, conventional peanut-shaped structure forms early via the buckling of the bar, and experiences little evolution once it stabilises. This feature is consistent in terms of size, strength and morphology, with peanut structures observed in nearby galaxies. The outer structure, however, displays a strong X, or "bow-tie", morphology. It forms just after the inner peanut, and gradually extends in time (within 1 to 1.5 Gyr) to almost the end of the bar, a radial scale where ansae occur. We conclude that, although both structures form, and are dynamically coupled to, the same bar, they are supported by inherently different mechanisms.

In this paper we introduce a new method for analysing Milky Way phase-space which allows us to reveal the imprint left by the Milky Way bar and spiral arms on the stars with full phase-space data in Gaia Data Release 2. The unprecedented quality and extended spatial coverage of these data enable us to discover six prominent stellar density structures in the disc to a distance of 5 kpc from the Sun. Four of these structures correspond to the spiral arms detected previously in the gas and young stars (Scutum-Centaurus, Sagittarius, Local and Perseus). The remaining two are associated with the main resonances of the Milky Way bar where corotation is placed at around 6.2 kpc and the outer Lindblad resonance beyond the Solar radius, at around 9 kpc. For the first time we provide evidence of the imprint left by spiral arms and resonances in the stellar densities not relying on a specific tracer, through enhancing the signatures left by these asymmetries. Our method offers new avenues for studying how the stellar populations in our Galaxy are shaped.

Previous studies based on the analysis of Gaia DR2 data have revealed that accreted stars, possibly originating from a single progenitor satellite, are a significant component of the halo of our Galaxy, potentially constituting most of the halo stars at $\rm [Fe/H] < -1$ within a few kpc from the Sun and beyond. In this paper, we couple astrometric data from Gaia DR2 with elemental abundances from APOGEE DR14 to characterize the kinematics and chemistry of in-situ and accreted populations up to $\rm [Fe/H] \sim -2$. Accreted stars appear to significantly impact the Galactic chemo-kinematic relations, not only at $\rm [Fe/H] < -1$, but also at metallicities typical of the thick and metal-poor thin discs. They constitute about 60% of all stars at $\rm [Fe/H] < -1$, the remaining 40% being made of (metal-weak) thick disc stars. We find that the stellar kinematic fossil record shows the imprint left by this accretion event which heated the old Galactic disc. We are able to age-date this kinematic imprint, showing that the accretion occurred between 9 and 11 Gyr ago, and that it led to the last significant heating of the Galactic disc. An important fraction of stars with abundances typical of the (metal-rich) thick disc, and heated by this interaction, is now found in the Galactic halo. Indeed about half of the kinematically defined halo at few kpc from the Sun is composed of metal-rich thick disc stars. Moreover, we suggest that this metal-rich thick disc component dominates the stellar halo of the inner Galaxy. The new picture that emerges from this study is one where the standard non-rotating in-situ halo population, the collapsed halo, seems to be more elusive than ever.

We explore the origin of phase-space substructures revealed by the second Gaia data release in the disc of the Milky Way, such as the ridges in the $V_{\phi}$-$r$ plane, the undulations in the $V_{\phi}$-$r$-$V_r$ space and the streams in the $V_{\phi}$-$V_r$ plane. We use a collisionless N-body simulation with co-spatial thin and thick discs, along with orbit integration, to study the orbital structure close to the Outer Lindblad Resonance (OLR) of the bar. We find that a prominent, long-lived ridge is formed in the $V_{\phi}$-$r$ plane due to the OLR which translates to streams in the $V_{\phi}$-$V_r$ plane and examine which closed periodic and trapped librating orbits are responsible for these features. We find that orbits which carry out small librations around the $x_1(1)$ family are preferentially found at negative $V_r$, giving rise to a `horn'-like feature, while orbits with larger libration amplitudes, trapped around the $x_1(2)$ and $x_1(1)$ families, constitute the positive $V_r$ substructure, i.e. the Hercules-like feature. This changing libration amplitude of orbits will translate to a changing ratio of thin/thick disc stars, which could have implications on the metallicity distribution in this plane. We find that a scenario in which the Sun is placed close to the OLR gives rise to a strong asymmetry in $V_r$ in the $V_{\phi}$-$V_r$ plane (i.e. Hercules vs. `the horn') and subsequently to undulations in the $V_{\phi}$-$r$-$V_r$ space. We also explore a scenario in which the Sun is placed closer to the bar corotation and find that the bar perturbation $alone$ cannot give rise to the these features.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

(Abridged) The Maunakea Spectroscopic Explorer (MSE) is an end-to-end science platform for the design, execution and scientific exploitation of spectroscopic surveys. It will unveil the composition and dynamics of the faint Universe and impact nearly every field of astrophysics across all spatial scales, from individual stars to the largest scale structures in the Universe. Major pillars in the science program for MSE include (i) the ultimate Gaia follow-up facility for understanding the chemistry and dynamics of the distant Milky Way, including the outer disk and faint stellar halo at high spectral resolution (ii) galaxy formation and evolution at cosmic noon, via the type of revolutionary surveys that have occurred in the nearby Universe, but now conducted at the peak of the star formation history of the Universe (iii) derivation of the mass of the neutrino and insights into inflationary physics through a cosmological redshift survey that probes a large volume of the Universe with a high galaxy density. MSE is positioned to become a critical hub in the emerging international network of front-line astronomical facilities, with scientific capabilities that naturally complement and extend the scientific power of Gaia, the Large Synoptic Survey Telescope, the Square Kilometer Array, Euclid, WFIRST, the 30m telescopes and many more.

We present a scenario of the chemical enrichment of the solar neighborhood that solves the G-dwarf problem by taking into account constraints on a larger scale. We argue that the Milky Way disk within 10 kpc has been enriched to solar metallicity by a massive stellar population: the thick disk, which itself formed from a massive turbulent gaseous disk. The key new consideration is that the pre-enrichment provided by the thick disk is not related to the mass fraction of this stellar population at the solar radius, as is classically assumed in inside-out scenarios, but is actually related to the formation of the entire massive thick disk, due to the vigorous gas phase mixing that occurred during its formation. Hence, the fact that this population represents only 15-25% of the local stellar surface density today is irrelevant for `solving' the G-dwarf problem. The only condition for this scenario to work is that the thick disk was formed from a turbulent gaseous disk that permitted a homogeneous -- not radially dependent -- distribution of metals, allowing the solar ring to be enriched to solar metallicity. At the solar radius, the gas flowing from the outer disk combined with the solar metallicity gas left over from thick disk formation, providing the fuel necessary to form the thin disk at the correct metallicity to solve the G-dwarf problem. Chemical evolution at R$>$6 kpc, and in particular beyond the solar radius, can be reproduced with the same scheme. These results imply that the local metallicity distribution is not connected to the gas accretion history of the Milky Way. Finally, we argue that the Sun is the result of the evolution typical of stars in the disk beyond $\sim$6 kpc (i.e., also undergoing dilution), and has none of the characteristics of inner disk stars. [Abridged]

Using a single N-body simulation ($N=0.14\times 10^9$) we explore the formation, evolution and spatial variation of the phase-space spirals similar to those recently discovered by Antoja et al. in the Milky Way disk, with Gaia DR2. For the first time in the literature, we use a self-consistent N-body simulation of an isolated Milky Way-type galaxy to show that the phase-space spirals develop naturally from vertical oscillations driven by the buckling of the stellar bar. We claim that the physical mechanism standing behind the observed incomplete phase-space mixing process can be internal and not necessarily due to the perturbation induced by a massive satellite. In our model, the bending oscillations propagate outwards and produce axisymmetric variations of the mean vertical coordinate and of the vertical velocity component. As a consequence, the phase-space wrapping results in the formation of patterns with various morphology across the disk, depending on the bar orientation, distance to the galactic center and time elapsed since the bar buckling. Once bending waves appear, they are supported for a long time via disk self-gravity. The underlying physical mechanism implies the link between in-plane and vertical motion that leads directly to phase-space structures whose amplitude and shape are in remarkable agreement with those of the phase-space spirals observed in the Milky Way disk. In our isolated galaxy simulation, phase-space spirals are still distinguishable, at the solar neighbourhood, 3 Gyr after the buckling phase. The long-lived character of the phase-space spirals generated by the bar buckling instability cast doubts on the timing argument used so far to get back at the time of the onset of the perturbation: phase-space spirals may have been caused by perturbations originated several Gyrs ago, and not as recent as suggested so far.

In the MW bulge, metal-rich stars form a strong bar and are more peanut-shaped than metal-poor stars. It has been recently claimed that this behavior is driven by the initial in-plane radial velocity dispersion of these populations, rather than by their initial vertical random motions. This has led to the suggestion that a thick disc is not necessary to explain the characteristics of the MW bulge. We rediscuss this issue by analyzing two dissipationless N-body simulations of boxy/peanut (b/p)-shaped bulges formed from composite stellar discs, made of kinematically cold and hot stellar populations, and we conclude that initial vertical random motions are as important as in-plane random motions in determining the relative contribution of cold and hot disc populations with height above the plane, the metallicity and age trends. Previous statements emphasizing the dominant role of in-plane motions in determining these trends are not confirmed. However, differences exist in the morphology and strength of the resulting b/p-shaped bulges: a model where disc populations have initially only different in-plane random motions, but similar thickness, results into a b/p bulge where all populations have a similar peanut shape, independently on their initial kinematics, or metallicity. We discuss the reasons behind these differences, and also predict the signatures that these two extreme initial conditions would leave on the vertical age and metallicity gradients of disc stars, outside the bulge region. We conclude that a metal-poor, kinematically (radial and vertical) hot component, that is a thick disc, is necessary in the MW before bar formation, supporting the scenario traced in previous works. [abridged]

Globular clusters (GCs), the oldest stellar systems observed in the Milky Way, have for long been considered single stellar populations. As such, they provided an ideal laboratory to understand stellar dynamics and primordial star formation processes. However, during the last two decades, observations unveiled their real, complex nature. Beside their pristine stars, GCs host one or more helium enriched and possibly younger stellar populations whose formation mechanism is still unknown. Even more puzzling is the existence of GCs showing star by star iron spreads. Using detailed N-body simulations we explore the hypothesis that these anomalies in metallicity could be the result of mutual stripping and mergers between a primordial population of disc GCs. In the first paper of this series we proved, both with analytical arguments and short-term N-body simulations, that disc GCs have larger fly-by and close encounter rates with respect to halo clusters. These interactions lead to mass exchange and even mergers that form new GCs, possibly showing metallicity spreads. Here, by means of long-term direct N-body simulations, we provide predictions on the dynamical properties of GCs that underwent these processes. The comparison of our predictions with available and future observational data could provide insights on the origin of GCs and on the Milky Way build-up history as a whole.

We present the results of a self-consistent $N$-body simulation following the evolution of a primordial population of thick disc globular clusters (GCs). We study how the internal properties of such clusters evolve under the action of mutual interactions, while they orbit a Milky Way-like galaxy. For the first time, through analytical and numerical considerations, we find that physical encounters between disc GCs are a crucial factor that contributed to the shape of the current properties of the Galactic GC system. Close passages or motion on similar orbits may indeed have a significant impact on the internal structure of clusters, producing multiple gravitationally bound sub-populations through the exchange of mass and even mergers. Our model produces two major mergers and a few small mass exchanges between pairs of GCs. Two of our GCs accrete stars from two companions, ending up with three internal sub-populations. We propose these early interactions and mergers between thick disc GCs with slightly different initial chemical compositions as a possible explanation for the presence of the spreads in metallicity observed in some of the massive Milky Ways GCs.

We investigate the nature of the double color-magnitude sequence observed in the Gaia DR2 HR diagram of stars with high transverse velocities. The stars in the reddest-color sequence are likely dominated by the dynamically-hot tail of the thick disk population. Information from Nissen & Schuster (2010) and from the APOGEE survey suggests that stars in the blue-color sequence have elemental abundance patterns that can be explained by this population having a relatively low star-formation efficiency during its formation. In dynamical and orbital spaces, such as the `Toomre diagram', the two sequences show a significant overlap, but with a tendency for stars on the blue-color sequence to dominate regions with no or retrograde rotation and high total orbital energy. In the plane defined by the maximal vertical excursion of the orbits versus their apocenters, stars of both sequences redistribute into discrete wedges. We conclude that stars which are typically assigned to the halo in the solar vicinity are actually both accreted stars lying along the blue sequence in the HR diagram, and the low velocity tail of the old Galactic disk, possibly dynamically heated by past accretion events. Our results imply that a halo population formed in situ and responsible for the early chemical enrichment prior to the formation of the thick disk is yet to be robustly identified, and that what has been defined as the stars of the in situ stellar halo of the Galaxy may be in fact fossil records of its last significant merger.

There is a long-standing debate on the origin of the metal-poor stellar populations of the Milky Way (MW) bulge, with the two leading scenarios being that these populations are either i) part of a classical metal-poor spheroid or ii) the same population as the chemically defined thick disc seen at the Solar neighbourhood. Here we test whether the latter scenario can reproduce the observed chemical properties of the MW bulge. To do so we compare an N-body simulation of a composite (thin+thick) stellar disc -- which evolves secularly to form a bar and a boxy/peanut (b/p) bulge -- to data from APOGEE DR13. This model, in which the thick disc is massive and centrally concentrated, can reproduce the morphology of the metal-rich and metal-poor stellar populations in the bulge, as well as the mean metallicity and [$\alpha$/Fe] maps as obtained from the APOGEE data. It also reproduces the trends, in both longitude and latitude, of the bulge metallicity distribution function (MDF). Additionally, we show that the model predicts small but measurable azimuthal metallicity variations in the inner disc due to the differential mapping of the thin and thick disc in the bar. We therefore see that the chemo-morphological relations of stellar populations in the MW bulge are naturally reproduced by mapping the thin and thick discs of the inner MW into a b/p.

To illustrate the potential of GDR2, we provide a first look at the kinematics of the Milky Way disc, within a radius of several kiloparsecs around the Sun. We benefit for the first time from a sample of 6.4 million F-G-K stars with full 6D phase-space coordinates, precise parallaxes, and precise Galactic cylindrical velocities . From this sample, we extracted a sub-sample of 3.2 million giant stars to map the velocity field of the Galactic disc from $\sim$5~kpc to $\sim$13~kpc from the Galactic centre and up to 2~kpc above and below the plane. We also study the distribution of 0.3 million solar neighbourhood stars ($r < 200$~pc), with median velocity uncertainties of 0.4~km/s, in velocity space and use the full sample to examine how the over-densities evolve in more distant regions. GDR2 allows us to draw 3D maps of the Galactocentric median velocities and velocity dispersions with unprecedented accuracy, precision, and spatial resolution. The maps show the complexity and richness of the velocity field of the galactic disc. We observe streaming motions in all the components of the velocities as well as patterns in the velocity dispersions. For example, we confirm the previously reported negative and positive galactocentric radial velocity gradients in the inner and outer disc, respectively. Here, we see them as part of a non-axisymmetric kinematic oscillation, and we map its azimuthal and vertical behaviour. We also witness a new global arrangement of stars in the velocity plane of the solar neighbourhood and in distant regions in which stars are organised in thin substructures with the shape of circular arches that are oriented approximately along the horizontal direction in the $U-V$ plane. Moreover, in distant regions, we see variations in the velocity substructures more clearly than ever before, in particular, variations in the velocity of the Hercules stream. (abridged)

The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.

Despite the huge amount of photometric and spectroscopic efforts targetting the Galactic bulge over the past few years, its age distribution remains controversial owing to both the complexity of determining the age of individual stars and the difficult observing conditions. Taking advantage of the recent release of very deep, proper-motion-cleaned colour--magnitude diagrams (CMDs) of four low reddening windows obtained with the Hubble Space Telescope (HST), we used the CMD-fitting technique to calculate the star formation history (SFH) of the bulge at -2deg > b > -4deg along the minor axis. We find that over 80 percent of the stars formed before 8 Gyr ago, but that a significant fraction of the super-solar metallicity stars are younger than this age. Considering only the stars that are within reach of the current generation of spectrographs (i.e. V < 21), we find that 10 percent of the bulge stars are younger than 5 Gyr, while this fraction rises to 20-25 percent in the metal-rich peak. The age-metallicity relation is well parametrized by a linear fit implying an enrichment rate of dZ/dt ~ 0.005 Gyr$^{-1}$. Our metallicity distribution function accurately reproduces that observed by several spectroscopic surveys of Baade's window, with the bulk of stars having metal-content in the range [Fe/H] ~ -0.7 to ~0.6, along with a sparse tail to much lower metallicities.

Although there is consensus that metal-rich stars in the Milky Way bulge are formed via secular evolution of the thin disc, the origin of their metal-poor counterparts is still under debate. Two different origins have been invoked for metal-poor stars: they might be classical bulge stars or stars formed via internal evolution of a massive thick disc. We use N-body simulations to calculate the kinematic signature given by the difference in the mean Galactocentric radial velocity ($\Delta\rm V_{\rm GC}$) between metal-rich stars ([Fe/H] $\ge$ 0) and moderately metal-poor stars (-1.0 $\le$ [Fe/H] $<$ 0) in two models, one containing a thin disc and a small classical bulge (B/D=0.1), and the other containing a thin disc and a massive centrally concentrated thick disc. We reasonably assume that thin-disk stars in each model may be considered as a proxy of metal-rich stars. Similarly, bulge stars and thick-disc stars may be considered as a proxy of metal-poor stars. We calculate $\Delta\rm V_{\rm GC}$ at different latitudes ($b=0^\circ$, $-2^\circ$, $-4^\circ$,$-6^\circ$, $-8^\circ$ and $-10^\circ$) and longitudes ($l=0^\circ$, $\pm5^\circ$, $\pm10^\circ$ and $\pm15^\circ$) and show that the $\Delta\rm V_{\rm GC}$ trends predicted by the two models are different. We compare the predicted results with ARGOS data and APOGEE DR13 data and show that moderately metal-poor stars are well reproduced with the co-spatial stellar discs model, which has a massive thick disc. Our results give more evidence against the scenario that most of the metal-poor stars are classical bulge stars. If classical bulge stars exists, most of them probably have metallicities [Fe/H] $<$ -1 dex, and their contribution to the mass of the bulge should be a small percentage of the total bulge mass.

We study how migration affects stars of a galaxy with a thin stellar disc and thicker stellar components. The simulated galaxy has a strong bar and lasting spiral arms. We find that the amplitude of the churning (change in angular momentum) is similar for thin and thick components, and of limited amplitude, and that stars of all components can be trapped at the corotation of the bar. At the exception of those stars trapped at the corotation, we find that stars that are far from their initial guiding radius are more likely so due to blurring rather than churning effects. We compare the simulation to orbits integration with a fixed gravitational potential rotating at a constant speed. In the latter case, stars trapped at corotation are churned periodically outside and inside the corotation radius, with a zero net average. However, as the bar speed of the simulated galaxy decreases and its corotation radius increases, stars trapped at corotation for several Gyrs can be churned outwards on average. We study the location of extreme migrators (stars experimenting the largest churning) and find that extreme migrators come from regions on the leading side of the effective potential local maxima.

We show that the bulge and the disk of the Milky Way (MW) at R$\lesssim$7~kpc are well described by a unique chemical evolution and a two-phase star-formation history (SFH). We argue that the populations within this inner disk, not the entire disk, are the same, and that the outer Lindblad resonance (OLR) of the bar plays a key role in explaining this uniformity. In our model of a two-phase star formation history, the metallicity, [$\alpha$/Fe] and [$\alpha$/H] distributions, and age-metallicity relation are all compatible with the observations of both the inner disk and bulge. The dip at [Fe/H]$\sim$0 dex seen in the metallicity distributions of the bulge and inner disk reflects the quenching episode in the SFH of the inner MW at age $\sim$8 Gyr, and the common evolution of the bulge and inner disk stars. We show that at z$\le$1.5, when the MW was starting to quench, transitioning between the end of the $\alpha$-enhanced thick disk formation to the start of the thin disk, and yet was still gas rich, the gas accretion rate could not have been significant. The [$\alpha$/Fe] abundance ratio before and after this quenching phase would be different, which is not observed. The present analysis suggests that Milky Way history, and in particular at the transition from the thick to the thin disk -- the epoch of the quenching -- must have been driven by a decrease of the star formation efficiency. We argue that the decline in the intensity of gas accretion, the formation of the bar, and the quenching of the SFR at the same epoch may be causally connected thus explaining their temporal coincidence (abridged).

This letter studies the formation of azimuthal metallicity variations in the disks of spiral galaxies in the absence of initial radial metallicity gradients. Using high-resolution $N$-body simulations, we model composite stellar discs, made of kinematically cold and hot stellar populations, and study their response to spiral arm perturbations. We find that, as expected, disk populations with different kinematics respond differently to a spiral perturbation, with the tendency for dynamically cooler populations to show a larger fractional contribution to spiral arms than dynamically hotter populations. By assuming a relation between kinematics and metallicity, namely the hotter the population, the more metal-poor it is, this differential response to the spiral arm perturbations naturally leads to azimuthal variations in the mean metallicity of stars in the simulated disk. Thus, azimuthal variations in the mean metallicity of stars across a spiral galaxy are not necessarily a consequence of the reshaping, by radial migration, of an initial radial metallicity gradient. They indeed arise naturally also in stellar disks which have initially only a negative vertical metallicity gradient.

We examine the metallicity trends in the Milky Way (MW) bulge - using APOGEE DR13 data - and explore their origin by comparing two N-body models of isolated galaxies which develop a bar and a boxy/peanut (b/p) bulge. Both models have been proposed as scenarios for reconciling a disc origin of the MW bulge with a negative vertical metallicity gradient. The first model is a superposition of co-spatial, i.e. overlapping, disc populations with different scaleheights, kinematics and metallicities. In this model the thick, metal-poor, and centrally concentrated disc populations contribute significantly to the stellar mass budget in the inner galaxy. The second model is a single disc with an initial steep radial metallicity gradient, which is mapped by the bar into the b/p bulge in such a way that the vertical metallicity gradient of the MW bulge is reproduced -- as shown already in previous works in the literature. However, as we show here, the latter model does not reproduce the positive longitudinal metallicity gradient of the inner disc, nor the metal-poor innermost regions of the Bulge seen in the data. On the other hand, the model with co-spatial thin and thick disc populations reproduces all the aforementioned trends. We therefore see that it is possible to reconcile a (primarily) disc origin for the MW bulge with the observed trends in metallicity by mapping the inner thin and thick discs of the MW into a b/p. For this scenario to reproduce the observations, the $\alpha$-enhanced, metal-poor, thick disc populations must have a significant mass contribution in the inner regions -- as has been suggested for the Milky Way.

Galaxy surveys have suggested that rapid and sustained decrease in the star-formation rate, "quenching", in massive disk galaxies is frequently related to the presence of a bar. Optical and near-IR observations reveal that nearly 60% of disk galaxies in the local universe are barred, thus it is important to understand the relationship between bars and star formation in disk galaxies. Recent observational results imply that the Milky Way quenched about 9-10 Gyr ago, at the transition between the cessation of the growth of the kinematically hot, old, metal-poor thick disk and the kinematically colder, younger, and more metal-rich thin disk. Although perhaps coincidental, the quenching episode could also be related to the formation of the bar. To explore the relation between bar formation and star formation in gas rich galaxies quantitatively, we simulated gas-rich disk isolated galaxies. Our simulations include prescriptions for star formation, stellar feedback, and for regulating the multi-phase interstellar medium. We find that the action of stellar bar efficiently quenches star formation, reducing the star-formation rate by a factor of 10 in less than 1 Gyr. Analytical and self-consistent galaxy simulations with bars suggest that the action of the stellar bar increases the gas random motions within the co-rotation radius of the bar. Indeed, we detect an increase in the gas velocity dispersion up to 20-35 km/s at the end of the bar formation phase. The star-formation efficiency decreases rapidly, and in all of our models, the bar quenches the star formation in the galaxy. The star-formation efficiency is much lower in simulated barred compared to unbarred galaxies and more rapid bar formation implies more rapid quenching.

We present the chemical distribution of the Milky Way, based on 2,900$\, {\rm deg^2}$ of $u$-band photometry taken as part of the Canada-France Imaging Survey. When complete, this survey will cover 10,000$\, {\rm deg^2}$ of the Northern sky. By combing the CFHT $u$-band photometry together with SDSS and Pan-STARRS $g,r,$ and $i$, we demonstrate that we are able to measure reliably the metallicities of individual stars to $\sim 0.2$ dex, and hence additionally obtain good photometric distance estimates. This survey thus permits the measurement of metallicities and distances of the dominant main-sequence population out to approximately 30 kpc, and provides much higher number of stars at large extraplanar distances than have been available from previous surveys. We develop a non-parametric distance-metallicity decomposition algorithm and apply it to the sky at $30\deg < |b| < 70\deg$ and to the North Galactic Cap. We find that the metallicity-distance distribution is well-represented by three populations whose metallicity distributions do not vary significantly with vertical height above the disk. As traced in main-sequence stars, the stellar halo component shows a vertical density profile that is close to exponential, with a scale height of around 3 kpc. This may indicate that the inner halo was formed partly from disk stars ejected in an ancient minor merger.

The Canada-France Imaging Survey (CFIS) will map the northern high Galactic latitude sky in the $u$-band ("CFIS-u", 10,000$\, {\rm deg^2}$) and in the $r$-band ("CFIS-r", 5,000$\, {\rm deg^2}$), enabling a host of stand-alone science investigations, and providing some of the ground-based data necessary for photometric redshift determination for the Euclid mission. In this first contribution we present the $u$-band component of the survey, describe the observational strategy, and discuss some first highlight results, based on approximately one third of the final area. We show that the Galactic anticenter structure is distributed continuously along the line of sight, out to beyond 20 kpc, and possesses a metallicity distribution that is essentially identical to that of the outer disk sampled by APOGEE. This suggests that it is probably a buckled disk of old metal-rich stars, rather than a stream or a flare. We also discuss the future potential for CFIS-u in discovering star-forming dwarf galaxies around the Local Group, the characterization of the white dwarf and blue straggler population of the Milky Way, as well as its sensitivity to low-surface brightness structures in external galaxies.

This paper presents the spatially resolved star formation history (2D-SFH) of a small sample of four local mergers: the early-stage mergers IC1623, NGC6090, and the Mice, and the more advanced merger NGC2623, by analyzing IFS data from the CALIFA survey and PMAS in LArr mode. Full spectral fitting techniques are applied to the datacubes to obtain the spatially resolved mass growth histories, the time evolution of the star formation rate intensity ($\Sigma_{SFR}$), and the local specific star formation rate, over three different timescales (30 Myr, 300 Myr, and 1 Gyr). The results are compared with non-interacting Sbc--Sc galaxies, to quantify if there is an enhancement of the star formation and to trace its time scale and spatial extent. Our results for the three LIRGs (IC1623W, NGC6090, and NGC2623) show that a major phase of star formation is occurring in time scales of 10$^{7}$ yr to few 10$^{8}$ yr, with global SFR enhancements of $\sim$2--6 with respect to main-sequence star forming (MSSF) galaxies. In the two early-stage mergers IC1623W and NGC6090, which are between first pericenter passage and coalescence, the most remarkable increase of the SFR with respect to non-interacting spirals occurred in the last 30 Myr, and it is spatially extended, with enhancements of factors 2--7 both in the centres ($r <$ 0.5 half light radius, HLR), and in the disks ($r >$ 1 HLR). In the more advanced merger NGC 2623 an extended phase of star formation occurred on a longer time-scale of $\sim$1 Gyr, with a SFR enhancement of a factor $\sim$2--3 larger than the one in Sbc--Sc MSSF galaxies over the same period, probably relic of the first pericenter passage epoch. A SFR enhancement in the last 30 Myr is also present, but only in NGC2623 centre, by a factor 3. In general, the spatially resolved SFHs of the LIRG-mergers are consistent with the predictions from high spatial resolution simulations.

We report on a detailed study of the stellar populations and ionized gas properties in the merger LIRG NGC 2623, analysing optical Integral Field Spectroscopy from the CALIFA survey and PMAS LArr, multiwavelength HST imaging, and OSIRIS narrow band H$\alpha$ and [NII]$\lambda$6584 imaging. The spectra were processed with the STARLIGHT full spectral fitting code, and the results compared with those for two early-stage merger LIRGs (IC 1623 W and NGC 6090), together with CALIFA Sbc/Sc galaxies. We find that NGC 2623 went through two periods of increased star formation (SF), a first and widespread episode, traced by intermediate-age stellar populations ISP (140 Myr-1.4 Gyr), and a second one, traced by young stellar populations YSP ($<$140 Myr), which is concentrated in the central regions ($<$1.4 kpc). Our results are in agreement with the epochs of the first peri-center passage ($\sim$200 Myr ago) and coalescence ($<$100 Myr ago) predicted by dynamical models, and with high resolution merger simulations in the literature, consistent with NGC 2623 representing an evolved version of the early-stage mergers. Most ionized gas is concentrated within $<$2.8 kpc, where LINER-like ionization and high velocity dispersion ($\sim$220 km/s) are found, consistent with the previously reported outflow. As revealed by the highest resolution OSIRIS and HST data, a collection of HII regions is also present in the plane of the galaxy, which explains the mixture of ionization mechanisms in this system. It is unlikely that the outflow in NGC 2623 will escape from the galaxy, given the low SFR intensity ($\sim$0.5 M$_{\odot}$yr$^{-1}$kpc$^{-2}$), the fact that the outflow rate is 3 times lower than the current SFR, and the escape velocity in the central areas higher than the outflow velocity.

We explore trends in the morphology and line-of-sight (los) velocity of stellar populations in the inner regions of disc galaxies, using N-body simulations with both a thin (kinematically cold) and a thick (kinematically hot) disc which form a bar and boxy/peanut (b/p) bulge. The bar in the thin disc component is $\sim$50\% stronger than the thick disc bar and is more elongated, with an axis ratio almost half that of the thick disc bar. The thin disc b/p bulge has a pronounced X-shape, while the thick disc b/p is weaker with a rather boxy shape. This leads to the signature of the b/p bulge in the thick disc to be weaker and further away from the plane than in the thin disc. Regarding the kinematics, we find that the los velocity of thick disc stars in the outer parts of the b/p bulge can be \emph{larger} than that of thin disc stars, by up to 40\% and 20\% for side-on and Milky Way-like orientations of the bar respectively. This is due to the different orbits followed by thin and thick disc stars in the bar-b/p region, which are affected by the fact that: i) thin disc stars are trapped more efficiently in the bar - b/p instability and thus lose more angular momentum than their thick disc counterparts and ii) thick disc stars have large radial excursions and therefore stars from large radii with high angular momenta can be found in the bar region. We also find that the difference between the los velocities of the thin and thick disc in the b/p bulge ($\Delta v_{los}$) correlates with the initial difference between the radial velocity dispersions of the two discs ($\Delta \sigma$) . We therefore conclude that stars in the bar - b/p bulge will have considerably different morphologies and kinematics depending on the kinematic properties of the disc population they originate from.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The role of major mergers in galaxy evolution is investigated through a detailed characterization of the stellar populations, ionized gas properties, and star formation rates (SFR) in the early-stage merger LIRGs IC 1623 W and NGC 6090, by analysing optical Integral Field Spectroscopy (IFS) and high resolution HST imaging. The spectra were processed with the Starlight full spectral fitting code, and the emission lines measured in the residual spectra. The results are compared with control non-interacting spiral galaxies from the CALIFA survey. Merger-induced star formation is extended and recent, as revealed by the young ages (50-80 Myr) and high contributions to light of young stellar populations (50-90$\%$), in agreement with merger simulations in the literature. These early-stage mergers have positive central gradients of the stellar metallicity, with an average $\sim$0.6 Z$_{\odot}$. Compared to non-interacting spirals, they have lower central nebular metallicity, and flatter profiles, in agreement with the gas inflow scenario. We find that they are dominated by star formation, although shock excitation cannot be discarded in some regions, where high velocity dispersion is found (170-200 km s$^{-1}$). The average SFR in these early-stage mergers ($\sim$23-32 M$_{\odot}$ yr$^{-1}$) is enhanced with respect to main-sequence Sbc galaxies by factors of 6-9, slightly above the predictions from classical merger simulations, but still possible in about 15$\%$ of major galaxy mergers, where U/LIRGs belong.

In this work, two new axisymmetric models for the Galactic mass distribution are presented. Motivated by recent results, these two models include the contribution of a stellar thin disc and of a thick disc, as massive as the thin counterpart but with a shorter scale-length. Both models satisfy a number of observational constraints: stellar densities at the solar vicinity, thin and thick disc scale lengths and heights, rotation curve(s), and the absolute value of the perpendicular force Kz as a function of distance to the Galactic centre. We numerically integrate into these new models the motion of all Galactic globular clusters for which distances, proper motions, and radial velocities are available, and the orbits of about one thousand stars in the solar vicinity. The retrieved orbital characteristics are compared to those obtained by integrating the clusters and stellar orbits in pure thin disc models. We find that, due to the possible presence of a thick disc, the computed orbital parameters of disc stars can vary by as much as 30-40%. We also show that the systematic uncertainties that affect the rotation curve still plague computed orbital parameters of globular clusters by similar amounts.

The $\Lambda$CDM cosmological scenario predicts that our Galaxy should contain hundreds of stellar streams at the solar vicinity, fossil relics of the merging history of the Milky Way and more generally of the hierarchical growth of galaxies. Because of the mixing time scales in the inner Galaxy, it has been claimed that these streams should be difficult to detect in configuration space but can still be identifiable in kinematic-related spaces like the energy/angular momenta spaces, E-Lz and Lperp-Lz, or spaces of orbital/velocity parameters. By means of high-resolution, dissipationless N-body simulations, containing between 25$\times10^6$ and 35$\times10^6$ particles, we model the accretion of a series of up to four 1:10 mass ratio satellites then up to eight 1:100 satellites and we search systematically for the signature of these accretions in these spaces. In all spaces considered (1) each satellite gives origin to several independent overdensities; (2) overdensities of multiple satellites overlap; (3) satellites of different masses can produce similar substructures; (4) the overlap between the in-situ and the accreted population is considerable everywhere; (5) in-situ stars also form substructures in response to the satellite(s) accretion. These points are valid even if the search is restricted to kinematically-selected halo stars only. As we are now entering the 'Gaia era', our results warn that an extreme caution must be employed before interpreting overdensities in any of those spaces as evidence of relics of accreted satellites. Reconstructing the accretion history of our Galaxy will require a substantial amount of accurate spectroscopic data, that, complemented by the kinematic information, will possibly allow us to (chemically) identify accreted streams and measure their orbital properties. (abridged)

Galaxy mergers are believed to trigger strong starbursts. This is well assessed by observations in the local Universe. However the efficiency of this mechanism has poorly been tested so far for high redshift, actively star forming, galaxies. We present a suite of pc-resolution hydrodynamical numerical simulations to compare the star formation process along a merging sequence of high and low z galaxies, by varying the gas mass fraction between the two models. We show that, for the same orbit, high-redshift gas-rich mergers are less efficient than low-redshift ones at producing starbursts: the star formation rate excess induced by the merger and its duration are both around 10 times lower than in the low gas fraction case. The mechanisms that account for the star formation triggering at low redshift - the increased compressive turbulence, gas fragmentation, and central gas inflows - are only mildly, if not at all, enhanced for high gas fraction galaxy encounters. Furthermore, we show that the strong stellar feedback from the initially high star formation rate in high redshift galaxies does not prevent an increase of the star formation during the merger. Our results are consistent with the observed increase of the number of major mergers with increasing redshift being faster than the respective increase in the number of starburst galaxies.

The determination of the age of the bulge has led to two contradictory results. On the one side, the color-magnitude diagrams in different bulge fields seem to indicate a uniformly old ($>$10 Gyr) population. On the other side, individual ages derived from dwarfs observed through microlensing events seem to indicate a large spread, from $\sim$ 2 to $\sim$ 13 Gyr. Because the bulge is now recognised as being mainly a boxy peanut-shaped bar, it is suggested that disk stars are one of its main constituents, and therefore also stars with ages significantly younger than 10 Gyr. Other arguments as well point to the fact that the bulge cannot be exclusively old, and in particular cannot be a burst population, as it is usually expected if the bulge was the fossil remnant of a merger phase in the early Galaxy. In the present study, we show that given the range of metallicities observed in the bulge, a uniformly old population would be reflected into a significant spread in color at the turn-off which is not observed. Inversely, we demonstrate that the correlation between age and metallicity expected to hold for the inner disk would conspire to form a color-magnitude diagram with a remarkably small spread in color, thus mimicking the color-magnitude diagram of a uniformly old population. If stars younger than 10 Gyr are part of the bulge, as must be the case if the bulge has been mainly formed through dynamical instabilities in the disk, then a very small spread at the turn-off is expected, as seen in the observations.

Quenching, the cessation of star formation, is one of the most significant events in the life cycle of galaxies. We show here the first evidence that the Milky Way experienced a generalised quenching of its star formation at the end of its thick disk formation $\sim$9 Gyr ago. Elemental abundances of stars studied as part of the APOGEE survey reveal indeed that in less than $\sim$2 Gyr the star formation rate in our Galaxy dropped by an order-of-magnitude. Because of the tight correlation between age and alpha abundance, this event reflects in the dearth of stars along the inner disk sequence in the [Fe/H]-[$\alpha$/Fe] plane. Before this phase, which lasted about 1.5 Gyr, the Milky Way was actively forming stars. Afterwards, the star formation resumed at a much lower level to form the thin disk. These events are very well matched by the latest observation of MW-type progenitors at high redshifts. In late type galaxies, quenching is believed to be related to a long and secular exhaustion of gas. In our Galaxy, it occurred on a much shorter time scale, while the chemical continuity before and after the quenching indicates that it was not due to the exhaustion of the gas. While quenching is generally associated with spheroids, our results show that it also occurs in galaxies like the Milky Way, possibly when they are undergoing a morphological transition from thick to thin disks. Given the demographics of late type galaxies in the local universe, in which classical bulges are rare, we suggest further that this may hold true generally in galaxies with mass lower than or approximately $M^*$, where quenching could be directly a consequence of thick disk formation. We emphasize that the quenching phase in the Milky Way could be contemporaneous with, and related to, the formation of the bar. We sketch a scenario on how a strong bar may inhibit star formation.

MSE is an 11.25m aperture observatory with a 1.5 square degree field of view that will be fully dedicated to multi-object spectroscopy. More than 3200 fibres will feed spectrographs operating at low (R ~ 2000 - 3500) and moderate (R ~ 6000) spectral resolution, and approximately 1000 fibers will feed spectrographs operating at high (R ~ 40000) resolution. MSE is designed to enable transformational science in areas as diverse as tomographic mapping of the interstellar and intergalactic media; the in-situ chemical tagging of thick disk and halo stars; connecting galaxies to their large scale structure; measuring the mass functions of cold dark matter sub-halos in galaxy and cluster-scale hosts; reverberation mapping of supermassive black holes in quasars; next generation cosmological surveys using redshift space distortions and peculiar velocities. MSE is an essential follow-up facility to current and next generations of multi-wavelength imaging surveys, including LSST, Gaia, Euclid, WFIRST, PLATO, and the SKA, and is designed to complement and go beyond the science goals of other planned and current spectroscopic capabilities like VISTA/4MOST, WHT/WEAVE, AAT/HERMES and Subaru/PFS. It is an ideal feeder facility for E-ELT, TMT and GMT, and provides the missing link between wide field imaging and small field precision astronomy. MSE is optimized for high throughput, high signal-to-noise observations of the faintest sources in the Universe with high quality calibration and stability being ensured through the dedicated operational mode of the observatory. (abridged)

The Galactic bulge, that is the prominent out-of-plane over-density present in the inner few kiloparsecs of the Galaxy, is a complex structure, as the morphology, kinematics, chemistry and ages of its stars indicate. To understand the nature of its main components -- those at [Fe/H] >~ -1 dex -- it is necessary to make an inventory of the stellar populations of the Galactic disc(s), and of their borders : the chemistry of the disc at the solar vicinity, well known from detailed studies of stars over many years, is not representative of the whole disc. This finding, together with the recent revisions of the mass and sizes of the thin and thick discs, constitutes a major step in understanding the bulge complexity. N-body models of a boxy/peanut-shaped bulge formed from a thin disc through the intermediary of a bar have been successful in interpreting a number of global properties of the Galactic bulge, but they fail in reproducing the detailed chemo-kinematic relations satisfied by its components and their morphology. It is only by adding the thick disc to the picture that we can understand the nature of the Galactic bulge.

(Abridged) By using an N-body simulation of a bulge that was formed via a bar instability mechanism, we analyse the imprints of the initial (i.e. before bar formation) location of stars on the bulge kinematics, in particular on the heliocentric radial velocity distribution of bulge stars. Four different latitudes were considered: $b=-4^\circ$, $-6^\circ$, $-8^\circ$, and $-10^\circ$, along the bulge minor axis as well as outside it, at $l=\pm5^\circ$ and $l=\pm10^\circ$. The bulge X-shaped structure comprises stars that formed in the disk at different locations. Stars formed in the outer disk, beyond the end of the bar, which are part of the boxy peanut-bulge structure may show peaks in the velocity distributions at positive and negative heliocentric radial velocities with high absolute values that can be larger than 100 $\rm km$ $\rm s^{-1}$, depending on the observed direction. In some cases the structure of the velocity field is more complex and several peaks are observed. Stars formed in the inner disk, the most numerous, contribute predominantly to the X-shaped structure and present different kinematic characteristics. Our results may enable us to interpret the cold high-velocity peak observed in the APOGEE commissioning data, as well as the excess of high-velocity stars in the near and far arms of the X-shaped structure at $l$=$0^\circ$ and $b$=$-6^\circ$. When compared with real data, the kinematic picture becomes more complex due to the possible presence in the observed samples of classical bulge and/or thick disk stars. Overall, our results point to the existence of complex patterns and structures in the bulge velocity fields, which are generated by the bar. This suggests that caution should be used when interpreting the bulge kinematics: the presence of substructures, peaks and clumps in the velocity fields is not necessarily a sign of past accretion events.

Radial stellar migration in galactic discs has received much attention in studies of galactic dynamics and chemical evolution, but remains a dynamical phenomenon that needs to be fully quantified. In this work, using a Tree-SPH simulation of an Sb-type disc galaxy, we quantify the effects of blurring (epicyclic excursions) and churning (change of guiding radius). We quantify migration (either blurring or churning) both in terms of flux (the number of migrators passing at a given radius), and by estimating the population of migrators at a given radius at the end of the simulation compared to non-migrators, but also by giving the distance over which the migration is effective at all radii. We confirm that the corotation of the bar is the main source of migrators by churning in a bar-dominated galaxy, its intensity being directly linked to the episode of a strong bar, in the first 1-3 Gyr of the simulation. We show that within the outer Lindblad resonance (OLR), migration is strongly dominated by churning, while blurring gains progressively more importance towards the outer disc and at later times. Most importantly, we show that the OLR limits the exchange of angular momentum, separating the disc in two distinct parts with minimal or null exchange, except in the transition zone, which is delimited by the position of the OLR at the epoch of the formation of the bar, and at the final epoch. We discuss the consequences of these findings for our understanding of the structure of the Milky Way disc. Because the Sun is situated slightly outside the OLR, we suggest that the solar vicinity may have experienced very limited churning from the inner disc.

We analyse the chemical properties of a set of solar vicinity stars, and show that the small dispersion in abundances of \alpha-elements at all ages provides evidence that the SFH has been uniform throughout the thick disk. In the context of long time scale infall models, we suggest that this result points either to a limited dependence of the gas accretion on the Galactic radius in the inner disk (R<10 kpc), or to a decoupling of the accretion history and star formation history due to other processes governing the ISM in the early disk, suggesting that infall cannot be a determining parameter of the chemical evolution at these epochs. We argue however that these results and other recent observational constraints -- namely the lack of radial metallicity gradient and the non-evolving scale length of the thick disk -- are better explained if the early disk is viewed as a pre-assembled gaseous system, with most of the gas settled before significant star formation took place -- formally the equivalent of a closed-box model. In any case, these results point to a weak, or non-existent inside-out formation history in the thick disk, or in the first 3-5 Gyr of the formation of the Galaxy. We argue however that the growing importance of an external disk whose chemical properties are distinct from those of the inner disk would give the impression of an inside-out growth process when seen through snapshots at different epochs. However, the progressive, continuous process usually invoked may not have actually existed in the Milky Way.

The apparent correlation between the specific star formation rate (sSFR) and total stellar mass (M_star) of galaxies is a fundamental relationship indicating how they formed their stellar populations. To attempt to understand this relation, we hypothesize that the relation and its evolution is regulated by the increase in the stellar and gas mass surface density in galaxies with redshift, which is itself governed by the angular momentum of the accreted gas, the amount of available gas, and by self-regulation of star formation. With our model, we can reproduce the specific SFR-M_star relations at z~1-2 by assuming gas fractions and gas mass surface densities similar to those observed for z=1-2 galaxies. We further argue that it is the increasing angular momentum with cosmic time that causes a decrease in the surface density of accreted gas. The gas mass surface densities in galaxies are controlled by the centrifugal support (i.e., angular momentum), and the sSFR is predicted to increase as, sSFR(z)=(1+z)^3/t_H0, as observed (where t_H0 is the Hubble time and no free parameters are necessary). At z>~2, we argue that star formation is self-regulated by high pressures generated by the intense star formation itself. The star formation intensity must be high enough to either balance the hydrostatic pressure (a rather extreme assumption) or to generate high turbulent pressure in the molecular medium which maintains galaxies near the line of instability (i.e. Toomre Q~1). The most important factor is the increase in stellar and gas mass surface density with redshift, which allows distant galaxies to maintain high levels of sSFR. Without a strong feedback from massive stars, such galaxies would likely reach very high sSFR levels, have high star formation efficiencies, and because strong feedback drives outflows, ultimately have an excess of stellar baryons (abridged).

We develop a chemical evolution model in order to study the star formation history of the Milky Way. Our model assumes that the Milky Way is formed from a closed box-like system in the inner regions, while the outer parts of the disc experience some accretion. Unlike the usual procedure, we do not fix the star formation prescription (e.g. Kennicutt law) in order to reproduce the chemical abundance trends. Instead, we fit the abundance trends with age in order to recover the star formation history of the Galaxy. Our method enables one to recover with unprecedented accuracy the star formation history of the Milky Way in the first Gyrs, in both the inner (R<7-8kpc) and outer (R>9-10kpc) discs as sampled in the solar vicinity. We show that, in the inner disc, half of the stellar mass formed during the thick disc phase, in the first 4-5 Gyr. This phase was followed by a significant dip in the star formation activity (at 8-9 Gyr) and a period of roughly constant lower level star formation for the remaining 8 Gyr. The thick disc phase has produced as many metals in 4 Gyr as the thin disc in the remaining 8 Gyr. Our results suggest that a closed box model is able to fit all the available constraints in the inner disc. A closed box system is qualitatively equivalent to a regime where the accretion rate, at high redshift, maintains a high gas fraction in the inner disc. In such conditions, the SFR is mainly governed by the high turbulence of the ISM. By z~1 it is possible that most of the accretion takes place in the outer disc, while the star formation activity in the inner disc is mostly sustained by the gas not consumed during the thick disc phase, and the continuous ejecta from earlier generations of stars. The outer disc follows a star formation history very similar to that of the inner disc, although initiated at z~2, about 2 Gyr before the onset of the thin disc formation in the inner disc.

By analyzing a N-body simulation of a bulge formed simply via a bar instability mechanism operating on a kinematically cold stellar disk, and by comparing the results of this analysis with the structural and kinematic properties of the main stellar populations of the Milky Way bulge, we conclude that the bulge of our Galaxy is not a pure stellar bar formed from a pre-existing thin stellar disk, as some studies have recently suggested. On the basis of several arguments emphasized in this paper, we propose that the bulge population which, in the Milky Way, is observed not to be part of the peanut structure corresponds to the old galactic thick disk, thus implying that the Milky Way is a pure thin+thick disk galaxy, with only a possible limited contribution of a classical bulge.

We compare the star-formation history and dynamics of the Milky Way (MW) with the properties of distant disk galaxies. During the first ~4 Gyr of its evolution, the MW formed stars with a high star-formation intensity (SFI), Sigma_SFR~0.6 Msun/yr/kpc2 and as a result, generated outflows and high turbulence in its interstellar medium. This intense phase of star formation corresponds to the formation of the thick disk. The formation of the thick disk is a crucial phase which enables the MW to have formed approximately half of its total stellar mass by z~1 which is similar to "MW progenitor galaxies" selected by abundance matching. This agreement suggests that the formation of the thick disk may be a generic evolutionary phase in disk galaxies. Using a simple energy injection-kinetic energy relationship between the 1-D velocity dispersion and SFI, we can reproduce the average perpendicular dispersion in stellar velocities of the MW with age. This relationship, its inferred evolution, and required efficiency are consistent with observations of galaxies from z~0-3. The high turbulence generated by intense star formation naturally resulted in a thick disk, a chemically well-mixed ISM, and is the mechanism that links the evolution of MW to the observed characteristics of distant disk galaxies.

By means of idealized, dissipationless N-body simulations which follow the formation and subsequent buckling of a stellar bar, we study the characteristics of boxy/peanut-shaped bulges and compare them with the properties of the stellar populations in the Milky Way bulge. The main results of our modeling, valid for the general family of boxy/peanut shaped bulges, are the following: (i) because of the spatial redistribution in the disk initiated at the epoch of bar formation, stars from the innermost regions to the outer Lindblad resonance of the stellar bar are mapped into a boxy bulge; (ii) the contribution of stars to the local bulge density depends on their birth radius: stars born in the innermost disk tend to dominate the innermost regions of the boxy bulge, while stars originating closer to the OLR are preferentially found in the outer regions of the boxy/peanut structure; (iii) stellar birth radii are imprinted in the bulge kinematics, the larger the birth radii of stars ending up in the bulge, the greater their rotational support and the higher their line-of- sight velocity dispersions (but note that this last trend depends on the bar viewing angle); (iv) the higher the classical bulge-over-disk ratio, the larger its fractional contribution of stars at large vertical distance from the galaxy mid-plane. (ABRIDGED) On the basis of their chemical and kinematic characteristics, the results of our modeling suggests that the populations A, B and C, as defined by the ARGOS survey, can be associated, respectively, with the inner thin disk, to the young thick and to the old thick disk, following the nomenclature recently suggested for stars in the solar neighborhood by Haywood et al. (2013).

Constraints on the Galactic bulge/bar structure and formation history from stellar kinematics and metallicities mainly come from relatively high-latitude fields (|b|>4) where a complex mix of stellar population is seen. We aim here to constrain the formation history of the Galactic bar by studying the radial velocity and metallicity distributions of stars in-situ (|b|<1). We observed red clump stars in four fields along the bar's major axis (l=10,6,-6 and b=0 plus a field at l=0,b=1) with low-resolution spectroscopy from VLT/FLAMES, observing around the CaII triplet. We developed robust methods for extracting radial velocity and metallicity estimates from these low signal-to-noise spectra. We derived distance probability distributions using Bayesian methods rigorously handling the extinction law. We present radial velocities and metallicity distributions, as well as radial velocity trends with distance. We observe an increase in the radial velocity dispersion near the Galactic plane. We detect the streaming motion of the stars induced by the bar in fields at l=+/-6, the highest velocity components of this bar stream being metal-rich ([Fe/H]~0.2 dex). Our data is consistent with a bar inclined at 26+/-3 from the Sun-Galactic centre line. We observe a significant fraction of metal-poor stars, in particular in the field at l=0,b=1. We confirm the flattening of the metallicity gradient along the minor axis when getting closer to the plane, with a hint that it could actually be inverted. Our stellar kinematics corresponds to the expected behaviour of a bar issued from the secular evolution of the Galactic disc. The mix of several populations, seen further away from the plane, is also seen in the bar in-situ since our metallicity distributions highlight a different spatial distribution between metal-poor and metal-rich stars, the more metal-poor stars being more centrally concentrated.

We report the first robust measurement of the Milky Way star formation history using the imprint left on chemical abundances of long-lived stars. The formation of the Galactic thick disc occurs during an intense star formation phase between 9.0 (z~1.5) and 12.5 Gyr (z~4.5) ago and is followed by a dip (at z~1.1) lasting about 1 Gyr. Our results imply that the thick disc is as massive as the Milky Way's thin disc, suggesting a fundamental role of this component in the genesis of our Galaxy, something that had been largely unrecognized. This new picture implies that huge quantities of gas necessary to feed the building of the thick disc must have been present at these epochs, in contradiction with the long-term infall assumed by chemical evolution models in the last two decades. These results allow us to fit the Milky Way within the emerging features of the evolution of disc galaxies in the early Universe.

We analyze high quality abundances data of solar neighborhood stars and show that there are two distinct regimes of [alpha/Fe] versus age which we identify as the epochs of the thick and thin disk formation. A tight correlation between metallicity and [alpha/Fe] versus age is clearly identifiable on thick disk stars, implying that this population formed from a well mixed ISM, over a time scale of 4-5 Gyr. Thick disk stars vertical velocity dispersion correlate with age, with the youngest objects having as small scale heights as those of thin disk stars. A natural consequence of these two results is that a vertical metallicity gradient is expected in this population. We suggest that the thick disk set the initial conditions for the formation of the inner thin disk. This provides also an explanation of the apparent coincidence between the step in metallicity at 7-10 kpc in the thin disk and the confinment of the thick disk at about R<10 kpc. We suggest that the outer thin disk developped outside the influence of the thick disk, but also that the high alpha-enrichment of the outer regions may originate from a primordial pollution by the gas expelled from the thick disk. Local metal-poor thin disk stars, whose properties are best explained by an origin in the outer disk, are shown to be as old as the youngest thick disk (9-10 Gyr), implying that the outer thin disk started to form while the thick disk formation was still on-going in the inner Galaxy. We point out that, given the tight age-abundance relations in the thick disk, an inside-out process would give rise to a radial gradient in abundances in this population which is not observed. Finally, we argue that the data discussed here leave little room for radial migration, either to have contaminated the solar vicinity, or, to have redistributed stars in significant proportion across the solar annulus.

We explore a second order Hamiltonian vertical resonance model for X-shaped or peanut-shaped galactic bulges. The X-shape is caused by the 2:1 vertical Lindblad resonance with the bar, with two vertical oscillation periods per orbital period in the bar frame. We examine N-body simulations and find that due to the bar slowing down and disk thickening during bar buckling, the resonance and associated peanut-shape moves outward. The peanut-shape is consistent with the location of the vertical resonance, independent of whether the bar buckled or not. We estimate the resonance width from the potential m=4 Fourier component and find that the resonance is narrow, affecting orbits over a narrow range in the angular momentum distribution, dL/L ~ 0.05. As the resonance moves outward, stars originally in the mid plane are forced out of the mid plane into orbits just within the resonance separatrix. The height of the separatrix orbits, estimated from the Hamiltonian model, is approximately consistent with the peanut-shape height. The X-shape is comprised of stars in the vicinity of the resonance separatrix. The velocity distributions from the simulations illustrate that low inclination orbits are depleted within resonance. Within resonance, the vertical velocity distribution is broad, consistent with resonant heating caused by the passage of the resonance through the disk. In the Milky Way bulge we relate the azimuthally averaged mid-plane mass density near the vertical resonance to the rotation curve and bar pattern speed. At an estimated vertical resonance galactocentric radius of ~1.3 kpc, we confirm a mid-plane density of ~5x10^8 Msol/kpc^3, consistent with recently estimated mass distributions. We find that the rotation curve, bar pattern speed, 2:1 vertical resonance location, X-shape tips, and mid-plane mass density, are all self-consistent in the Milky Way galaxy bulge.

(abridged) We have analyzed the properties of the rest-frame optical emission lines of a sample of 53 intensely star forming galaxies at z=1.3 to 2.7 observed with SINFONI on the ESO-VLT. We find large velocity dispersions in the lines, sigma=30-250 km/s. Our data agree well with simulations where we applied beam-smearing and assumed a scaling relation of the form: velocity dispersion is proportional to the square root of the star-formation intensity (star-formation rate per unit area). We conclude that the dispersions are primarily driven by star formation. To explain the high surface brightness and optical line ratios, high thermal pressures in the warm ionized medium, WIM, are required (log P/k (K/cm^3)>~6-7). Such thermal pressures in the WIM are similar to those observed in nearby starburst galaxies, but occur over much larger physical scales. Moreover, the relatively low ionization parameters necessary to fit the high surface brightnesses and optical line ratios suggest that the gas is not only directly associated with regions of star formation, but is wide spread throughout the general ISM. Thus the optical emission line gas is a tracer of the large scale dynamics of the bulk of the ISM. We present a simple model for the energy input from young stars in an accreting galaxy, to argue that the intense star-formation is supporting high turbulent pressure, which roughly balances the gravitational pressure and thus enables distant gas accreting disks to maintain a Toomre disk instability parameter Q~1. For a star formation efficiency of 3%, only 5-15% of the mechanical energy from young stars that is deposited in the ISM is needed to support the level of turbulence required for maintaining this balance. Since this balance is maintained by energy injected into the ISM by the young stars themselves, this suggests that star formation in high redshift galaxies is self-regulating.

In this paper we report a study (see Mastrobuono-Battisti et al., 2012) about the formation and characteristics of the tidal tails around Palomar 5 along its orbit in the Milky Way potential, by means of direct N-body simulations and simplified numerical models. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested.

By means of N-body simulations, we show that radial migration in galaxy disks, induced by bar and spiral arms, leads to significant azimuthal variations in the metallicity distribution of old stars at a given distance from the galaxy center. Metals do not show an axisymmetric distribution during phases of strong migration. Azimuthal variations are visible during the whole phase of strong bar phase, and tend to disappear as the effect of radial migration diminishes, together with a reduction in the bar strength. These results suggest that the presence of inhomogeneities in the metallicity distribution of old stars in a galaxy disk can be a probe of ongoing strong migration. Such signatures may be detected in the Milky Way by Gaia (and complementary spectroscopic data), as well as in external galaxies, by IFU surveys like CALIFA and ATLAS3D. Mixing - defined as the tendency toward a homogeneous, azimuthally symmetric, stellar distribution in the disk - and migration turns out to be two distinct processes, the effects of mixing starting to be visible when strong migration is over.

We investigate the evolution of galactic disks in N-body Tree-SPH simulations. We find that disks, initially truncated at three scale-lengths, can triple their radial extent, solely driven by secular evolution. Both Type I (single exponential) and Type II (down-turning) observed disk surface-brightness profiles can be explained by our findings. We relate these results to the strong angular momentum outward transfer, resulting from torques and radial migration associated with multiple patterns, such as central bars and spiral waves of different multiplicity. We show that even for stars ending up on cold orbits, the changes in angular momentum exhibit complex structure as a function of radius, unlike the expected effect of transient spirals alone. Focussing on one of our models, we find evidence for non-linear coupling among m=1, 2, 3 and 4 density waves, where m is the pattern multiplicity. We suggest that the naturally occurring larger resonance widths at galactic radii beyond four scale-lengths may have profound consequences on the formation and location of breaks in disk density profiles, provided spirals are present at such large distances. We also consider the effect of gas inflow and show that when in-plane smooth gas accretion of ~5 M_sun/yr is included, the outer disks become more unstable, leading to a strong increase in the stellar velocity dispersion. This, in turn, causes the formation of a Type III (up-turning) profile in the old stellar population. We propose that observations of Type III surface brightness profiles, combined with an up-turn in the stellar velocity dispersions beyond the disk break, could be a signature of ongoing gas-accretion. The results of this study suggest that disk outskirts comprised of stars migrated from the inner disk would have relatively large radial velocity dispersions, and significant thickness when seen edge-on. [Abridged]

By means of direct N-body simulations and simplified numerical models, we study the formation and characteristics of the tidal tails around Palomar 5, along its orbit in the Milky Way potential. Unlike previous findings, we are able to reproduce the substructures observed in the stellar streams of this cluster, without including any lumpiness in the dark matter halo. We show that overdensities similar to those observed in Palomar 5 can be reproduced by the epicyclic motion of stars along its tails, i.e. a simple local accumulation of orbits of stars that escaped from the cluster with very similar positions and velocities. This process is able to form stellar clumps at distances of several kiloparsecs from the cluster, so it is not a phenomenon confined to the inner part of Palomar 5's tails, as previously suggested. Our models can reproduce the density contrast between the clumps and the surrounding tails found in the observed streams, without including any lumpiness in the dark halo, suggesting new upper limits on its granularity.

By means of N-body simulations we investigate the impact of minor mergers on the angular momentum and dynamical properties of the merger remnant. Our simulations cover a range of initial orbital characteristics and gas-to-stellar mass fractions (from 0 to 20%), and include star formation and supernova feedback. We confirm and extend previous results by showing that the specific angular momentum of the stellar component always decreases independently of the orbital parameters or morphology of the satellite, and that the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy of the orbits. However, the decrease affects only the old stellar population, and not the new population formed from gas during the merging process. This means that the merging process induces an increasing difference in the rotational support of the old and young stellar components, with the old one lagging with respect to the new. Even if our models are not intended specifically to reproduce the Milky Way and its accretion history, we find that, under certain conditions, the modeled rotational lag found is compatible with that observed in the Milky Way disk, thus indicating that minor mergers can be a viable way to produce it. The lag can increase with the vertical distance from the disk midplane, but only if the satellite is accreted along a direct orbit, and in all cases the main contribution to the lag comes from stars originally in the primary disk rather than from stars in the satellite galaxy. We also discuss the possibility of creating counter-rotating stars in the remnant disk, their fraction as a function of the vertical distance from the galaxy midplane, and the cumulative effect of multiple mergers on their creation.

An important property of star-forming galaxies at z~1-2 is the high local star-formation intensities they maintain over tens of kiloparsecs at levels that are only observed in the nearby Universe in the most powerful nuclear starbursts. To investigate how these high star-formation intensities affect the warm ionized medium, we present an analysis of the average spectra of about 50 such galaxies at z~1.2-2.6 and of subsamples selected according to their local and global star-formation intensity. Stacking allows us to probe relatively weak lines like [SII]\lambda \lambda 6716,6731 and [OI]\lambda 6300, which are tracers of the conditions of the ISM and are undetectable in most individual targets. We find higher gas densities (hence pressures) in intensely star-forming regions compared to fainter diffuse gas and, overall, values that are comparable to starburst regions and the diffuse ISM in nearby galaxies. By modeling the H\alpha\ surface brightnesses and [SII]/H\alpha\ line ratios with the Cloudy photoionization code, we find that our galaxies continue trends observed in local galaxies, where gas pressures scale with star-formation intensity. We discuss these results in the context of models of self-regulated star formation, where star formation determines the average thermal and turbulent pressure in the ISM, which in turn determines the rate at which stars can form, finding good agreement with our data. We also confirm the detection of broad, faint lines underlying H\alpha\ and [NII], which have previously been considered evidence of either outflows or active galactic nuclei. Finding that the broad component is only significantly detected in stacks with the highest average local and global star-formation intensities strongly supports the outflow interpretation, and further emphasizes the importance of star-formation feedback and self-regulation in the early Universe.

In previous studies, it has been shown that the large line widths observed in high surface brightness H-alpha emitters at low and high redshifts are likely due to the mechanical energy injected by intense star formation. Here we discuss the possibility that the high surface brightnesses observed are not due to star formation, but due to cosmological gas accretion. We assume that all of the accretion energy is dissipated as shocks from the accreting gas. We show that in order to explain the high surface brightnesses both the mass accretion rate and energy would have to be much higher than expected from simulations or from equating the star formation with the accretion rate. We also investigate scaling relations between the surface brightness expected from accretion and for star formation through mechanical heating and photo-ionization, trying to identify a regime where such accretion may become evident in galaxies. Unfortunately, the surface brightness necessary to detect the gas in optical line emission is about an order of magnitude lower than what has currently been achieved with near-infrared observations of distant galaxies.

By means of N-body/SPH simulations we investigate the morphological properties of thick disks formed through minor mergers. We show that the vertical surface density profile of the post-merger thick disk follows a sech function and has an excess in the regions far from the disk mid-plane (z>2kpc). This stellar excess also follows a sech function with a larger scale height than the main thick disk component, and it is usually dominated by stars from the primary galaxy. Stars in the excess have a rotational velocity lower than that of stars in the thick disk, and they may thus be confused with stars in the inner galactic halo. The thick disk scale height increases with radius and the rate of its increase is smaller for more gas rich primary galaxies. On the contrary, the scale height of the stellar excess is independent of both radius and gas fraction. We also find that the post-merger thick disk has a radial scale length which is 10-50% larger than that of the thin disk. Two consecutive mergers have basically the same effect on heating the stellar disk as a single merger of the same total mass. To investigate how thick disks produced through secular processes may differ from those produced by minor mergers, we also simulated gravitationally unstable gas-rich disks. These disks do not produce either a stellar excess or a ratio of thick to thin disk scale lengths greater than one. Our simulation results are consistent with observations of the ratio of thick to thin disk scale lengths of the Milky Way and nearby galaxies, and with the Toomre diagram of the Milky Way. We conclude that minor mergers are a viable mechanism for the creation of galactic thick disks and investigating stars at several kpc above the mid-plane of the Milky Way and other galaxies may provide a quantitative method for studying the minor merger history of galaxies.

It has been known since many decades that galaxy interactions can induce star formation (hereafter SF) enhancements and that one of the driving mechanisms of this enhancement is related to gas inflows into the central galaxy regions, induced by asymmetries in the stellar component, like bars. In the last years many evidences have been accumulating, showing that interacting pairs have central gas-phase metallicities lower than those of field galaxies, by {\sim} 0.2-0.3 dex on average. These diluted ISM metallicities have been explained as the result of inflows of metal-poor gas from the outer disk to the galaxy central regions. A number of questions arises: What's the timing and the duration of this dilution? How and when does the SF induced by the gas inflow enrich the circumnuclear gas with re-processed material? Is there any correlation between the timing and strength of the dilution and the timing and intensity of the SF? By means of Tree-SPH simulations of galaxy major interactions, we have studied the effect that gas inflows have on the ISM dilution, and the effect that the induced SF has, subsequently, in re-enriching the nuclear gas. In this contribution, we present the main results of this study.

We have recently identified a new radial migration mechanism resulting from the overlap of spiral and bar resonances in galactic disks. Here we confirm the efficiency of this mechanism in fully self-consistent, Tree-SPH simulations, as well as high-resolution pure N-body simulations. In all barred cases we clearly identify the effect of spiral-bar resonance overlap by measuring a bimodality in the changes of angular momentum in the disk, dL, whose maxima are near the bar's corotation and outer Lindblad resonance. This contrasts with the smooth distribution of dL for a simulation with no stable bar present, where strong radial migration is induced by multiple spirals. The presence of a disk gaseous component appears to increase the rate of angular momentum exchange by about 20%. The efficiency of this mechanism is such that galactic stellar disks can extend to over 10 scale-lengths within 1-3 Gyr in both Milky Way size and low-mass galaxies (circular velocity ~100 km/s). We also show that metallicity gradients can flatten in less than 1 Gyr rendering mixing in barred galaxies an order of magnitude more efficient than previously thought.

We study the distribution of orbital eccentricities of stars in thick disks generated by the heating of a pre-existing thin stellar disk through a minor merger (mass ratio 1:10), using N-body/SPH numerical simulations of interactions that span a range of gas fractions in the primary disk and initial orbital configurations. The resulting eccentricity distributions have an approximately triangular shape, with a peak at 0.2-0.35, and a relatively smooth decline towards higher values. Stars originally in the satellite galaxy tend to have higher eccentricities (on average from e = 0.45 to e = 0.75), which is in general agreement with the models of Sales and collaborators, although in detail we find fewer stars with extreme values and no evidence of their secondary peak around e = 0.8. The absence of this high-eccentricity feature results in a distribution that qualitatively matches the observations. Moreover, the increase in the orbital eccentricities of stars in the solar neighborhood with vertical distance from the Galactic mid-plane recently found by Diericxk and collaborators can be qualitatively reproduced by our models, but only if the satellite is accreted onto a direct orbit. We thus speculate that if minor mergers were the dominant means of formating the Milky Way thick disk, the primary mechanism should be merging with satellite(s) on direct orbits.

We present simulations of the 21-cm signal during the epoch of reionization. We focus on properly modeling the absorption regime in the presence of inhomogeneous Wouthuysen-Field effect and X-ray heating. We ran radiative transfer simulations for three bands in the source spectrum (Lyman, UV, and X-ray) to fully account for these processes. We find that the brightness temperature fluctuation of the 21 cm signal has an amplitude greater than 100 mK during the early reionization, up to 10 times greater than the typical amplitude of a few 10 mK obtained during the later emission phase. More importantly, we find that even a rather high contribution from QSO-like sources only damps the absorption regime without erasing it. Heating the IGM with X-ray takes time. Our results show that observations of the early reionization will probably benefit from a higher signal-to-noise value than during later stages. After analyzing the statistical properties of the signal (power spectrum and PDF) we find three diagnostics to constrain the level of X-ray, hence the nature of the first sources.

The "Extreme starbursts in the local universe" workshop was held at the Insituto de Astrofisica de Andalucia in Granada, Spain on 21-25 June 2010. Bearing in mind the advent of a new generation of facilities such as JWST, Herschel, ALMA, eVLA and eMerlin, the aim of the workshop was to bring together observers and theorists to review the latest results. The purpose of the workshop was to address the following issues: what are the main modes of triggering extreme starbursts in the local Universe? How efficiently are stars formed in extreme starbursts? What are the star formation histories of local starburst galaxies? How well do the theoretical simulations model the observations? What can we learn about starbursts in the distant Universe through studies of their local counterparts? How important is the role of extreme starbursts in the hierarchical assembly of galaxies? How are extreme starbursts related to the triggering of AGN in the nuclei of galaxies? Overall, 41 talks and 4 posters with their corresponding 10 minutes short talks were presented during the workshop. In addition, the workshop was designed with emphasis on discussions, and therefore, there were 6 discussion sessions of up to one hour during the workshop. Here is presented a summary of the purposes of the workshop as well as a compilation of the abstracts corresponding to each of the presentations. The summary and conclusions of the workshop along with a description of the future prospects by Sylvain Veilleux can be found in the last section of this document. A photo of the assistants is included.

Strong inflows of gas from the outer disk to the inner kiloparsecs are induced during the interaction of disk galaxies. This inflow of relatively low-metallicity gas dilutes the metallicity of the circumnuclear gas. We have investigated several aspects of the process as the timing and duration of the dilution and its correlation with the induced star formation. We analysed major (1:1) gas-rich interactions and mergers, spanning a range of initial orbital characteristics. Star formation and metal enrichment from SNe are included in our model. Our results show that the strongest trend is between the star formation rate and the dilution of the metals in the nuclear region; i.e., the more intense the central burst of star formation, the more the gas is diluted. This trend comes from strong inflows of relatively metal-poor gas from the outer regions of both disks, which fuels the intense star formation and lowers the overall metallicity for a time. The strong inflows happen on timescales of about 10^8 years or less, and the most intense star formation and lowest gas phase metallicities are seen generally after the first pericentre passage. As the star formation proceeds and the merger advances, the dilution reduces and enrichment becomes dominant - ultimately increasing the metallicity of the circumnuclear gas to a level higher than the initial metallicities of the merging galaxies. The "fly-bys" - pairs that interact but do not merge - also cause some dilution. We even see some dilution early in the merger or in the "fly-bys" and thus do not observe a strong trend between the nuclear metallicities and separation in our simulations until the merger is well advanced. We also analyse the O and Fe enrichment of the ISM, and show that the evolution of the alpha/Fe ratios, as well as the dilution of the central gas metallicity, can be used as a clock for "dating" the interaction.

We present the GalMer database, a library of galaxy merger simulations, made available to users through tools compatible with the Virtual Observatory (VO) standards adapted specially for this theoretical database. To investigate the physics of galaxy formation through hierarchical merging, it is necessary to simulate galaxy interactions varying a large number of parameters: morphological types, mass ratios, orbital configurations, etc. On one side, these simulations have to be run in a cosmological context, able to provide a large number of galaxy pairs, with boundary conditions given by the large-scale simulations, on the other side the resolution has to be high enough at galaxy scales, to provide realistic physics. The GalMer database is a library of thousands simulations of galaxy mergers at moderate spatial resolution and it is a compromise between the diversity of initial conditions and the details of underlying physics. We provide all coordinates and data of simulated particles in FITS binary tables. The main advantages of the database are VO access interfaces and value-added services which allow users to compare the results of the simulations directly to observations: stellar population modelling, dust extinction, spectra, images, visualisation using dedicated VO tools. The GalMer value-added services can be used as virtual telescope producing broadband images, 1D spectra, 3D spectral datacubes, thus making our database oriented towards the usage by observers. We present several examples of the GalMer database scientific usage obtained from the analysis of simulations and modelling their stellar population properties, including: (1) studies of the star formation efficiency in interactions; (2) creation of old counter-rotating components; (3) reshaping metallicity profiles in elliptical galaxies; (4) orbital to internal angular momentum transfer; (5) reproducing observed colour bimodality of galaxies.

We have investigated the impact of dissipationless minor galaxy mergers on the angular momentum of the remnant. Our simulations cover a range of initial orbital characteristics and the system consists of a massive galaxy with a bulge and disk merging with a much less massive (one-tenth or one-twentieth) gasless companion which has a variety of morphologies (disk- or elliptical-like) and central baryonic mass concentrations. During the process of merging, the orbital angular momentum is redistributed into the internal angular momentum of the final system; the internal angular momentum of the primary galaxy can increase or decrease depending on the relative orientation of the orbital spin vectors (direct or retrograde), while the initially non-rotating dark matter halo always gains angular momentum. The specific angular momentum of the stellar component always decreases independent of the orbital parameters or morphology of the satellite, the decrease in the rotation velocity of the primary galaxy is accompanied by a change in the anisotropy of the orbits, and the ratio of rotation speed to velocity dispersion of the merger remnant is lower than the initial value, not only due to an increase in the dispersion but also to the slowing -down of the disk rotation. We briefly discuss several astrophysical implications of these results, suggesting that minor mergers do not cause a "random walk" process of the angular momentum of the stellar disk component of galaxies, but rather a steady decrease. Minor mergers may play a role in producing the large scatter observed in the Tully-Fisher relation for S0 galaxies, as well as in the increase of the velocity dispersion and the decrease in $v/\sigma$ at large radii as observed in S0 galaxies.

Many evidences exist for a connection between galaxy interactions and induced star formation. However, a large range of responses of galaxies to tidal interactions is found, both in observations and in numerical simulations. We will discuss some recent results obtained analysing a large sample (~ 1000) of simulations of interacting pairs and their agreement with the most recent observational works.

Through the analysis of a set of numerical simulations of major mergers between initially non-rotating, pressure supported progenitor galaxies with a range of central mass concentrations, we have shown that: (1) it is possible to generate elliptical-like galaxies, with v/sigma > 1 outside one effective radius, as a result of the conversion of orbital- into internal-angular momentum; (2) the outer regions acquire part of the angular momentum first; (3) both the baryonic and the dark matter components of the remnant galaxy acquire part of the angular momentum, the relative fractions depend on the initial concentration of the merging galaxies. For this conversion to occur the initial baryonic component must be sufficiently dense and/or the encounter should take place on a orbit with high angular momentum. Systems with these hybrid properties have been recently observed through a combination of stellar absorption lines and planetary nebulae for kinematic studies of early-type galaxies. Our results are in qualitative agreement with such observations and demonstrate that even mergers composed of non-rotating, pressure-supported progenitor galaxies can produce early-type galaxies with significant rotation at large radii.

Lenticular galaxies are believed to form by a combination of environmental effects and secular evolution. We study the nearby disc-dominated S0 galaxy NGC 6340 photometrically and spectroscopically to understand the mechanisms of S0 formation and evolution in groups. We use SDSS images to build colour maps and light profile of NGC 6340 which we decompose using a three-component model including Sersic and two exponential profiles. We also use Spitzer images to study the morphology of regions containing warm ISM and dust. Then, we re-process and re-analyse deep long-slit spectroscopic data for NGC 6340 and recover its stellar and gas kinematics, distribution of age and metallicity with the NBursts full spectral fitting. We obtain the profiles of internal kinematics, age, and metallicity out to >2 half-light radii. The three structural components of NGC 6340 are found to have distinct kinematical and stellar population properties. We see a kinematical misalignment between inner and outer regions of the galaxy. We confirm the old metal-rich centre and a wrapped inner gaseous polar disc (r~1 kpc) having weak ongoing star formation, counter-rotating in projection with respect to the stars. The central compact pseudo-bulge of NGC 6340 looks very similar to compact elliptical galaxies. In accordance with the results of numerical simulations, we conclude that properties of NGC 6340 can be explained as the result of a major merger of early-type and spiral galaxies which occurred about 12 Gyr ago. The intermediate exponential structure might be a triaxial pseudo-bulge formed by a past bar structure. The inner compact bulge could be the result of a nuclear starburst triggered by the merger. The inner polar disc appeared recently, 1/3-1/2 Gyr ago as a result of another minor merger or cold gas accretion.

Using a large suite of galaxies with a variety of concentrations and metallicity gradients, we study the evolution of non-dissipative ("dry") equal mass mergers. Our purpose in generating these simulations is to explore how the metallicity gradient in dry mergers depends on the structure and metallicity gradients of the galaxies involved in the merger. Specifically, we would like to answer: Could dry mergers lead to metallicity gradients as observed in elliptical galaxies in the local Universe? Do dry mergers always lead to a flattening of the initial (i.e., pre-merger) gradient? From this modeling, we conclude that: The ratio of the remnant and the initial galaxy slopes span a wide range of values, up to values greater than 1 (with values greater than one resulting only when companions have gradients twice the progenitor). For a merger between two ellipticals having identical initial metallicity slopes (i.e., equal companion and galaxy slopes), the metallicity profile of the remnant flattens, with a final gradient about 0.6 times the initial one. Ellipticals can maintain their original pre-merger metallicity gradient if the companion slope is sufficiently steep. The amount of flattening neither depends on the characteristics of the orbit of the progenitors or on their initial concentration. Given the diversity in outcomes of the mergers, we conclude that dry mergers do not violate any observational constraints on the systematic characteristics of metallicity gradients in local ellipticals. In fact, dry mergers may be important within the context of the results of our simulations and may explain the large scatter in gradients for massive ellipticals and the relative paucity of massive ellipticals with no or shallow metallicity gradients.

We analyze the physical conditions in the interstellar gas of 11 actively star-forming galaxies at z~2, based on integral-field spectroscopy from the ESO-VLT and HST/NICMOS imaging. We concentrate on the high H-alpha surface brightnesses, large line widths, line ratios and the clumpy nature of these galaxies. We show that photoionization calculations and emission line diagnostics imply gas pressures and densities that are similar to the most intense nearby star-forming regions at z=0 but over much larger scales (10-20 kpc). A relationship between surface brightness and velocity dispersion can be explained through simple energy injection arguments and a scaling set by nearby galaxies with no free parameters. The high velocity dispersions are a natural consequence of intense star formation thus regions of high velocity dispersion are not evidence for mass concentrations such as bulges or rings. External mechanisms like cosmological gas accretion generally do not have enough energy to sustain the high velocity dispersions. In some cases, the high pressures and low gas metallicites may make it difficult to robustly distinguish between AGN ionization cones and star formation, as we show for BzK-15504 at z=2.38. We construct a picture where the early stages of galaxy evolution are driven by self-gravity which powers strong turbulence until the velocity dispersion is high. Then massive, dense, gas-rich clumps collapse, triggering star formation with high efficiencies and intensities as observed. At this stage, the intense star formation is likely self-regulated by the mechanical energy output of massive stars.

The 21 cm emission of neutral hydrogen is the most promising probe of the epoch of reionization(EoR). In the next few years, the SKA pathfinders will provide statistical measurements of this signal. Numerical simulations predicting these observations are necessary to optimize the design of the instruments. The main difficulty is the computation of the spin temperature of neutral hydrogen which depends on the gas kinetic temperature and on the level of the local Lyman-alpha flux. A T_s >> T_cmb assumption is usual. However, this assumption does not apply early in the reionization history, or even later in the history as long as the sources of X-rays are too weak to heat the intergalactic medium significantly. This work presents the first EoR numerical simulations including, beside dynamics and ionizing continuum radiative transfer, a self-consistent treatment of the Ly-alpha radiative transfer. This allows us to compute the spin temperature more accurately. We use two different box sizes, 20 Mpc/h and 100 Mpc/h, and a star source model. Using the redshift dependence of average quantities, maps, and power spectra, we quantify the effect of using different assumptions to compute the spin temperature and the influence of the box size. The first effect comes from allowing for a signal in absorption. The magnitude of this effect depends on the amount of heating by hydrodynamic shocks and X-rays in the intergalactic medium(IGM). The second effects comes from using the real, local, Lyman-alpha flux. This effect is important for an average ionization fraction of less than 10%: it changes the overall amplitude of the 21 cm signal, and adds its own fluctuations to the power spectrum.

We investigate the intensity enhancement and the duration of starburst episodes, triggered by major galaxy interactions and mergers. To this aim, we analyze two large statistical datasets of numerical simulations. These have been obtained using two independent and different numerical techniques to model baryonic and dark matter evolution, that are extensively compared for the first time. One is a Tree-SPH code, the other one is a grid-based N-body sticky-particles code. We show that, at low redshift, galaxy interactions and mergers in general trigger only moderate star formation enhancements. Strong starbursts where the star formation rate is increased by a factor larger than 5 are rare and found only in about 15% of major galaxy interactions and mergers. Merger-driven starbursts are also rather short-lived, with a typical duration of the activity of a few 10^8 yr. These conclusions are found to be robust, independent from the numerical techniques and star formation models. At higher redshifts where galaxies contain more gas, gas inflow-induced starbursts are neither stronger neither longer than their local counterparts. In turn, the formation of massive gas clumps, results of local Jeans instability that can occur spontaneously in gas-rich disks or be indirectly favored by galaxy interactions, could play a more important role in determining the duration and intensity of star formation episodes.

We investigate the properties of a clump-cluster galaxy at redshift 1.57. The morphology of this galaxy is dominated by eight star-forming clumps in optical observations, and has photometric properties typical of most clump-cluster and chain galaxies. Its complex asymmetrical morphology has led to the suggestion that this system is a group merger of several initially separate proto-galaxies. We performed H_alpha integral field spectroscopy of this system using SINFONI on VLT UT4. These observations reveal a large-scale velocity gradient throughout the system, but with large local kinematic disturbances. Using a numerical model of gas-rich disk fragmentation, we find that clump interactions and migration can account for the observed disturbed rotation. On the other hand, the global rotation would not be expected for a multiply merging system. We further find that this system follows the stellar mass vs. metallicity, star formation rate and size relations expected for a disk at this redshift, and exhibits a disk-like radial metallicity gradient, so that the scenario of internal disk fragmentation is the most likely one. A red and metallic central concentration appears to be a bulge in this proto-spiral clumpy galaxy. A chain galaxy at redshift 2.07 in the same field also shows disk-like rotation. Such systems are likely progenitors of the present-day bright spiral galaxies, forming their exponential disks through clump migration and disruption and fueling their bulges. Our present results show that disturbed morphologies and kinematics are not necessarily signs of galaxy mergers and interactions, and can instead result from the internal evolution of primordial disks.

We investigate, by means of numerical simulations, the kinematics of elliptical-spiral merger remnants. Counterrotation can appear both in coplanar and in non-coplanar retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the external regions of the two interacting galaxies acquire part of the orbital angular momentum, due to the action of tidal forces.

We investigate, by means of numerical simulations, the possibility of forming counter-rotating old stellar components by major mergers between an elliptical and a spiral galaxy. We show that counter-rotation can appear both in dissipative and dissipationless retrograde mergers, and it is mostly associated to the presence of a disk component, which preserves part of its initial spin. In turn, the external regions of the two interacting galaxies acquire part of the orbital angular momentum, due to the action of tidal forces exerted on each galaxy by the companion.

We have undertaken a large set of simulations of galaxy interactions and mergers (GalMer Project) in order to study the physical processes related to galaxy encounters. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct and retrograde, according to the initial relative energy and impact parameter. The self gravity of stars, gas and dark matter is taken into account through a tree-code algorithm, the gas hydrodynamics through SPH. Star formation is included adopting a density-dependent Schmidt law. This wide library of galaxy interactions and mergers, containing, at present, about 900 simulations of major encounters, represents an unique tool to investigate statistically the chemodynamical evolution of interacting systems. In the following, we present and discuss some results obtained exploring the dataset, together with some future perspectives.

We present the results of detailed N-body simulations of clusters moving in a realistic Milky Way (MW) potential. The strong interaction with the bulge and the disk of the Galaxy leads to the formation of tidal tails, emanating from opposite sides of the cluster. Their orientation and morphology may be interpreted easily in terms of a comoving frame of coordinates.

We investigate the enhancement of star formation efficiency in galaxy interactions and mergers, by numerical simulations of several hundred galaxy collisions. All morphological types along the Hubble sequence are considered in the initial conditions of the two colliding galaxies, with varying bulge-to-disk ratios and gas mass fractions. Different types of orbits are simulated, direct and retrograde, according to the initial relative energy and impact parameter, and the resulting star formation history is compared to that occuring in the two galaxies when they are isolated. Our principal results are: (1) retrograde encounters have a larger star formation efficiency (SFE) than direct encounters; (2) the amount of gas available in the galaxy is not the main parameter governing the SFE in the burst phase; (3) there is an anticorrelation between the amplitude of the star forming burst and the tidal forces exerted per unit of time, which is due to the large amount of gas dragged outside the galaxy by tidal tails in strong interactions; (4) globally, the Kennicutt-Schmidt law is retrieved statistically for isolated galaxies, interacting pairs and mergers; (5) the enhanced star formation is essentially occurring in nuclear starbursts, triggered by inward gas flows driven by non-axisymmetries in the galaxy disks. Direct encounters develop more pronounced asymmetries than retrograde ones. Based on these statistical results, we derive general laws for the enhancement of star formation in galaxy interactions and mergers, as a function of the main parameters of the encounter.

We present the results of detailed N-body simulations of clusters moving in a realistic Milky Way (MW) potential. The strong interaction with the bulge and the disk of the Galaxy leads to the formation of tidal tails, emanating from opposite sides of the cluster. Some characteristic features in the morphology and orientation of these streams are recognized and intepreted. The tails have a complex morphology, in particular when the cluster approaches its apogalacticon, showing multiple ``arms'' in remarkable similarity to the structures observed around NGC 288 and Willman 1. Actually, the tails are generally good tracers of the cluster path quite far from the cluster center (>7--8 tidal radii), while on the smaller scale they are mainly pointing in the direction of the Galaxy center. In particular, the orientation of the inner part of the tails is highly correlated to the cluster orbital phase and to the local orbital angular acceleration. This implies that, in general, the orbital path cannot be estimated directly from the orientation of the tails, unless a sufficient large field around the cluster is available.

In the last decade, observational studies have shown the existence of tidal streams in the outer part of many galactic globular clusters. The most striking examples of clusters with well defined tidal tails are represented by Palomar 5 and NGC 5466 (both observed in the framework of the Sloan Digital Sky Survey), which show structures elongated for 4 kpc and 1 kpc in length, respectively. Unfortunately, most of the observational studies about globular clusters (GCs) do not cover such a large field of the sky as the SDSS does. In this framework, by mean of a parallel, adaptive tree-code, we performed detailed N-body simulations of GCs moving in a realistic three-components (bulge, disk and halo) Milky Way potential, in order to clarify whether and to what extent tails in the clusters outer regions (few tidal radii) are tracers of the local orbits and, also, if some kind of correlation exists among the cluster orbital phase and the orientation of such streams.

We present the results of detailed N-body simulations regarding the interaction of four massive globular clusters in the central region of a triaxial galaxy. The systems undergo a full merging event, producing a sort of 'Super Star Cluster' (SSC) whose features are close to those of a superposition of the individual initial mergers. In contrast with other similar simulations, the resulting SSC structural parameters are located along the observed scaling relations of globular clusters. These findings seem to support the idea that a massive SSC may have formed in early phases of the mother galaxy evolution and contributed to the growth of a massive nucleus.

The main topic of this paper is the investigation of the modes of interaction of globular clusters (GCs) moving in the inner part of a galaxy. This is tackled by means of high-resolution N-body simulations, whose first results are presented in this article. Our simulations dealt with primordial very massive (order of 10^7 solar masses) GCs that were able to decay, because of dynamical friction, into the inner regions of triaxial galaxies on a time much shorter than their internal relaxation time. To check the disruptive role of both tidal forces and GC-GC collisions, we maximised the tidal interaction considering GCs on quasi-radial orbits. The available CPU resources allowed us to follow 8 oscillations of the GCs along their orbits and the main findings are: i) clusters with an initial high enough King concentration parameter (c>=1.2), preserve up to 50% of their initial mass; ii) the inner density distribution of the survived clusters keep a King model profile; iii) GC-GC collisions have a negligible effect with respect to that caused by the passage through the galactic center; iv) the orbital energy dissipation due to the tidal interaction is of the same order of that caused by dynamical friction; v) complex sub-structures like "ripples" and "clumps" formed, as observed around real clusters. These findings support the validity of the hypothesis of merging of GCs in the galactic central region, with modes that deserve further careful investigations.

Astrocomp is a joint project, developed by the INAF-Astrophysical Observatory of Catania, University of Roma La Sapienza and Enea. The project has the goal of providing the scientific community of a web-based user-friendly interface which allows running parallel codes on a set of high-performance computing (HPC) resources, without any need for specific knowledge about parallel programming and Operating Systems commands. Astrocomp provides, also, computing time on a set of parallel computing systems, available to the authorized user. At present, the portal makes a few codes available, among which: FLY, a cosmological code for studying three-dimensional collisionless self-gravitating systems with periodic boundary conditions; ATD, a parallel tree-code for the simulation of the dynamics of boundary-free collisional and collisionless self-gravitating systems and MARA, a code for stellar light curves analysis. Other codes are going to be added to the portal.

We present some results of numerical simulations of a globular cluster orbiting in the central region of a triaxial galaxy on a set of 'loop' orbits. Tails start forming after about a quarter of the globular cluster orbital period and develop, in most cases, along the cluster orbit, showing clumpy substructures as observed, for example, in Palomar 5. If completely detectable, clumps can contain about 7,000 solar masses each, i.e. about 10% of the cluster mass at that epoch. The morphology of tails and clumps and the kinematical properties of stars in the tails are studied and compared with available observational data. Our finding is that the stellar velocity dispersion tends to level off at large radii, in agreement to that found for M15 and Omega Centauri.

The modern study of the dynamics of stellar systems requires the use of high-performance computers. Indeed, an accurate modelization of the structure and evolution of self-gravitating systems like planetary systems, open clusters, globular clusters and galaxies imply the evaluation of body-body interaction over the whole size of the structure, a task that is computationally very expensive, in particular when it is performed over long intervals of time. In this report we give a concise overview of the main problems of stellar systems simulations and present some exciting results we obtained about the interaction of globular clusters with the parent galaxy.