Chemical evolution with radial mixing

Models of the chemical evolution of our Galaxy are extended to include radial migration of stars and flow of gas through the disc. The models track the production of both iron and α-elements. A model is chosen that provides an excellent fit to the metallicity distribution of stars in the Geneva–Copenhagen survey (GCS) of the solar neighbourhood and a good fit to the local Hess diagram. The model provides a good fit to the distribution of GCS stars in the age–metallicity plane, although this plane was not used in the fitting process. Although this model's star formation rate is monotonically declining, its disc naturally splits into an α-enhanced thick disc and a normal thin disc. In particular, the model's distribution of stars in the ([O/Fe], [Fe/H]) plane resembles that of Galactic stars in displaying a ridge line for each disc. The thin-disc's ridge line is entirely due to stellar migration, and there is the characteristic variation of stellar angular momentum along it that has been noted by Haywood in survey data. Radial mixing of stellar populations with high  σ _z   from inner regions of the disc to the solar neighbourhood provides a natural explanation of why measurements yield a steeper increase of  σ _z   with age than predicted by theory. The metallicity gradient in the interstellar medium is predicted to be steeper than in earlier models, but appears to be in good agreement with data for both our Galaxy and external galaxies. The models are inconsistent with a cut-off in the star formation rate at low gas surface densities. The absolute magnitude of the disc is given as a function of time in several photometric bands, and radial colour profiles are plotted for representative times.     

Disc–halo interaction – I. Three-dimensional evolution of the Galactic disc

The results of a three-dimensional model for disc–halo interaction are presented here. The model considers explicitly the input of energy and mass by isolated and correlated supernovae in the disc. Once disrupted by the explosions, the disc never returns to its initial state. Instead it approaches a state where a thin H  i disc is formed in the Galactic plane, overlaid by thick H  i and H  ii gas discs with scaleheights of 500 pc and 1–1.5 kpc, respectively. The upper parts of the thick H  ii disc (the diffuse ionized medium) act as a disc–halo interface, and its formation and stability are directly correlated to the supernova rate per unit area in the simulated disc.     

Origin and structure of the Galactic disc(s)

We examine the chemical and dynamical structure in the solar neighbourhood of a model Galaxy that is the endpoint of a simulation of the chemical evolution of the Milky Way in the presence of radial mixing of stars and gas. Although the simulation's star formation rate declines monotonically from its unique peak and no merger or tidal event ever takes place, the model replicates all known properties of a thick disc, as well as matching special features of the local stellar population such as a metal-poor extension of the thin disc that has high rotational velocity. We divide the disc by chemistry and relate this dissection to observationally more convenient kinematic selection criteria. We conclude that the observed chemistry of the Galactic disc does not provide convincing evidence for a violent origin of the thick disc, as has been widely claimed.     

Simulations of minor mergers – II. The phase-space structure of thick discs

We analyse the phase-space structure of simulated thick discs that are the result of a 5:1 mass-ratio merger between a disc galaxy and a satellite. Our main goal is to establish what would be the imprints of a merger origin for the Galactic thick disc. We find that the spatial distribution predicted for thick-disc stars is asymmetric, seemingly in agreement with recent observations of the Milky Way thick disc. Near the Sun, the accreted stars are expected to rotate more slowly, to have broad velocity distributions and to occupy preferentially the wings of the line-of-sight velocity distributions. The majority of the stars in our model thick discs have low eccentricity orbits (in clear reference to the pre-existing heated disc) which give rise to a characteristic (sinusoidal) pattern for their line-of-sight velocities as a function of galactic longitude. The z -component of the angular momentum of thick-disc stars provides a clear discriminant between stars from the pre-existing disc and those from the satellite, particularly at large radii. These results are robust against the particular choices of initial conditions made in our simulations.     

Dynamics of an interacting luminous disc, dark halo and satellite companion

This paper describes a method, based on linear perturbation theory, to determine the dynamical interaction between extended halo and spheroid components and an environmental disturbance. One finds that resonant interaction between a galaxy and passing interlopers or satellite companions can carry the disturbance inward, deep inside the halo, where it can perturb the disc.    Applied to the Milky Way for example, the LMC and SMC appear to be sufficient to cause the observed Galactic warp and possibly seed other asymmetries. This is a multi-scale interaction in which the halo wake has a feature at roughly half the satellite orbital radius owing to a 2:1 orbital resonance. The rotating disturbance then excites an m  = 1 vertical disc mode which has the classic integral-sign morphology. A polar satellite orbit produces the largest warp and therefore the inferred LMC orbit is nearly optimal for maximum warp production.   Both the magnitude and morphology of the response depend on the details of the disc and halo models. Most critically, a change in the halo profile will shift the resonant frequencies and response location and consequently alter the coupling to the bending disc. Increasing the halo support relative to the disc, a sub-maximal disc model, decreases the warp amplitude.   Finally, the results and prognosis for N -body simulations are discussed. Discreteness noise in the halo, similar to that arising from a population of 10⁶-M⊙ black holes, can produce observable warping.     

