Modelling chemical abundance distributions for dwarf galaxies in the Local Group: the impact of turbulent metal diffusion

We investigate stellar metallicity distribution functions (MDFs), including Fe and ${\alpha}$-element abundances, in dwarf galaxies from the Feedback in Realistic Environments (FIRE) project. We examine both isolated dwarf galaxies and those that are satellites of a Milky Way-mass galaxy. In particular, we study the effects of including a sub-grid turbulent model for the diffusion of metals in gas. Simulations that include diffusion have narrower MDFs and abundance ratio distributions, because diffusion drives individual gas and star particles toward the average metallicity. This effect provides significantly better agreement with observed abundance distributions of dwarf galaxies in the Local Group, including the small intrinsic scatter in [${\alpha}$/Fe] vs. [Fe/H] (less than 0.1 dex). This small intrinsic scatter arises in our simulations because the interstellar medium (ISM) in dwarf galaxies is well-mixed at nearly all cosmic times, such that stars that form at a given time have similar abundances to within 0.1 dex. Thus, most of the scatter in abundances at z = 0 arises from redshift evolution and not from instantaneous scatter in the ISM. We find similar MDF widths and intrinsic scatter for satellite and isolated dwarf galaxies, which suggests that environmental effects play a minor role compared with internal chemical evolution in our simulations. Overall, with the inclusion of metal diffusion, our simulations reproduce abundance distribution widths of observed low-mass galaxies, enabling detailed studies of chemical evolution in galaxy formation.Comment: 19 pages, 13 figures, published in MNRA

Reconciling observed and simulated stellar halo masses

We use cosmological hydrodynamical simulations of Milky-Way-mass galaxies from the FIRE project to evaluate various strategies for estimating the mass of a galaxy's stellar halo from deep, integrated-light images. We find good agreement with integrated-light observations if we mimic observational methods to measure the mass of the stellar halo by selecting regions of an image via projected radius relative to the disk scale length or by their surface density in stellar mass . However, these observational methods systematically underestimate the accreted stellar component, defined in our (and most) simulations as the mass of stars formed outside of the host galaxy, by up to a factor of ten, since the accreted component is centrally concentrated and therefore substantially obscured by the galactic disk. Furthermore, these observational methods introduce spurious dependencies of the estimated accreted stellar component on the stellar mass and size of galaxies that can obscure the trends in accreted stellar mass predicted by cosmological simulations, since we find that in our simulations the size and shape of the central galaxy is not strongly correlated with the assembly history of the accreted stellar halo. This effect persists whether galaxies are viewed edge-on or face-on. We show that metallicity or color information may provide a way to more cleanly delineate in observations the regions dominated by accreted stars. Absent additional data, we caution that estimates of the mass of the accreted stellar component from single-band images alone should be taken as lower limits.Comment: Version accepted by Ap

Simulating galaxies in the reionization era with FIRE-2: morphologies and sizes

We study the morphologies and sizes of galaxies at z>5 using high-resolution cosmological zoom-in simulations from the Feedback In Realistic Environments project. The galaxies show a variety of morphologies, from compact to clumpy to irregular. The simulated galaxies have more extended morphologies and larger sizes when measured using rest-frame optical B-band light than rest-frame UV light; sizes measured from stellar mass surface density are even larger. The UV morphologies are usually dominated by several small, bright young stellar clumps that are not always associated with significant stellar mass. The B-band light traces stellar mass better than the UV, but it can also be biased by the bright clumps. At all redshifts, galaxy size correlates with stellar mass/luminosity with large scatter. The half-light radii range from 0.01 to 0.2 arcsec (0.05-1 kpc physical) at fixed magnitude. At z>5, the size of galaxies at fixed stellar mass/luminosity evolves as (1+z)^{-m}, with m~1-2. For galaxies less massive than M_star~10^8 M_sun, the ratio of the half-mass radius to the halo virial radius is ~10% and does not evolve significantly at z=5-10; this ratio is typically 1-5% for more massive galaxies. A galaxy's "observed" size decreases dramatically at shallower surface brightness limits. This effect may account for the extremely small sizes of z>5 galaxies measured in the Hubble Frontier Fields. We provide predictions for the cumulative light distribution as a function of surface brightness for typical galaxies at z=6.Comment: 11 pages, 11 figures, resubmitted to MNRAS after revision for referee's comment

