Formation Of Disk Galaxies Research Articles

Formation of disc galaxies around z ∼ 2

We present combined evolution of morphological and stellar properties of galaxies on the two sides of z=2 (2.0<z<4.0 and 1.5<z<2.0) in CDFS, with ground-based spectroscopic redshifts. We perform bulge-disc decomposition on their images in J and H filters, from the 3DHST Legacy Survey obtained using HST/WFC3. Combining morphological information with stellar properties, we provide a detailed account of the formation/growth of discs and spheroids around z~2. The fraction of 2-component (bulge+disc) systems increases from ~46% for z>2 to ~70% for z<2, compensating for the fall in population of pure discs and pure spheroids. All quiescent outliers of our full sample on the main-sequence, are 2-component systems, belonging to the lower redshift range (z<2). The doubling of stellar mass of 2-component systems and decrease in their SFR by the same factor, suggests that mechanisms involved in morphological transformations are also responsible for the quenching of their star formation activity. Interestingly, while there is substantial increase in the size (~2.5 times) and mass (~5 times) of pure discs, from z>2 to z<2, pure spheroids maintain roughly the same values. Additionally, while bulge hosting discs witness an expansion in scale length (~1.3 times), their bulge sizes as well as bulge to total light ratio see no evolution, suggesting that z~2 is pre-dominantly a disc formation period.

Rapid Formation of Black Holes in Galaxies: A Self-limiting Growth Mechanism

Abstract We present high-quality fluid dynamical simulations of isothermal gas flows in a rotating barred potential. We show that a large quantity of gas is driven right into the nucleus of a galaxy when the model lacks a central mass concentration, but the inflow stalls at a nuclear ring in comparison simulations that include a central massive object. The radius of the nuclear gas ring increases linearly with the mass of the central object. We argue that bars drive gas right into the nucleus in the early stages of disk galaxy formation, where a nuclear star cluster and perhaps a massive black hole could be created. The process is self-limiting, however, because inflow stalls at a nuclear ring once the mass of gas and stars in the nucleus exceeds ∼1% of the disk mass, which shuts off rapid growth of the black hole. We briefly discuss the relevance of these results to the seeding of massive black holes in galaxies, the merger model for quasar evolution, and the existence of massive black holes in disk galaxies that lack a significant classical bulge.

Perspectives on Galactic Structure

The Milky Way is currently the subject of great observational effort. This includes both ESA's unique Gaia mission, as well as a multitude of ground-based surveys. Several of these are already returning data of unprecedented depth and quality for large numbers of Milky Way stars. These new data are likely to lead to a quantum step in our understanding of Milky Way structure and evolution. Because the new data will allow us to study our Galaxy at much greater resolution than possible in other galaxies, we also expect to greatly improve our understanding of disk galaxy formation in general.

Star Formation of Merging Disk Galaxies with AGN Feedback Effects

Abstract Using a numerical hydrodynamics code, we perform various idealized galaxy merger simulations to study the star formation (SF) of two merging disk galaxies. Our simulations include gas accretion onto supermassive black holes and active galactic nucleus (AGN) feedback. By comparing AGN simulations with those without AGNs, we attempt to understand when the AGN feedback effect is significant. Using ∼70 simulations, we investigate SF with the AGN effect in mergers with a variety of mass ratios, inclinations, orbits, galaxy structures, and morphologies. Using these merger simulations with AGN feedback, we measure merger-driven SF using the burst efficiency parameter introduced by Cox et al. We confirm previous studies which demonstrated that, in galaxy mergers, AGN suppresses SF more efficiently than in isolated galaxies. However, we also find that the effect of AGNs on SF is larger in major than in minor mergers. In minor merger simulations with different primary bulge-to-total ratios, the effect of bulge fraction on the merger-driven SF decreases due to AGN feedback. We create models of Sa-, Sb-, and Sc-type galaxies and compare their SF properties while undergoing mergers. With the current AGN prescriptions, the difference in merger-driven SF is not as pronounced as in the recent observational study of Kaviraj. We discuss the implications of this discrepancy.

