Evolution Of Galactic Disks Research Articles

Comparing radial migration in dark matter and MOND regimes

Multiple studies on radial migration in disc galaxies have proven the importance of the effect of resonances with non-axisymmetric components on the evolution of galactic discs. However, the dynamical effects of classic Newtonian dynamics with dark matter (DM) differ from MOdified Newtonian Dynamics (MOND) and might trigger different radial migration. A thorough analysis of radial migration considering these two gravitational regimes might shed some light on different predictions of DM and MOND theories. We aim to quantitatively and qualitatively compare the effects of resonances and stellar radial migration (churning) in a Milky Way-like (MW-like) galaxy in the DM and MOND regimes. We performed numerical simulations of a MW-like galaxy to analyse the effect of non-axisymmetric structures (galactic bar and spiral arms) considering various parameters of the spiral structure. We conducted a two-dimensional numerical simulation consisting of the integration of $2 10^6$ stars in a static rotating galactic potential for $6 Gyr $. We analysed the changes in the stellar disc caused by stellar radial migration. To study the effects of radial migration and resonances in detail, we analysed the change in the star's position, the guiding radius, as well as the frequency phase space. We investigated DM and MOND approaches. The outcome of the simulation shows that the radial migration is much more pronounced in the MOND regime compared to the DM one. Increasing both the spiral amplitude and the spiral pattern speed results in more prominent radial migration. Compared to the DM approach, in the MOND regime, we observe up to five times as many stars with a maximum change in the guiding radius of more than $1.5 kpc $ during the time interval from $2-6 Gyr $. Increasing the amplitude of the spiral structure amplifies this effect. On the other hand, increasing the spiral pattern speed reduces the difference between DM and MOND. Analysis of the frequency phase space reveals that the most prominent resonances in all DM and MOND configurations are the co-rotation resonance with the spiral arms ($m=p=1$), outer Lindblad resonance with the galactic bar and spiral arms, and the co-rotation resonance ($m=2$, $p=1$) with the superposition of the galactic bar and spiral arms, $2 sp

Large-scale stellar age-velocity spiral pattern in NGC 4030

The processes driving the formation and evolution of late-type galaxies continue to be a debated subject in extragalactic astronomy. Investigating stellar kinematics, especially when combined with age estimates, provides crucial insights into the formation and subsequent development of galactic discs. Post-processing of exceptionally high-quality integral field spectroscopy data of NGC 4030 acquired with the Multi Unit Spectroscopic Explorer (MUSE) has revealed a striking grand design spiral pattern in the velocity dispersion map, that has not been detected in other galaxies. This pattern spatially correlates with HII regions, suggesting that stars currently being born exhibit lower velocity dispersion as compared to surrounding areas where star-formation is less active. We examined the age-velocity relation (AVR) and propose that its configuration might be shaped by a combination of heating mechanisms, seemingly consistent with findings from recent high-resolution cosmological zoom-in simulations. The complex structure of the uncovered AVR of NGC 4030 supports the hypothesis that stellar populations initially inherit the velocity dispersion sigma of the progenitor cold molecular gas, which depends on formation time and galactocentric distance, subsequently experiencing kinematic heating due to cumulative gravitational interactions during their lifetime. While advancing our understanding of the AVR, these findings also offer a new framework for investigating disc heating mechanisms and their role in the evolution of galactic discs.

The impact of spurious collisional heating on the morphological evolution of simulated galactic discs

ABSTRACT We use a suite of idealized N-body simulations to study the impact of spurious heating of star particles by dark matter (DM) particles on the kinematics and morphology of simulated galactic discs. We find that spurious collisional heating leads to a systematic increase of the azimuthal velocity dispersion (σϕ) of stellar particles and a corresponding decrease in their mean azimuthal velocities ($\overline{v}_\phi$). The rate of heating is dictated primarily by the number of DM halo particles (or equivalently, by the DM particle mass at fixed halo mass) and by radial gradients in the local DM density along the disc; it is largely insensitive to the stellar particle mass. Galaxies within haloes resolved with fewer than ≈106 DM particles are particularly susceptible to spurious morphological evolution, irrespective of the total halo mass (with even more particles required to prevent heating of the galactic centre). Collisional heating transforms galactic discs from flattened structures into rounder spheroidal systems, causing them to lose rotational support in the process. It also affects the locations of galaxies in standard scaling relations that link their various properties: at fixed stellar mass, it increases the sizes of galaxies, and reduces their mean stellar rotation velocities and specific angular momenta. Our results urge caution when extrapolating simulated galaxy scaling relations to low masses where spurious collisional effects can bias their normalization, slope, and scatter.