The structure of the Galactic bar

We present a deep near-infrared wide-angle photometric analysis of the structure of the inner Galactic bar and central disc. The presence of a triaxial structure at the centre of the Galaxy is confirmed, consistent with a bar inclined at  22°± 55  from the Sun—Galactic Centre line, extending to approximately 2.5 kpc from the Galactic Centre and with a rather small axis ratio. A feature at  ℓ=−98  not aligned with this triaxiality suggests the existence of a second structure in the inner Galaxy, a double triaxiality or an inner ring. We argue that this is likely to be the signature of the end of the Galactic bar, at approximately 2.5–3 kpc, which is circumscribed by an inner pseudo-ring. No thick dust lane preceding the bar is detected and a hole in the dust distribution of the disc inside the bar radius is inferred.     

The orientations of molecular clouds in the outer Galaxy: evidence for the scale of the turbulence driver?

Supernova (SN) explosions inject a considerable amount of energy into the interstellar medium (ISM) in regions with high-to-moderate star formation rates. In order to assess whether the driving of turbulence by supernovae is also important in the outer Galactic disc, where the star formation rates are lower, we study the spatial distribution of molecular cloud (MC) inclinations with respect to the Galactic plane. The latter contains important information on the nature of the mechanism of energy injection into the ISM. We analyse the spatial correlations between the position angles (PAs) of a selected sample of MCs (the largest clouds in the catalogue of the outer Galaxy published by Heyer et al). Our results show that when the PAs of the clouds are all mapped to values into the  [0°, 90°]  interval, there is a significant degree of spatial correlation between the PAs on spatial scales in the range of 100–800 pc. These scales are of the order of the sizes of individual SN shells in low-density environments such as those prevailing in the outer Galaxy and where the metallicity of the ambient gas is of the order of the solar value or smaller. These findings suggest that individual SN explosions, occurring in the outer regions of the Galaxy and in likewise spiral galaxies, albeit at lower rates, continue to play an important role in shaping the structure and dynamics of the ISM in those regions. The SN explosions we postulate here are likely associated with the existence of young stellar clusters in the far outer regions of the Galaxy and the ultraviolet emission and low levels of star formation observed with the Galaxy Evolution Explorer (GALEX) satellite in the outer regions of local galaxies.     

Storm fronts over galaxy discs: models of how waves generate extraplanar gas and its anomalous kinematics

The existence of partially ionized, diffuse gas and dust clouds at kiloparsec scale distances above the central planes of edge-on, galaxy discs was an unexpected discovery about 20 years ago. Subsequent observations showed that this extended or extraplanar diffuse interstellar gas (EDIG) has rotation velocities approximately 10–20 per cent lower than those in the central plane, and has been hard to account for. Here, we present results of hydrodynamic models, with radiative cooling and heating from star formation. We find that in models with star formation generated stochastically across the disc, an extraplanar gas layer is generated as long as the star formation is sufficiently strong. However, this gas rotates at nearly the same speed as the midplane gas. We then studied a range of models with imposed spiral or bar waves in the disc. EDIG layers were also generated in these models, but primarily over the wave regions, not over the entire disc. Because of this partial coverage, the EDIG clouds move radially, as well as vertically, with the result that observed kinematic anomalies are reproduced. The implication is that the kinematic anomalies are the result of three-dimensional motions when the cylindrical symmetry of the disc is broken. Thus, the kinematic anomalies are the result of bars or strong waves, and more face-on galaxies with such waves should have an asymmetric EDIG component. The models also indicate that the EDIG can contain a significant fraction of cool gas, and that some star formation can be triggered at considerable heights above the disc mid-plane. We expect all of these effects to be more prominent in young, forming discs, to play a role in rapidly smoothing disc asymmetries and in working to self-regulate disc structure.     

Hydrodynamical simulations of Galactic fountains – I. Evolution of single fountains

