Source author record

Andreas Burkert

Andreas Burkert appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph astro-ph.SR astro-ph.HE astro-ph.IM physics.flu-dyn physics.plasm-ph

Catalog footprint

What is connected

96works

8topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2025arXiv

NOEMA$^{\rm 3D}$: A first kpc resolution study of a $z\sim1.5$ main sequence barred galaxy channeling gas into a growing bulge

We present a very deep CO(3-2) observation of a massive, gas-rich, main sequence, barred spiral galaxy at $z\approx1.52$. Our data were taken with the IRAM-NOEMA interferometer for a 12-antenna equivalent on-source integration time of $\sim$ 50 hours. We fit the major axis kinematics using forward modelling of a rotating disk, and then subtract the two-dimensional beam convolved best-fit model revealing signatures of planar non-circular motions in the residuals. The inferred in-plane radial velocities are remarkably large, of the order of $\approx60$ km/s. Direct comparisons with a high-resolution, simulated, gas-rich, barred galaxy, obtained with the moving mesh code AREPO and the TNG sub-grid model, show that the observed non-circular gas flows can be explained as radial flows driven by the central bar, with an inferred net inflow rate of the order of the SFR. Given the recent evidence for a higher-than-expected fraction of barred disk galaxies at cosmic noon, our results suggest that rapid gas inflows due to bars could be important evolutionary drivers for the dominant population of star-forming galaxies at the peak epoch of star and galaxy formation.

preprint2023arXiv

Taking off the edge -- simultaneous filament and end core formation

Simulations of idealised star-forming filaments of finite length typically show core growth which is dominated by two cores forming at its respective end. The end cores form due to a strong increasing acceleration at the filament ends which leads to a sweep-up of material during the filament collapse along its axis. As this growth mode is typically faster than any other core formation mode in a filament, the end cores usually dominate in mass and density compared to other cores forming inside a filament. However, observations of star-forming filaments do not show this prevalence of cores at the filament ends. We explore a possible mechanism to slow the growth of the end cores using numerical simulations of simultaneous filament and embedded core formation, in our case a radially accreting filament forming in a finite converging flow. While such a setup still leads to end cores, they soon begin to move inwards and a density gradient is formed outside of the cores by the continued accumulation of material. As a result, the outermost cores are not longer located at the exact ends of the filament and the density gradient softens the inward gravitational acceleration of the cores. Therefore, the two end cores do not grow as fast as expected and thus do not dominate over other core formation modes in the filament.

preprint2022arXiv

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Spectral lines of ammonia, NH$_3$, are useful probes of the physical conditions in dense molecular cloud cores. In addition to advantages in spectroscopy, ammonia has also been suggested to be resistant to freezing onto grain surfaces, which should make it a superior tool for studying the interior parts of cold, dense cores. Here we present high-resolution NH$_3$ observations with the Very Large Array (VLA) and Green Bank Telescope (GBT) towards a prestellar core. These observations show an outer region with a fractional NH$_3$ abundance of X(NH$_3$) = (1.975$\pm$0.005)$\times 10^{-8}$ ($\pm 10\%$ systematic), but it also reveals that after all, the X(NH$_3$) starts to decrease above a H$_2$ column density of $\approx 2.6 \times 10^{22}$ cm$^{-2}$. We derive a density model for the core and find that the break-point in the fractional abundance occurs at the density n(H$_2$) $\sim 2\times10^5$ cm$^{-3}$, and beyond this point the fractional abundance decreases with increasing density, following the power law $n^{-1.1}$. This power-law behavior is well reproduced by chemical models where adsorption onto grains dominates the removal of ammonia and related species from the gas at high densities. We suggest that the break-point density changes from core to core depending on the temperature and the grain properties, but that the depletion power law is anyway likely to be close to $n^{-1}$ owing to the dominance of accretion in the central parts of starless cores.

preprint2022arXiv

Characterizing the 3D Kinematics of Young Stars in the Radcliffe Wave

