Search

We present a suite of high-resolution simulations of an isolated dwarf galaxy using four different hydrodynamical codes: {\sc Gizmo}, {\sc Arepo}, {\sc Gadget}, and {\sc Ramses}. All codes adopt the same physical model which includes radiative cooling, photoelectric heating, star formation, and supernova (SN) feedback. Individual SN explosions are directly resolved without resorting to sub-grid models, eliminating one of the major uncertainties in cosmological simulations. We find reasonable agreement on the time-averaged star formation rates as well as the joint density-temperature distributions between all codes. However, the Lagrangian codes show significantly burstier star formation, larger supernova-driven bubbles, and stronger galactic outflows compared to the Eulerian code. This is caused by the behavior in the dense, collapsing gas clouds when the Jeans length becomes unresolved: gas in Lagrangian codes collapses to much higher densities than in Eulerian codes, as the latter is stabilized by the minimal cell size. Therefore, more of the gas cloud is converted to stars and SNe are much more clustered in the Lagrangian models, amplifying their dynamical impact. The differences between Lagrangian and Eulerian codes can be reduced by adopting a higher star formation efficiency in Eulerian codes, which significantly enhances SN clustering in the latter. Adopting a zero SN delay time reduces burstiness in all codes, resulting in vanishing outflows as SN clustering is suppressed., Comment: accepted version by ApJ (including a new simulation in Appendix B suggested by the referee)

Apparent deviations between properties of dwarf galaxies from observations and simulations are known to exist, such as the "Missing Dwarfs" problem, the too-big-to-fail problem, and the cusp-core problem, to name a few. Recent studies have shown that these issues can at least be partially resolved by taking into account the systematic differences between simulations and observations. This work aims to investigate and address any systematic differences affecting the comparison of simulations with observations. To this aim, we analyzed a set of 24 realistically simulated MoRIA (Models of Realistic dwarfs In Action) dwarf galaxies in an observationally motivated way. We first constructed "observed" color-magnitude diagrams (CMDs) of the simulated dwarf galaxies in the typically used V- and I-bands. Then we used the CMD-fitting method to recover their star-formation histories (SFHs) from their observed CMDs. These solved SFHs were then directly compared to the true SFHs from the simulation star-particle data, mainly in terms of the star-formation rate(SFR) and the age-metallicity relation (AMR). We applied a dust extinction prescription to the simulation data to produce observed CMDs affected by dust in star-formation regions. Since future facilities, such as the JWST and E-ELT will focus on the near IR rather than the optical, we also constructed and analyzed CMDs using the I- and H-bands. We find a very good agreement between the recovered and the true SFHs of all the simulated dwarf galaxies in our sample, from the synthetic CMD analysis of their V-I versus I as well as the I-H versus H CMDs. Dust leads to an underestimation of the SFR during the last few hundred million years. Overall, our analysis indicates that quantities like SFR and AMR derived from the photometric observations of galaxies are directly comparable to their simulated counterparts., Comment: 19 pages, 15 figures, to appear in Section 4. Extragalactic astronomy of Astronomy and Astrophysics

We introduce the "Engineering Dwarfs at Galaxy formation's Edge" (EDGE) project to study the cosmological formation and evolution of the smallest galaxies in the Universe. In this first paper, we explore the effects of resolution and sub-grid physics on a single low mass halo ($M_{\rm halo}=10^{9}~M_\odot$), simulated to redshift $z=0$ at a mass and spatial resolution of $\sim 20~M_\odot$ and $\sim 3$ pc. We consider different star formation prescriptions, supernova feedback strengths and on-the-fly radiative transfer (RT). We show that RT changes the mode of galactic self-regulation at this halo mass, suppressing star formation by causing the interstellar and circumgalactic gas to remain predominantly warm ($\sim 10^4$ K) even before cosmic reionisation. By contrast, without RT, star formation regulation occurs only through starbursts and their associated vigorous galactic outflows. In spite of this difference, the entire simulation suite (with the exception of models without any feedback) matches observed dwarf galaxy sizes, velocity dispersions, $V$-band magnitudes and dynamical mass-to-light-ratios. This is because such structural scaling relations are predominantly set by the host dark matter halo, with the remaining model-to-model variation being smaller than the observational scatter. We find that only the stellar mass-metallicity relation differentiates the galaxy formation models. Explosive feedback ejects more metals from the dwarf, leading to a lower metallicity at a fixed stellar mass. We conclude that the stellar mass-metallicity relation of the very smallest galaxies provides a unique constraint on galaxy formation physics., Comment: 17 pages, 6 figures, Minor changes to mach published MNRAS version

We test whether advanced galaxy models and analysis techniques of simulations can alleviate the Too Big To Fail problem (TBTF) for late-type galaxies, which states that isolated dwarf galaxy kinematics imply that dwarfs live in lower-mass halos than is expected in a {\Lambda}CDM universe. Furthermore, we want to explain this apparent tension between theory and observations. To do this, we use the MoRIA suite of dwarf galaxy simulations to investigate whether observational effects are involved in TBTF for late-type field dwarf galaxies. To this end, we create synthetic radio data cubes of the simulated MoRIA galaxies and analyse their HI kinematics as if they were real, observed galaxies. We find that for low-mass galaxies, the circular velocity profile inferred from the HI kinematics often underestimates the true circular velocity profile, as derived directly from the enclosed mass. Fitting the HI kinematics of MoRIA dwarfs with a theoretical halo profile results in a systematic underestimate of the mass of their host halos. We attribute this effect to the fact that the interstellar medium of a low-mass late-type dwarf is continuously stirred by supernova explosions into a vertically puffed-up, turbulent state to the extent that the rotation velocity of the gas is simply no longer a good tracer of the underlying gravitational force field. If this holds true for real dwarf galaxies as well, it implies that they inhabit more massive dark matter halos than would be inferred from their kinematics, solving TBTF for late-type field dwarf galaxies., Comment: 21 pages, 21 figures. Accepted for publication in A&A. Corrected certain values in Table 1