The evolution of turbulent galactic discs: gravitational instability, feedback, and accretion

ABSTRACT We study the driving of turbulence in star-forming disc galaxies of different masses at different epochs, using an analytic ‘bathtub’ model. The disc of gas and stars is assumed to be in marginal Toomre instability. Turbulence is assumed to be sustained via an energy balance between its dissipation and three simultaneous energy sources. These are stellar feedback, inward transport due to disc instability and clumpy accretion via streams. The transport rate is computed with two different formalisms, with similar results. To achieve the energy balance, the disc self-regulates either the mass fraction in clumps or the turbulent viscous torque parameter. In this version of the model, the efficiency by which the stream kinetic energy is converted into turbulence is a free parameter, ξa. We find that the contributions of the three energy sources are in the same ball park, within a factor of ∼2 in all discs at all times. In haloes that evolve to a mass $\le 10^{12}\, {\rm M_{\odot }}$ by z = 0 ($\le 10^{11.5}\, {\rm M_{\odot }}$ at z ∼ 2), feedback is the main driver throughout their lifetimes. Above this mass, the main driver is either transport or accretion for very low or very high values of ξa, respectively. For an assumed ξa(t) that declines in time, galaxies in haloes with present-day mass >1012 M⊙ make a transition from accretion to transport dominance at intermediate redshifts, z ∼ 3, when their mass was ${\ge }10^{11.5}\, {\rm M_{\odot }}$. The predicted relation between star formation rate and gas velocity dispersion is consistent with observations.

The AMBRE Project: Solar neighbourhood chemodynamical constraints on Galactic disc evolution

Context. The abundance of α-elements relative to iron ([α/Fe]) is an important fossil signature in Galactic archaeology for tracing the chemical evolution of disc stellar populations. High-precision chemical abundances, together with accurate stellar ages, distances, and dynamical data, are crucial to infer the Milky Way formation history. Aims. The aim of this paper is to analyse the chemodynamical properties of the Galactic disc using precise magnesium abundance estimates for solar neighbourhood stars with accurate Gaia astrometric measurements. Methods. We estimated ages and dynamical properties for 366 main sequence turn-off stars from the AMBRE Project using PARSEC isochrones together with astrometric and photometric values from Gaia DR2. We use precise global metallicities [M/H] and [Mg/Fe] abundances from a previous study in order to estimate gradients and temporal chemodynamic relations for these stars. Results. We find a radial gradient of −0.099 ± 0.031 dex kpc−1 for [M/H] and +0.023 ± 0.009 dex kpc−1 for the [Mg/Fe] abundance. The steeper [Mg/Fe] gradient than that found in the literature is a result of the improvement of the AMBRE [Mg/Fe] estimates in the metal-rich regime. In addition, we find a significant spread of stellar age at any given [Mg/Fe] value, and observe a clear correlated dispersion of the [Mg/Fe] abundance with metallicity at a given age. While for [M/H] ≤ − 0.2, a clear age–[Mg/Fe] trend is observed, more metal-rich stars display ages from 3 up to 12 Gyr, describing an almost flat trend in the [Mg/Fe]–age relation. Moreover, we report the presence of radially migrated and/or churned stars for a wide range of stellar ages, although we note the large uncertainties of the amplitude of the inferred change in orbital guiding radii. Finally, we observe the appearance of a second chemical sequence in the outer disc, 10–12 Gyr ago, populating the metal-poor, low-[Mg/Fe] tail. These stars are more metal-poor than the coexisting stellar population in the inner parts of the disc, and show lower [Mg/Fe] abundances than prior disc stars of the same metallicity, leading to a chemical discontinuity. Our data favour the rapid formation of an early disc that settled in the inner regions, followed by the accretion of external metal-poor gas –probably related to a major accretion event such as the Gaia-Enceladus/Sausage one– that may have triggered the formation of the thin disc population and steepened the abundance gradient in the early disc.

