Recent integral field spectroscopic (IFS) surveys have revealed radial gradients in the optical spectral indices of post-starburst galaxies, which can be used to constrain their formation histories. We study the spectral indices of post-processed mock IFS datacubes of binary merger simulations, carefully matched to the properties of the MaNGA IFS survey, with a variety of black hole feedback models, progenitor galaxies, orbits and mass ratios. Based on our simulation sample, we find that only major mergers on prograde-prograde or retrograde-prograde orbits in combination with a mechanical black hole feedback model can form galaxies with weak enough ongoing star formation, and therefore absent H$\alpha$ emission, to be selected by traditional PSB selection methods. We find strong fluctuations in nebular emission line strengths, even within the PSB phase, suggesting that H$\alpha$ selected PSBs are only a subsample of the underlying population. The global PSB population can be more robustly identified using stellar continuum-based approaches. The difficulty in reproducing the very young PSBs in simulations potentially indicates that new sub-resolution star formation recipes are required to properly model the process of star formation quenching. In our simulations, we find that the starburst peaks at the same time at all radii, but is stronger and more prolonged in the inner regions. This results in a strong time evolution in the radial gradients of the spectral indices which can be used to estimate the age of the starburst without reliance on detailed star formation histories from spectral synthesis models.

We investigate the connection between galaxy--galaxy mergers and enhanced black hole (BH) growth using the cosmological hydrodynamical EAGLE simulation. We do this via three methods of analysis, investigating: the merger fraction of AGN, the AGN fraction of merging systems and the AGN fraction of galaxies with close companions. In each case, we find an increased abundance of AGN within merging systems relative to control samples of inactive or isolated galaxies (by up to a factor of $\approx 3$ depending on the analysis method used), confirming that mergers are enhancing BH accretion rates for at least a subset of the galaxy population. The greatest excess of AGN triggered via a merger are found in lower mass ($M_* \sim 10^{10}$~\Msol) gas rich ($f_{\mathrm{gas}} > 0.2$) central galaxies with lower mass BHs ($M_{\mathrm{BH}} \sim 10^{7}$~\Msol) at lower redshifts ($z<1$). We find no enhancement of AGN triggered via mergers in more massive galaxies ($M_* \gtrsim 10^{11}$~\Msol). The enhancement of AGN is not uniform throughout the phases of a merger, and instead peaks within the early \emph{remnants} of merging systems (typically lagging $\approx 300$~Myr post-coalescence of the two galaxies at $z=0.5$). We argue that neither major ($M_{\mathrm{*,1}} / M_{\mathrm{*,2}} \geq \frac{1}{4}$) nor minor mergers ($\frac{1}{10} < M_{\mathrm{*,1}} / M_{\mathrm{*,2}} < \frac{1}{4}$) are statistically relevant for enhancing BH masses globally. Whilst at all redshifts the galaxies experiencing a merger have accretion rates that are on average 2--3 times that of isolated galaxies, the majority of mass that is accreted onto BHs occurs outside the periods of a merger. We compute that on average no more than 15\% of a BHs final day mass comes from the enhanced accretion rates triggered via a merger

We study structure formation in a set of cosmological simulations to uncover the scales in the initial density field that gave rise to the formation of present-day structures. Our simulations share a common primordial power spectrum (here Lambda-CDM), but the introduction of hierarchical variations of the phase information allows us to systematically study the scales that determine the formation of structure at later times. We consider the variance in z=0 statistics such as the matter power spectrum and halo mass function. We also define a criterion for the existence of individual haloes across simulations, and determine what scales in the initial density field contain sufficient information for the non-linear formation of unique haloes. We study how the characteristics of individual haloes such as the mass and concentration, as well as the position and velocity, are affected by variations on different scales, and give scaling relations for haloes of different mass. Finally, we use the example of a cluster-mass halo to show how our hierarchical parametrisation of the initial density field can be used to create variants of particular objects. With properties such as mass, concentration, kinematics and substructure of haloes set on distinct and well-determined scales, and its unique ability to introduce variations localised in real space, our method is a powerful tool to study structure formation in cosmological simulations.

We describe a population of young star clusters (SCs) formed in a hydrodynamical simulation of a gas-rich dwarf galaxy merger resolved with individual massive stars at sub-parsec spatial resolution. The simulation is part of the GRIFFIN (Galaxy Realizations Including Feedback From INdividual massive stars) project. The star formation environment during the simulation spans seven orders of magnitude in gas surface density and thermal pressure, and the global star formation rate surface density ($\Sigma_\mathrm{SFR}$) varies by more than three orders of magnitude during the simulation. Young SCs more massive than $M_{\mathrm{*,cl}}\sim 10^{2.5}\,M_{\odot}$ form along a mass function with a power-law index $\alpha\sim-1.7$ ($\alpha\sim-2$ for $M_{\mathrm{*,cl}}\gtrsim10^{3}\,M_{\odot}$) at all merger phases, while the normalization and the highest SC masses (up to $\sim 10^6 M_{\odot}$) correlate with $\Sigma_\mathrm{SFR}$. The cluster formation efficiency varies from $\Gamma\sim20\%$ in early merger phases to $\Gamma\sim80\%$ at the peak of the starburst and is compared to observations and model predictions. The massive SCs ($\gtrsim10^4\,M_{\odot}$) have sizes and mean surface densities similar to observed young massive SCs. Simulated lower mass clusters appear slightly more concentrated than observed. All SCs form on timescales of a few Myr and lose their gas rapidly resulting in typical stellar age spreads between $\sigma\sim0.1-2$ Myr ($1\sigma$), consistent with observations. The age spreads increase with cluster mass, with the most massive cluster ($\sim10^6\, M_{\odot}$) reaching a spread of $5\, \mathrm{Myr}$ once its hierarchical formation finishes. Our study shows that it is now feasible to investigate the SC population of entire galaxies with novel high-resolution numerical simulations.

We present the novel algorithmically regularised integration method MSTAR for high accuracy ($|\Delta E/E| \gtrsim 10^{-14}$) integrations of N-body systems using minimum spanning tree coordinates. The two-fold parallelisation of the $\mathcal{O}(N_\mathrm{part}^2)$ force loops and the substep divisions of the extrapolation method allows for a parallel scaling up to $N_\mathrm{CPU} = 0.2 \times N_\mathrm{part}$. The efficient parallel scaling of MSTAR makes the accurate integration of much larger particle numbers possible compared to the traditional algorithmic regularisation chain (AR-CHAIN) methods, e.g. $N_\mathrm{part} = 5000$ particles on $400$ CPUs for $1$ Gyr in a few weeks of wall-clock time. We present applications of MSTAR on few particle systems, studying the Kozai mechanism and N-body systems like star clusters with up to $N_\mathrm{part} =10^4$ particles. Combined with a tree or a fast multipole based integrator the high performance of MSTAR removes a major computational bottleneck in simulations with regularised subsystems. It will enable the next generation galactic-scale simulations with up to $10^9$ stellar particles (e.g. $m_\star = 100 M_\odot$ for a $M_\star = 10^{11} M_\odot$ galaxy) including accurate collisional dynamics in the vicinity of nuclear supermassive black holes.

We study the orbital evolution and gravitational wave (GW) emission of supermassive black hole (SMBH) binaries formed in gas-free mergers of massive early-type galaxies using the hybrid tree-regularized N-body code KETJU. The evolution of the SMBHs and the surrounding galaxies is followed self-consistently from the large-scale merger down to the final few orbits before the black holes coalesce. Post-Newtonian corrections are included up to PN3.5-level for the binary dynamics, and the GW calculations include the corresponding corrections up to PN1.0-level. We analyze the significance of the stellar environment on the evolution of the binary and the emitted GW signal during the final GW emission dominated phase of the binary hardening and inspiral. Our simulations are compared to semi-analytic models that have often been used for making predictions for the stochastic GW background emitted by SMBHs. We find that the commonly used semi-analytic parameter values produce large differences in merger timescales and eccentricity evolution, but result in only $\sim 10\%$ differences in the GW spectrum emitted by a single binary at frequencies $f\gtrsim 10^{-1} \, \rm yr^{-1}$, which are accessible by current pulsar timing arrays. These differences are in part caused by the strong effects of the SMBH binaries on the surrounding stellar population, which are not included in the semi-analytic models.

Interpolation of data represented in curvilinear coordinates and possibly having some non-trivial, typically Riemannian or semi-Riemannian geometry is an ubiquitous task in all of physics. In this work we present a covariant generalization of the barycentric coordinates and the barycentric interpolation method for Riemannian and semi-Riemannian spaces of arbitrary dimension. We show that our new method preserves the linear accuracy property of barycentric interpolation in a coordinate-invariant sense. In addition, we show how the method can be used to interpolate constrained quantities so that the given constraint is automatically respected. We showcase the method with two astrophysics related examples situated in the curved Kerr spacetime. The first problem is interpolating a locally constant vector field, in which case curvature effects are expected to be maximally important. The second example is a General Relativistic Magnetohydrodynamics simulation of a turbulent accretion flow around a black hole, wherein high intrinsic variability is expected to be at least as important as curvature effects.

We propose a space-based interferometer surveying the gravitational wave (GW) sky in the milli-Hz to $\mu$-Hz frequency range. By the 2040s', the $\mu$-Hz frequency band, bracketed in between the Laser Interferometer Space Antenna (LISA) and pulsar timing arrays, will constitute the largest gap in the coverage of the astrophysically relevant GW spectrum. Yet many outstanding questions related to astrophysics and cosmology are best answered by GW observations in this band. We show that a $\mu$-Hz GW detector will be a truly overarching observatory for the scientific community at large, greatly extending the potential of LISA. Conceived to detect massive black hole binaries from their early inspiral with high signal-to-noise ratio, and low-frequency stellar binaries in the Galaxy, this instrument will be a cornerstone for multimessenger astronomy from the solar neighbourhood to the high-redshift Universe.

We present a hydro-dynamical simulation at sub-parsec and few solar mass resolution of a merger between two gas-rich dwarf galaxies. Our simulation includes a detailed model for the multi-phase interstellar medium (ISM) and is able to follow the entire formation history of spatially resolved star clusters, including feedback from individual massive stars. Shortly after the merger we find a population of $\sim 900$ stellar clusters with masses above $10^{2.5}\; \rm{M_\odot}$ and a cluster mass function (CMF), which is well fitted with a power-law with a slope of $\alpha=-1.70\pm0.08$. We describe here in detail the formation of the three most massive clusters (M$_{*} \gtrsim 10^5$ M$_\odot$), which populate the high-mass end of the CMF. The simulated clusters form rapidly on a timescale of $6$-$8$ Myr in converging flows of dense gas. The embedded merger phase has extremely high star formation rate surface densities of $\Sigma_\mathrm{SFR}>10\; \mathrm{M}_\odot\; \mathrm{yr}^{-1}\; \mathrm{kpc}^{-2}$ and thermal gas pressures in excess of $P_{\rm th}\sim10^7 \; \mathrm{k}_{\rm B}\;(\rm K\;\mathrm{cm}^{-3})^{-1}$. The formation process is terminated by rapid gas expulsion driven by the first generation of supernovae, after which the cluster centers relax and both their structure and kinematics become indistinguishable from observed local globular clusters. The simulation presented here provides a general model for the formation of metal-poor globular clusters in chemically unevolved starbursting environments of low-mass dwarf galaxies, which are common at high redshifts.

