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.

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.

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.

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.

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.

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.

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.

Blue compact dwarf galaxies (BCDs) form stars at, for their sizes, extraordinarily high rates. In this paper, we study what triggers this starburst and what is the fate of the galaxy once its gas fuel is exhausted. We select four BCDs with smooth outer regions, indicating them as possible progenitors of dwarf elliptical galaxies. We have obtained photometric and spectroscopic data with the FORS and ISAAC instruments on the VLT. We analyse their infra-red spectra using a full spectrum fitting technique which yields the kinematics of their stars and ionized gas together with their stellar population characteristics. We find that the_stellar_ velocity to velocity dispersion ratio (v/sigma) of our BCDs is of the order of 1.5, similar to that of dwarf elliptical galaxies. Thus, those objects do not require significant (if any) loss of angular momentum to fade into early type dwarfs. This finding is in discordance with previous studies, which however compared the stellar kinematics of dwarf elliptical galaxies with the gaseous kinematics of star forming dwarfs. The stellar velocity fields of our objects are very disturbed and the star-formation regions are often kinematically decoupled from the rest of the galaxy. These regions can be more or less metal rich with respect to the galactic body, and sometimes they are long lived. These characteristics prevent us from pinpointing a unique trigger of the star formation, even within the same galaxy. Gas impacts, mergers, and in-spiraling gas clumps are all possible star-formation ignitors for our targets.