The role of gas fraction and feedback in the stability and evolution of galactic discs: implications for cosmological galaxy formation models

ABSTRACT High-redshift star-forming galaxies often have irregular morphologies with giant clumps containing up to 108−109 solar masses of gas and stars. The origin and evolution of giant clumps are debated both theoretically and observationally. In most cosmological simulations, high-redshift galaxies have regular spiral structures or short-lived clumps, in contradiction with many idealized high-redshift disc models. Here, we test whether this discrepancy can be explained by the low gas fractions of galaxies in cosmological simulations. We present a series of simulations with varying gas fractions, from 25 per cent, typical of galaxies in most cosmological simulations, to 50 per cent, typical of observed galaxies at 1.5 < z < 3. We find that gas-poor models have short-lived clumps, that are unbound and mostly destroyed by galactic shear, even with weak stellar feedback. In contrast, gas-rich models form long-lived clumps even with boosted stellar feedback. This shows that the gas mass fraction is the primary physical parameter driving violent disc instabilities and the evolution of giant clumps on ∼108 yr time-scales, with lower impact from the calibration of the stellar feedback. Many cosmological simulations of galaxy formation have relatively gas-poor galactic discs, which could explain why giant clumps are absent or short-lived in such models. Similar baryonic and dark matter mass distribution could produce clumpy galaxies with long-lived clumps at z ∼ 2 if the gas fraction was in better agreement with observations.

Exploring the Origin of Thick Disks Using the NewHorizon and Galactica Simulations

Abstract Ever since a thick disk was proposed to explain the vertical distribution of the Milky Way disk stars, its origin has been a recurrent question. We aim to answer this question by inspecting 19 disk galaxies with stellar mass greater than 1010 M ⊙ in recent cosmological high-resolution zoom-in simulations: galactica and NewHorizon. The thin and thick disks are reasonably reproduced by the simulations with scale heights and luminosity ratios as observed. We then spatially classify the thin and thick disks and find that the thick disk stars are older, metal-poorer, kinematically hotter, and higher in accreted star fraction, while both disks are dominated by the stars formed in situ. Half of the in situ stars in the thick disks are formed before the galaxies develop their disks, and the rest are formed in spatially and kinematically thinner disks and then thickened with time by heating. However, the 19 galaxies have various properties and evolutionary routes, highlighting the need for statistically large samples to draw general conclusions. We conclude from our simulations that the thin and thick disk components are not entirely distinct in terms of formation processes but rather markers of the evolution of galactic disks. Moreover, as the combined result of the thickening of the existing disk stars and the continued formation of young thin disk stars, the vertical distribution of stars does not change much after the disks settle, pointing to the modulation of both orbital diffusion and star formation by the same confounding factor: the proximity of galaxies to marginal stability.

Virialization of the Inner CGM in the FIRE Simulations and Implications for Galaxy Disks, Star Formation, and Feedback

Abstract We use the FIRE-2 cosmological simulations to study the formation of a quasi-static, virial-temperature gas phase in the circumgalactic medium (CGM) at redshifts 0 < z < 5 and how the formation of this virialized phase affects the evolution of galactic disks. We demonstrate that when the halo mass crosses ∼1012 M ⊙, the cooling time of shocked gas in the inner CGM (∼0.1R vir, where R vir is the virial radius) exceeds the local free-fall time. The inner CGM then experiences a transition from on average subvirial temperatures (T ≪ T vir), large pressure fluctuations, and supersonic inflow/outflow velocities to virial temperatures (T ∼ T vir), uniform pressures, and subsonic velocities. This transition occurs when the outer CGM (∼0.5R vir) is already subsonic and has a temperature ∼T vir, indicating that the longer cooling times at large radii allow the outer CGM to virialize at lower halo masses than the inner CGM. This outside-in CGM virialization scenario is in contrast with inside-out scenarios commonly envisioned based on more idealized simulations. We demonstrate that inner CGM virialization coincides with abrupt changes in the central galaxy and its stellar feedback: the galaxy settles into a stable rotating disk, star formation transitions from “bursty” to “steady,” and stellar-driven galaxy-scale outflows are suppressed. Our results thus suggest that CGM virialization is initially associated with the formation of rotation-dominated thin galactic disks, rather than with the quenching of star formation as often assumed.

