We present an updated version of SKIRT, a 3D Monte Carlo radiative transfer code developed to simulate dusty galaxies. The main novel characteristics of the SKIRT code are the use of a stellar foam to generate random positions, an efficient combination of eternal forced scattering and continuous absorption, and a new library approach that links the radiative transfer code to the DustEM dust emission library. This approach enables a fast, accurate and self-consistent calculation of the dust emission of arbitrary mixtures of transiently heated dust grains and polycyclic aromatic hydrocarbons, even for full 3D models containing millions of dust cells. We have demonstrated the accuracy of the SKIRT code through a set of simulations based on the edge-on spiral galaxy UGC 4754. The models we ran were gradually refined from a smooth, 2D, LTE model to a fully 3D model that includes NLTE dust emission and a clumpy structure of the dusty ISM. We find that clumpy models absorb UV and optical radiation less efficiently than smooth models with the same amount of dust, and that the dust in clumpy models is on average both cooler and less luminous. Our simulations demonstrate that, given the appropriate use of optimization techniques, it is possible to efficiently and accurately run Monte Carlo radiative transfer simulations of arbitrary 3D structures of several million dust cells, including a full calculation of the NLTE emission by arbitrary dust mixtures.

We show results based on a large suite of N-Body/SPH simulations of isolated, flat dwarf galaxies, both rotating and non-rotating. The main goal is to investigate possible mechanisms to explain the observed dichotomy in radial stellar metallicity profiles of dwarf galaxies: dwarf irregulars (dIrr) and flat, rotating dwarf ellipticals (dE) generally possess flat metallicity profiles, while rounder and non-rotating dEs show strong negative metallicity gradients. These simulations show that flattening by rotation is key to reproducing the observed characteristics of flat dwarf galaxies, proving particularly efficient in erasing metallicity gradients. We propose a "centrifugal barrier mechanism" as an alternative to the previously suggested "fountain mechanism" for explaining the flat metallicity profiles of dIrrs and flat, rotating dEs. While only flattening the dark-matter halo has little influence, the addition of angular momentum slows down the infall of gas, so that star formation (SF) and the ensuing feedback are less centrally concentrated, occurring galaxy-wide. Additionally, this leads to more continuous SFHs by preventing large-scale oscillations in the SFR ("breathing"), and creates low density holes in the ISM, in agreement with observations of dIrrs. Our general conclusion is that rotation has a significant influence on the evolution and appearance of dwarf galaxies, and we suggest angular momentum as a second parameter (after galaxy mass as the dominant parameter) in dwarf galaxy evolution. Angular momentum differentiates between SF modes, making our fast rotating models qualitatively resemble dIrrs, which does not seem possible without rotation.

In this paper we investigate whether Smoothed Particle Hydrodynamics (SPH), equipped with artificial conductivity, is able to capture the physics of density/energy discontinuities in the case of the so-called shearing layers test, a test for examining Kelvin-Helmholtz (KH) instabilities. We can trace back each failure of SPH to show KH rolls to two causes: i) shock waves travelling in the simulation box and ii) particle clumping, or more generally, particle noise. The probable cause of shock waves is the Local Mixing Instability (LMI), previously identified in the literature. Particle noise on the other hand is a problem because it introduces a large error in the SPH momentum equation. We also investigate the role of artificial conductivity (AC). Including AC is necessary for the long-term behavior of the simulation (e.g. to get $\lambda=1/2, 1$ KH rolls). In sensitive hydrodynamical simulations great care is however needed in selecting the AC signal velocity, with the default formulation leading to too much energy diffusion. We present new signal velocities that lead to less diffusion. The effects of the shock waves and of particle disorder become less important as the time-scale of the physical problem (for the shearing layers problem: lower density contrast and higher Mach numbers) decreases. At the resolution of current galaxy formation simulations mixing is probably not important. However, mixing could become crucial for next-generation simulations.

We investigate the global photometric scaling relations traced by early-type galaxies in different environments, ranging from dwarf spheroidals, over dwarf elliptical galaxies, up to giant ellipticals (-8 mag > M_V > -24 mag). These results are based in part on our new HST/ACS F555W and F814W imagery of dwarf spheroidal galaxies in the Perseus Cluster. These scaling relations are almost independent of environment, with Local Group and cluster galaxies coinciding in the various diagrams. We show that at M_V ~ -14 mag, the slopes of the photometric scaling relations involving the Sersic parameters change significantly. We argue that these changes in slope reflect the different physical processes that dominate the evolution of early-type galaxies in different mass regimes. As such, these scaling relations contain a wealth of information that can be used to test models for the formation of early-type galaxies.

We present new fully self-consistent models of the formation and evolution of isolated dwarf galaxies. We have used the publicly available N-body/SPH code HYDRA, to which we have added a set of star formation criteria, and prescriptions for chemical enrichment (taking into account contributions from both SNIa and SNII), supernova feedback, and gas cooling. The models follow the evolution of an initially homogeneous gas cloud collapsing in a pre-existing dark-matter halo. These simplified initial conditions are supported by the merger trees of isolated dwarf galaxies extracted from the milli-Millennium Simulation. The star-formation histories of the model galaxies exhibit burst-like behaviour. These bursts are a consequence of the blow-out and subsequent in-fall of gas. The amount of gas that leaves the galaxy for good is found to be small, in absolute numbers, ranging between 3x10^7 Msol and 6x10^7 Msol . For the least massive models, however, this is over 80 per cent of their initial gas mass. The local fluctuations in gas density are strong enough to trigger star-bursts in the massive models, or to inhibit anything more than small residual star formation for the less massive models. Between these star-bursts there can be time intervals of several Gyrs. We have compared model predictions with available data for the relations between luminosity and surface brightness profile, half-light radius, central velocity dispersion, broad band colour (B-V) and metallicity, as well as the location relative to the fundamental plane. The properties of the model dwarf galaxies agree quite well with those of observed dwarf galaxies.