We present an analysis of the kinematics of the Radcliffe Wave, a 2.7-kpc-long sinusoidal band of molecular clouds in the solar neighborhood recently detected via 3D dust mapping. With Gaia DR2 astrometry and spectroscopy, we analyze the 3D space velocities of $\sim 1500$ young stars along the Radcliffe Wave in action-angle space, using the motion of the wave's newly born stars as a proxy for its gas motion. We find that the vertical angle of young stars -- corresponding to their orbital phase perpendicular to the Galactic plane -- varies significantly as a function of position along the structure, in a pattern potentially consistent with a wave-like oscillation. This kind of oscillation is not seen in a control sample of older stars from Gaia occupying the same volume, disfavouring formation channels caused by long-lived physical processes. We use a ``wavy midplane'' model to try to account for the trend in vertical angles seen in young stars, and find that while the best-fit parameters for the wave's spatial period and amplitude are qualitatively consistent with the existing morphology defined by 3D dust, there is no evidence for additional velocity structure. These results support more recent and/or transitory processes in the formation of the Radcliffe Wave, which would primarily affect the motion of the wave's gaseous material. Comparisons of our results with new and upcoming simulations, in conjunction with new stellar radial velocity measurements in Gaia DR3, should allow us to further discriminate between various competing hypotheses.

preprint2022arXiv

Star formation near the Sun is driven by expansion of the Local Bubble

For decades we have known that the Sun lies within the Local Bubble, a cavity of low-density, high-temperature plasma surrounded by a shell of cold, neutral gas and dust. However, the precise shape and extent of this shell, the impetus and timescale for its formation, and its relationship to nearby star formation have remained uncertain, largely due to low-resolution models of the local interstellar medium. Leveraging new spatial and dynamical constraints from the Gaia space mission, here we report an analysis of the 3D positions, shapes, and motions of dense gas and young stars within 200 pc of the Sun. We find that nearly all the star-forming complexes in the solar vicinity lie on the surface of the Local Bubble and that their young stars show outward expansion mainly perpendicular to the bubble's surface. Tracebacks of these young stars' motions support a scenario where the origin of the Local Bubble was a burst of stellar birth and then death (supernovae) taking place near the bubble's center beginning 14 Myr ago. The expansion of the Local Bubble created by the supernovae swept up the ambient interstellar medium into an extended shell that has now fragmented and collapsed into the most prominent nearby molecular clouds, in turn providing robust observational support for the theory of supernova-driven star formation.

preprint2021arXiv

Hierarchical fragmentation in high redshift galaxies revealed by hydrodynamical simulations

High-redshift star-forming galaxies have very different morphologies compared to nearby ones. Indeed, they are often dominated by bright star-forming structures of masses up to $10^{8-9}$ $\mathrm{M}_\odot$ dubbed «giant clumps». However, recent observations questioned this result by showing only low-mass structures or no structure at all. We use Adaptative Mesh Refinement hydrodynamical simulations of galaxies with parsec-scale resolution to study the formation of structures inside clumpy high-redshift galaxies. We show that in very gas-rich galaxies star formation occurs in small gas clusters with masses below $10^{7-8}$ $\mathrm{M}_\odot$ that are themselves located inside giant complexes with masses up to $10^8$ and sometimes $10^9$ $\mathrm{M}_\odot$. Those massive structures are similar in mass and size to the giant clumps observed in imaging surveys, in particular with the Hubble Space Telescope. Using mock observations of simulated galaxies, we show that at very high resolution with instruments like the Atacama Large Millimeter Array or through gravitational lensing, only low-mass structures are likely to be detected, and their gathering into giant complexes might be missed. This leads to the non-detection of the giant clumps and therefore introduces a bias in the detection of these structures. We show that the simulated giant clumps can be gravitationally bound even when undetected in mocks representative for ALMA observations and HST observations of lensed galaxies. We then compare the top-down fragmentation of an initially warm disc and the bottom-up fragmentation of an initially cold disc to show that the process of formation of the clumps does not impact their physical properties.

preprint2020arXiv

3D simulations of clump formation in stellar wind collisions

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.