We study the impact of merging supermassive black holes (SMBHs) on the central regions of massive early-type galaxies (ETGs) using a series of merger simulations with varying mass ratios. The ETG models include realistic stellar and dark matter components and are evolved with the GADGET based regularized tree code KETJU. We show that observed key properties of the nuclear stellar populations of massive ETGs, namely flat stellar density distributions (cores), tangentially biased velocity distributions and kinematically decoupled (counter-)rotation can naturally result from a single process $-$ the scouring by SMBHs. Major mergers with mass ratios of $q>1/3$ produce flat, tangentially biased cores with kinematically distinct components. These features originate from spin reversals during the decay of the SMBH orbits caused by gravitational torques from tidally expelled nuclear material. Minor mergers ($q\lesssim1/3$) on the other hand, form non-rotating cores and the tidal reversal becomes less important. Low-density stellar cores scoured in (multiple) minor mergers are less tangentially biased. This implies that the nuclear stellar properties of massive ETGs can be solely explained by stellar dynamical processes during their final assembly without any need for `feedback' from accreting black holes. We predict a strong correlation between decoupled cores, central anisotropy and merger history: decoupled cores form in binary mergers and we predict them to occur in elliptical galaxies with the strongest central anisotropy. Measurements of the central orbital structure are the key to understanding the number of mergers a given galaxy has experienced.

Given its velocity dispersion, the early-type galaxy NGC 1600 has an unusually massive ($M_\bullet = 1.7 \times 10^{10} M_\odot$) central supermassive black hole (SMBH), surrounded by a large core ($r_\mathrm{b} = 0.7$ kpc) with a tangentially biased stellar distribution. We present high-resolution equal-mass merger simulations including SMBHs to study the formation of such systems. The structural parameters of the progenitor ellipticals were chosen to produce merger remnants resembling NGC 1600. We test initial stellar density slopes of $\rho \propto r^{-1}$ and $\rho \propto r^{-3/2}$ and vary the initial SMBH masses from $8.5 \times 10^8$ to $8.5 \times 10^9$ $M_\odot$. With increasing SMBH mass the merger remnants show a systematic decrease in central surface brightness, an increasing core size, and an increasingly tangentially biased central velocity anisotropy. Two-dimensional kinematic maps reveal decoupled, rotating core regions for the most massive SMBHs. The stellar cores form rapidly as the SMBHs become bound, while the velocity anisotropy develops more slowly after the SMBH binaries become hard. The simulated merger remnants follow distinct relations between the core radius and the sphere-of-influence, and the SMBH mass, similar to observed systems. We find a systematic change in the relations as a function of the progenitor density slope, and present a simple scouring model reproducing this behavior. Finally, we find the best agreement with NGC 1600 using SMBH masses totaling the observed value of $M_\bullet = 1.7 \times 10^{10} M_\odot$. In general, density slopes of $\rho \propto r^{-3/2}$ for the progenitor galaxies are strongly favored for the equal-mass merger scenario.

Ray-tracing is a central tool for constructing mock observations of compact object emission and for comparing physical emission models with observations. We present Arcmancer, a publicly available general ray-tracing and tensor algebra library, written in C++ and providing a Python interface. Arcmancer supports Riemannian and semi-Riemannian spaces of any dimension and metric, and has novel features such as support for multiple simultaneous coordinate charts, embedded geometric shapes, local coordinate systems and automatic parallel propagation. The Arcmancer interface is extensively documented and user-friendly. While these capabilities make the library well suited for a large variety of problems in numerical geometry, the main focus of this paper is in general relativistic polarized radiative transfer. The accuracy of the code is demonstrated in several code tests and in a comparison with GRTRANS, an existing ray-tracing code. We then use the library in several scenarios as a way to showcase the wide applicability of the code. We study a thin variable-geometry accretion disk model, and find that polarization carries information of the inner disk opening angle. Next, we study rotating neutron stars and determine that to obtain polarized light curves at better than $\sim1\%$ level of accuracy, the rotation needs to be taken into account both in the space-time metric as well as in the shape of the star. Finally, we investigate the observational signatures of an accreting black hole lensed by an orbiting black hole. We find that these systems exhibit a characteristic asymmetric twin-peak profile both in flux and polarization properties.

Post-starburst galaxies can be identified via the presence of prominent Hydrogen Balmer absorption lines in their spectra. We present a comprehensive study of the origin of strong Balmer lines in a volume-limited sample of 189 galaxies with $0.01<z<0.05$, $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>9.5$ and projected axis ratio $b/a>0.32$. We explore their structural properties, environments, emission lines and star formation histories, and compare them to control samples of star-forming and quiescent galaxies, and simulated galaxy mergers. Excluding contaminants, in which the strong Balmer lines are most likely caused by dust-star geometry, we find evidence for three different pathways through the post-starburst phase, with most events occurring in intermediate-density environments: (1) a significant disruptive event, such as a gas-rich major merger, causing a starburst and growth of a spheroidal component, followed by quenching of the star formation (70% of post-starburst galaxies at $9.5<\log(\mbox{M}_{\star}/\mbox{M}_{\odot})<10.5$ and 60% at $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>10.5$); (2) at $9.5<\log(\mbox{M}_{\star}/\mbox{M}_{\odot})<10.5$, stochastic star formation in blue-sequence galaxies, causing a weak burst and subsequent return to the blue sequence (30%); (3) at $\log(\mbox{M}_{\star}/\mbox{M}_{\odot})>10.5$, cyclic evolution of quiescent galaxies which gradually move towards the high-mass end of the red sequence through weak starbursts, possibly as a result of a merger with a smaller gas-rich companion (40%). Our analysis suggests that AGN are `on' for $50%$ of the duration of the post-starburst phase, meaning that traditional samples of post-starburst galaxies with strict emission line cuts will be at least $50%$ incomplete due to the exclusion of narrow-line AGN.

We present a high-resolution smoothed particle hydrodynamics simulation of the Antennae galaxies (NGC 4038/4039) and follow the evolution $3$ Gyrs beyond the final coalescence. The simulation includes metallicity dependent cooling, star formation, and both stellar feedback and chemical enrichment. The simulated best-match Antennae reproduces well both the observed morphology and the off-nuclear starburst. We also produce for the first time a simulated two-dimensional metallicity map of the Antennae and find good agreement with the observed metallicity of off-nuclear stellar clusters, however the nuclear metallicities are overproduced by $\sim 0.5$ dex. Using the radiative transfer code SKIRT we produce multi-wavelength observations of both the Antennae and the merger remnant. The $1$ Gyr old remnant is well fitted with a S\'ersic profile of $n=7.07$, and with an $r$-band effective radius of $r_{\mathrm{e}}= 1.6$ kpc and velocity dispersion of $\sigma_{\mathrm{e}}=180$ km$/$s the remnant is located on the fundamental plane of early-type galaxies (ETGs). The initially blue Antennae remnant evolves onto the red sequence after $\sim 2.5$ Gyr of secular evolution. The remnant would be classified as a fast rotator, as the specific angular momentum evolves from $\lambda_{\mathrm{Re}}\approx0.11$ to $\lambda_{\mathrm{Re}}\approx0.14$ during its evolution. The remnant shows ordered rotation and a double peaked maximum in the mean 2D line-of-sight velocity. These kinematical features are relatively common among local ETGs and we specifically identify three local ETGs (NGC 3226, NGC 3379 and NGC 4494) in the ATLAS$^\mathrm{3D}$ sample, whose photometric and kinematic properties most resemble the Antennae remnant.

We present evidence from cosmological hydrodynamical simulations for a co-evolution of the slope of the total (dark and stellar) mass density profile, gamma_tot, and the dark matter fraction within the half-mass radius, f_DM, in early-type galaxies. The relation can be described as gamma_tot = A f_DM + B for all systems at all redshifts. The trend is set by the decreasing importance of gas dissipation towards lower redshifts and for more massive systems. Early-type galaxies are smaller, more concentrated, have lower f_DM and steeper gamma_tot at high redshifts and at lower masses for a given redshift; f_DM and gamma_tot are good indicators for growth by "dry" merging. The values for A and B change distinctively for different feedback models, and this relation can be used as a test for such models. A similar correlation exists between gamma_tot and the stellar mass surface density Sigma_*. A model with weak stellar feedback and feedback from black holes is in best agreement with observations. All simulations, independent of the assumed feedback model, predict steeper gamma_tot and lower f_DM at higher redshifts. While the latter is in agreement with the observed trends, the former is in conflict with lensing observations, which indicate constant or decreasing gamma_tot. This discrepancy is shown to be artificial: the observed trends can be reproduced from the simulations using observational methodology to calculate the total density slopes.

We present KETJU, a new extension of the widely-used smoothed particle hydrodynamics simulation code GADGET-3. The key feature of the code is the inclusion of algorithmically regularized regions around every supermassive black hole (SMBH). This allows for simultaneously following global galactic-scale dynamical and astrophysical processes, while solving the dynamics of SMBHs, SMBH binaries and surrounding stellar systems at sub-parsec scales. The KETJU code includes Post-Newtonian terms in the equations of motions of the SMBHs which enables a new SMBH merger criterion based on the gravitational wave coalescence timescale pushing the merger separation of SMBHs down to $\sim 0.005$ pc. We test the performance of our code by comparison to NBODY7 and rVINE. We set up dynamically stable multi-component merger progenitor galaxies to study the SMBH binary evolution during galaxy mergers. In our simulation sample the SMBH binaries do not suffer from the final-parsec problem, which we attribute to the non-spherical shape of the merger remnants. For bulge-only models, the hardening rate decreases with increasing resolution, whereas for models which in addition include massive dark matter halos the SMBH binary hardening rate becomes practically independent of the mass resolution of the stellar bulge. The SMBHs coalesce on average 200 Myr after the formation of the SMBH binary. However, small differences in the initial SMBH binary eccentricities can result in large differences in the SMBH coalescence times. Finally, we discuss the future prospects of KETJU, which allows for a straightforward inclusion of gas physics in the simulations.

