Constraining the subgrid physics in simulations of isolated dwarf galaxies

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

How the first stars shaped the faintest gas-dominated dwarf galaxies

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 Ap

The dynamics of general relativistic isotropic stellar cluster models -- Do relativistic extensions of the Plummer model exist?

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 alpha3.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 MNRA

Numerical simulations of dwarf galaxy merger trees

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 table

Gaseous infall triggering starbursts in simulated dwarf galaxies

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 figure

A new astrophysical solution to the Too Big To Fail problem - Insights from the MoRIA simulations

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

How the first stars shaped the faintest gas-dominated dwarf galaxies

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 has existed 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

Code Comparison in Galaxy Scale Simulations with Resolved Supernova Feedback: Lagrangian vs. Eulerian Methods

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