preprint2020arXiv

Accretion driven turbulence in filaments II: Effects of self-gravity

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.

preprint2020arXiv

Dwarfs in the Milky Way halo outer rim: first in-fall or backsplash satellites?

Leo T is a gas-rich dwarf located at 414kpc $(1.4R_{\rm vir})$ distance from the Milky Way (MW) and it is currently assumed to be on its first approach. Here, we present an analysis of orbits calculated backward in time for the dwarf with our new code {\sc delorean}, exploring a range of systematic uncertainties, e.g. MW virial mass and accretion, M31 potential, and cosmic expansion. We discover that orbits with tangential velocities in the Galactic Standard-of-Rest frame lower than $|\vec{u}_{\rm t}^{\rm GSR}|\!\leq\! 63^{+47}_{-39}{\rm\, km\, s^{-1}}$ result in backsplash solutions, i.e. orbits that entered and left the MW dark matter halo in the past, and that velocities above $|\vec{u}_{\rm t}^{\rm GSR}|\!\geq\!21^{+33}_{-21}{\rm\, km\, s^{-1}}$ result in wide orbit backsplash solutions with a minimum pericenter range of $D_{\rm min}\!\geq\!38^{+26}_{-16}{\rm \,kpc}$, which would allow this satellite to survive gas stripping and tidal disruption. Moreover, new proper motion estimates match with our region of backsplash solutions. We applied our method to other distant MW satellites, finding a range of gas stripped backsplash solutions for the gas-less Cetus and Eridanus II, providing a possible explanation for their lack of cold gas, while only first in-fall solutions are found for the HI rich Phoenix I. We also find that the cosmic expansion can delay their first pericenter passage when compared to the non-expanding scenario. This study explores the provenance of these distant dwarfs and provides constraints on the environmental and internal processes that shaped their evolution and current properties.

preprint2020arXiv

Efficient methanol production on the dark side of a prestellar core

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.

preprint2020arXiv

Globular Cluster Formation from Colliding Substructure

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.

preprint2020arXiv

Kinematics of Simulated Galaxies II: Probing the Stellar Kinematics of Galaxies out to Large Radii

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 $λ_\mathrm{R}$ and the closely related $(V/σ)$ 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.

preprint2020arXiv

Stellar angular momentum distribution linked to galaxy morphology

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.

preprint2019arXiv

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

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, $β$, of $\sim 0.046$ for the droplets, and when including both droplets and larger-scale dense cores, we find $β\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.

preprint2019arXiv

High-Precision Dark Halo Virial Masses from Globular Cluster Numbers: Implications for Globular Cluster Formation and Galaxy Assembly

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

preprint2019arXiv

On the origin of magnetic driven winds and the structure of the galactic dynamo in isolated 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_α$, $R_ω$ and $D$ we present strong evidence for an $α^2$-$Ω$ 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.

preprint2019arXiv

Stellar winds pump the heart of the Milky Way

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.

preprint2016arXiv

A Census of Large-Scale ($\ge$ 10 pc), Velocity-Coherent, Dense Filaments in the Northern Galactic Plane: Automated Identification Using Minimum Spanning Tree

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.

preprint2016arXiv