The ejection of the gas out of the disc in late-type galaxies is related to star formation and is due mainly to Type II supernovae. In this paper, we studied in detail the development of the Galactic fountains in order to understand their dynamical evolution and their influence on the redistribution of the freshly delivered metals over the disc. To this aim, we performed a number of 3D hydrodynamical radiative cooling simulations of the gas in the Milky Way where the whole Galaxy structure, the Galactic differential rotation and the supernova explosions generated by a single OB association are considered. A typical fountain powered by 100 Type II supernovae may eject material up to ∼2 kpc which than collapses back mostly in the form of dense, cold clouds and filaments. The majority of the gas lifted up by the fountains falls back on the disc remaining within a radial distance  Δ R = 0.5 kpc  from the place where the fountain originated. This localized circulation of disc gas does not influence the radial chemical gradients on large scale, as required by the chemical models of the Milky Way which reproduce the metallicity distribution without invoking large fluxes of metals. Simulations of multiple fountains fuelled by Type II supernovae of different OB associations will be presented in a companion paper.     

Stellar populations in the Canis Major overdensity

We performed a photometric multicolour survey of the core of the Canis Major overdensity at     , reaching   V ∼ 22  and covering  0.3 × 1.0  arcmin². The main aim is to unravel the complex mixture of stellar populations toward this Galactic direction, where in the recent past important signatures of an accretion event have been claimed to be detected. While our previous investigations were based on disjointed pointings aimed at revealing the large-scale structure of the third Galactic Quadrant, we now focus on a complete coverage of a smaller field centred on the Canis Major overdensity. A large wavelength baseline, in the UBVRI bands, allows us to build up a suite of colour–colour and colour–magnitude diagrams, providing a much better diagnostic tool to disentangle the stellar populations of the region. In fact, the simple use of one colour–magnitude diagram, widely employed in all the previous studies defending the existence of the Canis Major galaxy, does not allow one to separate the effects of the different parameters (reddening, age, metallicity and distance) involved in the interpretation of data, forcing to rely on heavy modelling. In agreement with our previous studies, in the same general region of the Milky Way, we recognize a young stellar population compatible with the expected structure and extension of the Local (Orion) and Outer (Norma–Cygnus) spiral arms in the Third Galactic Quadrant. Moreover, we interpret the conspicuous intermediate-age metal-poor population as belonging to the Galactic thick disc, distorted by the effect of strong disc warping at this latitude, and to the Galactic halo.     

Could the Galactic disc heating be due to globular cluster impacts?

So far, six mechanisms have been proposed to account for the Galactic disc heating. Of these, the most important appear to be a combination of scattering of stars by molecular clouds and by spiral arms. We study a further mechanism, namely the repeated disc impact of the original Galactic globular cluster population up to the present. We find that globular clusters could have contributed at most a small fraction of the current vertical energy of the disc, as they could heat the whole disc to  σ _z = 5.5 km s⁻¹  (c.f. the observed 18 and 39 km s⁻¹ for the thick and thin discs, respectively). We find that the rate of rise of disc heat (  α= 0.22  in  σ _z ∼ t ^α  with t being time) is close to that found for scattering by molecular clouds.     

Simulations of the formation of stellar discs in the Galactic Centre via cloud–cloud collisions

Young massive stars in the central parsec of our Galaxy are best explained by star formation within at least one, and possibly two, massive self-gravitating gaseous discs. With help of numerical simulations, we here consider whether the observed population of young stars could have originated from a large angle collision of two massive gaseous clouds at   R ≃ 1 pc  from Sgr A*. In all the simulations performed, the post-collision gas flow forms an inner, nearly circular gaseous disc and one or two eccentric outer filaments, consistent with the observations. Furthermore, the radial stellar mass distribution is always very steep,  Σ_*∝ R ⁻²  , again consistent with the observations. All of our simulations produce discs that are warped by between 30° and 60°, in accordance with the most recent observations. The three-dimensional velocity structure of the stellar distribution is sensitive to initial conditions (e.g. the impact parameter of the clouds) and gas cooling details. For example, the runs in which the inner disc is fed intermittently with material possessing fluctuating angular momentum result in multiple stellar discs with different orbital orientations, contradicting the observed data. In all the cases the amount of gas accreted by our inner boundary condition is large, enough to allow Sgr A* to radiate near its Eddington limit over ∼10⁵ yr. This suggests that a refined model would have physically larger clouds (or a cloud and a disc such as the circumnuclear disc) colliding at a distance of a few parsecs rather than 1 pc as in our simulations.     

Forming a large disc galaxy from a z < 1 major merger

