Normalized to: Schartmann, M.
[1]
oai:arXiv.org:2003.06153 [pdf] - 2063856
Ionized outflows in local luminous AGN: what are the real densities and
outflow rates?
Davies, R.;
Baron, D.;
Shimizu, T.;
Netzer, H.;
Burtscher, L.;
de Zeeuw, P. T.;
Genzel, R.;
Hicks, E. K. S.;
Koss, M.;
Lin, M. -Y.;
Lutz, D.;
Maciejewski, W.;
Müller-Sánchez, F.;
de Xivry, G. Orban;
Ricci, C.;
Riffel, R.;
Riffel, R. A.;
Rosario, D.;
Schartmann, M.;
Schnorr-Müller, A.;
Shangguan, J.;
Sternberg, A.;
Sturm, E.;
Storchi-Bergmann, T.;
Tacconi, L.;
Veilleux, S.
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]
oai:arXiv.org:1906.04181 [pdf] - 2034364
3D simulations of clump formation in stellar wind collisions
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]
oai:arXiv.org:1912.07734 [pdf] - 2050306
LLAMA: The $M_{BH}$ - $\sigma_{\star}$ Relation of the most luminous
local AGNs
Caglar, Turgay;
Burtscher, Leonard;
Brand, Bernhard;
Brinchmann, Jarle;
Davies, Richard I.;
Hicks, Erin K. S.;
Koss, Michael;
Lin, Ming-Yi;
Maciejewski, Witold;
Müller-Sánchez, Francisco;
Riffel, Rogemar A.;
Riffel, Rogério;
Rosario, David J.;
Schartmann, Marc;
Schnorr-Müller, Allan;
Shimizu, T. Taro;
Storchi-Bergmann, Thaisa;
Veilleux, Sylvain;
de Xivry, Gilles O.;
Bennert, Vardha N.
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]
oai:arXiv.org:1910.06976 [pdf] - 2026019
Stellar winds pump the heart of the Milky Way
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]
oai:arXiv.org:1912.01361 [pdf] - 2042294
An image of the dust sublimation region in the nucleus of NGC 1068
GRAVITY Collaboration;
Pfuhl, O.;
Davies, R.;
Dexter, J.;
Netzer, H.;
Hoenig, S.;
Lutz, D.;
Schartmann, M.;
Sturm, E.;
Amorim, A.;
Brandner, W.;
Clenet, Y.;
de Zeeuw, P. T.;
Eckart, A.;
Eisenhauer, F.;
Schreiber, N. M. Foerster;
Gao, F.;
Garcia, P. J. V.;
Genzel, R.;
Gillessen, S.;
Gratadour, D.;
Kishimoto, M.;
Lacour, S.;
Millour, F.;
Ott, T.;
Paumard, T.;
Perraut, K.;
Perrin, G.;
Peterson, B. M.;
Petrucci, P. O.;
Prieto, M. A.;
Rouan, D.;
Shangguan, J.;
Shimizu, T.;
Sternberg, A.;
Straub, O.;
Straubmeier, C.;
Tacconi, L. J.;
Tristram, K. R. W.;
Vermot, P.;
Waisberg, I.;
Widmann, F.;
Woillez, J.
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]
oai:arXiv.org:1907.03801 [pdf] - 1983882
The multiphase gas structure and kinematics in the circumnuclear region
of NGC 5728
Shimizu, T. Taro;
Davies, R. I.;
Lutz, D.;
Burtscher, L.;
Lin, M.;
Baron, D.;
Davies, R. L.;
Genzel, R.;
Hicks, E. K. S.;
Koss, M.;
Maciejewski, W.;
ller-Sanchez, F. Mu;
de Xivry, G. Orban;
Price, S. H.;
Ricci, C.;
Riffel, R.;
Riffel, R. A.;
Rosario, D.;
Schartmann, M.;
Schnorr-Muller, A.;
Sternberg, A.;
Sturm, E.;
Storchi-Bergmann, T.;
Tacconi, L.;
Veilleux, S.
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]
oai:arXiv.org:1910.00593 [pdf] - 2065272
The resolved size and structure of hot dust in the immediate vicinity of
AGN
GRAVITY Collaboration;
Dexter, J.;
Shangguan, J.;
Hönig, S.;
Kishimoto, M.;
Lutz, D.;
Netzer, H.;
Davies, R.;
Sturm, E.;
Pfuhl, O.;
Amorim, A.;
Bauböck, M.;
Brandner, W.;
Clénet, Y.;
de Zeeuw, P. T.;
Eckart, A.;
Eisenhauer, F.;
Schreiber, N. M. Förster;
Gao, F.;
Garcia, P. J. V.;
Genzel, R.;
Gillessen, S.;
Gratadour, D.;
Jiménez-Rosales, A.;
Lacour, S.;
Millour, F.;
Ott, T.;
Paumard, T.;
Perraut, K.;
Perrin, G.;
Peterson, B. M.;
Petrucci, P. O.;
Prieto, M. A.;
Rouan, D.;
Schartmann, M.;
Shimizu, T.;
Sternberg, A.;
Straub, O.;
Straubmeier, C.;
Tacconi, L. J.;
Tristram, K.;
Vermot, P.;
Waisberg, I.;
Widmann, F.;
Woillez, J.
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]
oai:arXiv.org:1906.06889 [pdf] - 1938474
A nuclear molecular outflow in the Seyfert galaxy NGC3227
Alonso-Herrero, A.;
Garcia-Burillo, S.;
Pereira-Santaella, M.;
Davies, R. I.;
Combes, F.;
Vestergaard, M.;
Raimundo, S. I.;
Bunker, A.;
Diaz-Santos, T.;
Gandhi, P.;
Garcia-Bernete, I.;
Hicks, E. K. S.;
Hönig, S. F.;
Hunt, L. K.;
Imanishi, M.;
Izumi, T.;
Levenson, N. A.;
Maciejewski, W.;
Packham, C.;
Almeida, C. Ramos;
Ricci, C.;
Rigopoulou, D.;
Roche, P. F.;
Rosario, D.;
Schartmann, M.;
Usero, A.;
Ward, M. J.
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]
oai:arXiv.org:1902.09370 [pdf] - 1846941
From kpcs to the central parsec of NGC 1097: feeding star formation and
a black hole at the same time
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]
oai:arXiv.org:1812.01416 [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
Gillessen, S.;
Plewa, P. M.;
Widmann, F.;
von Fellenberg, S.;
Schartmann, M.;
Habibi, M.;
Rosales, A. Jimenez;
Bauböck, M.;
Dexter, J.;
Gao, F.;
Waisberg, I.;
Eisenhauer, F.;
Pfuhl, O.;
Ott, T.;
Burkert, A.;
de Zeeuw, P. T.;
Genzel, R.
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]
oai:arXiv.org:1804.02860 [pdf] - 1736339
Simulating the pericentre passage of the Galactic centre star S2
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]
oai:arXiv.org:1805.05341 [pdf] - 1694163
The Galactic Centre source G2 was unlikely born in any of the known
massive binaries
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]
oai:arXiv.org:1803.06182 [pdf] - 1732608
Thick turbulent gas disks with magnetocentrifugal winds in active
galactic nuclei - Model infrared emission and optical polarization
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]
oai:arXiv.org:1710.04224 [pdf] - 1604978
LLAMA: Normal star formation efficiencies of molecular gas in the
centres of luminous Seyfert galaxies
Rosario, D. J.;
Burtscher, L.;
Davies, R. I.;
Koss, M.;
Ricci, C.;
Lutz, D.;
Riffel, R.;
Alexander, D. M.;
Genzel, R.;
Hicks, E. H.;
Lin, M. -Y.;
Maciejewski, W.;
Sanchez, F. Mueller-;
de Xivry, G. Orban;
Riffel, R. A.;
Schartmann, M.;
Schawinski, K.;
Schnorr-Mueller, A.;
Saintonge, A.;
Shimizu, T. T.;
Sternberg, A.;
Storchi-Bergmann, T.;
Sturm, E.;
Tacconi, L.;
Treister, E.;
Veilleux, S.
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]
oai:arXiv.org:1710.04098 [pdf] - 1589584
LLAMA: Nuclear stellar properties of Swift BAT AGN and matched inactive
galaxies
Lin, Ming-Yi;
Davies, R. I.;
Hicks, E. K. S.;
Burtscher, L.;
Contursi, A.;
Genzel, R.;
Koss, M.;
Lutz, D.;
Maciejewski, W.;
Müller-Sánchez, F.;
de Xivry, G. Orban;
Ricci, C.;
Riffel, R.;
Riffel, R. A.;
Rosario, D.;
Schartmann, M.;
Schnorr-Müller, A.;
Shimizu, T.;
Sternberg, A.;
Sturm, E.;
Storchi-Bergmann, T.;
Tacconi, L.;
Veilleux, S.
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]
oai:arXiv.org:1709.03367 [pdf] - 1588193
The life cycle of starbursting circumnuclear gas discs
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]
oai:arXiv.org:1708.06528 [pdf] - 1587347
Polarized radiative transfer modeling of warped and clumpy dusty tori
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]
oai:arXiv.org:1706.08547 [pdf] - 1697890
3D AMR hydrosimulations of a compact source scenario for the Galactic
Centre cloud G2
Ballone, A.;
Schartmann, M.;
Burkert, A.;
Gillessen, S.;
Plewa, P. M.;
Genzel, R.;
Pfuhl, O.;
Eisenhauer, F.;
Habibi, M.;
Ott, T.;
George, E. M.
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]
oai:arXiv.org:1703.03894 [pdf] - 1557938
IC 630: Piercing the Veil of the Nuclear Gas
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]
oai:arXiv.org:1610.09890 [pdf] - 1532342
The Role of Host Galaxy for the Environmental Dependence of Active
Nuclei in Local Galaxies
Davies, R. I.;
Hicks, E. K. S.;
Erwin, P.;
Burtscher, L.;
Contursi, A.;
Genzel, R.;
Janssen, A.;
Koss, M.;
Lin, M. -Y.;
Lutz, D.;
Maciejewski, W.;
Mueller-Sanchez, F.;
de Xivry, G. Orban;
Ricci, C.;
Riffel, R.;
Riffel, R. A.;
Rosario, D.;
Schartmann, M.;
Schnorr-Mueller, A.;
Shimizu, T.;
Sternberg, A.;
Sturm, E.;
Storchi-Bergmann, T.;
Tacconi, L.;
Veilleux, S.
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]
oai:arXiv.org:1610.07247 [pdf] - 1532200
Radiative transfer modelling of parsec-scale dusty warped discs
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]
oai:arXiv.org:1609.07318 [pdf] - 1531712
3D AMR simulations of G2 as an outflow
Ballone, A.;
Schartmann, M.;
Burkert, A.;
Gillessen, S.;
Plewa, P. M.;
Pfuhl, O.;
Genzel, R.;
Eisenhauer, F.;
Ott, T.;
George, E. M.;
Habibi, M.
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]
oai:arXiv.org:1609.07308 [pdf] - 1531711
3D AMR simulations of the evolution of the diffuse gas cloud G2 in the
Galactic Centre
Schartmann, M.;
Ballone, A.;
Burkert, A.;
Gillessen, S.;
Genzel, R.;
Pfuhl, O.;
Eisenhauer, F.;
Plewa, P. M.;
Ott, T.;
George, E. M.;
Habibi, M.
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]
oai:arXiv.org:1608.06995 [pdf] - 1483503
Multi-phase Nature of a Radiation-Driven Fountain with Nuclear Starburst
in a Low-mass Active Galactic Nucleus
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]
oai:arXiv.org:1607.07308 [pdf] - 1475452
Constraints on the Broad Line Region Properties and Extinction in Local
Seyferts
Schnorr-Müller, Allan;
Davies, R. I.;
Korista, K. T.;
Burtscher, L.;
Rosario, D.;
Storchi-Bergmann, T.;
Contursi, A.;
Genzel, R.;
Graciá-Carpio, J.;
Hicks, E. K. S.;
Janssen, A.;
Koss, M.;
Lin, M. -Y.;
Lutz, D.;
Maciejewski, W.;
Müller-Sánchez, F.;
de Xivry, G. Orban;
Riffel, R.;
Riffel, R. A.;
Schartmann, M.;
Sternberg, A.;
Sturm, E.;
Tacconi, L.;
Veilleux, S.;
Ulrich, O. A.
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]
oai:arXiv.org:1602.07305 [pdf] - 1378969
The G2+G2t complex as a fast and massive outflow?
Ballone, A.;
Schartmann, M.;
Burkert, A.;
Gillessen, S.;
Plewa, P. M.;
Genzel, R.;
Pfuhl, O.;
Eisenhauer, F.;
Ott, T.;
George, E. M.;
Habibi, M.
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]
oai:arXiv.org:1602.05592 [pdf] - 1422232
Mid-infrared interferometry of 23 AGN tori: On the significance of
polar-elongated emission
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]
oai:arXiv.org:1512.03430 [pdf] - 1374211
Clusters of Small Clumps Can Explain The Peculiar Properties of Giant
Clumps in High-Redshift Galaxies
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]
oai:arXiv.org:1511.05151 [pdf] - 1331142
Stellar feedback efficiencies: supernovae versus stellar winds
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]
oai:arXiv.org:1507.07012 [pdf] - 1327446
Clump formation through colliding stellar winds in the Galactic Centre
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]
oai:arXiv.org:1508.07060 [pdf] - 1292296
3D adaptive mesh refinement simulations of the gas cloud G2 born within
the disks of young stars in the Galactic Center
Schartmann, M.;
Ballone, A.;
Burkert, A.;
Gillessen, S.;
Genzel, R.;
Pfuhl, O.;
Eisenhauer, F.;
Plewa, P. M.;
Ott, T.;
George, E. M.;
Habibi, M.
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]
oai:arXiv.org:1505.00536 [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
Davies, R.;
Burtscher, L.;
Rosario, D.;
Storchi-Bergmann, T.;
Contursi, A.;
Genzel, R.;
Graciá-Carpio, J.;
Hicks, E.;
Janssen, A.;
Koss, M.;
Lin, M. -Y.;
Lutz, D.;
Maciejewski, W.;
Müller-Sánchez, F.;
de Xivry, G. Orban;
Ricci, C.;
Riffel, R.;
Riffel, R. A.;
Schartmann, M.;
Schnorr-Müller, A.;
Sternberg, A.;
Sturm, E.;
Tacconi, L.;
Veilleux, S.
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]
oai:arXiv.org:1408.5902 [pdf] - 929700
Structure Formation in Gas-Rich Galactic Discs with Finite Thickness:
From Discs to Rings
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]
oai:arXiv.org:1407.4354 [pdf] - 1215710
The Galactic Center cloud G2 and its gas streamer
Pfuhl, Oliver;
Gillessen, Stefan;
Eisenhauer, Frank;
Genzel, Reinhard;
Plewa, Philipp M.;
Ott, Thomas;
Ballone, Alessandro;
Schartmann, Marc;
Burkert, Andreas;
Fritz, Tobias K.;
Sari, Re'em;
Steinberg, Elad;
Madigan, Ann-Marie
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]
oai:arXiv.org:1409.7404 [pdf] - 1222255
Time-resolved infrared emission from radiation-driven central obscuring
structures in Active Galactic Nuclei
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]
oai:arXiv.org:1405.5653 [pdf] - 1209674
The central parsecs of active galactic nuclei: challenges to the torus
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]
oai:arXiv.org:1312.4386 [pdf] - 759676
Observations of the gas cloud G2 in the Galactic Center
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]
oai:arXiv.org:1312.4534 [pdf] - 760174
The dusty torus in the Circinus galaxy: a dense disk and the torus
funnel
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]
oai:arXiv.org:1311.5462 [pdf] - 1201789
Hydrodynamical simulations of a compact source scenario for G2
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]
oai:arXiv.org:1305.7238 [pdf] - 1171695
Hydrodynamical simulations of a compact source scenario for the Galactic
Center cloud G2
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]
oai:arXiv.org:1307.2068 [pdf] - 925555
A diversity of dusty AGN tori: Data release for the VLTI/MIDI AGN Large
Program and first results for 23 galaxies
Burtscher, Leonard;
Meisenheimer, Klaus;
Tristram, Konrad R. W.;
Jaffe, Walter;
Hönig, Sebastian F.;
Davies, Richard I.;
Kishimoto, Makoto;
Pott, Jörg-Uwe;
Röttgering, Huub;
Schartmann, Marc;
Weigelt, Gerd;
Wolf, Sebastian
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]
oai:arXiv.org:1306.1374 [pdf] - 1171851
Pericenter passage of the gas cloud G2 in the Galactic Center
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]
oai:arXiv.org:1305.2953 [pdf] - 1166603
Numerical simulations of the possible origin of the two sub-parsec scale
and counter-rotating stellar disks around SgrA*
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]
oai:arXiv.org:1302.4287 [pdf] - 646411
Dynamics of gas and dust clouds in active galactic nuclei
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]
oai:arXiv.org:1302.0124 [pdf] - 621981
Simulations of the origin and fate of the Galactic Center cloud G2
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]
oai:arXiv.org:1209.2272 [pdf] - 1151307
New observations of the gas cloud G2 in the Galactic Center
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]
oai:arXiv.org:1207.0785 [pdf] - 1124575
Magnetohydrodynamic stability of broad line region clouds
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]
oai:arXiv.org:1203.6356 [pdf] - 1117622
Simulations of the Origin and Fate of the Galactic Center Cloud G2
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]
oai:arXiv.org:1204.2823 [pdf] - 1117932
The complexity of parsec-scaled dusty tori in AGN
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]
oai:arXiv.org:1201.1414 [pdf] - 1092775
Physics of the Galactic Center Cloud G2, on its Way towards the
Super-Massive Black Hole
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]
oai:arXiv.org:1104.5023 [pdf] - 1076243
The Role of Secular Evolution in the Black Hole Growth of Narrow-Line
Seyfert 1 Galaxies
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]
oai:arXiv.org:1105.4875 [pdf] - 1076851
On the size-luminosity relation of AGN dust tori in the mid-infrared
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]
oai:arXiv.org:1103.4368 [pdf] - 1052892
Radiation feedback on dusty clouds during Seyfert activity
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]
oai:arXiv.org:0908.1100 [pdf] - 351010
Simulations of Direct Collisions of Gas Clouds with the Central Black
Hole
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]
oai:arXiv.org:1007.0112 [pdf] - 1033442
Stability of Cloud Orbits in the Broad Line Region of Active Galactic
Nuclei
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]
oai:arXiv.org:0911.4971 [pdf] - 197345
Radiatively enhanced elasticity and turbulence in clumpy tori of Active
Galactic Nuclei
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]
oai:arXiv.org:1004.1957 [pdf] - 1026250
Modelling Shear Flows with SPH and Grid Based Methods
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]
oai:arXiv.org:0912.4677 [pdf] - 1018899
Gas dynamics of the central few parsec region of NGC 1068 fuelled by the
evolving nuclear star cluster
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]
oai:arXiv.org:0911.5641 [pdf] - 1018556
Obscuring and feeding supermassive black holes with evolving nuclear
star clusters
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]
oai:arXiv.org:0911.4355 [pdf] - 1018470
[The Impact of Nuclear Star Formation on Gas Inflow to AGN
Davies, R. I.;
Hicks, E.;
Schartmann, M.;
Genzel, R.;
Tacconi, L. J.;
Engel, H.;
Burkert, A.;
Krause, M.;
Sternberg, A.;
Sanchez, F. Mueller;
Maciejewski, W.
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]
oai:arXiv.org:0811.1425 [pdf] - 18346
The effect of stellar feedback on the formation and evolution of gas and
dust tori in AGN
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]
oai:arXiv.org:0808.0655 [pdf] - 15163
Turbulent AGN tori
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]
oai:arXiv.org:0802.2604 [pdf] - 10213
Three-dimensional radiative transfer models of clumpy tori in Seyfert
galaxies
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]
oai:arXiv.org:0709.0209 [pdf] - 4494
Resolving the complex structure of the dust torus in the active nucleus
of the Circinus galaxy
Tristram, K. R. W.;
Meisenheimer, K.;
Jaffe, W.;
Schartmann, M.;
Rix, H. -W.;
Leinert, Ch.;
Morel, S.;
Wittkowski, M.;
Röttgering, H.;
Perrin, G.;
Lopez, B.;
Raban, D.;
Cotton, W. D.;
Graser, U.;
Paresce, F.;
Henning, Th.
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]
oai:arXiv.org:astro-ph/0504105 [pdf] - 72202
Towards a physical model of dust tori in Active Galactic Nuclei -
Radiative transfer calculations for a hydrostatic torus model
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.