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, espec
ially 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 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, st
ar 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.