Barred spiral galaxies in modified gravity theories

ABSTRACT When bars form within galaxy formation simulations in the standard cosmological context, dynamical friction with dark matter (DM) causes them to rotate rather slowly. However, almost all observed galactic bars are fast in terms of the ratio between corotation radius and bar length. Here, we explicitly display an 8σ tension between the observed distribution of this ratio and that in the EAGLE simulation at redshift 0. We also compare the evolution of Newtonian galactic discs embedded in DM haloes to their evolution in three extended gravity theories: Milgromian Dynamics (MOND), a model of non-local gravity, and a scalar–tensor–vector gravity theory (MOG). Although our models start with the same initial baryonic distribution and rotation curve, the long-term evolution is different. The bar instability happens more violently in MOND compared to the other models. There are some common features between the extended gravity models, in particular the negligible role played by dynamical friction − which plays a key role in the DM model. Partly for this reason, all extended gravity models predict weaker bars and faster bar pattern speeds compared to the DM case. Although the absence of strong bars in our idealized, isolated extended gravity simulations is in tension with observations, they reproduce the strong observational preference for ‘fast’ bar pattern speeds, which we could not do with DM. We confirm previous findings that apparently ‘ultrafast’ bars can be due to bar-spiral arm alignment leading to an overestimated bar length, especially in extended gravity scenarios where the bar is already fast.

The VVV open cluster project. Near-infrared sequences of NGC 6067, NGC 6259, NGC 4815, Pismis 18, Trumpler 23, and Trumpler 20

ABSTRACT Open clusters are central elements of our understanding of the Galactic disc evolution, as an accurate determination of their parameters leads to an unbiased picture of our Galaxy’s structure. Extending the analysis towards fainter magnitudes in cluster sequences has a significant impact on the derived fundamental parameters, such as extinction and total mass. We perform a homogeneous analysis of six open stellar clusters in the Galactic disc using kinematic and photometric information from the Gaia DR2 and VVV surveys: NGC 6067, NGC 6259, NGC 4815, Pismis 18, Trumpler 23, and Trumpler 20. We implement two coarse-to-fine characterization methods: first, we employ Gaussian mixture models to tag fields around each open cluster in the proper motion space, and then we apply an unsupervised machine learning method to make the membership assignment to each cluster. For the studied clusters, with ages in the ∼120–1900 Myr range, we report an increase of ∼45 per cent new member candidates on average in our sample. The data-driven selection approach of cluster members makes our catalogue a valuable resource for testing stellar evolutionary models and for assessing the cluster low-to-intermediate mass populations. This study is the first of a series intended to homogeneously reveal open cluster near-infrared sequences.

High-resolution spectroscopic study of dwarf stars in the northern sky