Using high-resolution SPH simulations in a fully cosmological Λ cold dark matter context, we study the formation of a bright disc-dominated galaxy that originates from a 'wet' major merger at   z = 0.8  . The progenitors of the disc galaxy are themselves disc galaxies that formed from early major mergers between galaxies with blue colours. A substantial thin stellar disc grows rapidly following the last major merger and the present-day properties of the final remnant are typical of early-type spiral galaxies, with an i -band bulge-to-disc ratio ∼0.65, a disc scalelength of 7.2 kpc,   g − r = 0.5 mag  , an H  i linewidth ( W ₂₀/2) of 238 km s⁻¹ and total magnitude   i =−22.4  . The key ingredients for the formation of a dominant stellar disc component after a major merger are (i) substantial and rapid accretion of gas through cold flows followed at late times by cooling of gas from the hot phase, (ii) supernova feedback that is able to partially suppress star formation during mergers and (iii) relative fading of the spheroidal component. The gas fraction of the progenitors' discs does not exceed 25 per cent at   z < 3  , emphasizing that the continuous supply of gas from the local environment plays a major role in the regrowth of discs and in keeping the galaxies blue. The results of this simulation alleviate the problem posed for the existence of disc galaxies by the high likelihood of interactions and mergers for galaxy-sized haloes at relatively low z .     

Forming stars on a viscous time-scale: the key to exponential stellar profiles in disc galaxies?

We argue for implementing star formation on a viscous time-scale in hydrodynamical simulations of disc galaxy formation and evolution. Modelling two-dimensional isolated disc galaxies with the Bhatnagar–Gross–Krook (BGK) hydrocode, we verify the analytic claim of various authors that if the characteristic time-scale for star formation is equal to the viscous time-scale in discs, the resulting stellar profile is exponential on several scalelengths whatever the initial gas and dark matter profile. This casts new light on both numerical and semi-analytical disc formation simulations that either (a) commence star formation in an already exponential gaseous disc, (b) begin a disc simulation with conditions known to lead to an exponential, i.e. the collapse of a spherically symmetric nearly uniform sphere of gas in solid-body rotation under the assumption of specific angular momentum conservation, or (c) in simulations performed in a hierarchical context, tune their feedback processes to delay disc formation until the dark matter haloes are slowly evolving and without much substructure so that the gas has the chance to collapse under conditions known to give exponentials. In such models, star formation follows a Schmidt-like law, which for lack of a suitable time-scale, resorts to an efficiency parameter. With star formation prescribed on a viscous time-scale, however, we find gas and star fractions after ∼12 Gyr that are consistent with observations without having to invoke a 'fudge factor' for star formation. Our results strongly suggest that despite our gap in understanding the exact link between star formation and viscosity, the viscous time-scale is indeed the natural time-scale for star formation.     

Simulating Gaia performances on white dwarfs

One of the most promising space missions of the European Space Agency is the astrometric satellite Gaia , which will provide very precise astrometry and multicolour photometry, for all 1.3 billion objects to   V ∼ 20  , and radial velocities with accuracies of a few km s⁻¹ for most stars brighter than   V ∼ 17  . Consequently, full homogeneous six-dimensional phase-space information for a huge number of stars will become available. Our Monte Carlo simulator has been used to estimate the number of white dwarfs potentially observable by Gaia . From this we assess the white dwarf luminosity functions that Gaia will obtain and discuss in depth the scientific returns of Gaia in the specific field of white dwarf populations. Scientifically attainable goals include, among others, a reliable determination of the age of the Galactic disc, a better knowledge of the halo of the Milky Way and the reconstruction of the star formation history of the Galactic disc. Our results also demonstrate the potential impact of a mission such as Gaia within the context of current understanding of white dwarf cooling theory.     

Hydrostatic equilibrium in a magnetized, warped Galactic disc

Hydrostatic equilibrium of the multiphase interstellar medium in the solar vicinity is reconsidered, with the regular and turbulent magnetic fields treated separately. The regular magnetic field strength required to support the gas is consistent with independent estimates, provided that energy equipartition is maintained between turbulence and random magnetic fields. Our results indicate that a mid-plane value of B ₀=4 μG for the regular magnetic field near the Sun leads to more attractive models than B ₀=2 μG . The vertical profiles of both the regular and random magnetic fields contain disc and halo components, the parameters of which we have determined. The layer at 1≲| z |≲4 kpc can be overpressured and an outflow at a speed of about 50 km s⁻¹ may occur there, presumably associated with a Galactic fountain flow, if B ₀≃2 μG .
We show that hydrostatic equilibrium in a warped disc must produce asymmetric density distributions in z , in rough agreement with H  i observations in the outer Galaxy. This asymmetry may be a useful diagnostic of the details of the warping mechanism in the Milky Way and other galaxies. We find indications that gas and magnetic field pressures are different above and below the warped midplane in the outer Galaxy, and quantify the difference in terms of turbulent velocity and/or magnetic field strength.     