Simulating dwarf galaxy halos in a reionizing Universe puts severe constraints on the sub-grid model employed in the simulations. Using the same sub-grid model that works for simulations without a UV-background (UVB) results in gas poor galaxies that stop forming stars very early on, except for halos with high masses. This is in strong disagreement with observed galaxies, which are gas rich and star forming down to a much lower mass range. To resolve this discrepancy, we ran a large suite of isolated dwarf galaxy simulations to explore a wide variety of sub-grid models and parameters, including timing and strength of the UVB, strength of the stellar feedback, and metallicity dependent Pop III feedback. We compared these simulations to observed dwarf galaxies by means of the baryonic Tully-Fisher relation (BTFR), which links the baryonic content of a galaxy to the observationally determined strength of its gravitational potential. We found that the results are robust to changes in the UVB. The strength of the stellar feedback shifts the results on the BTFR, but does not help to form gas rich galaxies at late redshifts. Only by including Pop III feedback are we able to produce galaxies that lie on the observational BTFR and that have neutral gas and ongoing star formation at redshift zero., Comment: Accepted for publication in MNRAS. 25 pages, 2 tables and 36 figures. Interactive plots can be found on http://www.dwarfs.ugent.be/btfr/

Low-mass dwarf galaxies are very sensitive test-beds for theories of cosmic structure formation since their weak gravitational fields allow the effects of the relevant physical processes to clearly stand out. Up to now, no unified account exists of the sometimes seemingly conflicting properties of the faintest isolated dwarfs in and around the Local Group, such as Leo T and the recently discovered Leo P and Pisces A systems. Using new numerical simulations, we show that this serious challenge to our understanding of galaxy formation can be effectively resolved by taking into account the regulating influence of the ultraviolet radiation of the first population of stars on a dwarf's star formation rate while otherwise staying within the standard cosmological paradigm for structure formation. These simulations produce faint, gas-dominated, star-forming dwarf galaxies that lie on the baryonic Tully-Fisher relation and that successfully reproduce a broad range of chemical, kinematical, and structural observables of real late-type dwarf galaxies. Furthermore, we stress the importance of obtaining properties of simulated galaxies in a manner as close as possible to the typically employed observational techniques., Comment: 13 pages, 2 tables, 12 figures. Accepted for publication in ApJ

We show that the general relativistic theory of the dynamics of isotropic stellar clusters can be developed essentially along the same lines as the Newtonian theory. We prove that the distribution function can be derived from any isotropic momentum moment and that every higher-order moment of the distribution can be written as an integral over a zeroth-order moment. We propose a mathematically simple expression for the distribution function of a family of isotropic general relativistic cluster models and investigate their dynamical properties. In the Newtonian limit, these models obtain a distribution function of the form F(E) ~ (E-E_0)^alpha, with E binding energy and E_0 a constant that determines the model's outer radius. The slope alpha sets the steepness of the distribution function and the corresponding radial density and pressure profiles. We show that the field equations only yield solutions with finite mass for alpha<3.5. Moreover, in the limit alpha->3.5, only Newtonian models exist. In other words: within the context of this family of models, no general relativistic version of the Plummer model exists. The most strongly bound model within the family is characterized by alpha=2.75 and a central redshift z_c~0.55., Comment: 10 pages, 5 figures, accepted for publication by MNRAS

We investigate the evolution of dwarf galaxies using N -body/SPH simulations that incorporate their formation histories through merger trees constructed using the ex- tended Press-Schechter formalism. The simulations are computationally cheap and have high spatial resolution. We compare the properties of galaxies with equal final mass but with different merger histories with each other and with those of observed dwarf spheroidals and irregulars. We show that the merger history influences many observable dwarf galaxy proper- ties. We identify two extreme cases that make this influence stand out most clearly: (i) merger trees with one massive progenitor that grows through relatively few mergers and (ii) merger trees with many small progenitors that merge only quite late. At a fixed halo mass, a type (i) tree tends to produce galaxies with larger stellar masses, larger half-light radii, lower central surface brightness, and since fewer potentially an- gular momentum cancelling mergers are required to build up the final galaxy, a higher specific angular momentum, compared with a type (ii) tree. We do not perform full-fledged cosmological simulations and therefore cannot hope to reproduce all observed properties of dwarf galaxies. However, we show that the simulated dwarfs are not unsimilar to real ones., Comment: 19 pages, 17 figures, 3 tables

Using computer simulations, we explored gaseous infall as a possible explanation for the starburst phase in Blue Compact Dwarf galaxies. We simulate a cloud impact by merging a spherical gas cloud into an isolated dwarf galaxy. We investigated which conditions were favourable for triggering a burst and found that the orbit and the mass of the gas cloud play an important role. We discuss the metallicity, the kinematical properties, the internal dynamics and the gas, stellar and dark matter distribution of the simulations during a starburst. We find that these are in good agreement with observations and depending on the set-up (e.g. rotation of the host galaxy, radius of the gas cloud), our bursting galaxies can have qualitatively very different properties. Our simulations offer insight in how starbursts start and evolve. Based on this, we propose what postburst dwarf galaxies will look like., Comment: Accepted for publication in MNRAS | 16 pages, 16 figures