The Direct Collapse Black Hole (DCBH) scenario provides a solution for forming the massive black holes powering bright quasars observed in the early Universe. A prerequisite for forming a DCBH is that the formation of (much less massive) Population III stars be avoided - this can be achieved by destroying H$_2$ via Lyman-Werner (LW) radiation (E$_{\rm{LW}}$ = 12.6 eV). We find that two conditions must be met in the proto-galaxy that will host the DCBH. First, prior star formation must be delayed; this can be achieved with a background LW flux of J$_{\rm BG} \gtrsim 100\ J_{21}$. Second, an intense burst of LW radiation from a neighbouring star-bursting proto-galaxy is required, just before the gas cloud undergoes gravitational collapse, to finally suppress star formation completely. We show here for the first time using high-resolution hydrodynamical simulations, including full radiative transfer, that this low-level background, combined with tight synchronisation and irradiation of a secondary proto-galaxy by a primary proto-galaxy, inevitably moves the secondary proto-galaxy onto the isothermal atomic cooling track, without the deleterious effects of either photo-evaporating the gas or polluting it by heavy elements. These, atomically cooled, massive proto-galaxies are expected to ultimately form a DCBH of mass $10^4 - 10^5 M_{\odot}$.

We present new SINFONI near-infrared integral field unit (IFU) spectroscopy and SALT optical long-slit spectroscopy characterising the history of a nearby merging luminous infrared galaxy, dubbed the Bird (IRAS19115-2114). The NIR line-ratio maps of the IFU data-cubes and stellar population fitting of the SALT spectra now allow dating of the star formation (SF) over the triple system uncovered from our previous adaptive optics data. The distinct components separate very clearly in a line-ratio diagnostic diagram. An off-nuclear pure starburst dominates the current SF of the Bird with 60-70% of the total, with a 4-7 Myr age, and signs of a fairly constant long-term star formation of the underlying stellar population. The most massive nucleus, in contrast, is quenched with a starburst age of >40 Myr and shows hints of budding AGN activity. The secondary massive nucleus is at an intermediate stage. The two major components have a population of older stars, consistent with a starburst triggered 1 Gyr ago in a first encounter. The simplest explanation of the history is that of a triple merger, where the strongly star forming component has joined later. We detect multiple gas flows in different phases. The Bird offers an opportunity to witness multiple stages of galaxy evolution in the same system; triggering as well as quenching of SF, and the early appearance of AGN activity. It also serves as a cautionary note on interpretations of observations with lower spatial resolution and/or without infrared data. At high-redshift the system would look like a clumpy starburst with crucial pieces of its puzzle hidden, in danger of misinterpretations.

We discuss the problem of polarized radiative transfer in general relativity. We present a set of equations suitable for solving the problem numerically for the case of an arbitrary space-time metric, and show numerical solutions to example problems. The solutions are computed with a new ray-tracing code, Arcmancer, developed by the authors.

The predicted abundance and properties of the low-mass substructures embedded inside larger dark matter haloes differ sharply among alternative dark matter models. Too small to host galaxies themselves, these subhaloes may still be detected via gravitational lensing, or via perturbations of the Milky Way's globular cluster streams and its stellar disk. Here we use the Apostle cosmological simulations to predict the abundance and the spatial and velocity distributions of subhaloes in the range 10^6.5-10^8.5 solar masses inside haloes of mass ~ 10^12 solar masses in LCDM. Although these subhaloes are themselves devoid of baryons, we find that baryonic effects are important. Compared to corresponding dark matter only simulations, the loss of baryons from subhaloes and stronger tidal disruption due to the presence of baryons near the centre of the main halo, reduce the number of subhaloes by ~ 1/4 to 1/2, independently of subhalo mass, but increasingly towards the host halo centre. We also find that subhaloes have non-Maxwellian orbital velocity distributions, with centrally rising velocity anisotropy and positive velocity bias which reduces the number of low-velocity subhaloes, particularly near the halo centre. We parameterise the predicted population of subhaloes in terms of mass, galactocentric distance, and velocities. We discuss implications of our results for the prospects of detecting dark matter substructures and for possible inferences about the nature of dark matter.

A nearby source of Lyman-Werner (LW) photons is thought to be a central component in dissociating H$_2$ and allowing for the formation of a direct collapse black hole seed. Nearby sources are also expected to produce copious amounts of hydrogen ionising photons and X-ray photons. We study here the feedback effects of the X-ray photons by including a spectrum due to high-mass X-ray binaries on top of a galaxy with a stellar spectrum. We explicitly trace photon packages emerging from the nearby source and track the radiative and chemical effects of the multi-frequency source $(E_{\rm photon} = \rm{0.76\ eV \rightarrow 7500\ eV}$). We find that X-rays have a strongly negative feedback effect, compared to a stellar only source, when the radiative source is placed at a separation greater than $\gtrsim 1 \ \rm kpc$. The X-rays heat the low and medium density gas in the envelope surrounding the collapsing halo suppressing the mass inflow. The result is a smaller enclosed mass compared to the stellar only case. However, for separations of $\lesssim 1 \ \rm kpc$, the feedback effects of the X-rays becomes somewhat neutral. The enhanced LW intensity at close separations dissociates more H$_2$ and this gas is heated due to stellar photons alone, the addition of X-rays is then not significant. This distance dependence of X-ray feedback suggests that a Goldilocks zone exists close to a forming galaxy where X-ray photons have a much smaller negative feedback effect and ideal conditions exist for creating massive black hole seeds.

The photo-dissociation of H$_2$ by a nearby anisotropic source of radiation is seen as a critical component in creating an environment in which a direct collapse black hole may form. Employing radiative transfer we model the effect of multi-frequency (0.76 eV - 60 eV) radiation on a collapsing halo at high redshift. We vary both the shape of the spectrum which emits the radiation and the distance to the emitting galaxy. We use blackbody spectra with temperatures of $\rm{T = 10^4\ K}$ and $\rm{T = 10^5\ K}$ and a realistic stellar spectrum. We find that an optimal zone exists between 1 kpc and 4 kpc from the emitting galaxy. If the halo resides too close to the emitting galaxy the photo-ionising radiation creates a large HII region which effectively disrupts the collapsing halo, too far from the source and the radiation flux drops below the level of the expected background and the H$_2$ fraction remains too high. When the emitting galaxy is initially placed between 1 kpc and 2 kpc from the collapsing halo, with a spectral shape consistent with a star-forming high redshift galaxy, then a large central core forms. The mass of the central core is between 5000 and 10000 $\rm{M_{\odot}}$ at a temperature of approximately 1000 K. This core is however surrounded by a reservoir of hotter gas at approximately 8000 K which leads to mass inflow rates of the order of $\sim 0.1$ $\rm{M_{\odot}}$ yr$^{-1}$. This environment has the potential to form a massive primordial star which can then lead to the formation of a direct collapse black hole.

We present a new morphological indicator designed for automated recognition of galaxies with faint asymmetric tidal features suggestive of an ongoing or past merger. We use this new indicator, together with preexisting diagnostics of galaxy structure to study the role of galaxy mergers in inducing (post-)starburst spectral signatures in local galaxies, and investigate whether (post-)starburst galaxies play a role in the build up of the `red sequence'. Our morphological and structural analysis of an evolutionary sample of 335 (post-)starburst galaxies in the SDSS DR7 with starburst ages 0<tSB<0.6 Gyr, shows that 45% of galaxies with young starbursts (tSB<0.1 Gyr) show signatures of an ongoing or past merger. This fraction declines with starburst age, and we find a good agreement between automated and visual classifications. The majority of the oldest (post-)starburst galaxies in our sample (tSB~0.6Gyr) have structural properties characteristic of early-type disks and are not as highly concentrated as the fully quenched galaxies commonly found on the `red sequence' in the present day Universe. This suggests that, if (post-)starburst galaxies are a transition phase between active star-formation and quiescence, they do not attain the structure of presently quenched galaxies within the first 0.6 Gyr after the starburst.

We examine the impact of dark matter particle resolution on the formation of a baryonic core in high resolution adaptive mesh refinement simulations. We test the effect that both particle smoothing and particle splitting have on the hydrodynamic properties of a collapsing halo at high redshift (z > 20). Furthermore, we vary the background field intensity, with energy below the Lyman limit (< 13.6 eV), as may be relevant for the case of metal-free star formation and super-massive black hole seed formation. We find that using particle splitting methods greatly increases our particle resolution without introducing any numerical noise and allows us to achieve converged results over a wide range of external background fields. Additionally, we find that for lower values of the background field a lower dark matter particle mass is required. We define the radius of the core as the point at which the enclosed baryonic mass dominates over the enclosed dark matter mass. For our simulations this results in $\rm{R_{core} \sim 5\ pc}$. We find that in order to produce converged results which are not affected by dark matter particles requires that the relationship ${M_{\rm{core}} / M_{\rm{DM}}} > 100.0$ be satisfied, where ${M_{\rm{core}}}$ is the enclosed baryon mass within the core and $M_{\rm{DM}}$ is the minimum dark matter particle mass. This ratio should provide a very useful starting point for conducting convergence tests before any production run simulations. We find that dark matter particle smoothing is a useful adjunct to already highly resolved simulations.

We study the evolution of the gaseous components in massive simulated galaxies and show that their early formation is fuelled by cold, low entropy gas streams. At lower redshifts of z<3 the simulated galaxies are massive enough to support stable virial shocks resulting in a transition from cold to hot gas accretion. The gas accretion history of early-type galaxies is directly linked to the formation of their stellar component in the two phased formation scenario, in which the central parts of the galaxy assemble rapidly through in situ star formation and the later assembly is dominated primarily by minor stellar mergers.

The direct collapse model of supermassive black hole seed formation provides an attractive solution to the origin of the quasars now routinely observed at $z \gtrsim 6$. We use the adaptive mesh refinement code Enzo to simulate the collapse of gas at high redshift, including a nine species chemical model of H, He, and H$_2$. The direct collapse model requires that the gas cools predominantly via atomic hydrogen. To this end we simulate the effect of an anisotropic radiation source on the collapse of a halo at high redshift. The radiation source is placed at a distance of 3 kpc (physical) from the collapsing object. The source is set to emit monochromatically in the center of the Lyman-Werner (LW) band only at $12.8 \ \rm{eV}$. The LW radiation emitted from the high redshift source is followed self-consistently using ray tracing techniques. We find that, due to self-shielding, a small amount of H$_2$ is able to form at the very center of the collapsing halo even under very strong LW radiation. Furthermore, we find that a radiation source, emitting $> 10^{54}\ (\sim10^3\ \rm{J_{21}})$ photons per second is required to cause the collapse of a clump of $\rm{M \sim 10^5}$ M$_{\odot}$. The resulting accretion rate onto the collapsing object is $\sim 0.25$ M$_{\odot}$ $\rm{yr^{-1}}$. Our results display significant differences, compared to the isotropic radiation field case, in terms of H$_2$ fraction at an equivalent radius. These differences will significantly effect the dynamics of the collapse. With the inclusion of a strong anisotropic radiation source, the final mass of the collapsing object is found to be $\rm{M \sim 10^5}$ M$_{\odot}$. This is consistent with predictions for the formation of a supermassive star or quasi-star leading to a supermassive black hole.