Formation of Globular Cluster Candidates in Merging Proto-galaxies at High Redshift: A View from the FIRE Cosmological Simulations

Using a state-of-the-art cosmological simulation of merging proto-galaxies at high redshift from the FIRE project, with explicit treatments of star formation and stellar feedback in the interstellar medium, we investigate the formation of star clusters and examine one of the formation hypothesis of present-day metal-poor globular clusters. We find that frequent mergers in high-redshift proto-galaxies could provide a fertile environment to produce long-lasting bound star clusters. The violent merger event disturbs the gravitational potential and pushes a large gas mass of ~> 1e5-6 Msun collectively to high density, at which point it rapidly turns into stars before stellar feedback can stop star formation. The high dynamic range of the reported simulation is critical in realizing such dense star-forming clouds with a small dynamical timescale, t_ff <~ 3 Myr, shorter than most stellar feedback timescales. Our simulation then allows us to trace how clusters could become virialized and tightly-bound to survive for up to ~420 Myr till the end of the simulation. Because the cluster's tightly-bound core was formed in one short burst, and the nearby older stars originally grouped with the cluster tend to be preferentially removed, at the end of the simulation the cluster has a small age spread.Comment: 14 pages, 14 figures, Accepted for publication in the Monthly Notices of the Royal Astronomical Society, High-resolution version of this article also available at http://www.jihoonkim.org/index/research.html#g

Dark Matter Halo Growth II: Diffuse Accretion and its Environmental Dependence

Dark matter haloes in Lambda CDM simulations grow by mergers with other haloes as well as accretion of "diffuse" non-halo material. We quantify the mass growth rates via these two processes, dM_mer/dt and dM_dif/dt, and their dependence on halo environment using the ~500,000 haloes in the Millennium simulation. Adopting a local mass density parameter as a measure of halo environment, we find the two rates show strong but opposite environmental dependence, with mergers playing an increasingly important role for halo growths in overdense regions and diffuse accretion dominating growth in voids. This behaviour is independent of the mass cuts used to define haloes vs non-haloes. For galaxy-scale haloes, these two opposite correlations largely cancel out, but a weak environmental dependence remains that results in a slightly lower mean total growth rate, and hence an earlier mean formation redshift, for haloes in denser regions. The mean formation redshift of the ~5000 cluster-mass haloes, on the other hand, appears to have no correlation with halo environment. The origin of the positive correlation of dM_mer/dt with local density can be traced to the surrounding mass reservoir outside the haloes, where more progenitor haloes are available in denser regions. The negative correlation of dM_dif/dt with density, however, is not explained by the available diffuse mass in the reservoir, which is in fact larger in denser regions. The non-halo component may therefore be partially comprised of truly diffuse dark matter particles that are dynamically hotter and are accreted at a suppressed rate in denser regions. We also discuss the implications of these results for how to modify the Extended Press-Schechter model of halo growth, which in its original form does not predict environmental dependence.Comment: 15 pages, 9 figures, accepted in MNRA

Fibril polymorphism affects immobilized non-amyloid flanking domains of huntingtin exon1 rather than its polyglutamine core

