Source author record

Marc Schartmann

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

18works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

18 published item(s)

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

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.

preprint2019arXiv

LLAMA: The $M_{BH}$ - $σ_{\star}$ Relation of the most luminous local AGNs

The $M_{BH}$ - $σ_{\star}$ relation is considered a result of co-evolution between the host galaxies and their super-massive black holes. For elliptical-bulge hosting inactive galaxies, this relation is well established, but there is still a debate whether active galaxies follow the same relation. In this paper, we estimate black hole masses for a sample of 19 local luminous AGNs (LLAMA) in order to test their location on the $M_{BH}$ - $σ_{\star}$ relation. Super-massive black hole masses ($M_{BH}$) were derived from the broad-line based relations for H$α$, H$β$ and Pa$β$ emission line profiles for the Type 1 AGNs. We compare the bulge stellar velocity dispersion ($σ_{\star}$) as determined from the Ca II triplet (CaT) with the dispersion measured from the near-infrared CO (2-0) absorption features for each AGN and find them to be consistent with each other. We apply an extinction correction to the observed broad line fluxes and we correct the stellar velocity dispersion by an average rotation contribution as determined from spatially resolved stellar kinematic maps. The H$α$-based black hole masses of our sample of AGNs were estimated in the range 6.34 $\leq$ $\log{M_{BH}}$ $\leq$ 7.75 M$_\odot$ and the $σ_{\star CaT}$ estimates range between 73 $\leq$ $σ_{\star CaT}$ $\leq$ 227 km s$^{-1}$. From the so-constructed $M_{BH}$ - $σ_{\star}$ relation for our Type 1 AGNs, we estimate the black hole masses for the Type 2 AGNs and the inactive galaxies in our sample. In conclusion, we find that our sample of local luminous AGNs is consistent with the $M_{BH}$ - $σ_{\star}$ relation of lower luminosity AGNs and inactive galaxies, after correcting for dust extinction and the rotational contribution to the stellar velocity dispersion.

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

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

Multi-phase Nature of a Radiation-Driven Fountain with Nuclear Starburst in a Low-mass Active Galactic Nucleus

The structures and dynamics of molecular, atomic, and ionized gases are studied around a low-luminosity active galactic nucleus (AGN) with a small ($2\times 10^6 M_\odot$) black hole using 3D radiation hydrodynamic simulations. We studied, for the first time, the non-equilibrium chemistry for the X-ray dominated region in the "radiation-driven fountain" (Wada 2012) with supernova feedback. A double hollow cone structure is naturally formed without postulating a thick "torus" around a central source. The cone is occupied with an inhomogeneous, diffuse ionized gas and surrounded by a geometrically thick ($h/r \gtrsim 1$) atomic gas. Dense molecular gases are distributed near the equatorial plane, and energy feedback from supernovae enhances their scale height. Molecular hydrogen exists in a hot phase ( > 1000 K) as well as in a cold ( < 100 K), dense ( >$10^3$ cm$^{-3}$) phase. The velocity dispersion of H$_2$ in the vertical direction is comparable to the rotational velocity, which is consistent with near infrared observations of nearby Seyfert galaxies. Using 3D radiation transfer calculations for the dust emission, we find polar emission in the mid-infrared band (12$μm$), which is associated with bipolar outflows, as suggested in recent interferometric observations of nearby AGNs. If the viewing angle for the nucleus is larger than 75 deg, the spectral energy distribution (~ 2 -- 60 $μm$) of this model is consistent with that of the Circinus galaxy. The multi-phase interstellar medium observed in optical/infrared and X-ray observations is also discussed.

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.

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

A diversity of dusty AGN tori: Data release for the VLTI/MIDI AGN Large Program and first results for 23 galaxies

The AGN-heated dust distribution (the "torus") is increasingly recognized not only as the absorber required in unifying models, but as a tracer for the reservoir that feeds the nuclear Super-Massive Black Hole. Yet, even its most basic structural properties (such as its extent, geometry and elongation) are unknown for all but a few archetypal objects. Since most AGNs are unresolved in the mid-infrared, we utilize the MID-infrared interferometric Instrument (MIDI) at the Very Large Telescope Interferometer (VLTI) that is sensitive to structures as small as a few milli-arcseconds (mas). We present here an extensive amount of new interferometric observations from the MIDI AGN Large Program (2009 - 2011) and add data from the archive to give a complete view of the existing MIDI observations of AGNs. Additionally, we have obtained high-quality mid-infrared spectra from VLT/VISIR. We present correlated and total flux spectra for 23 AGNs and derive flux and size estimates at 12 micron using simple axisymmetric geometrical models. Perhaps the most surprising result is the relatively high level of unresolved flux and its large scatter: The median "point source fraction" is 70 % for type 1 and 47 % for type 2 AGNs meaning that a large part of the flux is concentrated on scales smaller than about 5 mas (0.1 - 10 pc). Among sources observed with similar spatial resolution, it varies from 20 % - 100 %. For 18 of the sources, two nuclear components can be distinguished in the radial fits. While these models provide good fits to all but the brightest sources, significant elongations are detected in eight sources. The half-light radii of the fainter sources are smaller than expected from the size ~ L^0.5 scaling of the bright sources and show a large scatter, especially when compared to the relatively tight size--luminosity relation in the near-infrared.

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

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.

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.

preprint2011arXiv

On the size-luminosity relation of AGN dust tori in the mid-infrared

Interferometric measurements in the mid-infrared have shown that the sizes of the warm dust distributions in active galactic nuclei are consistent with their scaling with the square root of their luminosity. We carry out a more detailed analysis of this size-luminosity relation to investigate which of the general properties of the dusty tori in active galactic nuclei can be derived from this relation. We improve the accuracy of the size-luminosity relation by adding a few additional size measurements from more recent interferometric observations and compare the measured sizes to those derived from hydrodynamical and radiative transfer models of AGN tori. We find that a Gaussian approximation yields a reasonable estimate of the size of the brightness distribution, as long as the visibilities are within 0.2 {\leg} V {\leg} 0.9. The size estimates derived from the models are consistent with those determined from the measurements. However, the models predict a significant offset between the sizes derived for face-on and edge-on tori. This offset is not observed in the current data, probably because of the large uncertainties and low statistics of the present interferometric measurements. Furthermore, we find a ratio of the mid- to near-infrared sizes of approximately 30, whereas the first probes the body of the torus and the second is an estimate of the inner rim. The size-luminosity relation of AGN tori using Gaussian size estimates is a very simple and effective tool to investigate the internal structure and geometry of AGN tori and obtain constraints on the differences between type 1 and type 2 AGN. However, to fully exploit the possibilities of investigating the nuclear distributions of gas and dust in AGN using this size-luminosity relation, more accurate interferometric measurements of a larger sample of AGN are needed.

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.

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

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.