We present a catalogue and images of visually detected features, such as asymmetries, extensions, warps, shells, tidal tails, polar rings, and obvious signs of mergers or interactions, in the faint outer regions (at and outside of R_25) of nearby galaxies. This catalogue can be used in future quantitative studies that examine galaxy evolution due to internal and external factors. We are able to reliably detect outer region features down to a brightness level of 0.03 MJy/sr per pixel at 3.6 microns in the Spitzer Survey of Stellar Structure in Galaxies (S4G). We also tabulate companion galaxies. We find asymmetries in the outer isophotes in 22+/-1 per cent of the sample. The asymmetry fraction does not correlate with galaxy classification as an interacting galaxy or merger remnant, or with the presence of companions. We also compare the detected features to similar features in galaxies taken from cosmological zoom re-simulations. The simulated images have a higher fraction (33 per cent) of outer disc asymmetries, which may be due to selection effects and an uncertain star formation threshold in the models. The asymmetries may have either an internal (e.g., lopsidedness due to dark halo asymmetry) or external origin.

We use a large suite of hydrodynamical simulations of binary galaxy mergers to construct and calibrate a physical prescription for computing the effective radii and velocity dispersions of spheroids. We implement this prescription within a semi-analytic model embedded in merger trees extracted from the Bolshoi Lambda-CDM N-body simulation, accounting for spheroid growth via major and minor mergers as well as disk instabilities. We find that without disk instabilities, our model does not predict sufficient numbers of intermediate mass early-type galaxies in the local universe. Spheroids also form earlier in models with spheroid growth via disk instabilities. Our model correctly predicts the normalization, slope, and scatter of the low-redshift size-mass and Fundamental Plane relations for early type galaxies. It predicts a degree of curvature in the Faber-Jackson relation that is not seen in local observations, but this could be alleviated if higher mass spheroids have more bottom-heavy initial mass functions. The model also correctly predicts the observed strong evolution of the size-mass relation for spheroids out to higher redshifts, as well as the slower evolution in the normalization of the Faber-Jackson relation. We emphasize that these are genuine predictions of the model since it was tuned to match hydrodynamical simulations and not these observations.

We present optical integral field spectroscopy (IFS) observations of the Mice, a major merger between two massive (>10^11Msol) gas-rich spirals NGC4676A and B, observed between first passage and final coalescence. The spectra provide stellar and gas kinematics, ionised gas properties and stellar population diagnostics, over the full optical extent of both galaxies. The Mice provide a perfect case study highlighting the importance of IFS data for improving our understanding of local galaxies. The impact of first passage on the kinematics of the stars and gas has been significant, with strong bars likely induced in both galaxies. The barred spiral NGC4676B exhibits a strong twist in both its stellar and ionised gas disk. On the other hand, the impact of the merger on the stellar populations has been minimal thus far: star formation induced by the recent close passage has not contributed significantly to the global star formation rate or stellar mass of the galaxies. Both galaxies show bicones of high ionisation gas extending along their minor axes. In NGC4676A the high gas velocity dispersion and Seyfert-like line ratios at large scaleheight indicate a powerful outflow. Fast shocks extend to ~6.6kpc above the disk plane. The measured ram pressure and mass outflow rate (~8-20Msol/yr) are similar to superwinds from local ULIRGs, although NGC4676A has only a moderate infrared luminosity of 3x10^10Lsol. Energy beyond that provided by the mechanical energy of the starburst appears to be required to drive the outflow. We compare the observations to mock kinematic and stellar population maps from a merger simulation. The models show little enhancement in star formation during and following first passage, in agreement with the observations. We highlight areas where IFS data could help further constrain the models.

We study the effect of AGN mechanical and radiation feedback on the formation of bulge dominated galaxies via mergers of disc galaxies. The merging galaxies have mass-ratios of 1:1 to 6:1 and include pre-existing hot gaseous halos to properly account for the global impact of AGN feedback. Using smoothed particle hydrodynamics simulation code (GADGET-3) we compare three models with different AGN feedback models: (1) no black hole and no AGN feedback; (2) thermal AGN feedback; and (3) mechanical and radiative AGN feedback. The last model is motivated by observations of broad line quasars which show winds with initial velocities of $v_w \ge$ 10,000 km/s and also heating associated with the central AGN X-ray radiation. The primary changes in gas properties due to mechanical AGN feedback are lower thermal X-ray luminosity from the final galaxy - in better agreement with observations - and galactic outflows with higher velocity $\sim 1000$ km/s similar to recent direct observations of nearby merger remnants. The kinetic energy of the outflowing gas is a factor of $\sim$ 20 higher than in the thermal feedback case. All merger remnants with momentum-based AGN feedback with $v_w \sim 10,000$ km/s and $\epsilon_w=2 \times 10^{-3}$, independent of their progenitor mass-ratios, reproduce the observed relations between stellar velocity dispersion and black hole mass ($M_{\rm bh} - \sigma$) as well as X-ray luminosity ($L_X - \sigma$) with $10^{37.5}$ erg/s $\lesssim L_X (0.3-8~{\rm keV}) \lesssim 10^{39.2}$ erg/s for velocity dispersions in the range of 120 km/s $\lesssim \sigma \lesssim$ 190 km/s. In addition, the mechanical feedback produces a much greater AGN variability. We also show that gas is more rapidly and impulsively stripped from the galactic centres driving a moderate increase in galaxy size and decrease in central density with the mechanical AGN feedback model.

We have performed high-resolution numerical simulations with the hydrodynamical AMR code Enzo to investigate the formation of massive seed black holes in a sample of six dark matter haloes above the atomic cooling threshold. The aim of this study is to illustrate the effects of varying the maximum refinement level on the final object formed. The virial temperatures of the simulated haloes range from $\rm{T} \sim 10000\ \rm{K} - 16000\ \rm{K}$ and they have virial masses in the range $\rm{M} \sim 2 \times 10^7 \rm{M_{\odot}}$ to $\rm{M} \sim 7 \times 10^7 \rm{M_{\odot}}$ at $z \sim 15$. The outcome of our six fiducial simulations is both generic and robust. A rotationally supported, marginally gravitationally stable, disk forms with an exponential profile. The mass and scale length of this disk depends strongly on the maximum refinement level used. Varying the maximum refinement level by factors between 1 / 64 to 256 times the fiducial level illustrates the care that must be taken in interpreting the results. The lower resolution simulations show tentative evidence that the gas may become rotationally supported out to 20 pc while the highest resolution simulations show only weak evidence of rotational support due to the shorter dynamical times for which the simulation runs. The higher resolution simulations do, however, point to fragmentation at small scales of the order of $\sim 100$ AU. In the highest resolution simulations a central object of a few times $10^2\ \rm{M_{\odot}}$ forms with multiple strongly bound, Jeans unstable, clumps of $\sim 10\ \rm{M_{\odot}}$ and radii of 10 - 20 AU suggesting the formation of dense star clusters in these haloes.

The stellar velocity dispersion as a function of the galactocentric radius of an early-type galaxy can generally be well approximated by a power law $\sigma \propto r^{\beta}$. However, some observed dispersion profiles show a deviation from this fit at intermediate radii, usually between one and three $R_{\mathrm{eff}}$, where the velocity dispersion remains constant with radius, showing a bump-like behavior, which we term the "$\sigma$-"bump. To understand the origin of this $\sigma$-bump, we study a set of simulated early-type galaxies formed in major mergers. We find the $\sigma$-bump in all of our simulated early-type galaxies, with the size and position of the bump slightly varying from galaxy to galaxy, leading to the assumption that the bump is a characteristic of the major merger formation scenario. The feature can be seen both in the intrinsic and projected stellar velocity dispersions. In contrast to shells that form during the merger event but evolve with time and finally disappear, the $\sigma$-bump stays nearly constant with radius and is a permanent feature that is preserved until the end of the simulation. The $\sigma$-bump is not seen in the dark matter and gas components and we therefore conclude that it is a purely stellar feature of merger remnants.

We present a detailed two-dimensional stellar dynamical analysis of a sample of 44 cosmological hydrodynamical simulations of individual central galaxies and their satellites. Kinematic maps of the stellar line-of-sight velocity, velocity dispersion, and higher-order Gauss-Hermite moments $h_3$ and $h_4$ are constructed for each central galaxy and for the most massive satellites. The amount of rotation is quantified using the $\lambda_{\mathrm{R}}$-parameter. The velocity, velocity dispersion, $h_3$, and $h_4$ fields of the simulated galaxies show a diversity similar to observed kinematic maps of early-type galaxies in the ATLAS$^{\rm{3D}}$ survey. This includes fast (regular), slow, and misaligned rotation, hot spheroids with embedded cold disk components as well as galaxies with counter-rotating cores or central depressions in the velocity dispersion. We link the present day kinematic properties to the individual cosmological formation histories of the galaxies. In general, major galaxy mergers have a significant influence on the rotation properties resulting in both a spin-down as well as a spin-up of the merger remnant. Lower mass galaxies with significant in-situ formation of stars, or with additional gas-rich major mergers - resulting in a spin-up - in their formation history, form elongated fast rotators with a clear anti-correlation of $h_3$ and $v/\sigma$. An additional formation path for fast rotators includes gas poor major mergers leading to a spin-up of the remnants. This formation path does not result in anti-correlated $h_3$ and $v/\sigma$. The galaxies most consistent with the rare class of non-rotating round early-type galaxies grow by gas-poor minor mergers alone. In general, more massive galaxies have less in-situ star formation since $z \sim 2$, rotate slower and have older stellar populations. (shortened)