A Universal Angular Momentum Profile for Dark Matter Halos

Abstract The angular momentum distribution in dark matter halos and galaxies is a key ingredient in understanding their formation. Specifically, the internal distribution of angular momenta is closely related to the formation of disk galaxies. In this article, we use halos identified from a high-resolution simulation, the Bolshoi simulation, to study the spatial distribution of specific angular momenta, . We show that by stacking halos with similar masses to increase the signal-to-noise ratio, the profile can be fitted as a simple function, , with three free parameters, , and . Specifically, j s correlates with the halo mass M vir as , r s has a weak dependence on the halo mass as , and is independent of M vir. This profile agrees with that from a rigid shell model, though its origin is unclear. Our universal specific angular momentum profile is useful in modeling the angular momenta of halos. Furthermore, by using an empirical stellar mass–halo mass relation, we can infer the average angular momentum distribution of a dark matter halo. The specific angular momentum–stellar mass relation within a halo computed from our profile is shown to share a similar shape as that from the observed disk galaxies.

Barred galaxies in the EAGLE cosmological hydrodynamical simulation

We examine the properties of barred disc galaxies in a LCDM cosmological hydrodynamical simulation from the EAGLE project. Our study follows the formation of 269 discs identified at z = 0 in the stellar mass range 10.6 < log Mstr /M < 11. These discs show a wide range of bar strengths, from unbarred discs to weak bars to strongly barred systems (= 20%). Bars in these systems develop after redshift = 1.3, on timescales that depend sen- sitively on the strength of the pattern. Strong bars develop relatively quickly (in a few Gyr, = 10 disc rotation periods) in systems that are disc dominated, gas poor, and have declining rotation curves. Weak bars develop more slowly in systems where the disc is less gravitation- ally important, and are still growing at z = 0. Unbarred galaxies are comparatively gas-rich discs whose rotation speeds do not exceed the maximum circular velocity of the halos they inhabit. Bar lengths compare favourably with observations, ranging from 0.2 to 0.8 times the radius containing 90% of the stars. Bars slow down remarkably quickly as they grow, causing the inner regions of the surrounding dark halo to expand. At z = 0 strong bars have corota- tion radii roughly ten times the bar length. Such slow bars are inconsistent with the few cases where pattern speeds have been measured or inferred observationally, a discrepancy that, if confirmed, might prove a challenge for disc galaxy formation in LCDM.

Galaxy simulations in the Gaia era

AbstractWe live in an age where an enormous amount of astrometric, photometric, asteroseismic, and spectroscopic data of Milky Way stars are being acquired, many orders of magnitude larger than about a decade ago. Thanks to the Gaia astrometric mission and followup ground-based spectroscopic surveys in the next 5-10 years about 10-20 Million stars will have accurate 6D kinematics and chemical composition measurements. KEPLER-2, PLATO, and TESS will provide asteroseismic ages for a good fraction of those. In this article we outline some outstanding problems concerning the formation and evolution of the Milky Way and argue that, due to the complexity of physical processes involved in the formation of disk galaxies, numerical simulations in the cosmological context are needed for the interpretation of Milky Way observations. We also discuss in some detail the formation of the Milky Way thick disk, chemodynamical models, and the effects of radial migration.

Forming disk galaxies in major mergers