On the origin of warps and the role of the intergalactic medium

There is still no consensus as to what causes galactic discs to become warped. Successful models should account for the frequent occurrence of warps in quite isolated galaxies, their amplitude as well as the observed azimuthal and vertical distributions of the H  i layer. Intergalactic accretion flows and intergalactic magnetic fields may bend the outer parts of spiral galaxies. In this paper we consider the viability of these non-gravitational torques to take the gas off the plane. We show that magnetically generated warps are clearly flawed because they would wrap up into a spiral in less than two or three galactic rotations. The inclusion of any magnetic diffusivity to dilute the wrapping effect causes the amplitude of the warp to damp. We also consider the observational consequences of the accretion of an intergalactic plane-parallel flow at infinity. We have computed the amplitude and warp asymmetry in the accretion model, for a disc embedded in a flattened dark matter halo, including self-consistently the contribution of the modes with azimuthal wavenumbers   m = 0  and   m = 1  . Since the m = 0 component, giving a U-shaped profile, is not negligible compared to the m = 1 component, this model predicts quite asymmetric warps, maximum gas displacements on the two sides in the ratio 3 : 2 for the preferred Galactic parameters, and the presence of a fraction ∼3.5 per cent of U-shaped warps, at least. The azimuthal dependence of the moment transfer by the ram pressure would produce a strong asymmetry in the thickness of the H  i layer and asymmetric density distributions in z , in conflict with observational data for the warp in our Galaxy and in external galaxies. The amount of accretion that is required to explain the Galactic warp would give gas scaleheights in the far outer disc that are too small. We conclude that accretion of a flow with no net angular momentum cannot be the main and only cause of warps.     

Disentangling the Galaxy at low Galactic latitudes

We have used the field stars from the open cluster survey BOCCE (Bologna Open Clusters Chemical Evolution), to study three low-latitude fields imaged with the Canada–France–Hawaii telescope (CFHT), with the aim of better understanding the Galactic structure in those directions. Because of the deep and accurate photometry in these fields, they provide a powerful discriminant among Galactic structure models. In the present paper we discuss if a canonical star count model, expressed in terms of thin and thick disc radial scales, thick disc normalization and reddening distribution, can explain the observed colour–magnitude diagrams (CMDs). Disc and thick disc are described with double exponentials, the spheroid is represented with a De Vaucouleurs density law. In order to assess the fitting quality of a particular set of parameters, the colour distribution and luminosity function of synthetic photometry is compared to that of target stars selected from the blue sequence of the observed CMDs. Through a Kolmogorov–Smirnov test, we find that the classical decomposition halo-thin/thick disc is sufficient to reproduce the observations – no additional population is strictly necessary. In terms of solutions common to all three fields, we have found a thick disc scalelength that is equal to (or slightly longer than) the thin disc scale.     

Radial mixing and the transition between the thick and thin Galactic discs

The analysis of the kinematics of solar neighbourhood stars shows that the low- and high-metallicity tails of the thin disc are populated by objects which orbital properties suggest an origin in the outer and inner Galactic disc, respectively. Signatures of radial migration are identified in various recent samples, and are shown to be responsible for the high-metallicity dispersion in the age–metallicity distribution. Most importantly, it is shown that the population of low-metallicity wanderers of the thin disc (−0.7 < [Fe/H] < −0.3 dex) is also responsible for the apparent hiatus in metallicity with the thick disc (which terminal metallicity is about −0.2 dex). It implies that the thin disc at the solar circle has started to form stars at about this same metallicity. This is also consistent with the fact that 'transition' objects, which have α-element abundance intermediate between that of the thick and thin discs, are found in the range [−0.4, −0.2] dex. Once the metal-poor thin disc stars are recognized for what they are – wanderers from the outer thin disc – the parenthood between the two discs can be identified on stars genuinely formed at the solar circle through an evolutionary sequence in [α/Fe] and [Fe/H]. Another consequence is that stars that can be considered as truly resulting of the chemical evolution at the solar circle have a metallicity restricted to about [−0.2, +0.2] dex, confirming an old idea that most chemical evolution in the Milky Way have preceded the thin disc formation.