We present a set of hydro-dynamical numerical simulations of the Antennae galaxies in order to understand the origin of the central overlap starburst. Our dynamical model provides a good match to the observed nuclear and overlap star formation, especially when using a range of rather inefficient stellar feedback efficiencies (0.01 < q_EoS < 0.1). In this case a simple conversion of local star formation to molecular hydrogen surface density motivated by observations accounts well for the observed distribution of CO. Using radiative transfer post-processing we model synthetic far-infrared spectral energy distributions (SEDs) and two-dimensional emission maps for direct comparison with Herschel-PACS observations. For a gas-to-dust ratio of 62:1 and the best matching range of stellar feedback efficiencies the synthetic far-infrared SEDs of the central star forming region peak at values of ~65 - 81 Jy at 99 - 116 um, similar to a three-component modified black body fit to infrared observations. Also the spatial distribution of the far-infrared emission at 70 um, 100 um, and 160 um compares well with the observations: >50% (> 35%) of the emission in each band is concentrated in the overlap region while only < 30% (< 15%) is distributed to the combined emission from the two galactic nuclei in the simulations (observations). As a proof of principle we show that parameter variations in the feedback model result in unambiguous changes both in the global and in the spatially resolved observable far-infrared properties of Antennae galaxy models. Our results strengthen the importance of direct, spatially resolved comparative studies of matched galaxy merger simulations as a valuable tool to constrain the fundamental star formation and feedback physics.

Dynamical modeling and strong lensing data indicate that the total density profiles of early-type galaxies are close to isothermal, i.e., rho_tot ~ r^gamma with gamma approx -2. To understand the origin of this universal slope we study a set of simulated spheroids formed in isolated binary mergers as well as the formation within the cosmological framework. The total stellar plus dark matter density profiles can always be described by a power law with an index of gamma approx -2.1 with a tendency toward steeper slopes for more compact, lower-mass ellipticals. In the binary mergers the amount of gas involved in the merger determines the precise steepness of the slope. This agrees with results from the cosmological simulations where ellipticals with steeper slopes have a higher fraction of stars formed in situ. Each gas-poor merger event evolves the slope toward gamma ~ -2, once this slope is reached further merger events do not change it anymore. All our ellipticals have flat intrinsic combined stellar and dark matter velocity dispersion profiles. We conclude that flat velocity dispersion profiles and total density distributions with a slope of gamma ~ -2 for the combined system of stars and dark matter act as a natural attractor. The variety of complex formation histories as present in cosmological simulations, including major as well as minor merger events, is essential to generate the full range of observed density slopes seen for present-day elliptical galaxies.

We demonstrate that massive simulated galaxies assemble in two phases, with the initial growth dominated by compact in situ star formation, whereas the late growth is dominated by accretion of old stars formed in subunits outside the main galaxy. We also show that 1) gravitational feedback strongly suppresses late star formation in massive galaxies contributing to the observed galaxy colour bimodality that 2) the observed galaxy downsizing can be explained naturally in the two-phased model and finally that 3) the details of the assembly histories of massive galaxies are directly connected to their observed kinematic properties.

We study the growth of black holes (BHs) in galaxies using three-dimensional smoothed particle hydrodynamic simulations with new implementations of the momentum mechanical feedback, and restriction of accreted elements to those that are gravitationally bound to the BH. We also include the feedback from the X-ray radiation emitted by the BH, which heats the surrounding gas in the host galaxies, and adds radial momentum to the fluid. We perform simulations of isolated galaxies and merging galaxies and test various feedback models with the new treatment of the Bondi radius criterion. We find that overall the BH growth is similar to what has been obtained by earlier workers using the Springel, Di Matteo, & Hernquist algorithms. However, the outflowing wind velocities and mechanical energy emitted by winds are considerably higher (v_w ~ 1000-3000 km/s) compared to the standard thermal feedback model (v_w ~ 50-100 km/s). While the thermal feedback model emits only 0.1 % of BH released energy in winds, the momentum feedback model emits more than 30 % of the total energy released by the BH in winds. In the momentum feedback model, the degree of fluctuation in both radiant and wind output is considerably larger than in the standard treatments. We check that the new model of the BH mass accretion agrees with analytic results for the standard Bondi problem.