Context: In a series of papers, we study the major merger of two disk galaxies in order to establish whether or not such a merger can produce a disc galaxy. Aims: Our aim here is to describe in detail the technical aspects of our numerical experiments. Methods: We discuss the initial conditions of our major merger, which consist of two protogalaxies on a collision orbit. We show that such merger simulations can produce a non-realistic central mass concentration, and we propose simple, parametric, AGN-like feedback as a solution to this problem. Our AGN-like feedback algorithm is very simple: at each time-step we take all particles whose local volume density is above a given threshold value and increase their temperature to a preset value. We also compare the GADGET3 and GIZMO codes, by applying both of them to the same initial conditions. Results: We show that the evolution of isolated protogalaxies resembles the evolution of disk galaxies, thus arguing that our protogalaxies are well suited for our merger simulations. We demonstrate that the problem with the unphysical central mass concentration in our merger simulations is further aggravated when we increase the resolution. We show that our AGN-like feedback removes this non-physical central mass concentration, and thus allows the formation of realistic bars. Note that our AGN-like feedback mainly affects the central region of a model, without significantly modifying the rest of the galaxy. We demonstrate that, in the context of our kind of simulation, GADGET3 gives results which are very similar to those obtained with the PSPH (density independent SPH) flavor of GIZMO. Moreover, in the examples we tried, the differences between the results of the two flavors of GIZMO, namely PSPH, and MFM (mesh-less algorithm) are similar to and, in some comparisons, larger than the differences between the results of GADGET3 and PSPH.

The Diversity of Assembly Histories Leading to Disc Galaxy Formation in a ΛCDM Model

AbstractDisc galaxies forming in a LambdaCDM cosmology often experience violent mergers. The fact that disc galaxies are ubiquitous suggests that quiescent histories are not necessary. Modern cosmological simulations can now obtain realistic populations of disc galaxies, but it is still unclear how discs manage to survive massive mergers. Here we use a suite of hydrodynamical cosmological simulations to elucidate the fate of discs encountering massive mergers. We follow the changes in the post-merger disc-to-total ratios (D/T) of simulated galaxies and examine the relations between their present-day morphology, assembly history and gas fractions. We find that approximately half of present-day disc galaxies underwent at least one merger with a satellite more massive the host's stellar component and a third had mergers with satellites three times as massive. These mergers lead to a sharp, but often temporary, decrease in the D/T of the hosts, implying that discs are usually disrupted but then quickly re-grow. To do so, high cold gas fractions are required post-merger, as well as a relatively quiescent recent history (over a few Gyrs before z = 0). Our results show that discs can form via diverse merger pathways and that quiescent histories are not the dominant mode of disc formation.

The role of gas infall in the evolution of disc galaxies

Spiral galaxies are thought to acquire their gas through a protracted infall phase resulting in the inside-out growth of their associated discs. For field spirals, this infall occurs in the lower density environments of the cosmic web. The overall infall rate, as well as the galactocentric radius at which this infall is incorporated into the star-forming disc, plays a pivotal role in shaping the characteristics observed today. Indeed, characterising the functional form of this spatio-temporal infall in-situ is exceedingly difficult, and one is forced to constrain these forms using the present day state of galaxies with model or simulation predictions. We present the infall rates used as input to a grid of chemical evolution models spanning the mass spectrum of discs observed today. We provide a systematic comparison with alternate analytical infall schemes in the literature, including a first comparison with cosmological simulations. Identifying the degeneracies associated with the adopted infall rate prescriptions in galaxy models is an important step in the development of a consistent picture of disc galaxy formation and evolution.

The spatially-resolved star formation history of the M31 outer disc

Abstract We present deep Hubble Space Telescope Advanced Camera for Surveys observations of the stellar populations in two fields lying at 20 and 23 kpc from the centre of M31 along the south-west semimajor axis. These data enable the construction of colour–magnitude diagrams reaching the oldest main-sequence turn-offs (∼13 Gyr) which, when combined with another field at 25 kpc from our previous work, we use to derive the first precision constraints on the spatially-resolved star formation history of the M31 disc. The star formation rates exhibit temporal as well as field-to-field variations, but are generally always within a factor of 2 of their time average. There is no evidence of inside-out growth over the radial range probed. We find a median age of ∼7.5 Gyr, indicating that roughly half of the stellar mass in the M31 outer disc was formed before z ∼ 1. We also find that the age–metallicity relations (AMRs) are smoothly increasing from [Fe/H] ≃ −0.4 to solar metallicity between 10 and 3 Gyr ago, contrary to the flat AMR of the Milky Way disc at a similar number of scalelengths. Our findings provide insight on the roles of stellar feedback and radial migration in the formation and evolution of large disc galaxies.