Context. New space missions, such as NASA TESS or ESA PLATO, will focus on bright stars, which have been largely ignored by modern large surveys, especially in the northern sky. Spectroscopic information is of paramount importance in characterising the stars and analysing planets possibly orbiting them, and in studying the Galactic disc evolution. Aims. The aim of this work was to analyse all bright (V < 8 mag) F, G, and K dwarf stars using high-resolution spectra in the selected sky fields near the northern celestial pole. Methods. The observations were carried out with the 1.65 m diameter telescope at the Molėtai Astronomical Observatory and a fibre-fed high-resolution spectrograph covering a full visible wavelength range (4000–8500 Å). The atmospheric parameters were derived using the classical equivalent width approach while the individual chemical element abundances were determined from spectral synthesis. For both tasks the one-dimensional plane-parallel LTE MARCS stellar model atmospheres were applied. The NLTE effects for the majority of elemental abundances in our sample were negligible; however, we did calculate the NLTE corrections for the potassium abundances, as they were determined from the large 7698.9 Å line. For manganese and copper we have accounted for a hyperfine splitting. Results. We determined the main atmospheric parameters, kinematic properties, orbital parameters, and stellar ages for 109 newly observed stars and chemical abundances of Na I, Mg I, Al I, Si I, Si II, S I, K I, Ca I, Ca II, Sc I, Sc II, Ti I, Ti II, V I, Cr I, Cr II, Mn I, Fe I, Fe II, Co I, Ni I, Cu I, and Zn I for 249 F, G, and K dwarf stars observed in the present study and in our previous study. The [Mg I/Fe I] ratio was adopted to define the thin-disc (α-poor) and thick-disc (α-rich) stars in our sample. We explored the behaviour of 21 chemical species in the [El/Fe I] versus [Fe I/H] and [El/Fe I] versus age planes, and compared the results with the latest Galactic chemical evolution models. We also explored [El/Fe I] gradients according to the mean Galactocentric distances and maximum height above the Galactic plane. Conclusions. We found that in the Galactic thin-disc [El/Fe I] ratios of α-elements and aluminium have a positive trend with respect to age while the trend of Mn is clearly negative. Abundances of other species do not display significant trends. While the current theoretical models are able to reproduce the generic trends of the elements, they often seem to overestimate or underestimate the observational abundances. We found that the α-element and zinc abundances have slightly positive or flat radial and vertical gradients, while gradients for the odd-Z element Na, K, V, and Mn abundances are negative.

Towards a radially resolved semi-analytic model for the evolution of disc galaxies tuned with machine learning

ABSTRACT We present a flexible, detailed model for the evolution of galactic discs in a cosmological context since z ≈ 4, including a physically motivated model for radial transport of gas and stars within galactic discs. This expansion beyond traditional semi-analytic models that do not include radial structure, or include only a prescribed radial structure, enables us to study the internal structure of disc galaxies and the processes that drive it. In order to efficiently explore the large parameter space allowed by this model, we construct a neural-network-based emulator that can quickly return a reasonable approximation for many observables we can extract from the model, e.g. the star formation rate or the half-mass stellar radius, at different redshifts. We employ the emulator to constrain the model parameters with Bayesian inference by comparing its predictions to 11 observed galaxy scaling relations at a variety of redshifts. The constrained models agree well with observations, both those used to fit the data and those not included in the fitting procedure. These models will be useful theoretical tools for understanding the increasingly detailed observational data sets from Integral Field Units (IFUs).

Effects of coplanar satellite bands on galactic disc evolution

ABSTRACT Small dwarf companions have been long thought to have minimal influence on their host galaxy’s evolution without undergoing direct impacts to the host’s disc. However, in light of recent discoveries of coplanar, corotating satellite structures around the Milky Way, Andromeda, and Centaurus A, we use an N-body/test particle simulation to show that low-mass dwarf satellites within such structures are able to exert significant influence on their host’s disc, driving spiral waves and inducing stellar scattering. This is accomplished through quasi-periodic alignments of multiple small satellites within the structure that emulate the gravitational influence of a single, larger satellite such as Sagittarius Dwarf or the Large Magellanic Cloud. We find that the coplanar, corotational nature of such structures allows for repeated alignments on short enough time-scales to overcome damping within the disc, and in a consistent enough fashion to continually drive spiral waves over the course of 2 Gyr of simulation time. The spirals driven by this phenomenon tend to be flocculent and many armed due to the irregular intervals over which alignments occur. We additionally find that while the aligned satellites are able to induce noticeable thickening of the disc, their ability to drive surface density profile evolution is secondary to other effects.

Self-gravitational Force Calculation of High-order Accuracy for Infinitesimally Thin Gaseous Disks