Polyglutamine expansion in the huntingtin protein is the primary genetic cause of Huntington's disease (HD). Fragments coinciding with mutant huntingtin exon1 aggregate in vivo and induce HD-like pathology in mouse models. The resulting aggregates can have different structures that affect their biochemical behaviour and cytotoxic activity. Here we report our studies of the structure and functional characteristics of multiple mutant htt exon1 fibrils by complementary techniques, including infrared and solid-state NMR spectroscopies. Magic-angle-spinning NMR reveals that fibrillar exon1 has a partly mobile α-helix in its aggregation-accelerating N terminus, and semi-rigid polyproline II helices in the proline-rich flanking domain (PRD). The polyglutamine-proximal portions of these domains are immobilized and clustered, limiting access to aggregation-modulating antibodies. The polymorphic fibrils differ in their flanking domains rather than the polyglutamine amyloid structure. They are effective at seeding polyglutamine aggregation and exhibit cytotoxic effects when applied to neuronal cells

The Mass Spectrum of Light and Heavy Hadrons from Improved Lattice Actions

We use improved lattice actions for glue, light quarks and heavy quarks for which we use lattice NRQCD to compute hadron masses. Our results are in good agreement with experiment, except for charmed hadrons. It seems that charmed quar ks are not well approximated as heavy quarks nor as light quarks.Comment: 14 pages +6 pages figures, plain-tex fil

The Merger Rates and Mass Assembly Histories of Dark Matter Haloes in the Two Millennium Simulations

We construct merger trees of dark matter haloes and quantify their merger rates and mass growth rates using the joint dataset from the Millennium and Millennium-II simulations. The finer resolution of the Millennium-II Simulation has allowed us to extend our earlier analysis of halo merger statistics to an unprecedentedly wide range of descendant halo mass (10^10 < M0 < 10^15 Msun), progenitor mass ratio (10^-5 < xi < 1), and redshift (0 < z < 15). We update our earlier fitting form for the mean merger rate per halo as a function of M_0, xi, and z. The overall behavior of this quantity is unchanged: the rate per unit redshift is nearly independent of z out to z~15; the dependence on halo mass is weak (M0^0.13); and it is nearly a power law in the progenitor mass ratio (xi^-2). We also present a simple and accurate fitting formula for the mean mass growth rate of haloes as a function of mass and redshift. This mean rate is 46 Msun/yr for 10^12 Msun haloes at z=0, and it increases with mass as M^{1.1} and with redshift as (1+z)^2.5 (for z > 1). When the fit for the mean mass growth rate is integrated over a halo's history, we find excellent match to the mean mass assembly histories of the simulated haloes. By combining merger rates and mass assembly histories, we present results for the number of mergers over a halo's history and the statistics of the redshift of the last major merger.Comment: 12 pages, 9 figures, accepted in MNRA

FIRE-2 Simulations: Physics versus Numerics in Galaxy Formation

The Feedback In Realistic Environments (FIRE) project explores feedback in cosmological galaxy formation simulations. Previous FIRE simulations used an identical source code (“FIRE-1”) for consistency. Motivated by the development of more accurate numerics – including hydrodynamic solvers, gravitational softening, and supernova coupling algorithms – and exploration of new physics (e.g. magnetic fields), we introduce “FIRE-2”, an updated numerical implementation of FIRE physics for the GIZMO code. We run a suite of simulations and compare against FIRE-1: overall, FIRE-2 improvements do not qualitatively change galaxy-scale properties. We pursue an extensive study of numerics versus physics. Details of the star-formation algorithm, cooling physics, and chemistry have weak effects, provided that we include metal-line cooling and star formation occurs at higher-than-mean densities. We present new resolution criteria for high-resolution galaxy simulations. Most galaxy-scale properties are robust to numerics we test, provided: (1) Toomre masses are resolved; (2) feedback coupling ensures conservation, and (3) individual supernovae are time-resolved. Stellar masses and profiles are most robust to resolution, followed by metal abundances and morphologies, followed by properties of winds and circum-galactic media (CGM). Central (∼kpc) mass concentrations in massive (>L*) galaxies are sensitive to numerics (via trapping/recycling of winds in hot halos). Multiple feedback mechanisms play key roles: supernovae regulate stellar masses/winds; stellar mass-loss fuels late star formation; radiative feedback suppresses accretion onto dwarfs and instantaneous star formation in disks. We provide all initial conditions and numerical algorithms used