CO excitation of normal star-forming galaxies out toz= 1.5 as regulated by the properties of their interstellar medium

We investigate the CO excitation of normal star forming disk galaxies at z=1.5 using IRAM PdBI observations of the CO[2-1], CO[3-2] and CO[5-4] transitions for 4 galaxies, including VLA observations of CO[1-0] for 3 of them, with the aim of constraining the average state of H2 gas. Exploiting prior knowledge of the velocity range, spatial extent and size of the CO emission we measure reliable line fluxes with S/N>4-7 for individual transitions. While the average CO Spectral Line Energy Distribution (SLED) has a sub-thermal excitation similar to the Milky Way (MW) up to CO[3-2], we show that the average CO[5-4] emission is four times stronger than assuming MW excitation. This demonstrates the presence of an additional component of more excited, denser and possibly warmer molecular gas. The ratio of CO[5-4] to lower-J CO emission is lower than in local (U)LIRGs and high-redshift SMGs, correlating closely with the average intensity of the radiation field <U> and with the star formation surface density, but not with the SF efficiency (SFE). The CO[5-4] luminosity correlates linearly with LIR over 4 orders of magnitudes, with z=1.5 BzK galaxies following the same trend as local spirals and (U)LIRGs and high redshift star bursting SMGs. The CO[5-4] luminosity is thus empirically related to the dense gas, and might be a more convenient way to probe it than standard high--density tracers that are much fainter than CO. We see excitation variations among our sample galaxies, that can be linked to their evolutionary state and clumpiness in optical rest frame images. In one galaxy we see spatially resolved excitation variations, where the more highly excited part of the galaxy corresponds to the location of massive SF clumps. This supports to models that suggest that giant clumps are the main source of the high excitation CO emission in high redshift disk-like galaxies.

Stellar populations of the bulges of four spiral galaxies

AbstractKey information to understand the formation and evolution of disk galaxies are imprinted in the stellar populations of their bulges. This paper has the purpose to make available new measurements of the stellar population properties of the bulges of four spiral galaxies. Both the central values and radial profiles of the line strength of some of the most common Lick indices are measured along the major‐ and minor‐ axis of the bulge‐dominated region of the sample galaxies. The corresponding age, metallicity, and α /Fe ratio are derived by using the simple stellar population synthesis model predictions. The central values and the gradients of the stellar population properties of ESO‐LV 1890070, ESO‐LV 4460170, and ESO‐LV 5140100 are consistent with previous findings for bulges of spiral galaxies. On the contrary, the bulge of ESO‐LV 4500200 shows peculiar chemical properties possibly due to the presence of a central kinematically‐decoupled component. The negative metallicity gradient found in our bulges sample indicates a relevant role for the dissipative collapse in bulge formation. However, the shallow gradients found for the age and α /Fe ratio suggests that merging can not be completely ruled out for the sample bulges. This is confirmed by the properties of ESO‐LV 4500200 which can hardly be explained without invoking the capture of external material. (© 2015 WILEY‐VCH Verlag GmbH &amp; Co. KGaA, Weinheim)

Giant disc galaxies: where environment trumps mass in galaxy evolution

We identify some of the most HI massive and fastest rotating disk galaxies in the local universe with the aim of probing the processes that drive the formation of these extreme disk galaxies. By combining data from the Cosmic Flows project, which has consistently reanalyzed archival galaxy HI profiles, and 3.6$\mu$m photometry obtained with the Spitzer Space Telescope, with which we can measure stellar mass, we use the baryonic Tully-Fisher (BTF) relationship to explore whether these massive galaxies are distinct. We discuss several results, but the most striking is the systematic offset of the HI-massive sample above the BTF. These galaxies have both more gas and more stars in their disks than the typical disk galaxy of similar rotational velocity. The "condensed" baryon fraction, $f_C$, the fraction of the baryons in a dark matter halo that settle either as cold gas or stars into the disk, is twice as high in the HI-massive sample than typical, and almost reaches the universal baryon fraction in some cases, suggesting that the most extreme of these galaxies have little in the way of a hot baryonic component or cold baryons distributed well outside the disk. In contrast, the star formation efficiency, measured as the ratio of the mass in stars to that in both stars and gas, shows no difference between the HI-massive sample and the typical disk galaxies. We conclude that the star formation efficiency is driven by an internal, self-regulating process, while $f_C$ is affected by external factors. We also found that the most massive HI detected galaxies are located preferentially in filaments. We present the first evidence of an environmental effect on galaxy evolution using a dynamical definition of a filament.