Abstract Self-gravitational force calculation for infinitesimally thin disks is important for studies on the evolution of galactic and protoplanetary disks. Although high-order methods have been developed for hydrodynamic and magnetohydrodynamic equations, high-order improvement is desirable for solving self-gravitational forces for thin disks. In this work, we present a new numerical algorithm that is of linear complexity and of high-order accuracy. This approach is fast since the force calculation is associated with a convolution form, and the fast calculation can be achieved using Fast Fourier Transform. The nice properties, such as the finite supports and smoothness, of basis spline functions are exploited to stably interpolate a surface density and to achieve a high-order accuracy in forces. Moreover, if the mass distribution of interest is exclusively confined within a calculation domain, the method does not require artificial boundary values to be specified before the force calculation. To validate the proposed algorithm, a series of numerical tests, ranging from first- to third-order implementations, are performed, and the results are compared with analytic expressions derived for third- and fourth-order generalized Maclaurin disks. We conclude that the improvement on the numerical accuracy is significant with the order of the method, with only little increase of the complexity of the method.

Thorium in solar twins: implications for habitability in rocky planets

We have investigated the thorium (Th) abundance in a sample of 53 thin disc solar twins covering a wide range of ages. These data provide constrains on the mantle energy budget of terrestrial planets that can be formed over the evolution of the Galaxy's thin disc. We have estimated Th abundances with an average precision of 0.025\,dex (in both [Th/H] and [Th/Fe]) through comprehensive spectral synthesis of a Th\,II line present at 4019.1290\,{\AA}, using very high resolution (R\,=\,115,000) high quality HARPS spectra obtained at the ESO La Silla Observatory. We have confirmed that there is a large energy budget from Th decay for maintaining mantle convection inside potential rocky planets around solar twins, from the Galactic thin disc formation until now, because the pristine [Th/H]$_{\rm ZAMS}$ is super-solar on average under a uniform dispersion of 0.056\,dex (varying from +0.037 up to +0.138\,dex based on linear fits against isochrone stellar age). Comparing to neodymium (Nd) and europium (Eu), two others neutron-capture elements, the stellar pristine abundance of Th follows Eu along the Galactic thin disc evolution, but it does not follow Nd, probably because neodymium has a significant contribution from the $s$-process (about 60\,per\,cent).

The Formation and Evolution of Galactic Disks with APOGEE and the Gaia Survey

Abstract We explore the structure and evolutionary history of Galactic disks with Apache Point Observatory Galactic Evolution Experiment data release 13 (DR13 hereafter) and Gaia Tycho-Gaia Astrometric Solution data. We use the [α/M] ratio to allocate stars into particular Galactic components to elucidate the chemical and dynamical properties of the thin and thick disks. The spatial motions of the sample stars are obtained in Galactic Cartesian and cylindrical coordinates. We analyze the abundance trends and metallicity and [α/M] gradients of the thick and thin disks. We confirm the existence of metal-weak thick-disk stars in Galactic disks. A kinematical method is used to select the thin- and thick-disk stars for comparison. We calculate the scale length and scale height of the kinematically and chemically selected thick and thin disks based on the axisymmetric Jeans equation. We conclude that the scale length of the thick disk is approximately equal to that of the thin disk via a kinematical approach. For the chemical selection, this disparity is about 1 kpc. Finally, we get the stellar orbital parameters and try to unveil the formation scenario of the thick disk. We conclude that the gas-rich merger and radial migration are more reasonable formation scenarios for the thick disk.

New Features of Parallel Implementation of N-Body Problems on GPU

This paper focuses on the parallel implementation of a direct $N$-body method~(particle-particle algorithm) and the application of multiple GPUs for galactic dynamics simulations. Application of a hybrid OpenMP-CUDA technology is considered for models with a number of particles $N \sim 10^5 \div 10^7$. By means of $N$-body simulations of gravitationally unstable stellar galactic we have investigated the algorithms parallelization efficiency for various Nvidia Tesla graphics processors~(K20, K40, K80). Particular attention was paid to the parallel performance of simulations and accuracy of the numerical solution by comparing single and double floating-point precisions~(SP and DP). We showed that the double-precision simulations are slower by a factor of~$1.7$ than the single-precision runs performed on Nvidia Tesla K-Series processors. We also claim that application of the single-precision operations leads to incorrect result in the evolution of the non-axisymmetric gravitating $N$-body systems. In particular, it leads to significant quantitative and even qualitative distortions in the galactic disk evolution. For instance, after $10^4$ integration time steps for the single-precision numbers the total energy, momentum, and angular momentum of a system with $N = 2^{20}$ conserve with accuracy of $10^{-3}$, $10^{-2}$ and $10^{-3}$ respectively, in comparison to the double-precision simulations these values are $10^{-5}$, $10^{-15}$ and $10^{-13}$, respectively. Our estimations evidence in favour of usage of the second-order accuracy schemes with double-precision numbers since it is more efficient than in the fourth-order schemes with single-precision numbers.

