Source author record

R. I. Davies

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

32works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

32 published item(s)

preprint2022arXiv

From the Circumnuclear Disk in the Galactic Center to thick, obscuring tori of AGNs -- Modelling the molecular emission of a parsec-scale torus as found in NGC1068

The accretion rates needed to fuel the central black hole in a galaxy can be achieved via viscous torques in thick disks and rings, which can be resolved by millimetre interferometry within the inner ~20pc of the active galaxy NGC1068 at comparable scales and sensitivity to single dish observations of the Circumnuclear Disk (CND) in the Galactic Center. To interpret observations of these regions and determine the physical properties of their gas distribution, we present a modelling effort that includes (i) a simple dynamical simulations involving partially inelastic collisions between disk gas clouds, (ii) an analytical model of a turbulent clumpy gas disk calibrated by the dynamical model and observations, (iii) local turbulent and cosmic ray gas heating and cooling via H2O, H2, and CO emission, and (iv) determination of the molecular abundances. We also consider photodissociation regions (PDR) where gas is directly illuminated by the central engine. We compare the resulting model datacubes of the CO, HCN, HCO+, and CS brightness temperatures to available observations. In both cases the kinematics can be explained by one or two clouds colliding with a pre-existing ring, in a prograde sense for the CND and retrograde for NGC1068. And, with only dense disk clouds, the line fluxes can be reproduced to within a factor of about two. To avoid self-absorption of the intercloud medium, turbulent heating at the largest scales, comparable to the disk height, has to be decreased by a factor of 50-200. Our models indicate that turbulent mechanical energy input is the dominant gas heating mechanism within the thick gas disks. In N1068, while the bulk of the AGN X-ray radiation is absorbed in a layer of Compton-thick gas inside the dust sublimation radius, the optical/UV radiation may enhance the molecular line emission from photodissociation regions by ~50% at the inner edge of the gas ring.

preprint2020arXiv

Cold molecular gas and PAH emission in Seyfert galaxies

