sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Schartmann, M.

Normalized to: Schartmann, M.

65 article(s) in total. 166 co-authors, from 1 to 37 common article(s). Median position in authors list is 3,0.

[1]  [pdf] - 2063856
Ionized outflows in local luminous AGN: what are the real densities and outflow rates?
Comments: Submitted to MNRAS; 25 pages (the last 5 are extra figures in an appendix)
Submitted: 2020-03-13
We report on the determination of electron densities, and their impact on the outflow masses and rates, measured in the central few hundred parsecs of 11 local luminous active galaxies. We show that the peak of the integrated line emission in the AGN is significantly offset from the systemic velocity as traced by the stellar absorption features, indicating that the profiles are dominated by outflow. In contrast, matched inactive galaxies are characterised by a systemic peak and weaker outflow wing. We present three independent estimates of the electron density in these AGN, discussing the merits of the different methods. The electron density derived from the [SII] doublet is significantly lower than than that found with a method developed in the last decade using auroral and transauroral lines, as well as a recently introduced method based on the ionization parameter. The reason is that, for gas photoionized by an AGN, much of the [SII] emission arises in an extended partially ionized zone where the implicit assumption that the electron density traces the hydrogen density is invalid. We propose ways to deal with this situation and we derive the associated outflow rates for ionized gas, which are in the range 0.001--0.5 M$_{\odot}$yr$^{-1}$ for our AGN sample. We compare these outflow rates to the relation between $\dot{M}_{out}$ and $L_{AGN}$ in the literature, and argue that it may need to be modified and rescaled towards lower mass outflow rates.
[2]  [pdf] - 2034364
3D simulations of clump formation in stellar wind collisions
Comments: 21 pages, 20 Figures, 2 Tables. Accepted for publication in MNRAS
Submitted: 2019-06-10, last modified: 2020-01-09
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.
[3]  [pdf] - 2050306
LLAMA: The $M_{BH}$ - $\sigma_{\star}$ Relation of the most luminous local AGNs
Comments: Accepted by A&A, 30 pages, 19 figures
Submitted: 2019-12-16
The $M_{BH}$ - $\sigma_{\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}$ - $\sigma_{\star}$ relation. Super-massive black hole masses ($M_{BH}$) were derived from the broad-line based relations for H$\alpha$, H$\beta$ and Pa$\beta$ emission line profiles for the Type 1 AGNs. We compare the bulge stellar velocity dispersion ($\sigma_{\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$\alpha$-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 $\sigma_{\star CaT}$ estimates range between 73 $\leq$ $\sigma_{\star CaT}$ $\leq$ 227 km s$^{-1}$. From the so-constructed $M_{BH}$ - $\sigma_{\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}$ - $\sigma_{\star}$ relation of lower luminosity AGNs and inactive galaxies, after correcting for dust extinction and the rotational contribution to the stellar velocity dispersion.
[4]  [pdf] - 2026019
Stellar winds pump the heart of the Milky Way
Comments: 6 pages and 5 figures. Accepted for publication in ApJL
Submitted: 2019-10-15, last modified: 2019-12-04
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.
[5]  [pdf] - 2042294
An image of the dust sublimation region in the nucleus of NGC 1068
Comments: Accepted for publication in Astronomy and Astrophysics. 13 pages with 7 figures
Submitted: 2019-12-03
We present near-infrared interferometric data on the Seyfert 2 galaxy NGC 1068, obtained with the GRAVITY instrument on the European Southern Observatory Very Large Telescope Interferometer. The extensive baseline coverage from 5 to 60 M\lambda allowed us to reconstruct a continuum image of the nucleus with an unrivaled 0.2 pc resolution in the K-band. We find a thin ring-like structure of emission with a radius r = 0.24+/-0.03 pc, inclination i = 70+/-5 deg, position angle PA = -50+/-4 deg, and h/r < 0.14, which we associate with the dust sublimation region. The observed morphology is inconsistent with the expected signatures of a geometrically and optically thick torus. Instead, the infrared emission shows a striking resemblance to the 22 GHz maser disc, which suggests they share a common region of origin. The near-infrared spectral energy distribution indicates a bolometric luminosity of (0.4-4.7) x 10^45 erg/s, behind a large A_K ~ 5.5 (A_V ~ 90) screen of extinction that also appears to contribute significantly to obscuring the broad line region.
[6]  [pdf] - 1983882
The multiphase gas structure and kinematics in the circumnuclear region of NGC 5728
Comments: accepted for publication in MNRAS; 22 pages + Appendix, 16 figures
Submitted: 2019-07-08, last modified: 2019-10-02
We report on our combined analysis of HST, VLT/MUSE, VLT/SINFONI, and ALMA observations of the local Seyfert 2 galaxy, NGC 5728 to investigate in detail the feeding and feedback of the AGN. The datasets simultaneously probe the morphology, excitation, and kinematics of the stars, ionized gas, and molecular gas over a large range of spatial scales (10 pc--10 kpc). NGC 5728 contains a large stellar bar which is driving gas along prominent dust lanes to the inner 1 kpc where the gas settles into a circumnuclear ring. The ring is strongly star forming and contains a substantial population of young stars as indicated by the lowered stellar velocity dispersion and gas excitation consistent with HII regions. We model the kinematics of the ring using the velocity field of the CO (2--1) emission and stars and find it is consistent with a rotating disk. The outer regions of the disk, where the dust lanes meet the ring, show signatures of inflow at a rate of 1 M$_{\sun}$ yr$^{-1}$. Inside the ring, we observe three molecular gas components corresponding to the circular rotation of the outer ring, a warped disk, and the nuclear stellar bar. The AGN is driving an ionized gas outflow that reaches a radius of 250 pc with a mass outflow rate of 0.08 M$_{\sun}$ yr$^{-1}$ consistent with its luminosity and scaling relations from previous studies. While we observe distinct holes in CO emission which could be signs of molecular gas removal, we find that largely the AGN is not disrupting the structure of the circumnuclear region.
[7]  [pdf] - 2065272
The resolved size and structure of hot dust in the immediate vicinity of AGN
Comments: 13 pages, 9 figures, submitted to A&A
Submitted: 2019-10-01
We use VLTI/GRAVITY near-infrared interferometry measurements of 8 bright, Type 1 AGN to study the size and structure of hot dust heated by the central engine. We partially resolve each source, and report Gaussian FWHM sizes in the range 0.3-0.8 milliarcseconds. In all but one object, we find no evidence for significant elongation or asymmetry (closure phases < 1 deg). The effective physical radius increases with bolometric luminosity as found from past reverberation and interferometry measurements. The measured sizes for Seyfert galaxies are systematically larger than for the two quasars in our sample when measured relative to the previously reported R ~ L^1/2 relationship explained by emission at the sublimation radius. This could be evidence of evolving near-infrared emission region structure as a function of central luminosity.
[8]  [pdf] - 1938474
A nuclear molecular outflow in the Seyfert galaxy NGC3227
Comments: Accepted for publication in Astronomy and Astrophysics, 17 pages, 12 figures
Submitted: 2019-06-17
We present ALMA observations of the CO(2-1) and CO(3-2) molecular gas transitions and associated (sub)-mm continua of the nearby Seyfert 1.5 galaxy NGC3227 with angular resolutions 0.085-0.21" (7-15pc). On large scales the cold molecular gas shows circular motions as well as streaming motions on scales of a few hundred parsecs associated with a large scale bar. We fitted the nuclear ALMA 1.3mm emission with an unresolved component and an extended component. The 850$\mu$m emission shows at least two extended components, one along the major axis of the nuclear disk and the other along the axis of the ionization cone. The molecular gas in the central region (1" ~73pc) shows several CO clumps with complex kinematics which appears to be dominated by non-circular motions. While we cannot demonstrate conclusively the presence of a warped nuclear disk, we also detected non-circular motions along the kinematic minor axis. They reach line-of-sight velocities of v-vsys =150-200km/s. Assuming that the radial motions are in the plane of the galaxy, then we interpret them as a nuclear molecular outflow due to molecular gas in the host galaxy being entrained by the AGN wind. We derive molecular outflow rates of $5\,M_\odot\,{\rm yr}^{-1}$ and $0.6\,M_\odot\,{\rm yr}^{-1}$ at projected distances of up to 30pc to the northeast and southwest of the AGN, respectively. At the AGN location we estimate a mass in molecular gas of $5\times 10^{5}\,M_\odot$ and an average column density $N({\rm H}_2) = 2-3\times 10^{23}\,{\rm cm}^{-2}$ in the inner 15pc. The nuclear molecular gas and sub-mm continuum emission of NGC3227 do not resemble the classical compact torus. Rather, these emissions extend for several tens of parsecs and appear connected with the circumnuclear ring in the host galaxy disk, as found in other local AGN. (Abridged)
[9]  [pdf] - 1846941
From kpcs to the central parsec of NGC 1097: feeding star formation and a black hole at the same time
Comments: 16 pages, 9 color figures. To appear in MNRAS
Submitted: 2019-02-25
A panchromatic view of the star forming ring and feeding process in the central kpc of the galaxy NGC 1097 is presented. The assembled IR to UV images at $\sim$10 pc resolution allow us to characterise the population of circa 250 clusters in the ring and disentangle the network of filaments of dust and gas that enshroud and feed them. The ring is a place of intermittent star bursts over the last 100 Myr. Four major episodes covering a proto-cluster phase of eleven mid-IR sources at the molecular clouds core, and two (three) previous bursts with a time separation of 20 - 30 Myr are identified. The extinction map of the inner few kpc resolves NGC1097's two major dust lanes in bundles of narrow, $<$25 pc width, filaments running along the galaxy's bar. As they approach the ring, some circularise along it, others curve to the centre to produce a nuclear spiral. We believe these are kpc-scale dust-gas streamers feeding the ring and the black-hole. The total mass in clusters formed in the ring in the last 100 Myr is $< 10^7\, \rm{M_\odot}$, i.e. $< 1\% $ of the $10^{9} M_\odot$ of molecular gas in the ring; yet, at its current star formation rate, $\sim1.8\, \rm{M_\odot \, yr^{-1}}$, an order of magnitude more in stellar mass should have been produced over that period. This means that the availability of gas in the ring is not the sole star formation driver, perhaps the rate at which dense gas accumulates in the ring is key.
[10]  [pdf] - 1791986
Detection of a drag force in G2's orbit: Measuring the density of the accretion flow onto Sgr A* at 1000 Schwarzschild radii
Comments: accepted by ApJ, 15 pages, 6 figures
Submitted: 2018-12-04
The Galactic Center black hole Sgr A* is the archetypical example of an underfed massive black hole. The extremely low accretion rate can be understood in radiatively inefficient accretion flow models. Testing those models has proven to be difficult due to the lack of suitable probes. Radio and submm polarization measurements constrain the flow very close to the event horizon. X-ray observations resolving the Bondi radius yield an estimate roughly four orders of magnitude further out. Here, we present a new, indirect measurement of the accretion flow density at intermediate radii. We use the dynamics of the gas cloud G2 to probe the ambient density. We detect the presence of a drag force slowing down G2 with a statistical significance of approx 9 {\sigma}. This probes the accretion flow density at around 1000 Schwarzschild radii and yields a number density of approx. 4 x 10^3 cm^-3. Self-similar accretion models where the density follows a power law radial profile between the inner zone and the Bondi radius have predicted similar values.
[11]  [pdf] - 1736339
Simulating the pericentre passage of the Galactic centre star S2
Comments: accepted by A&A Letters
Submitted: 2018-04-09, last modified: 2018-08-01
The so-called S2 star reached its closest approach to the massive black hole (BH) at around 1500 $R_\mathrm{s}$ in May 2018. It has been proposed that the interaction of its stellar wind with the high-density accretion flow at this distance from Sgr A* will lead to a detectable, month-long X-ray flare. Our goal is to verify whether or not the S2 star wind can be used as a diagnostic tool to infer the properties of the accretion flow towards Sgr A* at its pericentre (an unprobed distance regime), putting important constraints on BH accretion flow models. We run a series of three-dimensional adaptive mesh refinement simulations with the help of the Ramses code which include the realistic treatment of the interaction of S2's stellar wind with the accretion flow along its orbit and - apart from hydrodynamical and thermodynamical effects - include the tidal interaction with the massive BH. These are post-processed to derive the X-ray emission in the observable 2-10 keV window. No significant excess of X-ray emission from Sgr A* is found for typical accretion flow models. A measurable excess is produced for a significantly increased density of the accretion flow. This can, however, be ruled out for standard power-law accretion flow models as in this case the thermal X-ray emission without the S2 wind interaction would already exceed the observed quiescent luminosity. Only a significant change of the wind parameters (increased mass loss rate and decreased wind velocity) might lead to an (marginally) observable X-ray flaring event. Even the detection of an (month-long) X-ray flare during the pericentre passage of the S2 star would not allow for strict constraints to be put on the accretion flow around Sgr A* due to the degeneracy caused by the dependence on multiple parameters (of the accretion flow model as well as the stellar wind).
[12]  [pdf] - 1694163
The Galactic Centre source G2 was unlikely born in any of the known massive binaries
Comments: 14 pages, 10 figures. Accepted for publication in MNRAS
Submitted: 2018-05-14, last modified: 2018-05-17
The source G2 has already completed its pericentre passage around Sgr A*, the super-massive black hole in the centre of our Galaxy. Although it has been monitored for 15 years, its astrophysical nature and origin still remain unknown. In this work, we aim to test the hypothesis of G2 being the result of a stellar wind collision. To do so, we study the motion and final fate of gas clumps formed as a result of collisions of stellar winds in massive binaries. Our approach is based on a test-particle model in order to describe the trajectories of such clumps. The model takes into account the gravitational field of Sgr A*, the interaction of the clumps with the interstellar medium as well as their finite lifetimes. Our analysis allows us to reject the hypothesis based on four arguments: i) if G2 has followed a purely Keplerian orbit since its formation, it cannot have been produced in any of the known massive binaries since their motions are not consistent; ii) in general, gas clumps are evaporated through thermal conduction on very short timescale (< 100yr) before getting close enough to Sgr A*; iii) IRS 16SW, the best candidate for the origin of G2, cannot generate clumps as massive as G2; and iv) clumps ejected from IRS 16SW describe trajectories significantly different to the observed motion of G2.
[13]  [pdf] - 1732608
Thick turbulent gas disks with magnetocentrifugal winds in active galactic nuclei - Model infrared emission and optical polarization
Comments: accepted for publication in A&A
Submitted: 2018-03-16
(Abridged) Infrared high-resolution imaging and interferometry have shown that the dust distribution is frequently elongated along the polar direction of an AGN. To explain these findings, we developed a model scenario for the inner ~30 pc of an AGN. We assume a rotating thick gas disk between about one and ten parsec. External gas accretion adds mass and injects energy via gas compression into this gas disk and drives turbulence. Our disks are assumed to be strongly magnetized via equipartition between the turbulent gas pressure and the energy density of the magnetic field. In a second step, we built three dimensional density cubes based on the analytical model, illuminated them with a central source, and made radiative transfer calculations. In a third step, we calculated MIR visibility amplitudes and compared them to available interferometric observations. We show that magnetocentrifugal winds starting from a thin and thick gas disk are viable in active galaxy centers. Once the wind is launched, it is responsible for the transport of angular momentum and the gas disk can become thin. The outflow scenario can account for the elongated dust structures, outer edges of the thin maser disks, and molecular outflows observed in local AGN. The models reproduce the observed terminal wind velocities, the scatter of the MIR/intrinsic X-ray correlation, and point source fractions. An application of the model to the Circinus Galaxy and NGC 1068 shows that the IR SED, available MIR interferometric observations, and optical polarization can be reproduced in a satisfactory way, provided that (i) a puff-up at the inner edge of the thin disk is present and (ii) a local screen with an optical depth of tau_V 20 in form of a local gas filament and/or a warp of the thick disk hide a significant fraction of both nuclei.
[14]  [pdf] - 1604978
LLAMA: Normal star formation efficiencies of molecular gas in the centres of luminous Seyfert galaxies
Comments: 24 pages, 18 figures, 2 part Appendix; accepted for publication in MNRAS
Submitted: 2017-10-11
Using new APEX and JCMT spectroscopy of the CO 2-1 line, we undertake a controlled study of cold molecular gas in moderately luminous Active Galactic Nuclei (AGN) and inactive galaxies from the Luminous Local AGN with Matched Analogs (LLAMA) survey. We use spatially resolved infrared photometry of the LLAMA galaxies from 2MASS, WISE, IRAS & Herschel, corrected for nuclear emission using multi-component spectral energy distribution (SED) fits, to examine the dust-reprocessed star-formation rates (SFRs), molecular gas fractions and star formation efficiencies (SFEs) over their central 1 - 3 kpc. We find that the gas fractions and central SFEs of both active and inactive galaxies are similar when controlling for host stellar mass and morphology (Hubble type). The equivalent central molecular gas depletion times are consistent with the discs of normal spiral galaxies in the local Universe. Despite energetic arguments that the AGN in LLAMA should be capable of disrupting the observable cold molecular gas in their central environments, our results indicate that nuclear radiation only couples weakly with this phase. We find a mild preference for obscured AGN to contain higher amounts of central molecular gas, which suggests a connection between AGN obscuration and the gaseous environment of the nucleus. Systems with depressed SFEs are not found among the LLAMA AGN. We speculate that the processes that sustain the collapse of molecular gas into dense pre-stellar cores may also be a prerequisite for the inflow of material on to AGN accretion disks.
[15]  [pdf] - 1589584
LLAMA: Nuclear stellar properties of Swift BAT AGN and matched inactive galaxies
Comments: 34 pages, 25 figures. Accepted for publication in MNRAS
Submitted: 2017-10-11
In a complete sample of local 14-195 keV selected AGNs and inactive galaxies, matched by their host galaxy properties, we study the spatially resolved stellar kinematics and luminosity distributions at near-infrared wavelengths on scales of 10-150 pc, using SINFONI on the VLT. In this paper, we present the first half of the sample, which comprises 13 galaxies, 8 AGNs and 5 inactive galaxies. The stellar velocity fields show a disk-like rotating pattern, for which the kinematic position angle is in agreement with the photometric position angle obtained from large scale images. For this set of galaxies, the stellar surface brightness of the inactive galaxy sample is generally comparable to the matched sample of AGN but extends to lower surface brightness. After removal of the bulge contribution, we find a nuclear stellar light excess with an extended nuclear disk structure, and which exhibits a size-luminosity relation. While we expect the excess luminosity to be associated with a dynamically cooler young stellar population, we do not typically see a matching drop in dispersion. This may be because these galaxies have pseudo-bulges in which the intrinsic dispersion increases towards the centre. And although the young stars may have an impact in the observed kinematics, their fraction is too small to dominate over the bulge and compensate the increase in dispersion at small radii, so no dispersion drop is seen. Finally, we find no evidence for a difference in the stellar kinematics and nuclear stellar luminosity excess between these active and inactive galaxies.
[16]  [pdf] - 1588193
The life cycle of starbursting circumnuclear gas discs
Comments: accepted by MNRAS
Submitted: 2017-09-11
High-resolution observations from the sub-mm to the optical wavelength regime resolve the central few 100pc region of nearby galaxies in great detail. They reveal a large diversity of features: thick gas and stellar discs, nuclear starbursts, in- and outflows, central activity, jet interaction, etc. Concentrating on the role circumnuclear discs play in the life cycles of galactic nuclei, we employ 3D adaptive mesh refinement hydrodynamical simulations with the RAMSES code to self-consistently trace the evolution from a quasi-stable gas disc, undergoing gravitational (Toomre) instability, the formation of clumps and stars and the disc's subsequent, partial dispersal via stellar feedback. Our approach builds upon the observational finding that many nearby Seyfert galaxies have undergone intense nuclear starbursts in their recent past and in many nearby sources star formation is concentrated in a handful of clumps on a few 100pc distant from the galactic centre. We show that such observations can be understood as the result of gravitational instabilities in dense circumnuclear discs. By comparing these simulations to available integral field unit observations of a sample of nearby galactic nuclei, we find consistent gas and stellar masses, kinematics, star formation and outflow properties. Important ingredients in the simulations are the self-consistent treatment of star formation and the dynamical evolution of the stellar distribution as well as the modelling of a delay time distribution for the supernova feedback. The knowledge of the resulting simulated density structure and kinematics on pc scale is vital for understanding inflow and feedback processes towards galactic scales.
[17]  [pdf] - 1587347
Polarized radiative transfer modeling of warped and clumpy dusty tori
Comments: 15 pages, 13 Figures, accepted for publication in A&A
Submitted: 2017-08-22
Context: Active galactic nuclei (AGN) are anisotropic objects surrounded by an optically thick equatorial medium whose true geometry still defy observers. Aims: In this paper, we aim to explore the optical, scattering-induced, polarization emerging from clumpy and warped dusty tori to check whether they can fit with the unified model predictions. Methods: We run polarized radiative transfer simulations in a set of warped and non warped clumpy tori to explore the differences induced by distorted dust distributions. We then include warped tori in a more complex model representative of an AGN to check, using polarimetry and imaging methods, if warps can reproduce the expected polarization dichotomy between Seyfert-1 and Seyfert-2 AGN. Results: The main results from our simulations highlight that isolated warped structures imprint the polarization degree and angle with distinctive signatures at Seyfert-1 orientations. Included in an AGN model, the signatures of warps are easily (but not always) washed out by multiple scattering in a clumpy environment. Imaging polarimetry may help to detect warped tori but we prove that warps can exist in AGN circumnuclear regions without contradicting observations. Conclusions: Two warped tori with a non significant difference in geometry in terms of photometry or spectroscopy can have totally different signatures in polarimetry. Testing the geometry of any alternative model to the usual dusty torus using polarized radiative transfer is a necessary approach to verify or reject a hypothesis.
[18]  [pdf] - 1697890
3D AMR hydrosimulations of a compact source scenario for the Galactic Centre cloud G2
Comments: resubmitted to MNRAS after referee report, 16 pages, 11 figures
Submitted: 2017-06-26
The nature of the gaseous and dusty cloud G2 in the Galactic Centre is still under debate. We present three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations of G2, modeled as an outflow from a "compact source" moving on the observed orbit. The construction of mock position-velocity (PV) diagrams enables a direct comparison with observations and allow us to conclude that the observational properties of the gaseous component of G2 could be matched by a massive ($\dot{M}_\mathrm{w}=5\times 10^{-7} \;M_{\odot} \mathrm{yr^{-1}}$) and slow ($50 \;\mathrm{km \;s^{-1}}$) outflow, as observed for T Tauri stars. In order for this to be true, only the material at larger ($>100 \;\mathrm{AU}$) distances from the source must be actually emitting, otherwise G2 would appear too compact compared to the observed PV diagrams. On the other hand, the presence of a central dusty source might be able to explain the compactness of G2's dust component. In the present scenario, 5-10 years after pericentre the compact source should decouple from the previously ejected material, due to the hydrodynamic interaction of the latter with the surrounding hot and dense atmosphere. In this case, a new outflow should form, ahead of the previous one, which would be the smoking gun evidence for an outflow scenario.
[19]  [pdf] - 1557938
IC 630: Piercing the Veil of the Nuclear Gas
Comments: 32 pages, 19 figures Accepted for publication in ApJ
Submitted: 2017-03-10
IC 630 is a nearby early-type galaxy with a mass of $6 \times 10^{10} M_{\odot}$ with an intense burst of recent (6 Myr) star formation. It shows strong nebular emission lines, with radio and X-ray emission, which classifies it as an AGN. With VLT-SINFONI and Gemini North-NIFS adaptive optics observations (plus supplementary ANU 2.3m WiFeS optical IFU observations), the excitation diagnostics of the nebular emission species show no sign of standard AGN engine excitation; the stellar velocity dispersion also indicate that a super-massive black hole (if one is present) is small ($M_{\bullet} = 2.25 \times 10^{5}~M_{\odot}$). The luminosity at all wavelengths is consistent with star formation at a rate of about $1-2 M_{\odot}$/yr. We measure gas outflows driven by star formation at a rate of $0.18 M_{\odot}$/yr in a face-on truncated cone geometry. We also observe a nuclear cluster or disk and other clusters. Photo-ionization from young, hot stars is the main excitation mechanism for [Fe II] and hydrogen, whereas shocks are responsible for the H$_2$ excitation. Our observations are broadly comparable with simulations where a Toomre-unstable, self-gravitating gas disk triggers a burst of star formation, peaking after about 30 Myr and possibly cycling with a period of about 200 Myr.
[20]  [pdf] - 1532342
The Role of Host Galaxy for the Environmental Dependence of Active Nuclei in Local Galaxies
Comments: accepted for publication in MNRAS. 12 pages with 10 figures
Submitted: 2016-10-31, last modified: 2017-01-06
We discuss the environment of local hard X-ray selected active galaxies, with reference to two independent group catalogues. We find that the fraction of these AGN in S0 host galaxies decreases strongly as a function of galaxy group size (halo mass) - which contrasts with the increasing fraction of galaxies of S0 type in denser environments. However, there is no evidence for an environmental dependence of AGN in spiral galaxies. Because most AGN are found in spiral galaxies, this dilutes the signature of environmental dependence for the population as a whole. We argue that the differing results for AGN in disk-dominated and bulge-dominated galaxies is related to the source of the gas fuelling the AGN, and so may also impact the luminosity function, duty cycle, and obscuration. We find that there is a significant difference in the luminosity function for AGN in spiral and S0 galaxies, and tentative evidence for some difference in the fraction of obscured AGN.
[21]  [pdf] - 1532200
Radiative transfer modelling of parsec-scale dusty warped discs
Comments: accepted by MNRAS
Submitted: 2016-10-23
Warped discs have been found on (sub-)parsec scale in some nearby Seyfert nuclei, identified by their maser emission. Using dust radiative transfer simulations we explore their observational signatures in the infrared in order to find out whether they can partly replace the molecular torus. Strong variations of the brightness distributions are found, depending on the orientation of the warp with respect to the line of sight. Whereas images at short wavelengths typically show a disc-like and a point source component, the warp itself only becomes visible at far-infrared wavelengths. A similar variety is visible in the shapes of the spectral energy distributions. Especially for close to edge-on views, the models show silicate feature strengths ranging from deep absorption to strong emission for variations of the lines of sight towards the warp. To test the applicability of our model, we use the case of the Circinus galaxy, where infrared interferometry has revealed a highly elongated emission component matching a warped maser disc in orientation and size. Our model is for the first time able to present a physical explanation for the observed dust morphology as coming from the AGN heated dust. As opposed to available torus models, a warped disc morphology produces a variety of silicate feature shapes for grazing lines of sight, close to an edge-on view. This could be an attractive alternative to a claimed change of the dust composition for the case of the nearby Seyfert 2 galaxy NGC 1068, which harbours a warped maser disc as well.
[22]  [pdf] - 1531712
3D AMR simulations of G2 as an outflow
Comments: 2 pages, 1 figure, Proceedings of IAU Symposium 322: The Multi-Messenger Astrophysics of the Galactic Centre
Submitted: 2016-09-23
We study the evolution of G2 in a \textit{Compact Source Scenario}, where G2 is the outflow from a low-mass central star moving on the observed orbit. This is done through 3D AMR simulations of the hydrodynamic interaction of G2 with the surrounding hot accretion flow. A comparison with observations is done by means of mock position-velocity (PV) diagrams. We found that a massive ($\dot{M}_\mathrm{w}=5\times 10^{-7} \;M_{\odot} \; \mathrm{yr^{-1}}$) and slow ($v_\mathrm{w}=50 \;\mathrm{km\; s^{-1}}$) outflow can reproduce G2's properties. A faster outflow ($v_\mathrm{w}=400 \;\mathrm{km\; s^{-1}}$) might also be able to explain the material that seems to follow G2 on the same orbit.
[23]  [pdf] - 1531711
3D AMR simulations of the evolution of the diffuse gas cloud G2 in the Galactic Centre
Comments: Proceedings IAU Symposium No. 322, The Multi-Messenger Astrophysics of the Galactic Centre, 2016, S. Longmore, G. Bicknell & R. Crocker, eds
Submitted: 2016-09-23
With the help of 3D AMR hydrodynamical simulations we aim at understanding G2's nature, recent evolution and fate in the coming years. By exploring the possible parameter space of the diffuse cloud scenario, we find that a starting point within the disc of young stars is favoured by the observations, which may hint at G2 being the result of stellar wind interactions.
[24]  [pdf] - 1483503
Multi-phase Nature of a Radiation-Driven Fountain with Nuclear Starburst in a Low-mass Active Galactic Nucleus
Comments: 9 pages, 5 figures. Accepted for ApJL. A movie file for Fig.5b can be downloaded from
Submitted: 2016-08-24
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$\mu 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 $\mu 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.
[25]  [pdf] - 1475452
Constraints on the Broad Line Region Properties and Extinction in Local Seyferts
Comments: 19 pages, 14 figures, accepted for publication in Monthly Notices of the Royal Astronomical Society
Submitted: 2016-07-25
We use high spectral resolution (R > 8000) data covering 3800-13000\r{A} to study the physical conditions of the broad line region (BLR) of nine nearby Seyfert 1 galaxies. Up to six broad HI lines are present in each spectrum. A comparison - for the first time using simultaneous optical to near-infrared observations - to photoionisation calculations with our devised simple scheme yields the extinction to the BLR at the same time as determining the density and photon flux, and hence distance from the nucleus, of the emitting gas. This points to a typical density for the HI emitting gas of 10$^{11}$cm$^{-3}$ and shows that a significant amount of this gas lies at regions near the dust sublimation radius, consistent with theoretical predictions. We also confirm that in many objects the line ratios are far from case B, the best-fit intrinsic broad-line H$\alpha$/H$\beta$ ratios being in the range 2.5-6.6 as derived with our photoionization modeling scheme. The extinction to the BLR, based on independent estimates from HI and HeII lines, is A$_V$ $\le$ 3 for Seyfert 1-1.5s, while Seyfert 1.8-1.9s have A$_V$ in the range 4-8. A comparison of the extinction towards the BLR and narrow line region (NLR) indicates that the structure obscuring the BLR exists on scales smaller than the NLR. This could be the dusty torus, but dusty nuclear spirals or filaments could also be responsible. The ratios between the X-ray absorbing column N$_H$ and the extinction to the BLR are consistent with the Galactic gas-to-dust ratio if N$_H$ variations are considered.
[26]  [pdf] - 1378969
The G2+G2t complex as a fast and massive outflow?
Comments: 6 pages, 4 figures. Accepted for publication in ApJL
Submitted: 2016-02-23
Observations of the gas component of the cloud G2 in the Galactic Center have revealed its connection to a tail (G2t) lying on the same orbit. More recent studies indicate a connection between G2 and G1, another cloud detected on the blueshifted side of G2's orbit, suggesting a scenario in which G2 is a denser clump in a stream of gas. In this Letter we show that a simulation of an outflow by a central source (possibly a T Tauri star) moving on G2's orbit and interacting with a hot atmosphere surrounding SgrA* can have G2 and G2t as a byproduct. G2 would be the bow-shock formed in the head of the source, while G2t might be the result of the stripping of the rest of the shocked material by the ram pressure of the surrounding hot gas and of its successive accumulation in the trailing region. Mock position-velocity diagrams for the Br$\gamma$ emission for this simulation can indeed reproduce the correct position and velocity of G2t, as well as the more tenuous material in between. Though some tension between the observations and the simulated model remains, we argue that this might be due to issues in the construction of observed position-velocity (PV) diagrams and/or to a poor treatment of some physical processes - like hydrodynamic mixing - in our simulation.
[27]  [pdf] - 1422232
Mid-infrared interferometry of 23 AGN tori: On the significance of polar-elongated emission
Comments: 12 pages, 6 figures. Accepted for publication on A&A
Submitted: 2016-02-17
Context. Detailed high resolution studies of AGN with mid-infrared (MIR) interferometry have revealed parsec-sized dust emission elongated in the polar direction in four sources. Aims. Using a larger, coherently analyzed sample of AGN observed with MIR interferometry, we aim to identify elongated mid-infrared emission in a statistical sample of sources. More specifically we wish to determine if there is indeed a preferred direction of the elongation and whether this direction is consistent with a torus-like structure or with a polar emission. Methods. We investigate the significance of the detection of an elongated shape in the MIR emission by fitting elongated Gaussian models to the interferometric data at 12 um. We pay special attention to 1) the uncertainties caused by an inhomogeneous (u,v) coverage, 2) the typical errors in the measurements and 3) the spatial resolution achieved for each object. Results. From our sample of 23 sources we are able to find elongated parsec-scale MIR emission in five sources: three type 2s, one type 1i and one type 1. Elongated emission in four of these sources has been published before; NGC5506 is a new detection. The observed axis ratios are typically around 2 and the position angle of the 12 um emission for all the elongated sources seems to be always closer to the polar axis of the system than to the equatorial axis. Two other objects, NGC4507 and MCG-5-23-16 with a reasonably well mapped (u,v) coverage and good signal-to-noise ratios, appear to have a less elongated 12 um emission. Conclusions. Our finding that sources showing elongated mid-infrared emission are preferentially extended in polar direction sets strong constraints on torus models or implies that both the torus and the NLR/outflow region have to be modeled together. Especially also models used for SED fitting will have to be revised to include emission from polar dust.
[28]  [pdf] - 1374211
Clusters of Small Clumps Can Explain The Peculiar Properties of Giant Clumps in High-Redshift Galaxies
Comments: Modifications for more clarity! ApJL in press! Brings fresh blood into the giant clump paradigm!
Submitted: 2015-12-10, last modified: 2016-02-05
Giant clumps are a characteristic feature of observed high-redshift disk galaxies. We propose that these kpc-sized clumps have a complex substructure and are the result of many smaller clumps self-organizing themselves into clump clusters (CC). This bottom-up scenario is in contrast to the common top-down view that these giant clumps form first and then sub fragment. Using a high resolution hydrodynamical simulation of an isolated, fragmented massive gas disk and mimicking the observations from Genzel et al. (2011) at $z \sim 2$, we find remarkable agreement in many details. The CCs appear as single entities of sizes $R_{\mathrm{HWHM}} \simeq 0.9-1.4$ kpc and masses $\sim 1.5-3 \times 10^9 \ \mathrm{M_{\odot}}$ representative of high-z observations. They are organized in a ring around the center of the galaxy. The origin of the observed clumps' high intrinsic velocity dispersion $\sigma_{\mathrm{intrinsic}} \simeq 50 - 100 \ \mathrm{km \ s^{-1}}$ is fully explained by the internal irregular motions of their substructure in our simulation. No additional energy input, e.g. via stellar feedback, is necessary. Furthermore, in agreement with observations, we find a small velocity gradient $V_{\mathrm{grad}} \simeq 8 - 27 \ \mathrm{km \ s^{-1} \ kpc^{-1}}$ along the CCs in the beam smeared velocity residual maps which corresponds to net prograde and retrograde rotation with respect to the rotation of the galactic disk. The CC scenario could have strong implications for the internal evolution, lifetimes and the migration timescales of the observed giant clumps, bulge growth and AGN activity, stellar feedback and the chemical enrichment history of galactic disks.
[29]  [pdf] - 1331142
Stellar feedback efficiencies: supernovae versus stellar winds
Comments: The manuscript MN-15-2374-MJ.R2 has been accepted for publication in Monthly Notices of the Royal Astronomical Society Main Journal on 2015 November 16
Submitted: 2015-11-16, last modified: 2015-11-18
Stellar winds and supernova (SN) explosions of massive stars ("stellar feedback") create bubbles in the interstellar medium (ISM) and insert newly produced heavy elements and kinetic energy into their surroundings, possibly driving turbulence. Most of this energy is thermalized and immediately removed from the ISM by radiative cooling. The rest is available for driving ISM dynamics. In this work we estimate the amount of feedback energy retained as kinetic energy when the bubble walls have decelerated to the sound speed of the ambient medium. We show that the feedback of the most massive star outweighs the feedback from less massive stars. For a giant molecular cloud (GMC) mass of 1e5 solar masses (as e.g. found in the Orion GMCs) and a star formation efficiency of 8% the initial mass function predicts a most massive star of approximately 60 solar masses. For this stellar evolution model we test the dependence of the retained kinetic energy of the cold GMC gas on the inclusion of stellar winds. In our model winds insert 2.34 times the energy of a SN and create stellar wind bubbles serving as pressure reservoirs. We find that during the pressure driven phases of the bubble evolution radiative losses peak near the contact discontinuity (CD), and thus, the retained energy depends critically on the scales of the mixing processes across the CD. Taking into account the winds of massive stars increases the amount of kinetic energy deposited in the cold ISM from 0.1% to a few percent of the feedback energy.
[30]  [pdf] - 1327446
Clump formation through colliding stellar winds in the Galactic Centre
Comments: Accepted for publication in MNRAS, 11 pages, 6 figures
Submitted: 2015-07-24, last modified: 2015-11-11
The gas cloud G2 is currently being tidally disrupted by the Galactic Centre super-massive black hole, Sgr A*. The region around the black hole is populated by $\sim 30$ Wolf-Rayet stars, which produce strong outflows. We explore the possibility that gas clumps, such as G2, originate from the collision of stellar winds via the non-linear thin shell instability. Following an analytical approach, we study the thermal evolution of slabs formed in the symmetric collision of winds, evaluating whether instabilities occur, and estimating possible clump masses. We find that the collision of relatively slow ($< 750$ km s$^{-1}$) and strong ($\sim 10^{-5}$ Msun yr$^{-1}$) stellar winds from stars at short separations ($<10$ mpc) is a process that indeed could produce clumps of G2's mass and above. Such short separation encounters of single stars along their known orbits are not common in the Galactic Centre, making this process a possible but unlikely origin for G2. We also discuss clump formation in close binaries such as IRS 16SW and in asymmetric encounters as promising alternatives that deserve further numerical study.
[31]  [pdf] - 1292296
3D adaptive mesh refinement simulations of the gas cloud G2 born within the disks of young stars in the Galactic Center
Comments: resubmitted to ApJ after referee report, 13 pages, 11 figures
Submitted: 2015-08-27
The dusty, ionized gas cloud G2 is currently passing the massive black hole in the Galactic Center at a distance of roughly 2400 Schwarzschild radii. We explore the possibility of a starting point of the cloud within the disks of young stars. We make use of the large amount of new observations in order to put constraints on G2's origin. Interpreting the observations as a diffuse cloud of gas, we employ three-dimensional hydrodynamical adaptive mesh refinement (AMR) simulations with the PLUTO code and do a detailed comparison with observational data. The simulations presented in this work update our previously obtained results in multiple ways: (1) high resolution three-dimensional hydrodynamical AMR simulations are used, (2) the cloud follows the updated orbit based on the Brackett-$\gamma$ data, (3) a detailed comparison to the observed high-quality position-velocity diagrams and the evolution of the total Brackett-$\gamma$ luminosity is done. We concentrate on two unsolved problems of the diffuse cloud scenario: the unphysical formation epoch only shortly before the first detection and the too steep Brackett-$\gamma$ light curve obtained in simulations, whereas the observations indicate a constant Brackett-$\gamma$ luminosity between 2004 and 2013. For a given atmosphere and cloud mass, we find a consistent model that can explain both, the observed Brackett-$\gamma$ light curve and the position-velocity diagrams of all epochs. Assuming initial pressure equilibrium with the atmosphere, this can be reached for a starting date earlier than roughly 1900, which is close to apo-center and well within the disks of young stars.
[32]  [pdf] - 987922
Insights on the Dusty Torus and Neutral Torus from Optical and X-ray Obscuration in a Complete Volume Limited Hard X-ray AGN Sample
Comments: Accepted for ApJ; 15 pages with 10 figures
Submitted: 2015-05-04
We describe a complete volume limited sample of nearby active galaxies selected by their 14-195keV luminosity, and outline its rationale for studying the mechanisms regulating gas inflow and outflow. We describe also a complementary sample of inactive galaxies, selected to match the AGN host galaxy properties. The active sample appears to have no bias in terms of AGN type, the only difference being the neutral absorbing column which is two orders of magnitude greater for the Seyfert 2s. In the luminosity range spanned by the sample, log L_{14-195keV} [erg/s] = 42.4-43.7, the optically obscured and X-ray absorbed fractions are 50-65%. The similarity of these fractions to more distant spectroscopic AGN samples, although over a limited luminosity range, suggests that the torus does not strongly evolve with redshift. Our sample confirms that X-ray unabsorbed Seyfert 2s are rare, comprising not more than a few percent of the Seyfert 2 population. At higher luminosities, the optically obscured fraction decreases (as expected for the increasing dust sublimation radius), but the X-ray absorbed fraction changes little. We argue that the cold X-ray absorption in these Seyfert 1s can be accounted for by neutral gas in clouds that also contribute to the broad line region (BLR) emission; and suggest that a geometrically thick neutral gas torus co-exists with the BLR and bridges the gap to the dusty torus.
[33]  [pdf] - 929700
Structure Formation in Gas-Rich Galactic Discs with Finite Thickness: From Discs to Rings
Comments: Accepted to MNRAS
Submitted: 2014-08-25, last modified: 2015-02-02
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.
[34]  [pdf] - 1215710
The Galactic Center cloud G2 and its gas streamer
Comments: 28 pages, paper as accepted by ApJ
Submitted: 2014-07-16, last modified: 2014-11-05
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.
[35]  [pdf] - 1222255
Time-resolved infrared emission from radiation-driven central obscuring structures in Active Galactic Nuclei
Comments: 16 pages, 12 figures, accepted for publication in MNRAS
Submitted: 2014-09-25
The central engines of Seyfert galaxies are thought to be enshrouded by geometrically thick gas and dust structures. In this article, we derive observable properties for a self-consistent model of such toroidal gas and dust distributions, where the geometrical thickness is achieved and maintained with the help of X-ray heating and radiation pressure due to the central engine. Spectral energy distributions (SEDs) and images are obtained with the help of dust continuum radiative transfer calculations with RADMC-3D. For the first time, we are able to present time-resolved SEDs and images for a physical model of the central obscurer. Temporal changes are mostly visible at shorter wavelengths, close to the combined peak of the dust opacity as well as the central source spectrum and are caused by variations in the column densities of the generated outflow. Due to the three-component morphology of the hydrodynamical models -- a thin disc with high density filaments, a surrounding fluffy component (the obscurer) and a low density outflow along the rotation axis -- we find dramatic differences depending on wavelength: whereas the mid-infrared images are dominated by the elongated appearance of the outflow cone, the long wavelength emission is mainly given by the cold and dense disc component. Overall, we find good agreement with observed characteristics, especially for those models, which show clear outflow cones in combination with a geometrically thick distribution of gas and dust, as well as a geometrically thin, but high column density disc in the equatorial plane.
[36]  [pdf] - 1209674
The central parsecs of active galactic nuclei: challenges to the torus
Comments: 23 pages, 10 figures, accepted for publication in MNRAS
Submitted: 2014-05-22
Type 2 AGN are by definition nuclei in which the broad-line region and continuum light are hidden at optical/UV wavelengths by dust. Via accurate registration of infrared (IR) Very Large Telescope adaptive optics images with optical \textit{Hubble Space Telescope} images we unambiguously identify the precise location of the nucleus of a sample of nearby, type 2 AGN. Dust extinction maps of the central few kpc of these galaxies are constructed from optical-IR colour images, which allow tracing the dust morphology at scales of few pc. In almost all cases, the IR nucleus is shifted by several tens of pc from the optical peak and its location is behind a dust filament, prompting to this being a major, if not the only, cause of the nucleus obscuration. These nuclear dust lanes have extinctions $A_V \geq 3-6$ mag, sufficient to at least hide the low-luminosity AGN class, and in some cases are observed to connect with kpc-scale dust structures, suggesting that these are the nuclear fueling channels. A precise location of the ionised gas H$\alpha$ and [\textsc{Si\,vii}] 2.48 $\mu$m coronal emission lines relative to those of the IR nucleus and dust is determined. The H$\alpha$ peak emission is often shifted from the nucleus location and its sometimes conical morphology appears not to be caused by a nuclear --torus-- collimation but to be strictly defined by the morphology of the nuclear dust lanes. Conversely, [\textsc{Si\,vii}] 2.48 $\mu$m emission, less subjected to dust extinction, reflects the truly, rather isotropic, distribution of the ionised gas. All together, the precise location of the dust, ionised gas and nucleus is found compelling enough to cast doubts on the universality of the pc-scale torus and supports its vanishing in low-luminosity AGN. Finally, we provide the most accurate position of the NGC 1068 nucleus, located at the South vertex of cloud B.
[37]  [pdf] - 759676
Observations of the gas cloud G2 in the Galactic Center
Comments: Proceedings of IAU Symposium #303, "The Galactic Center: Feeding and Feedback in a Normal Galactic Nucleus"; 30 Sep - 4 Oct 2013, Santa Fe / NM (USA)
Submitted: 2013-12-16
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.
[38]  [pdf] - 760174
The dusty torus in the Circinus galaxy: a dense disk and the torus funnel
Comments: 30 pages, 12 figures; A&A in press
Submitted: 2013-12-16
(Abridged) With infrared interferometry it is possible to resolve the nuclear dust distributions that are commonly associated with the dusty torus in active galactic nuclei (AGN). The Circinus galaxy hosts the closest Seyfert 2 nucleus and previous interferometric observations have shown that its nuclear dust emission is well resolved. To better constrain the dust morphology in this active nucleus, extensive new observations were carried out with MIDI at the Very Large Telescope Interferometer. The emission is distributed in two distinct components: a disk-like emission component with a size of ~ 0.2 $\times$ 1.1 pc and an extended component with a size of ~ 0.8 $\times$ 1.9 pc. The disk-like component is elongated along PA ~ 46{\deg} and oriented perpendicular to the ionisation cone and outflow. The extended component is elongated along PA ~ 107{\deg}, roughly perpendicular to the disk component and thus in polar direction. It is interpreted as emission from the inner funnel of an extended dust distribution and shows a strong increase in the extinction towards the south-east. We find no evidence of an increase in the temperature of the dust towards the centre. From this we infer that most of the near-infrared emission probably comes from parsec scales as well. We further argue that the disk component alone is not sufficient to provide the necessary obscuration and collimation of the ionising radiation and outflow. The material responsible for this must instead be located on scales of ~ 1 pc, surrounding the disk. The clear separation of the dust emission into a disk-like emitter and a polar elongated source will require an adaptation of our current understanding of the dust emission in AGN. The lack of any evidence of an increase in the dust temperature towards the centre poses a challenge for the picture of a centrally heated dust distribution.
[39]  [pdf] - 1201789
Hydrodynamical simulations of a compact source scenario for G2
Comments: 4 pages, 3 figures; Proceeding of the IAU 303: "The GC: Feeding and Feedback in a Normal Galactic Nucleus" / September 30 - October 4, 2013, Santa Fe, New Mexico (USA)
Submitted: 2013-11-21
The origin of the dense gas cloud G2 discovered in the Galactic Center (Gillessen et al. 2012) is still a debated puzzle. G2 might be a diffuse cloud or the result of an outflow from an invisible star embedded in it. We present here detailed simulations of the evolution of winds on G2's orbit. We find that the hydrodynamic interaction with the hot atmosphere present in the Galactic Center and the extreme gravitational field of the supermassive black hole must be taken in account when modeling such a source scenario. We find that the hydrodynamic interaction with the hot atmosphere present in the Galactic Center and the extreme gravitational field of the supermassive black hole must be taken in account when modeling such a source scenario. We also find that in this scenario most of the Br\gamma\ luminosity is expected to come from the highly filamentary densest shocked wind material. G2's observational properties can be used to constrain the properties of the outflow and our best model has a mass outflow rate of Mdot,w=8.8 x 10^{-8} Msun/yr and a wind velocity of vw = 50 km/s. These values are compatible with those of a young TTauri star wind, as already suggested by Scoville & Burkert (2013).
[40]  [pdf] - 1171695
Hydrodynamical simulations of a compact source scenario for the Galactic Center cloud G2
Comments: 25 pages, 14 figures. Accepted by ApJ
Submitted: 2013-05-30, last modified: 2013-09-19
The origin of the dense gas cloud G2 discovered in the Galactic Center (Gillessen et al. 2012) is still a debated puzzle. G2 might be a diffuse cloud or the result of an outflow from an invisible star embedded in it. We present hydrodynamical simulations of the evolution of different spherically symmetric winds of a stellar object embedded in G2. We find that the interaction with the ambient medium and with the extreme gravitational field of the supermassive black hole in the Galactic Center must be taken into account for such a source scenario. The thermal pressure of the hot and dense atmosphere confines the wind, while its ram pressure shapes it via stripping along the orbit, with the details depending on the wind parameters. Tidal forces squeeze the wind near pericenter, reducing it to a thin and elongated filament. We also find that in this scenario most of the Br\gamma\ luminosity is expected to come from the densest part of the wind, which has a highly filamentary structure with low filling factor. The observations can be best matched by a mass outflow rate of Mdot,w=8.8 x 10^{-8} Msun/yr and a wind velocity of vw = 50 km/s. These values are compatible with those of a young T Tauri star wind, as already suggested by (Scoville & Burkert 2013).
[41]  [pdf] - 925555
A diversity of dusty AGN tori: Data release for the VLTI/MIDI AGN Large Program and first results for 23 galaxies
Comments: A&A in press; 93 pages, 63 figures, 39 tables; data available only via CDS
Submitted: 2013-07-08
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.
[42]  [pdf] - 1171851
Pericenter passage of the gas cloud G2 in the Galactic Center
Comments: 10 figures, submitted to ApJ
Submitted: 2013-06-06
We have further followed the evolution of the orbital and physical properties of G2, the object currently falling toward the massive black hole in the Galactic Center on a near-radial orbit. New, very sensitive data were taken in April 2013 with NACO and SINFONI at the ESO VLT . The 'head' of G2 continues to be stretched ever further along the orbit in position-velocity space. A fraction of its emission appears to be already emerging on the blue-shifted side of the orbit, past pericenter approach. Ionized gas in the head is now stretched over more than 15,000 Schwarzschild radii RS around the pericenter of the orbit, at ~ 2000 RS ~ 20 light hours from the black hole. The pericenter passage of G2 will be a process stretching over a period of at least one year. The Brackett-{\gamma} luminosity of the head has been constant over the past 9 years, to within +- 25%, as have the line ratios Brackett-{\gamma} / Paschen-{\alpha} and Brackett-{\gamma} / Helium-I. We do not see any significant evidence for deviations of G2's dynamical evolution, due to hydrodynamical interactions with the hot gas around the black hole, from a ballistic orbit of an initially compact cloud with moderate velocity dispersion. The constant luminosity and the increasingly stretched appearance of the head of G2 in the position-velocity plane, without a central peak, is not consistent with several proposed models with continuous gas release from an initially bound zone around a faint star on the same orbit as G2.
[43]  [pdf] - 1166603
Numerical simulations of the possible origin of the two sub-parsec scale and counter-rotating stellar disks around SgrA*
Comments: 14 pages, 12 figures, submitted to ApJ
Submitted: 2013-05-13
We present a high resolution simulation of an idealized model to explain the origin of the two young, counter-rotating, sub-parsec scale stellar disks around the supermassive black hole SgrA* at the Center of the Milky Way. In our model, the collision of a single molecular cloud with a circum-nuclear gas disk (similar to the one observed presently) leads to multiple streams of gas flowing towards the black hole and creating accretion disks with angular momentum depending on the ratio of cloud and circum-nuclear disk material. The infalling gas creates two inclined, counter-rotating sub-parsec scale accretion disks around the supermassive black hole with the first disk forming roughly 1 Myr earlier, allowing it to fragment into stars and get dispersed before the second, counter-rotating disk forms. Fragmentation of the second disk would lead to the two inclined, counter-rotating stellar disks which are observed at the Galactic Center. A similar event might be happening again right now at the Milky Way Galactic Center. Our model predicts that the collision event generates spiral-like filaments of gas, feeding the Galactic Center prior to disk formation with a geometry and inflow pattern that is in agreement with the structure of the so called mini-spiral that has been detected in the Galactic Center.
[44]  [pdf] - 646411
Dynamics of gas and dust clouds in active galactic nuclei
Comments: 4 pages, 3 figures; language corrections, minor formatting changes
Submitted: 2013-02-18, last modified: 2013-03-12
We analyse the motion of single optically thick clouds in the potential of a central mass under the influence of an anisotropic radiation field ~|cos(\theta)|, a model applicable to the inner region of active galactic nuclei. Resulting orbits are analytically soluble for constant cloud column densities. All stable orbits are closed, although they have non-trivial shapes. Furthermore, there exists a stability criterion in the form of a critical inclination, which depends on the luminosity of the central source and the column density of the cloud.
[45]  [pdf] - 621981
Simulations of the origin and fate of the Galactic Center cloud G2
Comments: 6 pages, 3 figures, Proceedings of Science, "Nuclei of Seyfert galaxies and QSOs - Central engine & conditions of star formation", November 6-8, 2012, MPIfR, Bonn
Submitted: 2013-02-01
We investigate the origin and fate of the recently discovered gas cloud G2 close to the Galactic Center. Our hydrodynamical simulations focussing on the dynamical evolution of the cloud in combination with currently available observations favour two scenarios: a Compact Cloud which started around the year 1995 and an extended Spherical Shell of gas, with an apocenter distance within the disk(s) of young stars. The former is able to explain the detected signal of G2 in the position-velocity-diagram of the year 2008.5 and 2011.5 data. The latter can account for both, G2's signal as well as the fainter extended tail-like structure G2t seen at larger distances to the black hole and smaller velocities. From these first idealised simulations we expect a rise of the current activity of Sgr A* shortly after the closest approach and a constant feeding through a nozzle-like structure over a long period. The near future evolution of the cloud will be a sensitive probe of the conditions of the gas distribution in the milli-parsec environment of the massive black hole in the Galactic Center and will also give us invaluable information of the feeding of black holes and the activation of the central source.
[46]  [pdf] - 1151307
New observations of the gas cloud G2 in the Galactic Center
Comments: submitted to ApJ
Submitted: 2012-09-11
We present new observations of the recently discovered gas cloud G2 currently falling towards the massive black hole in the Galactic Center. The new data confirm that G2 is on a highly elliptical orbit with a predicted pericenter passage mid 2013. The updated orbit has an even larger eccentricity of 0.966, an epoch of pericenter two months later than estimated before, and a nominal minimum distance of 2200 Schwarzschild radii only. The velocity gradient of G2 has developed further to 600 km/s FWHM in summer 2012. We also detect the tail of similar total flux and on the same orbit as G2 along the trajectory at high significance. No hydrodynamic effects are detected yet, since the simple model of a tidally shearing gas cloud still describes the data very well. The flux of G2 has not changed by more than 10% between 2008 and 2012, disfavoring models where additional gas from a reservoir is released to the disrupting diffuse gas component.
[47]  [pdf] - 1124575
Magnetohydrodynamic stability of broad line region clouds
Comments: 16 pages, 11 figures, accepted by MNRAS, Figure 6 updated due to inconsistent linestyles, corrected mistake in Alfven speed formula, some very minor language corrections. The definitive version is available at www.blackwell-
Submitted: 2012-07-03, last modified: 2012-08-16
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.
[48]  [pdf] - 1117622
Simulations of the Origin and Fate of the Galactic Center Cloud G2
Comments: 16 pages, 16 figures, accepted by ApJ
Submitted: 2012-03-28, last modified: 2012-06-19
We investigate the origin and fate of the recently discovered gas cloud G2 close to the Galactic Center. Our hydrodynamical simulations focussing on the dynamical evolution of the cloud in combination with currently available observations favor two scenarios: a Compact Cloud which started around the year 1995 and a Spherical Shell of gas, with an apocenter distance within the disk(s) of young stars and a radius of a few times the size of the Compact Cloud. The former is able to explain the detected signal of G2 in the position-velocity diagram of the Br gamma emission of the year 2008.5 and 2011.5 data. The latter can account for both, G2's signal as well as the fainter extended tail-like structure G2t seen at larger distances from the black hole and smaller velocities. In contrast, gas stripped from a compact cloud by hydrodynamical interactions is not able to explain the location of the detected G2t emission in the observed position-velocity diagrams. This favors the Spherical Shell Scenario and might be a severe problem for the Compact Cloud as well as the so-called Compact Source Scenario. From these first idealized simulations we expect a roughly constant feeding of the supermassive black hole through a nozzle-like structure over a long period, starting shortly after the closest approach in 2013.51 for the Compact Cloud. If the matter accretes in the hot accretion mode, we do not expect a significant boost of the current activity of Sgr A* for the Compact Cloud model, but a boost of the average infrared and X-ray luminosity by roughly a factor of 80 for the Spherical Shell scenario with order of magnitude variations on a timescale of a few months. The near-future evolution of the cloud will be a sensitive probe of the conditions of the gas distribution in the milli-parsec environment of the massive black hole in the Galactic Center.
[49]  [pdf] - 1117932
The complexity of parsec-scaled dusty tori in AGN
Comments: 8 pages, 3 figures, to appear in the proceedings of the conference "The central kiloparsec in Galactic Nuclei: Astronomy at High Angular Resolution 2011", open access Journal of Physics: Conference Series (JPCS), published by IOP Publishing
Submitted: 2012-04-12
Warm gas and dust surround the innermost regions of active galactic nuclei (AGN). They provide the material for accretion onto the super-massive black hole and they are held responsible for the orientation-dependent obscuration of the central engine. The AGN-heated dust distributions turn out to be very compact with sizes on scales of about a parsec in the mid-infrared. Only infrared interferometry currently provides the necessary angular resolution to directly study the physical properties of this dust. Size estimates for the dust distributions derived from interferometric observations can be used to construct a size--luminosity relation for the dust distributions. The large scatter about this relation suggests significant differences between the dust tori in the individual galaxies, even for nuclei of the same class of objects and with similar luminosities. This questions the simple picture of the same dusty doughnut in all AGN. The Circinus galaxy is the closest Seyfert 2 galaxy. Because its mid-infrared emission is well resolved interferometrically, it is a prime target for detailed studies of its nuclear dust distribution. An extensive new interferometric data set was obtained for this galaxy. It shows that the dust emission comes from a very dense, disk-like structure which is surrounded by a geometrically thick, similarly warm dust distribution as well as significant amounts of warm dust within the ionisation cone.
[50]  [pdf] - 1092775
Physics of the Galactic Center Cloud G2, on its Way towards the Super-Massive Black Hole
Comments: 22 pages, 13 figures, submitted to ApJ
Submitted: 2012-01-06
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.
[51]  [pdf] - 1076243
The Role of Secular Evolution in the Black Hole Growth of Narrow-Line Seyfert 1 Galaxies
Comments: 18 pages, 10 figures, 2 tables, accepted by MNRAS
Submitted: 2011-04-26, last modified: 2011-07-15
Narrow-Line Seyfert 1 (NLS1) galaxies show extreme properties with respect to the other Seyfert galaxies. Indeed, they are thought to be accreting at Eddington rates and to possess low mass black holes. Therefore, they may represent a key class of objects for understanding the co-evolution of black holes and their host galaxies. We propose that NLS1s represent a class of AGN in which the black hole growth is, and has always been, dominated by secular evolution. Firstly, by looking at the NLS1 host galaxy properties in the literature, we show that the evolution of NLS1s is presently driven by secular processes, much more so than for Broad-Line Seyfert 1s (BLS1s). Secondly, we study the bulges of NLS1 and BLS1 galaxies. Our results demonstrate that NLS1 host bulges are pseudo-bulges and are statistically different from BLS1 bulges. This difference points to the particular importance of secular processes in the past evolution of their hosts. We build on this result to understand the implications on their evolution and the duration of their duty cycle. We show that NLS1s are not necessarily in a special phase of black hole growth and that several Gyr are required for their black hole masses to become similar to BLS1s. Finally, in the light of our results, we discuss the location of NLS1 galaxies on the MBH-\sigma plane and speculate about the connection between the NLS1 galaxy properties and their black hole spin.
[52]  [pdf] - 1076851
On the size-luminosity relation of AGN dust tori in the mid-infrared
Comments: 10 pages, 4 figures, accepted for publication in A&A
Submitted: 2011-05-24
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.
[53]  [pdf] - 1052892
Radiation feedback on dusty clouds during Seyfert activity
Comments: 13 pages, 18 figures, accepted by MNRAS
Submitted: 2011-03-22
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.
[54]  [pdf] - 351010
Simulations of Direct Collisions of Gas Clouds with the Central Black Hole
Comments: 21 pages, 18 figures, Accepted for publication in MNRAS
Submitted: 2009-08-07, last modified: 2010-10-27
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.
[55]  [pdf] - 1033442
Stability of Cloud Orbits in the Broad Line Region of Active Galactic Nuclei
Comments: 7 pages, 4 fig.s, accepted by MNRAS, some text was missing - corrected
Submitted: 2010-07-01, last modified: 2010-10-11
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.
[56]  [pdf] - 197345
Radiatively enhanced elasticity and turbulence in clumpy tori of Active Galactic Nuclei
Comments: 10 pages, 5 figures, accepted by MNRAS, this paper has been withdrawn by the authors due to an essential flaw in the assumptions
Submitted: 2009-11-25, last modified: 2010-07-01
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.
[57]  [pdf] - 1026250
Modelling Shear Flows with SPH and Grid Based Methods
Comments: 16 pages, 17 figures
Submitted: 2010-04-12
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.
[58]  [pdf] - 1018899
Gas dynamics of the central few parsec region of NGC 1068 fuelled by the evolving nuclear star cluster
Comments: 13 pages, 11 figures, accepted for publication in MNRAS, updated author list and references
Submitted: 2009-12-23, last modified: 2010-01-19
High resolution observations with the NIR adaptive optics integral field spectrograph SINFONI at the VLT proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of three-dimensional high resolution hydrodynamical simulations with the Pluto code, we follow the evolution of such clusters, focusing on stellar mass loss. This leads to clumpy or filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. In order to capture the relevant physics in the inner region, we treat this disc separately by viscously evolving the radial surface density distribution. This enables us to link the tens of parsec scale region (accessible via SINFONI observations) to the (sub-)parsec scale region (observable with the MIDI instrument and via water maser emission). In this work, we concentrate on the effects of a parametrised turbulent viscosity to generate angular momentum and mass transfer in the disc and additionally take star formation into account. Input parameters are constrained by observations of the nearby Seyfert 2 galaxy NGC 1068. At the current age of its nuclear starburst of 250 Myr, our simulations yield disc sizes of the order of 0.8 to 0.9 pc, gas masses of 1.0e6 solar masses and mass transfer rates of 0.025 solar masses per year through the inner rim of the disc. This shows that our large scale torus model is able to approximately account for the disc size as inferred from interferometric observations in the mid-infrared and compares well to the extent and mass of a rotating disc structure as inferred from water maser observations. Several other observational constraints are discussed as well.
[59]  [pdf] - 1018556
Obscuring and feeding supermassive black holes with evolving nuclear star clusters
Comments: 6 pages, 3 figures, to appear in the proceedings of the IAU General Assembly 2009, Rio de Janeiro, S267 Co-evolution of Central Black Holes and Galaxies
Submitted: 2009-11-30
Recently, high resolution observations with the help of the near-infrared adaptive optics integral field spectrograph SINFONI at the VLT proved the existence of massive and young nuclear star clusters in the centres of a sample of Seyfert galaxies. With the help of high resolution hydrodynamical simulations with the PLUTO-code, we follow the evolution of such clusters, especially focusing on mass and energy feedback from young stars. This leads to a filamentary inflow of gas on large scales (tens of parsec), whereas a turbulent and very dense disc builds up on the parsec scale. Here, we concentrate on the long-term evolution of the nuclear disc in NGC 1068 with the help of an effective viscous disc model, using the mass input from the large scale simulations and accounting for star formation in the disc. This two-stage modelling enables us to connect the tens of parsec scale region (observable with SINFONI) with the parsec scale environment (MIDI observations). At the current age of the nuclear star cluster, our simulations predict disc sizes of the order of 0.8 to 0.9 pc, gas masses of 1.0e6 Msun and mass transfer rates through the inner boundary of 0.025 Msun/yr in good agreement with values derived from observations.
[60]  [pdf] - 1018470
[The Impact of Nuclear Star Formation on Gas Inflow to AGN
Comments: to appear in: Co-Evolution of Central Black Holes and Galaxies; 7 pages
Submitted: 2009-11-23
Our adaptive optics observations of nearby AGN at spatial resolutions as small as 0.085arcsec show strong evidence for recent, but no longer active, nuclear star formation. We begin by describing observations that highlight two contrasting methods by which gas can flow into the central tens of parsecs. Gas accumulation in this region will inevitably lead to a starburst, and we discuss the evidence for such events. We then turn to the impact of stellar evolution on the further inflow of gas by combining a phenomenological approach with analytical modelling and hydrodynamic simulations. These complementary perspectives paint a picture in which all the processes are ultimately regulated by the mass accretion rate into the central hundred parsecs, and the ensuing starburst that occurs there. The resulting supernovae delay accretion by generating a starburst wind, which leaves behind a clumpy interstellar medium. This provides an ideal environment for slower stellar outflows to accrete inwards and form a dense turbulent disk on scales of a few parsecs. Such a scenario may resolve the discrepancy between the larger scale structure seen with adaptive optics and the small scale structure seen with VLTI.
[61]  [pdf] - 18346
The effect of stellar feedback on the formation and evolution of gas and dust tori in AGN
Comments: 16 pages, 11 figures, accepted by MNRAS, high resolution version can be downloaded from:
Submitted: 2008-11-10
Recently, the existence of geometrically thick dust structures in Active Galactic Nuclei (AGN) has been directly proven with the help of mid-infrared interferometry. The observations are consistent with a two-component model made up of a geometrically thin and warm central disk, surrounded by a colder, fluffy torus component. In an exploratory study, we investigate one possible physical mechanism, which could produce such a structure, namely the effect of stellar feedback from a young nuclear star cluster on the interstellar medium in centres of AGN. The model is realised with the help of the hydrodynamics code TRAMP. We follow the evolution of the interstellar medium by taking discrete mass loss and energy ejection due to stellar processes, as well as optically thin radiative cooling into account. In a post-processing step, we calculate observable quantities (spectral energy distributions and images) with the help of the radiative transfer code MC3D. The interplay between injection of mass, supernova explosions and radiative cooling leads to a two-component structure made up of a cold geometrically thin, but optically thick and very turbulent disk residing in the vicinity of the angular momentum barrier, surrounded by a filamentary structure. The latter consists of cold long radial filaments flowing towards the disk and a hot tenuous medium in between, which shows both inwards and outwards directed motions. This modelling is able to reproduce the range of observed neutral hydrogen column densities of a sample of Seyfert galaxies as well as the relation between them and the strength of the silicate 10 micron spectral feature. Despite being quite crude, our mean Seyfert galaxy model is even able to describe the SEDs of two intermediate type Seyfert galaxies observed with the Spitzer Space Telescope.
[62]  [pdf] - 15163
Turbulent AGN tori
Comments: 4 pages, 2 figures; Talk given at "The Central Kiloparsec: Active Galactic Nuclei and Their Hosts", Ierapetra, Greece, 4-6 June, 2008; to appear in Volume 79 of the Memorie della Societa Astronomica Italiana
Submitted: 2008-08-05
Recently, the MID-infrared Interferometric instrument (MIDI) at the VLTI has shown that dust tori in the two nearby Seyfert galaxies NGC 1068 and the Circinus galaxy are geometrically thick and can be well described by a thin, warm central disk, surrounded by a colder and fluffy torus component. By carrying out hydrodynamical simulations with the help of the TRAMP code (Klahr et al. 1999), we follow the evolution of a young nuclear star cluster in terms of discrete mass-loss and energy injection from stellar processes. This naturally leads to a filamentary large scale torus component, where cold gas is able to flow radially inwards. The filaments open out into a dense and very turbulent disk structure. In a post-processing step, we calculate observable quantities like spectral energy distributions or images with the help of the 3D radiative transfer code MC3D (Wolf 2003). Good agreement is found in comparisons with data due to the existence of almost dust-free lines of sight through the large scale component and the large column densities caused by the dense disk.
[63]  [pdf] - 10213
Three-dimensional radiative transfer models of clumpy tori in Seyfert galaxies
Comments: 15 pages, 23 figures, accepted by A&A
Submitted: 2008-02-19
Tori of Active Galactic Nuclei are made up of a mixture of hot and cold gas, as well as dust. In order to protect the dust grains from destruction by the hot gas as well as by the energetic radiation of the accretion disk, the dust is often assumed to be distributed in clouds. In our new 3D model of AGN dust tori, the torus is modelled as a wedge-shaped disk in which dusty clouds are randomly distributed, by taking the dust density distribution of the corresponding continuous model into account. We especially concentrate on the differences between clumpy and continuous models in terms of the temperature distributions, the surface brightness distributions and interferometric visibilities, as well as spectral energy distributions. To this end, we employ radiative transfer calculations with the help of the 3D Monte Carlo code MC3D. In a second step, interferometric visibilities are calculated from the simulated surface brightness distributions, which can be directly compared to observations with the MIDI instrument. The radial temperature distributions of clumpy models possess significantly enhanced scatter compared to the continuous cases. Even at large distances, clouds can be heated directly by the central accretion disk. The existence of the silicate 10 micron-feature in absorption or in emission depends sensitively on the distribution, the size and optical depth of clouds in the innermost part of the torus, due to shadowing effects of clouds there. This explains failure and success of previous modelling efforts of clumpy tori. After adapting the parameters of our clumpy standard model to the circumstances of the Seyfert 2 Circinus galaxy, it can qualitatively explain recent mid-infrared interferometric observations performed with MIDI, as well as high resolution spectral data.
[64]  [pdf] - 4494
Resolving the complex structure of the dust torus in the active nucleus of the Circinus galaxy
Comments: 20 pages, 16 figures, accepted for publication by A&A
Submitted: 2007-09-03
To test the dust torus model for active galactic nuclei directly, we study the extent and morphology of the nuclear dust distribution in the Circinus galaxy using high resolution interferometric observations in the mid-infrared with the MIDI instrument at the Very Large Telescope Interferometer. We find that the dust distribution in the nucleus of Circinus can be explained by two components, a dense and warm disk-like component of 0.4 pc size and a slightly cooler, geometrically thick torus component with a size of 2.0 pc. The disk component is oriented perpendicular to the ionisation cone and outflow and seems to show the silicate feature at 10 micron in emission. It coincides with a nuclear maser disk in orientation and size. From the energy needed to heat the dust, we infer a luminosity of the accretion disk corresponding to 20% of the Eddington luminosity of the nuclear black hole. We find that the interferometric data are inconsistent with a simple, smooth and axisymmetric dust emission. The irregular behaviour of the visibilities and the shallow decrease of the dust temperature with radius provide strong evidence for a clumpy or filamentary dust structure. We see no evidence for dust reprocessing, as the silicate absorption profile is consistent with that of standard galactic dust. We argue that the collimation of the ionising radiation must originate in the geometrically thick torus component. Our findings confirm the presence of a geometrically thick, torus-like dust distribution in the nucleus of Circinus, as required in unified schemes of Seyfert galaxies. Several aspects of our data require that this torus is irregular, or "clumpy".
[65]  [pdf] - 72202
Towards a physical model of dust tori in Active Galactic Nuclei - Radiative transfer calculations for a hydrostatic torus model
Comments: 23 pages, 22 figures, accepted by A&A, figure 2 corrected
Submitted: 2005-04-05, last modified: 2005-05-20
We explore physically self-consistent models of dusty molecular tori in Active Galactic Nuclei (AGN) with the goal of interpreting VLTI observations and fitting high resolution mid-IR spectral energy distributions (SEDs). The input dust distribution is analytically calculated by assuming hydrostatic equilibrium between pressure forces - due to the turbulent motion of the gas clouds - and gravitational and centrifugal forces as a result of the contribution of the nuclear stellar distribution and the central black hole. For a fully three-dimensional treatment of the radiative transfer problem through the tori we employ the Monte Carlo code MC3D. We find that in homogeneous dust distributions the observed mid-infrared emission is dominated by the inner funnel of the torus, even when observing along the equatorial plane. Therefore, the stratification of the distribution of dust grains - both in terms of size and composition - cannot be neglected. In the current study we only include the effect of different sublimation radii which significantly alters the SED in comparison to models that assume an average dust grain property with a common sublimation radius, and suppresses the silicate emission feature at 9.7 micron. In this way we are able to fit the mean SED of both type I and type II AGN very well. Our fit of special objects for which high angular resolution observations (less than 0.3 arcseconds) are available indicates that the hottest dust in NGC 1068 reaches the sublimation temperature while the maximum dust temperature in the low-luminosity AGN Circinus falls short of 1000 K.