Structure formation in gas-rich galactic discs with finite thickness: from discs to rings

Gravitational instabilities play an important role in structure formation of gas-rich high-redshift disc galaxies. In this paper, we revisit the axisymmetric perturbation theory and the resulting growth of structure by taking the realistic thickness of the disc into account. In the unstable regime, which corresponds for thick discs to a Toomre parameter below the critical value Q0, crit = 0.696, we find a fastest growing perturbation wavelength that is always a factor 1.93 times larger than in the classical razor-thin disc approximation. This result is independent of the adopted disc scaleheight and by this independent of temperature and surface density. In order to test the analytical theory, we compare it with a high-resolution hydrodynamical simulation of an isothermal gravitationally unstable gas disc with the typical vertical sech2 density profile and study its break up into rings that subsequently fragment into dense clumps. In the first phase, rings form, that organize themselves discretely, with distances corresponding to the local fastest growing perturbation wavelength. We find that the disc scaleheight has to be resolved initially with five or more grid cells in order to guarantee proper growth of the ring structures, which follow the analytical prediction. These rings later on contract to a thin and dense line, while at the same time accreting more gas from the inter-ring region. It is these dense, circular filaments, that subsequently fragment into a large number of clumps. Contrary to what is typically assumed, the clump sizes are therefore not directly determined by the fastest growing wavelength.

How to bend galaxy disc profiles: the role of halo spin

Abstract The radial density profiles of stellar galaxy discs can be well approximated as an exponential. Compared to this canonical form, however, the profiles in the majority of disc galaxies show downward or upward breaks at large radii. Currently, there is no coherent explanation in a galaxy formation context of the radial profile per se, along with the two types of profile breaks. Using a set of controlled hydrodynamic simulations of disc galaxy formation, we find a correlation between the host halo's initial angular momentum and the resulting radial profile of the stellar disc: galaxies that live in haloes with a low spin parameter λ ≲ 0.03 show an up-bending break in their disc density profiles, while galaxies in haloes of higher angular momentum show a down-bending break. We find that the case of pure exponential profiles (λ ≈ 0.035) coincides with the peak of the spin parameter distribution from cosmological simulations. Our simulations not only imply an explanation of the observed behaviours, but also suggest that the physical origin of this effect is related to the amount of radial redistribution of stellar mass, which is anticorrelated with λ.

Formation of disc galaxies in preheated media: a preventative feedback model

We introduce a semi-analytic galaxy formation model implementing a self-consistent treatment for the hot halo gas configuration and the assembly of central disks. Using the model, we explore a preventative feedback model, in which the circum-halo medium is assumed to be preheated up to a certain entropy level by early starbursts or other processes, and compare it with an ejective feedback model, in which baryons are first accreted into dark matter halos and subsequently ejected out by feedback. The model demonstrates that when the medium is preheated to an entropy comparable to the halo virial entropy the baryon accretion can be largely reduced and delayed. In addition, the preheated medium can establish an extended low density gaseous halo when it accretes into the dark matter halos, and result in a specific angular momentum of the cooling gas large enough to form central disks as extended as those observed. Combined with simulated halo assembly histories, the preventative feedback model can reproduce remarkably well a number of observational scaling relations. These include the cold baryon (stellar plus cold gas) mass fraction-halo mass relations, star formation histories, disk size-stellar mass relation and its evolution, and the number density of low-mass galaxies as a function of redshift. In contrast, the conventional ejective feedback model fails to reproduce these observational trends. Using the model, we demonstrate that the properties of disk galaxies are closely tied to the thermal state of hot halo gas and even possibly the circum-halo medium, which suggests that observational data for the disk properties and circum-galactic hot/warm medium may jointly provide interesting constraints for galaxy formation models.