We investigate the relation between the detection of the $11.3\,μ$m PAH feature in the nuclear ($\sim 24-230\,$pc) regions of 22 nearby Seyfert galaxies and the properties of the cold molecular gas. For the former we use ground-based (0.3-0.6" resolution) mid-infrared (mid-IR) spectroscopy. The cold molecular gas is traced by ALMA and NOEMA high (0.2-1.1") angular resolution observations of the CO(2-1) transition. Galaxies with a nuclear detection of the $11.3\,μ$m PAH feature contain more cold molecular gas (median $1.6\times 10^7\,M_\odot$) and have higher column densities ($N({\rm H}_2) = 2 \times 10^{23}\,{\rm cm}^{-2}$) over the regions sampled by the mid-IR slits than those without a detection. This suggests that molecular gas plays a role in shielding the PAH molecules in the harsh environments of Seyfert nuclei. Choosing the PAH molecule naphthalene as an illustration, we compute its half-life in the nuclear regions of our sample when exposed to 2.5keV hard X-ray photons. We estimate shorter half-lives for naphthalene in nuclei without a $11.3\,μ$m PAH detection than in those with a detection. The Spitzer/IRS PAH ratios on circumnuclear scales ($\sim$ 4" $\sim$ 0.25-1.3kpc) are in between model predictions for neutral and partly ionized PAHs. However, Seyfert galaxies in our sample with the highest nuclear H$_2$ column densities are not generally closer to the neutral PAH tracks. This is because in the majority of our sample galaxies, the CO(2-1) emission in the inner $\sim$ 4" is not centrally peaked and in some galaxies traces circumnuclear sites of strong star formation activity. Spatially resolved observations with the MIRI medium-resolution spectrograph (MRS) on the James Webb Space Telescope will be able to distinguish the effects of an active galactic nucleus (AGN) and star formation on the PAH emission in nearby AGN.

preprint2020arXiv

From Nuclear to Circumgalactic: Zooming in on AGN-Driven Outflows at z~2.2 with SINFONI

We use deep adaptive optics assisted integral field spectroscopy from SINFONI on the VLT to study the spatially resolved properties of ionized gas outflows driven by active galactic nuclei (AGN) in three galaxies at z~2.2 -- K20-ID5, COS4-11337 and J0901+1814. These systems probe AGN feedback from nuclear to circumgalactic scales, and provide unique insights into the different mechanisms by which AGN-driven outflows interact with their host galaxies. K20-ID5 and COS4-11337 are compact star forming galaxies with powerful $\sim$1500 km s$^{-1}$ AGN-driven outflows that dominate their nuclear H$α$ emission. The outflows do not appear to have any impact on the instantaneous star formation activity of the host galaxies, but they carry a significant amount of kinetic energy which could heat the halo gas and potentially lead to a reduction in the rate of cold gas accretion onto the galaxies. The outflow from COS4-11337 is propagating directly towards its companion galaxy COS4-11363, at a projected separation of 5.4 kpc. COS4-11363 shows signs of shock excitation and recent truncation of star formation activity, which could plausibly have been induced by the outflow from COS4-11337. J0901+1814 is gravitationally lensed, giving us a unique view of a compact (R = 470 $\pm$ 70 pc), relatively low velocity ($\sim$650 km s$^{-1}$) AGN-driven outflow. J0901+1814 has a similar AGN luminosity to COS4-11337, suggesting that the difference in outflow properties is not related to the current AGN luminosity, and may instead reflect a difference in the evolutionary stage of the outflow and/or the coupling efficiency between the AGN ionizing radiation field and the gas in the nuclear regions.

preprint2020arXiv

Local AGN Survey (LASr): I. Galaxy sample, infrared colour selection and predictions for AGN within 100 Mpc

In order to answer some of the major open questions in the fields of supermassive black hole (SMBH) and galaxy evolution, a complete census of SMBH growth, i.e., active galactic nuclei (AGN), is required. Thanks to deep all-sky surveys, such as those by the Wide-field Infrared Survey Explorer (WISE) and the Spectrum-Roentgen-Gamma (SRG) missions, this task is now becoming feasible in the nearby Universe. Here, we present a new survey, the Local AGN Survey (LASr), with the goal of identifying AGN unbiased against obscuration and determining the intrinsic Compton-thick (CT) fraction. First, we construct the most complete all-sky sample of galaxies within 100 Mpc from astronomical databases (90% completeness for log M_star/M_sun ~9.4), four times deeper than the current local galaxy reference, the Two Micron All-Sky Survey Redshift Survey (2MRS), which turns out to miss 20% of known luminous AGN. These 49k galaxies serve as parent sample for LASr, called LASr-GPS. It contains 4.3k already known AGN, >=82% of these are estimated to have L(12um) < 10^42.3 erg/s, i.e., are low-luminosity AGN. As a first method for identifying Seyfert-like AGN, we use WISE-based infrared colours, finding 221 galaxies at L(12um) >= 10^42.3 erg/s to host an AGN at 90% reliability, This includes 61 new AGN candidates and implies and optical type 2 fraction of 50 to 71%. We quantify the efficiency of this technique and estimate the total number of AGN with L_int(2-10 keV) >= 10^42 erg/s in the volume to be 362^+145_-116 (8.6^+3.5_-2.8 10^-5 Mpc^-3). X-ray brightness estimates indicate the CT fraction to be 40-55% to explain the Swift non-detections of the infrared selected objects. One third of the AGN within 100 Mpc remain to be identified and we discuss the prospects for the eROSITA all-sky survey.

preprint2020arXiv

The regulation of galaxy growth along the size-mass relation by star-formation, as traced by H$α$ in KMOS$^{3D}$ galaxies at 0.7 < z < 2.7

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

preprint2019arXiv

Molecular outflows in local galaxies: Method comparison and a role of intermittent AGN driving

We report new detections and limits from a NOEMA and ALMA CO(1-0) search for molecular outflows in 13 local galaxies with high FIR surface brightness, and combine with results from the literature. CO line ratios and outflow structure provide some constraints on the conversion from observables to quantities such as molecular mass outflow rates. Ratios between outflow emission in higher J CO transitions and in CO(1-0) typically are consistent with excitation Ri1<~1. For IRAS 13120-5453, however, R31=2.10 indicates optically thin CO in the outflow. Like much of the outflow literature, we use alpha(CO) = 0.8, and we present arguments for using C=1 in deriving molecular mass outflow rates Mdot = C*M*v/R. We compare the two main methods for molecular outflow detection: CO mm interferometry and Herschel OH spectroscopy. For 26 sources studied with both methods, we find 80% agreement in detecting vout>~150km/s outflows, and non-matches can be plausibly ascribed to outflow geometry and SNR. For 12 bright ULIRGs with detailed OH-based outflow modeling, CO outflows are detected in all but one. Outflow masses, velocities, and sizes for these 11 sources agree well between the two methods, and modest remaining differences may relate to the different but overlapping regions sampled by CO emission and OH absorption. Outflow properties correlate better with AGN luminosity and with bolometric luminosity than with FIR surface brightness. The most massive outflows are found for systems with current AGN activity, but significant outflows in non-AGN systems must relate to star formation or to AGN activity in the recent past. We report scaling relations for the increase of outflow mass, rate, momentum rate, and kinetic power with bolometric luminosity. Short ~10^6yr flow times and some sources with resolved multiple outflow episodes support a role of intermittent driving, likely by AGN. (abridged)

preprint2016arXiv

A SINFONI view of the nuclear activity and circum-nuclear star formation in NGC 4303

We present new maps of emission-line flux distributions and kinematics in both ionized (traced by HI and [FeII] lines) and molecular (H2) gas of the inner 0.7x0.7kpc2 of the galaxy NGC4303, with a spatial resolution 40-80pc and velocity resolution 90-150 km/s obtained from near-IR integral field specroscopy using the VLT instrument SINFONI. The most promiment feature is a 200-250pc ring of circum-nuclear star-forming regions. The emission from ionized and molecular gas shows distinct flux distributions: while the strongest HI and [FeII] emission comes from regions in the west side of the ring (ages~4Myr), the H2 emission is strongest at the nucleus and in the east side of the ring (ages>10Myr). We find that regions of enhanced hot H2 emission are anti-correlated with those of enhanced [FeII] and HI emission, which can be attributed to post starburst regions that do not have ionizing photons anymore but still are hot enough (~2000K) to excite the H2 molecule. The line ratios are consistent with the presence of an AGN at the nucleus. The youngest regions have stellar masses in the range 0.3-1.5E5 MSun and ionized and hot molecular gas masses of ~0.25-1.2E4 Msun and 2.5-5 Msun, respectively. The stellar and gas velocity fields show a rotation pattern, with the gas presenting larger velocity amplitudes than the stars, with a deviation observed for the H2 along the nuclear bar, where increased velocity dispersion is also observed, possibly associated with non circular motions along the bar. The stars in the ring show smaller velocity dispersion than the surroundings, that can be attributed to a cooler dynamics due to their recent formation from cool gas.

preprint2016arXiv

Constraints on the Broad Line Region Properties and Extinction in Local Seyferts

We use high spectral resolution (R > 8000) data covering 3800-13000Å 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$α$/H$β$ 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.

preprint2016arXiv

On the relation of optical obscuration and X-ray absorption in Seyfert galaxies

The optical classification of a Seyfert galaxy and whether it is considered X-ray absorbed are often used interchangeably. But there are many borderline cases and also numerous examples where the optical and X-ray classifications appear to be in conflict. In this article we re-visit the relation between optical obscuration and X-ray absorption in AGNs. We make use of our "dust color" method (Burtscher et al. 2015) to derive the optical obscuration A_V and consistently estimated X-ray absorbing columns using 0.3--150 keV spectral energy distributions. We also take into account the variable nature of the neutral gas column N_H and derive the Seyfert sub-classes of all our objects in a consistent way. We show in a sample of 25 local, hard-X-ray detected Seyfert galaxies (log L_X / (erg/s) ~ 41.5 - 43.5) that there can actually be a good agreement between optical and X-ray classification. If Seyfert types 1.8 and 1.9 are considered unobscured, the threshold between X-ray unabsorbed and absorbed should be chosen at a column N_H = 10^22.3 / cm^2 to be consistent with the optical classification. We find that N_H is related to A_V and that the N_H/A_V ratio is approximately Galactic or higher in all sources, as indicated previously. But in several objects we also see that deviations from the Galactic ratio are only due to a variable X-ray column, showing that (1) deviations from the Galactic N_H/A_V can simply be explained by dust-free neutral gas within the broad line region in some sources, that (2) the dust properties in AGNs can be similar to Galactic dust and that (3) the dust color method is a robust way to estimate the optical extinction towards the sublimation radius in all but the most obscured AGNs.

preprint2016arXiv

Thick Disks, and an Outflow, of Dense Gas in the Nuclei of Nearby Seyfert Galaxies

We discuss the dense molecular gas in central regions of nearby Seyfert galaxies, and report new arcsec resolution observations of HCN(1-0) and HCO$^+$(1-0) for 3 objects. In NGC 3079 the lines show complex profiles as a result of self-absorption and saturated continuum absorption. H$^{13}$CN reveals the continuum absorption profile, with a peak close to the galaxy's systemic velocity that traces disk rotation, and a second feature with a blue wing extending to $-350$km s$^{-1}$ that most likely traces a nuclear outflow. The morphological and spectral properties of the emission lines allow us to constrain the dense gas dynamics. We combine our kinematic analysis for these 3 objects, as well as another with archival data, with a previous comparable analysis of 4 other objects, to create a sample of 8 Seyferts. In 7 of these, the emission line kinematics imply thick disk structures on radial scales of $\sim$100pc, suggesting such structures are a common occurrence. We find a relation between the circumnuclear LHCN and Mdyn that can be explained by a gas fraction of 10% and a conversion factor αHCN $\sim$ 10 between gas mass and HCN luminosity. Finally, adopting a different perspective to probe the physical properties of the gas around AGN, we report on an analysis of molecular line ratios which indicates that the clouds in this region are not self-gravitating.

preprint2015arXiv

Obscuration in AGNs: near-infrared luminosity relations and dust colors

We combine two approaches to isolate the AGN luminosity at near-infrared wavelengths and relate the near-IR pure AGN luminosity to other tracers of the AGN. Using integral-field spectroscopic data of an archival sample of 51 local AGNs, we estimate the fraction of non-stellar light by comparing the nuclear equivalent width of the stellar 2.3 micron CO absorption feature with the intrinsic value for each galaxy. We compare this fraction to that derived from a spectral decomposition of the integrated light in the central arc second and find them to be consistent with each other. Using our estimates of the near-IR AGN light, we find a strong correlation with presumably isotropic AGN tracers. We show that a significant offset exists between type 1 and type 2 sources in the sense that type 1 sources are 7 (10) times brighter in the near-IR at log L_MIR = 42.5 (log L_X = 42.5). These offsets only becomes clear when treating infrared type 1 sources as type 1 AGNs. All AGNs have very red near-to-mid-IR dust colors. This, as well as the range of observed near-IR temperatures, can be explained with a simple model with only two free parameters: the obscuration to the hot dust and the ratio between the warm and hot dust areas. We find obscurations of A_V (hot) = 5 - 15 mag for infrared type 1 sources and A_V (hot) = 15 - 35 mag for type 2 sources. The ratio of hot dust to warm dust areas of about 1000 is nicely consistent with the ratio of radii of the respective regions as found by infrared interferometry.

preprint2015arXiv

The SINS/zC-SINF Survey of z~2 Galaxy Kinematics: Rest-frame Morphology, Structure, and Colors from Near-infrared Hubble Space Telescope Imaging

We present the analysis of HST $J$- and $H$-band imaging for 29 galaxies on the star-forming main sequence at $z\sim2$, which have Adaptive Optics VLT SINFONI integral field spectroscopy from our SINS/zC-SINF program. The SINFONI H$α$ data resolve the on-going star-formation and the ionized gas kinematics on scales of $1-2$ kpc; the near-IR images trace the galaxies' rest-frame optical morphologies and distributions of stellar mass in old stellar populations at a similar resolution. The global light profiles of most galaxies show disk-like properties well described by a single Sérsic profile with $n\sim1$, with only $\sim15%$ requiring a high $n>3$ Sérsic index, all more massive than $10^{10}M_\odot$. In bulge+disk fits, about $40%$ of galaxies have a measurable bulge component in the light profiles, with $\sim15%$ showing a substantial bulge-to-total ratio $B/T\ge0.3$. This is a lower limit to the frequency of $z\sim2$ massive galaxies with a developed bulge component in stellar mass because it could be hidden by dust and/or outshined by a thick actively star-forming disk component. The galaxies' rest-optical half-light radii range between $1-7$ kpc, with a median of 2.1 kpc, and lie slightly above the size-mass relation at these epochs reported in the literature. This is attributed to differences in sample selection and definitions of size and/or mass measurements. The $(u-g)_{rest}$ color gradient and scatter within individual $z\sim2$ massive galaxies with $\ge10^{11}M_\odot$ are as high as in $z=0$ low-mass, late-type galaxies, and are consistent with the high star-formation rates of massive $z\sim2$ galaxies being sustained at large galactocentric distances.

preprint2014arXiv

Fueling AGN II: Spatially Resolved Molecular Inflows and Outflows

We analyse the 2-dimensional distribution and kinematics of the stars as well as molecular and ionised gas in the central few hundred parsecs of 5 active and 5 matched inactive galaxies. The equivalent widths of the Br-gamma line indicate there is no on-going star formation in their nuclei, although recent (terminated) starbursts are possible in the active galaxies. The stellar velocity fields show no signs of non-circular motions, while the 1-0S(1) H_2 kinematics exhibit significant deviations from simple circular rotation. In the active galaxies the H_2 kinematics reveal inflow and outflow superimposed on disk rotation. Steady-state circumnuclear inflow is seen in three AGN, and hydrodynamical models indicate it can be driven by a large scale bar. In three of the five AGN, molecular outflows are spatially resolved. The outflows are oriented such that they intersect, or have an edge close to, the disk - which may be the source of molecular gas in the outflow. The relatively low speeds imply the gas will fall back onto the disk; and with moderate outflow rates, they will have only a local impact on the host galaxy. H_2 was detected in two inactive galaxies. These exhibit chaotic circumnuclear dust morphologies and have molecular structures that are counter-rotating with respect to the main gas component, which could lead to gas inflow in the near future. In our sample, all four galaxies with chaotic dust morphology in the circumnuclear region exist in moderately dense groups with 10-15 members where accretion of stripped gas can easily occur.

preprint2013arXiv

The Black Hole and Central Stellar Population of MCG--6-30-15

We present the first near-infrared integral field spectroscopy observations of the galaxy MCG--6-30-15. The H-band data studied in this paper cover the central 500 pc of the galaxy at the best resolution (0".1) so far. The spectra of the innermost regions are dominated by broad brackett series emission lines and non-stellar continuum, under which we are able to trace the distribution and kinematics of the stars and also the [Fe II] line emission. We find that there is a counter-rotating stellar core extending out to 125 pc, which appears to be associated with the [Fe II] emission. Based on the mass-to-light ratio, and the presence of this emission line, we estimate the age of the central stellar population to be of order of 65 Myr. We show that the gas needed to fuel the black hole is, at most, only 1 per cent of that needed to form these stars. We derive independent constraints on the black hole mass using the dynamical information and determine an upper limit for the black hole mass, Mbh < 6e7 Msun, that is consistent with other estimates.

preprint2013arXiv

The Mid-infrared Emission of Narrow-Line Active Galactic Nuclei: Star-Formation, Nuclear Activity and two populations revealed by WISE

We explore the nature of the long-wavelength mid-infrared (MIR) emission of a sample of 13000 local Type II (narrow-line) Active Galactic Nuclei (AGNs) from the Sloan Digital Sky Survey (SDSS) using 12 and 22 micron photometry from the WISE all-sky survey. In combination with FIRST 1.4 GHz measurements, we show that AGNs divide into two relatively distinct populations or "branches" in the plane of MIR and radio luminosity. Seyfert galaxies lie almost exclusively on a MIR-bright branch (Branch A), while low-ionization nuclear emission line galaxies (LINERs) are split evenly into Branch A and the MIR-faint Branch B. We devise various tests to constrain the processes that define the branches, including a comparison to the properties of pure star-forming (SF) inactive galaxies on the MIR-Radio plane. We demonstrate that the total MIR emission of objects on Branch A, including most Seyfert galaxies, is governed primarily by host star-formation, with about 15% of the 22 micron luminosity coming from AGN-heated dust. This implies that on-going dusty star-formation is a general property of Seyfert host galaxies. We show that the 12 micron broad-band luminosity of AGNs on Branch A is suppressed with respect to star-forming galaxies, possibly due to the destruction of PAHs or deeper 10 microns Si absorption in AGNs. We uncover a correlation between the MIR luminosity and [O III] luminosity in AGNs. This suggests a relationship between the SFR and nuclear luminosity in the AGN population, but we caution on the importance of selection effects inherent to such AGN-dominated emission-line galaxies in driving such a correlation. We highlight the MIR-radio plane as a useful tool in comparative studies of SF and nuclear activity in AGN.

preprint2013arXiv

The Mrk 231 molecular outflow as seen in OH

We report on the Herschel/PACS observations of OH in Mrk 231, with detections in 9 doublets observed within the PACS range, and present radiative transfer models for the outflowing OH. Signatures of outflowing gas are found in up to 6 OH doublets with different excitation requirements. At least two outflowing components are identified, one with OH radiatively excited, and the other with low excitation, presumably spatially extended. Particularly prominent, the blue wing of the absorption detected in the in-ladder 2Pi_{3/2} J=9/2-7/2 OH doublet at 65 um, with E_lower=290 K, indicates that the excited outflowing gas is generated in a compact and warm (circum)nuclear region. Because the excited, outflowing OH gas in Mrk 231 is associated with the warm, far-IR continuum source, it is likely more compact (diameter of 200-300 pc) than that probed by CO and HCN. Nevertheless, its mass-outflow rate per unit of solid angle as inferred from OH is similar to that previously derived from CO, >~70x(2.5x10^{-6}/X_{OH}) Msun yr^{-1} sr^{-1}, where X_{OH} is the OH abundance relative to H nuclei. In spherical symmetry, this would correspond to >~850x(2.5x10^{-6}/X_{OH}) Msun yr^{-1}, though significant collimation is inferred from the line profiles. The momentum flux of the excited component attains ~15 L_{AGN}/c, with an OH column density of (1.5-3)x10^{17} cm^-2 and a mechanical luminosity of ~10^{11} Lsun. The detection of very excited OH peaking at central velocities indicates the presence of a nuclear reservoir of gas rich in OH, plausibly the 130-pc scale circumnuclear torus previously detected in OH megamaser emission, that may be feeding the outflow. An exceptional ^{18}OH enhancement, with OH/^{18}OH<~30 at both central and blueshifted velocities, is likely the result of interstellar-medium processing by recent starburst/SNe activity.

preprint2013arXiv

The quenching of star formation in accretion-driven clumpy turbulent tori of active galactic nuclei

Galactic gas-gas collisions involving a turbulent multiphase ISM share common ISM properties: dense extraplanar gas visible in CO, large linewidths (>= 50 km/s), strong mid-infrared H_2 line emission, low star formation activity, and strong radio continuum emission. Gas-gas collisions can occur in the form of ICM ram pressure stripping, galaxy head-on collisions, compression of the intragroup gas and/or galaxy ISM by an intruder galaxy which flies through the galaxy group at a high velocity, or external gas accretion on an existing gas torus in a galactic center. We suggest that the common theme of all these gas-gas interactions is adiabatic compression of the ISM leading to an increase of the turbulent velocity dispersion of the gas. The turbulent gas clouds are then overpressured and star formation is quenched. Within this scenario we developed a model for turbulent clumpy gas disks where the energy to drive turbulence is supplied by external infall or the gain of potential energy by radial gas accretion within the disk. The cloud size is determined by the size of a C-type shock propagating in dense molecular clouds with a low ionization fraction at a given velocity dispersion. We give expressions for the expected volume and area filling factors, mass, density, column density, and velocity dispersion of the clouds. The latter is based on scaling relations of intermittent turbulence whose open parameters are estimated for the CND in the Galactic Center. The properties of the model gas clouds and the external mass accretion rate necessary for the quenching of the star formation rate due to adiabatic compression are consistent with those derived from high-resolution H_2 line observations. Based on these findings, a scenario for the evolution of gas tori in galactic centers is proposed and the implications for star formation in the Galactic Center are discussed.

preprint2012arXiv

Physical properties of dense molecular gas in centres of Seyfert galaxies

We present new ~1" resolution data of the dense molecular gas in the central 50-100 pc of four nearby Seyfert galaxies. PdBI observations of HCN and, in 2 of the 4 sources, simultaneously HCO+ allow us to carefully constrain the dynamical state of the dense gas surrounding the AGN. Analysis of the kinematics shows large line widths of 100-200 km/s FWHM that can only partially arise from beam smearing of the velocity gradient. The observed morphological and kinematic parameters (dimensions, major axis position angle, red and blue channel separation, and integrated line width) are well reproduced by a thick disk, where the emitting dense gas has a large intrinsic dispersion (20-40 km/s), implying that it exists at significant scale heights (25-30% of the disk radius). To put the observed kinematics in the context of the starburst and AGN evolution, we estimate the Toomre Q parameter. We find this is always greater than the critical value, i.e. Q is above the limit such that the gas is stable against rapid star formation. This is supported by the lack of direct evidence, in these 4 Seyfert galaxies, for on-going star formation close around the AGN. Instead, any current star formation tends to be located in a circumnuclear ring. We conclude that the physical conditions are indeed not suited to star formation within the central ~100 pc.

preprint2011arXiv

An E-ELT Case Study: Colour-Magnitude Diagrams of an Old Galaxy in the Virgo Cluster

One of the key science goals for a diffraction limited imager on an Extremely Large Telescope (ELT) is the resolution of individual stars down to faint limits in distant galaxies. The aim of this study is to test the proposed capabilities of a multi-conjugate adaptive optics (MCAO) assisted imager working at the diffraction limit, in IJHK$_s$ filters, on a 42m diameter ELT to carry out accurate stellar photometry in crowded images in an Elliptical-like galaxy at the distance of the Virgo cluster. As the basis for realistic simulations we have used the phase A studies of the European-ELT project, including the MICADO imager (Davies & Genzel 2010) and the MAORY MCAO module (Diolaiti 2010). We convolved a complex resolved stellar population with the telescope and instrument performance expectations to create realistic images. We then tested the ability of the currently available photometric packages STARFINDER and DAOPHOT to handle the simulated images. Our results show that deep Colour-Magnitude Diagrams (photometric error, $\pm$0.25 at I$\ge$27.2; H$\ge$25. and K$_s\ge$24.6) of old stellar populations in galaxies, at the distance of Virgo, are feasible at a maximum surface brightness, $μ_V \sim$ 17 mag/arcsec$^2$ (down to M$_I > -4$ and M$_H \sim$ M$_K > -6$), and significantly deeper (photometric error, $\pm$0.25 at I$\ge$29.3; H$\ge$26.6 and K$_s\ge$26.2) for $μ_V \sim$ 21 mag/arcsec$^2$ (down to M$_I \ge -2$ and M$_H \sim$ M$_K \ge -4.5$). The photometric errors, and thus also the depth of the photometry should be improved with photometry packages specifically designed to adapt to an ELT MCAO Point Spread Function. We also make a simple comparison between these simulations and what can be expected from a Single Conjugate Adaptive Optics feed to MICADO and also the James Webb Space Telescope.

preprint2011arXiv

Arp220: extinction and merger-induced star formation

We analyse new spatially resolved integral field spectroscopic H- and K-band data at a resolution of 0.3arcsec (100pc) and re-analyse interferometric CO(2-1) line observations of the prototypical merging system Arp220. We find that the majority of the K-band luminosity is due to a 10Myr old starburst, with a significant contribution from an underlying ca 1Gyr old stellar population, and a small contribution from stars less than 8Myr old. The Calzetti (2000) reddening law provides the best fit to photometric datapoints spanning 0.45um -- 2.12um. Furthermore, estimates of the bolometric luminosity from IRAS fluxes in conjunction with our stellar population analysis indicate that we observe less than 10% of the emitted K-band light. The stellar and CO(2-1) kinematic centre of the western nucleus coincides with the compact hot dust emission, indicating that the latter marks the centre of the gravitational potential. In the eastern nucleus, the CO(2-1) data are well matched by a model in which the gas orbits around the peak of the dust emission. This, and the similarity of the K-band tracer kinematics, shows that despite the irregular morphology, the eastern nucleus is also a kinematically coherent structure. Comparison of the extinction map and EW_CO and EW_Br-gamma maps indicates that the lower half of the eastern nucleus is significantly more extincted than the upper half, suggesting that the lower half is buried in the larger-scale gas disk.

preprint2011arXiv

Outflows from AGN: Kinematics of the Narrow-Line and Coronal-Line Regions in Seyfert Galaxies

As part of an extensive study of the physical properties of active galactic nuclei (AGN) we report high spatial resolution near-IR integral-field spectroscopy of the narrow-line region (NLR) and coronal-line region (CLR) of seven Seyfert galaxies. These measurements elucidate for the first time the two-dimensional spatial distribution and kinematics of the recombination line Brγ and high-ionization lines [Sivi], [Alix] and [Caviii] on scales <300 pc from the AGN. The observations reveal kinematic signatures of rotation and outflow in the NLR and CLR. The spatially resolved kinematics can be modeled as a combination of an outflow bicone and a rotating disk coincident with the molecular gas. High-excitation emission is seen in both components, suggesting it is leaking out of a clumpy torus. While NGC 1068 (Seyfert 2) is viewed nearly edge-on, intermediate-type Seyferts are viewed at intermediate angles, consistent with unified schemes. A correlation between the outflow velocity and the molecular gas mass in r<30 pc indicates that the accumulation of gas around the AGN increases the collimation and velocity of the outflow. The outflow rate is 2-3 orders of magnitude greater than the accretion rate, implying that the outflow is mass-loaded by the surrounding interstellar medium (ISM). In half of the observed AGN the kinetic power of the outflow is of the order of the power required by two-stage feedback models to be thermally coupled to the ISM and match the M-σ* relation. In these objects the radio jet is clearly interacting with the ISM, indicative of a link between jet power and outflow power.

preprint2010arXiv

Adaptive optics near infrared integral field spectroscopy of NGC 2992

NGC 2992 is an intermediate Seyfert 1 galaxy showing outflows on kilo parsec scales which might be due either to AGN or starburst activity. We therefore aim at investigating its central region for a putative starburst in the past and its connection to the AGN and the outflows. Observations were performed with the adaptive optics near infrared integral field spectrograph SINFONI on the VLT, complemented by longslit observations with ISAAC on the VLT, as well as N- and Q-band data from the Spitzer archive. The spatial and spectral resolutions of the SINFONI data are 50 pc and 83 km/s, respectively. The field of view of 3" x 3" corresponds to 450 pc x 450 pc. Br_gamma equivalent width and line fluxes from PAHs were compared to stellar population models to constrain the age of the putative recent star formation. A simple geometric model of two mutually inclined disks and an additional cone to describe an outflow was developed to explain the observed complex velocity field in H_2 1-0S(1). The morphologies of the Br_gamma and the stellar continuum are different suggesting that at least part of the Br_gamma emission comes from the AGN. This is confirmed by PAH emission lines at 6.2 micron and 11.2 micron and the strength of the silicon absorption feature at 9.7 micron, which point to dominant AGN activity with a relatively minor starburst contribution. We find a starburst age of 40 Myr - 50 Myr from Br_gamma line diagnostics and the radio continuum; ongoing star formation can be excluded. Both the energetics and the timescales indicate that the outflows are driven by the AGN rather than the starburst. The complex velocity field observed in H_2 1-0S(1) in the central 450 pc can be explained by the superposition of the galaxy rotation and an outflow.

preprint2010arXiv

Gas dynamics of the central few parsec region of NGC 1068 fuelled by the evolving nuclear star cluster

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.

preprint2010arXiv

Most Submillimetre Galaxies are Major Mergers

We analyse subarcsecond resolution interferometric CO line data for twelve sub-millimetre-luminous (S850um > 5mJy) galaxies with redshifts between 1 and 3, presenting new data for four of them. Morphologically and kinematically most of the twelve systems appear to be major mergers. Five of them are well-resolved binary systems, and seven are compact or poorly resolved. Of the four binary systems for which mass measurements for both separate components can be made, all have mass ratios of 1:3 or closer. Furthermore, comparison of the ratio of compact to binary systems with that observed in local ULIRGs indicates that at least a significant fraction of the compact SMGs must also be late-stage mergers. In addition, the dynamical and gas masses we derive are most consistent with the lower end of the range of stellar masses published for these systems, favouring cosmological models in which SMGs result from mergers. These results all point to the same conclusion, that likely most of the bright SMGs with L_IR > 5x10e12L_sun are major mergers.

preprint2010arXiv

NGC6240: Merger-Induced Star Formation & Gas Dynamics

We present spatially resolved integral field spectroscopic K-band data at a resolution of 0.13" (60pc) and interferometric CO(2-1) line observations of the prototypical merging system NGC6240. Despite the clear rotational signature, the stellar kinematics in the two nuclei are dominated by dispersion. We use Jeans modelling to derive the masses and the mass-to-light ratios of the nuclei. Combining the luminosities with the spatially resolved Br-gamma equivalent width shows that only 1/3 of the K-band continuum from the nuclei is associated with the most recent star forming episode; and that less than 30% of the system's bolometric luminosity and only 9% of its stellar mass is due to this starburst. The star formation properties, calculated from typical merger star formation histories, demonstrate the impact of different assumptions about the star formation history. The properties of the nuclei, and the existence of a prominent old stellar population, indicate that the nuclei are remnants of the progenitor galaxies' bulges.

preprint2010arXiv

The Impact of cold gas accretion above a mass floor on galaxy scaling relations

Using the cosmological baryonic accretion rate and normal star formation efficiencies, we present a very simple model for star-forming galaxies (SFGs) that accounts for the mass and redshift dependencies of the SFR-Mass and Tully-Fisher relations from z=2 to the present. The time evolution follows from the fact that each modelled galaxy approaches a steady state where the SFR follows the (net) cold gas accretion rate. The key feature of the model is a halo mass floor M_{min}~10^{11} below which accretion is quenched in order to simultaneously account for the observed slopes of the SFR-Mass and Tully-Fischer relations. The same successes cannot be achieved via a star-formation threshold (or delay) nor by varying the SF efficiency or the feedback efficiency. Combined with the mass ceiling for cold accretion due to virial shock heating, the mass floor M_{min} explains galaxy "downsizing", where more massive galaxies formed earlier and over a shorter period of time. It turns out that the model also accounts for the observed galactic baryon and gas fractions as a function of mass and time, and the cosmic SFR density from z~6 to z=0, which are all resulting from the mass floor M_{min}. The model helps to understand that it is the cosmological decline of accretion rate that drives the decrease of cosmic SFR density between z~2 and z=0 and the rise of the cosmic SFR density allows us to put a constraint on our main parameter M_{min}~10^{11} solar masses. Among the physical mechanisms that could be responsible for the mass floor, we view that photo-ionization feedback (from first in-situ hot stars) lowering the cooling efficiency is likely to play a large role.

preprint2009arXiv

[The Impact of Nuclear Star Formation on Gas Inflow to AGN

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.

preprint2009arXiv

Obscuring and feeding supermassive black holes with evolving nuclear star clusters

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.

preprint2009arXiv

The Role of Molecular Gas in Obscuring Seyfert Active Galactic Nuclei

In a sample of local active galactic nuclei studied at a spatial resolution on the order of 10 pc we show that the interstellar medium traced by the molecular hydrogen v=1-0 S(1) 2.1um line forms a geometrically thick, clumpy disk. The kinematics of the molecular gas reveals general rotation, although an additional significant component of random bulk motion is required by the high local velocity dispersion. The size scale of the typical gas disk is found to have a radius of ~30 pc with a comparable vertical height. Within this radius the average gas mass is estimated to be ~10^7 Msun based on a typical gas mass fraction of 10%, which suggests column densities of Nh ~ 5x10^23 cm^-2. Extinction of the stellar continuum within this same region suggest lower column densities of Nh ~ 2x10^22 cm^-2, indicating that the gas distribution on these scales is dominated by dense clumps. In half of the observed Seyfert galaxies this lower column density is still great enough to obscure the AGN at optical/infrared wavelengths. We conclude, based on the spatial distribution, kinematics, and column densities that the molecular gas observed is spatially mixed with the nuclear stellar population and is likely to be associated with the outer extent of any smaller scale nuclear obscuring structure. Furthermore, we find that the velocity dispersion of the molecular gas is correlated with the star formation rate per unit area, suggesting a link between the two phenomena, and that the gas surface density follows known "Schmidt-Kennicutt" relations. The molecular/dusty structure on these scales may be dynamic since it is possible that the velocity dispersion of the gas, and hence the vertical disk height, is maintained by a short, massive inflow of material into the nuclear region and/or by intense, short-lived nuclear star formation.

preprint2008arXiv

Mergers and Mass Accretion Rates in Galaxy Assembly: The Millennium Simulation Compared to Observations of z~2 Galaxies

Recent observations of UV-/optically selected, massive star forming galaxies at z~2 indicate that the baryonic mass assembly and star formation history is dominated by continuous rapid accretion of gas and internal secular evolution, rather than by major mergers. We use the Millennium Simulation to build new halo merger trees, and extract halo merger fractions and mass accretion rates. We find that even for halos not undergoing major mergers the mass accretion rates are plausibly sufficient to account for the high star formation rates observed in z~2 disks. On the other hand, the fraction of major mergers in the Millennium Simulation is sufficient to account for the number counts of submillimeter galaxies (SMGs), in support of observational evidence that these are major mergers. When following the fate of these two populations in the Millennium Simulation to z=0, we find that subsequent mergers are not frequent enough to convert all z~2 turbulent disks into elliptical galaxies at z=0. Similarly, mergers cannot transform the compact SMGs/red sequence galaxies at z~2 into observed massive cluster ellipticals at z=0. We argue therefore, that secular and internal evolution must play an important role in the evolution of a significant fraction of z~2 UV-/optically and submillimeter selected galaxy populations.

preprint2008arXiv

Molecular gas streamers feeding and obscuring the active nucleus of NGC1068

We report the first direct observations of neutral, molecular gas streaming in the nucleus of NGC1068 on scales of <30 pc using SINFONI near-infrared integral field spectroscopy. At a resolution of 0.075", the flux map of 2.12 $μ$m 1-0 S(1) molecular hydrogen emission around the nucleus in the central arcsec reveals two prominent linear structures leading to the AGN from the north and south. The kinematics of the gas in these features are dominated by non-circular motions and indicate that material is streaming towards the nucleus on highly elliptical or parabolic trajectories whose orientations are compatible with that of the disk plane of the galaxy. We interpret the data as evidence for fueling of gas to the central region. The radial transport rate from ~30 pc to a few parsec from the nucleus is ~15 M$_\sun$ yr$^{-1}$. One of the infalling clouds lies directly in front of the central engine. We interpret it as a tidally disrupted streamer that forms the optically thick outerpart of an amorphous clumpy molecular/dusty structure which contributes to the nuclear obscuration.

preprint1998arXiv

ALFA: First Operational Experience of the MPE/MPIA Laser Guide Star System for Adaptive Optics

The sodium laser guide star adaptive optics system ALFA has been constructed at the Calar Alto 3.5-m telescope. Following the first detection of the laser beacon on the wavefront sensor in 1997 the system is now being optimized for best performance. In this contribution we discuss the current status of the launch beam and the planned improvements and upgrades. We report on the performance level achieved when it is used with the adaptive optics system, and relate various aspects of our experience during operation of the system. We have begun to produce scientific results and mention two of these.