We present a sample of nine high resolution cosmological simulations in the mass range of M_vir=7x10^11-4x10^12 M_sun starting from LambdaCDM initial conditions. Our simulations include primordial radiative cooling, photoionization, star formation, supernova II feedback, but exclude supernova driven winds and AGN feedback. The simulated galaxies assemble in two phases, with the initial growth dominated by compact (r<r_eff) in situ star formation fueled by cold, low entropy gas streams resulting in a very similar mean assembly redshift of z_{f,ins}~2.5 for the in situ stellar component in all galaxies. The late growth is dominated by accretion of old stars formed in subunits outside the main galaxy (r>r_eff) resulting in an assembly redshift of z_{f,acc}~0.5-1.5 with much larger scatter. We find a positive correlation between the fraction of accreted stars and the final mass of our galaxies. We show that gravitational feedback strongly suppresses late star formation in massive galaxies contributing to the observed galaxy color bimodality. The accretion of stellar material is also responsible for the observed size growth of early-type galaxies. In addition, we find that the dark matter fractions within the stellar half-mass radii continuously increase towards lower redshift from about f_DM~0.05 at z~3 to f_DM~0.1-0.3 at z=0. Furthermore, the logarithmic slope of the total density profile is nearly isothermal at the present-day (gamma'~1.9-2.2). Finally, the input of gravitational heating lowers the central dark matter densities in the galaxies, with the effect being smaller compared to simulations without supernova feedback.

We analyze 40 cosmological re-simulations of individual massive galaxies with present-day stellar masses of $M_{*} > 6.3 \times 10^{10} M_{\odot}$ in order to investigate the physical origin of the observed strong increase in galaxy sizes and the decrease of the stellar velocity dispersions since redshift $z \approx 2$. At present 25 out of 40 galaxies are quiescent with structural parameters (sizes and velocity dispersions) in agreement with local early type galaxies. At z=2 all simulated galaxies with $M_* \gtrsim 10^{11}M_{\odot}$ (11 out of 40) at z=2 are compact with projected half-mass radii of $\approx$ 0.77 ($\pm$0.24) kpc and line-of-sight velocity dispersions within the projected half-mass radius of $\approx$ 262 ($\pm$28) kms$^{-1}$ (3 out of 11 are already quiescent). Similar to observed compact early-type galaxies at high redshift the simulated galaxies are clearly offset from the local mass-size and mass-velocity dispersion relations. Towards redshift zero the sizes increase by a factor of $\sim 5-6$, following $R_{1/2} \propto (1+z)^{\alpha}$ with $\alpha = -1.44$ for quiescent galaxies ($\alpha = -1.12$ for all galaxies). The velocity dispersions drop by about one-third since $z \approx 2$, following $\sigma_{1/2} \propto (1+z)^{\beta}$ with $\beta = 0.44$ for the quiescent galaxies ($\beta = 0.37$ for all galaxies). The simulated size and dispersion evolution is in good agreement with observations and results from the subsequent accretion and merging of stellar systems at $z\lesssim 2$ which is a natural consequence of the hierarchical structure formation. A significant number of the simulated massive galaxies (7 out of 40) experience no merger more massive than 1:4 (usually considered as major mergers). On average, the dominant accretion mode is stellar minor mergers with a mass-weighted mass-ratio of 1:5. (abridged)

Many observed massive star-forming z\approx2 galaxies are large disks that exhibit irregular morphologies, with \sim1kpc, \sim10^(8-10)Msun clumps. We present the largest sample to date of high-resolution cosmological SPH simulations that zoom-in on the formation of individual M*\sim10^(10.5)Msun galaxies in \sim10^(12)Msun halos at z\approx2. Our code includes strong stellar feedback parameterized as momentum-driven galactic winds. This model reproduces many characteristic features of this observed class of galaxies, such as their clumpy morphologies, smooth and monotonic velocity gradients, high gas fractions (f_g\sim50%) and high specific star-formation rates (\gtrsim1Gyr^(-1)). In accord with recent models, giant clumps (Mclump\sim(5x10^8-10^9)Msun) form in-situ via gravitational instabilities. However, the galactic winds are critical for their subsequent evolution. The giant clumps we obtain are short-lived and are disrupted by wind-driven mass loss. They do not virialise or migrate to the galaxy centers as suggested in recent work neglecting strong winds. By phenomenologically implementing the winds that are observed from high-redshift galaxies and in particular from individual clumps, our simulations reproduce well new observational constraints on clump kinematics and clump ages. In particular, the observation that older clumps appear closer to their galaxy centers is reproduced in our simulations, as a result of inside-out formation of the disks rather than inward clump migration.

Hydrodynamic simulations of galaxies with active galactic nuclei (AGN) have typically employed feedback that is purely local: i.e., an injection of energy to the immediate neighborhood of the black hole. We perform GADGET-2 simulations of massive elliptical galaxies with an additional feedback component: an observationally calibrated X-ray radiation field which emanates from the black hole and heats gas out to large radii from the galaxy center. We find that including the heating and radiation pressure associated with this X-ray flux in our simulations enhances the effects which are commonly reported from AGN feedback. This new feedback model is twice as effective as traditional feedback at suppressing star formation, produces 3 times less star formation in the last 6 Gyr, and modestly lowers the final BH mass (30%). It is also significantly more effective than an X-ray background in reducing the number of satellite galaxies.

We present high resolution simulations of a multiple merger of three disk galaxies including the evolution of magnetic fields performed with the N-body/SPH code Gadget. For the first time, we embed the galaxies in a magnetized, low-density medium, thus modeling an ambient IGM. The simulations include radiative cooling and a model for star formation and supernova feedback. The progenitor disks have initial magnetic seed fields in the range of 10e-9 to 10e-6 G and the IGM has initial fields of 10e-12 to 10e-9 G. The simulations are compared to a run excluding magnetic fields. We show that the propagation of interaction-driven shocks depends significantly on the initial magnetic field strength. The shocks propagate faster in simulations with stronger initial field, suggesting that the shocks are supported by magnetic pressure. The Mach numbers of the shocks range from approximately M=1.5 for the non-magnetized case up to M=6 for the highest initial magnetization, resulting in higher temperatures of the shock heated IGM gas. The magnetic field in the system saturates rapidly after the mergers at ~ 10e-6 G within the galaxies and ~ 10e-8 G in the IGM independent of the initial value. These field strengths agree with observed values and correspond to the equipartition value of the magnetic pressure with the turbulent pressure in the system. We also present synthetic radio and polarization maps for different phases of the evolution showing that shocks driven by the interaction produce a high amount of polarized emission. These idealized simulations indicate that magnetic fields play an important role for the hydrodynamics of the IGM during galactic interactions. We also show that even weak seed fields are efficiently strengthened during multiple galactic mergers. This interaction driven amplification might have been a key process for the magnetization of the Universe.

Cosmological simulations of galaxy formation appear to show a two-phase character with a rapid early phase at z>2 during which in-situ stars are formed within the galaxy from infalling cold gas followed by an extended phase since z<3 during which ex-situ stars are primarily accreted. In the latter phase massive systems grow considerably in mass and radius by accretion of smaller satellite stellar systems formed at quite early times (z>3) outside of the virial radius of the forming central galaxy. These tentative conclusions are obtained from high resolution re-simulations of 39 individual galaxies in a full cosmological context with present-day virial halo masses ranging from 7e11 M_sun h^-1 < M_vir < 2.7e13 M_sun h^-1 and central galaxy masses between 4.5e10 M_sun h^-1 < M_* < 3.6e11 M_sun h^-1. The simulations include the effects of a uniform UV background, radiative cooling, star formation and energetic feedback from SNII. The importance of stellar accretion increases with galaxy mass and towards lower redshift. In our simulations lower mass galaxies ($M_* < 9e10 M_sun h^-1) accrete about 60 per cent of their present-day stellar mass. High mass galaxy ($M_* > 1.7e11 M_sun h^-1) assembly is dominated by accretion and merging with about 80 per cent of the stars added by the present-day. In general the simulated galaxies approximately double their mass since z=1. For massive systems this mass growth is not accompanied by significant star formation. The majority of the in-situ created stars is formed at z>2, primarily out of cold gas flows. We recover the observational result of archaeological downsizing, where the most massive galaxies harbor the oldest stars. We find that this is not in contradiction with hierarchical structure formation. Most stars in the massive galaxies are formed early on in smaller structures, the galaxies themselves are assembled late.

We perform numerical simulations of clouds in the Galactic Centre (GC) engulfing the nuclear super-massive black hole and show that this mechanism leads to the formation of gaseous accretion discs with properties that are similar to the expected gaseous progenitor discs that fragmented into the observed stellar disc in the GC. As soon as the cloud hits the black hole, gas with opposite angular momentum relative to the black hole collides downstream. This process leads to redistribution of angular momentum and dissipation of kinetic energy, resulting in a compact gaseous accretion disc. A parameter study using thirteen high resolution simulations of homogeneous clouds falling onto the black hole and engulfing it in parts demonstrates that this mechanism is able to produce gaseous accretion discs that could potentially be the progenitor of the observed stellar disc in the GC. A comparison of simulations with different equations of state (adiabatic, isothermal and full cooling) demonstrates the importance of including a detailed thermodynamical description. However the simple isothermal approach already yields good results on the radial mass transfer and accretion rates, as well as disc eccentricities and sizes. We find that the cloud impact parameter strongly influences the accretion rate whereas the impact velocity has a small affect on the accretion rate.

Even though the dark-matter power spectrum in the absence of biasing predicts a number density of halos n(M) ~ M^-2 (i.e. a Schechter alpha value of -2) at the low-mass end (M < 10^10 M_solar), hydrodynamic simulations have typically produced values for stellar systems in good agreement with the observed value alpha ~ -1. We explain this with a simple physical argument and show that an efficient external gas-heating mechanism (such as the UV background included in all hydro codes) will produce a critical halo mass below which halos cannot retain their gas and form stars. We test this conclusion with GADGET-2-based simulations using various UV backgrounds, and for the first time we also investigate the effect of an X-ray background. We show that at the present epoch alpha is depends primarily on the mean gas temperature at the star-formation epoch for low-mass systems (z <~ 3): with no background we find alpha ~ -1.5, with UV only alpha ~ -1.0, and with UV and X-rays alpha ~ -0.75. We find the critical final halo mass for star formation to be ~4x10^8 M_solar with a UV background and ~7x10^8 M_solar with UV and X-rays.

We present spatially resolved integral field spectroscopic K-band data at a resolution of 0.13" (60pc) and interferometric CO(2-1) line observations of the prototypical merging system NGC6240. Despite the clear rotational signature, the stellar kinematics in the two nuclei are dominated by dispersion. We use Jeans modelling to derive the masses and the mass-to-light ratios of the nuclei. Combining the luminosities with the spatially resolved Br-gamma equivalent width shows that only 1/3 of the K-band continuum from the nuclei is associated with the most recent star forming episode; and that less than 30% of the system's bolometric luminosity and only 9% of its stellar mass is due to this starburst. The star formation properties, calculated from typical merger star formation histories, demonstrate the impact of different assumptions about the star formation history. The properties of the nuclei, and the existence of a prominent old stellar population, indicate that the nuclei are remnants of the progenitor galaxies' bulges.

We use a series of idealized, numerical SPH simulations to study the formation and evolution of galactic, gas-rich disks forming from gas infall within dark matter halos. The temperature and density structure of the gas is varied in order to differentiate between (i) simultaneous gas infall at a large range of radii and (ii) the inside-out build-up of a disk. In all cases, the disks go through phases of ring formation, gravitational instability and break-up into massive clumps. Ring formation can be enhanced by a focal point effect. The position of the ring is determined by the angular momentum distribution of the material it forms from. We study the ring and clump morphologies, the characteristic properties of the resulting velocity dispersion field and the effect of star formation. In the early phases, gas accretion leads to a high vertical velocity dispersion. We find that the disk fragmentation by gravitational instability and the subsequent clump-clump interactions drive high velocity dispersions mainly in the plane of the disk while at the same time the vertical velocity dispersion dissipates. The result is a strong variation of the line-of-sight velocity dispersion with inclination angle. For a face-on view, clumps appear as minima in the (vertical) dispersion, whereas for a more edge-on view, they tend to correspond to maxima. There exists observational evidence of a systematic variation of the velocity dispersion with inclination angle in high-redshift disks, which could be partly explained by our simulation results. Additional energetic sources to drive velocity dispersion that are not included in our models are also expected to contribute to the observational results.

We present self-consistent high-resolution simulations of NGC4038/4039 (the "Antennae galaxies") including star formation, supernova feedback and magnetic fields performed with the N-body/SPH code Gadget, in which magnetohydrodynamics are followed with the SPH method. We vary the initial magnetic field in the progenitor disks from 1 nG to 100 muG. At the time of the best match with the central region of the Antennae system the magnetic field has been amplified by compression and shear flows to an equilibrium field of approximately 10 muG, independent of the initial seed field. These simulations are a proof of the principle that galaxy mergers are efficient drivers for the cosmic evolution of magnetic fields. We present a detailed analysis of the magnetic field structure in the central overlap region. Simulated radio and polarization maps are in good morphological and quantitative agreement with the observations. In particular, the two cores with the highest synchrotron intensity and ridges of regular magnetic fields between the cores and at the root of the southern tidal arm develop naturally in our simulations. This indicates that the simulations are capable of realistically following the evolution of the magnetic fields in a highly non-linear environment. We also discuss the relevance of the amplification effect for present day magnetic fields in the context of hierarchical structure formation.

We present a new high-resolution N-body/SPH simulation of an encounter of two gas-rich disk galaxies which closely matches the morphology and kinematics of the interacting Antennae galaxies (NGC 4038/39). The simulation includes radiative cooling, star formation and feedback from SNII. The large-scale morphology and kinematics are determined by the internal structure and the orbit of the progenitor disks. The properties of the central region, in particular the starburst in the overlap region, only match the observations for a very short time interval of ~20 Myr after the second encounter. This indicates that the Antennae galaxies are in a special phase only about 40 Myr after the second encounter and 50 Myr before their final collision. This is the only phase in the simulation when a gas-rich overlap region between the nuclei is forming accompanied by enhanced star formation. The star formation rate as well as the recent star formation history in the central region agree well with observational estimates. For the first time this new model explains the distributed extra-nuclear star formation in the Antennae galaxies as a consequence of the recent second encounter. The proposed model predicts that the Antennae are in a later merger stage than the Mice (NGC 4676) and would therefore lose their first place in the classical Toomre sequence.

We present a first-principles' prediction that two charged particles of masses M_1 and M_2 separated R apart in a dielectric vacuum act on each other always an attractive force in addition to other known forces in between. This component attractive force on one charge results as the Lorentz force in the radiation depolarization- and magnetic- fields of the other charge, being an attractive radiation force, and is in addition to the ordinary repulsive radiation force. The exact solution for the attractive radiation force is F_g=G' M_1M_2/R^2, an identical formula to Newton's law of gravitation. G'=\chi_{0^*}e^4/4\pi\epsilon_0^2\hbar^2\rho_l is identifiable with Newton's gravitational constant, \chi_{0^*} being the susceptibility and \rho_l the linear mass density of the vacuum, and the remaining fundamental constants of the usual meaning. The F_g force is conveyed by a transverse vacuuonic dipole-moment wave traveling at the velocity of light and can penetrate matter freely. In all of respects, the F_g force represents a viable cause of Newton's universal gravity.

We study the effect of dissipational gas physics on the vertical heating and thickening of disc galaxies during minor mergers. We produce a suite of minor merger simulations for Milky Way-like galaxies. This suite consists of collisionless simulations as well as hydrodynamical runs including a gaseous component in the galactic disc. We find that in dissipationless simulations minor mergers cause the scale height of the disc to increase by up to a factor of ~2. When the presence of gas in the disc is taken into account this thickening is reduced by 25% (50%) for an initial disc gas fraction of 20% (40%), leading to a final scale height z0 between 0.6 and 0.7 kpc, regardless of the initial scale height. We argue that the presence of gas reduces disc heating via two mechanisms: absorption of kinetic impact energy by the gas and/or formation of a new thin stellar disc that can cause heated stars to recontract towards the disc plane. We show that in our simulations most of the gas is consumed during the merger and thus the regrowth of a new thin disc has a negligible impact on the z0 of the post merger galaxy. Final disc scale heights found in our simulations are in good agreement with studies of the vertical structure of spiral galaxies where the majority of the systems are found to have scale heights of 0.4 kpc < z0 < 0.8 kpc. We also found no tension between recent measurements of the scale height of the Milky Way thin disc and results coming from our hydrodynamical simulations. We conclude that the existence of a thin disc in the Milky Way and in external galaxies is not in obvious conflict with the predictions of the CDM model.

We study the evolution of black holes (BHs) on the M_BH-sigma and M_BH-M_bulge planes as a function of time in disk galaxies undergoing mergers. We begin the simulations with the progenitor black hole masses being initially below (Delta log M_BH=-2), on (Delta log M_BH=0) and above (Delta log M_BH=0.5) the observed local relations. The final relations are rapidly established after the final coalescense of the galaxies and their BHs. Progenitors with low initial gas fractions (f_gas=0.2) starting below the relations evolve onto the relations (Delta log M_BH=-0.18), progenitors on the relations stay there (Delta log M_BH=0) and finally progenitors above the relations evolve towards the relations, but still remaining above them (Delta log M_BH=0.35). Mergers in which the progenitors have high initial gas fractions (f_gas=0.8) evolve above the relations in all cases (Delta log M_BH=0.5). We find that the initial gas fraction is the prime source of scatter in the observed relations, dominating over the scatter arising from the evolutionary stage of the merger remnants. The fact that BHs starting above the relations do not evolve onto the relations, indicates that our simulations rule out the scenario in which overmassive BHs evolve onto the relations through gas-rich mergers. By implication our simulations thus disfavor the picture in which supermassive BHs develop significantly before their parent bulges.

Understanding the details of how the red sequence is built is a key question in galaxy evolution. What are the relative roles of gas-rich vs. dry mergers, major vs. minor mergers or galaxy mergers vs. gas accretion? In Wild et al. 2009 we compare hydrodynamic simulations with observations to show how gas-rich major mergers result in galaxies with strong post-starburst spectral features, a population of galaxies easily identified in the real Universe using optical spectra. Using spectra from the VVDS deep survey with z~0.7, and a principal component analysis technique to provide indices with high enough SNR, we find that 40% of the mass flux onto the red-sequence could enter through a strong post-starburst phase, and thus through gas-rich major mergers. The deeper samples provided by next generation galaxy redshift surveys will allow us to observe the primary physical processes responsible for the shut-down in starformation and build-up of the red sequence.

We find that the amount and nature of the assumed ionizing background can strongly affect galaxy formation and evolution. Galaxy evolution simulations typically incorporate an ultraviolet background which falls off rapidly above z=3; e.g., that of Haardt & Madau (1996). However, this decline may be too steep to fit the WMAP constraints on electron scattering optical depth or observations of intermediate redshift (z ~ 2-4) Ly-alpha forest transmission. As an alternative, we present simulations of the cosmological formation of individual galaxies with UV backgrounds that decline more slowly at high redshift: both a simple intensity rescaling and the background recently derived by Faucher-Giguere (2009), which softens the spectrum at higher redshifts. We also test an approximation of the X-ray background with a similar z-dependence. We find for the test galaxies that an increase in either the intensity or hardness of ionizing radiation generically pushes star formation towards lower redshifts: although overall star formation in the simulation boxes is reduced by 10-25%, the galaxies show a factor of ~2 increase in the fraction of stars within a 30 kpc radius that are formed after z=1. Other effects include late gas inflows enhanced up to 30 times, stellar half-mass radii decreased by up to 30%, central velocity dispersions increased up to 40%, and a strong reduction in substructure. The magnitude of the effects depends on the environmental/accretion properties of the particular galaxy.

Using a high resolution hydrodynamical cosmological simulation of the formation of a massive spheroidal galaxy we show that elliptical galaxies can be very compact and massive at high redshift in agreement with recent observations. Accretion of stripped in-falling stellar material increases the size of the system with time and the central concentration is reduced by dynamical friction of the surviving stellar cores. In a specific case of a spheroidal galaxy with a final stellar mass of $1.5 \times 10^{11} M_{\odot}$ we find that the effective radius $r_e$ increases from $0.7 \pm 0.2 \rm kpc$ at z = 3 to $r_e = 2.4 \pm 0.4 \rm kpc$ at z = 0 with a concomitant decrease in the effective density of an order of magnitude and a decrease of the central velocity dispersion by approximately 20% over this time interval. A simple argument based on the virial theorem shows that during the accretion of weakly bound material (minor mergers) the radius can increase as the square of the mass in contrast to the usual linear rate of increase for major mergers. By undergoing minor mergers compact high redshift spheroids can evolve into present-day systems with sizes and concentrations similar to observed local ellipticals. This indicates that minor mergers may be the main driver for the late evolution of sizes and densities of early-type galaxies.

We present a set of global, self-consistent N-body/SPH simulations of the dynamic evolution of galactic discs with gas and including magnetic fields. We have implemented a description to follow the evolution of magnetic fields with the ideal induction equation in the SPH part of the Vine code. Results from a direct implementation of the field equations are compared to a representation by Euler potentials, which pose a div(B)-free description, an constraint not fulfilled for the direct implementation. All simulations are compared to an implementation of magnetic fields in the Gadget code which includes also cleaning methods for div(B). Starting with a homogeneous seed field we find that by differential rotation and spiral structure formation of the disc the field is amplified by one order of magnitude within five rotation periods of the disc. The amplification is stronger for higher numerical resolution. Moreover, we find a tight connection of the magnetic field structure to the density pattern of the galaxy in our simulations, with the magnetic field lines being aligned with the developing spiral pattern of the gas. Our simulations clearly show the importance of non-axisymmetry for the evolution of the magnetic field.

We study the thermal formation history of four simulated galaxies that were shown in Naab et al. (2007) to reproduce a number of observed properties of elliptical galaxies. The temperature of the gas in the galaxies is steadily increasing with decreasing redshift, although much of the gas has a cooling time shorter than the Hubble time. The gas is being heated and kept hot by gravitational heating processes through the release of potential energy from infalling stellar clumps. The energy is dissipated in supersonic collisions of infalling gas lumps with the ambient gas and through the dynamical capturing of satellite systems causing gravitational wakes that transfer energy to the surrounding gas. Furthermore dynamical friction from the infalling clumps pushes out dark matter, lowering the central dark matter density by up to a factor of two from z=3 to z=0. In galaxies in which the late formation history (z<2) is dominated by minor merging and accretion the energy released (E~5x10^{59} ergs) from gravitational feedback is sufficient to form red and dead elliptical galaxies by z~1 even in the absence of supernova and AGN feedback.

From the VIMOS VLT DEEP Survey (VVDS) we select a sample of 16 galaxies with spectra which identify them as having recently undergone a strong starburst and subsequent fast quenching of star formation. These post-starburst galaxies lie in the redshift range 0.5<z<1.0 with masses >10^9.75Msun. They have a number density of 1x10^-4 per Mpc^3, almost two orders of magnitude sparser than the full galaxy population with the same mass limit. We compare with simulations to show that the galaxies are consistent with being the descendants of gas rich major mergers. Starburst mass fractions must be larger than ~5-10% and decay times shorter than ~10^8 years for post-starburst spectral signatures to be observed in the simulations. We find that the presence of black hole feedback does not greatly affect the evolution of the simulated merger remnants through the post-starburst phase. The multiwavelength spectral energy distributions of the post-starburst galaxies show that 5/16 have completely ceased the formation of new stars. These 5 galaxies correspond to a mass flux entering the red-sequence of rhodot(A->Q, PSB) = 0.0038Msun/Mpc^3/yr, assuming the defining spectroscopic features are detectable for 0.35Gyr. If the galaxies subsequently remain on the red sequence, this accounts for 38(+4/-11)% of the growth rate of the red sequence. Finally, we compare our high redshift results with a sample of galaxies with 0.05<z<0.1 observed in the SDSS and UKIDSS surveys. We find a very strong redshift evolution: the mass density of strong post-starburst galaxies is 230 times lower at z~0.07 than at z~0.7.

In the last decade we have seen an enormous increase in the size and quality of spectroscopic galaxy surveys, both at low and high redshift. New statistical techniques to analyse large portions of galaxy spectra are now finding favour over traditional index based methods. Here we will review a new robust and iterative Principal Component Analysis (PCA) algorithm, which solves several common issues with classic PCA. Application to the 4000AA break region of galaxies in the VIMOS VLT Deep Survey (VVDS) and Sloan Digital Sky Survey (SDSS) gives new high signal-to-noise ratio spectral indices easily interpretable in terms of recent star formation history. In particular, we identify a sample of post-starburst galaxies at z~0.7 and z~0.07. We quantify for the first time the importance of post-starburst galaxies, consistent with being descendants of gas-rich major mergers, for building the red sequence. Finally, we present a comparison with new low and high redshift "mock spectroscopic surveys" derived from a Millennium Run semi-analytic model.

In the framework of hierarchical structure formation ellipticals can form from merging of smaller disk galaxies. The nearby interacting 'Antennae' galaxy pair (NGC 4038/39) is one of the best-studied local systems of merging spirals, thus presenting us with an ideal laboratory for the study of galaxy evolution models. The Antennae are believed to be in a state prior to their final encounter with rapid subsequent merging, which puts them in the first position of the Toomre (1977) merger sequence. Here we present first numerical high-resolution, self-consistent, smoothed particle hydrodynamics (SPH) simulations of the Antennae system, including star formation and stellar feedback, and compare our results to VLA HI observations by Hibbard et al. (2001). We are able to obtain a close, but not yet perfect match to the observed morphology and kinematics of the system.

We present binary galaxy merger simulations of gas-rich disks (Sp-Sp), of early-type galaxies and disks (E-Sp, mixed mergers), and mergers of early-type galaxies (E-E, dry mergers) with varying mass ratios and different progenitor morphologies. The simulations include radiative cooling, star formation and black hole (BH) accretion and the associated feedback processes. We find for Sp-Sp mergers, that the peak star formation rate and BH accretion rate decrease and the growth timescales of the central black holes and newly formed stars increase with higher progenitor mass ratios. The termination of star formation by BH feedback in disk mergers is significantly less important for higher progenitor mass ratios (e.g. 3:1 and higher). In addition, the inclusion of BH feedback suppresses efficiently star formation in dry E-E mergers and mixed E-Sp mergers.

We present binary galaxy merger simulations with varying mass ratios and different progenitor morphologies. The simulations include mergers of gas-rich disks (Sp-Sp), of early-type galaxies and disks (E-Sp, mixed mergers), and mergers of early-type galaxies (E-E, dry mergers). We follow the dynamics of gas, stars and dark matter, and include radiative cooling, star formation and black hole (BH) accretion. For Sp-Sp mergers, the peak star formation rate and BH accretion rate decrease and the growth timescales of the central black holes and newly formed stars increase with higher progenitor mass ratios. The peak BH accretion rate typically occurs shortly after the time of BH merging for low progenitor mass ratios (e.g. 3:1 and lower), whereas for higher progenitor mass ratios there is no clear correlation between the peak BH accretion rate and BH merging time. The termination of star formation by BH feedback in disk mergers is significantly less important for higher progenitor mass ratios (e.g. 3:1 and higher). In addition, the inclusion of BH feedback suppresses efficiently star formation in dry E-E mergers and mixed E-Sp mergers. All merger remnants, independent of their progenitors, follow the observed relations between the central BH mass and the stellar velocity dispersion M_BH-sigma, the bulge mass M_BH-M_* and the bulge binding energy M_BH-M_{*}sigma^2, with the dominant source of scatter arising from variations in the initial gas mass fraction. The normalizations for all relations and the simulated slope of the M_BH-sigma and M_BH-M_{*}sigma^2 relations are in good agreement with the observations, whereas the simulated slope of the M_BH-M_* relation is slightly steeper compared to the observations. (abridged)

We present a simple set of kinematic criteria that can distinguish between galaxies dominated by ordered rotational motion and those involved in major merger events. Our criteria are based on the dynamics of the warm ionized gas (as traced by H-alpha) within galaxies, making this analysis accessible to high-redshift systems, whose kinematics are primarily traceable through emission features. Using the method of kinemetry (developed by Krajnovic and co-workers), we quantify asymmetries in both the velocity and velocity dispersion maps of the warm gas, and the resulting criteria enable us to empirically differentiate between non-merging and merging systems at high redshift. We apply these criteria to 11 of our best-studied rest-frame UV/optical-selected z~2 galaxies for which we have near infrared integral field spectroscopic data from SINFONI on the VLT. Of these 11 systems, we find that >50% have kinematics consistent with a single rotating disk interpretation, while the remaining systems are more likely undergoing major mergers. This result, combined with the short formation timescales of these systems, provides evidence that rapid, smooth accretion of gas plays a significant role in galaxy formation at high redshift.

(abridged) We present near-infrared (NIR) adaptive optics imaging obtained with VLT/NACO and optical spectroscopy from the Southern African Large Telescope (SALT) of a luminous infrared galaxy (LIRG) IRAS 19115-2124. These data are combined with archival HST imaging and Spitzer imaging and spectroscopy, allowing us to study this disturbed interacting/merging galaxy, dubbed the Bird, in extraordinary detail. In particular, the data reveal a triple system where the LIRG phenomenon is dominated by the smallest of the components. One nucleus is a regular barred spiral with significant rotation, while another is highly disturbed with a surface brightness distribution intermediate to that of disk and bulge systems, and hints of remaining arm/bar structure. We derive dynamical masses in the range 3-7x10^10 M_solar for both. The third component appears to be a 1-2x10^10 M_solar mass irregular galaxy. The total system exhibits HII galaxy-like optical line ratios and strengths, and no evidence for AGN activity is found from optical or mid-infrared data. The star formation rate is estimated to be 190 M_solar/yr. We search for SNe, super star clusters, and detect 100-300 km/s outflowing gas from the Bird. Overall, the Bird shows kinematic, dynamical, and emission line properties typical for cool ultra luminous IR galaxies. However, the interesting features setting it apart for future studies are its triple merger nature, and the location of its star formation peak - the strongest star formation does not come from the two major K-band nuclei, but from the third irregular component. Aided by simulations, we discuss scenarios where the irregular component is on its first high-speed encounter with the more massive components.

The intrinsic anisotropy delta and flattening epsilon of simulated merger remnants is compared with elliptical galaxies that have been observed by the SAURON collaboration, and that were analysed using axisymmetric Schwarzschild models. Collisionless binary mergers of stellar disks and disk mergers with an additional isothermal gas component, neglecting star formation cannot reproduce the observed trend delta = 0.55 epsilon (SAURON relationship). An excellent fit of the SAURON relationship for flattened ellipticals with epsilon >= 0.25 is however found for merger simulations of disks with gas fractions >= 20%, including star formation and stellar energy feedback. Massive black hole feedback does not strongly affect this result. Subsequent dry merging of merger remnants however does not generate the slowly-rotating SAURON ellipticals which are characterized by low ellipticities epsilon < 0.25 and low anisotropies. This indicates that at least some ellipticals on the red galaxy sequence did not form by binary mergers of disks or early-type galaxies. We show that stellar spheroids resulting from multiple, hierarchical mergers of star-bursting subunits in a cosmological context are in excellent agreement with the low ellipticities and anisotropies of the slowly rotating SAURON ellipticals and their observed trend of delta with epsilon. The numerical simulations indicate that the SAURON relation might be a result of strong violent relaxation and phase mixing of multiple, kinematically cold stellar subunits with the angular momentum of the system determining its location on the relation.

Galactic spheroids can form as a result of galaxy interactions and mergers of disks. Detailed analyses of the photometric properties, the intrinsic orbital structure, the line-of-sight velocity distributions and the kinemetry of simulated merger remnants, which depend critically on the geometry and the gas content of the interacting progenitors, indicate that low and intermediate mass rotating ellipticals can form from mergers of disks. The masses and metallicities of all massive ellipticals and the kinematics of some massive non-rotating ellipticals cannot be explained by binary mergers. Thus these galaxies might have formed in a different way.

We describe high resolution Smoothed Particle Hydrodynamics (SPH) simulations of three approximately $M_*$ field galaxies starting from \LCDM initial conditions. The simulations are made intentionally simple, and include photoionization, cooling of the intergalactic medium, and star formation but not feedback from AGN or supernovae. All of the galaxies undergo an initial burst of star formation at $z \approx 5$, accompanied by the formation of a bubble of heated gas. Two out of three galaxies show early-type properties at present whereas only one of them experienced a major merger. Heating from shocks and -PdV work dominates over cooling so that for most of the gas the temperature is an increasing function of time. By $z \approx 1$ a significant fraction of the final stellar mass is in place and the spectral energy distribution resembles those of observed massive red galaxies. The galaxies have grown from $z=1 \to 0$ on average by 25% in mass and in size by gas poor (dry) stellar mergers. By the present day, the simulated galaxies are old ($\approx 10 {\rm Gyrs}$), kinematically hot stellar systems surrounded by hot gaseous haloes. Stars dominate the mass of the galaxies up to $\approx 4$ effective radii ($\approx 10$ kpc). Kinematic and most photometric properties are in good agreement with those of observed elliptical galaxies. The galaxy with a major merger develops a counter-rotating core. Our simulations show that realistic intermediate mass giant elliptical galaxies with plausible formation histories can be formed from \LCDM initial conditions even without requiring recent major mergers or feedback from supernovae or AGN.

We apply a physically motivated stellar feedback model to analyse the statistical properties of damped Lyman-alpha systems (DLAs) expected in the concordance cold dark matter (CDM) model. Our feedback model produces extended low-metallicity cold gaseous discs around small galaxies. Since the space density of galaxies with low circular speeds is high, these discs dominate the cross-section for the identification of DLAs at all redshifts. The combined effects of star formation, outflows and infall in our models result in mild evolution of the N_{HI}-weighted metallicity content in DLAs with redshift, consistent with observations. According to our model, DLAs contribute only a small fraction of the volume averaged star formation rate at redshifts z \simlt 5. Our model predicts weak evolution in Omega_{HI} over the redshift range z=0-5. Furthermore, we show that the cosmological evolution of Omega_{HI} and the cosmic star formation rate are largely disconnected and conclude that the evolution of Omega_{HI} as a function of redshift is more likely to tell us about feedback processes and the evolution of the outer gaseous components of small galaxies than about the cosmic history of star formation.

Using N-body simulations with a large set of massless test particles we compare the predictions of two theories of violent relaxation, the well known Lynden-Bell theory and the more recent theory by Nakamura. We derive ``weaken'' versions of both theories in which we use the whole equilibrium coarse-grained distribution function as a constraint instead of the total energy constraint. We use these weaken theories to construct expressions for the conditional probability $K_i(\tau)$ that a test particle initially at the phase-space coordinate $\tau$ would end-up in the $i$'th macro-cell at equilibrium. We show that the logarithm of the ratio $R_{ij}(\tau) \equiv K_i(\tau)/K_j(\tau)$ is directly proportional to the initial phase-space density $f_0(\tau)$ for the Lynden-Bell theory and inversely proportional to $f_0(\tau)$ for the Nakamura theory. We then measure $R_{ij}(\tau)$ using a set of N-body simulations of a system undergoing a gravitational collapse to check the validity of the two theories of violent relaxation. We find that both theories are at odds with the numerical results, qualitatively and quantitatively.

We present a sample of mid-infrared detected sources from the European Large Area ISO Survey (ELAIS) regions characterised by strong mid-IR radiation with faint near-IR and optical counterparts. These extreme mid-to-near-IR objects (EMNOs) are defined here by a flux ratio of f_15um / f_2.2um > 25. This population is not obvious in deeper small area ISO surveys, though it produces more than 20% of the observed cosmic IR background radiation (CIRB) at 15um above 1 mJy. Near-future large area deep mid-IR surveys with the Spitzer Space Telescope, however, are bound to uncover large amounts of these objects, which we argue to most likely be obscured AGN, based on SED shapes and X-ray data. Very strong dusty starbursts at z>1 may also have high mid-to-near-IR flux ratios, but using the MIR/NIR and FIR/MIR ratios these may be separated. Most of our EMNOs appear to be ULIRGs, half are also extremely red objects (ERO). A curious case of a low redshift, less luminous object with a very young stellar population is also found. We predict that the simple broad band selection method makes EMNOs a useful window into high-redshift obscured nuclear activity and its sought after relation to star-formation, in a similar way that EROs have been used to define samples of high-redshift early type galaxies.

We have searched for Extremely Red Objects (EROs) around faint mid-IR selected galaxies in ELAIS fields. We find a significant overdensity, by factors of 2 to 5, of these EROs compared to field EROs in the same region and literature random field ERO counts. The excess is similar to that found previously in the fields of known high redshift quasars and AGN. While with the present data it cannot be definitely shown whether the overdensity is physically connected to the mid-IR source, a derived radial distribution does suggest this. The fraction of EROs among K-selected galaxies is high in the mid-IR fields in agreement with the picture that the EROs responsible for the overdensity are members of high redshift clusters of galaxies. We find R-K>5 selected EROs to be more clustered around the mid-IR galaxies than I-K>4 EROs, though statistics are weak because of small numbers. However, this would be consistent with a cluster/galaxy group scenario if, as we argue, the former colour selection finds preferentially more strongly clustered early type galaxies, whereas the latter selection includes a larger fraction of dusty EROs. Finally, using the mid-IR data, we are able to limit the fraction of ULIRG type very dusty EROs at K<18 magnitude to less than 10% of the total ERO population.

We present results on number counts of Extremely Red Objects (EROs) in a 2850 arcmin^2 near-infrared survey performed in European Large Area ISO Survey (ELAIS) fields at K<17.5. Counts of EROs are extended to brighter levels than available previously, giving 0.002 +/- 0.001 arcmin^-2 at K<16.5 and consistent numbers with literature values at fainter magnitudes. Photometric redshifts from HYPERZ as well as GRASIL model SEDs of galaxies imply that our EROs are located in the range z=0.7-1.5, with the bulk of the population at z~1. Taking advantage of the ISO data in the fields, we use mid-IR detections to constrain the number of dusty EROs, and also discuss the superior capabilities of Spitzer Space Telescope to detect dusty EROs. Both the mid-IR data and the use of colour-colour diagrammes indicate that at most 10-20% of the EROs in this bright regime are dusty starbursting systems. The space density of our EROs, interpreted to be counterparts of local >2-3L^star massive galaxies at around z~1, is estimated to be approx. 2x10^-5 Mpc^-3, which is consistent with local values. Furthermore, the cumulative number counts at our bright magnitudes are remarkably well fitted by pure luminosity evolution models.