Clusters of Small Clumps Can Explain The Peculiar Properties of Giant Clumps in High-Redshift Galaxies

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 $σ_{\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.

preprint2016arXiv

Constructing multi-scale gravitational energy spectra from molecular cloud surface density PDF -- Interplay between turbulence and gravity

(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 $ρ(r)$. For a region with $N_{\rm col} \sim N^{-γ_{\rm N}}$, the gravitational energy spectra is $E_{\rm p}(k)\sim k^{-4(1 - 1/γ_{\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 $ρ\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.

preprint2016arXiv

Dynamical friction and scratches of orbiting satellite galaxies on host systems

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.

preprint2016arXiv

Instability of Supersonic Cold Streams Feeding Galaxies I: Linear Kelvin-Helmholtz Instability with Body Modes

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.

preprint2016arXiv

Probing the multi-scale interplay between gravity and turbulence - Power-law like gravitational energy spectra of the Orion Complex

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.

preprint2016arXiv

The evolution of metallicity and metallicity gradients from z=2.7-0.6 with KMOS3D

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.

preprint2015arXiv

26Al kinematics: superbubbles following the spiral arms? Constraints from the statistics of star clusters and HI supershells

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

preprint2015arXiv

Squeezed between shells? On the origin of the Lupus I molecular cloud. APEX/LABOCA, Herschel, and Planck observations

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.

preprint2015arXiv

Structure Formation in Gas-Rich Galactic Discs with Finite Thickness: From Discs to Rings

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.

preprint2015arXiv

The chemical evolution of local star forming galaxies: Radial profiles of ISM metallicity, gas mass, and stellar mass and constraints on galactic accretion and winds

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.

preprint2015arXiv

The co-planarity of satellite galaxies delivered by randomly aligned cold mode accretion streams

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.

preprint2015arXiv

The Structure and Dark Halo Core Properties of Dwarf Spheroidal Galaxies

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.

preprint2014arXiv

A Consistent Study of Metallicity Evolution at 0.8 < z < 2.6

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.

preprint2014arXiv

Cosmological simulations of black hole growth: AGN luminosities and downsizing

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)

preprint2014arXiv

Development of multiple tidal tails around globular clusters and dwarf satellite galaxies

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.

preprint2014arXiv

First Investigation of the Combined Impact of Ionizing Radiation and Momentum Winds from a Massive Star on a Self-Gravitating Core

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.

preprint2014arXiv

From Filaments to Oscillating Starless Cores

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.

preprint2014arXiv

Re-examining the Too-Big-To-Fail Problem for Dark Matter Haloes with Central Density Cores

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

preprint2014arXiv

The Bones of the Milky Way

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.

preprint2014arXiv

The connection between the cusp-to-core transformation and observational universalities of DM halos

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 $μ_{\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 $μ_{\rm 0D}$ relation and that Strigari relation follows from the $μ_{\rm 0D}$ relation for dwarf galaxies. We also show that the central densities of cored dark halos provide valuable information about their formation redshifts.

preprint2014arXiv

The formation of filamentary bundles in turbulent molecular clouds

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.

preprint2014arXiv

The Galactic Center cloud G2 and its gas streamer

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.

preprint2013arXiv

Balance among gravitational instability, star formation, and accretion determines the structure and evolution of disk galaxies

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.

preprint2013arXiv

Galactic winds - How to launch galactic outflows in typical Lyman-break galaxies

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.

preprint2013arXiv

Kinematics of massive star ejecta in the Milky Way as traced by 26Al

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.

preprint2013arXiv

Observations of the gas cloud G2 in the Galactic Center

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.

preprint2013arXiv

On the origin of the fundamental metallicity relation and the scatter in galaxy scaling relations

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.

preprint2013arXiv

Pericenter passage of the gas cloud G2 in the Galactic Center

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-γ luminosity of the head has been constant over the past 9 years, to within +- 25%, as have the line ratios Brackett-γ / Paschen-α and Brackett-γ / 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.

preprint2013arXiv

The Dark Halo - Spheroid Conspiracy and the Origin of Elliptical Galaxies

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.

preprint2013arXiv

The Dependence of Star Formation Efficiency on Gas Surface Density

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.

preprint2013arXiv

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

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.

preprint2013arXiv

The Mysterious Sickle Object in the Carina Nebula: A stellar wind induced bow shock grazing a clump?

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.

preprint2013arXiv

The Mystery of the $σ$-Bump -- A new Signature for Major Mergers in Early-type Galaxies?

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 $σ\propto r^β$. 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 "$σ$-"bump. To understand the origin of this $σ$-bump, we study a set of simulated early-type galaxies formed in major mergers. We find the $σ$-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 $σ$-bump stays nearly constant with radius and is a permanent feature that is preserved until the end of the simulation. The $σ$-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.

preprint2013arXiv

The Return of the Proplyds - Understanding the Dynamics of Ionization Triggered Stars

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.

preprint2013arXiv

The SINS/zC-SINF Survey of z~2 Galaxy Kinematics: The Nature of Dispersion Dominated Galaxies

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.

preprint2013arXiv

Wet Disc Contraction to Galactic Blue Nuggets and Quenching to Red Nuggets

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.

preprint2012arXiv

Evolving Gravitationally Unstable Disks Over Cosmic Time: Implications For Thick Disk Formation

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.

preprint2012arXiv

Feedback by massive stars and the emergence of superbubbles I. Energy efficiency & Vishniac instabilities

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.

preprint2012arXiv

Galacto-forensic of LMC's orbital history as a probe for the dark matter potential in the outskirt of the Galaxy

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.

preprint2012arXiv

Magnetohydrodynamic stability of broad line region clouds

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.

preprint2012arXiv

Origin of the anti-hierarchical growth of black holes

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.

preprint2012arXiv

Physics of the Galactic Center Cloud G2, on its Way towards the Super-Massive Black Hole

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.

preprint2012arXiv

Relaxation and Stripping: The Evolution of Sizes, Dispersions and Dark Matter Fractions in Major and Minor Mergers of Elliptical 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^α$, reaches values of $α\approx 2.4$, significantly exceeding the simple theoretical limit of $α= 2$. For major mergers the sizes grow with $α\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)

preprint2012arXiv

Shocked Superwinds from the z~2 Clumpy Star-forming Galaxy, ZC406690

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.

preprint2012arXiv

The SINS/zC-SINF survey of z~2 galaxy kinematics: Outflow properties

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

preprint2012arXiv

Threshing in Action - The tidal disruption of a dwarf galaxy by the Hydra I Cluster

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.

preprint2011arXiv

Delayed metal recycling in galaxies: the inefficiency of cold gas enrichment in colliding supershell simulations

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.

preprint2011arXiv

Radiation feedback on dusty clouds during Seyfert activity

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.

preprint2011arXiv

Rotational Support of Giant Clumps in High-z Disc Galaxies

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.

preprint2011arXiv

Short-lived star-forming giant clumps in cosmological simulations of z~2 disks

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.

preprint2010arXiv

A correlation between central supermassive black holes and the globular cluster systems of early-type galaxies

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.

preprint2010arXiv

A Toomre-like stability criterion for the clumpy and turbulent interstellar medium

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.

preprint2010arXiv

Accretion-Driven Turbulence and the Transition to Global Instability in Young Galaxy Disks

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.

preprint2010arXiv

Detailed Numerical Simulations on the Formation of Pillars around HII-regions

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.

preprint2010arXiv

Formation of Cold Filamentary Structure from Wind Blown Superbubbles

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.

preprint2010arXiv

Modelling Shear Flows with SPH and Grid Based Methods

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.

preprint2010arXiv

On the Dynamics and Evolution of Gravitational Instability-Dominated Disks

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.

preprint2010arXiv

Radiatively enhanced elasticity and turbulence in clumpy tori of Active Galactic Nuclei

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.

preprint2010arXiv

Simulations of Direct Collisions of Gas Clouds with the Central Black Hole

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.

preprint2010arXiv

Stability of Cloud Orbits in the Broad Line Region of Active Galactic Nuclei

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.

preprint2010arXiv

The structure of gravitationally unstable gas-rich disk galaxies

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.

preprint2010arXiv

The Two Phases of Galaxy Formation

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.

preprint2009arXiv

Bulge n and B/T in High Mass Galaxies: Constraints on the Origin of Bulges in Hierarchical Models

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.

preprint2009arXiv

Effects of ram pressure on the gas distribution and star formation in the Large Magellanic Cloud

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.

preprint2009arXiv

The evolution of Black Hole scaling relations in galaxy mergers

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.

preprint2008arXiv

Specific Angular Momentum of Disc Merger Remnants and the $λ_R$-Parameter

We use two-dimensional kinematic maps of simulated binary disc mergers to investigate the $λ_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 $λ_R > 0.1$ and slow rotators $λ_R < 0.1$. Tests on merger remnants reveal that $λ_R$ is a robust indicator of the true angular momentum content in elliptical galaxies. We find the same range of $λ_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 $λ_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 $λ_R$ as high as our progenitor disc galaxies.

preprint2008arXiv

The Challenge of Modelling Galactic Disks

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.

preprint2008arXiv

The Connection between Orbits and Isophotal Shape in Elliptical Galaxies

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.

preprint2008arXiv

The Inevitable Future of the Starless Core Barnard 68

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.

preprint2007arXiv

The Virial Balance of Clumps and Cores in Molecular Clouds

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

preprint2006arXiv

The Separation/Period Gap in the Distribution of Extrasolar Planets around Stars with Masses M >= 1.2 Msolar

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.

preprint2004arXiv

Galaxy Formation and the Cosmological Angular Momentum Problem

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.

preprint2003arXiv

Major Mergers and the Origin of Elliptical Galaxies

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.

preprint2003arXiv

On the Surface Heating of Synchronously-Spinning Short-Period Jovian Planets

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.

preprint2001arXiv

The Formation of the Milky Way in the Cosmological Context

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.

preprint2000arXiv

Star Formation-Regulated Growth of Black Holes in Protogalactic Spheroids

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.

preprint2000arXiv

The Effect of Gas Loss on the Formation of Bound Stellar Clusters

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

preprint1995arXiv

Constraints on massive black holes as dark matter candidates using Galactic globular clusters

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

Andreas Burkert

What is connected

Connect this record

See the researcher in context

Building this map preview

96 published item(s)

NOEMA$^{\rm 3D}$: A first kpc resolution study of a $z\sim1.5$ main sequence barred galaxy channeling gas into a growing bulge

Taking off the edge -- simultaneous filament and end core formation

An Interferometric View of H-MM1. I. Direct Observation of NH3 Depletion

Characterizing the 3D Kinematics of Young Stars in the Radcliffe Wave

Star formation near the Sun is driven by expansion of the Local Bubble

Hierarchical fragmentation in high redshift galaxies revealed by hydrodynamical simulations

3D simulations of clump formation in stellar wind collisions

Accretion driven turbulence in filaments II: Effects of self-gravity

Dwarfs in the Milky Way halo outer rim: first in-fall or backsplash satellites?

Efficient methanol production on the dark side of a prestellar core

Globular Cluster Formation from Colliding Substructure

Kinematics of Simulated Galaxies II: Probing the Stellar Kinematics of Galaxies out to Large Radii

Stellar angular momentum distribution linked to galaxy morphology

Droplets II: Internal Velocity Structures and Potential Rotational Motions in Pressure-dominated Coherent Structures

High-Precision Dark Halo Virial Masses from Globular Cluster Numbers: Implications for Globular Cluster Formation and Galaxy Assembly

On the origin of magnetic driven winds and the structure of the galactic dynamo in isolated galaxies

Stellar winds pump the heart of the Milky Way

A Census of Large-Scale ($\ge$ 10 pc), Velocity-Coherent, Dense Filaments in the Northern Galactic Plane: Automated Identification Using Minimum Spanning Tree

Clusters of Small Clumps Can Explain The Peculiar Properties of Giant Clumps in High-Redshift Galaxies

Constructing multi-scale gravitational energy spectra from molecular cloud surface density PDF -- Interplay between turbulence and gravity

Dynamical friction and scratches of orbiting satellite galaxies on host systems

Instability of Supersonic Cold Streams Feeding Galaxies I: Linear Kelvin-Helmholtz Instability with Body Modes

Probing the multi-scale interplay between gravity and turbulence - Power-law like gravitational energy spectra of the Orion Complex

The evolution of metallicity and metallicity gradients from z=2.7-0.6 with KMOS3D

26Al kinematics: superbubbles following the spiral arms? Constraints from the statistics of star clusters and HI supershells

Squeezed between shells? On the origin of the Lupus I molecular cloud. APEX/LABOCA, Herschel, and Planck observations

Structure Formation in Gas-Rich Galactic Discs with Finite Thickness: From Discs to Rings

The chemical evolution of local star forming galaxies: Radial profiles of ISM metallicity, gas mass, and stellar mass and constraints on galactic accretion and winds

The co-planarity of satellite galaxies delivered by randomly aligned cold mode accretion streams

The Structure and Dark Halo Core Properties of Dwarf Spheroidal Galaxies

A Consistent Study of Metallicity Evolution at 0.8 < z < 2.6

Cosmological simulations of black hole growth: AGN luminosities and downsizing

Development of multiple tidal tails around globular clusters and dwarf satellite galaxies

First Investigation of the Combined Impact of Ionizing Radiation and Momentum Winds from a Massive Star on a Self-Gravitating Core

From Filaments to Oscillating Starless Cores

Re-examining the Too-Big-To-Fail Problem for Dark Matter Haloes with Central Density Cores

The Bones of the Milky Way

The connection between the cusp-to-core transformation and observational universalities of DM halos

The formation of filamentary bundles in turbulent molecular clouds

The Galactic Center cloud G2 and its gas streamer

Balance among gravitational instability, star formation, and accretion determines the structure and evolution of disk galaxies

Galactic winds - How to launch galactic outflows in typical Lyman-break galaxies

Kinematics of massive star ejecta in the Milky Way as traced by 26Al

Observations of the gas cloud G2 in the Galactic Center

On the origin of the fundamental metallicity relation and the scatter in galaxy scaling relations

Pericenter passage of the gas cloud G2 in the Galactic Center

The Dark Halo - Spheroid Conspiracy and the Origin of Elliptical Galaxies

The Dependence of Star Formation Efficiency on Gas Surface Density

The effect of metal enrichment and galactic winds on galaxy formation in cosmological zoom simulations

The Mysterious Sickle Object in the Carina Nebula: A stellar wind induced bow shock grazing a clump?

The Mystery of the $σ$-Bump -- A new Signature for Major Mergers in Early-type Galaxies?

The Return of the Proplyds - Understanding the Dynamics of Ionization Triggered Stars

The SINS/zC-SINF Survey of z~2 Galaxy Kinematics: The Nature of Dispersion Dominated Galaxies

Wet Disc Contraction to Galactic Blue Nuggets and Quenching to Red Nuggets

Evolving Gravitationally Unstable Disks Over Cosmic Time: Implications For Thick Disk Formation

Feedback by massive stars and the emergence of superbubbles I. Energy efficiency & Vishniac instabilities

Galacto-forensic of LMC's orbital history as a probe for the dark matter potential in the outskirt of the Galaxy

Magnetohydrodynamic stability of broad line region clouds

Origin of the anti-hierarchical growth of black holes

Physics of the Galactic Center Cloud G2, on its Way towards the Super-Massive Black Hole

Relaxation and Stripping: The Evolution of Sizes, Dispersions and Dark Matter Fractions in Major and Minor Mergers of Elliptical Galaxies

Shocked Superwinds from the z~2 Clumpy Star-forming Galaxy, ZC406690

The SINS/zC-SINF survey of z~2 galaxy kinematics: Outflow properties

Threshing in Action - The tidal disruption of a dwarf galaxy by the Hydra I Cluster

Delayed metal recycling in galaxies: the inefficiency of cold gas enrichment in colliding supershell simulations

Radiation feedback on dusty clouds during Seyfert activity

Rotational Support of Giant Clumps in High-z Disc Galaxies

Short-lived star-forming giant clumps in cosmological simulations of z~2 disks

A correlation between central supermassive black holes and the globular cluster systems of early-type galaxies

A Toomre-like stability criterion for the clumpy and turbulent interstellar medium

Accretion-Driven Turbulence and the Transition to Global Instability in Young Galaxy Disks

Detailed Numerical Simulations on the Formation of Pillars around HII-regions

Formation of Cold Filamentary Structure from Wind Blown Superbubbles

Modelling Shear Flows with SPH and Grid Based Methods