Building disc structure and galaxy properties through angular momentum: the Dark Sage semi-analytic model

We present the new semi-analytic model of galaxy evolution, DARK SAGE, a heavily modified version of the publicly available SAGE code. The model is designed for detailed evolution of galactic discs. We evolve discs in a series of annuli with fixed specific angular momentum, which allows us to make predictions for the radial and angular-momentum structure of galaxies. Most physical processes, including all channels of star formation and associated feedback, are performed in these annuli. We present the surface density profiles of our model spiral galaxies, both as a function of radius and specific angular momentum, and find the discs naturally build a pseduobulge-like component. Our main results are focussed on predictions relating to the integrated mass--specific angular momentum relation of stellar discs. The model produces a distinct sequence between these properties in remarkable agreement with recent observational literature. We investigate the impact Toomre disc instabilities have on shaping this sequence and find they are crucial for regulating both the mass and spin of discs. Without instabilities, high-mass discs would be systematically deficient in specific angular momentum by a factor of ~2.5, with increased scatter. Instabilities also appear to drive the direction in which the mass--spin sequence of spiral galaxy discs evolves. With them, we find galaxies of fixed mass have higher specific angular momentum at later epochs.

Instability of counter-rotating stellar disks

We use an N-body simulation, constructed using GADGET-2, to investigate an accretion flow onto an astrophysical disk that is in the opposite sense to the disk's rotation. In order to separate dynamics intrinsic to the counter-rotating flow from the impact of the flow onto the disk, we consider an initial condition in which the counter-rotating flow is in an annular region immediately exterior the main portion of the astrophysical disk. Such counter-rotating flows are seen in systems such as NGC 4826 (known as the "Evil Eye Galaxy"). Interaction between the rotating and counter-rotating components is due to two-stream instability in the boundary region. A multi-armed spiral density wave is excited in the astrophysical disk and a density distribution with high azimuthal mode number is excited in the counter-rotating flow. Density fluctuations in the counter-rotating flow aggregate into larger clumps and some of the material in the counter-rotating flow is scattered to large radii. Accretion flow processes such as this are increasingly seen to be of importance in the evolution of multi-component galactic disks.

Simulations of magnetized multiphase galactic disc regulated by supernovae explosions

What exactly controls star formation in the Galaxy remains controversial. In particular, the role of feedback and magnetic field are still partially understood. We investigate the role played by supernovae feedback and magnetic field onto the star formation and the structure of the Galactic disk. We perform numerical simulations of the turbulent, magnetized, self-gravitating, multi-phase, supernovae regulated ISM within a 1 kpc stratified box. We implemented various schemes for the supernovae. This goes from a random distribution at a fixed rate to distributions for which the supernovae are spatially and temporally correlated to the formation of stars. To study the influence of magnetic field on star formation, we perform both hydrodynamical and magneto-hydrodynamical simulations. We find that supernovae feedback has a drastic influence on the galactic evolution. The supernovae distribution is playing a very significant role. When the supernovae are not correlated with star formation events, they do not modify significantly the very high star formation rate obtained without feedback. When the supernovae follow the accretion, the star formation rate can be reduced by a factor up to 30. Magnetic field is also playing a significant role. It reduces the star formation rate by a factor up to 2-3 and reduces the number of collapse sites by a factor of about 2. The exact correlation between the supernovae and the dense gas appears to have significant consequences on the galactic disk evolution and the star formation. This implies that small scale studies are necessary to understand and quantify the feedback efficiency. Magnetic field does influence the star formation at galactic scales by reducing the star formation rate and the number of star formation sites.