skip to main content


Title: Emergence of galactic morphologies at cosmic dawn: input from numerical modelling
ABSTRACT

We employ a series of high-resolution zoom-in cosmological simulations to analyse the emerging morphology of main galaxies in dark matter haloes at z ≳ 2. We choose haloes of similar masses, ${\rm log}\, M_{\rm vir}/{\rm M_\odot }\sim 11.65\pm 0.05$, at the target zf = 6, 4, and 2. The rationale for this choice allows us to analyse how the different growth rate in these haloes propagates down to galaxy scales, affecting their basic parameters. Halos were embedded in high/low overdensity regions, and two versions of a galactic wind feedback were employed. Our main results are: (1) Although our galaxies evolve in different epochs, their global parameters remain within narrow range. Their morphology, kinematics, and stellar populations differ substantially, yet all host sub-kpc stellar bars; (2) The star formation rates appear higher for larger zf; (3) Bulges and stellar spheroids were separated by stellar kinematics, discy bulges were revealed using the Sersic method and photometry.The bulge-to-total mass ratios appear independent of the last merger time for all zf. The spheroid-to-total mas ratios lie within ∼0.5–0.8; (4) The synthetic redshifted, pixelized, and PSF-degraded JWST images allow detection of stellar discs at all zf. (5) Based on the kinematic decomposition, rotational support in discs depends on the feedback type, but increases with decreasing zf; (6) Finally, the ALMA images detect discs at all zf, but spiral structure is detectable in zf = 2 galaxies. Moreover, galaxies follow the Tully–Fisher relation, being separated only by the galactic wind feedback.

 
more » « less
NSF-PAR ID:
10366598
Author(s) / Creator(s):
; ;
Publisher / Repository:
Oxford University Press
Date Published:
Journal Name:
Monthly Notices of the Royal Astronomical Society
Volume:
513
Issue:
1
ISSN:
0035-8711
Format(s):
Medium: X Size: p. 693-712
Size(s):
["p. 693-712"]
Sponsoring Org:
National Science Foundation
More Like this
  1. ABSTRACT

    Recent observations indicate that galactic outflows are ubiquitous in high-redshift (high-z) galaxies, including normal star-forming galaxies, quasar hosts, and dusty star-forming galaxies (DSFGs). However, the impact of outflows on the evolution of their hosts is still an open question. Here, we analyse the star-formation histories and galactic outflow properties of galaxies in massive haloes ($10^{12}\, {\rm M}_{\odot }\ \lt\ M_{\rm vir}\ \lt\ 5\times 10^{12}\, {\rm M}_{\odot }$) at z ≳ 5.5 in three zoom-in cosmological simulations from the MassiveFIRE suite, as part of the Feedback In Realistic Environments (FIRE) project. The simulations were run with the FIRE-2 model, which does not include feedback from active galactic nuclei. The simulated galaxies resemble z > 4 DSFGs, with star-formation rates of $\sim\!{1000}\ {\rm M}_{\odot }\, \rm yr^{-1}$ and molecular gas masses of Mmol ∼ 1010 M⊙. However, the simulated galaxies are characterized by higher circular velocities than those observed in high-z DSFGs. The mass loading factors from stellar feedback are of the order of ∼0.1, implying that stellar feedback is inefficient in driving galactic outflows and gas is consumed by star formation on much shorter time-scales than it is expelled from the interstellar medium. We also find that stellar feedback is highly inefficient in self-regulating star formation in this regime, with an average integrated star formation efficiency (SFE) per dynamical time of 30 per cent. Finally, compared with FIRE-2 galaxies hosted in similarly massive haloes at lower redshift, we find lower mass loading factors and higher SFEs in the high-z sample. We argue that both effects originate from the higher total and gas surface densities that characterize high-z massive systems.

     
    more » « less
  2. ABSTRACT We analyse the cold dark matter density profiles of 54 galaxy haloes simulated with Feedback In Realistic Environments (FIRE)-2 galaxy formation physics, each resolved within $0.5{{\ \rm per\ cent}}$ of the halo virial radius. These haloes contain galaxies with masses that range from ultrafaint dwarfs ($M_\star \simeq 10^{4.5}\, \mathrm{M}_{\odot }$) to the largest spirals ($M_\star \simeq 10^{11}\, \mathrm{M}_{\odot }$) and have density profiles that are both cored and cuspy. We characterize our results using a new, analytic density profile that extends the standard two-parameter Einasto form to allow for a pronounced constant density core in the resolved innermost radius. With one additional core-radius parameter, rc, this three-parameter core-Einasto profile is able to characterize our feedback-impacted dark matter haloes more accurately than other three-parameter profiles proposed in the literature. To enable comparisons with observations, we provide fitting functions for rc and other profile parameters as a function of both M⋆ and M⋆/Mhalo. In agreement with past studies, we find that dark matter core formation is most efficient at the characteristic stellar-to-halo mass ratio M⋆/Mhalo ≃ 5 × 10−3, or $M_{\star } \sim 10^9 \, \mathrm{M}_{\odot }$, with cores that are roughly the size of the galaxy half-light radius, rc ≃ 1−5 kpc. Furthermore, we find no evidence for core formation at radii $\gtrsim 100\ \rm pc$ in galaxies with M⋆/Mhalo < 5 × 10−4 or $M_\star \lesssim 10^6 \, \mathrm{M}_{\odot }$. For Milky Way-size galaxies, baryonic contraction often makes haloes significantly more concentrated and dense at the stellar half-light radius than DMO runs. However, even at the Milky Way scale, FIRE-2 galaxy formation still produces small dark matter cores of ≃ 0.5−2 kpc in size. Recent evidence for a ∼2 kpc core in the Milky Way’s dark matter halo is consistent with this expectation. 
    more » « less
  3. ABSTRACT

    We use numerical simulations to follow evolution of barred galaxies in a suite of models with progressively more massive stellar bulges, with bulge-to-total (disc+bulge) mass ratios of B/T ∼ 0–0.25, embedded in dark matter (DM) haloes with spin $\lambda\sim 0\!-\!0.09$. We focus on models with a sequence of initial rotational support for bulges, and analyse their spinup and spindown. We find that (1) the presence of a bulge affects evolution of bars, i.e. the time-scale of bar instability, bar pattern speed, and its decay, and the vertical buckling instability. Bar strength is nearly independent of B/T in haloes with spin $\lambda=0$, and is suppressed by a factor ∼2 for haloes with $\lambda=0.09$; (2) The main effect of the bulge is the destruction of the harmonic core which affects the buckling; (3) The bulge plays a minor role in the exchange of angular momentum between the barred disc and the DM halo, during its spinup and spindown; (4) Buckling process triggers different response above/below the disc mid-plane, which anticorrelates with the bulge mass; (5) In spinning haloes, the buckling process has a prolonged amplitude tail, extending by few Gyr, as verified by measuring distortions in the Laplace plane; (6) Furthermore, as verified by orbital spectral analysis, the bulge gains its spin from the bar mainly via the inner Lindblad resonance, while losing it via a number of resonances lying between the outer and inner Lindblad resonance. The corollary is that we do not expect to find non-rotating bulges in barred galaxies.

     
    more » « less
  4. ABSTRACT

    We quantify the impact of galaxy formation on dark matter halo shapes using cosmological simulations at redshift z = 0. Using magnetohydrodynamic simulations from the IllustrisTNG project, we focus on haloes of mass $10^{10\!-\!14} \, \rm M_{\odot }$ from the 50 Mpc (TNG50) and 100 Mpc (TNG100) boxes and compare them to dark matter-only (DMO) analogues and other simulations, e.g. Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) and Evolution and Assembly of GaLaxies and their Environments (EAGLE). We further quantify the prediction uncertainty by varying the feedback models using smaller 25 ${\rm Mpc}\, h^{-1}$ boxes. We find that (i) galaxy formation results in rounder haloes compared to DMO simulations, in qualitative agreement with past results. Haloes of mass ${\approx }2\times 10^{12} \, \rm M_{\odot }$ are most spherical, with an average minor-to-major axial ratio of $\langle s \rangle$ ≈ 0.75 in the inner halo, an increase of 40 per cent compared to their DMO counterparts. No significant difference is present for low-mass $10^{10} \, \rm M_{\odot }$ haloes; (ii) stronger feedback, e.g. increasing galactic wind speed, reduces the impact of baryons; (iii) the inner halo shape correlates with the stellar mass fraction, explaining the dependence of halo shapes on feedback models; and (iv) the fiducial and weaker feedback models are most consistent with observational estimates of the Milky Way halo shape. At fixed halo mass, very diverse and possibly unrealistic feedback models all predict inner shapes closer to one another than to the DMO results. Because of the large halo-to-halo variation in halo shape, a larger observational sample is required to statistically distinguish different baryonic prescriptions.

     
    more » « less
  5. null (Ed.)
    ABSTRACT Galaxy internal structure growth has long been accused of inhibiting star formation in disc galaxies. We investigate the potential physical connection between the growth of dispersion-supported stellar structures (e.g. classical bulges) and the position of galaxies on the star-forming main sequence at z ∼ 0. Combining the might of the SAMI and MaNGA galaxy surveys, we measure the λRe spin parameter for 3289 galaxies over $9.5 \lt \log M_{\star } [\rm {M}_{\odot }] \lt 12$. At all stellar masses, galaxies at the locus of the main sequence possess λRe values indicative of intrinsically flattened discs. However, above $\log M_{\star }[\rm {M}_{\odot }]\sim 10.5$ where the main sequence starts bending, we find tantalizing evidence for an increase in the number of galaxies with dispersion-supported structures, perhaps suggesting a connection between bulges and the bending of the main sequence. Moving above the main sequence, we see no evidence of any change in the typical spin parameter in galaxies once gravitationally interacting systems are excluded from the sample. Similarly, up to 1 dex below the main sequence, λRe remains roughly constant and only at very high stellar masses ($\log M_{\star }[\rm {M}_{\odot }]\gt 11$), do we see a rapid decrease in λRe once galaxies decline in star formation activity. If this trend is confirmed, it would be indicative of different quenching mechanisms acting on high- and low-mass galaxies. The results suggest that whilst a population of galaxies possessing some dispersion-supported structure is already present on the star-forming main sequence, further growth would be required after the galaxy has quenched to match the kinematic properties observed in passive galaxies at z ∼ 0. 
    more » « less