We present a hybrid code combining the OpenMP-parallel tree code VINE with an algorithmic chain regularization scheme. The new code, called "rVINE", aims to significantly improve the accuracy of close encounters of massive bodies with supermassive black holes in galaxy-scale numerical simulations. We demonstrate the capabilities of the code by studying two test problems, the sinking of a single massive black hole to the centre of a gas-free galaxy due to dynamical friction and the hardening of a supermassive black hole binary due to close stellar encounters. We show that results obtained with rVINE compare well with NBODY7 for problems with particle numbers that can be simulated with NBODY7. In particular, in both NBODY7 and rVINE we find a clear N-dependence of the binary hardening rate, a low binary eccentricity and moderate eccentricity evolution, as well as the conversion of the galaxy's inner density profile from a cusp to a a core via the ejection of stars at high velocity. The much larger number of particles that can be handled by rVINE will open up exciting opportunities to model stellar dynamics close to SMBHs much more accurately in a realistic galactic context. This will help to remedy the inherent limitations of commonly used tree solvers to follow the correct dynamical evolution of black holes in galaxy scale simulations.

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.

We reexamine the age distribution of star clusters in the Antennae in the context of N-body+hydrodynamical simulations of these interacting galaxies. All of the simulations that account for the observed morphology and other properties of the Antennae have star formation rates that vary relatively slowly with time, by factors of only 1.3 - 2.5 in the past 10^8 yr. In contrast, the observed age distribution of the clusters declines approximately as a power law, dN/dt \propto t^{gamma} with gamma = -1.0, for ages 10^6 yr \la t \la 10^9 yr. These two facts can only be reconciled if the clusters are disrupted progressively for at least 10^8 yr and possibly 10^9 yr. When we combine the simulated formation rates with a power-law model, f_surv \propto t^{delta}, for the fraction of clusters that survive to each age t, we match the observed age distribution with exponents in the range -0.9 \la delta \la -0.6 (with a slightly different delta for each simulation). The similarity between delta and gamma indicates that dN/dt is shaped mainly by the disruption of clusters rather than variations in their formation rate. Thus, the situation in the interacting Antennae resembles that in relatively quiescent galaxies such as the Milky Way and the Magellanic Clouds.

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

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.