A systematic look at the effects of radiative feedback on disc galaxy formation

Galaxy formation models and simulations rely on various feedback mechanisms to reproduce the observed baryonic scaling relations and galaxy morphologies. Although dwarf galaxy and giant elliptical properties can be explained using feedback from supernova and active galactic nuclei, Milky Way-sized galaxies still represent a challenge to current theories of galaxy formation. In this paper, we explore the possible role of feedback from stellar radiation in regulating the main properties of disc galaxies such as our own Milky Way. We have performed a suite of cosmological simulations of the same ∼1012 M⊙ halo selected based on its rather typical mass accretion history. We have implemented radiative feedback from young stars using a crude model of radiative transfer for ultraviolet and infrared radiation. However, the model is realistic enough such that the dust opacity plays a direct role in regulating the efficiency of our feedback mechanism. We have explored various models for the dust opacity, assuming different constant dust temperatures, as well as a varying dust temperature model. We find that while strong radiative feedback appears as a viable mechanism to regulate the stellar mass fraction in massive galaxies, it also prevents the formation of discs with reasonable morphologies. In models with strong stellar radiation feedback, stellar discs are systematically too thick while the gas disc morphology is completely destroyed due to the efficient mixing between the feedback-affected gas and its surroundings. At the resolution of our simulation suite, we find it impossible to preserve spiral disc morphology while at the same time expelling enough baryons to satisfy the abundance matching constraints.

GAUSSIAN RANDOM FIELD: PHYSICAL ORIGIN OF SERSIC PROFILES

While the Sersic profile family provide adequate fits for the surface brightness profiles of observed galaxies, the physical origin is unknown. We show that, if the cosmological density field are seeded by random gaussian fluctuations, as in the standard cold dark matter model, galaxies with steep central profiles have simultaneously extended envelopes of shallow profiles in the outskirts, whereas galaxies with shallow central profiles are accompanied by steep density profiles in the outskirts. These properties are in accord with those of the Sersic profile family. Moreover, galaxies with steep central profiles form their central regions in smaller denser subunits that possibly merge subsequently, which naturally leads to formation of bulges. In contrast, galaxies with shallow central profiles form their central regions in a coherent fashion without significant substructure, a necessary condition for disk galaxy formation. Thus, the scenario is self-consistent with respect to the correlation between observed galaxy morphology and Sersic index. We predict further that clusters of galaxies should display a similar trend, which should be verifiable observationally.

Diffuse gas properties and stellar metallicities in cosmological simulations of disc galaxy formation

We analyse the properties of the circum-galactic medium and the metal content of the stars comprising the central galaxy in eight hydrodynamical `zoom-in' simulations of disc galaxy formation. We use these properties as a benchmark for our model of galaxy formation physics implemented in the moving-mesh code AREPO, which succeeds in forming quite realistic late-type spirals in the set of `Aquarius' initial conditions of Milky Way-sized haloes. Galactic winds significantly influence the morphology of the circum-galactic medium and induce bipolar features in the distribution of heavy elements. They also affect the thermodynamic properties of the circum-galactic gas by supplying an energy input that sustains its radiative losses. Although a significant fraction of the heavy elements are transferred from the central galaxy to the halo, and even beyond the virial radius, enough metals are retained by stars to yield a peak in their metallicity distributions at about $Z_{\odot}$. All our default runs overestimate the stellar [O/Fe] ratio, an effect that we demonstrate can be rectified by an increase of the adopted SN type Ia rate. Nevertheless, the models have difficulty in producing stellar metallicity gradients of the same strength as observed in the Milky Way.