Whereas cold dark matter (CDM) simulations predict central dark matter cusps with densities that diverge as $\rho$(r)$\sim$ 1/r observations often indicate constant density cores with finite central densities $\rho_0$ and a flat density distribution within a core radius r$_0$. This paper investigates whether this core-cusp problem can be solved by fuzzy dark matter (FDM), a hypothetical particle with a mass of order m$\approx$10$^{-22}$eV and a corresponding de Broglie wavelength on astrophysical scales. We show that galaxies with CDM halo virial masses M$_{vir} \leq 10^{11}$M$_{\odot}$ follow two core scaling relations. In addition to the well known universal core column density $\Sigma_0 \equiv \rho_0 \times$r$_0$ = 75 M$_{\odot}$pc$^{-2}$ core radii increase with virial masses as r$_0 \sim$ M$_{vir}^{\gamma}$ with $\gamma$ of order unity. Using the simulations by Schive et al. (2014) we demonstrate that FDM can explain the r$_0$-M$_{vir}$ scaling relation if the virial masses of the observed galaxy sample scale with formation redshift z as M$_{vir}\sim$(1+z)$^{-0.4}$. The observed constant $\Sigma_0$ is however in complete disagreement with FDM cores which are characterised by a steep dependence $\Sigma_0 \sim$r$_0^{-3}$, independent of z. More high-resolution simulations are now required to confirm the Schive et al. simulations and explore especially the transition region between the soliton core and the surrounding halo. If these results hold, FDM can be ruled out as the origin of observed dark matter cores and other physical processes are required to account for their formation.

We study the formation, evolution and collapse of dense cores by tracking density structures in a magnetohydrodynamic (MHD) simulation. We identify cores using the dendrogram algorithm and utilize machine learning techniques, including principal component analysis (PCA) and the k-means clustering algorithm to analyze the full density and velocity dispersion profiles of these cores. We find that there exists an evolutionary sequence consisting of three distinct phases: i) the formation of turbulent density structures (Phase I), ii) the dissipation of turbulence and the formation of coherent cores (Phase II), and iii) the transition to protostellar cores through gravitational collapse (Phase III). In dynamically evolving molecular clouds, the existence of these three phases corresponds to the coexistence of three populations of cores with distinct physical properties. The prestellar and protostellar cores frequently analyzed in previous studies of observations and simulations belong to the last phase in this evolutionary picture. We derive typical lifetimes of 1.4$\pm$1.0$\times$10$^5$ yr, 3.3$\pm$1.4$\times$10$^5$ yr and 3.3$\pm$1.4$\times$10$^5$ yr, respectively for Phase I, II and III. We find that cores can form from both converging flows and filament fragmentation and that cores may form both inside and outside the filaments. We then compare our results to previous observations of coherent cores and provide suggestions for future observations to study cores belonging to the three phases.

We report high quality, H{\alpha} or CO rotation curves (RCs) to several Re for 41 large, massive, star-forming disk galaxies (SFGs), across the peak of cosmic galaxy evolution (z~0.67-2.45), taken with the ESO-VLT, the LBT and IRAM-NOEMA. Most RC41 SFGs have reflection symmetric RCs plausibly described by equilibrium dynamics. We fit the major axis position-velocity cuts with beam-convolved, forward modeling with a bulge, a turbulent rotating disk, and a dark matter (DM) halo. We include priors for stellar and molecular gas masses, optical light effective radii and inclinations, and DM masses from abundance matching scaling relations. Two-thirds or more of the z>1.2 SFGs are baryon dominated within a few Re of typically 5.5 kpc, and have DM fractions less than maximal disks (<fDM (Re)>=0.12). At lower redshift (z<1.2) that fraction is less than one-third. DM fractions correlate inversely with the baryonic angular momentum parameter, baryonic surface density and bulge mass. Inferred low DM fractions cannot apply to the entire disk & halo but more plausibly reflect a flattened, or cored, inner DM density distribution. The typical central 'DM deficit' in these cores relative to NFW distributions is ~30% of the bulge mass. The observations are consistent with rapid radial transport of baryons in the first generation massive gas rich halos forming globally gravitationally unstable disks, and leading to efficient build-up of massive bulges and central black holes. A combination of heating due to dynamical friction and AGN feedback may drive DM out of the initial cusps.

We present ALMA maps of the starless molecular cloud core Ophiuchus/H-MM1 in the lines of deuterated ammonia (ortho-NH2D), methanol (CH3OH), and sulphur monoxide (SO). The dense core is seen in NH2D emission, whereas the CH3OH and SO distributions form a halo surrounding the core. Because methanol is formed on grain surfaces, its emission highlights regions where desorption from grains is particularly efficient. Methanol and sulphur monoxide are most abundant in a narrow zone that follows the eastern side of the core. This side is sheltered from the stronger external radiation field coming from the west. We show that photodissociation on the illuminated side can give rise to an asymmetric methanol distribution, but that the stark contrast observed in H-MM1 is hard to explain without assuming enhanced desorption on the shaded side. The region of the brightest emission has a wavy structure that rolls up at one end. This is the signature of Kelvin-Helmholtz instability occurring in sheared flows. We suggest that in this zone, methanol and sulphur are released as a result of grain-grain collisions induced by shear vorticity.

We extend our previous work on simulations with the code RAMSES on accretion driven turbulence by including self-gravity and study the effects of core formation and collapse. We show that radial accretion onto filaments drives turbulent motions which are not isotropic but radially dominated. In contrast to filaments without gravity, the velocity dispersion of self-gravitating filaments does not settle in an equilibrium. Despite showing similar amounts of driven turbulence, they continually dissipate their velocity dispersion until the onset of core formation. This difference is connected to the evolution of the radius as it determines the dissipation rate. In the non-gravitational case filament growth is not limited and its radius grows linearly with time. In contrast, there is a maximum extent in the self-gravitational case resulting in an increased dissipation rate. Furthermore, accretion driven turbulence shows a radial profile which is anti-correlated with density. This leads to a constant turbulent pressure throughout the filament. As the additional turbulent pressure does not have a radial gradient it does not contribute to the stability of filaments and does not increase the critical line-mass. However, this radial turbulence does affect the radius of a filament, adding to the extent and setting its maximum value. Moreover, the radius evolution also affects the growth timescale of cores which compared to the timescale of collapse of an accreting filament limits core formation to high line-masses.

We study the spatially-resolved stellar specific angular momentum $j_*$ in a high-quality sample of 24 CALIFA galaxies covering a broad range of visual morphology, accounting for stellar velocity and velocity dispersion. The shape of the spaxel-wise probability density function of normalised $s=j_*/j_{*mean}$, PDF($s$), deviates significantly from the near-universal initial distribution expected of baryons in a dark matter halo and can be explained by the expected baryonic effects in galaxy formation that remove and redistribute angular momentum. Further we find that the observed shape of the PDF($s$) correlates significantly with photometric morphology, where late-type galaxies have PDF($s$) that is similar to a normal distribution, whereas early types have a strongly-skewed PDF($s$) resulting from an excess of low-angular momentum material. Galaxies that are known to host pseudobulges (bulge Sersic index $n_b <2.2$) tend to have less skewed bulge PDF($s$), with skewness $(b_{1rb})\lesssim0.8$. The PDF($s$) encodes both kinematic and photometric information and appears to be a robust tracer of morphology. Its use is motivated by the desire to move away from traditional component-based classifications which are subject to observer bias, to classification on a galaxy's fundamental (stellar mass, angular momentum) properties. In future, PDF($s$) may also be useful as a kinematic decomposition tool.

We present half-light sizes measured from H${\alpha}$ emission tracing star-formation in 281 star-forming galaxies from the KMOS3D survey at 0.7 < z < 2.7. Sizes are derived by fitting 2D exponential disk models, with bootstrap errors averaging 20%. H${\alpha}$ sizes are a median (mean) of 1.19 (1.26) times larger than those of the stellar continuum, which due to radial dust gradients places an upper limit on the growth in stellar size via star formation, with just 43% intrinsic scatter. At fixed continuum size the H${\alpha}$ size shows no residual trend with stellar mass, star formation rate, redshift or morphology. The only significant residual trend is with the excess obscuration of H${\alpha}$ by dust, at fixed continuum obscuration. The scatter in continuum size at fixed stellar mass is likely driven by the scatter in halo spin parameters. The stability of the ratio of H${\alpha}$ size to continuum size demonstrates a high degree of stability in halo spin and in the transfer of angular momentum to the disk over a wide range of physical conditions and cosmic time. This may require local regulation by feedback processes. The implication of our results, as we demonstrate using a toy model, is that our upper limit on star-formation driven growth is sufficient only to evolve star-forming galaxies approximately along the observed size-mass relation, consistent with the size growth of galaxies at constant cumulative co-moving number density. To explain the observed evolution of the size-mass relation of star-forming disk galaxies other processes, such as the preferential quenching of compact galaxies or galaxy mergers, may be required.

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.

The inner parsec of our Galaxy contains tens of Wolf-Rayet stars whose powerful outflows are constantly interacting while filling the region with hot, diffuse plasma. Theoretical models have shown that, in some cases, the collision of stellar winds can generate cold, dense material in the form of clumps. However, their formation process and properties are not well understood yet. In this work we present, for the first time, a statistical study of the clump formation process in unstable wind collisions. We study systems with dense outflows (${\sim}10^{-5}\rm\ M_{\odot}\ yr^{-1}$), wind speeds of $500$-$1500\rm\ km\ s^{-1}$, and stellar separations of ${\sim}20$-$200\rm\ au$. We develop 3D high resolution hydrodynamical simulations of stellar wind collisions with the adaptive-mesh refinement grid-based code Ramses. We aim to characterise the initial properties of clumps that form through hydrodynamic instabilities, mostly via the non-linear thin shell instability (NTSI). Our results confirm that more massive clumps are formed in systems whose winds are close to the transition between the radiative and adiabatic regimes. Increasing either the wind speed or the degree of asymmetry increases the dispersion of the clump mass and ejection speed distributions. Nevertheless, the most massive clumps are very light (${\sim}10^{-3}$-$10^{-2}\rm\ M_{\oplus}$), about three orders of magnitude less massive than theoretical upper limits. Applying these results to the Galactic Centre we find that clumps formed through the NTSI should not be heavy enough either to affect the thermodynamic state of the region or to survive for long enough to fall onto the central super-massive black hole.

We investigate the stellar kinematics of a sample of galaxies extracted from the hydrodynamic cosmological Magneticum Pathfinder simulations out to $5$ half-mass radii. We construct differential radial stellar spin profiles quantified by the observationally widely used $\lambda_\mathrm{R}$ and the closely related $(V/\sigma)$ parameters. We find three characteristic profile shapes: profiles exhibiting a (i) peak within $2.5$ half-mass radii and a subsequent decrease (ii) continuous increase that plateaus at larger radii typically with a high amplitude (iii) completely flat behaviour typically with low amplitude, in agreement with observations. This shows that the kinematic state of the stellar component can vary significantly with radius, suggesting a distinct interplay between in-situ star formation and ex-situ accretion of stars. Following the evolution of our sample through time, we provide evidence that the accretion history of galaxies with decreasing profiles is dominated by the anisotropic accretion of low mass satellites that get disrupted beyond $ \sim 2.0$ half-mass radii, building up a stellar halo with non-ordered motion while maintaining the central rotation already present at $z=2$. In fact, at $z=2$ decreasing profiles are the predominant profile class. Hence, we can predict a distinct formation pathway for galaxies with a decreasing profile and show that the centre resembles an old embedded disk. Furthermore, we show that the radius of the kinematic transition provides a good estimation for the transition radius from in-situ stars in the centre to accreted stars in the halo.

The central super-massive black hole of the Milky Way, Sgr A*, accretes at a very low rate making it a very underluminous galactic nucleus. Despite the tens of Wolf-Rayet stars present within the inner parsec supplying ${\sim}10^{-3}\rm\ M_{\odot}\ yr^{-1}$ in stellar winds, only a negligible fraction of this material ($<10^{-4}$) ends up being accreted onto Sgr A*. The recent discovery of cold gas (${\sim}10^4\rm\ K$) in its vicinity raised questions about how such material could settle in the hostile (${\sim}10^7\rm\ K$) environment near Sgr A*. In this work we show that the system of mass-losing stars blowing winds can naturally account for both the hot, inefficient accretion flow, as well as the formation of a cold disk-like structure. We run hydrodynamical simulations using the grid-based code Ramses starting as early in the past as possible to observe the state of the system at the present time. Our results show that the system reaches a quasi-steady state in about ${\sim}500\rm\ yr$ with material being captured at a rate of ${\sim}10^{-6}\rm\ M_{\odot}\ yr^{-1}$ at scales of ${\sim}10^{-4}\rm\ pc$, consistent with the observations and previous models. However, on longer timescales ($\gtrsim3000\rm\ yr$) the material accumulates close to the black hole in the form of a disk. Considering the duration of the Wolf-Rayet phase (${\sim}10^5\rm\ yr$), we conclude that this scenario likely has already happened, and could be responsible for the more active past of Sgr A*, and/or its current outflow. We argue that the hypothesis of the mass-losing stars being the main regulator of the activity of the black hole deserves further consideration.

We confirm that the number of globular clusters (GCs), N$_{GC}$, is an excellent tracer of their host galaxy's halo virial mass M$_{vir}$. The simple linear relation M$_{vir} = 5 \times 10^9$ M$_{\odot} \times$ N$_{GC}$ fits the data perfectly from M$_{vir} = 10^{10}$ M$_{\odot}$ to M$_{vir} = 2 \times 10^{15}$ M$_{\odot}$. This result is independent of galaxy morphology and extends statistically into the dwarf galaxy regime with M$_{vir} = 10^8 - 10^{10}$ M$_{\odot}$, including the extreme ultra diffuse galaxy DF44. As this correlation does not depend on GC mass it is ideally suited for high-precision determinations of M$_{vir}$. The linearity is most simply explained by cosmological merging of a high-redshift halo seed population that hosted on average one GC per $5 \times 10^8$ M$_{\odot}$ of dark matter. We show that hierarchical merging is also extremely powerful in restoring a linear correlation and erasing signatures of even a strong secular evolution of GC systems. The cosmological merging scenario also implies a strong decline of the scatter in $N_{GC}$ with increasing virial mass $\delta N_{GC}/N_{GC} \sim M_{vir}^{-1/2}$ in contrast with the observations that show a roughly constant scatter, independent of virial mass. This discrepancy can be explained if errors in determining virial masses from kinematical tracers and gravitational lensing are on the order of a factor of 2. GCs in dwarf satellite galaxies pose a serious problem for high-redshift GC formation scenarios; the dark halo masses of dwarf galaxies hosting GCs therefore might need to be an order of magnitude larger than currently estimated.

We present an analysis of the internal velocity structures of the newly identified sub-0.1 pc coherent structures, droplets, in L1688 and B18. By fitting 2D linear velocity fields to the observed maps of velocity centroids, we determine the magnitudes of linear velocity gradients and examine the potential rotational motions that could lead to the observed velocity gradients. The results show that the droplets follow the same power-law relation between the velocity gradient and size found for larger-scale dense cores. Assuming that rotational motion giving rise to the observed velocity gradient in each core is a solid-body rotation of a rotating body with a uniform density, we derive the "net rotational motions" of the droplets. We find a ratio between rotational and gravitational energies, $\beta$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $\beta \sim 0.039$. We then examine the alignment between the velocity gradient and the major axis of each droplet, using methods adapted from the histogram of relative orientations (HRO) introduced by Soler et al. (2013). We find no definitive correlation between the directions of velocity gradients and the elongations of the cores. Lastly, we discuss physical processes other than rotation that may give rise to the observed velocity field.

We use a dedicated 0.7-m telescope to image the halos of 119 galaxies in the Local Volume to $\mu_r \sim 28-30$ mag/arcsec$^2$. The sample is primarily from the 2MASS Large Galaxy Atlas and extended to include nearby dwarf galaxies and more distant giant ellipticals, and spans fully the galaxy colour-magnitude diagram including the blue cloud and red sequence. We present an initial overview, including deep images of our galaxies. Our observations reproduce previously reported low surface brightness structures, including extended plumes in M51, and a newly discovered tidally extended dwarf galaxy in NGC7331. Low surface brightness structures, or "envelopes", exceeding 50 kpc in diameter are found mostly in galaxies with $M_V<-20.5$, and classic interaction signatures are infrequent. Defining a halo diameter at the surface brightness 28 mag/arcsec$^2$, we find that halo diameter is correlated with total galaxy luminosity. Extended signatures of interaction are found throughout the galaxy colour-magnitude diagram without preference for the red or blue sequences, or the green valley. Large envelopes may be found throughout the colour-magnitude diagram with some preference for the bright end of the red sequence. Spiral and S0 galaxies have broadly similar sizes, but ellipticals extend to notably greater diameters, reaching 150 kpc. We propose that the extended envelopes of disk galaxies are dominated by an extension of the disk population rather than by a classical population II halo.

The effect of galactic orbits on a galaxy's internal evolution within a galaxy cluster environment has been the focus of heated debate in recent years. To understand this connection, we use both the $(0.5 \,$Gpc)$^3$ and the Gpc$^3$ boxes from the cosmological hydrodynamical simulation set Magneticum Pathfinder. We investigate the velocity-anisotropy, phase space, and the orbital evolution of up to $\sim 5 \cdot 10^{5}$ resolved satellite galaxies within our sample of 6776 clusters with $M_{\mathrm{vir}} > 10^{14} \, \mathrm{M_{\odot}}$ at low redshift, which we also trace back in time. In agreement with observations, we find that star-forming satellite galaxies inside galaxy clusters are characterised by more radially dominated orbits, independent of cluster mass. Furthermore, the vast majority of star-forming satellite galaxies stop forming stars during their first passage. We find a strong dichotomy both in line-of-sight and radial phase space between star-forming and quiescent galaxies, in line with observations. The tracking of individual orbits shows that the star-formation of almost all satellite galaxies drops to zero within $1 \, \mathrm{Gyr}$ after in-fall. Satellite galaxies that are able to remain star-forming longer are characterised by tangential orbits and high stellar mass. All this indicates that in galaxy clusters the dominant quenching mechanism is ram-pressure stripping.

We investigate the contribution of clumps and satellites to the galaxy mass assembly. We analyzed spatially-resolved Hubble Space Telescope observations (imaging and slitless spectroscopy) of 53 star-forming galaxies at z ~ 1 - 3. We created continuum and emission line maps and pinpointed residual "blobs" detected after subtracting the galaxy disk. Those were separated into compact (unresolved) and extended (resolved) components. Extended components have sizes ~ 2 kpc and comparable stellar mass and age as the galaxy disks, whereas the compact components are 1.5 dex less massive and 0.4 dex younger than the disks. Furthermore the extended blobs are typically found at larger distances from the galaxy barycenter than the compact ones. Prompted by these observations and by the comparison with simulations, we suggest that compact blobs are in-situ formed clumps, whereas the extended ones are accreting satellites. Clumps and satellites enclose respectively ~ 20% and ~< 80% of the galaxy stellar mass, ~ 30% and ~ 20% of its star formation rate. Considering the compact blobs, we statistically estimated that massive clumps (Mstar >~ 10^9 Msun) have lifetimes of ~ 650 Myr, and the less massive ones (10^8 < Mstar < 10^9 Msun) of ~ 145 Myr. This supports simulations predicting long-lived clumps (lifetime > 100 Myr). Finally, ~< 30% (13%) of our sample galaxies are undergoing single (multiple) merger(s), they have a projected separation ~< 10 kpc, and the typical mass ratio of our satellites is 1:5 (but ranges between 1:10 and 1:1), in agreement with literature results for close pair galaxies.

We investigate the build-up of the galactic dynamo and subsequently the origin of a magnetic driven outflow. We use a setup of an isolated disc galaxy with a realistic circum-galactic medium (CGM). We find good agreement of the galactic dynamo with theoretical and observational predictions from the radial and toroidal components of the magnetic field as function of radius and disc scale height. We find several field reversals indicating dipole structure at early times and quadrupole structure at late times. Together with the magnetic pitch angle and the dynamo control parameters $R_{\alpha}$, $R_{\omega}$ and $D$ we present strong evidence for an $\alpha^2$-$\Omega$ dynamo. The formation of a bar in the centre leads to further amplification of the magnetic field via adiabatic compression which subsequently drives an outflow. Due to the Parker-Instability the magnetic field lines rise to the edge of the disc, break out and expand freely in the CGM driven by the magnetic pressure. Finally, we investigate the correlation between magnetic field and star formation rate. Globally, we find that the magnetic field is increasing as function of the star formation rate surface density with a slope between $0.3$ and $0.45$ in good agreement with predictions from theory and observations. Locally, we find that the magnetic field can decrease while star formation increases. We find that this effect is correlated with the diffusion of magnetic field from the spiral arms to the inter-arm regions which we explicitly include by solving the induction equation and accounting for non-linear terms.

Galaxy mergers are a fundamental part of galaxy evolution. To study the resulting mass distributions of different kinds of galaxy mergers, we present a simulation suite of 36 high-resolution isolated merger simulations, exploring a wide range of parameter space in terms of mass ratios (mu = 1:5, 1:10, 1:50, 1:100) and orbital parameters. We find that mini mergers deposit a higher fraction of their mass in the outer halo compared to minor mergers, while their contribution to the central mass distribution is highly dependent on the orbital impact parameter: for larger pericentric distances we find that the centre of the host galaxy is almost not contaminated by merger particles. We also find that the median of the resulting radial mass distribution for mini mergers differs significantly from the predictions of simple theoretical tidal-force models. Furthermore, we find that mini mergers can increase the size of the host disc significantly without changing the global shape of the galaxy, if the impact occurs in the disc plane, thus providing a possible explanation for extended low-surface brightness disks reported in observations. Finally, we find clear evidence that streams are a strong indication of nearly circular infall of a satellite (with large angular momentum), whereas the appearance of shells clearly points to (nearly) radial satellite infall.

Giant clumps on ~kpc scales and with masses of 10^8-10^9 Msol are ubiquitous in observed high-redshift disc galaxies. Recent simulations and observations with high spatial resolution indicate the existence of substructure within these clumps. We perform high-resolution simulations of a massive galaxy to study the substructure formation within the framework of gravitational disc instability. We focus on an isolated and pure gas disk with an isothermal equation of state with T=10^4 K that allows capturing the effects of self-gravity and hydrodynamics robustly. The main mass of the galaxy resides in rotationally supported clumps which grow by merging to a maximum clump mass of 10^8 Msol with diameter ~120 pc for the dense gas. They group to clump clusters (CCs) within relatively short times (<< 50 Myr), which are present over the whole simulation time. We identify several mass and size scales on which the clusters appear as single objects at the corresponding observational resolution between ~10^8 - 10^9 Msol. Most of the clusters emerge as dense groups and for larger beams they are more likely to be open structures represented by a single object. In the high resolution runs higher densities can be reached, and the initial structures can collapse further and fragment to many clumps smaller than the initial Toomre length. In our low resolution runs, the clumps directly form on larger scales 0.3-1 kpc with 10^8-10^9 Msol. Here, the artificial pressure floor which is typically used to prevent spurious fragmentation strongly influences the initial formation of clumps and their properties at very low densities.

We present the $0.6<z<2.6$ evolution of the ionized gas velocity dispersion in 175 star-forming disk galaxies based on data from the full KMOS$^{\rm 3D}$ integral field spectroscopic survey. In a forward-modelling Bayesian framework including instrumental effects and beam-smearing, we fit simultaneously the observed galaxy velocity and velocity dispersion along the kinematic major axis to derive the intrinsic velocity dispersion $\sigma_0$. We find a reduction of the average intrinsic velocity dispersion of disk galaxies as a function of cosmic time, from $\sigma_0\sim45$ km s$^{-1}$ at $z\sim2.3$ to $\sigma_0\sim30$ km s$^{-1}$ at $z\sim0.9$. There is substantial intrinsic scatter ($\sigma_{\sigma_0, {\rm int}}\approx10$ km s$^{-1}$) around the best-fit $\sigma_0-z$-relation beyond what can be accounted for from the typical measurement uncertainties ($\delta\sigma_0\approx12$ km s$^{-1}$), independent of other identifiable galaxy parameters. This potentially suggests a dynamic mechanism such as minor mergers or variation in accretion being responsible for the scatter. Putting our data into the broader literature context, we find that ionized and atomic+molecular velocity dispersions evolve similarly with redshift, with the ionized gas dispersion being $\sim10-15$ km s$^{-1}$ higher on average. We investigate the physical driver of the on average elevated velocity dispersions at higher redshift, and find that our galaxies are at most marginally Toomre-stable, suggesting that their turbulent velocities are powered by gravitational instabilities, while stellar feedback as a driver alone is insufficient. This picture is supported through comparison with a state-of-the-art analytical model of galaxy evolution.

We present the observation and analysis of newly discovered coherent structures in the L1688 region of Ophiuchus and the B18 region of Taurus. Using data from the Green Bank Ammonia Survey (GAS), we identify regions of high density and near-constant, almost-thermal, velocity dispersion. Eighteen coherent structures are revealed, twelve in L1688 and six in B18, each of which shows a sharp "transition to coherence" in velocity dispersion around its periphery. The identification of these structures provides a chance to study the coherent structures in molecular clouds statistically. The identified coherent structures have a typical radius of 0.04 pc and a typical mass of 0.4 Msun, generally smaller than previously known coherent cores identified by Goodman et al. (1998), Caselli et al. (2002), and Pineda et al. (2010). We call these structures "droplets." We find that unlike previously known coherent cores, these structures are not virially bound by self-gravity and are instead predominantly confined by ambient pressure. The droplets have density profiles shallower than a critical Bonnor-Ebert sphere, and they have a velocity (VLSR) distribution consistent with the dense gas motions traced by NH3 emission. These results point to a potential formation mechanism through pressure compression and turbulent processes in the dense gas. We present a comparison with a magnetohydrodynamic simulation of a star-forming region, and we speculate on the relationship of droplets with larger, gravitationally bound coherent cores, as well as on the role that droplets and other coherent structures play in the star formation process.

We present a census of ionized gas outflows in 599 normal galaxies at redshift 0.6<z<2.7, mostly based on integral field spectroscopy of Ha, [NII], and [SII] line emission. The sample fairly homogeneously covers the main sequence of star-forming galaxies with masses 9.0<log(M*/Msun)<11.7, and probes into the regimes of quiescent galaxies and starburst outliers. About 1/3 exhibits the high-velocity component indicative of outflows, roughly equally split into winds driven by star formation (SF) and active galactic nuclei (AGN). The incidence of SF-driven winds correlates mainly with star formation properties. These outflows have typical velocities of ~450 km/s, local electron densities of n_e~380 cm^-3, modest mass loading factors of ~0.1-0.2 at all galaxy masses, and energetics compatible with momentum driving by young stellar populations. The SF-driven winds may escape from log(M*/Msun)<10.3 galaxies but substantial mass, momentum, and energy in hotter and colder outflow phases seem required to account for low galaxy formation efficiencies in the low-mass regime. Faster AGN-driven outflows (~1000-2000 km/s) are commonly detected above log(M*/Msun)~10.7, in up to ~75% of log(M*/Msun)>11.2 galaxies. The incidence, strength, and velocity of AGN-driven winds strongly correlates with stellar mass and central concentration. Their outflowing ionized gas appears denser (n_e~1000 cm^-3), and possibly compressed and shock-excited. These winds have comparable mass loading factors as the SF-driven winds but carry ~10 (~50) times more momentum (energy). The results confirm our previous findings of high duty cycle, energy-driven outflows powered by AGN above the Schechter mass, which may contribute to star formation quenching.

A panchromatic view of the star forming ring and feeding process in the central kpc of the galaxy NGC 1097 is presented. The assembled IR to UV images at $\sim$10 pc resolution allow us to characterise the population of circa 250 clusters in the ring and disentangle the network of filaments of dust and gas that enshroud and feed them. The ring is a place of intermittent star bursts over the last 100 Myr. Four major episodes covering a proto-cluster phase of eleven mid-IR sources at the molecular clouds core, and two (three) previous bursts with a time separation of 20 - 30 Myr are identified. The extinction map of the inner few kpc resolves NGC1097's two major dust lanes in bundles of narrow, $<$25 pc width, filaments running along the galaxy's bar. As they approach the ring, some circularise along it, others curve to the centre to produce a nuclear spiral. We believe these are kpc-scale dust-gas streamers feeding the ring and the black-hole. The total mass in clusters formed in the ring in the last 100 Myr is $< 10^7\, \rm{M_\odot}$, i.e. $< 1\% $ of the $10^{9} M_\odot$ of molecular gas in the ring; yet, at its current star formation rate, $\sim1.8\, \rm{M_\odot \, yr^{-1}}$, an order of magnitude more in stellar mass should have been produced over that period. This means that the availability of gas in the ring is not the sole star formation driver, perhaps the rate at which dense gas accumulates in the ring is key.

We investigate the relationship between star formation activity and outflow properties on kiloparsec scales in a sample of 28 star forming galaxies at $z\sim$ 2-2.6, using adaptive optics assisted integral field observations from SINFONI on the VLT. The narrow and broad components of the H$\alpha$ emission are used to simultaneously determine the local star formation rate surface density ($\Sigma_{\rm SFR}$), and the outflow velocity $v_{\rm out}$ and mass outflow rate $\dot{M}_{\rm out}$, respectively. We find clear evidence for faster outflows with larger mass loading factors at higher $\Sigma_{\rm SFR}$. The outflow velocities scale as $v_{\rm out}$ $\propto$ $\Sigma_{\rm SFR}^{0.34 \pm 0.10}$, which suggests that the outflows may be driven by a combination of mechanical energy released by supernova explosions and stellar winds, as well as radiation pressure acting on dust grains. The majority of the outflowing material does not have sufficient velocity to escape from the galaxy halos, but will likely be re-accreted and contribute to the chemical enrichment of the galaxies. In the highest $\Sigma_{\rm SFR}$ regions the outflow component contains an average of $\sim$45% of the H$\alpha$ flux, while in the lower $\Sigma_{\rm SFR}$ regions only $\sim$10% of the H$\alpha$ flux is associated with outflows. The mass loading factor, $\eta$ = $\dot{M}_{\rm out}$/SFR, is positively correlated with $\Sigma_{\rm SFR}$ but is relatively low even at the highest $\Sigma_{\rm SFR}$: $\eta \lesssim$ 0.5 $\times$ (380 cm$^{-3}$/n$_e$). This may be in tension with the $\eta$ $\gtrsim$ 1 required by cosmological simulations, unless a significant fraction of the outflowing mass is in other gas phases and has sufficient velocity to escape the galaxy halos.

The abundance and demographics of dark matter substructure is important for many areas in astrophysics and cosmological $N$-body simulations have been the primary tool used to investigate them. However, it has recently become clear that the simulations are subject to numerical artefacts, which hampers a proper treatment of the tidal evolution of subhaloes. Unfortunately, no analytical models that accurately describe subhalo evolution exist either. We therefore present a library of idealized, high resolution $N$-body simulations of the tidal evolution of individual subhaloes that can be used to calibrate semi-analytical models and to complement cosmological simulations. The simulations focus on minor mergers, i.e., the mass of the subhalo is much smaller than that of the host halo, such that the impact of dynamical friction is negligible. This setup allows the adoption of a fixed analytical potential for modelling the host halo. The dynamical evolution of subhaloes is followed with $N$-body computations. In the library, four parameters, two of which characterize the subhalo orbit with respect to the host halo, and the two concentrations of the host- and subhalo, are varied over the ranges encountered in cosmological simulations. We show several representative examples from the library that illustrate the evolution of the subhalo mass and velocity dispersion profiles. Additionally, we make publicly available a pre-trained non-parametric model of the subhalo mass evolution based on random forest regression. This model is able to interpolate the simulation data at the 0.1\,dex level and provides efficient access to the data for further use in modelling.

Following the success of the Plateau de Bure high-z Blue Sequence Survey (PHIBSS), we present the PHIBSS2 legacy program, a survey of the molecular gas properties of star-forming galaxies on and around the star formation main sequence (MS) at different redshifts using NOEMA. This survey significantly extends the existing sample of star-forming galaxies with CO molecular gas measurements, probing the peak epoch of star formation (z=1-1.6) as well as its building-up (z=2-3) and winding-down (z=0.5-0.8) phases. The targets are drawn from the GOODS, COSMOS, and AEGIS deep fields and uniformly sample the MS in the stellar mass (M*) - star formation rate (SFR) plane with log(M*/Msun) = 10-11.8. We describe the survey strategy and sample selection before focusing on the results obtained at z=0.5-0.8, where we report 60 CO(2-1) detections out of 61 targets. We determine their molecular gas masses and separately obtain disc sizes and bulge-to-total (B/T) luminosity ratios from HST I-band images. The median molecular gas-to-stellar mass ratio, gas fraction, and depletion time as well as their dependence with M* and offset from the MS follow published scaling relations for a much larger sample of galaxies spanning a wider range of redshifts. The galaxy-averaged Kennicutt-Schmidt (KS) relation between molecular gas and SFR surface densities is strikingly linear, pointing towards similar star formation timescales within galaxies at any given epoch. In terms of morphology, the molecular gas content, the SFR, the disc stellar mass, and the disc molecular gas fraction do not seem to correlate with B/T and the stellar surface density, which suggests an ongoing supply of fresh molecular gas to compensate for the build-up of the bulge. Our measurements do not yield any significant variation of the depletion time with B/T and hence no strong evidence for morphological quenching within the scatter of the MS.

We investigate the future evolution of the Smith Cloud by performing hydrodynamical simulations of the cloud impact onto the gaseous Milky Way Galactic disk. We assume a local origin for the cloud and thus do not include a dark matter component to stabilize it. Our main focus is the cloud's influence on the local and global star formation rate (SFR) of the Galaxy and whether or not it leads to an observable event in the far future. Our model assumes two extremes for the mass of the Smith Cloud, an upper mass limit of 10$^7$ M$_{\odot}$ and a lower mass limit of 10$^6$ M$_{\odot}$, compared to the observational value of a few 10$^6$ M$_{\odot}$. In addition, we also make the conservative assumption that the entirety of the cloud mass of the extended Smith Cloud is concentrated within the tip of the cloud. We find that the impact of the low-mass cloud produces no noticeable change in neither the global SFR nor the local SFR at the cloud impact site within the galactic disk. For the high-mass cloud we find a short-term (roughly 5 Myr) increase of the global SFR of up to 1 M$_{\odot}$ yr$^{-1}$, which nearly doubles the normal Milky Way SFR. This highly localized starburst should be observable.

The Galactic Center black hole Sgr A* is the archetypical example of an underfed massive black hole. The extremely low accretion rate can be understood in radiatively inefficient accretion flow models. Testing those models has proven to be difficult due to the lack of suitable probes. Radio and submm polarization measurements constrain the flow very close to the event horizon. X-ray observations resolving the Bondi radius yield an estimate roughly four orders of magnitude further out. Here, we present a new, indirect measurement of the accretion flow density at intermediate radii. We use the dynamics of the gas cloud G2 to probe the ambient density. We detect the presence of a drag force slowing down G2 with a statistical significance of approx 9 {\sigma}. This probes the accretion flow density at around 1000 Schwarzschild radii and yields a number density of approx. 4 x 10^3 cm^-3. Self-similar accretion models where the density follows a power law radial profile between the inner zone and the Bondi radius have predicted similar values.

State-of-the-art integral field surveys like $\mathrm{ATLAS^{3D}}$, SLUGGS, CALIFA, SAMI, and MaNGA provide large data sets of kinematical observations of early-type galaxies (ETGs), yielding constraints on the formation of ETGs. Using the cosmological hydrodynamical \textit{Magneticum Pathfinder} simulations, we investigate the paradigm of fast and slow rotating ETGs in a fully cosmological context. We show that the ETGs within the \textit{Magneticum} simulation are in remarkable agreement with the observations, revealing fast and slow rotators quantified by the angular momentum proxy $\lambda_{\mathrm{R}}$ and the flattening $\epsilon$ with the observed prevalence. Taking full advantage of the three-dimensional data, we demonstrate that the dichotomy between fast and slow rotating galaxies gets enhanced, showing an upper and lower population separated by an underpopulated region in the edge-on $\lambda_{\mathrm{R}}$-$\epsilon$ plane. Following the evolution of the $\lambda_{\mathrm{R}}$-$\epsilon$ plane through cosmic time, we find that, while the upper population is already in place at $z=2$, the lower population gets statistically significant below $z=1$ with a gradual increase. At least $50\%$ of the galaxies transition from fast to slow rotators on a short timescale, in most cases associated to a significant merger event. Furthermore, we connect the $M_{*}$-$j_{*}$ plane, quantified by the $b$-value, with the $\lambda_{\mathrm{R}}$-$\epsilon$ plane, revealing a strong correlation between the position of a galaxy in the $\lambda_{\mathrm{R}}$-$\epsilon$ plane and the $b$-value. Going one step further, we classify our sample based on features in their velocity map, finding all five common kinematic groups, also including the recently observed group of prolate rotators, populating distinct regions in the $\lambda_{\mathrm{R}}$-$b$ plane.

We use the $\mathit{Gaia}$ DR2 distances of about 700 mid-infrared selected young stellar objects in the benchmark giant molecular cloud Orion A to infer its 3D shape and orientation. We find that Orion A is not the fairly straight filamentary cloud that we see in (2D) projection, but instead a cometary-like cloud oriented toward the Galactic plane, with two distinct components: a denser and enhanced star-forming (bent) Head, and a lower density and star-formation quieter $\sim$75 pc long Tail. The true extent of Orion A is not the projected $\sim$40 pc but $\sim$90 pc, making it by far the largest molecular cloud in the local neighborhood. Its aspect ratio ($\sim$30:1) and high column-density fraction ($\sim45\%$) make it similar to large-scale Milky Way filaments ("bones"), despite its distance to the galactic mid-plane being an order of magnitude larger than typically found for these structures.

We exploited observational constraints on stars, gas and nucleosynthesis ashes for the closest region of recent massive-star formation, Scorpius-Centaurus OB2, and combined them with 3D hydrodynamical simulations, in order to address physics and history for the case of the Scorpius-Centaurus superbubble. We used published cold gas observations through PLANCK survey data processing, HERSCHEL and APEX, continuum and molecular line observations. We analysed the Galactic All Sky Survey (GASS) to investigate shell structures in atomic hydrogen, and used HIPPARCOS and Gaia data in combination with interstellar absorption against stars to obtain new constraints for the distance to the Hi features. Hot gas is traced in soft X-rays via the ROSAT all sky survey. Nucleosynthesis ejecta from massive stars were traced with new INTEGRAL spectrometer observations via 26Al radioactivity. We also performed 3D hydrodynamical simulations for the Sco-Cen superbubble. To investigate the impact of massive star feedback on extended clouds, we simulate the interaction of a turbulent cloud with the hot, pressurised gas in a superbubble. The hot gas fills the tenuous regions of the cloud and compresses the denser parts. Stars formed in these dense clumps would have distinct spatial and kinematic distributions. The combined results from observations and simulations are consistent with a scenario where dense gas was initially distributed in a band elongated in the direction now occupied by the OB association. Superbubbles powered by massive stars would then repeatedly break out of the elongated parent cloud, surround and squash the denser parts of the gas sheet and thus induce more star formation. The expected spatial and kinematic distribution of stars is consistent with observations of Sco-Cen. The scenario might apply to many similar regions in the Galaxy and also to AGN-related superbubbles.

Observations of prestellar cores in star-forming filaments show two distinct morphologies. While molecular line measurements often show broad cores, submillimeter continuum observations predominantly display pinched cores compared to the bulk of the filament gas. In order to explain how different morphologies arise, we use the gravitational instability model where prestellar cores form by growing density perturbations. The radial extent at each position is set by the local line-mass. We show that the ratio of core radius to filament radius is determined by the initial line-mass of the filament. Additionally, the core morphology is independent of perturbation length scale and inclination, which makes it an ideal diagnostic for observations. Filaments with a line-mass of less than half its critical value should form broad cores, whereas filaments with more than half its critical line-mass value should form pinched cores. For filaments embedded in a constant background pressure, the dominant perturbation growth times significantly differ for low and high line-mass filaments. Therefore, we predict that only one population of cores is present if all filaments within a region begin with similar initial perturbations.

We study accretion driven turbulence for different inflow velocities in star forming filaments using the code ramses. Filaments are rarely isolated objects and their gravitational potential will lead to radially dominated accretion. In the non-gravitational case, accretion by itself can already provoke non-isotropic, radially dominated turbulent motions responsible for the complex structure and non-thermal line widths observed in filaments. We find that there is a direct linear relation between the absolute value of the total density weighted velocity dispersion and the infall velocity. The turbulent velocity dispersion in the filaments is independent of sound speed or any net flow along the filament. We show that the density weighted velocity dispersion acts as an additional pressure term supporting the filament in hydrostatic equilibrium. Comparing to observations, we find that the projected non-thermal line width variation is generally subsonic independent of inflow velocity.

The so-called S2 star reached its closest approach to the massive black hole (BH) at around 1500 $R_\mathrm{s}$ in May 2018. It has been proposed that the interaction of its stellar wind with the high-density accretion flow at this distance from Sgr A* will lead to a detectable, month-long X-ray flare. Our goal is to verify whether or not the S2 star wind can be used as a diagnostic tool to infer the properties of the accretion flow towards Sgr A* at its pericentre (an unprobed distance regime), putting important constraints on BH accretion flow models. We run a series of three-dimensional adaptive mesh refinement simulations with the help of the Ramses code which include the realistic treatment of the interaction of S2's stellar wind with the accretion flow along its orbit and - apart from hydrodynamical and thermodynamical effects - include the tidal interaction with the massive BH. These are post-processed to derive the X-ray emission in the observable 2-10 keV window. No significant excess of X-ray emission from Sgr A* is found for typical accretion flow models. A measurable excess is produced for a significantly increased density of the accretion flow. This can, however, be ruled out for standard power-law accretion flow models as in this case the thermal X-ray emission without the S2 wind interaction would already exceed the observed quiescent luminosity. Only a significant change of the wind parameters (increased mass loss rate and decreased wind velocity) might lead to an (marginally) observable X-ray flaring event. Even the detection of an (month-long) X-ray flare during the pericentre passage of the S2 star would not allow for strict constraints to be put on the accretion flow around Sgr A* due to the degeneracy caused by the dependence on multiple parameters (of the accretion flow model as well as the stellar wind).

Nuclear outflows driven by accreting massive black holes are one of the main feedback mechanisms invoked at high-z to reproduce the distinct separation between star-forming, disk galaxies and quiescent spheroidal systems. Yet, our knowledge of feedback at high-z remains limited by the lack of observations of the multiple gas phases in galaxy outflows. In this work we use new deep, high-spatial resolution ALMA CO(3-2) and archival VLT/SINFONI H$\alpha$ observations to study the molecular and ionized components of the AGN-driven outflow in zC400528 ---a massive, main sequence galaxy at z=2.3 in the process of quenching. We detect a powerful molecular outflow that shows a positive velocity gradient and extends for at least ~10 kpc from the nuclear region, about three times the projected size of the ionized wind. The molecular gas in the outflow does not reach velocities high enough to escape the galaxy and is therefore expected to be reaccreted. Keeping in mind the various assumptions involved in the analysis, we find that the mass and energetics of the outflow are dominated by the molecular phase. The AGN-driven outflow in zC400528 is powerful enough to deplete the molecular gas reservoir on a timescale at least twice shorter than that needed to exhaust it by star formation. This suggests that the nuclear outflow is one of the main quenching engines at work in the observed suppression of the central star-formation activity in zC400528.

In the hierarchical two-phase formation scenario, the extended halos of early type galaxies (ETGs) are expected to have different physical properties from those of the galaxies' central regions. This work aims at characterizing the kinematic properties of ETG halos using planetary nebulae (PNe) as tracers, which allow us to overcome the limitations of absorption line spectroscopy of continuum at low surface brightness. We present two-dimensional velocity and velocity dispersion fields for 33 ETGs, including both fast (FRs) and slow rotators (SRs), making this the largest kinematic survey to-date of extragalactic PNe. The velocity fields extend out to a median 5.6 effective radii (Re), with a range [3Re-13Re] for the PN.S ETGs. We complemented the PN kinematics with absorption line data from the literature. We find that ETGs typically show a kinematic transition between inner regions and halo. Estimated transition radii in units of Re anticorrelate with stellar mass. SRs have increased but still modest rotational support at large radii, while most of the FRs show a decrease in rotation, due to the fading of the stellar disk in the outer, more slowly rotating spheroid. 30% of the FRs are dominated by rotation also at large radii. Most ETGs have flat or slightly falling halo velocity dispersion profiles, but 15% of the sample have steeply falling profiles. 40% of the FRs show kinematic twists, misalignments, or rotation along two axes, indicating that they turn from oblate at the center to triaxial in the halo, consistently with the distribution of isophote twists from extended photometry. ETGs have more diverse kinematic properties in their halos than in the central regions. The observed kinematic transition to the halo and its dependence on stellar mass is consistent with LambdaCDM simulations and supports a two-phase formation scenario. [abridged]

The formation of compact dusty clouds in diffuse interstellar medium (ISM) has been recently proposed and studied by Tsytovich et al. (2014). In the present paper, an effect of the clouds on the rate of H$\to$H$_2$ transition in the ISM is examined. We discuss the mechanisms leading to the formation of the clouds and the creation of gaseous clumps around them, analyze the main processes determining the efficiency of the recombination of atomic hydrogen on dust grains, and estimate the expected enhancement of the global H$_2$ formation due to the presence of the clouds. In conclusion, we argue that the compact dusty clouds provide optimum conditions for the atomic recombination in diffuse ISM, and point out some astrophysical implications of the resulting H$_2$ formation enhancement.

The source G2 has already completed its pericentre passage around Sgr A*, the super-massive black hole in the centre of our Galaxy. Although it has been monitored for 15 years, its astrophysical nature and origin still remain unknown. In this work, we aim to test the hypothesis of G2 being the result of a stellar wind collision. To do so, we study the motion and final fate of gas clumps formed as a result of collisions of stellar winds in massive binaries. Our approach is based on a test-particle model in order to describe the trajectories of such clumps. The model takes into account the gravitational field of Sgr A*, the interaction of the clumps with the interstellar medium as well as their finite lifetimes. Our analysis allows us to reject the hypothesis based on four arguments: i) if G2 has followed a purely Keplerian orbit since its formation, it cannot have been produced in any of the known massive binaries since their motions are not consistent; ii) in general, gas clumps are evaporated through thermal conduction on very short timescale (< 100yr) before getting close enough to Sgr A*; iii) IRS 16SW, the best candidate for the origin of G2, cannot generate clumps as massive as G2; and iv) clumps ejected from IRS 16SW describe trajectories significantly different to the observed motion of G2.

We present the "SINS/zC-SINF AO survey" of 35 star-forming galaxies, the largest sample with deep adaptive optics-assisted (AO) near-infrared integral field spectroscopy at z~2. The observations, taken with SINFONI at the Very Large Telescope, resolve the Ha and [NII] line emission and kinematics on scales of ~1.5 kpc. In stellar mass, star formation rate, rest-optical colors and size, the AO sample is representative of its parent seeing-limited sample and probes the massive (M* ~ 2x10^9 - 3x10^11 Msun), actively star-forming (SFR ~ 10-600 Msun/yr) part of the z~2 galaxy population over a wide range in colors ((U-V)_rest ~ 0.15-1.5 mag) and half-light radii (R_e,H ~ 1-8.5 kpc). The sample overlaps largely with the "main sequence" of star-forming galaxies in the same redshift range to a similar K_AB = 23 magnitude limit; it has ~0.3 dex higher median specific SFR, ~0.1 mag bluer median (U-V)_rest color, and ~10% larger median rest-optical size. We describe the observations, data reduction, and extraction of basic flux and kinematic properties. With typically 3-4 times higher resolution and 4-5 times longer integrations (up to 23hr) than the seeing-limited datasets of the same objects, the AO data reveal much more detail in morphology and kinematics. The now complete AO observations confirm the majority of kinematically-classified disks and the typically elevated disk velocity dispersions previously reported based on subsets of the data. We derive typically flat or slightly negative radial [NII]/Ha gradients, with no significant trend with global galaxy properties, kinematic nature, or the presence of an AGN. Azimuthal variations in [NII]/Ha are seen in several sources and are associated with ionized gas outflows, and possible more metal-poor star-forming clumps or small companions. [Abridged]

Selecting disk galaxies from the cosmological, hydrodynamical simulation Magneticum Pathfinder we show that almost half of our poster child disk galaxies at $z=2$ show significantly declining rotation curves and low dark matter fractions, very similar to recently reported observations. These galaxies do not show any anomalous behavior, reside in standard dark matter halos and typically grow significantly in mass until $z = 0$, where they span all morphological classes, including disk galaxies matching present day rotation curves and observed dark matter fractions. Our findings demonstrate that declining rotation curves and low dark matter fractions in rotation dominated galaxies at $z=2$ appear naturally within the $\Lambda$CDM paradigm and reflect the complex baryonic physics, which plays a role at the peak epoch of star-formation. In addition, we find some dispersion dominated galaxies at $z=2$ which host a significant gas disk and exhibit similar shaped rotation curves as the disk galaxy population, rendering it difficult to differentiate between these two populations with currently available observation techniques.

This paper provides an update of our previous scaling relations (Genzel et al.2015) between galaxy integrated molecular gas masses, stellar masses and star formation rates, in the framework of the star formation main-sequence (MS), with the main goal to test for possible systematic effects. For this purpose our new study combines three independent methods of determining molecular gas masses from CO line fluxes, far-infrared dust spectral energy distributions, and ~1mm dust photometry, in a large sample of 1444 star forming galaxies (SFGs) between z=0 and 4. The sample covers the stellar mass range log(M*/M_solar)=9.0-11.8, and star formation rates relative to that on the MS, delta_MS=SFR/SFR(MS), from 10^{-1.3} to 10^{2.2}. Our most important finding is that all data sets, despite the different techniques and analysis methods used, follow the same scaling trends, once method-to-method zero point offsets are minimized and uncertainties are properly taken into account. The molecular gas depletion time t_depl, defined as the ratio of molecular gas mass to star formation rate, scales as (1+z)^{-0.6}x(delta_MS)^{-0.44}, and is only weakly dependent on stellar mass. The ratio of molecular-to-stellar mass mu_gas depends on (1+z)^{2.5}x (delta_MS)^{0.52}x(M*)^{-0.36}, which tracks the evolution of the specific star formation rate. The redshift dependence of mu_gas requires a curvature term, as may the mass-dependences of t_depl and mu_gas. We find no or only weak correlations of t_depl and mu_gas with optical size R or surface density once one removes the above scalings, but we caution that optical sizes may not be appropriate for the high gas and dust columns at high-z.

The interplay between gravity, turbulence and the magnetic field determines the evolution of the molecular ISM and the formation of the stars. In spite of growing interests, there remains a lack of understanding of the importance of magnetic field over multiple scales. We derive the magnetic energy spectrum -- a measure that constraints the multi-scale distribution of the magnetic energy, and compare it with the gravitational energy spectrum derived in Li \& Burkert (2016). In our formalism, the gravitational energy spectrum is purely determined by the surface density PDF, and the magnetic energy spectrum is determined by both the surface density PDF and the magnetic-field-density relation. If regions have density PDFs close to $P(\Sigma)\sim \Sigma^{-2}$ and a universal magnetic field-density relation $B\sim \rho^{1/2}$, we expect a multi-scale near equipartition between gravity and the magnetic fields. This equipartition is found to be true in NGC6334 where estimates of magnetic fields over multiple scales (from 0.1 pc to a few parsec) are available. However, the current observations are still limited in sample size. In the future, it is necessary to obtain multi-scale measurements of magnetic fields from different clouds with different surface density PDFs and apply our formalism to further study the gravity-magnetic field interplay.

Accurately predicting the demographics of dark matter (DM) substructure is of paramount importance for many fields of astrophysics, including gravitational lensing, galaxy evolution, halo occupation modeling, and constraining the nature of dark matter. Because of its strongly non-linear nature, DM substructure is typically modeled using N-body simulations, which reveal that large fractions of DM subhaloes undergo complete disruption. In this paper we use both analytical estimates and idealized numerical simulations to investigate whether this disruption is mainly physical, due to tidal heating and stripping, or numerical (i.e., artificial). We show that, contrary to naive expectation, subhaloes that experience a tidal shock $\Delta E$ that exceeds the subhalo's binding energy, $|E_{\rm b}|$, do not undergo disruption, even when $\Delta E/|E_{\rm b}|$ is as large as 100. Along the same line, and contrary to existing claims in the literature, instantaneously stripping matter from the outskirts of a DM subhalo also does not result in its complete disruption, even when the instantaneous remnant has positive binding energy. In addition, we show that tidal heating due to high-speed (impulsive) encounters with other subhaloes (`harassment'), is negligible compared to the tidal effects due to the host halo. Hence, we conclude that, in the absence of baryonic processes, the complete, physical disruption of CDM substructure is extremely rare, and that most disruption in numerical simulations therefore must be artificial. We discuss various processes that have been associated with numerical overmerging, and conclude that inadequate force-softening is the most likely culprit.

The star-forming ability of a molecular cloud depends on the fraction of gas it can cycle into the dense-phase. Consequently, one of the crucial questions in reconciling star-formation in clouds is to understand the factors that control this process. While it is widely accepted that the variation in ambient conditions can alter significantly the ability of a cloud to spawn stars, the observed variation in the star-formation rate in nearby clouds that experience similar ambient conditions, presents an interesting question. In this work we attempted to reconcile this variation within the paradigm of colliding flows. To this end we develop self-gravitating, hydrodynamic realisations of identical flows, but allowed to collide off-centre. Typical observational diagnostics such as the gas-velocity dispersion, the fraction of dense-gas, the column density distribution ({\small N-PDF}), the distribution of gas mass as a function of $K$-band extinction and the strength of compressional/solenoidal modes in the post-collision cloud were deduced for different choices of the impact parameter of collision. We find that a strongly sheared cloud is terribly inefficient in cycling gas into the dense phase and that such a cloud can possibly reconcile the sluggish nature of star-formation reported for some clouds. Within the paradigm of cloud-formation via colliding flows this is possible in case of flows colliding with a relatively large impact parameter. We conclude that compressional modes - though probably essential - are insufficient to ensure a relatively higher star-formation efficiency in a cloud.

[Abridged] The Lupus I cloud is found between the Upper-Scorpius and the Upper-Centaurus-Lupus sub-groups, where the expanding USco HI shell appears to interact with a bubble currently driven by the winds of the remaining B-stars of UCL. We investigate if the Lupus I molecular could have formed in a colliding flow, and how the kinematics of the cloud might have been influenced by the larger scale gas dynamics. We performed APEX 13CO and C18O observations of three parts of Lupus. We compare these results to the atomic hydrogen data from the GASS HI survey and our dust emission results presented in the previous paper. Based on the velocity information, we present a geometric model for the interaction zone between the USco shell and the UCL wind bubble. We present evidence that the molecular gas of Lupus I is tightly linked to the atomic material of the USco shell. The CO emission in Lupus I is found mainly at velocities in the same range as the HI velocities. Thus, the molecular cloud is co-moving with the expanding USco atomic Hi shell. The gas in the cloud shows a complex kinematic structure with several line-of-sight components that overlay each other. The non-thermal velocity dispersion is in the transonic regime in all parts of the cloud and could be injected by external compression. Our observations and the derived geometric model agree with a scenario where Lupus I is located in the interaction zone between the USco shell and the UCL wind bubble. The kinematics observations are consistent with a scenario where the Lupus I cloud formed via shell instabilities. The particular location of Lupus I between USco and UCL suggests that counter-pressure from the UCL wind bubble and pre-existing density enhancements, perhaps left over from the gas stream that formed the stellar subgroups, may have played a role in its formation.

The bathtub model of the star forming interstellar medium is based on the powerful constraint that mass has to be conserved when gas flows through its various thermal and density phases, ending up eventually in a young star or being blown away by stellar feedback. It predicts that the star formation rate of a molecular cloud is not determined by the cloud's mass or its internal collapse timescale, but rather by the accretion rate of new gas. For the most simple case of a constant accretion flow an equilibrium state is reached quickly where the star formation rate equals the accretion rate and where the dense gas mass is constant and independent of time. The mass of the young star cluster, on the other hand, increases linearly with time. The stellar mass fraction therefore represents a sensitive clock to measure the age of the star-forming region. The bathtub model predicts that the efficiency of star formation is small, of order 1%, even in the dense filamentary phases of molecular clouds. It provides a simple explanation for the dense gas fraction of order 10% in molecular clouds and for the large gas depletion timescales of star-forming galaxies of order $5 \times 10^8 - 10^9$ yrs.

High-resolution observations from the sub-mm to the optical wavelength regime resolve the central few 100pc region of nearby galaxies in great detail. They reveal a large diversity of features: thick gas and stellar discs, nuclear starbursts, in- and outflows, central activity, jet interaction, etc. Concentrating on the role circumnuclear discs play in the life cycles of galactic nuclei, we employ 3D adaptive mesh refinement hydrodynamical simulations with the RAMSES code to self-consistently trace the evolution from a quasi-stable gas disc, undergoing gravitational (Toomre) instability, the formation of clumps and stars and the disc's subsequent, partial dispersal via stellar feedback. Our approach builds upon the observational finding that many nearby Seyfert galaxies have undergone intense nuclear starbursts in their recent past and in many nearby sources star formation is concentrated in a handful of clumps on a few 100pc distant from the galactic centre. We show that such observations can be understood as the result of gravitational instabilities in dense circumnuclear discs. By comparing these simulations to available integral field unit observations of a sample of nearby galactic nuclei, we find consistent gas and stellar masses, kinematics, star formation and outflow properties. Important ingredients in the simulations are the self-consistent treatment of star formation and the dynamical evolution of the stellar distribution as well as the modelling of a delay time distribution for the supernova feedback. The knowledge of the resulting simulated density structure and kinematics on pc scale is vital for understanding inflow and feedback processes towards galactic scales.

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.

A large fraction of the gas in the Galaxy is cold, dense, and molecular. If all this gas collapsed under the influence of gravity and formed stars in a local free-fall time, the star formation rate in the Galaxy would exceed that observed by more than an order of magnitude. Other star-forming galaxies behave similarly. Yet observations and simulations both suggest that the molecular gas is indeed gravitationally collapsing, albeit hierarchically. Prompt stellar feedback offers a potential solution to the low observed star formation rate if it quickly disrupts star-forming clouds during gravitational collapse. However, this requires that molecular clouds must be short-lived objects, raising the question of how so much gas can be observed in the molecular phase. This can occur only if molecular clouds form as quickly as they are destroyed, maintaining a global equilibrium fraction of dense gas. We therefore examine cloud formation timescales. We first demonstrate that supernova and superbubble sweeping cannot produce dense gas at the rate required to match the cloud destruction rate. On the other hand, Toomre gravitational instability can reach the required production rate. We thus argue that, although dense, star-forming gas may last only around a single global free-fall time, the dense gas in star-forming galaxies can globally exist in a state of dynamic equilibrium between formation by gravitational instability, and disruption by stellar feedback. At redshift z >~ 2, the Toomre instability timescale decreases, resulting in a prediction of higher molecular gas fractions at early times, in agreement with observations.

The nature of the gaseous and dusty cloud G2 in the Galactic Centre is still under debate. We present three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of G2, modeled as an outflow from a "compact source" moving on the observed orbit. The construction of mock position-velocity (PV) diagrams enables a direct comparison with observations and allow us to conclude that the observational properties of the gaseous component of G2 could be matched by a massive ($\dot{M}_\mathrm{w}=5\times 10^{-7} \;M_{\odot} \mathrm{yr^{-1}}$) and slow ($50 \;\mathrm{km \;s^{-1}}$) outflow, as observed for T Tauri stars. In order for this to be true, only the material at larger ($>100 \;\mathrm{AU}$) distances from the source must be actually emitting, otherwise G2 would appear too compact compared to the observed PV diagrams. On the other hand, the presence of a central dusty source might be able to explain the compactness of G2's dust component. In the present scenario, 5-10 years after pericentre the compact source should decouple from the previously ejected material, due to the hydrodynamic interaction of the latter with the surrounding hot and dense atmosphere. In this case, a new outflow should form, ahead of the previous one, which would be the smoking gun evidence for an outflow scenario.

We investigate the stellar mass and baryonic mass Tully-Fisher relations (TFRs) of massive star-forming disk galaxies at redshift z~2.3 and z~0.9 as part of the KMOS^3D integral field spectroscopy survey. Our spatially resolved data allow reliable modelling of individual galaxies, including the effect of pressure support on the inferred gravitational potential. At fixed circular velocity, we find higher baryonic masses and similar stellar masses at z~2.3 as compared to z~0.9. Together with the decreasing gas-to-stellar mass ratios with decreasing redshift, this implies that the contribution of dark matter to the dynamical mass at the galaxy scale increases towards lower redshift. A comparison to local relations reveals a negative evolution of the stellar and baryonic TFR zero-points from z=0 to z~0.9, no evolution of the stellar TFR zero-point from z~0.9 to z~2.3, and a positive evolution of the baryonic TFR zero-point from z~0.9 to z~2.3. We discuss a toy model of disk galaxy evolution to explain the observed, non-monotonic TFR evolution, taking into account the empirically motivated redshift dependencies of galactic gas fractions, and of the relative amount of baryons to dark matter on galaxy and halo scales.

We explore the H-alpha emission in the massive quiescent galaxies observed by the KMOS-3D survey at 0.7 < z < 2.7. The H-alpha line is robustly detected in 20 out of 120 UVJ-selected quiescent galaxies, and we classify the emission mechanism using the H-alpha line width and the [NII]/H-alpha line ratio. We find that AGN are likely to be responsible for the line emission in more than half of the cases. We also find robust evidence for star formation activity in nine quiescent galaxies, which we explore in detail. The H-alpha kinematics reveal rotating disks in five of the nine galaxies. The dust-corrected H-alpha star formation rates are low (0.2 - 7 Msun/yr), and place these systems significantly below the main sequence. The 24micron-based infrared luminosities, instead, overestimate the star formation rates. These galaxies present a lower gas-phase metallicity compared to star-forming objects with similar stellar mass, and many of them have close companions. We therefore conclude that the low-level star formation activity in these nine quiescent galaxies is likely to be fueled by inflowing gas or minor mergers, and could be a sign of rejuvenation events.

We present a solution for the observed core fragmentation of filaments in the Taurus L1517 dark cloud which previously could not be explained (Hacar et. al 2011). Core fragmentation is a vital step for the formation of stars. Observations suggest a connection to the filamentary structure of the cloud gas but it remains unclear which process is responsible. We show that the gravitational instability process of an isothermal cylinder can account for the exhibited fragmentation under the assumption that the perturbation grows on the dominant wavelength. We use numerical simulations with the code RAMSES, estimate observed column densities and line-of-sight velocities and compare them to the observations. A critical factor for the observed fragmentation is that cores grow by redistributing mass within the filament and thus the density between the cores decreases over the fragmentation process. This often leads to wrong dominant wavelength estimates as it is strongly dependent on the initial central density. We argue that non-linear effects also play an important role on the evolution of the fragmentation. Once the density perturbation grows above the critical line-mass, non-linearity leads to an enhancement of the central core density in comparison to the analytical prediction. Choosing the correct initial conditions with perturbation strengths of around 20\%, leads to inclination corrected line-of-sight velocities and central core densities within the observational measurement error in a realistic evolution time.

In early 2014 the fast-moving near-infrared source G2 reached its closest approach to the supermassive black hole Sgr A* in the Galactic Center. We report on the evolution of the ionized gaseous component and the dusty component of G2 immediately after this event, revealed by new observations obtained in 2015 and 2016 with the SINFONI integral field spectrograph and the NACO imager at the ESO VLT. The spatially resolved dynamics of the Br$\gamma$ line emission can be accounted for by the ballistic motion and tidal shearing of a test-particle cloud that has followed a highly eccentric Keplerian orbit around the black hole for the last 12 years. The non-detection of a drag force or any strong hydrodynamic interaction with the hot gas in the inner accretion zone limits the ambient density to less than a few 10$^3$ cm$^{-3}$ at the distance of closest approach (1500 $R_s$), assuming G2 is a spherical cloud moving through a stationary and homogeneous atmosphere. The dust continuum emission is unresolved in L'-band, but stays consistent with the location of the Br$\gamma$ emission. The total luminosity of the Br$\gamma$ and L' emission has remained constant to within the measurement uncertainty. The nature and origin of G2 are likely related to that of the precursor source G1, since their orbital evolution is similar, though not identical. Both object are also likely related to a trailing tail structure, which is continuously connected to G2 over a large range in position and radial velocity.

We exploit the deep resolved Halpha kinematic data from the KMOS^3D and SINS/zC-SINF surveys to examine the largely unexplored outer disk kinematics of star-forming galaxies (SFGs) out to the peak of cosmic star formation. Our sample contains 101 SFGs representative of the more massive (9.3 < log(M*/Msun) < 11.5) main sequence population at 0.6<z<2.6. Through a novel stacking approach we are able to constrain a representative rotation curve extending out to ~4 effective radii. This average rotation curve exhibits a significant drop in rotation velocity beyond the turnover, with a slope of Delta(V)/Delta(R) = $-0.26^{+0.10}_{-0.09}$ in units of normalized coordinates V/V_max and R/R_turn. This result confirms that the fall-off seen previously in some individual galaxies is a common feature of our sample of high-z disks. We show that this outer fall-off strikingly deviates from the flat or mildly rising rotation curves of local spiral galaxies of similar masses. We furthermore compare our data with models including baryons and dark matter demonstrating that the falling stacked rotation curve can be explained by a high mass fraction of baryons relative to the total dark matter halo (m_d>~0.05) in combination with a sizeable level of pressure support in the outer disk. These findings are in agreement with recent studies demonstrating that star-forming disks at high redshift are strongly baryon dominated within the disk scale, and furthermore suggest that pressure gradients caused by large turbulent gas motions are present even in their outer disks. We demonstrate that these results are largely independent of our model assumptions such as the presence of a central stellar bulge, the effect of adiabatic contraction at fixed m_d, and variations in the concentration parameter.

The geometry and intrinsic ellipticity distribution of ultra diffuse galaxies (UDGs) is determined from the line-of-sight distribution of axial ratios q of a large sample of UDGs, detected by Koda et al. (2015) in the Coma cluster. With high significance the data rules out an oblate, disk-like geometry, characterised by major axi a=b>c. The data is however in good agreement with prolate shapes, corresponding to a=b<c. This indicates that UDGs are not thickened, rotating, axisymmetric disks, puffed up by violent processes. Instead they are anisotropic elongated cigar- or bar-like structures, similar to the prolate dwarf spheroidal galaxy population of the Local Group. The intrinsic distribution of axial ratios of the Coma UDGs is flat in the range of 0.4 <= a/c <= 0.9 which will provide important constraints for theoretical models of their origin. Formation scenarios that could explain the extended prolate nature of UDGs are discussed.

In cold dark matter cosmology, the baryonic components of galaxies are thought to be mixed with and embedded in non-baryonic and non-relativistic dark matter, which dominates the total mass of the galaxy and its dark matter halo. In the local Universe, the mass of dark matter within a galactic disk increases with disk radius, becoming appreciable and then dominant in the outer, baryonic regions of the disks of star-forming galaxies. This results in rotation velocities of the visible matter within the disk that are constant or increasing with disk radius. Comparison between the dynamical mass and the sum of stellar and cold gas mass at the peak epoch of galaxy formation, inferred from ancillary data, suggest high baryon factions in the inner, star-forming regions of the disks. Although this implied baryon fraction may be larger than in the local Universe, the systematic uncertainties (stellar initial mass function, calibration of gas masses) render such comparisons inconclusive in terms of the mass of dark matter. Here we report rotation curves for the outer disks of six massive star-forming galaxies, and find that the rotation velocities are not constant, but decrease with radius. We propose that this trend arises because of two main factors: first, a large fraction of the massive, high-redshift galaxy population was strongly baryon dominated, with dark matter playing a smaller part than in the local Universe; and second, the large velocity dispersion in high-redshift disks introduces a substantial pressure term that leads to a decrease in rotation velocity with increasing radius. The effect of both factors appears to increase with redshift. Qualitatively, the observations suggest that baryons in the early Universe efficiently condensed at the centres of dark matter halos when gas fractions were high, and dark matter was less concentrated. [Abridged]

In the past years several authors studied the abundance of satellites around galaxies in order to better estimate the halo masses of host galaxies. To investigate this connection, we analyze galaxies with $M_\mathrm{star}\geq\,10^{10}\,M_{\odot}$ from the hydrodynamical cosmological simulation Magneticum. We find that the satellite fraction of centrals is independent of their morphology. With the exception of very massive galaxies at low redshift, our results do not support the assumption that the dark matter (DM) haloes of spheroidal galaxies are significantly more massive than those of disc galaxies at fixed $M_\mathrm{star}$. We show that the density-morphology-relation starts to build up at $z\sim2$ and is independent of the star-formation properties of central galaxies. We conclude that environmental quenching is more important for satellites than for centrals. Our simulations indicate that conformity is already in place at $z=2$, where formation redshift and current star-formation rate (SFR) of central and satellite galaxies correlate. Centrals with low SFRs have formed earlier (at fixed $M_\mathrm{star}$) while centrals with high SFR formed later, with typical formation redshifts well in agreement with observations. However, we confirm the recent observations that the apparent number of satellites of spheroidal galaxies is significantly larger than for disc galaxies. This difference completely originates from the inclusion of companion galaxies, i.e. galaxies that do not sit in the potential minimum of a DM halo. Thus, due to the density-morphological-relation the number of satellites is not a good tracer for the halo mass, unless samples are restricted to the central galaxies of DM haloes.

We discuss the possibility that gravitational focusing, is responsible for the power-law mass function of star clusters $N(\log M) \propto M^{-1}$. This power law can be produced asymptotically when the mass accretion rate of an object depends upon the mass of the accreting body as $\dot{M} \propto M^2$. While Bondi-Hoyle-Littleton accretion formally produces this dependence on mass in a uniform medium, realistic environments are much more complicated. However, numerical simulations in SPH allowing for sink formation yield such an asymptotic power-law mass function. We perform pure N-body simulations to isolate the effects of gravity from those of gas physics and to show that clusters naturally result with the power-law mass distribution. We also consider the physical conditions necessary to produce clusters on appropriate timescales. Our results help support the idea that gravitationally-dominated accretion is the most likely mechanism for producing the cluster mass function.

Recently reported variations in the typical physical properties of Galactic and extra-Galactic molecular clouds (MCs), and in their ability to form stars have been attributed to local variations in the magnitude of interstellar pressure. Inferences from these surveys have called into question two long-standing beliefs that the MCs : 1 are Virialised entities and (2) have approximately constant surface density i.e., the validity of the Larson's third law. In this work we invoke the framework of cloud-formation via collisions between warm gas flows. Post-collision clouds forming in these realisations cool rapidly and evolve primarily via the interplay between the Non-linear Thin Shell Instability (NTSI), and the self-gravity. Over the course of these simulations we traced the temporal evolution of the surface density of the assembled clouds, the fraction of dense gas, the distribution of gas column density (NPDF), and the Virial nature of the assembled clouds. We conclude, these physical properties of MCs not only exhibit temporal variation, but their respective peak-magnitude also increases in proportion with the magnitude of external pressure, $P_{ext}$. The velocity dispersion in assembled clouds appears to follow the power-law, $\sigma_{gas}\propto P_{ext}^{0.23}$. Also, the power-law tail at higher densities becomes shallower with increasing magnitude of external pressure, for magnitudes, $P_{ext}/k_{B}\lesssim 10^{7}$ K cm$^{-3}$, at higher magnitudes such as those typically found in the Galactic CMZ ($P_{ext}/k_{B} > 10^{7}$ K cm$^{-3}$), the power-law shows significant steepening. Thus while our results are broadly consistent with inferences from various recent observational surveys, it appears, MCs hardly exhibit a unique set of properties, but rather a wide variety, that can be reconciled with a range of magnitudes of pressure between 10$^{4}$ K cm$^{-3}$ - 10$^{8}$ K cm$^{-3}$.

We study the stability of filaments in equilibrium between gravity and internal as well as external pressure using the grid based AMR-code RAMSES. A homogeneous, straight cylinder below a critical line mass is marginally stable. However, if the cylinder is bent, e.g. with a slight sinusoidal perturbation, an otherwise stable configuration starts to oscillate, is triggered into fragmentation and collapses. This previously unstudied behavior allows a filament to fragment at any given scale, as long as it has slight bends. We call this process `geometrical fragmentation'. In our realization the spacing between the cores matches the wavelength of the sinusoidal perturbation, whereas up to now, filaments were thought to be only fragmenting on the characteristical scale set by the mass-to-line ratio. Using first principles, we derive the oscillation period as well as the collapse timescale analytically. To enable a direct comparison with observations, we study the line-of-sight velocity for different inclinations. We show that the overall oscillation pattern can hide the infall signature of cores.

We present 0".2-resolution Atacama Large Millimeter/submillimeter Array observations at 870 um for 25 Halpha-seleced star-forming galaxies (SFGs) around the main-sequence at z=2.2-2.5. We detect significant 870 um continuum emission in 16 (64%) of these SFGs. The high-resolution maps reveal that the dust emission is mostly radiated from a single region close to the galaxy center. Exploiting the visibility data taken over a wide $uv$ distance range, we measure the half-light radii of the rest-frame far-infrared emission for the best sample of 12 massive galaxies with logM*>11. We find nine galaxies to be associated with extremely compact dust emission with R_{1/2,870um}<1.5 kpc, which is more than a factor of 2 smaller than their rest-optical sizes, R_{1/2,1.6um}=3.2 kpc, and is comparable with optical sizes of massive quiescent galaxies at similar redshifts. As they have an exponential disk with Sersic index of n=1.2 in the rest-optical, they are likely to be in the transition phase from extended disks to compact spheroids. Given their high star formation rate surface densities within the central 1 kpc of Sigma SFR1kpc=40 Msol/yr/kpc^2, the intense circumnuclear starbursts can rapidly build up a central bulge with Sigma M*1kpc>1e10 Msol/kpc^2 in several hundred Myr, i.e. by z~2. Moreover, ionized gas kinematics reveal that they are rotation-supported with an angular momentum as large as that of typical SFGs at z=1-3. Our results suggest bulges are commonly formed in extended rotating disks by internal processes, not involving major mergers.

We have used dedicated 0.7m telescopes in California and Israel to image the halos of ~200 galaxies in the Local Volume to 29 mag/sq arcsec, the sample mainly drawn from the 2MASS Large Galaxy Atlas (LGA). We supplement the LGA sample with dwarf galaxies and more distant giant ellipticals. Low surface brightness halos exceeding 50 kpc in diameter are found only in galaxies more luminous than L* and classic interaction signatures are relatively infrequent. Halo diameter is correlated with total galaxy luminosity. Extended low surface brightness halos are present even in galaxies as faint as M_V=-18. Edge-on galaxies with boxy bulges tend to lack extended spheroidal halos, while those with large classical bulges exhibit extended round halos, supporting the notions that boxy or barlike bulges originate from disks. Most face-on spiral galaxies present features that appear to be irregular extensions of spiral arms, although rare cases show smooth boundaries with no sign of star formation. Although we serendipitously discovered a dwarf galaxy undergoing tidal disruption in the halo of NGC 4449, we found no comparable examples in our general survey. A search for similar examples in the Local Volume identified hcc087, a tidally disrupting dwarf galaxy in the Hercules Cluster, but we do not confirm an anomalously large half-light radius reported for the dwarf VCC 1661.

Gravity plays a determining role in the evolution of the molecular ISM. In \citet{2016arXiv160304342L}, we proposed a measure called gravitational energy spectrum to quantify the importance of gravity on multiple physical scales. In this work, using a wavelet-based decomposition technique, we derive the gravitational energy spectra of the Orion A and the Orion B molecular cloud from observational data. The gravitational energy spectra exhibit power-law-like behaviours. From a few pc down to $\sim 0.1 $ pc scale, the Orion A and Orion B molecular cloud have $E_{\rm p}(k)\sim k^{-1.88}$ and $E_{\rm p}(k)\sim k^{-2.09}$, respectively. These scaling exponents are close to the scaling exponents of the kinetic energy power spectrum of compressible turbulence (where $E\sim k^{-2}$), with a near-equipartition of turbulent versus gravitational energy on multiple scales. This provides a clear evidence that gravity is able to counteract effectively against turbulent motion for these length scales. The results confirm our earlier analytical estimates. For the Orion A molecular cloud, gravity inevitably dominates over turbulence inside the cloud. Our results provide a clear observational proof that gravity is playing a determining role in the evolution these molecular clouds from the cloud scale down to $\sim 0.1\;\rm pc$. { However, turbulence is likely to dominate in clouds like California. } The method is general and should be applicable to all the astrophysical problems where gravity plays a role.

We study the evolution of G2 in a \textit{Compact Source Scenario}, where G2 is the outflow from a low-mass central star moving on the observed orbit. This is done through 3D AMR simulations of the hydrodynamic interaction of G2 with the surrounding hot accretion flow. A comparison with observations is done by means of mock position-velocity (PV) diagrams. We found that a massive ($\dot{M}_\mathrm{w}=5\times 10^{-7} \;M_{\odot} \; \mathrm{yr^{-1}}$) and slow ($v_\mathrm{w}=50 \;\mathrm{km\; s^{-1}}$) outflow can reproduce G2's properties. A faster outflow ($v_\mathrm{w}=400 \;\mathrm{km\; s^{-1}}$) might also be able to explain the material that seems to follow G2 on the same orbit.

With the help of 3D AMR hydrodynamical simulations we aim at understanding G2's nature, recent evolution and fate in the coming years. By exploring the possible parameter space of the diffuse cloud scenario, we find that a starting point within the disc of young stars is favoured by the observations, which may hint at G2 being the result of stellar wind interactions.

Massive galaxies at high redshift are predicted to be fed from the cosmic web by narrow, dense, cold streams. These streams penetrate supersonically through the hot medium encompassed by a stable shock near the virial radius of the dark-matter halo. Our long-term goal is to explore the heating and dissipation rate of the streams and their fragmentation and possible breakup, in order to understand how galaxies are fed, and how this affects their star-formation rate and morphology. We present here the first step, where we analyze the linear Kelvin-Helmholtz instability (KHI) of a cold, dense slab or cylinder flowing through a hot, dilute medium in the transonic regime. The current analysis is limited to the adiabatic case with no gravity and assuming equal pressure in the stream and the medium. By analytically solving the linear dispersion relation, we find a transition from a dominance of the familiar rapidly growing surface modes in the subsonic regime to more slowly growing body modes in the supersonic regime. The system is parameterized by three parameters: the density contrast between the stream and the medium, the Mach number of stream velocity with respect to the medium, and the stream width with respect to the halo virial radius. We find that a realistic choice for these parameters places the streams near the mode transition, with the KHI exponential-growth time in the range 0.01-10 virial crossing times for a perturbation wavelength comparable to the stream width. We confirm our analytic predictions with idealized hydrodynamical simulations. Our linear-KHI estimates thus indicate that KHI may in principle be effective in the evolution of streams by the time they reach the galaxy. More definite conclusions await the extension of the analysis to the nonlinear regime and the inclusion of cooling, thermal conduction, the halo potential well, self-gravity and magnetic fields.

This paper continues our study of the foreground population to the Orion molecular clouds. The goal is to characterize the foreground population north of NGC 1981 and to investigate the star formation history in the large Orion star-forming region. We focus on a region covering about 25 square degrees, centered on the $\epsilon$ Orionis supergiant (HD 37128, B0\,Ia) and covering the Orion Belt asterism. We used a combination of optical (SDSS) and near-infrared (2MASS) data, informed by X-ray (\textit{XMM-Newton}) and mid-infrared (WISE) data, to construct a suite of color-color and color-magnitude diagrams for all available sources. We then applied a new statistical multiband technique to isolate a previously unknown stellar population in this region. We identify a rich and well-defined stellar population in the surveyed region that has about 2\,000 objects that are mostly M stars. We infer the age for this new population to be at least 5\, Myr and likely $\sim10$\,Myr and estimate a total of about 2\,500 members, assuming a normal IMF. This new population, which we call the Orion Belt population, is essentially extinction-free, disk-free, and its spatial distribution is roughly centered near $\epsilon$ Ori, although substructure is clearly present. The Orion Belt population is likely the low-mass counterpart to the Ori OB Ib subgroup. Although our results do not rule out Blaauw's sequential star formation scenario for Orion, we argue that the recently proposed blue streams scenario provides a better framework on which one can explain the Orion star formation region as a whole. We speculate that the Orion Belt population could represent the evolved counterpart of an Orion nebula-like cluster.

We exploit deep integral-field spectroscopic observations with KMOS/VLT of 240 star-forming disks at 0.6 < z < 2.6 to dynamically constrain their mass budget. Our sample consists of massive ($\gtrsim 10^{9.8} M_\odot$) galaxies with sizes $R_e \gtrsim 2$ kpc. By contrasting the observed velocity and dispersion profiles to dynamical models, we find that on average the stellar content contributes $32^{+8}_{-7}\%$ of the total dynamical mass, with a significant spread among galaxies (68th percentile range f_star ~ 18 - 62%). Including molecular gas as inferred from CO- and dust-based scaling relations, the estimated baryonic mass adds up to $56^{+17}_{-12}\%$ of total for the typical galaxy in our sample, reaching ~ 90% at z > 2. We conclude that baryons make up most of the mass within the disk regions of high-redshift star-forming disk galaxies, with typical disks at z > 2 being strongly baryon-dominated within $R_e$. Substantial object-to-object variations in both stellar and baryonic mass fractions are observed among the galaxies in our sample, larger than what can be accounted for by the formal uncertainties in their respective measurements. In both cases, the mass fractions correlate most strongly with measures of surface density. High $\Sigma_{star}$ galaxies feature stellar mass fractions closer to unity, and systems with high inferred gas or baryonic surface densities leave less room for additional mass components other than stars and molecular gas. Our findings can be interpreted as more extended disks probing further (and more compact disks probing less far) into the dark matter halos that host them.

Large-scale gaseous filaments with length up to the order of 100 pc are on the upper end of the filamentary hierarchy of the Galactic interstellar medium. Their association with respect to the Galactic structure and their role in Galactic star formation are of great interest from both observational and theoretical point of view. Previous "by-eye" searches, combined together, have started to uncover the Galactic distribution of large filaments, yet inherent bias and small sample size limit conclusive statistical results to be drawn. Here, we present (1) a new, automated method to identify large-scale velocity-coherent dense filaments, and (2) the first statistics and the Galactic distribution of these filaments. We use a customized minimum spanning tree algorithm to identify filaments by connecting voxels in the position-position-velocity space, using the Bolocam Galactic Plane Survey spectroscopic catalog. In the range of $7.^{\circ}5 \le l \le 194^{\circ}$, we have identified 54 large-scale filaments and derived mass ($\sim 10^3 - 10^5 \, M_\odot$), length (10-276 pc), linear mass density (54-8625 $M_\odot \, \rm{pc}^{-1}$), aspect ratio, linearity, velocity gradient, temperature, fragmentation, Galactic location and orientation angle. The filaments concentrate along major spiral arms. They are widely distributed across the Galactic disk, with 50% located within $\pm$20 pc from the Galactic mid-plane and 27% run in the center of spiral arms (aka "bones"). An order of 1% of the molecular ISM is confined in large filaments. Massive star formation is more favorable in large filaments compared to elsewhere. This is the first comprehensive catalog of large filaments useful for a quantitative comparison with spiral structures and numerical simulations.

(Abridged) We derive an analytical formula which provides estimates on multiscale gravitational energy distribution using the observed surface density PDF. Our analytical formalism also enables one to convert the observed column density PDF into an estimated volume density PDF, and to obtain average radial density profile $\rho(r)$. For a region with $N_{\rm col} \sim N^{-\gamma_{\rm N}}$, the gravitational energy spectra is $E_{\rm p}(k)\sim k^{-4(1 - 1/\gamma_{\rm N})}$. We apply the formula to observations of molecular clouds, and find that a scaling index of $-2$ of the surface density PDF implies that $\rho \sim r^{-2}$ and $E_{\rm p}(k) \sim k^{-2}$. The results are valid from the cloud scale (a few parsec) to around $\sim 0.1 \;\rm pc$. Because of the resemblance the scaling index of the gravitational energy spectrum and the that of the kinetic energy power spectrum of the Burgers turbulence (where $E\sim k^{-2}$), our result indicates that gravity can act effectively against turbulence over a multitude of physical scales. This is the critical scaling index which divides molecular clouds into two categories: clouds like Orion and Ophiuchus have shallower power laws, and the amount of gravitational energy is too large for turbulence to be effective inside the cloud. Because gravity dominates, we call this type of cloud g-type clouds. On the other hand, clouds like the California molecular cloud and the Pipe nebula have steeper power laws, and turbulence can overcome gravity if it can cascade effectively from the large scale. We call this type of cloud t-type clouds. The analytical formula can be used to determine if gravity is dominating cloud evolution when the column density probability distribution function (PDF) can be reliably determined.

The outer stellar halos of galaxies contain vital information about the formation history of galaxies, since the relaxation timescales in the outskirts are long enough to keep the memory, while the information about individual formation events in the central parts has long been lost due to mixing, star formation and relaxation. To unveil some of the information encoded in these faint outer halo regions, we study the stellar outskirts of galaxies selected from a fully hydrodynamical high resolution cosmological simulation, called Magneticum. We find that the density profiles of the outer stellar halos of galaxies over a broad mass range can be well described by an Einasto profile. For a fixed total mass range, the free parameters of the Einasto fits are closely correlated. Galaxies which had more (dry) merger events tend to have lesser curved outer stellar halos, however, we find no indication that the amount of curvature is correlated with galaxy morphology. The Einasto-like shape of the outer stellar halo densities can also explain the observed differences between the Milky Way and Andromeda outer stellar halos.

We analyze the angular momenta of massive star forming galaxies (SFGs) at the peak of the cosmic star formation epoch (z~0.8-2.6). Our sample of ~360 log(M*/Msun) ~ 9.3-11.8 SFGs is mainly based on the KMOS3D and SINS/zC-SINF surveys of H$\alpha$ kinematics, and collectively provides a representative subset of the massive star forming population. The inferred halo scale angular momentum distribution is broadly consistent with that theoretically predicted for their dark matter halos, in terms of mean spin parameter <$\lambda$> ~ 0.037 and its dispersion ($\sigma_{log(\lambda)}$~0.2). Spin parameters correlate with the disk radial scale, and with their stellar surface density, but do not depend significantly on halo mass, stellar mass, or redshift. Our data thus support the long-standing assumption that on average, even at high redshifts, the specific angular momentum of disk galaxies reflects that of their dark matter halos (j_d = j_DM). The lack of correlation between $\lambda$ x (j_d/j_DM) and the nuclear stellar density $\Sigma_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to "compaction" inside massive high-redshift disks. For our sample, the inferred average stellar-to-dark matter mass ratio is ~2%, consistent with abundance matching results. Including the molecular gas, the total baryonic disk-to-dark matter mass ratio is ~5% for halos near $10^{12}$ Msun, which corresponds to 31% of the cosmologically available baryons, implying that high-redshift disks are strongly baryon dominated.

(Abridged) Many of the observed CO line profiles exhibit broad linewidths that greatly exceed the thermal broadening expected within molecular clouds. These suprathermal CO linewidths are assumed to be originated from the presence of unresolved supersonic motions inside clouds. Typically overlooked in the literature, in this paper we aim to quantify the impact of the opacity broadening effects on the current interpretation of the CO suprathermal line profiles. Without any additional contributions to the gas velocity field, a large fraction of the apparently supersonic (${\cal M}\sim$2-3) linewidths measured in both $^{12}$CO and $^{13}$CO (J=1-0) lines can be explained by the saturation of their corresponding sonic-like, optically-thin C$^{18}$O counterparts assuming standard isotopic fractionation. Combined with the presence of multiple components detected in our C$^{18}$O spectra, these opacity effects seem to be also responsible of the highly supersonic linewidths (${\cal M}>$8-10) detected in the broadest $^{12}$CO and $^{13}$CO spectra in Taurus. Our results demonstrate that most of the suprathermal $^{12}$CO and $^{13}$CO linewidths could be primarily created by a combination of opacity broadening effects and multiple gas velocity components blended in these saturated emission lines. Once corrected by their corresponding optical depth, each of these gas components present transonic intrinsic linewidths consistently traced by the three CO isotopologues within a factor of 2. Highly correlated and velocity-coherent at large scales, the largest and highly supersonic velocity differences inside clouds are generated by the relative motions between individual gas components. This highly discretized structure of the molecular gas traced in CO suggest that the gas dynamics inside molecular clouds could be better described by the properties of a fully-resolved macroscopic turbulence.

Gravity plays important roles in the evolution of molecular clouds. We present an acceleration mapping method to estimate the acceleration induced by gravitational interactions in molecular clouds based on observational data. We find that the geometry of a region has a significant impact on the behavior of gravity. In the Pipe nebula which can be approximated as a gas filament, we find that gravitational acceleration can effectively compress the end of this filament, which may have triggered star formation. We identify this as the "gravitational focusing" effect proposed by Burkert & Hartman (2004). In the sheet-like IC348-B3 region, gravity can lead to collapse at its edge, while in the centrally condensed NGC1333 cluster-forming region gravity can drive accretion towards the center. In general, gravitational acceleration tends to be enhanced in the localized regions around the ends of the filaments and the edges of sheet-like structures. Neglecting magnetic fields, these "gravitational focusing" and "edge collapse" effects can promote the formation of dense gas in a timescale that is much shorter than the global dynamical time. Since the interstellar medium is in general structured, these edge effects should be prevalent.

We present measurements of the [NII]/Ha ratio as a probe of gas-phase oxygen abundance for a sample of 419 star-forming galaxies at z=0.6-2.7 from the KMOS3D near-IR multi-IFU survey. The mass-metallicity relation (MZR) is determined consistently with the same sample selection, metallicity tracer, and methodology over the wide redshift range probed by the survey. We find good agreement with long-slit surveys in the literature, except for the low-mass slope of the relation at z~2.3, where this sample is less biased than previous samples based on optical spectroscopic redshifts. In this regime we measure a steeper slope than some literature results. Excluding the AGN contribution from the MZR reduces sensitivity at the high mass end, but produces otherwise consistent results. There is no significant dependence of the [NII]/Ha ratio on SFR or environment at fixed redshift and stellar mass. The IFU data allow spatially resolved measurements of [NII]/Ha, from which we can infer abundance gradients for 180 galaxies, thus tripling the current sample in the literature. The observed gradients are on average flat, with only 15 gradients statistically offset from zero at >3sigma. We have modelled the effect of beam-smearing, assuming a smooth intrinsic radial gradient and known seeing, inclination and effective radius for each galaxy. Our seeing-limited observations can recover up to 70% of the intrinsic gradient for the largest, face-on disks, but only 30% for the smaller, more inclined galaxies. We do not find significant trends between observed or corrected gradients and any stellar population, dynamical or structural galaxy parameters, mostly in agreement with existing studies with much smaller sample sizes. In cosmological simulations, strong feedback is generally required to produce flat gradients at high redshift.

Observations of the gas component of the cloud G2 in the Galactic Center have revealed its connection to a tail (G2t) lying on the same orbit. More recent studies indicate a connection between G2 and G1, another cloud detected on the blueshifted side of G2's orbit, suggesting a scenario in which G2 is a denser clump in a stream of gas. In this Letter we show that a simulation of an outflow by a central source (possibly a T Tauri star) moving on G2's orbit and interacting with a hot atmosphere surrounding SgrA* can have G2 and G2t as a byproduct. G2 would be the bow-shock formed in the head of the source, while G2t might be the result of the stripping of the rest of the shocked material by the ram pressure of the surrounding hot gas and of its successive accumulation in the trailing region. Mock position-velocity diagrams for the Br$\gamma$ emission for this simulation can indeed reproduce the correct position and velocity of G2t, as well as the more tenuous material in between. Though some tension between the observations and the simulated model remains, we argue that this might be due to issues in the construction of observed position-velocity (PV) diagrams and/or to a poor treatment of some physical processes - like hydrodynamic mixing - in our simulation.

Giant clumps are a characteristic feature of observed high-redshift disk galaxies. We propose that these kpc-sized clumps have a complex substructure and are the result of many smaller clumps self-organizing themselves into clump clusters (CC). This bottom-up scenario is in contrast to the common top-down view that these giant clumps form first and then sub fragment. Using a high resolution hydrodynamical simulation of an isolated, fragmented massive gas disk and mimicking the observations from Genzel et al. (2011) at $z \sim 2$, we find remarkable agreement in many details. The CCs appear as single entities of sizes $R_{\mathrm{HWHM}} \simeq 0.9-1.4$ kpc and masses $\sim 1.5-3 \times 10^9 \ \mathrm{M_{\odot}}$ representative of high-z observations. They are organized in a ring around the center of the galaxy. The origin of the observed clumps' high intrinsic velocity dispersion $\sigma_{\mathrm{intrinsic}} \simeq 50 - 100 \ \mathrm{km \ s^{-1}}$ is fully explained by the internal irregular motions of their substructure in our simulation. No additional energy input, e.g. via stellar feedback, is necessary. Furthermore, in agreement with observations, we find a small velocity gradient $V_{\mathrm{grad}} \simeq 8 - 27 \ \mathrm{km \ s^{-1} \ kpc^{-1}}$ along the CCs in the beam smeared velocity residual maps which corresponds to net prograde and retrograde rotation with respect to the rotation of the galactic disk. The CC scenario could have strong implications for the internal evolution, lifetimes and the migration timescales of the observed giant clumps, bulge growth and AGN activity, stellar feedback and the chemical enrichment history of galactic disks.

We study the dynamical response of extended systems, hosts, to smaller systems, satellites, orbiting around the hosts using extremely high-resolution N-body simulations with up to one billion particles. This situation corresponds to minor mergers which are ubiquitous in the scenario of hierarchical structure formation in the universe. According to Chandrasekhar (1943), satellites create density wakes along the orbit and the wakes cause a deceleration force on satellites, i.e. dynamical friction. This study proposes an analytical model to predict the dynamical response of hosts as reflected in their density distribution and finds not only traditional wakes but also mirror images of over- and underdensities centered on the host. Our controlled N-body simulations with high resolutions verify the predictions of the analytical model. We apply our analytical model to the expected dynamical response of nearby interacting galaxy pairs, the Milky Way - Large Magellanic Cloud system and the M31 - M33 system.

Stellar winds and supernova (SN) explosions of massive stars ("stellar feedback") create bubbles in the interstellar medium (ISM) and insert newly produced heavy elements and kinetic energy into their surroundings, possibly driving turbulence. Most of this energy is thermalized and immediately removed from the ISM by radiative cooling. The rest is available for driving ISM dynamics. In this work we estimate the amount of feedback energy retained as kinetic energy when the bubble walls have decelerated to the sound speed of the ambient medium. We show that the feedback of the most massive star outweighs the feedback from less massive stars. For a giant molecular cloud (GMC) mass of 1e5 solar masses (as e.g. found in the Orion GMCs) and a star formation efficiency of 8% the initial mass function predicts a most massive star of approximately 60 solar masses. For this stellar evolution model we test the dependence of the retained kinetic energy of the cold GMC gas on the inclusion of stellar winds. In our model winds insert 2.34 times the energy of a SN and create stellar wind bubbles serving as pressure reservoirs. We find that during the pressure driven phases of the bubble evolution radiative losses peak near the contact discontinuity (CD), and thus, the retained energy depends critically on the scales of the mixing processes across the CD. Taking into account the winds of massive stars increases the amount of kinetic energy deposited in the cold ISM from 0.1% to a few percent of the feedback energy.

The gas cloud G2 is currently being tidally disrupted by the Galactic Centre super-massive black hole, Sgr A*. The region around the black hole is populated by $\sim 30$ Wolf-Rayet stars, which produce strong outflows. We explore the possibility that gas clumps, such as G2, originate from the collision of stellar winds via the non-linear thin shell instability. Following an analytical approach, we study the thermal evolution of slabs formed in the symmetric collision of winds, evaluating whether instabilities occur, and estimating possible clump masses. We find that the collision of relatively slow ($< 750$ km s$^{-1}$) and strong ($\sim 10^{-5}$ Msun yr$^{-1}$) stellar winds from stars at short separations ($<10$ mpc) is a process that indeed could produce clumps of G2's mass and above. Such short separation encounters of single stars along their known orbits are not common in the Galactic Centre, making this process a possible but unlikely origin for G2. We also discuss clump formation in close binaries such as IRS 16SW and in asymmetric encounters as promising alternatives that deserve further numerical study.

The evolution and distribution of the angular momentum of dark matter (DM) halos have been discussed in several studies over the past decades. In particular, the idea arose that angular momentum conservation should allow to infer the total angular momentum of the entire DM halo from measuring the angular momentum of the baryonic component, which is populating the center of the halo, especially for disk galaxies. To test this idea and to understand the connection between the angular momentum of the DM halo and its galaxy, we use the Magneticum simulations. We successfully produce populations of spheroidal and disk galaxies self-consistently. Thus, we are able to study the dependence of galactic properties on their morphology. We find that (1) the specific angular momentum of stars in disk and spheroidal galaxies as a function of their stellar mass compares well with observational results; (2) the specific angular momentum of the stars in disk galaxies is slightly smaller compared to the specific angular momentum of the cold gas, in good agreement with observations; (3) simulations including the baryonic component show a dichotomy in the specific stellar angular momentum distribution when splitting the galaxies according to their morphological type (this dichotomy can also be seen in the spin parameter, where disk galaxies populate halos with slightly larger spin compared to spheroidal galaxies); (4) disk galaxies preferentially populate halos in which the angular momentum vector of the DM component in the central part shows a better alignment to the angular momentum vector of the entire halo; and (5) the specific angular momentum of the cold gas in disk galaxies is approximately 40 percent smaller than the specific angular momentum of the total DM halo and shows a significant scatter.

The Lupus I cloud is found between the Upper-Scorpius (USco) and the Upper-Centaurus-Lupus (UCL) sub-groups of the Sco-Cen OB-association, where the expanding USco H I shell appears to interact with a bubble currently driven by the winds of the remaining B-stars of UCL.We want to study how collisions of large-scale interstellar gas flows form and influence new dense clouds in the ISM.We performed LABOCA continuum sub-mm observations of Lupus I that provide for the first time a direct view of the densest, coldest cloud clumps and cores at high angular resolution.We complemented those by Herschel and Planck data from which we constructed column density and temperature maps.We calculated PDFs to characterize the density structure of the cloud.The northern part of Lupus I is found to have on average lower densities and higher temperatures as well as no active star formation.The center-south part harbors dozens of pre-stellar cores where density and temperature reach their maximum and minimum, respectively.Our analysis of the column density PDFs from the Herschel data show double peak profiles for all parts of the cloud which we attribute to an external compression.In those parts with active star formation, the PDF shows a power-law tail at high densities.The PDFs we calculated from our LABOCA data trace the denser parts of the cloud showing one peak and a power-law tail.With LABOCA we find 15 cores with masses between 0.07 and 1.71 Msun and a total mass of ~8 Msun.The total gas and dust mass of the cloud is ~164 Msun and hence 5% of the mass is in cores.From the Herschel and Planck data we find a total mass of ~174 Msun and ~171 Msun, respectively.The position, orientation and elongated shape of Lupus I, the double peak PDFs and the population of pre-stellar and protostellar cores could be explained by the large-scale compression from the advancing USco H I shell and the UCL wind bubble.

Filaments represent a key structure during the early stages of the star formation process. Simulations show filamentary structure commonly formed before and during the formation of cores. Aims. The Serpens Core represents an ideal laboratory to test the state-of-the-art of simulations of turbulent Giant Molecular Clouds. We use Herschel observations of the Serpens Core to compute temperature and column density maps of the region. Among the simulations of Dale et al. (2012), we select the early stages of their Run I, before stellar feedback is initiated, with similar total mass and physical size as the Serpens Core. We derive temperature and column density maps also from the simulations. The observed distribution of column densities of the filaments has been analysed first including and then masking the cores. The same analysis has been performed on the simulations as well. A radial network of filaments has been detected in the Serpens Core. The analysed simulation shows a striking morphological resemblance to the observed structures. The column density distribution of simulated filaments without cores shows only a log-normal distribution, while the observed filaments show a power-law tail. The power-law tail becomes evident in the simulation if one focuses just on the column density distribution of the cores. In contrast, the observed cores show a flat distribution. Even though the simulated and observed filaments are subjectively similar-looking, we find that they behave in very different ways. The simulated filaments are turbulence-dominated regions, the observed filaments are instead self-gravitating structures that will probably fragment into cores.

The dusty, ionized gas cloud G2 is currently passing the massive black hole in the Galactic Center at a distance of roughly 2400 Schwarzschild radii. We explore the possibility of a starting point of the cloud within the disks of young stars. We make use of the large amount of new observations in order to put constraints on G2's origin. Interpreting the observations as a diffuse cloud of gas, we employ three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations with the PLUTO code and do a detailed comparison with observational data. The simulations presented in this work update our previously obtained results in multiple ways: (1) high resolution three-dimensional hydrodynamical AMR simulations are used, (2) the cloud follows the updated orbit based on the Brackett-$\gamma$ data, (3) a detailed comparison to the observed high-quality position-velocity diagrams and the evolution of the total Brackett-$\gamma$ luminosity is done. We concentrate on two unsolved problems of the diffuse cloud scenario: the unphysical formation epoch only shortly before the first detection and the too steep Brackett-$\gamma$ light curve obtained in simulations, whereas the observations indicate a constant Brackett-$\gamma$ luminosity between 2004 and 2013. For a given atmosphere and cloud mass, we find a consistent model that can explain both, the observed Brackett-$\gamma$ light curve and the position-velocity diagrams of all epochs. Assuming initial pressure equilibrium with the atmosphere, this can be reached for a starting date earlier than roughly 1900, which is close to apo-center and well within the disks of young stars.

We present Karl G. Jansky Very Large Array observations of the CO J=1-0 transition in a sample of four $z\sim2$ main sequence galaxies. These galaxies are in the blue sequence of star-forming galaxies at their redshift, and are part of the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) which imaged them in CO J=3-2. Two galaxies are imaged here at high signal-to-noise, allowing determinations of their disk sizes, line profiles, molecular surface densities, and excitation. Using these and published measurements, we show that the CO and optical disks have similar sizes in main-sequence galaxies, and in the galaxy where we can compare CO J=1-0 and J=3-2 sizes we find these are also very similar. Assuming a Galactic CO-to-H$_2$ conversion, we measure surface densities of $\Sigma_{mol}\sim1200$ M$_\odot$pc$^{-2}$ in projection and estimate $\Sigma_{mol}\sim500-900$ M$_\odot$pc$^{-2}$ deprojected. Finally, our data yields velocity-integrated Rayleigh-Jeans brightness temperature line ratios $r_{31}$ that are approximately unity. In addition to the similar disk sizes, the very similar line profiles in J=1-0 and J=3-2 indicate that both transitions sample the same kinematics, implying that their emission is coextensive. We conclude that in these two main sequence galaxies there is no evidence for significant excitation gradients or a large molecular reservoir that is diffuse or cold and not involved in active star-formation. We suggest that $r_{31}$ in very actively star-forming galaxies is likely an indicator of how well mixed the star formation activity and the molecular reservoir are.

The structure and dark matter halo core properties of dwarf spheroidal galaxies (dSphs) are investigated. A double-isothermal model of an isothermal stellar system, embedded in an isothermal dark halo core provides an excellent fit to the various observed stellar surface density distributions. The stellar system can be well characterised by King profiles with a broad distribution of concentration parameters c. The core scale length of the stellar system a_* is sensitive to the central dark matter density rho_0. In contrast to single-component systems, the cut-off radius of the stellar system, rs_t, however does not trace the tidal radius but the core radius r_c of its dark matter halo. c is therefore sensitive to the ratio of the stellar to the dark matter velocity dispersion, sigma_*/sigma_0. Simple empirical relationships are derived that allow to calculate the dark halo core parameters rho_0, r_c and sigma_0, given the observable quantities sigma_*, a_* and c. The DIS model is applied to the Milky Way's dSphs. Their halo velocity dispersions lie in a narrow range of 10km/s <= sigma_0 <= 18km/s with halo core radii of 280pc <= r_c <= 1.3kpc and r_c=2a_*. All dSphs follow closely the same universal dark halo core scaling relation rho_0*r_c=75 Msolar/pc^2 that characterises the cores of more massive galaxies over several orders of magnitude in mass. The dark matter core mass is a strong function of core radius. Inside a fixed radius r_u, with r_u the logarithmic mean of the dSph's core radii, the total enclosed mass M_u is however roughly constant, although outliers should exist. For our dSphs we find r_u=400pc and M_u=2.6*10^7 Msolar. The core densities of the Galaxy's dSphs are very high, with rho_0=0.2 Msolar/pc^3. They should therefore be tidally undisturbed. Observational evidence for tidal effects might then provide a serious challenge for the cold dark matter scenario.

Recent observations have shown that the majority of the Andromeda galaxy's satellites are aligned in a thin plane. On the theoretical side it has been proposed that galaxies acquire their gas via cold streams. In addition, numerical simulations show that the same streams also deliver gas clumps which could potentially develop into satellite galaxies. Assuming that cold streams are a major source of satellite systems around galaxies we calculate the probabilities in two different models to find a certain fraction of satellites within a thin plane around the central galaxy of the host halo with and without having the same sense of rotation. Using simple geometrical considerations and adopting a random orientation of the streams we demonstrate that the vast thin disk of satellites detected around Andromeda can naturally be explained within this framework. In fact, without any satellite scattering, two streams would lead to too many satellites in the thin plane, compared with the observations. Three streams reproduce the observations very well. Natural implications from our model are that all massive galaxies should have a thin plane of satellites and that the satellites should naturally distribute themselves not only into a single plane but into several inclined ones. We estimate the effect of additional satellites accreted from random directions and find it to be of minor relevance for a mild inflow of satellites from random directions.

High energy resolution spectroscopy of the 1.8 MeV radioactive decay line of 26Al with the SPI instrument on board the INTEGRAL satellite has recently revealed that diffuse 26Al has large velocities in comparison to other components of the interstellar medium in the Milky Way. 26Al shows Galactic rotation in the same sense as the stars and other gas tracers, but reaches excess velocities up to 300 km/s. We investigate if this result can be understood in the context of superbubbles, taking into account the statistics of young star clusters and H I supershells, as well as the association of young star clusters with spiral arms. We derive energy output and 26Al mass of star clusters as a function of the cluster mass via population synthesis from stellar evolutionary tracks of massive stars. [...] We link this to the size distribution of HI supershells and assess the properties of likely 26Al-carrying superbubbles. 26Al is produced by star clusters of all masses above about 200 solar masses, roughly equally contributed over a logarithmic star cluster mass scale, and strongly linked to the injection of feedback energy. The observed superbubble size distribution cannot be related to the star cluster mass function in a straight forward manner. In order to avoid that the added volume of all superbubbles exceeds the volume of the Milky Way, individual superbubbles have to merge frequently. If any two superbubbles merge, or if 26Al is injected off-centre in a bigger HI supershell we expect the hot 26Al-carrying gas to obtain velocities of the order of the typical sound speed in superbubbles, about 300 km/s before decay. [...]

The radially averaged metallicity distribution of the ISM and the young stellar population of a sample of 20 disk galaxies is investigated by means of an analytical chemical evolution model which assumes constant ratios of galactic wind mass loss and accretion mass gain to star formation rate. Based on this model the observed metallicities and their gradients can be described surprisingly well by the radially averaged distribution of the ratio of stellar mass to ISM gas mass. The comparison between observed and model predicted metallicity is used to constrain the rate of mass loss through galactic wind and accretion gain in units of the star formation rate. Three groups of galaxies are found: galaxies with either mostly winds and only weak accretion, or mostly accretion and only weak winds, and galaxies where winds are roughly balanced by accretion. The three groups are distinct in the properties of their gas disks. Galaxies with approximately equal rates of mass-loss and accretion gain have low metallicity, atomic hydrogen dominated gas disks with a flat spatial profile. The other two groups have gas disks dominated by molecular hydrogen out to 0.5 to 0.7 isophotal radii and show a radial exponential decline, which is on average steeper for the galaxies with small accretion rates. The rates of accretion (<1.0 x SFR) and outflow (<2.4 x SFR) are relatively low. The latter depend on the calibration of the zero point of the metallicity determination from the use of HII region strong emission lines.

Gravitational instabilities play an important role in structure formation of gas-rich high-redshift disc galaxies. In this paper, we revisit the axisymmetric perturbation theory and the resulting growth of structure by taking the realistic thickness of the disc into account. In the unstable regime, which corresponds for thick discs to a Toomre parameter below the critical value Q_0,crit = 0.696, we find a fastest growing perturbation wavelength that is always a factor 1.93 times larger than in the classical razor-thin disc approximation. This result is independent of the adopted disc scaleheight and by this independent of temperature and surface density. In order to test the analytical theory, we compare it with a high-resolution hydrodynamical simulation of an isothermal gravitationally unstable gas disc with the typical vertical sech^2 density profile and study its break up into rings that subsequently fragment into dense clumps. In the first phase, rings form, that organize themselves discretely, with distances corresponding to the local fastest growing perturbation wavelength. We find that the disc scaleheight has to be resolved initially with five or more grid cells in order to guarantee proper growth of the ring structures, which follow the analytical prediction. These rings later on contract to a thin and dense line, while at the same time accreting more gas from the inter-ring region. It is these dense, circular filaments, that subsequently fragment into a large number of clumps. Contrary to what is typically assumed, the clump sizes are therefore not directly determined by the fastest growing wavelength.

Massive stars shape the surrounding ISM by emitting ionizing photons and ejecting material through stellar winds. To study the impact of the momentum from the wind of a massive star on the surrounding neutral or ionized material, we implemented a new HEALPix-based momentum conserving wind scheme in the Smoothed Particle Hydrodynamics (SPH) code SEREN. A qualitative study of the impact of the feedback from an O7.5-like star on a self gravitating sphere shows that, on its own, the transfer of momentum from a wind onto cold surrounding gas has both a compressing and dispersing effect. It mostly affects gas at low and intermediate densities. When combined with a stellar source's ionizing UV radiation, we find the momentum driven wind to have little direct effect on the gas. We conclude that, during a massive star's main sequence, the UV ionizing radiation is the main feedback mechanism shaping and compressing the cold gas. Overall, the wind's effects on the dense gas dynamics and on the triggering of star formation are very modest. The structures formed in the ionization-only simulation and in the combined feedback simulation are remarkably similar. However, in the combined feedback case, different SPH particles end up being compressed. This indicates that the microphysics of gas mixing differ between the two feedback simulations and that the winds can contribute to the localized redistribution and reshuffling of gas.

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

We present new, deep near-infrared SINFONI @ VLT integral field spectroscopy of the gas cloud G2 in the Galactic Center, from late summer 2013 and spring 2014. G2 is visible in recombination line emission. The spatially resolved kinematic data track the ongoing tidal disruption. As expected for an observation near pericenter passage, roughly half of the gas in 2014 is found at the redshifted, pre-pericenter side of the orbit, while the other half is at the post-pericenter, blueshifted side. We also present an orbital solution for the gas cloud G1, which was discovered a decade ago in L'-band images when it was spatially almost coincident with Sgr A*. The orientation of the G1 orbit in the three angles is almost identical to the one of G2, but it has a lower eccentricity and smaller semi-major axis. We show that the observed astrometric positions and radial velocities of G1 are compatible with the G2 orbit, assuming that (i) G1 was originally on the G2 orbit preceding G2 by 13 years and (ii) a simple drag force acted on it during pericenter passage. Taken together with the previously described tail of G2, which we detect in recombination line emission and thermal broadband emission, we propose that G2 may be a bright knot in a much more extensive gas streamer. This matches purely gaseous models for G2, such as a stellar wind clump or the tidal debris from a partial disruption of a star.

The formation and evolution of tidal tails like those observed around some globular clusters and dwarf satellite galaxies is examined with an N-body simulation. In particular, we analyse in detail the evolving tidal features of a one-component satellite that is moving on a highly eccentric orbit in the external field of a host galaxy potential like our own. The results show that every time the satellite approaches apogalacticon, a fresh pair of tidal tails becomes notably prominent, and that eventually, the satellite possesses multiple tidal tails via repeating apocentre passages. Accordingly, the number of observed tidal arms can be used as a tracer of the number of orbital periods that such a system has completed around the centre of its host galaxy. By identifying the arm particles included in each of the first three consecutively formed pairs of tidal tails, we find that each pair of tidal tails is practically identical to one another regarding the energy and angular-momentum distributions. In addition, we demonstrate that the density profiles of these three pairs of tidal tails at their first apogalacticons after formation agree well with one another. It therefore follows that the multiplicity in tidal features originates from the repeated episode of tidal-arm formation in the course of the precessing motion of a satellite.

Recent studies found the densities of dark matter (DM) subhaloes which surround nearby dwarf spheroidal galaxies (dSphs) to be significantly lower than those of the most massive subhaloes expected around Milky Way sized galaxies in cosmological simulations, the so called "too-big-to-fail" (TBTF) problem. A caveat of previous work has been that dark substructures were assumed to contain steep density cusps in the center of DM haloes even though the central density structure of DM haloes is still under debate. In this study, we re-examine the TBTF problem for models of DM density structure with cores or shallowed cusps. Our analysis demonstrates that the TBTF problem is alleviated as the logarithmic slope of the central cusp becomes shallower. We find that the TBTF problem is avoided if the central cusps of DM haloes surrounding dSphs are shallower than $r^{-0.6}$.

The central engines of Seyfert galaxies are thought to be enshrouded by geometrically thick gas and dust structures. In this article, we derive observable properties for a self-consistent model of such toroidal gas and dust distributions, where the geometrical thickness is achieved and maintained with the help of X-ray heating and radiation pressure due to the central engine. Spectral energy distributions (SEDs) and images are obtained with the help of dust continuum radiative transfer calculations with RADMC-3D. For the first time, we are able to present time-resolved SEDs and images for a physical model of the central obscurer. Temporal changes are mostly visible at shorter wavelengths, close to the combined peak of the dust opacity as well as the central source spectrum and are caused by variations in the column densities of the generated outflow. Due to the three-component morphology of the hydrodynamical models -- a thin disc with high density filaments, a surrounding fluffy component (the obscurer) and a low density outflow along the rotation axis -- we find dramatic differences depending on wavelength: whereas the mid-infrared images are dominated by the elongated appearance of the outflow cone, the long wavelength emission is mainly given by the cold and dense disc component. Overall, we find good agreement with observed characteristics, especially for those models, which show clear outflow cones in combination with a geometrically thick distribution of gas and dust, as well as a geometrically thin, but high column density disc in the equatorial plane.

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

The very long and thin infrared dark cloud "Nessie" is even longer than had been previously claimed, and an analysis of its Galactic location suggests that it lies directly in the Milky Way's mid-plane, tracing out a highly elongated bone-like feature within the prominent Scutum-Centaurus spiral arm. Re-analysis of mid-infrared imagery from the Spitzer Space Telescope shows that this IRDC is at least 2, and possibly as many as 5 times longer than had originally been claimed by Nessie's discoverers (Jackson et al. 2010); its aspect ratio is therefore at least 300:1, and possibly as large as 800:1. A careful accounting for both the Sun's offset from the Galactic plane ($\sim 25$ pc) and the Galactic center's offset from the $(l^{II},b^{II})=(0,0)$ position shows that the latitude of the true Galactic mid-plane at the 3.1 kpc distance to the Scutum-Centaurus Arm is not $b=0$, but instead closer to $b=-0.4$, which is the latitude of Nessie to within a few pc. An analysis of the radial velocities of low-density (CO) and high-density (${\rm NH}_3$) gas associated with the Nessie dust feature suggests that Nessie runs along the Scutum-Centaurus Arm in position-position-velocity space, which means it likely forms a dense `spine' of the arm in real space as well. The Scutum-Centaurus arm is the closest major spiral arm to the Sun toward the inner Galaxy, and, at the longitude of Nessie, it is almost perpendicular to our line of sight, making Nessie the easiest feature to see as a shadow elongated along the Galactic Plane from our location. Future high-resolution dust mapping and molecular line observations of the harder-to-find Galactic "bones" should allow us to exploit the Sun's position above the plane to gain a (very foreshortened) view "from above" of the Milky Way's structure.

We present the correlations between stellar mass, star formation rate (SFR) and [NII]/Ha flux ratio as indicator of gas-phase metallicity for a sample of 222 galaxies at 0.8 < z < 2.6 and log(M*/Msun)=9.0-11.5 from the LUCI, SINS/zC-SINF and KMOS3D surveys. This sample provides a unique analysis of the mass-metallicity relation (MZR) over an extended redshift range using consistent data analysis techniques and strong-line metallicity indicator. We find a constant slope at the low-mass end of the relation and can fully describe its redshift evolution through the evolution of the characteristic turnover mass where the relation begins to flatten at the asymptotic metallicity. At fixed mass and redshift, our data do not show a correlation between the [NII]/Ha ratio and SFR, which disagrees with the 0.2-0.3dex offset in [NII]/Ha predicted by the "fundamental relation" between stellar mass, SFR and metallicity discussed in recent literature. However, the overall evolution towards lower [NII]/Ha at earlier times does broadly agree with these predictions.

In this study, we present a detailed, statistical analysis of black hole growth and the evolution of active galactic nuclei (AGN) using cosmological hydrodynamic simulations run down to $z=0$. The simulations self-consistently follow radiative cooling, star formation, metal enrichment, black hole growth and associated feedback processes from both supernovae typeII/Ia and AGN. We consider two simulation runs, one with a large co-moving volume of $(500\ \mathrm{Mpc})^3$ and one with a smaller volume of $(68\ \mathrm{Mpc})^3$ but with a by a factor of almost 20 higher mass resolution. Consistently with previous results, our simulations can widely match observed black hole properties of the local Universe. Furthermore, our simulations can successfully reproduce the evolution of the bolometric AGN luminosity function for both the low-luminosity and the high-luminosity end up to $z=3.0$. In addition, the smaller but higher resolution run is able to match the observational data of the low bolometric luminosity end at higher redshifts $z=3-4$. We also perform a direct comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for a torus-level obscuration, and find a similarly good agreement. These results nicely demonstrate that the observed "anti-hierarchical" trend in the AGN number density evolution (i.e. the number densities of luminous AGN peak at higher redshifts than those of faint AGN) is self-consistently predicted by our simulations. Implications of this downsizing behaviour on active black holes, their masses and Eddington-ratios are discussed. Overall, the downsizing behaviour in the AGN number density as a function of redshift can be mainly attributed to the evolution of the gas density in the resolved vicinity of a (massive) black hole. (shortened)

We report the detection of ubiquitous powerful nuclear outflows in massive (> 10^11 Msun) z~2 star-forming galaxies (SFGs), which are plausibly driven by an Active Galactic Nucleus (AGN). The sample consists of the eight most massive SFGs from our SINS/zC-SINF survey of galaxy kinematics with the imaging spectrometer SINFONI, six of which have sensitive high-resolution adaptive optics (AO) assisted observations. All of the objects are disks hosting a significant stellar bulge. The spectra in their central regions exhibit a broad component in Halpha and forbidden [NII] and [SII] line emission, with typical velocity FWHM ~ 1500 km/s, [NII]/Halpha ratio ~ 0.6, and intrinsic extent of 2 - 3 kpc. These properties are consistent with warm ionized gas outflows associated with Type 2 AGN, the presence of which is confirmed via independent diagnostics in half the galaxies. The data imply a median ionized gas mass outflow rate of ~ 60 Msun/yr and mass loading of ~ 3. At larger radii, a weaker broad component is detected but with lower FWHM ~ 485 km/s and [NII]/Halpha ~ 0.35, characteristic for star formation-driven outflows as found in the lower-mass SINS/zC-SINF galaxies. The high inferred mass outflow rates and frequent occurrence suggest the nuclear outflows efficiently expel gas out of the centers of the galaxies with high duty cycles, and may thus contribute to the process of star formation quenching in massive galaxies. Larger samples at high masses will be crucial to confirm the importance and energetics of the nuclear outflow phenomenon, and its connection to AGN activity and bulge growth.

Long wavelength sonic oscillations are observed or inferred in many of the small, dark molecular clouds, the starless cores, that are the precursors to protostars. The oscillations provide significant internal energy and the time scale for their dissipation may control the rate of star formation in the starless cores. Despite their importance, their origin is unknown. We explore one hypothesis that the oscillations develop as a starless core forms from a filament. We model this process with a numerical hydrodynamic simulation and generate synthetic molecular line observations with a radiative transfer simulation. Comparison with actual observations shows general agreement suggesting the proposed mechanism is viable.

The classical picture of a star-forming filament is a near-equilibrium structure, with collapse dependent on its gravitational criticality. Recent observations have complicated this picture, revealing filaments as a mess of apparently interacting subfilaments, with transsonic internal velocity dispersions and mildly supersonic intra-subfilament dispersions. How structures like this form is unresolved. Here we study the velocity structure of filamentary regions in a simulation of a turbulent molecular cloud. We present two main findings: first, the observed complex velocity features in filaments arise naturally in self gravitating hydrodynamic simulations of turbulent clouds without the need for magnetic or other effects. Second, a region that is filamentary only in projection and is in fact made of spatially distinct features can displays these same velocity characteristics. The fact that these disjoint structures can masquerade as coherent filaments in both projection and velocity diagnostics highlights the need to continue developing sophisticated filamentary analysis techniques for star formation observations.

Observations have revealed interesting universal properties of dark matter (DM) halos especially around low-mass galaxies. Strigari et al. (2008) showed that DM halos have common enclosed masses within 300pc (Strigari relation). Kormendy & Freeman (2004) reported DM halos having almost identical central surface densities (the $\mu_{\rm 0D}$ relation). In addition, there exists a core-cusp problem, a discrepancy of the central density distribution between simulated halos and observations. We investigate whether a scenario where cuspy halos transform into cores by some dynamical processes can also explain their universal structural properties. It is shown that a cusp-to-core transformation model naturally reproduces the $\mu_{\rm 0D}$ relation and that Strigari relation follows from the $\mu_{\rm 0D}$ relation for dwarf galaxies. We also show that the central densities of cored dark halos provide valuable information about their formation redshifts.

In 2011, we discovered a compact gas cloud ("G2") with roughly three Earth masses that is falling on a near-radial orbit toward the massive black hole in the Galactic Center. The orbit is well constrained and pericenter passage is predicted for early 2014. Our data beautifully show that G2 gets tidally sheared apart due to the massive black hole's force. During the next months, we expect that in addition to the tidal effects, hydrodynamics get important, when G2 collides with the hot ambient gas around Sgr A*. Simulations show that ultimately, the cloud's material might fall into the massive black hole. Predictions for the accretion rate and luminosity evolution, however, are very difficult due to the many unknowns. Nevertheless, this might be a unique opportunity in the next years to observe how gas feeds a massive black hole in a galactic nucleus.

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 study the origin of high-redshift, compact, quenched spheroids (red nuggets) through the dissipative shrinkage of gaseous discs into compact star-forming systems (blue nuggets). The discs, fed by cold streams, undergo violent disc instability (VDI) that drives gas into the centre (along with mergers). The inflow is dissipative when its timescale is shorter than the star formation timescale. This implies a threshold of 0.28 in the cold-to-total mass ratio within the disc radius. For the typical gas fraction 0.5 at z=2, this threshold is traced back to a maximum spin parameter of 0.05, implying that half the star-forming galaxies contract to blue nuggets, while the rest form extended stellar discs. Thus, the surface density of blue galaxies is expected to be bimodal about 10^9 Msun/kpc^2, slightly increasing with mass. Blue nuggets are expected to be rare at low z when the gas fraction is low. The blue nuggets quench to red nuggets by complementary internal and external mechanisms. Internal quenching by a compact bulge, operating as a fast mode and especially at high z, may involve starbursts, stellar and AGN feedback, or Q-quenching. Quenching due to hot-medium haloes above 10^{12} Msun provides maintenance and a slower mode at low redshift. These predictions are confirmed in simulations and are consistent with observations at z=0-3.

The origin of the dense gas cloud G2 discovered in the Galactic Center (Gillessen et al. 2012) is still a debated puzzle. G2 might be a diffuse cloud or the result of an outflow from an invisible star embedded in it. We present here detailed simulations of the evolution of winds on G2's orbit. We find that the hydrodynamic interaction with the hot atmosphere present in the Galactic Center and the extreme gravitational field of the supermassive black hole must be taken in account when modeling such a source scenario. We find that the hydrodynamic interaction with the hot atmosphere present in the Galactic Center and the extreme gravitational field of the supermassive black hole must be taken in account when modeling such a source scenario. We also find that in this scenario most of the Br\gamma\ luminosity is expected to come from the highly filamentary densest shocked wind material. G2's observational properties can be used to constrain the properties of the outflow and our best model has a mass outflow rate of Mdot,w=8.8 x 10^{-8} Msun/yr and a wind velocity of vw = 50 km/s. These values are compatible with those of a young TTauri star wind, as already suggested by Scoville & Burkert (2013).

We present a simple toy model to understand what sets the scatter in star formation and metallicity of galaxies at fixed mass. The scatter ultimately arises from the intrinsic scatter in the accretion rate, but may be substantially reduced depending on the timescale on which the accretion varies compared to the timescale on which the galaxy loses gas mass. This model naturally produces an anti-correlation between star formation and metallicity at a fixed mass, the basis of the fundamental metallicity relation. We show that observational constraints on the scatter in galaxy scaling relations can be translated into constraints on the galaxy-to-galaxy variation in the mass loading factor, and the timescales and magnitude of stochastic accretion onto star-forming galaxies. We find a remarkably small scatter in the mass loading factor, < 0.1 dex, and that the scatter in accretion rates is smaller than expected from N-body simulations.

A novel mechanism of the formation of compact dusty clouds in astrophysical environments is discussed. It is shown that the balance of collective forces operating in space dusty plasmas can result in the effect of dust self-confinement, generating equilibrium spherical clusters. The distribution of dust and plasma density inside such objects and their stability are investigated. Spherical dusty clouds can be formed in a broad range of plasma parameters, suggesting that this process of dust self-organization might be a generic phenomenon occurring in different astrophysical media. We argue that compact dusty clouds can represent condensation seeds for a population of small-scale, cold, gaseous clumps in the diffuse interstellar medium. They could play an important role in regulating its small-scale structure and its thermodynamical evolution.

Proplyds and stars inside HII-regions are a well studied phenomenon. It is possible that they were triggered by the expansion of the HII-region itself. Here, we present calculations on the dynamics of HII-regions. We show that the triggered stars that form in the expanding shell of swept up material around the HII region rarely return into the HII regions on timescales that are inferred for the proplyds and observed young stars. However, in very dense environments like Orion, the triggered stars return in time. Thus, our model can explain why proplyds are barely observed in other HII regions. We propose that the properties of young stellar objects in HII regions in general depend critically on the distance from the massive, ionizing central star cluster. Closest in, there are proplyds, where the disk of a young star interacts directly with the feedback of the massive star. Further out are Class II protostars, where the ionization already removed the envelope. Even further away, one should find Class I stars, which either have been triggered by the ionizing radiation or pre-existed and have not lost their envelope yet. This radial sequence is not necessarily an age sequence but rather a result of the dwindling importance of stellar winds and ionizing radiation with distance. We investigate the observational signature of triggered star formation and find that the stellar distribution for ionization triggered star formation shows a distinct feature, a peak at the current position of the ionization front. Therefore, it is generally possible to tell triggered and in situ distributions of stars apart.

As part of the SINS/zC-SINF surveys of high-z galaxy kinematics, we derive the radial distributions of H-alpha surface brightness, stellar mass surface density, and dynamical mass at ~2 kpc resolution in 19 z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. From these data we infer the radial distribution of the Toomre Q-parameter for these main-sequence star forming galaxies (SFGs), covering almost two decades of stellar mass (10^9.6 to 10^11.5 solar masses). In more than half of our SFGs, the H-alpha distributions cannot be fit by a centrally peaked distribution, such as an exponential, but are better described by a ring, or the combination of a ring and an exponential. At the same time the kinematic data indicate the presence of a mass distribution more centrally concentrated than a single exponential distribution for 5 of the 19 galaxies. The resulting Q-distributions are centrally peaked for all, and significantly exceed unity there for three quarters of the SFGs. The occurrence of H-alpha rings and of large nuclear Q-values is strongly correlated, and is more common for the more massive SFGs. While our sample is small and there remain substantial uncertainties and caveats, our observations are consistent with a scenario in which cloud fragmentation and global star formation are secularly suppressed in gas rich high-z disks from the inside out, as the central stellar mass density of the disks grows.

We investigate the differential effects of metal cooling and galactic stellar winds on the cosmological formation of individual galaxies with three sets of cosmological, hydrodynamical zoom simulations of 45 halos in the mass range 10^11<M_halo<10^13M_sun. Models including both galactic winds and metal cooling (i) suppress early star formation at z>1 and predict reasonable star formation histories, (ii) produce galaxies with high cold gas fractions (30-60 per cent) at high redshift, (iii) significantly reduce the galaxy formation efficiencies for halos (M_halo<10^12M_sun) at all redshifts in agreement with observational and abundance matching constraints, (iv) result in high-redshift galaxies with reduced circular velocities matching the observed Tully-Fisher relation at z~2, and (v) significantly increase the sizes of low-mass galaxies (M_stellar<3x10^10M_sun) at high redshift resulting in a weak size evolution - a trend in agreement with observations. However, the low redshift (z<0.5) star formation rates of massive galaxies are higher than observed (up to ten times). No tested model predicts the observed size evolution for low-mass and high-mass galaxies simultaneously. Due to the delayed onset of star formation in the wind models, the metal enrichment of gas and stars is delayed and agrees well with observational constraints. Metal cooling and stellar winds are both found to increase the ratio of in situ formed to accreted stars - the relative importance of dissipative vs. dissipationless assembly. For halo masses below ~10^12M_sun, this is mainly caused by less stellar accretion and compares well to predictions from semi-analytical models but still differs from abundance matching models. For higher masses, the fraction of in situ stars is over-predicted due to the unrealistically high star formation rates at low redshifts.

The origin of the dense gas cloud G2 discovered in the Galactic Center (Gillessen et al. 2012) is still a debated puzzle. G2 might be a diffuse cloud or the result of an outflow from an invisible star embedded in it. We present hydrodynamical simulations of the evolution of different spherically symmetric winds of a stellar object embedded in G2. We find that the interaction with the ambient medium and with the extreme gravitational field of the supermassive black hole in the Galactic Center must be taken into account for such a source scenario. The thermal pressure of the hot and dense atmosphere confines the wind, while its ram pressure shapes it via stripping along the orbit, with the details depending on the wind parameters. Tidal forces squeeze the wind near pericenter, reducing it to a thin and elongated filament. We also find that in this scenario most of the Br\gamma\ luminosity is expected to come from the densest part of the wind, which has a highly filamentary structure with low filling factor. The observations can be best matched by a mass outflow rate of Mdot,w=8.8 x 10^{-8} Msun/yr and a wind velocity of vw = 50 km/s. These values are compatible with those of a young T Tauri star wind, as already suggested by (Scoville & Burkert 2013).

Massive stars form in groups and their winds and supernova explosions create superbubbles up to kpc in size. Their ejecta are important for the dynamics of the interstellar medium and chemical evolution models. However, ejecta kinematics and the characteristic scales in space and time are rather unexplored beyond ~10 ka. Through measurement of radioactive 26Al with its decay time constant of ~10^6 years, we aim to trace the kinematics of massive-star and supernova ejecta over million-year time scales. Our goal is to identify the mixing time scale and the spatio-kinematics of such ejecta from the pc to kpc scale. We use the SPI spectrometer on INTEGRAL and its observations along the Galactic ridge to trace the line shape of the 1808.63 keV gamma-ray line from 26Al decay. We determine line centroids and compare these to Doppler shift expectations from large-scale systematic rotation around the Galaxy's center, as observed in other Galactic objects. We measure the radial velocities of gas traced by 26Al, averaged over the line of sight, as a function of Galactic longitude. We find substantially higher velocities than expected from Galactic rotation, the average bulk velocity being ~200 km s^-1 larger than the Galactic-rotation prediction. The observed radial velocity spread implies a Doppler-broadening of the gamma-ray line that is consistent with our measurements of the overall line width. We can reproduce the observed characteristics with 26Al sources located along the inner spiral arms, when we add a global blow-out into the forward direction away from arms, such as expected from massive stars offset towards the spiral-arm leading edge. The superbubbles-halo connection implies angular-momentum transfer in the disk-halo system and also radial gas flows.

Based on high-resolution, spatially resolved data of 10 z~2 star-forming galaxies from the SINS/zC-SINF survey and LUCI data for 12 additional galaxies, we probe the excitation properties of high-z galaxies and the impact of active galactic nuclei (AGN), shocks and photoionization. We explore how these spatially-resolved line ratios can inform our interpretation of integrated emission line ratios obtained at high redshift. Many of our galaxies fall in the `composite' region of the z~0 [NII]/Ha versus [OIII]/Hb diagnostic (BPT) diagram, between star-forming galaxies and those with AGN. Based on our resolved measurements, we find that some of these galaxies likely host an AGN, while others appear to be affected by the presence of shocks possibly caused by an outflow or from enhanced ionization parameter as compared with HII regions in normal local star-forming galaxies. We find that the Mass-Excitation (MEx) diagnostic, which separates purely star-forming and AGN hosting local galaxies in the [OIII]/Hb versus stellar mass plane, does not properly separate z~2 galaxies classified according to the BPT diagram. However, if we shift the galaxies based on the offset between the local and z~2 mass-metallicity relation (i.e. to the mass they would have at z~0 with the same metallicity), we find better agreement between the MEx and BPT diagnostics. Finally, we find that metallicity calibrations based on [NII]/Ha are more biased by shocks and AGN at high-z than the [OIII]/Hb/[NII]/Ha calibration.

Studies by Lada (2010) and Heiderman (2010) have suggested that star formation mostly occurs above a threshold in gas surface density Sigma of Sigma_c = 120 Msun pc^{-2} (A_K = 0.8). Heiderman infer a threshold by combining low-mass star-forming regions, which show a steep increase in the star formation rate per unit area Sigma_SFR with increasing Sigma, and massive cores forming luminous stars which show a linear relation. We argue that these observations do not require a particular density threshold. The steep dependence of Sigma_SFR, approaching unity at protostellar core densities, is a natural result of the increasing importance of self-gravity at high densities along with the corresponding decrease in evolutionary timescales. The linear behavior of Sigma_SFR vs. Sigma in massive cores is consistent with probing dense gas in gravitational collapse, forming stars at a characteristic free-fall timescale given by the use of a particular molecular tracer. The low-mass and high-mass regions show different correlations between gas surface density and the area A spanned at that density, with A=Sigma^{-3} for low-mass regions and A=Sigma^{-1} for the massive cores; this difference, along with the use of differing techniques to measure gas surface density and star formation, suggests that connecting the low-mass regions with massive cores is problematic. We show that the approximately linear relationship between dense gas mass and stellar mass used by Lada similarly does not demand a particular threshold for star formation, and requires continuing formation of dense gas. Our results are consistent with molecular clouds forming by galactic hydrodynamic flows with subsequent gravitational collapse.

We perform hydrodynamical simulations of a young galactic disc embedded in a hot gaseous halo using parameters typical for Lyman break galaxies (LBGs). We take into account the (static) gravitational potentials due to a dark matter halo, a stellar bulge and a disc of stars and gas. Star formation is treated by a local Kennicutt-Schmidt law. We simplify the structure of the interstellar medium by restricting the computational domain to a 25th of the full azimuthal angle, effectively assuming large-scale axisymmetry and neglecting any effects of spiral structure, and focus on the large-scale ISM drivers, the superbubbles. Supernovae are triggered randomly and have preset event sizes of several tens to hundreds. We further investigate different halo gas pressures and energy injection methods. Many of our simulated galaxies, but not all, develop bipolar outflows. We characterise the strength of the outflow by mass and energy outflow rates, and investigate the effect of changes to the details of the model. We find that supernovae are more effective if comprised into larger superbubbles. The weight and the pressure of the halo gas is able to quench galactic outflows. The wind emerges from a series of superbubbles in regions where a critical star formation density is exceeded. The superbubbles expand into the gaseous halo at slightly supersonic speed, producing radiative shock waves with similar characteristics as the absorptions systems observed around LBGs.

We have further followed the evolution of the orbital and physical properties of G2, the object currently falling toward the massive black hole in the Galactic Center on a near-radial orbit. New, very sensitive data were taken in April 2013 with NACO and SINFONI at the ESO VLT . The 'head' of G2 continues to be stretched ever further along the orbit in position-velocity space. A fraction of its emission appears to be already emerging on the blue-shifted side of the orbit, past pericenter approach. Ionized gas in the head is now stretched over more than 15,000 Schwarzschild radii RS around the pericenter of the orbit, at ~ 2000 RS ~ 20 light hours from the black hole. The pericenter passage of G2 will be a process stretching over a period of at least one year. The Brackett-{\gamma} luminosity of the head has been constant over the past 9 years, to within +- 25%, as have the line ratios Brackett-{\gamma} / Paschen-{\alpha} and Brackett-{\gamma} / Helium-I. We do not see any significant evidence for deviations of G2's dynamical evolution, due to hydrodynamical interactions with the hot gas around the black hole, from a ballistic orbit of an initially compact cloud with moderate velocity dispersion. The constant luminosity and the increasingly stretched appearance of the head of G2 in the position-velocity plane, without a central peak, is not consistent with several proposed models with continuous gas release from an initially bound zone around a faint star on the same orbit as G2.

Over the past 10 Gyr, star-forming galaxies have changed dramatically, from clumpy and gas rich, to rather quiescent stellar-dominated disks with specific star formation rates lower by factors of a few tens. We present a general theoretical model for how this transition occurs, and what physical processes drive it, making use of 1D axisymmetric thin disk simulations with an improved version of the Gravitational Instability-Dominated Galaxy Evolution Tool (GIDGET) code. We show that at every radius galaxies tend to be in a slowly evolving equilibrium state wherein new accretion is balanced by star formation, galactic winds, and radial transport of gas through the disk by gravitational instability (GI) -driven torques. The gas surface density profile is determined by which of these terms are in balance at a given radius - direct accretion is balanced by star formation and galactic winds near galactic centers, and by transport at larger radii. We predict that galaxies undergo a smooth transition from a violent disk instability phase to secular evolution. This model provides a natural explanation for the high velocity dispersions and large clumps in z~2 galaxies, the growth and subsequent quenching of bulges, and features of the neutral gas profiles of local spiral galaxies.

We present a high resolution simulation of an idealized model to explain the origin of the two young, counter-rotating, sub-parsec scale stellar disks around the supermassive black hole SgrA* at the Center of the Milky Way. In our model, the collision of a single molecular cloud with a circum-nuclear gas disk (similar to the one observed presently) leads to multiple streams of gas flowing towards the black hole and creating accretion disks with angular momentum depending on the ratio of cloud and circum-nuclear disk material. The infalling gas creates two inclined, counter-rotating sub-parsec scale accretion disks around the supermassive black hole with the first disk forming roughly 1 Myr earlier, allowing it to fragment into stars and get dispersed before the second, counter-rotating disk forms. Fragmentation of the second disk would lead to the two inclined, counter-rotating stellar disks which are observed at the Galactic Center. A similar event might be happening again right now at the Milky Way Galactic Center. Our model predicts that the collision event generates spiral-like filaments of gas, feeding the Galactic Center prior to disk formation with a geometry and inflow pattern that is in agreement with the structure of the so called mini-spiral that has been detected in the Galactic Center.

Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us to study the star formation efficiency on subgalactic scales around redshift z = 1.2. We present a method for obtaining the gas and star formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding scales are not resolved. This method is based on identifying these structures in position-velocity diagrams corresponding to slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star forming distant galaxies - EGS13003805, EGS13004291, EGS12007881, and EGS13019128 in the AEGIS terminology - to determine the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [OII] line. We identify the ensembles of clumps with the angular resolution available in both CO and [OII] spectroscopy; i.e., 1-1.5". SFR and gas surface densities are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation on a scale of ~8 kpc. The data generally indicates an average depletion time of 1.9 Gyr, but with significant variations from point to point within the galaxies.

Optical and near-infrared images of the Carina Nebula show a peculiar arc-shaped feature, which we call the "Sickle", next to the B-type star Trumpler 14 MJ 218. We use multi-wavelength observations to explore and constrain the nature and origin of the nebulosity. Using sub-mm data from APEX/LABOCA as well as Herschel far-infrared maps, we discovered a dense, compact clump with a mass of ~ 40 Msun located close to the apex of the Sickle. We investigate how the B-star MJ 218, the Sickle, and the clump are related. Our numerical simulations show that, in principle, a B-type star located near the edge of a clump can produce a crescent-shaped wind shock front, similar to the observed morphology. However, the observed proper motion of MJ 218 suggest that the star moves with high velocity (~ 100 km/s) through the ambient interstellar gas. We argue that the star is just about to graze along the surface of the clump, and the Sickle is a bow shock induced by the stellar wind, as the object moves supersonically through the density gradient in the envelope of the clump.

We report matched resolution, imaging spectroscopy of the CO J=3-2 line (with the IRAM Plateau de Bure millimeter interferometer) and of the H-alpha line (with LUCI at the Large Binocular Telescope)in the massive z=1.53 main-sequence galaxy EGS 13011166, as part of the "Plateau de Bure high-z, blue sequence survey (PHIBSS). We combine these data with HST V-J-J-H-band maps to derive spatially resolved distributions of stellar surface density, star formation rate, molecular gas surface density, optical extinction and gas kinematics. The spatial distribution and kinematics of the ionized and molecular gas are remarkably similar and are well modeled by a turbulent, globally Toomre unstable rotating disk. The stellar surface density distribution is smoother than the clumpy rest-frame UV/optical light distribution, and peaks in an obscured, star forming massive bulge near the dynamical center. The molecular gas surface density and the effective optical screen extinction track each other and are well modeled by a 'mixed' extinction model. The inferred slope of the spatially resolved molecular gas to star formation rate relation depends strongly on the adopted extinction model and can vary from 0.8 to 1.7. For the preferred mixed dust-gas model we find a near linear slope.

We analyze the spectra, spatial distributions and kinematics of Ha, [NII] and [SII] emission in a sample of 42, z~2.2 UV/optically selected star forming galaxies (SFGs) from the SINS & zC-SINF surveys, 35 of which were observed in the adaptive optics mode of SINFONI. This is supplemented by kinematic data from 48 z~1-2.5 galaxies from the literature. We find that the kinematic classification of the high-z SFGs as `dispersion dominated' or `rotation dominated' correlates most strongly with their intrinsic sizes. Smaller galaxies are more likely `dispersion-dominated' for two main reasons: 1) The rotation velocity scales linearly with galaxy size but intrinsic velocity dispersion does not depend on size, and as such, their ratio is systematically lower for smaller galaxies, and 2) Beam smearing strongly decreases large-scale velocity gradients and increases observed dispersion much more for galaxies with sizes at or below the resolution. Dispersion dominated SFGs may thus have intrinsic properties similar to `rotation dominated' SFGs, but are primarily more compact, lower mass, less metal enriched and may have higher gas fractions, plausibly because they represent an earlier evolutionary state.

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 explore the possibility that the G2 gas cloud falling in towards \sgra is the mass loss envelope of a young TTauri star. As the star plunges to smaller radius at 1000 to 6000 \kms, a strong bow shock forms where the stellar wind is impacted by the hot X-ray emitting gas in the vicinity of \sgra. For a stellar mass loss rate of $4\times10^{-8}$ \msun per yr and wind velocity 100 \kms, the bow shock will have an emission measure ($EM = n^2 vol$) at a distance $\sim10^{16}$ cm, similar to that inferred from the IR emission lines. The ionization of the dense bow shock gas is potentially provided by collisional ionization at the shock front and cooling radiation (X-ray and UV) from the post shock gas. The former would predict a constant line flux as a function of distance from \sgra, while the latter will have increasing emission at lesser distances. In this model, the star and its mass loss wind should survive pericenter passage since the wind is likely launched at 0.2 AU and this is much less than the Roche radius at pericenter ($\sim3$ AU for a stellar mass of 2\msun). In this model, the emission cloud will probably survive pericenter passage, discriminating this scenario from others.

We analyse the motion of single optically thick clouds in the potential of a central mass under the influence of an anisotropic radiation field ~|cos(\theta)|, a model applicable to the inner region of active galactic nuclei. Resulting orbits are analytically soluble for constant cloud column densities. All stable orbits are closed, although they have non-trivial shapes. Furthermore, there exists a stability criterion in the form of a critical inclination, which depends on the luminosity of the central source and the column density of the cloud.

We investigate the origin and fate of the recently discovered gas cloud G2 close to the Galactic Center. Our hydrodynamical simulations focussing on the dynamical evolution of the cloud in combination with currently available observations favour two scenarios: a Compact Cloud which started around the year 1995 and an extended Spherical Shell of gas, with an apocenter distance within the disk(s) of young stars. The former is able to explain the detected signal of G2 in the position-velocity-diagram of the year 2008.5 and 2011.5 data. The latter can account for both, G2's signal as well as the fainter extended tail-like structure G2t seen at larger distances to the black hole and smaller velocities. From these first idealised simulations we expect a rise of the current activity of Sgr A* shortly after the closest approach and a constant feeding through a nozzle-like structure over a long period. The near future evolution of the cloud will be a sensitive probe of the conditions of the gas distribution in the milli-parsec environment of the massive black hole in the Galactic Center and will also give us invaluable information of the feeding of black holes and the activation of the central source.

We study the evolution of isolated and merging interstellar bubbles of three stars (25, 32 and 60 M\odot) in a homogeneous background medium with a density of 10 mp /ccm via 3D-hydrodynamic simulations with standard ISM thermodynamics (optically thin radiative cooling and photo-electric heating) and time dependent energy and mass input according to stellar evolutionary tracks. We vary the position of the three stars relative to each other to compare the energy response for cases of isolated, merging and initially cospatial bubbles. Due to mainly the Vishniac instability, our simulated bubbles develop thick shells and filamentary internal structures in column density. The shell widths reach tens of per cent of the outer bubble radius, which compares favourably to observations. More energy is retained in the ISM for more closely packed groups, by up to a factor of three and typically a factor of two for intermediate times after the first supernova. Once the superbubble is established, different positions of the contained stars make only a minor difference to the energy tracks. For our case of three massive stars, the energy deposition varies only very little for distances up to about 30 pc between the stars. Energy injected by supernovae is entirely dissipated in a superbubble on a timescale of about 1 Myr, which increases slightly with the superbubble size at the time of the explosion. The Vishniac instability may be responsible for the broadening of the shells of interstellar bubbles. Massive star winds are significant energetically due to their - in the long run - more efficient, steady energy injection and because they evacuate the space around the massive stars. For larger scale simulations, the feedback effect of close groups of stars or clusters may be subsumed into one effective energy input with insignificant loss of energy accuracy.

Based on SINFONI Ha, [NII] and [SII] AO data of 30 z \sim 2 star-forming galaxies (SFGs) from the SINS and zcSINF surveys, we find a strong correlation of the Ha broad flux fraction with the star formation surface density of the galaxy, with an apparent threshold for strong outflows occurring at 1 Msun yr^-1 kpc^-2. Above this threshold, we find that SFGs with logm_\ast>10 have similar or perhaps greater wind mass loading factors (eta = Mdotout/SFR) and faster outflow velocities than lower mass SFGs. This trend suggests that the majority of outflowing gas at z \sim 2 may derive from high-mass SFGs, and that the z \sim 2 mass-metallicity relation is driven more by dilution of enriched gas in the galaxy gas reservoir than by the efficiency of outflows. The mass loading factor is also correlated with the SFR and inclination, such that more star-forming and face-on galaxies launch more powerful outflows. For galaxies that have evidence for strong outflows, we find that the broad emission is spatially extended to at least the half-light radius (\sim a few kpc). We propose that the observed threshold for strong outflows and the observed mass loading of these winds can be explained by a simple model wherein break-out of winds is governed by pressure balance in the disk. Using the ratio of the [SII] doublet in a broad and narrow component, we find that outflowing gas has a density of \sim10-100 cm^-3, significantly less than that of the star forming gas (600 cm^-3).

We present PHIBSS, the IRAM Plateau de Bure high-z blue sequence CO 3-2 survey of the molecular gas properties in normal star forming galaxies (SFGs) near the cosmic star formation peak. PHIBSS provides 52 CO detections in two redshift slices at z~1.2 and 2.2, with log(M*(M_solar))>10.4 and log(SFR(M_solar/yr))>1.5. Including a correction for the incomplete coverage of the M*-SFR plane, we infer average gas fractions of ~0.33 at z~1.2 and ~0.47 at z~2.2. Gas fractions drop with stellar mass, in agreement with cosmological simulations including strong star formation feedback. Most of the z~1-3 SFGs are rotationally supported turbulent disks. The sizes of CO and UV/optical emission are comparable. The molecular gas - star formation relation for the z=1-3 SFGs is near-linear, with a ~0.7 Gyrs gas depletion timescale; changes in depletion time are only a secondary effect. Since this timescale is much less than the Hubble time in all SFGs between z~0 and 2, fresh gas must be supplied with a fairly high duty cycle over several billion years. At given z and M*, gas fractions correlate strongly with the specific star formation rate. The variation of specific star formation rate between z~0 and 3 is mainly controlled by the fraction of baryonic mass that resides in cold gas.

Observational studies have revealed a "downsizing" trend in black hole (BH) growth: the number densities of luminous AGN peak at higher redshifts than those of faint AGN. This would seem to imply that massive black holes formed before low mass black holes, in apparent contradiction to hierarchical clustering scenarios. We investigate whether this observed "downsizing" in BH growth is reproduced in a semi-analytic model for the formation and evolution of galaxies and black holes, set within the hierarchical paradigm for structure formation (Somerville et al. 2008; S08). In this model, black holes evolve from light seeds (\sim100M\odot) and their growth is merger-driven. The original S08 model (baseline model) reproduces the number density of AGN at intermediate redshifts and luminosities, but underproduces luminous AGN at very high redshift (z > 3) and overproduces them at low redshift (z < 1). In addition, the baseline model underproduces low-luminosity AGN at low redshift (z < 1). To solve these problems we consider several modifications to the physical processes in the model: (1) a 'heavy' black hole seeding scenario (2) a sub-Eddington accretion rate ceiling that depends on the cold gas fraction, and (3) an additional black hole accretion mode due to disk instabilities. With these three modifications, the models can explain the observed downsizing, successfully reproduce the bolometric AGN luminosity function and simultaneously reproduce galaxy and black hole properties in the local Universe. We also perform a comparison with the observed soft and hard X-ray luminosity functions of AGN, including an empirical correction for torus-level obscuration, and reach similar conclusions. Our best-fit model suggests a scenario in which disk instabilities are the main driver for moderately luminous Seyfert galaxies at low redshift, while major mergers are the main trigger for luminous AGN.

The 3D observed velocities of the Large and Small Magellanic Clouds(LMC and SMC) provide an opportunity to probe the Galactic potential in the outskirt of the Galactic halo. Based on a canonical NFW model of the Galactic potential, Besla et al.(2007) reconstructed LMC and SMC's orbits and suggested that they are currently on their first perigalacticon passage about the Galaxy. Motivated by several recent revisions of the Sun's motion around the Galactic center, we re-examine the LMC's orbital history and show that it depends sensitively on the dark-matter's mass distribution beyond its present Galactic distance. We utilize results of numerical simulations to consider a range of possible structural and evolutionary models for the Galactic potentials. We find that within the theoretical and observational uncertainties, it is possible for the LMC to have had multiple perigalacticon passages on the Hubble time scale, especially if the Galactic circular velocity at the location of the Sun is greater than $\sim 228$km s$^{-1}$. Based on these models, a more accurate determination of the LMC's motion may be used to determine the dark matter distribution in the outskirt of the Galactic halo.

We present new observations of the recently discovered gas cloud G2 currently falling towards the massive black hole in the Galactic Center. The new data confirm that G2 is on a highly elliptical orbit with a predicted pericenter passage mid 2013. The updated orbit has an even larger eccentricity of 0.966, an epoch of pericenter two months later than estimated before, and a nominal minimum distance of 2200 Schwarzschild radii only. The velocity gradient of G2 has developed further to 600 km/s FWHM in summer 2012. We also detect the tail of similar total flux and on the same orbit as G2 along the trajectory at high significance. No hydrodynamic effects are detected yet, since the simple model of a tidally shearing gas cloud still describes the data very well. The flux of G2 has not changed by more than 10% between 2008 and 2012, disfavoring models where additional gas from a reservoir is released to the disrupting diffuse gas component.

Hydrodynamic stability has been a longstanding issue for the cloud model of the broad line region in active galactic nuclei. We argue that the clouds may be gravitationally bound to the supermassive black hole. If true, stabilisation by thermal pressure alone becomes even more difficult. We further argue that if magnetic fields should be present in such clouds at a level that could affect the stability properties, they need to be strong enough to compete with the radiation pressure on the cloud. This would imply magnetic field values of a few Gauss for a sample of Active Galactic Nuclei we draw from the literature. We then investigate the effect of several magnetic configurations on cloud stability in axi-symmetric magnetohydrodynamic simulations. For a purely azimuthal magnetic field which provides the dominant pressure support, the cloud first gets compressed by the opposing radiative and gravitational forces. The pressure inside the cloud then increases, and it expands vertically. Kelvin-Helmholtz and column density instability lead to a filamentary fragmentation of the cloud. This radiative dispersion continues until the cloud is shredded down to the resolution level. For a helical magnetic field configuration, a much more stable cloud core survives with a stationary density histogram which takes the form of a power law. Our simulated clouds develop sub-Alfvenic internal motions on the level of a few hundred km/s.

We report on the discovery of strong tidal features around a dwarf spheroidal galaxy in the Hydra I galaxy cluster, indicating its ongoing tidal disruption. This very low surface brightness object, HCC-087, was originally classified as an early-type dwarf in the Hydra Cluster Catalogue (HCC), but our re-analysis of the ESO-VLT/FORS images of the HCC unearthed a clear indication of an S-shaped morphology and a large spatial extent. Its shape, luminosity (M_V=-11.6 mag), and physical size (at a half-light radius of 3.1 kpc and a full length of ~5.9 kpc) are comparable to the recently discovered NGC 4449B and the Sagittarius dwarf spheroidal, all of which are undergoing clear tidal disruption. Aided by N-body simulations we argue that HCC-087 is currently at its first apocenter, at 150 kpc, around the cluster center and that it is being tidally disrupted by the galaxy cluster's potential itself. An interaction with the near-by (50 kpc) S0 cluster galaxy HCC-005, at M* ~ 3 x 10^10 M_sun is rather unlikely, as this constellation requires a significant amount of dynamical friction and thus low relative velocities. The S-shaped morphology and large spatial extent of the satellite would, however, also appear if HCC-087 would orbit the cluster center. These features appear to be characteristic properties of satellites that are seen in the process of being tidally disrupted, independent of the environment of the destruction. An important finding of our simulations is an orientation of the tidal tails perpendicular to the orbit.

We analyse the results of four simulations of isolated galaxies: two with a rigid spiral potential of fixed pattern speed, but with different degrees of star-formation induced feedback, one with an axisymmetric galactic potential and one with a `live' self-gravitating stellar component. Since we use a Lagrangian method we are able to select gas that lies within giant molecular clouds (GMCs) at a particular timeframe, and to then study the properties of this gas at earlier and later times. We find that gas which forms GMCs is not typical of the interstellar medium at least 50 Myr before the clouds form and reaches mean densities within an order of magnitude of mean cloud densities by around 10 Myr before. The gas in GMCs takes at least 50 Myr to return to typical ISM gas after dispersal by stellar feedback, and in some cases the gas is never fully recycled. We also present a study of the two-dimensional, vertically-averaged velocity fields within the ISM. We show that the velocity fields corresponding to the shortest timescales (that is, those timescales closest to the immediate formation and dissipation of the clouds) can be readily understood in terms of the various cloud formation and dissipation mechanisms. Properties of the flow patterns can be used to distinguish the processes which drive converging flows (e.g.\ spiral shocks, supernovae) and thus molecular cloud formation, and we note that such properties may be detectable with future observations of nearby galaxies.

We revisit collisionless major and minor mergers of spheroidal galaxies in the context of the size evolution of elliptical galaxies. The simulations are performed as a series of mergers with mass-ratios of 1:1 and 1:10 for models representing pure bulges as well as bulges embedded in dark matter halos. For major and minor mergers, respectively, we identify and analyze two different processes, violent relaxation and stripping, leading to size evolution and a change of the dark matter fraction within the observable effective radius. Violent relaxation - which is the dominant mixing process for major mergers but less important for minor mergers - scatters relatively more dark matter particles than bulge particles to small. Stripping in minor mergers assembles stellar satellite particles at large radii in halo dominated regions of the massive host. This strongly increases the size of the bulge into regions with higher dark matter fractions leaving the inner host structure almost unchanged. A factor of two mass increase by minor mergers increases the dark matter fraction by 20 per cent. We present analytic corrections to simple one-component virial estimates for the evolution of the gravitational radii. If such a two-component system grows by minor mergers alone its size growth, $r_{\mathrm{e}} \propto M^\alpha$, reaches values of $\alpha \approx 2.4$, significantly exceeding the simple theoretical limit of $\alpha = 2$. For major mergers the sizes grow with $\alpha \lesssim 1$. Our results indicate that minor mergers of galaxies embedded in massive dark matter halos provide a potential mechanism for explaining the rapid size growth and the build-up of massive elliptical systems predicting significant dark matter fractions and radially biased velocity dispersions at large radii (abbreviated)

We investigate the origin and fate of the recently discovered gas cloud G2 close to the Galactic Center. Our hydrodynamical simulations focussing on the dynamical evolution of the cloud in combination with currently available observations favor two scenarios: a Compact Cloud which started around the year 1995 and a Spherical Shell of gas, with an apocenter distance within the disk(s) of young stars and a radius of a few times the size of the Compact Cloud. The former is able to explain the detected signal of G2 in the position-velocity diagram of the Br gamma emission of the year 2008.5 and 2011.5 data. The latter can account for both, G2's signal as well as the fainter extended tail-like structure G2t seen at larger distances from the black hole and smaller velocities. In contrast, gas stripped from a compact cloud by hydrodynamical interactions is not able to explain the location of the detected G2t emission in the observed position-velocity diagrams. This favors the Spherical Shell Scenario and might be a severe problem for the Compact Cloud as well as the so-called Compact Source Scenario. From these first idealized simulations we expect a roughly constant feeding of the supermassive black hole through a nozzle-like structure over a long period, starting shortly after the closest approach in 2013.51 for the Compact Cloud. If the matter accretes in the hot accretion mode, we do not expect a significant boost of the current activity of Sgr A* for the Compact Cloud model, but a boost of the average infrared and X-ray luminosity by roughly a factor of 80 for the Spherical Shell scenario with order of magnitude variations on a timescale of a few months. The near-future evolution of the cloud will be a sensitive probe of the conditions of the gas distribution in the milli-parsec environment of the massive black hole in the Galactic Center.

Observations of disk galaxies at z~2 have demonstrated that turbulence driven by gravitational instability can dominate the energetics of the disk. We present a 1D simulation code, which we have made publicly available, that economically evolves these galaxies from z~2 to z~0 on a single CPU in a matter of minutes, tracking column density, metallicity, and velocity dispersions of gaseous and multiple stellar components. We include an H$_2$ regulated star formation law and the effects of stellar heating by transient spiral structure. We use this code to demonstrate a possible explanation for the existence of a thin and thick disk stellar population and the age-velocity dispersion correlation of stars in the solar neighborhood: the high velocity dispersion of gas in disks at z~2 decreases along with the cosmological accretion rate, while at lower redshift, the dynamically colder gas forms the low velocity dispersion stars of the thin disk.

We have obtained high-resolution data of the z 2 ring-like, clumpy star-forming galaxy (SFG) ZC406690 using the VLT/SINFONI with AO (in K-band) and in seeing-limited mode (in H- and J-band). Our data includes all of the main strong optical emission lines: [OII], [OIII], Ha, Hb, [NII] and [SII]. We find broad, blueshifted Ha and [OIII] emission line wings in the spectra of the galaxy's massive, star-forming clumps (sigma \sim 85 km s^-1) and even broader wings (up to 70% of the total Ha flux, with sigma \sim 290 km s^-1) in regions spatially offset from the clumps by \sim 2 kpc. The broad emission likely originates from large-scale outflows with mass outflow rates from individual clumps that are 1-8x the SFR of the clumps. Based on emission line ratio diagnostics ([NII]/Ha and [SII]/Ha) and photoionization and shock models, we find that the emission from the clumps is due to a combination of photoionization from the star-forming regions and shocks generated in the outflowing component, with 5-30% of the emission deriving from shocks. In terms of the ionization parameter (6x10^7-10^8 cm/s, based on both the SFR and the O32 ratio), density (local electron densities of 300-1800 cm^-3 in and around the clumps, and ionized gas column densities of 1200-8000 Msol/pc^2), and SFR (10-40 Msol/yr), these clumps more closely resemble nuclear starburst regions of local ULIRGs and dwarf irregulars than HII regions in local galaxies. However, the star-forming clumps are not located in the nucleus as in local starburst galaxies but instead are situated in a ring several kpc from the center of their high-redshift host galaxy, and have an overall disk-like morphology. The two brightest clumps are quite different in terms of their internal properties, energetics and relative ages, and thus we are given a glimpse at two different stages in the formation and evolution of rapidly star-forming giant clumps at high-z.

We investigate the structure of the Milky Way by determining how features in a spatial map correspond to CO features in a velocity map. We examine structures including logarithmic spiral arms, a ring and a bar. We explore the available parameter space, including the pitch angle of the spiral arms, radius of a ring, and rotation curve. We show that surprisingly, a spiral arm provides a better fit to the observed molecular ring than a true ring feature. This is because both a spiral arm, and the observed feature known as the molecular ring, are curved in velocity longitude space. We find that much of the CO emission in the velocity longitude map can be fitted by a nearly symmetric 2 armed spiral pattern. One of the arms corresponds to the molecular ring, whilst the opposite arm naturally reproduces the Perseus arm. Multiple arms also contribute to further emission in the vicinity of the molecular ring and match other observed spiral arms. Whether the Galactic structure consists primarily of two, or several spiral arms, the presence of 2 symmetric logarithmic spirals, which begin in the vicinity of the ends of the bar, suggest a spiral density wave associated with the bar.

The origin, structure and evolution of the small gas cloud, G2, is investigated, that is on an orbit almost straight into the Galactic central supermassive black hole (SMBH). G2 is a sensitive probe of the hot accretion zone of Sgr A*, requiring gas temperatures and densities that agree well with models of captured shock-heated stellar winds. Its mass is equal to the critical mass below which cold clumps would be destroyed quickly by evaporation. Its mass is also constrained by the fact that at apocenter its sound crossing timescale was equal to its orbital timescale. Our numerical simulations show that the observed structure and evolution of G2 can be well reproduced if it formed in pressure equilibrium with the surrounding in 1995 at a distance from the SMBH of 7.6e16 cm. If the cloud would have formed at apocenter in the 'clockwise' stellar disk as expected from its orbit, it would be torn into a very elongated spaghetti-like filament by 2011 which is not observed. This problem can be solved if G2 is the head of a larger, shell-like structure that formed at apocenter. Our numerical simulations show that this scenario explains not only G2's observed kinematical and geometrical properties but also the Br_gamma observations of a low surface brightness gas tail that trails the cloud. In 2013, while passing the SMBH G2 will break up into a string of droplets that within the next 30 years mix with the surrounding hot gas and trigger cycles of AGN activity.

Measurements of stellar orbits provide compelling evidence that the compact radio source Sagittarius A* at the Galactic Centre is a black hole four million times the mass of the Sun. With the exception of modest X-ray and infrared flares, Sgr A* is surprisingly faint, suggesting that the accretion rate and radiation efficiency near the event horizon are currently very low. Here we report the presence of a dense gas cloud approximately three times the mass of Earth that is falling into the accretion zone of Sgr A*. Our observations tightly constrain the cloud's orbit to be highly eccentric, with an innermost radius of approach of only ~3,100 times the event horizon that will be reached in 2013. Over the past three years the cloud has begun to disrupt, probably mainly through tidal shearing arising from the black hole's gravitational force. The cloud's dynamic evolution and radiation in the next few years will probe the properties of the accretion flow and the feeding processes of the super-massive black hole. The kilo-electronvolt X-ray emission of Sgr A* may brighten significantly when the cloud reaches pericentre. There may also be a giant radiation flare several years from now if the cloud breaks up and its fragments feed gas into the central accretion zone.

We address the internal support against total free-fall collapse of the giant clumps that form by violent gravitational instability in high-z disc galaxies. Guidance is provided by an analytic model, where the proto-clumps are cut from a rotating disc and collapse to equilibrium while preserving angular momentum. This model predicts prograde clump rotation. This is confirmed in hydro-AMR zoom-in simulations of galaxies in a cosmological context. In most high-z clumps, the centrifugal force dominates the support, R=Vrot^2/Vcirc^2 > 0.5, where Vrot is the rotation velocity and Vcirc is the circular velocity. The clump spin indeed tends to be in the sense of the global disc angular momentum, but substantial tilts are frequent. Most clumps are in Jeans equilibrium, with the rest of the support provided by turbulence. Simulations of isolated gas-rich discs that resolve the clump substructure reveal that the cosmological simulations may overestimate R by ~30%, but the dominance of rotational support at high-z is not a resolution artifact. In turn, isolated gas-poor disc simulations produce at z=0 smaller gaseous non-rotating transient clouds, indicating that the difference in rotational support is associated with the fraction of cold baryons in the disc. In our current cosmological simulations, the clump rotation velocity is typically Vrot~100 km/s, but when beam smearing of \geq 0.1 arcsec is imposed, the rotation signal is reduced to a small gradient of \leq 30 km/s/kpc across the clump. The velocity dispersion in the simulated clumps is comparable to the disc dispersion so it is expected to leave only a marginal signal. Retrograde minor-merging galaxies could lead to massive clumps that do not show rotation.Testable predictions of the scenario as simulated are that the mean stellar age of the clumps, and the stellar fraction, are declining linearly with distance from the disc center.

We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor {\alpha}CO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log({\alpha}CO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.

The fate of metals ejected by young OB associations into the Interstellar Medium (ISM) is investigated numerically. In particular, we study the enrichment of the cold gas phase, which is the material that forms molecular clouds. Following previous work, the expansion and collision of two supershells in a diffuse ISM is simulated, in this case also introducing an advected quantity which represents the metals expelled by the young stars. We adopt the simplest possible approach, not differentiating between metals coming from stellar winds and those coming from supernovae. Even though the hot, diffuse phase of the ISM receives a significant amount of metals from the stars, the cold phase is efficiently shielded, with very little metal enrichment. Significant enrichment of the cold ISM will therefore be delayed by at least the cooling time of this hot phase. No variations in cloud metallicity with distance from the OB association or with direction are found, which means that the shell collision does little to enhance the metallicity of the cold clumps. We conclude that the stellar generation that forms out of molecular structures, triggered by shell collisions cannot be significantly enriched.

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.

We present a detailed comparison between numerical cosmological hydrodynamic zoom simulations and semi-analytic models (SAMs) run within merger trees extracted from the simulations. The high-resolution simulations represent 48 individual halos with virial masses in the range 2.4*10^11M_sun < M_Halo < 3.3*10^13M_sun. They include radiative H & He cooling, photo-ionization, star formation and thermal SN feedback. We compare with different SAM versions including only this complement of physical processes, and also ones including supernova driven winds, metal cooling, and feedback from AGN. Our analysis is focused on the cosmic evolution of the baryon content in galaxies and its division into various components (stars, cold gas, and hot gas). Both the SAMs and simulations are compared with observational relations between halo mass and stellar mass, and between stellar mass and star formation rate, at low and high redshift. The simulations turn out to have much higher star formation efficiencies (by about a factor of ten) than the SAMs. Therefore the cold gas is consumed much more rapidly in the simulations and stars form much earlier. Also, simulations show a transition between stellar mass growth that is dominated by in situ formation of stars to growth that is predominantly through accretion of stars formed in external galaxies. In SAMs, stellar growth is always dominated by in situ star formation. In addition, SAMs overpredict the overall gas accretion rates relative to the simulations, and overestimate the fraction of "hot" relative to "cold" accretion. We discuss the reasons for these discrepancies, and identify several physical processes that are missing in our SAM. We also highlight physical processes that are neglected in the simulations studied here, but which appear to be crucial in order to understand the properties of real galaxies.

The kinematics of distant galaxies, from z=0.1 to z>2, play a key role in our understanding of galaxy evolution from early times to the present. One of the important parameters is the intrinsic, or local, velocity dispersion of a galaxy, which allows one to quantify the degree of non-circular motions such as pressure support. However, this is difficult to measure because the observed dispersion includes the effects of (often severe) beam smearing on the velocity gradient. Here we investigate four methods of measuring the dispersion that have been used in the literature, to assess their effectiveness at recovering the intrinsic dispersion. We discuss the biasses inherent in each method, and apply them to model disk galaxies in order to determine which methods yield meaningful quantities, and under what conditions. All the mean weighted dispersion estimators are affected by (residual) beam smearing. In contrast, the dispersion recovered by fitting a spatially and spectrally convolved disk model to the data is unbiassed by the beam smearing it is trying to compensate. Because of this, and because the bias it does exhibit depends only on the signal-to-noise, it can be considered reliable. However, at very low signal-to-noise, all methods should be used with caution.

We perform calculations of isolated disc galaxies to investigate how the properties of the ISM, the nature of molecular clouds, and the global star formation rate depend on the level of stellar feedback. We adopt a simple physical model, which includes a galactic potential, a standard cooling and heating prescription of the ISM, and self gravity of the gas. Stellar feedback is implemented by injecting energy into dense, gravitationally collapsing gas, but is independent of the Schmidt-Kennicutt relation. We obtain fractions of gas, and filling factors for different phases of the ISM in reasonable ageement with observations. Supernovae are found to be vital to reproduce the scale heights of the different components of the ISM, and velocity dispersions. The GMCs formed in the simulations display mass spectra similar to the observations, their normalisation dependent on the level of feedback. We find ~40 per cent of the clouds exhibit retrograde rotation, induced by cloud-cloud collisions. The star formation rates we obtain are in good agreement with the observed Schmidt-Kennicutt relation, and are not strongly dependent on the star formation efficiency we assume, being largely self regulated by the feedback. We also investigate the effect of spiral structure by comparing calculations with and without the spiral component of the potential. The main difference with a spiral potential is that more massive GMCs are able to accumulate in the spiral arms. Thus we are able to reproduce massive GMCs, and the spurs seen in many grand design galaxies, even with stellar feedback. The presence of the spiral potential does not have an explicit effect on the star formation rate, but can increase the star formation rate indirectly by enabling the formation of long-lived, strongly bound clouds.

We study the influence of gas metallicity, turbulence, and non-equilibrium chemistry on the evolution of the two-phase interstellar medium (warm and cold atomic phases), and thereby constrain the initial conditions for star formation prevailing in turbulent gas. We perform high-resolution simulations in three dimensions, including a realistic non-equilibrium treatment of the ionization state of the gas, and examine both driven and decaying turbulence. This allows us to explore variations in the metallicity Z. In this paper, we study solar metallicity, Z=Z_sun, and low metallicity, Z=1d-3 Z_sun, gas. For driven, large-scale turbulence, we find that the influence of the metallicity on the amount of mass in the cold gas component is small. However, in decaying turbulent conditions this picture is much changed. While cold regions survive in the case of solar metallicity, they are quickly heated and dispersed in low-metallicity gas. This result suggests that star formation can be suppressed in environments of low metallicity, unless a strong turbulent driver is acting on time scales shorter than a few turbulent crossing times. Inter alia this finding could explain the overall inefficient star formation as well as the burst-like mode of star formation found in metal-poor, gas-rich systems like dwarf galaxies.

We investigate the evolution of dusty gas clouds falling into the centre of an active Seyfert nucleus. Two-dimensional high-resolution radiation hydrodynamics simulations are performed to study the fate of single clouds and the interaction between two clouds approaching the Active Galactic Nucleus. We find three distinct phases of the evolution of the cloud: (i) formation of a lenticular shape with dense inner rim caused by the interaction of gravity and radiation pressure (the lense phase), (ii) formation of a clumpy sickle-shaped structure as the result of a converging flow (the clumpy sickle phase) and (iii) a filamentary phase caused by a rapidly varying optical depth along the sickle. Depending on the column density of the cloud, it will either be pushed outwards or its central (highest column density) parts move inwards, while there is always some material pushed outwards by radiation pressure effects. The general dynamical evolution of the cloud can approximately be described by a simple analytical model.

We have studied the properties of giant star forming clumps in five z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. The clumps reside in disk regions where the Toomre Q-parameter is below unity, consistent with their being bound and having formed from gravitational instability. Broad H{\alpha}/[NII] line wings demonstrate that the clumps are launching sites of powerful outflows. The inferred outflow rates are comparable to or exceed the star formation rates, in one case by a factor of eight. Typical clumps may lose a fraction of their original gas by feedback in a few hundred million years, allowing them to migrate into the center. The most active clumps may lose much of their mass and disrupt in the disk. The clumps leave a modest imprint on the gas kinematics. Velocity gradients across the clumps are 10-40 km/s/kpc, similar to the galactic rotation gradients. Given beam smearing and clump sizes, these gradients may be consistent with significant rotational support in typical clumps. Extreme clumps may not be rotationally supported; either they are not virialized, or they are predominantly pressure supported. The velocity dispersion is spatially rather constant and increases only weakly with star formation surface density. The large velocity dispersions may be driven by the release of gravitational energy, either at the outer disk/accreting streams interface, and/or by the clump migration within the disk. Spatial variations in the inferred gas phase oxygen abundance are broadly consistent with inside-out growing disks, and/or with inward migration of the clumps.

The most recent observational evidence seems to indicate that giant molecular clouds are predominantly gravitationally unbound objects. In this paper we show that this is a natural consequence of a scenario in which cloud-cloud collisions and stellar feedback regulate the internal velocity dispersion of the gas, and so prevent global gravitational forces from becoming dominant. Thus, while the molecular gas is for the most part gravitationally unbound, local regions within the denser parts of the gas (within the clouds) do become bound and are able to form stars. We find that the observations, in terms of distributions of virial parameters and cloud structures, can be well modelled provided that the star formation efficiency in these bound regions is of order 5 - 10 percent. We also find that in this picture the constituent gas of individual molecular clouds changes over relatively short time scales, typically a few Myr.

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.

The expansion and collision of two wind-blown superbubbles is investigated numerically. Our models go beyond previous simulations of molecular cloud formation from converging gas flows by exploring this process with realistic flow parameters, sizes and timescales. The superbubbles are blown by time-dependent winds and supernova explosions, calculated from population synthesis models. They expand into a uniform or turbulent diffuse medium. We find that dense, cold gas clumps and filaments form naturally in the compressed collision zone of the two superbubbles. Their shapes resemble the elongated, irregular structure of observed cold, molecular gas filaments and clumps. At the end of the simulations, between 65 and 80 percent of the total gas mass in our simulation box is contained in these structures. The clumps are found in a variety of physical states, ranging from pressure equilibrium with the surrounding medium to highly under-pressured clumps with large irregular internal motions and structures which are rotationally supported.

The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fraction. The galaxies segregate into dispersion-dominated systems with 1<vmax/sigma<3, maximum rotational velocities vmax<200 km/s and disk half-light radii rd=1-3 kpc and rotation-dominated systems with vmax>200 km/s, vmax/sigma>3 and rd=4-8 kpc. For the dispersion-dominated sample, radial pressure gradients partly compensate the gravitational force, reducing the rotational velocities. Including this pressure effect in the MMW model, dispersion-dominated galaxies can be fitted well with spin parameters lf lambda_d=0.03-0.05 for high disk mass fractions of m_d=0.2 and with lambda_d=0.01-0.03 for m_d=0.05. These values are in good agreement with cosmological expectations. For the rotation-dominated sample however pressure effects are small and better agreement with theoretically expected disk spin parameters can only be achieved if the dark halo mass contribution in the visible disk regime (2-3*rd) is smaller than predicted by the MMW model. We argue that these galaxies can still be embedded in standard cold dark matter halos if the halos did not contract adiabatically in response to disk formation. It is shown that the observed high turbulent gas motions of the galaxies are consistent with a Toomre instability parameter Q=1 which is equal to the critical value, expected for gravitational disk instability to be the major driver of turbulence. The dominant energy source of turbulence is then the potential energy of the gas in the disk.

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.

We investigate the global dynamic stability of spherical clouds in the Broad Line Region (BLR) of Active Galactic Nuclei (AGN), exposed to radial radiation pressure, gravity of the central black hole (BH), and centrifugal forces assuming the clouds adapt their size according to the local pressure. We consider both, isotropic and anisotropic light sources. In both cases, stable orbits exist also for very sub-Keplerian rotation for which the radiation pressure contributes substantially to the force budget. We demonstrate that highly eccentric, very sub-Keplerian stable orbits may be found. This gives further support for the model of Marconi et al. 2008, who pointed out that black hole masses might be significantly underestimated if radiation pressure is neglected. That model improved the agreement between black hole masses derived in certain active galaxies based on BLR dynamics, and black hole masses derived by other means in other galaxies by inclusion of a luminosity dependent term. For anisotropic illumination, energy is conserved for averages over long time intervals, only, but not for individual orbits. This leads to Rosetta orbits that are systematically less extended in the direction of maximum radiation force. Initially isotropic relatively low column density systems would therefore turn into a disk when an anisotropic AGN is switched on.

We derive the evolution equations describing a thin axisymmetric disk of gas and stars with an arbitrary rotation curve that is kept in a state of marginal gravitational instability and energy equilibrium due to the balance between energy released by accretion and energy lost due to decay of turbulence. Rather than adopt a parameterized alpha prescription, we instead use the condition of marginal gravitational instability to self-consistently determine the position- and time-dependent transport rates. We show that there is a steady-state configuration for disks dominated by gravitational instability, and that this steady state persists even when star formation is taken into account if the accretion rate is sufficiently large. For disks in this state we analytically determine the velocity dispersion, surface density, and rates of mass and angular momentum transport as a function of the gas mass fraction, the rotation curve, and the rate of external accretion onto the disk edge. We show that disks that are initially out of steady state will evolve into it on the viscous timescale of the disk, which is comparable to the orbital period if the accretion rate is high. Finally, we discuss the implications of these results for the structure of disks in a broad range of environments, including high redshift galaxies, the outer gaseous disks of local galaxies, and accretion disks around protostars.

We study the structural evolution of turbulent molecular clouds under the influence of ionizing radiation emitted from a nearby massive star by performing a high resolution parameter study with the iVINE code. The temperature is taken to be 10K or 100K, the mean number density is either 100cm^3 or 300cm^3. Besides, the turbulence is varied between Mach 1.5 and Mach 12.5 and the main driving scale between 1pc and 8pc. We vary the ionizing flux by an order of magnitude. In our simulations the ionizing radiation enhances the initial turbulent density distribution and thus leads to the formation of pillar-like structures observed adjacent to HII regions in a natural way. Gravitational collapse occurs regularly at the tips of the structures. We find a clear correlation between the initial state of the turbulent cold cloud and the final morphology and physical properties of the structures formed. The most favorable regime for the formation of pillars is Mach 4-10. Structures and therefore stars only form if the initial density contrast between the high density unionized gas and the gas that is going to be ionized is lower than the temperature contrast between the hot and the cold gas. The density of the resulting pillars is determined by a pressure equilibrium between the hot and the cold gas. A thorough analysis of the simulations shows that the complex kinematical and geometrical structure of the formed elongated filaments reflects that of observed pillars to an impressive level of detail. In addition, we find that the observed line-of sight velocities allow for a distinct determination of different formation mechanisms. Comparing the current simulations to previous results and recent observations we conclude that e.g. the pillars of creation in M16 formed by the mechanism proposed here and not by the radiation driven implosion of pre-existing clumps.

Elliptical, lenticular, and early-type spiral galaxies show a remarkably tight power-law correlation between the mass M_BH of their central supermassive black hole (SMBH) and the number N_GC of globular clusters: M_BH=m*N_GC^(1.08+/-0.04) with m=1.7*10^5 solar masses. Thus, to a good approximation the SMBH mass is the same as the total mass of the globular clusters. Based on a limited sample of 13 galaxies, this relation appears to be a better predictor of SMBH mass (rms scatter 0.2 dex) than the M_BH-sigma relation between SMBH mass and velocity dispersion sigma. The small scatter reflects the fact that galaxies with high globular cluster specific frequency S_N tend to harbor SMBHs that are more massive than expected from the M_BH-sigma relation.

The paper assumes radiation forces proportional to distance between equal temperature clouds. However, we assume there are clouds in any direction. The forces then cancel almost entirely, besides small velocity effects. Therefore, the presented theory is inadequate.

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 results on the evolution of the intrinsic scatter of black hole masses considering different implementations of a model in which black holes only grow via mergers. We demonstrate how merger driven growth affects the correlations between black hole mass and host bulge mass. The simple case of an initially log-normal distributed scatter in black hole and bulge masses combined with random merging within the galaxy population results in a decreasing scatter with merging generation/number as predicted by the Central-limit theorem. In general we find that the decrease in scatter {\sigma} is well approximated by {\sigma}merg(m) = {\sigma}ini \times (m + 1)^(-a/2) with a = 0.42 for a range of mean number of mergers m < 50. For a large mean number of mergers (m > 100) we find a convergence to a = 0.61. This is valid for a wide range of different initial distributions, refill-scenarios or merger mass-ratios. Growth scenarios based on halo merger trees of a (100 Mpc)^3 dark matter LambdaCDM-simulation show a similar behaviour with a scatter decrease of a = 0.30 with typical number of mergers m < 50 consistent with random merging (best matching model: a = 0.34). Assuming a present day scatter of 0.3 dex in black hole mass and a mean number of mergers not exceeding m = 50 our results imply a scatter of 0.6 dex at z = 3 and thus a possible scenario in which overmassive (and undermassive) black holes at high redshift are a consequence of a larger intrinsic scatter in black hole mass. A simple toy model connecting the growth of black holes to the growth of LambdaCDM dark matter halos via mergers, neglecting any contribution from accretion, yields a consistent M\cdot -MBulge relation at z = 0 - if we assume the correct initial relation.

We use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.

We explore the gravitational instability of clumpy and turbulent gas discs, taking into account the Larson-type scaling laws observed in giant molecular clouds (GMCs) and HI, as well as more general scaling relations. This degree of freedom is of special interest in view of the coming high-z ISM surveys, and is thus potentially important for understanding the dynamical effects of turbulence at all epochs of galaxy evolution. Our analysis shows that turbulence has a deep impact on the gravitational instability of the disc. It excites a rich variety of stability regimes, several of which have no classical counterpart. Among other diagnostics, we provide two useful tools for observers and simulators: (1) the stability map of turbulence, which illustrates our stability scenario and relates it to the phenomenology of interstellar turbulence: GMC/HI observations, simulations and models; (2) a Toomre-like stability criterion, $Q\geq\bar{Q}$, which applies to a large class of clumpy/turbulent discs. We make specific predictions about GMC and cold-HI turbulence, and point out the implications of our analysis for high-z galaxy surveys.

Given the importance of shear flows for astrophysical gas dynamics, we study the evolution of the Kelvin-Helmholtz instability (KHI) analytically and numerically. We derive the dispersion relation for the two-dimensional KHI including viscous dissipation. The resulting expression for the growth rate is then used to estimate the intrinsic viscosity of four numerical schemes depending on code-specific as well as on physical parameters. Our set of numerical schemes includes the Tree-SPH code VINE, an alternative SPH formulation developed by Price (2008), and the finite-volume grid codes FLASH and PLUTO. In the first part, we explicitly demonstrate the effect of dissipation-inhibiting mechanisms such as the Balsara viscosity on the evolution of the KHI. With VINE, increasing density contrasts lead to a continuously increasing suppression of the KHI (with complete suppression from a contrast of 6:1 or higher). The alternative SPH formulation including an artificial thermal conductivity reproduces the analytically expected growth rates up to a density contrast of 10:1. The second part addresses the shear flow evolution with FLASH and PLUTO. Both codes result in a consistent non-viscous evolution (in the equal as well as in the different density case) in agreement with the analytical prediction. The viscous evolution studied with FLASH shows minor deviations from the analytical prediction.

A simple model of gas accretion in young galaxy disks suggests that fast turbulent motions can be driven by accretion energy for a time t_acc~2(epsilon^{0.5} GM^2/xi V^3)^{0.5} where epsilon is the fraction of the accretion energy going into disk turbulence, M and V are the galaxy mass and rotation speed, and xi is the accretion rate. After t_acc, accretion is replaced by disk instabilities as a source of turbulence driving, and shortly after that, energetic feedback by young stars should become important. The star formation rate equilibrates at the accretion rate after 1 to 2 t_acc, depending on the star formation efficiency per dynamical time. The fast turbulence that is observed in high redshift starburst disks is not likely to be driven by accretion because the initial t_acc phase is over by the time the starburst is present. However, the high turbulent speeds that must have been present earlier, when the observed massive clumps first formed, could have been driven by accretion energy. The combined observations of a high relative velocity dispersion in the gas of z~2 clumpy galaxies and a gas mass comparable to the stellar mass suggests that either the star formation efficiency is fairly high, perhaps 10x higher than in local galaxies, or the observed turbulence is powered by young stars.

Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.

High resolution observations with the NIR adaptive optics integral field spectrograph SINFONI at the VLT proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of three-dimensional high resolution hydrodynamical simulations with the Pluto code, we follow the evolution of such clusters, focusing on stellar mass loss. This leads to clumpy or filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. In order to capture the relevant physics in the inner region, we treat this disc separately by viscously evolving the radial surface density distribution. This enables us to link the tens of parsec scale region (accessible via SINFONI observations) to the (sub-)parsec scale region (observable with the MIDI instrument and via water maser emission). In this work, we concentrate on the effects of a parametrised turbulent viscosity to generate angular momentum and mass transfer in the disc and additionally take star formation into account. Input parameters are constrained by observations of the nearby Seyfert 2 galaxy NGC 1068. At the current age of its nuclear starburst of 250 Myr, our simulations yield disc sizes of the order of 0.8 to 0.9 pc, gas masses of 1.0e6 solar masses and mass transfer rates of 0.025 solar masses per year through the inner rim of the disc. This shows that our large scale torus model is able to approximately account for the disc size as inferred from interferometric observations in the mid-infrared and compares well to the extent and mass of a rotating disc structure as inferred from water maser observations. Several other observational constraints are discussed as well.

Recently, high resolution observations with the help of the near-infrared adaptive optics integral field spectrograph SINFONI at the VLT proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of high resolution hydrodynamical simulations with the PLUTO-code, we follow the evolution of such clusters, especially focusing on mass and energy feedback from young stars. This leads to a filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. Here, we concentrate on the long-term evolution of the nuclear disc in NGC 1068 with the help of an effective viscous disc model, using the mass input from the large scale simulations and accounting for star formation in the disc. This two-stage modelling enables us to connect the tens of parsec scale region (observable with SINFONI) with the parsec scale environment (MIDI observations). At the current age of the nuclear star cluster, our simulations predict disc sizes of the order of 0.8 to 0.9 pc, gas masses of 1.0e6 Msun and mass transfer rates through the inner boundary of 0.025 Msun/yr in good agreement with values derived from observations.

Our adaptive optics observations of nearby AGN at spatial resolutions as small as 0.085arcsec show strong evidence for recent, but no longer active, nuclear star formation. We begin by describing observations that highlight two contrasting methods by which gas can flow into the central tens of parsecs. Gas accumulation in this region will inevitably lead to a starburst, and we discuss the evidence for such events. We then turn to the impact of stellar evolution on the further inflow of gas by combining a phenomenological approach with analytical modelling and hydrodynamic simulations. These complementary perspectives paint a picture in which all the processes are ultimately regulated by the mass accretion rate into the central hundred parsecs, and the ensuing starburst that occurs there. The resulting supernovae delay accretion by generating a starburst wind, which leaves behind a clumpy interstellar medium. This provides an ideal environment for slower stellar outflows to accrete inwards and form a dense turbulent disk on scales of a few parsecs. Such a scenario may resolve the discrepancy between the larger scale structure seen with adaptive optics and the small scale structure seen with VLTI.

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.

We investigate the collapse and fragmentation of low-mass, trans-sonically turbulent prestellar cores, using SPH simulations. The initial conditions are slightly supercritical Bonnor-Ebert spheres, all with the same density profile, the same mass (M_O=6.1 Msun) and the same radius (R_O=17,000 AU), but having different initial turbulent velocity fields. Four hundred turbulent velocity fields have been generated, all scaled so that the mean Mach number is M=1. Then a subset of these, having a range of net angular momenta, j, has been evolved. The evolution of these turbulent cores is not strongly correlated with j. Instead it is moderated by the formation of filamentary structures due to converging turbulent flows. A high fraction (~ 82%) of the protostars forming from turbulent cores are attended by protostellar accretion discs, but only a very small fraction (~16%) of these discs is sufficiently cool and extended to develop non-linear gravitational instabilities and fragment.

Galactic disk formation requires knowledge about the initial conditions under which disk galaxies form, the boundary conditions that affect their secular evolution and the micro-physical processes that drive the multi-phase interstellar medium and regulate their star formation history. Most of these ingredients are still poorly understood. Recent high-resolution observations of young high-redshift disk galaxies provide insight into early phases of galactic disk formation and evolution. Combined with low-redshift disk data these observations should eventually allow us to reconstruct the origin and evolution of late-type galaxies. I summarize some of the major problems that need to be addressed for a more consistent picture of galactic disk formation and evolution.

Abridged: We use three-dimensional SPH simulations to investigate the collapse of low-mass prestellar cores and the formation and early evolution of protostellar discs. The initial conditions are slightly supercritical Bonnor-Ebert spheres in rigid rotation. The core mass and initial radius are held fixed at M_O=6.1 M_sun and R_O=17,000 AU, and the only parameter that we vary is the initial angular speed \Omega_O. Protostellar discs forming from cores with \Omega_O<1.35 10d-13 1/s have radii between 100 and 300 AU and are quite centrally concentrated; due to heating by gas infall onto the disc and accretion onto the central object, they are also quite warm, T>100 K, and therefore stable against gravitational fragmentation. In contrast, more rapidly rotating cores form discs which are less concentrated and cooler, and have radii between 400 and 1000 AU; as a consequence they are prone to gravitational fragmentation and the formation of multiple systems. We derive a criterion that predicts whether a rigidly rotating core having given M_O, R_O and \Omega_O will produce a protostellar disc which fragments whilst material is still infalling from the core envelope. We then apply this criterion to core samples for which M_O, R_O and \Omega_O have been estimated observationally. We conclude that the observed cores are stable against fragmentation at this stage, due to their low angular speeds and the heat delivered at the accretion shock where the infalling material hits the disc.

We present a Fortran 95 code for simulating the evolution of astrophysical systems using particles to represent the underlying fluid flow. The code is designed to be versatile, flexible and extensible, with modular options that can be selected either at compile time or at run time. We include a number of general purpose modules describing a variety of physical processes commonly required in the astrophysical community. The code can be used as an N-body code to evolve a set of particles in two or three dimensions using either a Leapfrog or Runge-Kutta-Fehlberg integrator, with or without individual timesteps for each particle. Particles may interact gravitationally as $N$-body particles, and all or any subset may also interact hydrodynamically, using the Smoothed Particle Hydrodynamic (SPH) method. Massive point particles (`stars') which may accrete nearby SPH or $N$-body particles may also be included. The default free boundary conditions can be replaced by a module to include periodic boundaries. Cosmological expansion may also be included. An interface with special purpose `GRAPE' hardware may also be selected. If available, forces obtained from the GRAPE coprocessors may be transparently substituted for those obtained from the default tree based calculation. The code may be run without modification on single processors or in parallel using OpenMP compiler directives on large scale, shared memory parallel machines. In comparison to the Gadget-2 code of Springel 2005, the gravitational force calculation is $\approx 3.5 - 4.8$ times faster with VINE when run on 8 Itanium~2 processors in an SGI Altix, while producing nearly identical outcomes in our test problems. We present simulations of several test problems, including a merger simulation of two elliptical galaxies with 800000 particles.

We report the results from a deep HST NICMOS H-band imaging survey of 28 z < 0.3 QSOs from the Palomar-Green (PG) sample. This program is part of QUEST (Quasar / ULIRG Evolution STudy) and complements a similar set of data on 26 highly-nucleated ULIRGs presented in Paper I. Our analysis indicates that the fraction of QSOs with elliptical hosts is higher among QSOs with undetected far-infrared (FIR) emission, small infrared excess, and luminous hosts. The hosts of FIR-faint QSOs show a tendency to have less pronounced merger-induced morphological anomalies and larger QSO-to-host luminosity ratios on average than the hosts of FIR-bright QSOs, consistent with late-merger evolution from FIR-bright to FIR-faint QSOs. The spheroid sizes and total host luminosities of the radio-quiet PG QSOs in our sample are statistically indistinguishable from the ULIRG hosts presented in Paper I, while those of radio-loud PG QSOs are systematically larger and more luminous. ULIRGs and PG QSOs with elliptical hosts fall near, but not exactly on, the fundamental plane of inactive spheroids. We confirm the systematic trend noted in Paper I for objects with small (< 2 kpc) spheroids to be up to ~1 mag. brighter than inactive spheroids. The host colors and wavelength dependence of their sizes support the idea that these deviations are due at least in part to non-nuclear star formation. However, the amplitudes of these deviations does not depend on host R-H colors. Taken at face value (i.e., no correction for extinction or the presence of a young stellar population), the H-band spheroid-host luminosities imply BH masses ~5 -- 200 x 10^7 M_sun and sub-Eddington mass accretion rates for both QSOs and ULIRGs. These results are compared with published BH mass estimates derived from other methods. (abridged)

We use high resolution N-body/SPH simulations to study the hydrodynamical interaction between the Large Magellanic Cloud (LMC) and the hot halo of the Milky Way. We investigate whether ram-pressure acting on the satellite's ISM can explain the peculiarities observed in the HI distribution and the location of the recent star formation activity. Due to the present nearly edge-on orientation of the disk with respect to the orbital motion, compression at the leading edge can explain the high density region observed in HI at the south-east border. In the case of a face-on disk (according to Mastropietro et al. 2008 the LMC was moving almost face-on before the last perigalactic passage), ram-pressure directed perpendicularly to the disk produces a clumpy structure characterized by voids and high density filaments that resemble those observed by the Parkes HI survey. As a consequence of the very recent edge-on motion, the H-alpha emission is mainly concentrated on the eastern side where 30 Doradus and most of the supergiant shells are located, although some H-alpha complexes form a patchy distribution on the entire disk. In this scenario only the youngest stellar complexes show a progression in age along the leading border of the disk.

We use the bulge Sersic index n and bulge-to-total ratio (B/T) to explore the fundamental question of how bulges form. We perform 2D bulge-disk-bar decomposition on H-band images of 143 bright, high stellar mass (>1.0e10 solar masses) low-to-moderately inclined (i<70 degrees) spirals. Our results are: (1) Our H-band bar fraction (~58%) is consistent with that from ellipse fits. (2) 70% of the stellar mass is in disks, 10% in bars, and 20% in bulges. (3) A large fraction (~69%) of bright spirals have B/T<0.2, and ~76% have low n<2 bulges. These bulges exist in barred and unbarred galaxies across a wide range of Hubble types. (4) About 65% (68%) of bright spirals with n<2 (B/T<0.2) bulges host bars, suggesting a possible link between bars and bulges. (5) We compare the results with predictions from a set of LCDM models. In the models, a high mass spiral can have a bulge with a present-day low B/T<0.2 only if it did not undergo a major merger since z<2. The predicted fraction (~1.6%) of high mass spirals, which have undergone a major merger since z<4 and host a bulge with a present-day low B/T<0.2, is a factor of over thirty smaller than the observed fraction (~66%) of high mass spirals with B/T<0.2. Thus, contrary to common perception, bulges built via major mergers since z<4 seriously fail to account for the bulges present in ~66% of high mass spirals. Most of these present-day low B/T<0.2 bulges are likely to have been built by a combination of minor mergers and/or secular processes since z<4.

We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 from Halpha line emission. The galaxies are selected from the larger sample of the SINS survey, based on the prominence of ordered rotational motions with respect to more complex merger induced dynamics. The quality of the data allows us to carefully select systems with kinematics dominated by rotation, and to model the gas dynamics across the whole galaxy using suitable exponential disk models. We obtain a good correlation between the dynamical mass and the stellar mass, finding that large gas fractions Mgas~M*) are required to explain the difference between the two quantities. We use the derived stellar mass and maximum rotational velocity Vmax from the modeling to construct for the first time the stellar mass Tully-Fisher relation at z ~ 2.2. The relation obtained shows a slope similar to what is observed at lower redshift, but we detect an evolution of the zero point. We find that at z ~ 2.2 there is an offset in log(M*) for a given rotational velocity of 0.41+-0.11 with respect to the local Universe. This result is consistent with the predictions of the latest N-body/hydrodynamical simulations of disk formation and evolution, which invoke gas accretion onto the forming disk in filaments and cooling flows. This scenario is in agreement with other dynamical evidence from SINS, where gas accretion from the halo is required to reproduce the observed properties of a large fraction of the z ~ 2 galaxies.

The fact that dark matter (DM), thus far, has revealed itself only on scales of galaxies and larger, again thrusts onto astrophysics the opportunity and the responsibility to confront the age old mystery "What is the nature of matter?" By deriving basic data on the nature of DM - e.g., mass of its particle(s), present mean temperature, distribution in galaxies and other structures in the universe, and capacity for dissipational collapse - we will be uncovering the properties of the dominant species of matter in the universe and significantly extending the standard models of particle physics. Determining the mass of the DM particle to an order of magnitude would help to sort out the particle family to which it (or they) belongs. Beyond mass, there are issues of stability. The DM particle may be unstable with a measurable half-life, or it may become unstable after absorbing a certain amount of energy from collisions. In both cases it would contribute to the present hot dark matter component. Some key parameters of DM can most accurately be measured in the very nearby universe because DM dominates the mass in the outer Milky Way (MW), in other galaxies in the Local Group, and in the Local Group in its entirety. The presence and amount of DM can be quantified by study of dynamical processes observable in fine detail within these entities. Precise measurements of 3-D velocities for stars, coherent star streams, and stars in satellite stellar systems out to the edge of the Galaxy can reveal "what is the shape, orientation, density law, and lumpiness of the dark matter halo" as well as "what is the total mass of the Galaxy?"

The concordance Lambda Cold Dark Matter (Lambda-CDM) model for the formation of structure in the Universe, while remarkably successful at describing observations of structure on large scales, continues to be challenged by observations on galactic scales. Fortunately, CDM models and their various proposed alternatives make a rich variety of testable predictions that make the Local Group and our own Milky Way Galaxy key laboratories for exploring dark matter (DM) in this regime. However, some of the most definitive tests of local DM require microarcsecond astrometry of faint sources, an astrometric regime that is a unique niche of SIM Lite. This chapter explores the important and distinct contributions that can be made by SIM Lite in the exploration of galaxy dynamics and DM on galaxy scales and that have cosmological consequences. Key areas of potential SIM Lite exploration include (1) measuring the shape, orientation, density law, and lumpiness of the dark halo of the Milky Way and other nearby galaxies, (2) determining the orbits of Galactic satellites, which may be representatives of late infall from the hierarchical formation of the Milky Way, (3) ascertaining the distribution of angular momentum and orbital anisotropy of stars and globular clusters to the outer reaches of the Galactic halo, dynamical properties that hold clues to the early hierarchical formation of the Galaxy, (4) measuring the physical nature of DM by placing strong constraints on the phase space density in the cores of nearby dSph galaxies, and (5) reconstructing the dynamical history of the Local Group through the determination of orbits and masses of its constituent galaxies.

We propose to advance investigations of electromagnetic radiation originating in atomic nuclei beyond its current infancy to a true astronomy. This nuclear emission is independent from conditions of gas, thus complements more traditional stronomical methods used to probe the nearby universe. Radioactive gamma-rays arise from isotopes which are made in specific locations inside massive stars, their decay in interstellar space traces an otherwise not directly observable hot and tenuous phase of the ISM, which is crucial for feedback from massive stars. Its intrinsic clocks can measure characteristic times of processes within the ISM. Frontier questions that can be addressed with studies in this field are the complex interiors of massive stars and supernovae which are key agents in galactic dynamics and chemical evolution, the history of star-forming and supernova activity affecting our solar-system environment, and explorations of occulted and inaccessible regions of young stellar nurseries in our Galaxy.

We present high resolution simulations on the impact of ionizing radiation of massive O-stars on the surrounding turbulent interstellar medium (ISM). The simulations are performed with the newly developed software iVINE which combines ionization with smoothed particle hydrodynamics (SPH) and gravitational forces. We show that radiation from hot stars penetrates the ISM, efficiently heats cold low density gas and amplifies over-densities seeded by the initial turbulence. The formation of observed pillar-like structures in star forming regions (e.g. in M16) can be explained by this scenario. At the tip of the pillars gravitational collapse can be induced, eventually leading to the formation of low mass stars. Detailed analysis of the evolution of the turbulent spectra shows that UV-radiation of O-stars indeed provides an excellent mechanism to sustain and even drive turbulence in the parental molecular cloud.

Detailed models of galactic disk formation and evolution require knowledge about the initial conditions under which disk galaxies form, the boundary conditions that affect their secular evolution and the micro-physical processes that drive the multi-phase interstellar medium and regulate the star formation history. Unfortunately, up to now, most of these ingredients are still poorly understood. The challenge therefore is to, despite this caveat, construct realistic models of galactic disks with predictive power. This short review will summarize some problems related to numerical simulations of galactic disk formation and evolution.

We use oblate axisymmetric dynamical models including dark halos to determine the orbital structure of intermediate mass to massive Coma early-type galaxies. We find a large variety of orbital compositions. Averaged over all sample galaxies the unordered stellar kinetic energy in the azimuthal and the radial direction are of the same order, but they can differ by up to 40 percent in individual systems. In contrast, both for rotating and non-rotating galaxies the vertical kinetic energy is on average smaller than in the other two directions. This implies that even most of the rotating ellipticals are flattened by an anisotropy in the stellar velocity dispersions. Using three-integral axisymmetric toy models we show that flattening by stellar anisotropy maximises the entropy for a given density distribution. Collisionless disk merger remnants are radially anisotropic. The apparent lack of strong radial anisotropy in observed early-type galaxies implies that they may not have formed from mergers of disks unless the influence of dissipational processes was significant.

We present a three-dimensional, fully parallelized, efficient implementation of ionizing UV radiation for smoothed particle hydrodynamics (SPH) including self-gravity. Our method is based on the SPH/tree code VINE. We therefore call it iVINE (for Ionization + VINE). This approach allows detailed high-resolution studies of the effects of ionizing radiation from e.g. young massive stars on their turbulent parental molecular clouds. In this paper we describe the concept and the numerical implementation of the radiative transfer for a plain-parallel geometry and we discuss several test cases demonstrating the efficiency and accuracy of the new method. As a first application, we study the radiatively driven implosion of marginally stable molecular clouds at various distances of a strong UV source and show that they are driven into gravitational collapse. The resulting cores are very compact and dense exactly as it is observed in clustered environments. Our simulations indicate that the time of triggered collapse depends on the distance of the core from the UV source. Clouds closer to the source collapse several $10^5$ years earlier than more distant clouds. This effect can explain the observed age spread in OB associations where stars closer to the source are found to be younger. We discuss possible uncertainties in the observational derivation of shock front velocities due to early stripping of proto-stellar envelopes by ionizing radiation.

We use two-dimensional kinematic maps of simulated binary disc mergers to investigate the $\lambda_R$-parameter, which is a luminosity weighted measure of projected angular momentum per unit mass. This parameter was introduced to subdivide the SAURON sample of early type galaxies in so called fast $\lambda_R > 0.1$ and slow rotators $\lambda_R < 0.1$. Tests on merger remnants reveal that $\lambda_R$ is a robust indicator of the true angular momentum content in elliptical galaxies. We find the same range of $\lambda_R$ values in our merger remnants as in the SAURON galaxies. The merger mass ratio is decisive in creating a slow or a fast rotator in a single binary merger, the former being created mostly in an equal mass merger. Slow rotators have a $\lambda_R$ which does not vary with projection. The confusion rate with face-on fast rotators is very small. Merger with low gas fractions form slow rotators with smaller ellipticities and are in much better agreement with the SAURON slow rotators. Remergers of merger remnants are slow rotators but tend to have too high ellipticities. Fast rotators maintain the angular momentum content from the progenitor disc galaxy if merger mass ratio is high. Some SAURON galaxies have values of $\lambda_R$ as high as our progenitor disc galaxies.

We describe the implementation of slitless radial velocity measurements of extragalactic planetary nebulae (PNs) with the 8.2 m Subaru telescope and its Cassegrain imaging spectrograph, FOCAS. As a first application, we have extended a previous search for PNs in NGC 4697 to larger angular distances from its center. A total of 218 PNs were detected, and their radial velocities were measured. We have added 56 new PN detections to the existing sample of 535, observed previously with the ESO VLT + FORS imaging spectrograph; 36 of these new 56 PNs are located at angular distances larger than 230 arcsec from the center of NGC 4697. We compare the new FOCAS velocities with the earlier FORS velocities, for 158 of the 162 reobserved sources, finding good agreement. We now have kinematic information extending out to 5 effective radii from the center. The outer line-of-sight velocity dispersion is a bit lower than estimated earlier. This result is compatible with the existence of a dark matter halo plus some degree of radial anisotropy, but the dark matter halo is rather inconspicuous, and it is still unclear how massive it can be. A more detailed global dynamical study of the whole set of PN velocities will be required to decide if they permit to narrow down the range of possible dark matter distributions in NGC 4697. The new radial velocities reveal no evidence of rotation at 5 effective radii.

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)

Dense, small molecular cloud cores have been identified as the direct progenitors of stars. One of the best studied examples is Barnard 68 which is considered a prototype stable, spherical gas core, confined by a diffuse high-pressure environment. Observations of its radial density structure however indicate that Barnard 68 should be gravitationally unstable and collapsing which appears to be inconsistent with its inferred long lifetime and stability. We argue that Barnard 68 is currently experiencing a fatal collision with another small core which will lead to gravitational collapse. Despite the fact that this system is still in an early phase of interaction, our numerical simulations imply that the future gravitational collapse is already detectable in the outer surface density structure of the globule which mimicks the profile of a gravitationally unstable Bonnor-Ebert sphere. Within the next 200,000 years Barnard 68 will condense into a low-mass solar-type star(s), formed in isolation, and surrounded by diffuse, hot interstellar gas. As witnessed in situ for Barnard 68, core mergers might in general play an important role in triggering star formation and shaping the molecular core mass distribution and by that also the stellar initial mass function.

We examine the origin of photometrical properties of N-body merger remnants from the perpective of their orbital content. We show that disc mergers alone are unlikely to form a boxy-discy dichotomy in isophotal shape, because of the survival of disc-like orbits even in violent mergers. The shape of the different orbit families can vary strongly, depending on how violent the merger was. Minor axis tubes can become boxy and box orbits can be round. However, for edge-on projections isophotal shape, orbital content and line-of-sight velocity are well connected.

We present Ha integral field spectroscopy of well resolved, UV/optically selected z~2 star-forming galaxies as part of the SINS survey with SINFONI on the ESO VLT. Our laser guide star adaptive optics and good seeing data show the presence of turbulent rotating star forming rings/disks, plus central bulge/inner disk components, whose mass fractions relative to total dynamical mass appears to scale with [NII]/Ha flux ratio and star formation age. We propose that the buildup of the central disks and bulges of massive galaxies at z~2 can be driven by the early secular evolution of gas-rich proto-disks. High redshift disks exhibit large random motions. This turbulence may in part be stirred up by the release of gravitational energy in the rapid cold accretion flows along the filaments of the cosmic web. As a result dynamical friction and viscous processes proceed on a time scale of <1 Gyr, at least an order of magnitude faster than in z~0 disk galaxies. Early secular evolution thus drives gas and stars into the central regions and can build up exponential disks and massive bulges, even without major mergers. Secular evolution along with increased efficiency of star formation at high surface densities may also help to account for the short time scales of the stellar buildup observed in massive galaxies at z~2.

The observed rapid onset of star formation in molecular clouds requires rapid formation of dense fragments which can collapse individually before being overtaken by global gravitationally-driven flows. Many previous investigations have suggested that supersonic turbulence produces the necessary fragmentation, without addressing however the source of this turbulence. Motivated by our previous (numerical) work on the flow-driven formation of molecular clouds, we investigate the expected timescales of the dynamical and thermal instabilities leading to the rapid fragmentation of gas swept up by large-scale flows, and compare them with global gravitational collapse timescales. We identify parameter regimes in gas density, temperature and spatial scale within which a given instability will dominate. We find that dynamical instabilities disrupt large-scale coherent flows via generation of turbulence, while strong thermal fragmentation amplifies the resulting low-amplitude density perturbations, thus leading to small-scale, high-density fragments as seeds for {\em local} gravity to act upon. Global gravity dominates only on the largest scales; large-scale gravitationally-driven flows promote the formation of groups and clusters of stars formed by turbulence, thermal fragmentation, and rapid cooling.

We present high resolution N-body/SPH simulations of the interacting cluster 1E0657-56. The main and the sub-cluster are modeled using extended cuspy LCDM dark matter halos and isothermal beta-profiles for the collisional component. The hot gas is initially in hydrostatic equilibrium inside the global potential of the clusters. We investigate the X-ray morphology and derive the most likely impact parameters, mass ratios and initial relative velocities. We find that the observed displacement between the X-ray peaks and the associated mass distribution, the morphology of the bow shock, the surface brightness and projected temperature profiles across the shock discontinuity can be well reproduced by offset 1:6 encounters where the sub-cluster has initial velocity (in the rest frame of the main cluster) close to 2 times the virial velocity of the main cluster dark matter halo. A model with the same mass ratio and lower velocity (1.5 times the main cluster virial velocity) matches quite well most of the observations. However, it does not reproduce the morphology of the main cluster peak. Dynamical friction strongly affects the kinematics of the sub-cluster so that the low velocity bullet is actually bound to the main system at the end of the simulation. We find that a relatively high concentration (c=6) of the main cluster dark matter halo is necessary in order to prevent the disruption of the associated X-ray peak. For a selected sub-sample of runs we perform a detailed three dimensional analysis following the past, present and future evolution of the interacting systems. In particular, we investigate the kinematics of the gas and dark matter components as well as the changes in the density profiles and the motion of the system in the L_X-T diagram.

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 study numerically the formation of molecular clouds in large-scale colliding flows including self-gravity. The models emphasize the competition between the effects of gravity on global and local scales in an isolated cloud. Global gravity builds up large-scale filaments, while local gravity -- triggered by a combination of strong thermal and dynamical instabilities -- causes cores to form. The dynamical instabilities give rise to a local focusing of the colliding flows, facilitating the rapid formation of massive protostellar cores of a few 100 M$_\odot$. The forming clouds do not reach an equilibrium state, though the motions within the clouds appear comparable to ``virial''. The self-similar core mass distributions derived from models with and without self-gravity indicate that the core mass distribution is set very early on during the cloud formation process, predominantly by a combination of thermal and dynamical instabilities rather than by self-gravity.

We model mock observations of collisionless N-body disc-disc mergers with an axisymmetric orbit superposition program that has been used to model Coma ellipticals. The remnants sample representatively the shapes of disc-disc mergers including prolate, triaxial and oblate objects. The aim is to better understand how the assumption of axial symmetry affects reconstructed masses and stellar motions of systems which are intrinsically not axisymmetric, whether it leads to a bias and how such a potential bias can be recognised in models of real galaxies. The mass recovery at the half-light radius depends on viewing-angle and intrinsic shape: edge-on views allow to reconstruct total masses with an accuracy between 20% (triaxial/prolate remnants) and 3% (oblate remnant). Masses of highly flattened, face-on systems are underestimated by up to 50%. Deviations in local mass densities can be larger where remnants are strongly triaxial or prolate. Luminous M/L are sensitive to box orbits in the remnants. Box orbits cause the central kinematics to vary with viewing-angle. Reconstructed luminous M/L and central masses follow this variation. Luminous M/L are always underestimated (up to a factor of 2.5). Respective dark halos in the models can be overestimated by about the same amount, depending again on viewing angle. Reconstructed velocity anisotropies depend on viewing angle and on the orbital composition of the remnant. We construct N-body realisations of the Schwarzschild models to discuss chaotic orbits and the virial equilibrium in our models. Apparently flattened, rotating ellipticals of intermediate mass are likely close to both, axial symmetry and edge-on orientation. Our results imply that Schwarzschild models allow a reconstruction of their masses and stellar anisotropies with high accuracy. (abridged)

We investigate both analytically and numerically the motion of massless particles orbiting primary star in a close circular binary system with particular focus on the gas drag effects. These are the first calculations with particles ranging in size from 1m to 10km, which account for the presence of a tidally perturbed gaseous disk. We have found numerically that the radial mass transport by the tidal waves plays a crucial role in the orbital evolution of particles. Numerical results are confirmed analytical calculations that do not assume anything about origin of the radial flow in the disk. We demonstrate that the migration rate of a particle in a disk out of radial equilibrium is enhanced due to the enhanced mass flux of gas colliding with the particle and the migration is always directed inward regardless of the sign of the radial gas velocity. Within the framework of the perturbation theory we derive general, formulae for short-term variations of the particle semi-major axis, eccentricity and inclination in such disk. The formulae account for departures from axial symmetry by introducing effective components of the gas velocity. They agree with numerical results within several percent. We have also found in numerical simulations that the tidal waves introduce coherence in periastron longitude and eccentricity for particles on neighbouring orbits. The degree of the coherence depends on the particle size and on the distance from the primary star, being most prominent for particles with 10m radius. The results are important mainly in the context of planetesimal formation and, to a lesser degree, during the early planetesimal accretion stage.

We study the instantaneous virial balance of clumps and cores (CCs) in 3D simulations of driven, MHD, isothermal molecular clouds (MCs). The models represent a range of magnetic field strengths in MCs from subcritical to non-magnetic regimes. We identify CCs at different density thresholds, and for each object, we calculate all the terms that enter the Eulerian form of the virial theorem (EVT). A CC is considered gravitationally bound when the gravitational term in the EVT is larger than the amount for the system to be virialized, which is more stringent than the condition that it be large enough to make the total volume energy negative. We also calculate, quantities commonly used in the observations to indicate the state of gravitational boundedness of CCs such as the Jeans number J_c, the mass-to magnetic flux ratio mu_c, and the virial parameter alpha_vir. Our results show that: a) CCs are dynamical out-of-equilibrium structures. b) The surface energies are of the same order than their volume counterparts c) CCs are either in the process of being compressed or dispersed by the velocity field. Yet, not all CCs that have a compressive net kinetic energy are gravitationally bound. d) There is no 1-to-1 correspondence between the state of gravitational boundedness of a CC as described by the virial analysis or as implied by the classical indicators. In general, in the virial analysis, we observe that only the inner regions of the objects are gravitationally bound, whereas J_c, alpha_vir, and mu_c estimates tend to show that they are more gravitationally bound at the lowest threshold levels and more magnetically supercritical. g) We observe, in the non-magnetic simulation, the existence of a bound core with structural and dynamical properties that resemble those of the Bok globule Barnard 68 (B68).

We present a two-dimensional kinematic analysis for a sample of simulated binary disc merger remnants with mass ratios 1:1 and 3:1. For the progenitor discs we used pure stellar models as well as models with 10% of their mass in gas. A multitude of phenomena also observed in real galaxies are found in the simulations. These include misaligned rotation, embedded discs, gas rings, counter-rotating cores and kinematic misaligned discs. Using the 2D maps we illustrate projection effects and the change in properties of a merger remnant when gas is included in the merger. We find that kinematic peculiar subsystems are preferably formed in equal mass mergers. Equal-mass collisionless remnants can show almost no rotation, regular rotation or strong kinematic misalignment. The inclusion of gas makes the remnants appear more round(1:1) and axisymmetric(3:1). Counter-Rotating Cores (CRCs) are almost exclusively formed in equal-mass mergers with a dissipational component. 3:1 remnants show a much more regular structure. We quantify these properties by applying the kinemetric methods recently developed by Krajnovi\'c et al. This work will help to understand observations of elliptical galaxies with 2D field spectrographs, like SAURON.

The formation of tidal dwarf galaxies (TDG) inside tidal arms of interacting disk galaxies has been studied with N-body and N-body/SPH simulations at different resolutions. In pure N-body simulations no bound objects are formed at high resolution. At low resolution bound objects can form in tidal tails in agreement with previous work. We conclude that tidal dwarf galaxies are not likely to form by pure collisionless collapse in tidal tails. However, the presence of a sufficiently massive and extended gas component in the progenitor disk supports the formation of bound stellar objects in the tidal arms. Our results clearly favor a dissipation supported scenario in which the formation of TDGs is induced by the local collapse of gas which then triggers the collapse of the stellar component.

We revisit the analysis of the Non-linear Thin Shell Instability (NTSI) numerically, including magnetic fields. The magnetic tension force is expected to work against the main driver of the NTSI -- namely transverse momentum transport. However, depending on the field strength and orientation, the instability may grow. For fields aligned with the inflow, we find that the NTSI is suppressed only when the Alfv\'en speed surpasses the (supersonic) velocities generated along the collision interface. Even for fields perpendicular to the inflow, which are the most effective at preventing the NTSI from developing, internal structures form within the expanding slab interface, probably leading to fragmentation in the presence of self-gravity or thermal instabilities. High Reynolds numbers result in local turbulence within the perturbed slab, which in turn triggers reconnection and dissipation of the excess magnetic flux. We find that when the magnetic field is initially aligned with the flow, there exists a (weak) correlation between field strength and gas density. However, for transverse fields, this correlation essentially vanishes. In light of these results, our general conclusion is that instabilities are unlikely to be erased unless the magnetic energy in clouds is much larger than the turbulent energy. Finally, while our study is motivated by the scenario of molecular cloud formation in colliding flows, our results span a larger range of applicability, from supernovae shells to colliding stellar winds.

Cold clouds embedded in warm media are very common objects in astrophysics. Their disruption timescale depends strongly on the dynamical configuration. We discuss the evolution of an initially homogeneous cold cloud embedded in warm turbulent gas. Within a couple of dynamical timescales, the filling factor of the cold gas within the original cloud radius drops below 50%. Turbulent diffusivities estimated from the time evolution of radial filling factor profiles are not constant with time. Cold and warm gas are bodily transported by turbulence and mixed. This is only mildly indicated by column density maps. The radiation field within the cloud, however, increases by several orders of magnitudes due to the mixing, with possible consequences for cloud chemistry and evolution within a few dynamical timescales.

Recent observations with the Spitzer Space Telescope show clear evidence that star formation takes place in the surrounding of young massive O-type stars, which are shaping their environment due to their powerful radiation and stellar winds. In this work we investigate the effect of ionising radiation of massive stars on the ambient interstellar medium (ISM): In particular we want to examine whether the UV-radiation of O-type stars can lead to the observed pillar-like structures and can trigger star formation. We developed a new implementation, based on a parallel Smooth Particle Hydrodynamics code (called IVINE), that allows an efficient treatment of the effect of ionising radiation from massive stars on their turbulent gaseous environment. Here we present first results at very high resolution. We show that ionising radiation can trigger the collapse of an otherwise stable molecular cloud. The arising structures resemble observed structures (e.g. the pillars of creation in the Eagle Nebula (M16) or the Horsehead Nebula B33). Including the effect of gravitation we find small regions that can be identified as formation places of individual stars. We conclude that ionising radiation from massive stars alone can trigger substantial star formation in molecular clouds.

We suggest that the Orion A cloud is gravitationally collapsing on large scales, and is producing the Orion Nebula Cluster due to the focusing effects of gravity acting within a finite cloud geometry. In support of this suggestion, we show how an elliptical rotating sheet of gas with a modest density gradient along the major axis can collapse to produce a structure qualitatively resembling Orion A, with a fan-shaped structure at one end, ridges or filaments along the fan, and a narrow curved filament at the other end reminiscent of the famous integral-shaped filament. The model produces a local concentration of mass within the narrow filament which in principle could form a dense cluster of stars like that of the Orion Nebula. We suggest that global gravitational contraction might be a more common feature of molecular clouds than previously recognized, and that the formation of star clusters is a dynamic process resulting from the focusing effects of gravity acting upon the geometry of finite clouds.

The evidence for a shortage of exosolar planets with semimajor axes -1.1 <= log (a/AU) <= -0.2 is investigated. It is shown that this valley results from a gap in the radial distribution of planets, orbiting stars with masses M >= 1.2 Msolar (the high-mass sample, HMS). No underabundance is found for planets orbiting stars with smaller masses. The observational data also indicate that within the HMS population it is preferentially the more massive planets with M sin(i) >= 0.8 M_J that are missing. Monte-Carlo simulations of planet formation and migration are presented that reproduce the observed shortage of planets in the observed radius regime. A dependence on the disk depletion timescale tau_dep is found. The gap is more pronounced for tau_dep = 10^6 - 10^7 yrs than for tau_dep = 3*10^6 - 3*10^7 yrs. This might explain the observed trend with stellar mass if disks around stars with masses M* >= 1.2 Msolar have shorter depletion timescales than those around less massive stars. Possible reasons for such a dependence are a decrease of disk size and an increase of stellar EUV flux with stellar mass.

We present a large set of merger simulations of early-type disc galaxies with mass ratios of 1:1 and 3:1 and 10% of the total disc mass in gas. In contrast to the collisionless case equal-mass mergers with gas do not result in very boxy remnants which is caused by the suppression of box orbits and the change of the projected shape of minor-axis tube orbits in the more axisymmetric remnants. The isophotal shape of 3:1 remnants and the global kinematic properties of 1:1 and 3:1 remnants are only weakly affected by the presence of gas. 1:1 remnants are slowly rotating whereas 3:1 remnants are fast rotating and discy. The shape of the stellar LOSVD is strongly influenced by gas. The LOSVDs of collisionless remnants have broad leading wings while their gaseous counterparts show steep leading wings, more consistent with observations of elliptical galaxies. We show that this change is also caused by the suppressed populating of box orbits and it is amplified by the formation of extended gas discs in the merger remnants. If elliptical galaxies have formed from mergers our results indicate that massive, slowly rotating boxy elliptical galaxies can not have formed from dissipative mergers of discs. Pure stellar (dry) mergers are the more likely candidates. On the other hand lower mass, fast rotating and discy ellipticals can have formed from dissipative (wet) mergers of early-type discs. So far, only unequal-mass disc mergers with gas can successfully explain their observed substructure. This is consistent with the revised morphological classification scheme of increasing importance of gas dissipation when moving from boxy ellipticals to discy ellipticals and then to spiral galaxies, proposed by Kormendy & Bender (abbreviated).

We present results from our Very Large Telescope large program to study the dynamical evolution of local Ultraluminous Infrared Galaxies (ULIRGs) and QSOs. This paper is the second in a series presenting the stellar kinematics of 54 ULIRGs, derived from high resolution, long-slit H- and K-band spectroscopy. The data presented here, including observations of 17 new targets, are mainly focused on sources that have coalesced into a single nucleus. The stellar kinematics, extracted from the CO ro-vibrational bandheads in our spectra, indicate that ULIRG remnants are dynamically heated systems with a mean dispersion of 161 km/s. The combination of kinematic, structural, and photometric properties of the remnants indicate that they mostly originate from major mergers and that they result in the formation of systems supported by random motions, therefore, elliptical galaxies. The peak of the velocity dispersion distribution and the locus of ULIRGs on the fundamental plane of early-type galaxies indicate that the end products of ultraluminous mergers are typically moderate-mass ellipticals (of stellar mass ~10^10 - 10^11 M_sun). Converting the host dispersion into black hole mass with the aid of the M_BH-sigma relation yields black hole mass estimates of the order 10^7 - 10^8 M_sun and high accretion rates with Eddington efficiencies often >0.5.

Molecular Cloud Complexes (MCCs) are highly structured and ``turbulent''. Observational evidence suggests that MCCs are dynamically dominated systems, rather than quasi-equilibrium entities. The observed structure is more likely a consequence of the formation process rather than something that is imprinted after the formation of the MCC. Converging flows provide a natural mechanism to generate MCC structure. We present a detailed numerical analysis of this scenario. Our study addresses the evolution of a MCC from its birth in colliding atomic hydrogen flows up until the point when H$_2$ may begin to form. A combination of dynamical and thermal instabilities breaks up coherent flows efficiently, seeding the small-scale non-linear density perturbations necessary for local gravitational collapse and thus allowing (close to) instantaneous star formation. Many observed properties of MCCs come as a natural consequence of this formation scenario. Since converging flows are omnipresent in the ISM, we discuss the general applicability of this mechanism, from local star formation regions to galaxy mergers.

An overview is presented of the main properties of the interstellar medium. Evidence is summarized that the interstellar medium is highly turbulent, driven on different length scales by various energetic processes. Large-scale turbulence determines the formation of structures like filaments and shells in the diffuse interstellar medium. It also regulates the formation of dense, cold molecular clouds. Molecular clouds are now believed to be transient objects that form on timescales of order 1e7 yrs in regions where HI gas is compressed and cools. Supersonic turbulence in the compressed HI slab is generated by a combination of hydrodynamical instabilities, coupled with cooling. Turbulent dissipation is compensated by the kinetic energy input of the inflow. Molecular hydrogen eventually forms when the surface density in the slab reaches a threshold value of 1e21 cm^-2 at which point further cooling triggers the onset of star formation by gravitational collapse. A few Myrs later, the newly formed stars and resulting supernovae will disperse their molecular surrounding and generate new expanding shells that drive again turbulence in the diffuse gas and trigger the formation of a next generation of cold clouds. Although a consistent scenario of interstellar medium dynamics and star formation is emerging many details are still unclear and require more detailed work on microphysical processes as well as a better understanding of supersonic, compressible turbulence.

Almost all young stars are found in multiple systems. This suggests that protostellar cores almost always fragment into multiple objects. The observed properties of multiple systems such as their separation distribution and mass ratios provide strong constraints on star formation theories. We review the observed properties of young and old multiple systems and find that the multiplicity of stars changes. Such an evolution is probably due to (a) the dynamical decay of small-N systems and/or (b) the destruction of multiple systems within dense clusters. We review simulations of the fragmentation of rotating and turbulent molecular cores. Such models almost always produce multiple systems, however the properties of those systems do not match observations at all well. Magnetic fields appear to supress fragmentation, prehaps suggesting that they are not dynamically important in the formation of multiple systems. We finish by discussing possible reasons why theory fails to match observation, and the future prospects for this field.

The processes are investigated by which gas loses its angular momentum during the protogalactic collapse phase, leading to disk galaxies that are too compact with respect to the observations. High-resolution N-body/SPH simulations in a cosmological context are presented including cold gas and dark matter. A halo with quiet merging activity since z~3.8 and with a high spin parameter is analysed that should be an ideal candidate for the formation of an extended galactic disk. We show that the gas and the dark matter have similar specific angular momenta until a merger event occurs at z~2 with a mass ratio of 5:1. All the gas involved in the merger loses a substantial fraction of its specific angular momentum due to tidal torques and falls quickly into the center. Dynamical friction plays a minor role,in contrast to previous claims. In fact, after this event a new extended disk begins to form from gas that was not involved in the 5:1 merger event and that falls in subsequently. We argue that the angular momentum problem of disk galaxy formation is a merger problem: in cold dark matter cosmology substantial mergers with mass ratios of 1:1 to 6:1 are expected to occur in almost all galaxies. We suggest that energetic feedback processes could in principle solve this problem, however only if the heating occurs at the time or shortly before the last substantial merger event. Good candidates for such a coordinated feedback would be a merger-triggered star burst or central black hole heating. If a large fraction of the low angular momentum gas would be ejected as a result of these processes, late-type galaxies could form with a dominant extended disk component, resulting from late infall, a small bulge-to-disk ratio and a low baryon fraction, in agreement with observations.

The luminosity dependence of kinematical and isophotal properties of elliptical galaxies is investigated using numerical simulations of galaxy merging, combined with semi-analytical models of hierarchical structure formation. Mergers of spiral galaxies as the only formation mechanism for elliptical galaxies can neither reproduce the kinematical and photometric properties of very massive elliptical galaxies nor the change from rotationally flattened disky to anisotropic boxy systems with increasing luminosity. We present numerical simulations showing that binary mergers of early-type galaxies open an additional channel for the formation of anisotropic, slowly rotating and boxy ellipticals. Including this channel in a semi-analytical model we can successfully reproduce the observed trend that more luminous giant ellipticals are more boxy and less flattened by rotation. This trend can be strengthened by suppressing residual gas infall and star formation for galaxies with stellar bulge masses $M_* \geq 3 \times 10^{10} M_{\odot}$. Hence we propose that mergers of early-type galaxies play an important role for the assembly of massive elliptical galaxies.

We present first results from our Very Large Telescope large program to study the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs), which are the products of mergers of gas-rich galaxies. The full data set consists of high resolution, long-slit, H- and K-band spectra of 38 ULIRGs and 12 QSOs (between 0.042<z<0.268). In this paper, we present the sources that have not fully coalesced, and therefore have two distinct nuclei. This sub-sample consists of 21 ULIRGs, the nuclear separation of which varies between 1.6 and 23.3 kpc. From the CO bandheads that appear in our spectra, we extract the stellar velocity dispersion, sigma, and the rotational velocity, V_rot. The stellar dispersion equals 142 km/s on average, while V_rot is often of the same order. We combine our spectroscopic results with high-resolution infrared (IR) imaging data to study the conditions for ULIRG activity in interacting pairs. We find that the majority of ULIRGs are triggered by almost equal-mass major mergers of 1.5:1 average ratio. Less frequently, 3:1 encounters are also observed in our sample. However, less violent mergers of mass ratio >3:1 typically do not force enough gas into the center to generate ULIRG luminosities.

We present HI data of the dwarf galaxy DDO 47, aimed at testing the hypothesis that dark halo triaxiality might induce non-circular motions resulting in rotation curves best fitted by cored halos, even if the dark matter halo is intrinsically cuspy. We performed a harmonic decomposition of the velocity field in order to search for alleged non-circular motions needed to ``hide'' a cusp: in DDO 47 non-circular motions are globally at a level of 2-3 km s$^{-1}$, far from being sufficient to reconcile the observed rotation curve with the $\Lambda$CDM predictions. We conclude that the dark matter halo around DDO 47 is truly cored and that a cusp cannot be hidden by non-circular motions. More details are shown in Gentile et al. (2006, ApJL in press, astro-ph/0506538).

We present an analysis of high resolution HI data of the dwarf galaxy DDO 47, aimed at testing the hypothesis that dark halo triaxiality might induce non-circular motions resulting in rotation curves best fitted by cored halos, even if the dark matter halo is intrinsically cuspy. This hypothesis could be invoked in order to reconcile the predictions of the standard LambdaCDM theory with the rotation curves of disk galaxies. DDO 47 is an ideal case to test this hypothesis because it has a very regular velocity field, its rotation curve is best fitted by a cored halo and a NFW halo is inconsistent with the data. We analysed the velocity field through the higher-order harmonic terms in order to search for kinematical signatures of alleged non-circular motions needed to "hide" a cusp: the result is that globally non-circular motions are at a level of 2-3 km/s, and they are more likely to be associated to the presence of some spiral structure than to a global elongated potential (e.g. a triaxial halo). These non-circular motions are far from being sufficient to account for the discrepancy with the LambdaCDM predictions. We therefore conclude that the dark matter halo around the dwarf galaxy DDO 47 is truly cored and that a cusp cannot be hidden by non-circular motions.

We investigate the relationship between the velocity dispersion of the gas and the SN rate and feedback efficiency in the ISM. We explore the constancy of the velocity dispersion profiles in the outer parts of galactic disks at~6-8 km s^-1, and the transition to the starburst regime. Our results show that a) SN driving leads to constant velocity dispersions of sig~6 km s^-1 for the total gas and sigHI~3 km s^-1 for the HI gas, independent of the SN rate, for values of the rate between 0.01-0.5 the Galactic rate R_{G},b) the position of the transition to the starburst regime at SFR/Area~5*10^-3-10^-2 M_sol yr^-1 kpc^-2 observed in the simulations, is in good agreement with the transition to the starburst regime in the observations, c) for the high SN rates, no HI gas is present in the simulations box, however, for the total gas velocity dispersion, there is good agreement between the models and the observations,d) at the intermediate SN rates R/R_{G}~0.5-1, taking into account the thermal broadening of the HI line helps reach a good agreement in that regime between the models and the observations,e) for R/R_{G}<0.5, sig and sigHI fall below the observed values by a factor of~2. However, a set of simulation with different values of epsilon indicates that for larger values of the supernova feedback efficiencies, velocity dispersions of the HI gas of the order of 5-6 km s^{-1} can be obtained, in closer agreement with the observations. The fact that for R/R_{G}<0.5, the HI gas velocity dispersions are a factor ~2 smaller than the observed values could result from the fact that we might have underestimated the SN feedback efficiency. It might also be an indication that other physical processes couple to the stellar feedback in order to produce the observed level of turbulence in galactic disks.

We present results from our VLT Large Program to study the dynamical evolution of Ultraluminous Infrared Galaxies (ULIRGs). We have so far obtained near-infrared high-resolution ISAAC spectra of 53 local ULIRGs at several merger timescales and 12 Palomar-Green QSOs (more than half of which are IR-bright sources). We have extracted the stellar velocity dispersion and rotational velocity along our slits to derive the kinematics of the merging galaxies. These quantities enable us to answer the following questions about the evolution of ULIRGs: 1) What are the progenitor mass ratios?, 2) How do the stellar kinematics evolve?, and, 3) Is there a connection between ULIRGs and QSOs? We find that the Ultraluminous phase is mainly triggered by mergers of approximately equal mass galaxies, however, less violent minor mergers (of progenitor mass ratio ~3:1) also exist in our sample. Dynamical heating of the merging hosts is observed as the stellar systematic rotation decreases with time in favour of the increase of random motions. The merger remnants, being dispersion-dominated systems with non-negligible rotation, resemble elliptical galaxies. Placing ULIRGs on the fundamental plane of early-type galaxies shows that they resemble intermediate-mass ellipticals. After the nuclear coalescence, the black hole masses of ULIRGs, calculated from their relation to the host dispersions, are of the order 10^7-10^8 M_sun. To investigate whether ULIRGs go through a QSO phase during their evolution, we perform a similar (preliminary) analysis of the IR-bright-QSO kinematics. We find that the average dispersion of the IR-bright QSOs is similar to that of the ULIRG remnants, indicating that evolutionary links between the two populations may exist.

In order to specify cosmologically motivated initial conditions for major galaxy mergers (mass ratios $\leq$ 4:1) that are supposed to explain the formation of elliptical galaxies we study the orbital parameters of major mergers of cold dark matter halos using a high-resolution cosmological simulation. Almost half of all encounters are nearly parabolic with eccentricities $e \approx 1$ and no correlations between the halo spin planes or the orbital planes. The pericentric argument $\omega$ shows no correlation with the other orbital parameters and is distributed randomly. In addition we find that 50 % of typical pericenter distances are larger than half the halo's virial radii which is much larger than typically assumed in numerical simulations of galaxy mergers. In contrast to the usual assumption made in semi-analytic models of galaxy formation the circularities of major mergers are found to be not randomly distributed but to peak around a value of $\epsilon \approx 0.5$. Additionally all results are independent of the minimum progenitor mass and major merger definitions (i.e. mass ratios $\leq$ 4:1; 3:1; 2:1).

We use a set of twelve high-resolution N-body/hydrodynamical simulations in the $\Lambda$CDM cosmology to investigate the origin and formation rate of fossil groups (FGs), which are X-ray bright galaxy groups dominated by a large elliptical galaxy, with the second brightest galaxy being at least two magnitudes fainter. The simulations invoke star formation, chemical evolution with non-instantaneous recycling, metal dependent radiative cooling, strong star burst driven galactic super winds, effects of a meta-galactic UV field and full stellar population synthesis. We find an interesting correlation between the magnitude gap between the first and second brightest galaxy and the formation time of the group. It is found that FGs have assembled half of their final dark matter mass already at $z\ga1$, and subsequently typically grow by minor merging only, wheras non-FGs on average form later. The early assembly of FGs leaves sufficient time for galaxies of $L \sim L_*$ to merge into the central one by dynamical friction, resulting in the large magnitude gap at $z=0$. A fraction of 33$\pm$16% of the groups simulated are found to be fossil, whereas the observational estimate is $\sim$10-20%. The FGs are found to be X-ray over-luminous relative to non-FGs of the same optical luminosity, in qualitative agreement with observations. Finally, from a dynamical friction analysis is found that only because infall of $L \sim L_*$ galaxies happens along filaments with small impact parameters do FGs exist at all.

Molecular clouds (MCs) are highly structured and ``turbulent''. Colliding gas streams of atomic hydrogen have been suggested as a possible source of MCs, imprinting the filamentary structure as a consequence of dynamical and thermal instabilities. We present a 2D numerical analysis of MC formation via converging HI flows. Even with modest flow speeds and completely uniform inflows, non-linear density perturbations as possible precursors of MCs arise. Thus, we suggest that MCs are inevitably formed with substantial structure, e.g., strong density and velocity fluctuations, which provide the initial conditions for subsequent gravitational collapse and star formation in a variety of galactic and extragalactic environments.

We present a detailed investigation of the relation between the orbital content of merger remnants and observable properties of elliptical and S0 galaxies. Our analysis is based on the statistical sample of collisionless mergers of disk galaxies with different mass ratios and orbital parameters, published by Naab & Burkert. We use the spectral method by Carpintero & Aguilar to determine the orbital content of every remnant and correlate it with its intrinsic shape, and its projected kinematic and photometric properties. We discuss the influence of the bulge component and varying pericenter distances. The two most abundant orbit classes are the minor axis tubes and the box orbits. Their ratio seems to determine the basic properties of a remnant. On average, the fraction of minor axis tubes increases by a factor of two from a merger mass ratio of 1:1 to 4:1, whereas the fraction of box orbits decreases by 10%. At a given mass the central velocity dispersion of a remnant scales with the ratio of minor axis tubes to box orbits. Interestingly, the division line between rotational supported systems and pressure supported objects, $(v_{maj}/\sigma_0)^*=0.7$, turns out to coincide with a box to minor axis tube ratio of unity. The observed $h_3$-$v/\sigma$ anti-correlation for ellipticals can not be reproduced by collisionless merger remnants. We propose that this can only be reconciled by an additional physical process that significantly reduces the box orbit content. Remnants which are dominated by minor axis tube orbits have predominantly disky projections. Boxy remnants have always a box to minor axis tube ratio larger than one. This study will enable to identify observed ellipticals that could have formed, in the collisionless limit, by gas-poor disk mergers.

The evolution of the spin parameter of dark halos and the dependence on the halo merging history in a set of dissipationless cosmological LCDM simulations is investigated. Special focus is placed on the differences of the two commonly used versions of the spin parameter, namely lambda=J*E^1/2/(G*M^5/2) (Peebles 80) and lambda'=J/(sqrt(2)*M_vir*R_vir*V_vir) (Bullock et al. 01). Though the distribution of the spin transfer rate defined as the ratio of the spin parameters after and prior to a merger is similar to a high degree for both, lambda and lambda', we find considerable differences in the time evolution: while lambda' is roughly independent of redshift, lambda turns out to increase significantly with decreasing redshift. This distinct behaviour arises from small differences in the spin transfer during accretion events. The evolution of the spin parameter is strongly coupled with the virial ratio eta:=2*E_kin/|E_pot| of dark halos. Major mergers disturb halos and increase both their virial ratio and spin parameter for 1-2 Gyrs. At high redshifts (z=2-3) many halos are disturbed with an average virial ratio of eta = 1.3 which approaches unity until z=0. We find that the redshift evolution of the spin parameters is dominated by the huge number of minor mergers rather than the rare major merger events.

We suggest that large HI holes observed in the interstellar medium (ISM) of galaxies such as the Large Magellanic Cloud (LMC) and Holmberg II (Ho II, DDO 50, UGC 4305) can form as the combined result of turbulence coupled to thermal and gravitational instabilities. We investigate this problem with three dimensional hydrodynamical simulations, taking into account cooling and heating processes and the action of the self-gravity of the gas. We construct an algorithm for radiative transfer to post-process the simulated data and build emission maps in the 21 cm neutral hydrogen line. With this approach, we are able to reproduce the structure of the shells and holes as observed in regions of the ISM, where no stellar activity is detected. In order to quantify the comparison of our synthetic maps to the observations, we calculate the physical scale-autocorrelation length relation (L-L_{cr} relation) both on the synthetic HI maps and the HI map of Ho II. The L-L_{cr} relation shows a linear increase of the autocorrelation length with the physical scale up to the scale of energy injection and flattens for larger scales. The comparison of the L-L_{cr} relation between the observations and the synthetic maps suggests that turbulence is driven in the ISM of Ho II on large scales (~ 6 kpc). The slope of the L-L_{cr} relation in the linear regime in Ho II is better reproduced by models where turbulence is coupled with a weak efficiency cooling of the gas. These results demonstrate the importance of the interplay between turbulence and the thermodynamics of the gas for structure formation in the ISM. Our analysis can be used to determine the scale on which kinetic energy is injected in the ISM of dwarf irregular galaxies, and to derive, in a first approximation, the cooling rate of the gas.

The projected kinematical properties of unequal-mass merger remnants of disk galaxies are analysed and shown to agree well with observations of disky, fast rotating elliptical galaxies. This supports the major merger hypothesis of early-type galaxy formation. However, in contrast to previous claims, the merger remnants are very anisotropic with values of the anisotropy parameter that are similar to equal-mass merger remnants that form boxy, slowly rotating ellipticals. Including gas in the simulations does not change this result although the line-of-sight velocity profile and the intrinsic orbital structure are strongly affected by the presenceof gas. The kinematical difference between boxy and disky ellipticals appears not to be the amount of anisotropy but rather rotation and the shape of the velocity dispersion tensor. The apparent isotropy of observed disky ellipticals is shown to result from inclination effects. Even small inclination angles strongly reduce the measured anisotropy of fast rotating systems, seen in projection. A second problem is the limited amount of information that is available when measuring only the central velocity dispersion and a characteristic rotation and ellipticity. Methods are investigated that allow a better determination of the intrinsic anisotropy of fast rotating early-type galaxies with known inclination angles.

We investigate the dynamical state of a new class of extended star clusters, known as "Faint Fuzzies", that were discovered in two nearby S0 galaxies, NGC 1023 and NGC 3384. It is shown that the Faint Fuzzies of NGC 1023 lie in a fast rotating ring-like structure within the galactic disk with mean radius of 5 kpc, rotational velocity of 200 km/s and velocity dispersion of 115 km/s. We propose a scenario for the origin of Faint Fuzzies that is connected to the origin of S0 galaxies as a result of galaxy-galaxy interactions in dense environments. As is apparent in the Cartwheel galaxy, and is confirmed by numerical simulations, the passage of a small galaxy through, or close to, the center of a disk galaxy can form a ring of clumpy star formation with a radius comparable to the Faint Fuzzy ring radius in NGC 1023. In this case, the Faint Fuzzies are signposts for the transformation of spiral galaxies into lenticulars via such interactions.

Using semi-analytical models of galaxy formation, the origin of boxy and disky elliptical galaxies is investigated. We find that the simple scenario, motivated by N-body simulations, in which the isophotal shape is only dependent on the mass ratio of the last major merger, is not able to reproduce the observation that the fraction of boxy and disky ellipticals depends on galaxy luminosity. The observations can however be reproduce with the following reasonable assumptions: (i) equal-mass mergers lead to boxy ellipticals and unequal-mass mergers produce disky ellipticals (as motivated by N-body simulations) (ii) major mergers between bulge-dominated galaxies result always in boxy ellipticals, independent of the mass ratio, (iii) merger remnants that subsequently accrete gas leading to a secondary stellar disk with more than $20 %$ the total stellar fraction are always disky. Boxy and disky ellipticals can be divided into two subclasses, depending on their formation history. Boxy ellipticals either formed by equal mass mergers of disk galaxies or by major mergers of early-type galaxies. We find that disky ellipticals are indicators of unequal-mass mergers or late gas infall. Disky ellipticals with high luminosities preferentially are 1-component systems that result from unequal mass mergers whereas low-luminosity disky ellipticals are more likely to harbor secondary disk components. In addition, the fraction of disky ellipticals with secular disk components should increase in regions with higher galaxy densities. Taking into account the conversion of cuspy cores into flat low-density cores by black hole merging we find that disky ellipticals should contain central density cusps whereas boxy ellipticals should in general be characterized by flat cores.

We investigate the dynamical decay of non-hierarchical accreting triple systems and its implications on the ejection model as Brown Dwarf formation scenario. A modified chain-regularization scheme is used to integrate the equations of motion, that also allows for mass changes over time as well as for momentum transfer from the accreted gas mass onto the bodies. We integrate an ensemble of triple systems within a certain volume with different accretion rates, assuming several prescriptions of how momentum is transferred onto the bodies. We follow their evolution until the systems have decayed. We analyze the end states and decay times of these systems and determine the fraction of Brown Dwarfs formed, their escape speeds as well as the semi-major axis distribution of the formed Brown Dwarf binaries. We find that the formation probability of Brown Dwarfs depends strongly on the assumed momentum transfer which is related to the motion of the gas. Due to ongoing accretion and consequent shrinkage of the systems, the median escape velocity is increased by a factor of 2 and the binary separations are decreased by a factor of 5 compared with non-accreting systems. Furthermore, the obtained semi-major axis distribution drops off sharply to either side of the median, which is also supported by observations. We conclude that accretion and momentum transfer of accreted gas during the dynamical decay of triple systems is able to produce the observed distribution of close binary Brown Dwarfs, making the ejection model a viable option as Brown Dwarf formation scenario.

We analyze a three-dimensional smoothed particle hydrodynamics simulation of an evolving and later collapsing pre-stellar core. Using a three-dimensional continuum radiative transfer program, we generate images at 7 micron, 15 micron, 175 micron, and 1.3 mm for different evolutionary times and viewing angles. We discuss the observability of the properties of pre-stellar cores for the different wavelengths. For examples of non-symmetric fragments, it is shown that, misleadingly, the density profiles derived from a one-dimensional analysis of the corresponding images are consistent with one-dimensional core evolution models. We conclude that one-dimensional modeling based on column density interpretation of images does not produce reliable structural information and that multidimensional modeling is required.

We present two-dimensional simulations of finite, self-gravitating gaseous sheets. Unlike the case of infinite sheets, such configurations do not constitute equilibrium states but instead are subject to global collapse unless countered by pressure forces or rotation. The initial effect of finite geometry is to promote concentrations of material at the edges of the sheet. If the sheet is not perfectly circular, gravitational focussing results in enhanced concentrations of mass. In the second-most simple geometry, that of an elliptical outer boundary, the general result is collapse to a filamentary structure with the densest concentrations of mass at the ends of the filament. We suggest that these simple calculations have interesting implications for the gravitational evolution of overall molecular cloud structure, envisioning that such clouds might originate as roughly sheetlike sections of gas accumulated as a result of large-scale flows in the local interstellar medium. We show some examples of local clouds with overall filamentary shape and denser concentrations of mass and star clusters near the ends of the overall extended structure, suggestive of our simple ellipse collapse calculations. We suggest that cluster-forming gas is often concentrated as a result of gravity acting on irregular boundaries; this mechanism can result in very rapid infall of gas which may be of importance to the formation of massive stars. This picture suggests that much of the supersonic ``turbulence'' observed in molecular clouds might be gravitationally-generated. Our results may provide impetus for further theoretical explorations of global gravitational effects in molecular clouds and their implications for generating the substructure needed for fragmentation into stars and clusters.

The importance of angular momentum in regulating the sizes of galactic disks and by this their star formation history is highlighted. Tidal torques and accretion of satellites in principle provide enough angular momentum to form disks with sizes that are in agreement with observations. However three major problems have been identified that challenge cold dark matter theory and affect models of galaxy evolution: (1) too much angular momentum is transferred from the gas to the dark halos during infall, leading to disks with scale lengths that are too small, (2) bulgeless disks require more specific angular momentum than is generated cosmologically even if gas would not lose angular momentum during infall, (3) gravitational torques and hierarchical merging produce a specific angular momentum distribution that does not match the distribution required to form exponential disks naturally; some gas has exceptionally high angular momentum, leading to extended outer disks while another large gas fraction will contain very little specific angular momentum and is expected to fall into the galactic center, forming a massive and dominant bulge component. Any self-consistent theory of galaxy formation will require to provide solutions to these questions. Selective mass loss of low-angular-momentum gas in an early phase of galaxy evolution currently seems to be the most promising scenario. Such a process would have a strong affect on the early protogalactic evolution phase, the origin and evolution of galactic morphologies and link central properties of galaxies like the origin of central massive black holes with their global structure.

The specific angular momentum of Cold Dark Matter (CDM) halos in a $\Lambda$CDM universe is investigated. Their dimensionless specific angular momentum $\lambda'=\frac{j}{\sqrt{2}V_{vir} R{vir}}$ with $V_{vir}$ and $R_{vir}$ the virial velocity and virial radius, respectively depends strongly on their merging histories. We investigate a set of $\Lambda$CDM simulations and explore the specific angular momentum content of halos formed through various merging histories. Halos with a quiet merging history, dominated by minor mergers and accretion until the present epoch, acquire by tidal torques on average only 2% to 3% of the angular momentum required for their rotational support ($\lambda'=0.02$). This is in conflict with observational data for a sample of late-type bulgeless galaxies which indicates that those galaxies reside in dark halos with exceptionally high values of $\lambda' \approx 0.06-0.07$. Minor mergers and accretion preserve or slowly increase the specific angular momentum of dark halos with time. This mechanism is however not efficient enough in order to explain the observed spin values for late-type dwarf galaxies. Energetic feedback processes have been invoked to solve the problem that gas loses a large fraction of its specific angular momentum during infall. Under the assumption that dark halos hosting bulgeless galaxies acquire their mass via quiescent accretion, our results indicate yet another serious problem: the specific angular momentum gained during the formation of these objects is not large enough to explain their observed rotational properties,even if no angular momentum would be lost during gas infall.

The conditions that lead to self-regulated star formation, star bursts and the formation of massive stellar clusters are discussed. Massive stars have a strong impact on their environment, especially on the evolution of dwarf galaxies which are the building blocks of giant galaxies. Energy input by massive young clusters might help to solve some of the most important puzzles of galaxy formation: the cosmological substructure problem and the angular momentum problem.

The variation of the mass-to-light ratios M/L of early type galaxies as function of their luminosities L is investigated. It is shown that the tilt beta=0.27 (in the B--band) of the fundamental plane relation M/L ~ L^{beta} can be understood as a combination of two effects: about one-quarter (i.e. dbeta =0.07) is a result of systematic variations of the stellar population properties with increasing luminosity. The remaining three-quarters (i.e. dbeta =0.2) can be completely attributed to nonhomology effects that lead to a systematic change of the surface brightness profiles with increasing luminosity. Consequently, the observed tilt in the K-band (beta=0.17) where stellar population effects are negligible, is explained by nonhomology effects alone. After correcting for nonhomology, the mean value of the mass-to-light ratio of elliptical galaxies (M/L_B) is 7.1+-2.8 (1 sigma scatter).

The interpretation of microlensing results towards the Large Magellanic Cloud (LMC) still remains controversial. Whereas white dwarfs have been proposed to explain these results and, hence, to contribute significantly to the mass budget of our Galaxy, there are as well several constraints on the role played by white dwarfs. In this paper we analyze self-consistently and simultaneously four different results, namely, the local halo white dwarf luminosity function, the microlensing results reported by the MACHO team towards the LMC, the results of Hubble Deep Field (HDF) and the results of the EROS experiment, for several initial mass functions and halo ages. We find that the proposed log-normal initial mass functions do not contribute to solve the problem posed by the observed microlensing events and, moreover, they overproduce white dwarfs when compared to the results of the HDF and of the EROS survey. We also find that the contribution of hydrogen-rich white dwarfs to the dynamical mass of the halo of the Galaxy cannot be more than $\sim 4%$.

We consider the atmospheric flow on short-period extra-solar planets through two-dimensional numerical simulations of hydrodynamics with radiation transfer. One side is always exposed to the irradiation from the host star. The other is always in shadow. The temperature of the day side is determined by the equilibrium which the planetary atmosphere establishes with stellar radiation. Part of the thermal energy deposited on the day side is advected to the night side by a current. The radiation transfer, the night-side temperature distribution and by this the spectroscopic signature of the planet are sensitive functions of the atmospheric opacity. If the atmosphere contains grains with an abundanceand size distribution comparable to that of the interstellar medium, shallow heating occurs on the day side and the night side cools well below the day side. The temperature difference decreases as the abundance of grains is reduced. A simple analytic model of the dissipation of the circulation flow and associated kinetic heating is considered. This heating effect occurs mostly near the photosphere, not deep enough to significantly affect the size of planets. We show that the surface irradiation suppresses convection near the photospheric region on the day side. In some cases convection zones appear near the surface on the night side. This structural modification may influence the response and dissipation of tidal disturbances and alter the circularization and synchronization time scales.

We perform a large parameter survey of 112 collisionless N-body simulations of binary mergers of disk-galaxies with mass ratios of 1:1, 2:1, 3:1, and 4:1 using the special purpose hardware GRAPE. The basic properties of the remnants correlate with the mass ratio of the progenitor disks. We find that about 80% of the equal-mass merger simulations lead to slowly rotating remnants with significant minor-axis rotation. One half of all projected 1:1 remnants shows boxy and the other half shows disky isophotes. Remnants with mass ratios of 3:1 and 4:1 all rotate fast and show a small amount of minor-axis rotation. About 90% of the projected 3:1 and 4:1 remnants show disky isophotes. 2:1 remnants show intermediate properties. Projection effects lead to a large spread in the data in good agreement with observations. They do not change the fundamental kinematical differences between equal-mass and unequal-mass merger remnants. The results of this study weaken the disk merger scenario as the possible formation mechanism of massive boxy giant ellipticals. In general, the properties of equal-mass (and 2:1) merger remnants are consistent with the observed population of giant ellipticals in the intermediate mass regime between low mass fast rotating disky and bright massive boxy giant ellipticals. 3:1 and 4:1 merger remnants, however, are in very good agreement with the class of low luminosity, fast rotating giant elliptical galaxies. Binary mergers of disk galaxies are therefore still very good candidates for being the main formation mechanism for intermediate and low mass giant ellipticals. The homogeneous class of massive boxy ellipticals most likely formed by a different process.

The morphological properties of the Carina, Sculptor and Fornax dwarfs are investigated using new wide field data with a total area of 29 square degrees. The stellar density maps are derived, hinting that Sculptor possesses tidal tails indicating interaction with the Milky Way. Contrary to previous studies we cannot find any sign of breaks in the density profiles for the Carina and Fornax dwarfs. The possible existence of tidal tails in Sculptor and of King limiting radii in Fornax and Carina are used to derive global M/L ratios, without using kinematic data. By matching those M/L ratios to kinematically derived values we are able to constrain the orbital parameters of the three dwarfs. Fornax cannot have M/L smaller than 3 and must be close to its perigalacticon now. The other extreme is Sculptor that needs to be on an orbit with an eccentricity bigger than 0.5 to be able to form tidal tails despite its kinematic M/L.

Using the guidelines on SN element production provided by XMM-Newton, we summarize the results of ASCA observations on the element abundance in groups and clusters of galaxies. We show that while the metal production in groups could be described by a stellar population with a standard local IMF, clusters of galaxies require a more top-heavy IMF. We attribute an excess heavy element production to an IMF evolution with redshift. Dating the galaxy formation in clusters by observations of the star-formation rate, we conclude that the IMF variations have occurred preferentially at z>~4. We further combine our metallicity measurements with the mass function of clusters to estimate the role of clusters in the evolution of the metal content of the Universe. We argue that at no epoch stars are a major container of metals, unless groups of galaxies are not representative for the star-formation. This lends further support for the reduced (0.6 solar) mass-averaged oxygen abundance in the stellar population.

We summarize current models of the formation of spheroidal stellar systems. Whereas globular clusters form in an efficient mode of star formation inside turbulent molecular clouds, the origin of galactic spheroids, that is bulges, dwarf ellipticals, and giant ellipticals, is directly coupled with structure formation and merging of structures in the Universe. Disks are the fundamental building blocks of galaxies and the progenitors of galactic spheroids. The origin of the various types of spheroids and their global properties can be understood as a result of disk heating by external perturbations, internal disk instabilities, or minor and major mergers.

Based on theoretical arguments and quasi-stationary radiative MHD calculations, a model for an accretion-powered jet is presented. It is argued that accretion disks around BHs consist of 1) a cold, Keplerian-rotating and weakly magnetized mediumin the outer part, 2) a highly advective and turbulent-free plasma inside $r_\mm{tr}=10-20$ Schwarzschild radii, where magnetic fields are predominantly of large scale topology and in excess of thermal equipartition, and 3) an ion-dominated torus in the vicinity of the hole, where magnetic fields undergo a topological change into a monopole like-configuration. The action of magnetic fields interior to $r_\mm{tr}$ is to initiate torsional {\Alfven} waves that extract angular momentum from the disk-plasma and deposit it into the transition layer (-TL) above the disk, where the plasma is dissipative and tenuous. A significant fraction of the shear-generated toroidal magnetic field reconnects in the TL, thereby heating the plasma up to the virial-temperature and forming a super-Keplerian rotating, and hence centrifugally accelerated outflow. The strong magnetic field in the TL forces the electrons to cool rapidly yielding thereby an ion-dominated outflow. The toroidal magnetic field in the TL is in thermal equipartition with the ions, whereas the poloidal component is in equipartition with the electrons. Such a strong toroidal magnetic field is essential for increasing the jet-disk luminosity in the radio regime.

We use semi-analytical modelling techniques to investigate the progenitor morphologies of present day ellipticals. We find that, independent of the environment, the fraction of mergers of bulge dominated galaxies (early-types) increases with time. The last major merger of bright present day ellipticals with $M_{B} \lesssim -21$ is preferentially between bulge dominated galaxies, while those with $ M_{B} \sim -20$ have mainly experienced last major mergers between a bulge dominated and a disk dominated galaxy. Independent of specific model assumptions, more than 50% of present day elliptical in clusters with $M_B \lesssim -18$ had last major mergers which are not of spirals as usually expected within the standard merger scenario.

The formation of elliptical galaxies as a result of the merging of spiral galaxies is discussed. We analyse a large set of numerical N-Body merger simulations which show that major mergers can in principle explain the observed isophotal fine structure of ellipticals and its correlation with kinematical properties. Equal-mass mergers lead to boxy, slowly rotating systems, unequal-mass mergers produce fast rotating and disky ellipticals. However, several problems remain. Anisotropic equal mass mergers appear disky under certain projections which is not observed. The intrinsic ellipticities of remnants are often larger than observed. Finally, although unequal-mass mergers produce fast rotating ellipticals, the remnants are in general more anisotropic than expected from observations. Additional processes seem to play an important role which are not included in dissipationless mergers. Resolving these problems might provide interesting new information on the structure and gas content of the progenitors of early-type galaxies.

Self-interacting dark matter was proposed by Spergel & Steinhardt (2000) to alleviate two conflicts between Cold Dark Matter (CDM) models and observations. Firstly, CDM N-body simulations predict dark matter halo density profiles that diverge at the centre in disagreement with the constant density cores observed in late-type dwarf and Low Surface Brightness (LSB) galaxies. Secondly, N-body simulations predict an overabundance of subhalos in the Galactic halo. Using a simple semi-analytical argument we show that weakly self-interacting dark matter models, which can produce halo cores of the sizes observed in dark matter dominated galaxies, are unable to reconcile the number of satellites in the Galactic halo with the observed number of dwarf galaxies in the Local Group.

The interpretation of the old, cool white dwarfs recently found by Oppenheimer et al. (2001) is still controversial. Whereas these authors claim that they have finally found the elusive ancient halo white dwarf population that contributes significantly to the mass budget of the galactic halo, there have been several other contributions that argue that these white dwarfs are not genuine halo members but, instead, thick disk stars. We show here that the interpretation of this sample is based on the adopted distances, which are obtained from a color--magnitude calibration, and we demonstrate that when the correct distances are used a sizeable fraction of these putative halo white dwarfs belong indeed to the disk population. We also perform a maximum likelihood analysis of the remaining set of white dwarfs and we find that they most likely belong to the thick disk population. However, another possible explanation is that this sample of white dwarfs has been drawn from a 1:1 mixture of the halo and disk white dwarf populations.

Motivated by the presence of numerous dark matter clumps in the Milky Way's halo as expected from the cold dark matter cosmological model, we conduct numerical simulations to examine the heating of the disk. We construct an initial galaxy model in equilibrium, with a stable thin disk. The disk interacts with dark matter clumps for about 5 Gyr. Three physical effects are examined : first the mass spectrum of the dark matter clumps, second the initial thickness of the galactic disk, and third the spatial distribution of the clumps. We find that the massive end of the mass spectrum determines the amount of disk heating. Thicker disks suffer less heating. There is a certain thickness at which the heating owing to the interaction with the clumps becomes saturates. We also find that the heating produced by the model which mimics the distribution found in Standard CDM cosmology is significant and too high to explain the observational constraints. On the other hand, our model that corresponds to the clump distribution in a $\Lambda$CDM cosmology produces no significant heating. This result suggests that the $\Lambda$CDM cosmology is preferable with respect to the Standard CDM cosmology in explaining the thickness of the Milky Way.

We use high resolution numerical simulations to investigate the adiabatic contraction of dark matter halos with a Hernquist density profile. We test the response of the halos to the growth of additional axisymmetric disk potentials with various central concentrations and the spherically symmetric potential of a softened point mass. Adding the potentials on timescales that are long compared to the dynamical time scale of the halo, the contracted halos have density profiles that are in excellent agreement with analytical predictions based on the conservation of the adiabatic invariant $M(r)r$. This is surprising as this quantity is strictly conserved only for particles on circular orbits and in spherically symmetric potentials. If the same potentials are added on timescales that are short compared to the dynamical timescale, the result depends strongly on the adopted potential. The adiabatic approximation still works for disk potentials. It does, however, fail for the central potential.

Challenging stellar dynamical problems, such as the study of gravothermal oscillations in star clusters, have in the past initiated the very successful building of GRAPE special purpose computers. It is discussed, that present day tasks such as the formation and evolution of galactic nuclei with one or more massive black holes and the coupled stellar and gas dynamical processes in the formation of nuclei and star clusters, demand a new kind of hybrid architecture, using both GRAPE and a reconfigurable logics board called RACE. For such a system we have developed first implementations and floating point performance studies in the case of the SPH algorithm (smoothed particle hydrodynamics), which will be of great advantage for SPH modelling and also for direct $N$-body simulations using the more efficient Ahmad-Cohen neighbour schemes.

We present the first three-dimensional calculations of Radiation Driven Implosion of Molecular Clouds including self-gravity and ionization. We discuss the effects of initial density perturbations on the dynamics of ionizing globules and show that the onset of gravitational collapse can be significantly delayed for a multiple of the implosion timescale. We demonstrate that Radiation Driven Implosion could be an efficient process for injecting disordered kinetic energy into molecular clouds in the vicinity of massive stars.

The puzzling correlation between the spin parameter lambda of galactic disks and the disk-to-halo mass fraction fdisk is investigated. We show that such a correlation arises naturally from uncertainties in determining the virial masses of dark matter halos. This result leads to the conclusion that the halo properties derived from fits to observed rotation curves are still very uncertain which might explain part of the disagreements between cosmological models and observations. We analyse lambda and fdisk as function of the adopted halo virial mass. Reasonable halo concentrations require fdisk=0.01-0.07 which is significantly smaller than the universal baryon fraction. Most of the available gas either never settled into the galactic disks or was ejected subsequently. In both cases it is not very surprising that the specific angular momentum distribution of galactic disks does not agree with the cosmological predictions which neglect these effects.

We performed N-body/SPH simulations of merging gas-rich disk galaxies with mass ratios of 1:1 and 3:1. A stellar disk and bulge component and a dark halo was realized with collisionless particles, the gas was represented by SPH particles. Since we did not include star formation we focused on the gas dynamics and its influence on the structure of merger remnants. We find that equal-mass mergers are in general more effective in driving gas to the center. Around 50% of the gas resides in the central regions after the merger is complete. The gas in unequal-mass mergers keeps a large amount of its initial angular momentum and settles in a large scale disk while only 20% -30% of the gas is driven to the center. This general result is almost independent of the merger geometry. The gas in the outer regions accumulates in dense knots within tidal tails which could lead to the formation of open clusters or dwarf satellites. Later on, the gas knots loose angular momentum by dynamical friction and successively sink to the center of the remnant thereby increasing the total gas content of the disk. Due to the influence of the gas the simulated merger remnant becomes more oblate than its pure stellar counterpart. The shape of the stellar LOSVD changes and the third-order Gauss-Hermite coefficient $h_3$ is in good agreement with observations.

The one-point statistics of column density distributions of turbulent molecular cloud models are investigated and compared with observations. In agreement with the observations, the number N of pixels with surface density S is distributed exponentially N(S)=exp(-S/S0) in models of driven compressible supersonic turbulence. However, in contrast to the observations, the exponential slope defined by S0 is not universal but instead depends strongly on the adopted rms Mach number and on the smoothing of the data cube. We demonstrate that this problem can be solved if one restricts the analysis of the surface density distribution to subregions with sizes equal to the correlation length of the flow which turns out to be given by the driving scale. In this case, the column density distributions are universal with a slope that is in excellent agreement with the observations and independent of the Mach number or smoothing. The observed molecular clouds therefore are coherent structures with sizes of order their correlation lengths. Turbulence inside these clouds must be driven on the largest scales, if at all. Numerical models of turbulent molecular clouds have to be restricted to cubes with sizes similar to the correlation lengthscale in order to be compared with observations. Our results imply that turbulence is generated on scales that are much larger that the Jeans length. In this case, gravitational collapse cannot be suppressed and star formation should start in molecular clouds within a dynamical timescale after their formation.

We present a new correlation between the mass-to-light (M/L) ratio and the mean metallicity for the satellites of the Local Group. This relation together with their central surface brightness define a Fundamental Line where metal poor and low surface brightness dwarfs are dark matter dominated while metal rich high surface brightness systems will have a low M/L ratio. This Fundamental Line is independent of distance to the central galaxy (M31 and the Milky Way) and morphological type among other global parameters as their star formation history. The new M/L-metallicity relation indicates that dwarf spheroidal galaxies are dark matter confined. It can be interpreted if satellites experienced a continuous loss of metals with a final episode of gas ejection at the end of the star formation epoch. Only Globular Clusters are found to lie outside the Fundamental Line of dwarf satellites. Unlike satellite galaxies they are not dark matter confined.

We use semi-analytical modelling of galaxy formation to study the redshift evolution of the galaxy merger fractions and merger ratesin Lambda CDM and quintessence (QCDM) cosmologies, their dependence on physical parameters as the environment, the merger timescale, the way major mergers are defined, and the minimum mass of objects taken into account. We find that for a given final halo mass the redshift dependence of the merger fraction F_mg and the resulting merger rate can be fitted well by a power law for redshifts z <= 1. The normalization F_{mg}(0) and the slope m depend on the final halo mass. For a given merger timescale t_merg and an assumed maximum mass ratio R_major for major mergers, F_mg(0) and m depend exponentially on each other. The slope m depends logarithmically on the ratio of the final halo mass and the minimum halo mass taken into account. In addition, the local normalization F_{mg}(0) increases for larger R_{major} while m decreases. We compare the predicted merger fractionwith recent observations and find that the model cannot reproduce both the merger index and the normalization at the same time. In general the model underestimates F_{mg}(0) and m by a factor of 2.

We argue for implementing star formation on a viscous timescale in hydrodynamical simulations of disk galaxy formation and evolution. Modelling two-dimensional isolated disk galaxies with the Bhatnagar-Gross-Krook (BGK) hydrocode, we verify the analytic claim of various authors that if the characteristic timescale for star formation is equal to the viscous timescale in disks, the resulting stellar profile is exponential on several scale lengths whatever the initial gas and dark matter profile. This casts new light on both numerical and semi-analytical disk formation simulations which either (a) commence star formation in an already exponential gaseous disk, (b) begin a disk simulation with conditions known to lead to an exponential, i.e. the collapse of a spherically symmetric nearly uniform sphere of gas in solid body rotation under the assumption of specific angular momentum conservation, or (c) in simulations performed in a hierarchical context, tune their feedback processes to delay disk formation until the dark matter halos are slowly evolving and without much substructure so that the gas has the chance to collapse under conditions known to give exponentials. In such models, star formation follows a Schmidt-like law, which for lack of a suitable timescale, resorts to an efficiency parameter. With star formation prescribed on a viscous timescale however, we find gas and star fractions after $\sim$ 12 Gyr that are consistent with observations without having to invoke any ``fudge factor'' for star formation. Our results strongly suggest that despite our gap in understanding the exact link between star formation and viscosity, the viscous timescale is indeed the natural timescale for star formation.

The formation of massive stellar clusters in turbulent molecular clouds is investigated. We include artificial star formation and energy feedback of newly born stars. The obtained systems are not likely to survive. Case studies to determine conditions necessary for forming bound clusters will be done in the future.

Recent progress in the understanding of star formation is summarized. A consistent picture is emerging where molecular clouds form with turbulent velocity fields and clumpy substructure, imprinted already during their formation. The clouds are initially supported by supersonic turbulence which dissipates however within massive clumps on short timescales, of order their local dynamical timescales. As a result, the clumps collapse and fragment into stellar clusters. Subsequent energetic feedback processes by the newly formed stars will change the thermal and kinematical state of the surrounding molecular gas, triggering or suppressing star formation and leading eventually to the disruption of the molecular clouds. More work is required to understand these processes in greater details.

The formation of the Milky Way is discussed within the context of the cold dark matter scenario. Several problems arise which can be solved if the Galaxy experienced an early phase of gas heating and decoupling from the dark matter substructure. This model combines the Eggen, Lynden-Bell and Sandage picture of a monolithic protogalactic collapses with the Searle and Zinn picture of an early merging phase of substructures into one consistent scenario of Galactic formation.

We compute the specific angular momentum distributions of a sample of low mass disk galaxies observed by Swaters. We compare these distributions to those of dark matter haloes obtained by Bullock et al. from high resolution N-body simulations of structure formation in a LCDM Universe. We find that although the disk mass fractions are significantly smaller than the Universal baryon fraction, the total specific angular momenta of the disks are in good agreement with those of dark matter haloes. This suggests that disks form out of only a small fraction of the available baryons, but yet manage to draw most of the available angular momentum. In addition we find that the angular momentum distributions of disks are clearly distinct from those of the dark matter; disks lack predominantly both low and high specific angular momentum. Understanding these findings in terms of a coherent picture for disk formation is challenging. Cooling, feedback and stripping, which are the main mechanisms to explain the small disk mass fractions found, seem unable to simultaneously explain the disk's angular momentum distribution. In fact, it seems that the baryons that make up the disks must have been born out of angular momentum distributions that are clearly distinct from those of LCDM haloes. Therefore, understanding the angular momentum content of disk galaxies remains an important challenge for our picture of galaxy formation.

We investigate detailed kinematical properties of simulated collisionless merger remnants of disk galaxies with mass ratios of 1:1 and 3:1. The simulations were performed by direct summation using the new special hardware device GRAPE-5. In agreement with observations, the shape of the line-of-sight velocity distribution (LOSVD) is Gaussian with small deviations. For most cases we find that the retrograde wings of the LOSVD are steeper than the prograde ones. This is in contradiction with observations which show broad retrograde and steep prograde wings. This serious problem in the collisionless formation scenario of massive elliptical galaxies can be solved if all rotating ellipticals, even boxy ones, contain an additional stellar disk component with $\propto 15%$ of the total stellar mass and a scale length of order the effective radius of the spheroid. We propose that the progenitor galaxies of massive ellipticals must have contained a significant amount of gas that did not condense into stars during the merger process itself but formed an extended gaseous disk before the star formation epoch. The heating source that prevented the gas from forming stars early and the origin of the large specific angular momentum required for the gas component to form an extended disk are still unsolved problems.

Pre-main sequence and main-sequence binary systems are observed to have periods, P, ranging from one day to 10^(10) days and eccentricities, e, ranging from 0 to 1. We pose the problem if stellar-dynamical interactions in very young and compact star clusters may broaden an initially narrow period distribution to the observed width. N-body computations of extremely compact clusters containing 100 and 1000 stars initially in equilibrium and in cold collapse are preformed. In all cases the assumed initial period distribution is uniform in the narrow range 4.5 < log10(P) < 5.5 (P in days) which straddles the maximum in the observed period distribution of late-type Galactic-field dwarf systems. None of the models lead to the necessary broadening of the period distribution, despite our adopted extreme conditions that favour binary--binary interactions. Stellar-dynamical interactions in embedded clusters thus cannot, under any circumstances, widen the period distribution sufficiently. The wide range of orbital periods of very young and old binary systems is therefore a result of cloud fragmentation and immediate subsequent magneto-hydrodynamical processes operating within the multiple proto-stellar system.

We test the formation of low-luminosity elliptical galaxies through collisionless mergers of unequal-mass disk galaxies. The kinematic properties of a small survey of simulated merger end-products with initial disk mass ratios of 3:1 is compared to a sample of seven low-luminosity galaxies observed by Rix et al. that were chosen photometrically to be "ellipticals". In this paper, we go beyond a comparison in terms of global properties (using one number to characterize a model or a galaxy e.g. <a4>, ellipticity at some fixed radius or central velocity dispersion) and examine the detailed kinematics as function of galactocentric distance. The merger remnants are ``observed'' through a slit along the major and minor axis, using a pixel binning and slit width similar to the one used during the spectroscopic observations. We compare the rotational support (v/sigma) of the merger remnants, the normalized rotation on the minor axis and the major axis Gauss-Hermite moments h3, to that of the observed galaxies. The N-body remnants are very flattened when viewed edge-on (<epsilon> = 0.6) and should be inclined before making a fair comparison with the Rix et al. data set (which has <epsilon> = 0.3). When the merger remnants are appropriately inclined, their v/sigma profiles rise slower than the observed ones. Note that even when the remnants are viewed edge-on, the v/sigma profiles do not match the observations. The detailed comparison of the observations with purely collisionless 3:1-merger remnants suggests that this kind of dissipationless merger is not likely to be the dominant formation channel for low-luminosity elliptical galaxies. (abridged)

We performed N-body/SPH simulations of merging gas rich disk-galaxies with a mass ratio of 3:1. A stellar disk and bulge component and a dark halo is realized with collisionless particles, the gas is simulated using SPH particles. In this paper we focus on the gas dynamics in mergers and its influence on the structure of the final merger remnant, neglecting star formation. During the merger around $8\times 10^9 M_{\odot}$ of gas accumulate in the central region of the merger remnant (inside 300 pc). Part of the gas falls to the center, channeled by a tidally induced bar in the more massive galaxy. A peak accretion rate of 150 solar masses/year is reached when the two galaxy centers merge. It is likely that the gas will experience a central starburst and/or fuel a central black hole. The gas in the outer regions accumulates in dense knots within tidal tails which could lead to the formation of open clusters or dwarf satellites. Later on, the gas knots loose angular momentum by dynamical friction and also sink to the center of the remnant, thereby increasing the central gas content. In the end, the gas settles in a central power-law disk with surface density $\Sigma \propto r^{-2}$ surrounded by dense, orbiting gas clumps. The stellar profile is of de Vaucouleurs-type, like in simulations of pure stellar mergers

Collisionless N-body simulations of merging disk-galaxies with mass ratios ($\eta$) of 1:1, 2:1, 3:1, and 4:1 have been performed using direct summation with the special purpose hardware GRAPE. The simulations are used to examine whether the formation of elliptical galaxies can be explained in the context of the merger scenario. The photometric, kinematical and isophotal properties of the merger remnants are investigated and turn out to be in very good agreement with observations of giant elliptical galaxies. We conclude that equal mass mergers lead to slowly rotating, anisotropic remnants having predominantly boxy isophotes. Mergers with a mass ratio of 3:1 and 4:1, on the other hand, are fast isotropic rotators with disky isophotes. Projection effects can explain the observed scatter in the kinematical and isophotal properties of elliptical galaxies.

High-resolution numerical simulations of galaxy mergers are analysed. The global structure and isophotal shapes of the merger remnants are in good agreement with the observations. Whereas equal-mass mergers lead to anisotropic, boxy ellipticals, unequal-mass mergers result in disky and isotropic systems. The line-of-sight velocity distributions show small deviations from a Gaussian distribution. In all cases we find that the retrograde wings are steeper and that the prograde wings are broader than a Gaussian distribution. This is in contradiction with the observations which show broader retrograde and steeper prograde wings in all ellipticals. This fundamental difference between observation and theory can be explained if all ellipticals, even anistropic boxy ones, contain extended stellar disk components with luminosities of order 10 to 20 percent the total luminosity and scale radii of order the effective radii of the spheroids.

We model the star formation history (SFH) and the chemical evolution of the Galactic disk by combining an infall model and a limit-cycle model of the interstellar medium (ISM). Recent observations have shown that the SFH of the Galactic disk violently variates or oscillates. We model the oscillatory SFH based on the limit-cycle behavior of the fractional masses of three components of the ISM. The observed period of the oscillation ($\sim 1$ Gyr) is reproduced within the natural parameter range. This means that we can interpret the oscillatory SFH as the limit-cycle behavior of the ISM. We then test the chemical evolution of stars and gas in the framework of the limit-cycle model, since the oscillatory behavior of the SFH may cause an oscillatory evolution of the metallicity. We find however that the oscillatory behavior of metallicity is not prominent because the metallicity reflects the past integrated SFH. This indicates that the metallicity cannot be used to distinguish an oscillatory SFH from one without oscillations.

The evolution of halos consisting of weakly self-interacting dark matter particles is summarized. The halos initially contain a central density cusp as predicted by cosmological models. Weak self-interaction leads to the formation of an isothermal, low-density core which is in agreement with published data on the rotation curves of dwarf galaxies. However, subsequently, core collapse leads to an even steeper central density cusp. In order to explain the observed dark matter rotation curves by weak interaction, the cross section must be in a narrow range around sigma/mp= 0.1 cm*cm/g.

The observed relation between central black hole mass and spheroid velocity dispersion is interpreted in terms of a self-regulation model that incorporates a viscous Keplerian accretion disk to feed the black hole, embedded in a massive, self-gravitating star forming disk that eventually populates the spheroid. The model leads to a constant ratio between black hole mass and spheroid mass which is equal to the inverse of the critical Reynolds number for the onset of turbulence in the accretion disk surrounding the central black hole. Applying the fundamental plane correlation for spheroids, we find that the black hole mass has a power-law dependence on the spheroid velocity dispersion with a slope in the range of 4-5. We explain the larger scatter in the Magorrian relation with respect to the black hole mass-spheroid velocity dispersion relationship as a result of secular evolution of the spheroid that primarily affects its luminosity and to a much lesser extent its velocity dispersion.

Using hydrodynamic simulations we investigate the time evolution and fragmentation of regions within molecular clouds which have lost their turbulent support leading to gravitational contraction. The initial density distributions are described by random Gaussian fluctuations with varying slopes x of the power spectrum P(k) ~ k^-x, covering the range from flat (x=0) to very steep spectra (x=3). We consider molecular cloud volumes containing different masses relative to the average Jeans mass M_J, from 1 M_J to 222 M_J. This parameter study extends the detailed analysis of systems with initially P(k) ~ k^-2 and mass 222 M_J presented by Klessen & Burkert (2000). The dynamical evolution of the simulated molecular cloud regions is insensitive to the slope of the initial density fluctuation spectrum. The system evolves into a complicated network of intersecting filaments and collapsing clumps leading to the formation of a compact cluster of accreting and interacting embedded protostellar cores. The cluster builds up as bound entity, but dissolves later due to collisional effects. In all simulations, the mass spectrum of collapsed cores is very broad, has approximately log-normal shape and peaks roughly at the average Jeans mass. This supports the hypothesis that the average Jeans mass is the main parameter determining the peak in the stellar spectrum, and suggests that the interplay between self-gravity on the one side and thermal and turbulent pressure on the other side is the dominant process that regulates the formation of stellar clusters.

The Lin and Pringle model (1987) of galactic disk formation postulates that if star formation proceeds on the same timescale as the viscous redistribution of mass and angular momentum in disk galaxies, then the stars attain an exponential density profile. Their claim is that this result holds generally: regardless of the disk galaxy's initial gas and dark matter distribution and independent of the nature of the viscous processes acting in the disk. We present new results from a set of 2D hydro-simulations which investigate their analytic result.

We present a stochastic model based on the Monte Carlo technique to study the inhomogeneous chemical evolution of the Galactic halo. In particular, we consider the local enrichment and subsequent mixing of the interstellar gas resulting from bursts of star formation, following explicitly the fragmentation and coalescence of interstellar gas clouds. The model takes into account the mixing of halo gas through cloud-cloud collisions and the delayed mixing of supernova ejecta into the interstellar medium. The consequences of the infall of halo gas onto the Galactic disk are also discussed. The mass spectrum of clouds, the age-metallicity relation, and the G-dwarf distribution are investigated at different times. We analyze in detail the predictions of our model for the abundance of elements like Eu, Ba, and Sr in the Galactic halo, following their production by $r$-process nucleosynthesis for different assumptions on the supernova mass range. Finally, we compare our results with spectroscopic data for the chemical composition of metal-poor halo stars.

The effect of gas ejection on the structure and binding energy of newly formed stellar clusters is investigated. The star formation efficiency (SFE), necessary for forming a gravitationally bound stellar cluster, is determined. Two sets of numerical N-body simulations are presented: As a first simplified approach we treat the residual gas as an external potential. The gas expulsion is approximated by reducing the gas mass to zero on a given timescale, which is treated as a free parameter. In a second set of simulations we use smoothed particle hydrodynamics (SPH) to follow the dynamics of the outflowing residual gas self-consistently. We investigate cases where gas outflow is induced by an outwards propagating shock front and where the whole gas cloud is heated homogeneously, leading to ejection. If the stars are in virial equilibrium with the gaseous environment initially, bound clusters only form in regions where the local SFE is larger than 50% or where the gas expulsion timescale is long compared to the dynamical timescale. A small initial velocity dispersion of the stars leads to a compaction of the cluster during the expulsion phase and reduces the SFE needed to form bound clusters to less than 10%.

The rotational properties of the visible and dark components of low-mass disk galaxies (vrot<=100 km/s) are investigated using the Swaters sample. The rotational parameter lambda'=lambda_DM*(j_d/m_d) is determined, where lambda_DM is the dark halo spin parameter and j_d/m_d is the ratio between the specific angular momentum of the disk and that of the dark halo. The distribution of lambda' is in excellent agreement with cosmological predictions of hierarchical clustering if j_d/m_d=1, that is if the protogalactic gas did not loose a significant amount of specific angular momentum due to dynamical friction. This result is in disagreement with current cosmological Nbody/SPH simulations where the baryonic component looses 90% of its specific angular momentum while settling into the equatorial plane. The Swaters sample also shows a surprisingly strong correlation between lambda' and the baryonic mass fraction m_d: lambda'=0.4 m_d**0.6. This correlation can be explained if the total amount of gas in protogalaxies is a universal fraction of their dark matter mass, of order 10%, and if the variation in m_d is a result of the fact that only the inner parts of the primordial gas distribution managed to form a visible disk component. In this case the specific angular momentum of the gas out of which the disk formed is a factor of 2.75 larger than that of the dark halo, which would require a yet unknown spin-up process for the visible baryonic component.

The rotational properties of centrally condensed, turbulent molecular cloud cores with velocity fields that are characterized by Gaussian random fields are investigated. It is shown that the observed line width-size relationship can be reproduced if the velocity power spectrum is a power-law with P(k)=k**n and n=-3 to -4. The line-of-sight velocity maps of these cores show velocity gradients that can be interpreted as rotation. For n=-4, both, the deduced values of the angular velocity Omega=1.6 km/s/pc * (R/0.1 pc)**0.5 and the scaling relations between Omega and the core radius R are in very good agreement with the observations. As a result of the dominance of long wavelength modes, the cores also have a net specific angular momentum with an average value of j=7*(10**20)*(R/0.1 pc)**(1.5) cm**2/s with a large spread. Their internal dimensionless rotational parameter is beta=0.03, independent of the scale radius R. In general, the line-of-sight velocity gradient of an individual turbulent core does not provide a good estimate of its internal specific angular momentum. We find however that the distribution of the specific angular momenta of a large sample of cores which are described by the same power spectrum can be determined very accurately from the distribution of their line-of-sight velocity gradients Omega using the simple formula j=p*Omega*R*R where p depends on the density distribution of the core and has to be determined from a Monte-Carlo study. Our results show that for centrally condensed cores the intrinsic angular momentum is overestimated by a factor of 2-3 if p=0.4 is used.

We establish the presence of a dark matter core radius, for the first time in a very large number of spiral galaxies of all luminosities. Contrary to common opinion we find that the sizes of these cores and the " DM core problem" are bigger for more massive spirals. As a result the Burkert profile provides an excellent mass model for dark halos around disk galaxies. Moreover, we find that the spiral dark matter core densities $\rho_{0}$ and core radii $r_{0}$ lie in the same scaling relation $\rho_{0}=4.5\times 10^-2 (r_{0}/kpc)^{-2/3} M_{\odot}pc^{-3}$ of dwarf galaxies with core radii upto ten times more smaller.

The evolution of halos consisting of weakly self-interacting dark matter particles is investigated using a new numerical Monte-Carlo N-body method. The halos initially contain kinematically cold, dense 1/r-power-law cores. For interaction cross sections (sigma/mp) >= 10-100 cm**2/g weak self-interaction leads to the formation of isothermal, constant density cores within a Hubble time as a result of heat transfer into the cold inner regions. This core structure is in good agreement with the observations of dark matter rotation curves in dwarf galaxies. The isothermal core radii and core densities are a function of the halo scale radii and scale masses which depend on the cosmological model. Adopting the currently popular LCDM model, the predicted core radii and core densities are on good agreement with the observations. For large interaction cross sections, massive dark halos with scale radii rs >= 14000 cm**2/g/(sigma/mp) kpc could experience core collapse during their lifetime, leading to cores with singular isothermal density profiles.

The star formation process in molecular clouds usually leads to the formation of multiple stellar systems, mostly binaries. Remaining disks around those stars may be located around individual stars (circumstellar disks) or around the entire binary system (circumbinary disk). We shall briefly review the present observational evidence for both types of disks in binary stars, in particular the properties of circumbinary disks. We then present recent results of the theoretical modeling of the collapse and fragmentation of gravitationally unstable molecular cloud cores and their implications for binary and disk formation, and discuss the dynamical influence of the binary companions on disk truncation and gap formation. The presence of binaries may have profound influence on the process of planet formation as well. We present results on the stability and evolution of orbits of planets in disks around binaries.

The isothermal dynamical evolution of a clumpy and turbulent molecular cloud region and its fragmentation into a protostellar cluster is investigated numerically. The effect of different initial density and velocity distributions, generated from different realizations of a Gaussian random field with power spectrum P(k) ~ k^{-2} is investigated. As a result of the interplay between gas pressure and gravitational forces, a quasi-equilibrium clump mass spectrum emerges with a power-law distribution dN/dM ~ M^{-1.5}, independent of the initial conditions. Being part of a complex network of filaments, individual clumps are elongated, centrally condensed objects with 2:1 to 4:1 axis ratios with outer r^{-2} density distributions. Dense, Jeans-unstable gas clumps collapse and form protostellar cores which evolve through competitive accretion and $N$-body interactions with other cores. In contrast to the clumps, the core mass spectrum is best described by a log-normal distribution which peaks approximately at the average Jeans mass of the system. Scaled to physical conditions typical for star-forming molecular clouds, the mass function is in good agreement with the IMF of multiple stellar systems. The final dynamical state of the newly formed stellar cluster closely resembles observed young stellar clusters. It has a core/halo structure which is typical for collision dominated N-body systems. The 2-point correlation function of the spatial stellar distribution can be described by two power-laws with a break in the slope at the transition point from the binary to the large-scale clustering regime. The protostellar cluster is marginally bound.

A new method for the inclusion of ionizing radiation from uniform radiation fields into 3D Smoothed Particle Hydrodynamics (SPHI) simulations is presented. We calculate the optical depth for the Lyman continuum radiation from the source towards the SPHI particles by ray-tracing integration. The time-dependent ionization rate equation is then solved locally for the particles within the ionizing radiation field. Using test calculations, we explore the numerical behaviour of the code with respect to the implementation of the time-dependent ionization rate equation. We also test the coupling of the heating caused by the ionization to the hydrodynamical part of the SPHI code.

The radiative cooling of optically thin gaseous regions and the formation of a two-phase medium and of cold gas clouds with a clumpy substructure is investigated. In optically thin clouds, the growth rate of small isobaric density perturbations is independent of their length scale. However, the growth of a perturbation is limited by its transition from isobaric to isochoric cooling. The temperature at which this transition occurs decreases with the length scale of the perturbation. Consequently small scale perturbations have the potential to reach higher amplitudes than large scale perturbations. When the amplitude becomes nonlinear, advection overtakes the pressure gradient in promoting the compression resulting in an accelerated growth of the disturbance. The critical temperature for transition depends on the initial amplitude. The fluctuations which can first reach nonlinearity before their isobaric to isochoric transition will determine the characteristic size and mass of the cold dense clumps which would emerge from the cooling of an initially nearly homogeneous region of gas. Thermal conduction is in general very efficient in erasing isobaric, small-scale fluctuations, suppressing a cooling instability. A weak, tangled magnetic field can however reduce the conductive heat flux enough for low-amplitude fluctuations to grow isobarically and become non-linear if their length scales are of order 0.01 pc. Finally, we demonstrate how a 2-phase medium, with cold clumps being pressure confined in a diffuse hot residual background component, would be sustained if there is adequate heating to compensate the energy loss.

Many clusters of galaxies contain an appreciable amount of hot gas, the intracluster medium. As a consequence, gas will be stripped from galaxies that move through the inter cluster medium, if the ram pressure exceeds the internal gravitational force. We study the interaction between the intracluster medium and an extended gas component of dwarf galaxies confined by a surrounding cold dark matter halo analytically and numerically, using axisymmetric two-dimensional hydrodynamical simulations at high resolution. The results show that the gas within the dark matter halo is totally stripped in a typical galactic cluster. The process of ram pressure stripping therefore must have played an important role during the chemo-dynamical evolution of dwarf galaxies in galactic clusters. Our results predict that most of the chemical evolution of dwarf galaxies must have occurred at high redshift, before the virialized cluster had formed.

N-body simulations of major spiral-spiral mergers with different mass ratios q <= 4:1 are presented. All mergers lead to spheroidal stellar systems with de Vaucouleurs-like surface density profiles, resembling elliptical galaxies. The internal dynamical structure of the remnants does however depend critically on the adopted mass ratio. Equal-mass mergers with q <= 2:1 form radially anisotropic, slowly rotating systems with preferentially boxy isophotes. Unequal-mass mergers with q > 2:1 result in the formation of rotationally supported ellipticals with disky isophotes and small minor axis rotation. Projection effects lead in general to a large scatter in the kinematical and isophotal properties. If the observed dichotomy between boxy and disky ellipticals is a result of different mass ratios of the merger components massive ellipticals must have formed preferentially from equal-mass mergers in contrast to low-mass ellipticals which had to form from unequal-mass mergers.

The evolution of low mass (M < 2.5 Msun) binaries through the common envelope phase has been studied for systems in which one member is on its first ascent of the red giant branch. Three-dimensional hydrodynamical simulations have been carried out for a range of red giant masses (1 - 2 Msun) with degenerate helium cores (0.28 - 0.45 Msun) and companions (0.1 - 0.45 Msun) for initial orbital periods ranging from about 15 - 1000 days. The results suggest that these low mass binary systems can survive the common envelope phase provided that the helium degenerate core is more massive than about 0.2 - 0.25 Msun and that the mass of the red giant progenitor is less than about 2 Msun. For the observed short period double degenerate systems, it is found that evolutionary scenarios involving two phases of common envelope evolution are not likely and that a scenario involving an Algol-like phase of mass transfer followed by a common envelope phase is viable, suggesting that the first-formed white dwarf is often reheated by nuclear burning on its surface. A formation mechanism for two subdwarf B stars observed in eclipsing short period binaries with low mass main sequence stars is also described.

We present optical and near-infrared colour maps of the central regions of bulges of S0 and spiral galaxies obtained with WFPC2 and NICMOS on the HST. By combined use of HST and ground-based data the colour information spans from a few tens of pc to a few kpc. In almost all galaxies the colour profiles in the central 100-200 pc become more rapidly redder. We attribute the high central colour indices to a central concentration of dust. We infer an average extinction at the center of A_V = 0.6 - 1.0 mag. Several objects show central dust rings or disks at sub-kpc scales similar to those found by others in giant ellipticals. For galactic bulges of types S0 to Sb, the tightness of the B-I vs I-H relation suggests that the age spread among bulges of early type galaxies is small, of at most 2 Gyr. Colours at 1 R_eff, where we expect extinction to be negligible are similar to than those of elliptical galaxies in the Coma cluster, suggesting that these bulges formed at the same time as the bright galaxies in Coma. Furthermore the galaxy ages are found to be independent of their environment. Since it is likely that Coma was formed at redshift z > 3, our bulges, which are in groups and in the field, also must have been formed at this epoch. Bulges of early-type spirals cannot be formed by secular evolution of bars at recent epochs, since such bulges would be much younger. There are three galaxies of type Sbc and later, their bulges are younger and could perhaps arise from secular evolution of transient bars. Our results are in good agreement with semi-analytic predictions (Baugh et al. 1996, 1998, Kauffmann 1996), who also predict that bulges, in clusters and in the field, are as old as giant ellipticals in clusters.

We describe the basic features of a Monte Carlo model specifically designed to follow the inhomogenous chemical evolution of the Galactic halo, taking into account the effects of local enrichment and mixing of the halo gas, and with particular emphasis on elements like Eu produced by r-process nucleosynthesis. We compare our results with spectroscopic data for the chemical composition of metal-poor halo stars and globular clusters like M13, M5, M92, M4, and we infer some constraints on the star formation history of the halo and the rate and mass spectrum of supernovae during the first epoch of Galaxy evolution.

The origin of boxy and disky elliptical galaxies is investigated. The results of two collisionless N-body simulations of spiral-spiral mergers with mass ratios of 1:1 and 3:1 are discussed and the projected properties of the merger remnants are investigated. It is shown that the equal-mass merger leads to an anisotropic, slowly rotating system with preferentially boxy isophotes and significant minor axis rotation. The unequal-mass merger results in the formation of a rotationally supported elliptical with disky isophotes and small minor axis rotation. The observed scatter in the kinematical and isophotal properties of both classes of elliptical galaxies can be explained by projection effects.

We study a possible formation mechanism for elliptical/S0 galaxies using N-body simulations with GRAPE. A close galaxy encounter which does not lead to a merger can induce a strong bar in an initially axisymmetric disk. This bar is unstable against bending oscillations and loses angular momentum to the dark halo. After 6 Gyrs the remnant is spheroidal with a de Vaucouleurs-type surface profile and has ellipticities and rotation properties comparable to observed elliptical/S0 galaxies: the system rotates fast (v/sigma \approx 1) and has disky isophotes.

We have investigated effects of early stellar encounters on a protoplanetary disk (planetesimal disk) and found that they can explain the high eccentricities and inclinations observed in the outer part ($>42$AU) of the Edgeworth-Kuiper Belt (EKB). The proto-sun is considered as a member of a stellar aggregation that undergoes dissolution on a timescale $\sim 10^8$yrs, such that a planetesimal disk experiences a flyby encounter at pericenter distance ($q$) on the order of 100AU. We simulated dynamical evolution of a planetesimal (test particle) disk perturbed by a passing star. We show that the stellar encounter pumps the velocity dispersion in the disk in the outer parts ($> 0.25q$). Planet formation is forestalled in that region. We also find that a stellar encounter with pericenter distance $q \sim 100-200$AU could have pumped up the velocity dispersion of EKB objects outside 42AU to the observed magnitude while preserving that inside Neptune's 3:2 mean-motion resonance (located at 39.5AU), which allows for the efficient capture of objects by the sweeping of the 3:2 resonance during orbital migration by proto-Neptune.

It is shown that the cuspy density distributions observed in the cores of elliptical galaxies can be realized by dissipationless gravitational collapse. The initial models consist of power-law density spheres such as $\rho\propto r^{-1}$ with anisotropic velocity dispersions. Collapse simulations are carried out by integrating the collisionless Boltzmann equation directly, on the assumption of spherical symmetry. From the results obtained, the extent of constant density cores, formed through violent relaxation, decreases as the velocity anisotropy increases radially, and practically disappears for extremely radially anisotropic models. As a result, the relaxed density distributions become more cuspy with increasing radial velocity anisotropy. It is thus concluded that the velocity anisotropy could be a key ingredient for the formation of density cusps in a dissipationless collapse picture. The velocity dispersions increase with radius in the cores according to the nearly power-law density distributions. The power-law index, $n$, of the density profiles, defined as $\rho\propto r^{-n}$, changes from $n\approx 2.1$ at intermediate radii, to a shallower power than $n\approx 2.1$ toward the centre. This density bend can be explained from our postulated local phase-space constraint that the phase-space density accessible to the relaxed state is determined at each radius by the maximum phase-space density of the initial state.

The structure of dark matter halos as predicted from cosmological models is discussed and compared with observed rotation curves of dark matter-dominated dwarf galaxies. The theoretical models predict that dark matter halos represent a one-parameter family with a universal density profile. Observations of dark matter-dominated rotation curves indeed confirm the universal structure of dark halos. They are even in excellent agreement with the expected mass-radius scaling relations for the currently favoured cosmological model (standard cold dark matter with Omega_0=0.25 and Omega_{Lambda}=0.75). The rotation profiles however disagree with the predicted dark matter density distributions. Secular processes which might affect the inner halo structure do not seem to provide a good solution to this problem. We discuss, as an alternative, the possibility that dark halos consist of two separate components, a dark baryonic and a dark non-baryonic component.

In order to understand the dynamical and chemical evolution of our Galaxy it is of fundamental importance to study the local neighborhood. White dwarf stars are ideal candidates to probe the history of the solar neighborhood, since these ``fossil'' stars have very long evolutionary time-scales and, at the same time, their evolution is relatively well understood. In fact, the white dwarf luminosity function has been used for this purpose by several authors. However, a long standing problem arises from the relatively poor statistics of the samples, especially at low luminosities. In this paper we assess the statistical reliability of the white dwarf luminosity function by using a Monte Carlo approach.

The evolution of a giant planet within the stellar envelope of a main-sequence star is investigated as a possible mechanism for enhancing the stellar metallicities of the parent stars of extrasolar planetary systems. Three-dimensional hydrodynamical simulations of a planet subject to impacting stellar matter indicate that the envelope of a Jupiter-like giant planet can be completely stripped in the outer stellar convection zone of a solar-mass star. In contrast, Jupiter-like and less massive Saturn-like giant planets are able to survive through the base of the convection zone of a 1.22 solar-mass star. Although strongly dependent on details of planetary interior models, partial or total dissolution of giant planets can result in significant enhancements in the metallicity of host stars with masses between about 1.0 and 1.3 solar masses. The implications of these results with regard to planetary orbital migration are briefly discussed.

The isothermal gravitational collapse and fragmentation of a region within a molecular cloud and the subsequent formation of a protostellar cluster is investigated numerically. The clump mass spectrum which forms during the fragmentation phase can be well approximated by a power law distribution dN/dM ~ M^(-1.5). In contrast, the mass spectrum of protostellar cores that form in the centers of Jeans-unstable clumps and evolve through accretion and N-body interactions is described by a log-normal distribution with a width that is in excellent agreement with observations of multiple stellar systems.

The isothermal gravitational collapse and fragmentation of a molecular cloud region and the subsequent formation of a protostellar cluster is investigated numerically. The clump mass spectrum which forms during the fragmentation phase can be well approximated by a power law distribution dN/dM ~ M^(-1.5). In contrast, the mass spectrum of protostellar cores that form in the centers of Jeans unstable clumps and evolve through accretion and N-body interaction is best described by a log-normal distribution. Assuming a star formation efficiency of ~10%, it is in excellent agreement with the IMF of multiple stellar systems.

The evolution of binary systems consisting of an asymptotic giant branch star of mass equal to 3 M_sun or 5 M_sun, and a main sequence star of mass equal to 0.4 M_sun or 0.6 M_sun with orbital periods > 200 days has been followed from the onset through the late stages of the common envelope phase. Using a nested grid technique, the three-dimensional hydrodynamical simulations of an asymptotic giant branch star with radii approximately 1 A.U. indicate that a significant fraction of the envelope gas is unbound (about 31% and 23% for binaries of 3 M_sun and 0.4 M_sun, and 5 M_sun and 0.6 M_sun respectively) by the ends of the simulations, and that the efficiency of the mass ejection process is about 40%. While the original volume of the giant is virtually evacuated in the late stages, most of the envelope gas remains marginally bound on the grid. At the ends of our simulations, when the orbital decay timescale exceeds about 5 years, the giant core and companion orbit one another with a period of about 1 day (2.4 days for a binary involving a more evolved giant), although this is an upper limit to the final orbital period. For a binary of 5 M_sun and 0.4 M_sun, the common envelope may not be completely ejected.

We present numerical studies of compressible, decaying turbulence, with and without magnetic fields, with initial rms Alfven and Mach numbers ranging up to five, and apply the results to the question of the support of star-forming interstellar clouds of molecular gas. We find that, in 1D, magnetized turbulence actually decays faster than unmagnetized turbulence. In all the regimes that we have studied 3D turbulence-super-Alfvenic, supersonic, sub-Alfvenic, and subsonic-the kinetic energy decays as (t-t0)^(-x), with 0.85 < x < 1.2. We compared results from two entirely different algorithms in the unmagnetized case, and have performed extensive resolution studies in all cases, reaching resolutions of 256^3 zones or 350,000 particles. We conclude that the observed long lifetimes and supersonic motions in molecular clouds must be due to external driving, as undriven turbulence decays far too fast to explain the observations.

Using smoothed particle hydrodynamics in combination with the special-purpose hardware device GRAPE, we follow the fragmentation process and dynamical evolution in the interior of molecular clouds until most of the gas is converted to form a cluster of protostellar cores. The properties of this protostellar cluster are compared to observational data from star forming regions.

The best measured rotation curve for any galaxy is that of the dwarf spiralXXXX DDO 154, which extends out to about 20 disk scale lengths. It provides an ideal laboratory for testing the universal density profile prediction from high resolution numerical simulations of hierarchical clustering in cold dark matter dominated cosmological models. We find that the observed rotation curve cannot be fit either at small radii, as previously noted, or at large radii. We advocate a resolution of this dilemma by postulating the existence of a dark spheroid of baryons amounting to several times the mass of the observed disk component and comparable to that of the cold dark matter halo component. Such an additional mass component provides an excellent fit to the rotation curve provided that the outer halo is still cold dark matter-dominated with a density profile and mass-radius scaling relation as predicted by standard CDM-dominated models. The universal existence of such dark baryonic spheroidal components provides a natural explanation of the universal rotation curves observed in spiral galaxies, may have a similar origin and composition to the local counterpart that has been detected as MACHOs in our own galactic halo via gravitational microlensing, and is consistent with, and even motivated by, primordial nucleosynthesis estimates of the baryon fraction.

Hydrodynamical calculations in three space dimensions of the collapse of an isothermal, rotating 1 M\sol protostellar cloud are presented. The initial density stratification is a power law with density $\rho \propto r^{-p}$, with $p=1$. The case of the singular isothermal sphere ($p=2$) is not considered; however $p=1$ has been shown observationally to be a good representation of the density distribution in molecular cloud cores just before the beginning of collapse. The collapse is studied with two independent numerical methods, an SPH code with 200,000 particles, and a finite-difference code with nested grids which give high spatial resolution in the inner regions. Although previous numerical studies have indicated that such a power-law distribution would not result in fragmentation into a binary system, both codes show, in contrast, that multiple fragmentation does occur in the central regions of the protostar. Thus the process of binary formation by fragmentation is shown to be consistent with the fact that a large fraction of young stars are observed to be in binary or multiple systems.

The rotation curves of dark matter dominated dwarf galaxies are analysed. The observations show that dark matter halos represent a one-parameter family with self similar density profiles. The global halo parameters, like total mass and scale length are coupled by simple scaling relations. The inner halo regions resemble non-singular, isothermal spheres with constant density cores. The observations are compared with dark matter halos, resulting from cosmological cold dark matter simulations. The theoretical models predict that dark matter halos represent a one-parameter family in agreement with the observations. However, in contradiction to the observations, the calculations lead to dark matter halos with $r^{-1}$ density cusps in the center and non-isothermal velocity dispersion profiles. Processes which might affect the inner halo structure, resulting in isothermal, constant density cores are discussed.

The amount of dark matter in the four galactic dwarf spheroidals with large mass-to-light ratios is investigated. Sextans has a cut-off radius which is equal to the expected tidal radius, assuming a high mass-to-light ratio. This satellite very likely is dark matter dominated. Carina, Ursa Minor and Draco, on the other hand, cannot contain a dominating dark matter component if the observed 'extra-tidal' stars are located exterior to the tidal radii of these systems. The evidence for tidal stripping in the absence of dark matter is also supported by the fact that the observed cut-off radii of all three satellites are equal to their tidal radii, assuming a low, globular cluster like mass-to-light ratio. The large velocity dispersions of these galaxies, on the other hand, seem to provide strong evidence for a massive dark matter component. In this case, the 'extra-tidal' stars lie deeply embedded in the dark matter potential wells of the satellites. These stars then would represent a gravitationally bound, extended stellar component with unknown origin.

The kinematical and spatial properties of the metal-rich globular clusters in the Galaxy ([Fe/H] > -0.8) indicates that these objects do not comprise a homogeneous population. Three subsystems are identified. The highest-mass clusters exhibit a very slow net rotation with a speed of v_rot = 24 km/s and v_rot/sigma = 0.3, indicative of a centrally condensed, spheroidal subsystem. Roughly half of the lower-mass clusters appear to be located in an elongated bar-like structure which passes through the Galactic Center, and has similar properties to the central stellar bar of the Milky Way. The remaining lower-mass clusters exhibit very rapid net rotation, with a rotation speed of v_rot = 164 km/s and v_rot/sigma = 6. These clusters are located in the Galactic plane, within a ring of 4 to 6 kpc radial distance from the Galactic Center. The highest-mass clusters may have formed during relatively advanced stages of the dissipative evolution of the inner Galactic halo. The lower-mass bar clusters may be associated with the formation of the Galactic stellar bar or bulge. The lower-mass ring clusters appear to be real disk objects. They may represent a stage in cluster formation that was intermediate between that of the halo globular clusters and the oldest extant open clusters.

A review of type Ia supernovae (SNe Ia) statistics as expected from the different scenarios and explosion models so far suggested is presented.

A new solution is presented for the puzzling, observed universality of the exponential luminosity profiles, perpendicular to the disk plane of spiral and lenticular galaxies. It is shown that such exponential $z$-profiles result naturally from gaseous protodisks which settle into isothermal equilibrium prior to star formation. Subsequent cooling leads to a gravitational contraction of the gas towards the equatorial plane and to a stellar exponential $z$-profile if the star formation rate is assumed to be comparable to the cooling rate of the gas. The final stellar scale height depends only on the initial gas temperature and local surface density. This model therefore provides a new method to investigate the early energetic state of galactic protodisks with measured scale heights and surface densities along the disk plane.

Some dwarf spheroidal galaxies have experienced many Gyrs long periods without star formation. We show that Type Ia supernovae, formed from a first generation of stars, can delay a second epoch of star formation by several Gyrs, if the total gas mass in the dwarf spheroidals is smaller than $10^8 M_{\odot}$. The details depend on the assumed Type Ia supernova model and the mass fraction of stars formed in the first population. This scenario requires that dwarf spheroidals be confined by dark matter halos, as otherwise the heated gas would have been lost due to tidal effects and galactic winds. If Type Ia supernovae played an important role, heating the gas in dwarf spheroidals, they should also have increased their iron abundance during the epoch without star formation. The subsequent stellar generation should then start with a higher [Fe/H] than that of the most iron-rich stars found in the previous generation.

The occurrence and properties of Type Ia supernovae (SNe Ia) in single-degenerate binary systems (white dwarf [WD] + nondegenerate companion) is examined for galaxies of different types, and as a function of redshift. The rates and characteristics (peak luminosities, expansion velocities of the ejecta) expected from the explosion of mass-accreting WDs in symbiotic systems and "helium star cataclysmics" are found to be different from those arising in another class of candidate systems: cataclysmic-like (contact) systems (CLS), where a CO WD accretes hydrogen on a thermal time scale from a Roche-lobe filling main-sequence or subgiant companion. We derive the evolution of the SNe Ia rate and properties resulting from the thermonuclear explosion of sub-Chandrasekhar mass WDs in such systems when they detonate a helium layer accumulated from steady burning of hydrogen at the surface. A fraction of CLS are believed to form a subset of the observed luminous supersoft X-ray sources (SSS). Sub-Chandrasekhar explosions from CLS are disfavored in all types of galaxies at redshifts z > 1. On the other hand, CLS where the WD succeeds to grow to the Chandrasekhar mass are more likely found in spiral galaxies, and absent from early-type galaxies. SNe Ia statistics could (if the uncertainties still involved are reduced) help to discriminate among proposed SNe Ia scenarios. The range of variation of the characteristics of SNe Ia in the CLS scenario should be narrower than in symbiotics. The predicted correlation between peak luminosity and expansion velociy of the ejecta in SNe Ia coming from these systems is weak. For CLS, the distinction between the characteristics of SNe Ia respectively arising from sub-Chandrasekhar and from Chandrasekhar-mass explosions should be sharp, since all sub-Chandrasekhar

This work considers the idea of massive black holes being the constituents of the Galactic dark matter halo. It constrains the maximum black hole mass to $\mbh \sil 5 \times 10^4 \msun$ by examining their influence on the population of globular clusters in our Milky Way. In the adopted halo model, globular clusters are exposed to constant bombardment of halo objects on their orbits through the Galaxy and thus will steadily gain internal energy. Depending on the mass of these halo objects and the structural parameters of the globular clusters, they can be disrupted on time scales of a few billion years and below. These disruption time scales are calculated using a modification of the well known (classical) impulsive approximation and compared with direct N-body simulations of such encounter events to ensure the method works correctly. For a set of ten prototypical globular cluster models and black hole masses ranging from $10^3$ to $10^7\:\msun$, Monte-Carlo-simulations of $10\,000$ encounter histories over the period of 10 billion years were calculated each, at three different galactocentric distances \mbox{$R = 5$}, $10$ and $15\, \kpc$. These data were compared with the real globular cluster population in our Galaxy and used to obtain the above constraint of $\mbh \sil 5 \times 10^4 \msun$.

Hydrodynamical calculations in three space dimensions of the collapse of an isothermal, centrally condensed, rotating 1 M\sol protostellar cloud are presented. A numerical algorithm involving nested subgrids is used to resolve the region where fragmentation occurs in the central part of the protostar. A previous calculation by Boss, which produced a hierarchical multiple system, is evolved further, at comparable numerical resolution, and the end result is a binary, with more than half of the mass of the original cloud, whose orbital separation increases with time as a result of accretion of high-angular momentum material and as a result of merging with fragments that have formed farther out. Repeating the calculation with significantly higher resolution, we find that a sequence of binaries can be induced by fragmentation of circumbinary disks. The stability of the resulting multiple system is investigated Using $n$-body calculations, which indicate that such a system would transform on a short time scale into a more stable hierarchical structure. The outermost and most massive binary which forms in the high-resolution run has properties similar to that of the binary found in the low-resolution calculation. Thus the basic outcome is shown to be independent of the numerical spatial resolution. The high-resolution run, in addition, leads to the formation of a system of smaller fragments, which might be important for the understanding of the origin of close binaries with low-mass components and of low mass single stars.

The dependence of the Type Ia supernova (SN Ia) rate on galaxy type is examined for three currently proposed scenarios: merging of a Chandrasekhar--mass CO white dwarf (WD) with a CO WD companion, explosion of a sub--Chandrasekhar mass CO WD induced by accretion of material from a He star companion, and explosion of a sub--Chandrasekhar CO WD in a symbiotic system. The variation of the SNe Ia rate and explosion characteristics with time is derived, and its correlation with parent population age and galaxy redshift is discussed. Among current scenarios, CO + He star systems should be absent from E galaxies. Explosion of CO WDs in symbiotic systems could account for the SNe Ia rate in these galaxies. The same might be true for the CO + CO WD scenario, depending on the value of the common envelope parameter. A testable prediction of the sub--Chandrasekhar WD model is that the average brightness and kinetic energy of the SN Ia events should increase with redshift for a given Hubble type. Also for this scenario, going along the Hubble sequence from E to Sc galaxies SNe Ia events should be brighter on average and should show larger mean velocities of the ejecta. The observational correlations strongly suggest that the characteristics of the SNe Ia explosion are linked to parent population age. The scenario in which WDs with masses below the Chandrasekhar mass explode appears the most promising one to explain the observed variation of the SN Ia rate with galaxy type together with the luminosity--expansion velocity trend.

From spectrophotometric and chemical evolutionary modelling of interacting and merging galaxies we predict abundances of burst stars and of globular clusters (GCs) that may be formed in strong bursts. Observations of young GCs in merger remnants confirm the predicted metallicities. We investigate their formation as a result of a strong star burst with a high star formation efficiency over a large volume of the merger remnant. The evolution of GCs of various metallicities in terms of broad band colours and stellar metal indices are modeled. Comparing our results to the V-I, Mgb and Fe5270 observations of GCs in NGC 7252 largely removes the age metallicity degeneracy and allows for accurate age dating as well as for estimating the effects of self enrichment.

Recent observations indicate that dark matter haloes have flat central density profiles. Cosmological simulations with non-baryonic dark matter predict however self similar haloes with central density cusps. This contradiction has lead to the conclusion that dark matter must be baryonic. Here it is shown that the dark matter haloes of dwarf spiral galaxies represent a one parameter family with self similar density profiles. The observed global halo parameters are coupled with each other through simple scaling relations which can be explained by the standard cold dark matter model if one assumes that all the haloes formed from density fluctuations with the same primordial amplitude. We find that the finite central halo densities correlate with the other global parameters. This result rules out scenarios where the flat halo cores formed subsequently through violent dynamical processes in the baryonic component. These cores instead provide important information on the origin and nature of dark matter in dwarf galaxies.

It is shown that the violent relaxation of dissipationless stellar systems leads to universal de Vaucouleurs profiles only outside 1.5 effective radii $R_e$. Inside $1.5 R_e$ the surface density profiles depend strongly on the initial conditions and are in general not in agreement with the de Vaucouleurs law. This result is in contradiction to the observations which show that all elliptical galaxies have de Vaucouleurs profiles inside $1.5 R_e$ and show strong deviations outside $1.5 R_e$. The $r^{1/4}$-profiles of E-galaxies are therefore not a direct result of collisionless violent relaxation and angular momentum redistribution. It is rather a universal initial density profile which is required in order to explain the universal structure of E-galaxies. A new model for the formation of E-galaxies is presented which assumes that ellipticals and bulges formed from isothermal, gaseous spheres which accumulated in the core regions of dark matter halos and condensed into stars through a cooling instability. This scenario predicts that early type galaxies trace directly the core radii and core masses of their dark matter halos, providing important information on the nature of the dark matter component in galaxies.