Source author record

A. Sternberg

A. Sternberg 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

66works
7topics
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

66 published item(s)

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

Ionized outflows in local luminous AGN: what are the real densities and outflow rates?

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.

preprint2019arXiv

An image of the dust sublimation region in the nucleus of NGC 1068

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

preprint2019arXiv

Molecular Gas Inflows and Outflows in Ultraluminous Infrared Galaxies at $z\sim0.2$ and one QSO at $z=6.1$

Aims. We aim to search and characterize inflows and outflows of molecular gas in four ultraluminous infrared galaxies (ULIRGs) at $z\sim0.2-0.3$ and one distant QSO at $z=6.13$. Methods. We use Herschel PACS and ALMA Band 7 observations of the hydroxyl molecule (OH) line at rest-frame wavelength 119 $μ$m which in absorption can provide unambiguous evidence for inflows or outflows of molecular gas in nuclear regions of galaxies. Our study contributes to double the number of OH observations of luminous systems at $z\sim0.2-0.3$, and push the search for molecular outflows based on the OH transition to $z\sim6$. Results. We detect OH high-velocity absorption wings in three of the four ULIRGs. In two cases, IRAS F20036-1547 and IRAS F13352+6402, the blueshifted absorption profiles indicate the presence of powerful and fast molecular gas outflows. Consistent with an inside-out quenching scenario, these outflows are depleting the central reservoir of molecular gas at a similar rate than the intense star formation activity. In the case of the starburst-dominated system IRAS 10091+4704, we detect an inverted P-Cygni profile that is unique among ULIRGs and indicates the presence of a fast ($\sim400$ km s$^{-1}$) inflow of molecular gas at a rate of $\sim100~M_{\odot}~{\rm yr}^{-1}$ towards the central region. Finally, we tentatively detect ($\sim3σ$) the OH doublet in absorption in the $z=6.13$ QSO ULAS J131911+095051. The OH feature is blueshifted with a median velocity that suggests the presence of a molecular outflow, although characterized by a modest molecular mass loss rate of $\sim200~M_{\odot}~{\rm yr}^{-1}$. This value is comparable to the small mass outflow rates found in the stacking of the [CII] spectra of other $z\sim6$ QSOs and suggests that ejective feedback in this phase of the evolution of ULAS J131911+095051 has subsided.

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)

preprint2019arXiv

The resolved size and structure of hot dust in the immediate vicinity of AGN

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.

preprint2016arXiv

Broad [CII] line wings as tracer of molecular and multi-phase outflows in infrared bright galaxies

We report a tentative correlation between the outflow characteristics derived from OH absorption at $119\,μ\text{m}$ and [CII] emission at $158\,μ\text{m}$ in a sample of 22 local and bright ultraluminous infrared galaxies (ULIRGs). For this sample we investigate whether [CII] broad wings are a good tracer of molecular outflows, and how the two tracers are connected. Fourteen objects in our sample have a broad wing component as traced by [CII], and all of these also show OH119 absorption indicative of an outflow (in 1 case an inflow). The other eight cases, where no broad [CII] component was found, are predominantly objects with no OH outflow or a low-velocity ($\leq 100\,\text{km s}^{-1}$) OH outflow. The full width at half maximum (FWHM) of the broad [CII] component shows a trend with the OH119 blue-shifted velocity, although with significant scatter. Moreover, and despite large uncertainties, the outflow masses derived from OH and broad [CII] show a 1:1 relation. The main conclusion is therefore that broad [CII] wings can be used to trace molecular outflows. This may be particularly relevant at high redshift, where the usual tracers of molecular gas (like low-J CO lines) become hard to observe. Additionally, observations of blue-shifted Na I D $λλ5890,5896$ absorption are available for ten of our sources. Outflow velocities of Na I D show a trend with OH velocity and broad [CII] FWHM. These observations suggest that the atomic and molecular gas phases of the outflow are connected.

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

The angular momentum distribution and baryon content of star forming galaxies at z~1-3

We analyze the angular momenta of massive star forming galaxies (SFGs) at the peak of the cosmic star formation epoch (z~0.8-2.6). Our sample of ~360 log(M*/Msun) ~ 9.3-11.8 SFGs is mainly based on the KMOS3D and SINS/zC-SINF surveys of H$α$ kinematics, and collectively provides a representative subset of the massive star forming population. The inferred halo scale angular momentum distribution is broadly consistent with that theoretically predicted for their dark matter halos, in terms of mean spin parameter <$λ$> ~ 0.037 and its dispersion ($σ_{log(λ)}$~0.2). Spin parameters correlate with the disk radial scale, and with their stellar surface density, but do not depend significantly on halo mass, stellar mass, or redshift. Our data thus support the long-standing assumption that on average, even at high redshifts, the specific angular momentum of disk galaxies reflects that of their dark matter halos (j_d = j_DM). The lack of correlation between $λ$ x (j_d/j_DM) and the nuclear stellar density $Σ_{*}$(1kpc) favors a scenario where disk-internal angular momentum redistribution leads to "compaction" inside massive high-redshift disks. For our sample, the inferred average stellar-to-dark matter mass ratio is ~2%, consistent with abundance matching results. Including the molecular gas, the total baryonic disk-to-dark matter mass ratio is ~5% for halos near $10^{12}$ Msun, which corresponds to 31% of the cosmologically available baryons, implying that high-redshift disks are strongly baryon dominated.

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

500 Days of SN 2013dy: spectra and photometry from the ultraviolet to the infrared

SN 2013dy is a Type Ia supernova for which we have compiled an extraordinary dataset spanning from 0.1 to ~ 500 days after explosion. We present 10 epochs of ultraviolet (UV) through near-infrared (NIR) spectra with HST/STIS, 47 epochs of optical spectra (15 of them having high resolution), and more than 500 photometric observations in the BVrRiIZYJH bands. SN 2013dy has a broad and slowly declining light curve (delta m(B) = 0.92 mag), shallow Si II 6355 absorption, and a low velocity gradient. We detect strong C II in our earliest spectra, probing unburned progenitor material in the outermost layers of the SN ejecta, but this feature fades within a few days. The UV continuum of SN 2013dy, which is strongly affected by the metal abundance of the progenitor star, suggests that SN 2013dy had a relatively high-metallicity progenitor. Examining one of the largest single set of high-resolution spectra for a SN Ia, we find no evidence of variable absorption from circumstellar material. Combining our UV spectra, NIR photometry, and high-cadence optical photometry, we construct a bolometric light curve, showing that SN 2013dy had a maximum luminosity of 10.0^{+4.8}_{-3.8} * 10^{42} erg/s. We compare the synthetic light curves and spectra of several models to SN 2013dy, finding that SN 2013dy is in good agreement with a solar-metallicity W7 model.

preprint2015arXiv

A deep Herschel/PACS observation of CO(40-39) in NGC 1068: a search for the molecular torus

Emission from high-J CO lines in galaxies has long been proposed as a tracer of X-ray dominated regions (XDRs) produced by AGN. Of particular interest is the question of whether the obscuring torus, which is required by AGN unification models, can be observed via high-J CO cooling lines. Here we report on the analysis of a deep Herschel-PACS observation of an extremely high J CO transition (40-39) in the Seyfert 2 galaxy NGC 1068. The line was not detected, with a derived 3$σ$ upper limit of $2 \times 10^{-17}\,\text{W}\,\text{m}^{-2}$. We apply an XDR model in order to investigate whether the upper limit constrains the properties of a molecular torus in NGC 1068. The XDR model predicts the CO Spectral Line Energy Distributions for various gas densities and illuminating X-ray fluxes. In our model, the CO(40-39) upper limit is matched by gas with densities $\sim 10^{6}-10^{7}\,\text{cm}^{-3}$, located at $1.6-5\,\text{pc}$ from the AGN, with column densities of at least $10^{25}\,\text{cm}^{-2}$. At such high column densities, however, dust absorbs most of the CO(40-39) line emission at $λ= 65.69\, μ$m. Therefore, even if NGC 1068 has a molecular torus which radiates in the CO(40-39) line, the dust can attenuate the line emission to below the PACS detection limit. The upper limit is thus consistent with the existence of a molecular torus in NGC 1068. In general, we expect that the CO(40-39) is observable in only a few AGN nuclei (if at all), because of the required high gas column density, and absorption by dust.

preprint2015arXiv

High Resolution Imaging of PHIBSS z~2 Main Sequence Galaxies in CO J=1-0

We present Karl G. Jansky Very Large Array observations of the CO J=1-0 transition in a sample of four $z\sim2$ main sequence galaxies. These galaxies are in the blue sequence of star-forming galaxies at their redshift, and are part of the IRAM Plateau de Bure HIgh-$z$ Blue Sequence Survey (PHIBSS) which imaged them in CO J=3-2. Two galaxies are imaged here at high signal-to-noise, allowing determinations of their disk sizes, line profiles, molecular surface densities, and excitation. Using these and published measurements, we show that the CO and optical disks have similar sizes in main-sequence galaxies, and in the galaxy where we can compare CO J=1-0 and J=3-2 sizes we find these are also very similar. Assuming a Galactic CO-to-H$_2$ conversion, we measure surface densities of $Σ_{mol}\sim1200$ M$_\odot$pc$^{-2}$ in projection and estimate $Σ_{mol}\sim500-900$ M$_\odot$pc$^{-2}$ deprojected. Finally, our data yields velocity-integrated Rayleigh-Jeans brightness temperature line ratios $r_{31}$ that are approximately unity. In addition to the similar disk sizes, the very similar line profiles in J=1-0 and J=3-2 indicate that both transitions sample the same kinematics, implying that their emission is coextensive. We conclude that in these two main sequence galaxies there is no evidence for significant excitation gradients or a large molecular reservoir that is diffuse or cold and not involved in active star-formation. We suggest that $r_{31}$ in very actively star-forming galaxies is likely an indicator of how well mixed the star formation activity and the molecular reservoir are.

preprint2015arXiv

High-J CO SLEDs in nearby infrared bright galaxies observed by Herschel-PACS

We report the detection of far-infrared (FIR) CO rotational emission from nearby active galactic nuclei (AGN) and starburst galaxies, as well as several merging systems and Ultra-Luminous Infrared Galaxies (ULIRGs). Using Herschel-PACS, we have detected transitions in the J$_{upp}$ = 14 - 20 range ($λ\sim$ 130 - 185 $μ$m, $ν\sim$ 1612 - 2300 GHz) with upper limits on (and in two cases, detections of) CO line fluxes up to J$_{upp}$ = 30. The PACS CO data obtained here provide the first well-sampled FIR extragalactic CO SLEDs for this range, and will be an essential reference for future high redshift studies. We find a large range in the overall SLED shape, even amongst galaxies of similar type, demonstrating the uncertainties in relying solely on high-J CO diagnostics to characterize the excitation source of a galaxy. Combining our data with low-J line intensities taken from the literature, we present a CO ratio-ratio diagram and discuss its potential diagnostic value in distinguishing excitation sources and physical properties of the molecular gas. The position of a galaxy on such a diagram is less a signature of its excitation mechanism, than an indicator of the presence (or absence) of warm, dense molecular gas. We then quantitatively analyze the CO emission from a subset of the detected sources with Large Velocity Gradient (LVG) radiative transfer models to fit the CO SLEDs. Using both single-component and two-component LVG models to fit the kinetic temperature, velocity gradient, number density and column density of the gas, we derive the molecular gas mass and the corresponding CO-to-H$_2$ conversion factor, $α_{CO}$, for each respective source. For the ULIRGs we find $α$ values in the canonical range 0.4 - 5 M$_\odot$/(K kms$^{-1}$pc$^2$), while for the other objects, $α$ varies between 0.2 and 14.} Finally, we compare our best-fit LVG model ..

preprint2015arXiv

Insights on the Dusty Torus and Neutral Torus from Optical and X-ray Obscuration in a Complete Volume Limited Hard X-ray AGN Sample

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.

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.

preprint2014arXiv

A multi-wavelength investigation of the radio-loud supernova PTF11qcj and its circumstellar environment

We present the discovery, classification, and extensive panchromatic (from radio to X-ray) follow-up observations of PTF11qcj, a supernova discovered by the Palomar Transient Factory. PTF11qcj is located at a distance of dL ~ 124 Mpc. Our observations with the Karl G. Jansky Very Large Array show that this event is radio-loud: PTF11qcj reached a radio peak luminosity comparable to that of the famous gamma-ray-burst-associated supernova 1998bw (L_{5GHz} ~ 10^{29} erg/s/Hz). PTF11qcj is also detected in X-rays with the Chandra observatory, and in the infrared band with Spitzer. Our multi-wavelength analysis probes the supernova interaction with circumstellar material. The radio observations suggest a progenitor mass-loss rate of ~10^{-4} Msun/yr x (v_w/1000 km/s), and a velocity of ~(0.3-0.5)c for the fastest moving ejecta (at ~10d after explosion). However, these estimates are derived assuming the simplest model of supernova ejecta interacting with a smooth circumstellar material characterized by radial power-law density profile, and do not account for possible inhomogeneities in the medium and asphericity of the explosion. The radio light curve shows deviations from such a simple model, as well as a re-brightening at late times. The X-ray flux from PTF11qcj is compatible with the high-frequency extrapolation of the radio synchrotron emission (within the large uncertainties). An IR light echo from pre-existing dust is in agreement with our infrared data. Our analysis of pre-explosion data from the Palomar Transient Factory suggests that a precursor eruption of absolute magnitude M_r ~ -13 mag may have occurred ~ 2.5 yr prior to the supernova explosion. Based on our panchromatic follow-up campaign, we conclude that PTF11qcj fits the expectations from the explosion of a Wolf-Rayet star. Precursor eruptions may be a feature characterizing the final pre-explosion evolution of such stars.

preprint2014arXiv

AGC198606: A gas-bearing dark matter minihalo?

We present neutral hydrogen (HI) imaging observations with the Westerbork Synthesis Radio Telescope of AGC198606, an HI cloud discovered in the ALFALFA 21cm survey. This object is of particular note as it is located 16 km/s and 1.2 degrees from the gas-bearing ultra-faint dwarf galaxy Leo T while having a similar HI linewidth and approximately twice the flux density. The HI imaging observations reveal a smooth, undisturbed HI morphology with a full extent of 23'x16' at the 5x10^18 atoms cm^-2 level. The velocity field of AGC198606 shows ordered motion with a gradient of ~25 km/s across ~20'. The global velocity dispersion is 9.3 km/s with no evidence for a narrow spectral component. No optical counterpart to AGC198606 is detected. The distance to AGC198606 is unknown, and we consider several different scenarios: physical association with Leo T, a minihalo at a distance of ~150 kpc based on the models of Faerman et al. (2013), and a cloud in the Galactic halo. At a distance of 420 kpc, AGC198606 would have an HI mass of 6.2x10^5 Msun, an HI radius of 1.4 kpc, and a dynamical mass within the HI extent of 1.5x10^8 Msun.

preprint2014arXiv

Combined CO & Dust Scaling Relations of Depletion Time and Molecular Gas Fractions with Cosmic Time, Specific Star Formation Rate and Stellar Mass

We combine molecular gas masses inferred from CO emission in 500 star forming galaxies (SFGs) between z=0 and 3, from the IRAM-COLDGASS, PHIBSS1/2 and other surveys, with gas masses derived from Herschel far-IR dust measurements in 512 galaxy stacks over the same stellar mass/redshift range. We constrain the scaling relations of molecular gas depletion time scale (tdepl) and gas to stellar mass ratio (Mmolgas/M*) of SFGs near the star formation main-sequence with redshift, specific star formation rate (sSFR) and stellar mass (M*). The CO- and dust-based scaling relations agree remarkably well. This suggests that the CO-H2 mass conversion factor varies little within 0.6dex of the main sequence (sSFR(ms,z,M*)), and less than 0.3dex throughout this redshift range. This study builds on and strengthens the results of earlier work. We find that tdepl scales as (1+z)^-0.3 *(sSFR/sSFR(ms,z,M*))^-0.5, with little dependence on M*. The resulting steep redshift dependence of Mmolgas/M* ~(1+z)^3 mirrors that of the sSFR and probably reflects the gas supply rate. The decreasing gas fractions at high M* are driven by the flattening of the SFR-M* relation. Throughout the redshift range probed a larger sSFR at constant M* is due to a combination of an increasing gas fraction and a decreasing depletion time scale. As a result galaxy integrated samples of the Mmolgas-SFR rate relation exhibit a super-linear slope, which increases with the range of sSFR. With these new relations it is now possible to determine Mmolgas with an accuracy of 0.1dex in relative terms, and 0.2dex including systematic uncertainties.

preprint2014arXiv

Evidence for Wide-Spread AGN Driven Outflows in the Most Massive z~1-2 Star Forming Galaxies

In this paper we follow up on our previous detection of nuclear ionized outflows in the most massive (log(M*/Msun) >= 10.9) z~1-3 star-forming galaxies (Forster Schreiber et al.), by increasing the sample size by a factor of six (to 44 galaxies above log(M*/Msun) >= 10.9) from a combination of the SINS/zC-SINF, LUCI, GNIRS, and KMOS^3D spectroscopic surveys. We find a fairly sharp onset of the incidence of broad nuclear emission (FWHM in the Ha, [NII], and [SII] lines ~ 450-5300 km/s), with large [NII]/Ha ratios, above log(M*/Msun) ~ 10.9, with about two thirds of the galaxies in this mass range exhibiting this component. Broad nuclear components near and above the Schechter mass are similarly prevalent above and below the main sequence of star-forming galaxies, and at z~1 and ~2. The line ratios of the nuclear component are fit by excitation from active galactic nuclei (AGN), or by a combination of shocks and photoionization. The incidence of the most massive galaxies with broad nuclear components is at least as large as that of AGNs identified by X-ray, optical, infrared or radio indicators. The mass loading of the nuclear outflows is near unity. Our findings provide compelling evidence for powerful, high-duty cycle, AGN-driven outflows near the Schechter mass, and acting across the peak of cosmic galaxy formation.

preprint2014arXiv

The SINS/zC-SINF survey of z~2 galaxy kinematics: Evidence for powerful AGN-driven nuclear outflows in massive star-forming galaxies

We report the detection of ubiquitous powerful nuclear outflows in massive (> 10^11 Msun) z~2 star-forming galaxies (SFGs), which are plausibly driven by an Active Galactic Nucleus (AGN). The sample consists of the eight most massive SFGs from our SINS/zC-SINF survey of galaxy kinematics with the imaging spectrometer SINFONI, six of which have sensitive high-resolution adaptive optics (AO) assisted observations. All of the objects are disks hosting a significant stellar bulge. The spectra in their central regions exhibit a broad component in Halpha and forbidden [NII] and [SII] line emission, with typical velocity FWHM ~ 1500 km/s, [NII]/Halpha ratio ~ 0.6, and intrinsic extent of 2 - 3 kpc. These properties are consistent with warm ionized gas outflows associated with Type 2 AGN, the presence of which is confirmed via independent diagnostics in half the galaxies. The data imply a median ionized gas mass outflow rate of ~ 60 Msun/yr and mass loading of ~ 3. At larger radii, a weaker broad component is detected but with lower FWHM ~ 485 km/s and [NII]/Halpha ~ 0.35, characteristic for star formation-driven outflows as found in the lower-mass SINS/zC-SINF galaxies. The high inferred mass outflow rates and frequent occurrence suggest the nuclear outflows efficiently expel gas out of the centers of the galaxies with high duty cycles, and may thus contribute to the process of star formation quenching in massive galaxies. Larger samples at high masses will be crucial to confirm the importance and energetics of the nuclear outflow phenomenon, and its connection to AGN activity and bulge growth.

preprint2013arXiv

Dielectric properties of PLZT-x/65/35 (2 \leqslant x \leqslant 13) under mechanical stress, electric field and temperature loading

We investigated the effect of uniaxial pressure (0-1000 bars) applied parallely to the ac electric field on dielectric properties of PLZT-x/65/35 (2 \leqslant x \leqslant 13) ceramics. There was revealed a significant effect of the external stress on these properties. The application of uniaxial pressure leads to the change of the peak intensity of the electric permittivity (\varepsilon), of the frequency dispersion as well as of the dielectric hysteresis. The peak intensity $\varepsilon$ becomes diffused/sharpened and shifts to a higher/lower temperatures with increasing the pressure. It was concluded that the application of uniaxial pressure induces similar effects as increasing the Ti ion concentration in PZT system. We interpreted our results based on the domain switching processes under the action of combined electromechanical loading.

preprint2013arXiv

Fast Molecular Outflows in Luminous Galaxy Mergers: Evidence for Quasar Feedback from Herschel

We report the results from a systematic search for molecular (OH-119 um) outflows with Herschel-PACS in a sample of 43 nearby (z < 0.3) galaxy mergers, mostly ultraluminous infrared galaxies (ULIRGs) and QSOs. We find that the character of the OH feature (strength of the absorption relative to the emission) correlates with that of the 9.7-um silicate feature, a measure of obscuration in ULIRGs. Unambiguous evidence for molecular outflows, based on the detection of OH absorption profiles with median velocities more blueshifted than -50 km/sec, is seen in 26 (70%) of the 37 OH-detected targets, suggesting a wide-angle (~145 degrees) outflow geometry. Conversely, unambiguous evidence for molecular inflows, based on the detection of OH absorption profiles with median velocities more redshifted than +50 km/sec, is seen in only 4 objects, suggesting a planar or filamentary geometry for the inflowing gas. Terminal outflow velocities of ~-1000 km/sec are measured in several objects, but median outflow velocities are typically ~-200 km s^{-1}. While the outflow velocities show no statistically significant dependence on the star formation rate, they are distinctly more blueshifted among systems with large AGN fractions and luminosities [log (L_AGN / L_sun) > 11.8 +/- 0.3]. The quasars in these systems play a dominant role in driving the molecular outflows. In contrast, the most AGN dominated systems, where OH is seen purely in emission, show relatively modest OH line widths, despite their large AGN luminosities, perhaps indicating that molecular outflows subside once the quasar has cleared a path through the obscuring material.

preprint2013arXiv

Multi-Epoch High-Spectral-Resolution Observations of Neutral Sodium in 14 Type Ia Supernovae

One of the main questions concerning Type Ia supernovae is the nature of the binary companion of the exploding white dwarf. A major discriminant between different suggested models is the presence and physical properties of circumstellar material at the time of explosion. If present, this material will be ionized by the ultra-violet radiation of the explosion and later recombine. This ionization-recombination should manifest itself as time-variable absorption features that can be detected via multi-epoch high-spectral-resolution observations. Previous studies have shown that the strongest effect is seen in the neutral sodium D lines. We report on observations of neutral sodium absorption features observed in multi-epoch high-resolution spectra of 14 Type Ia supernova events. This is the first multi-epoch high-resolution study to include multiple SNe. No variability in line strength that can be associated with circumstellar material is detected. We find that ~18% of the events in the extended sample exhibit time-variable sodium features associated with circumstellar material. We explore the implication of this study on our understanding of the progenitor systems of Type Ia supernovae via the current Type Ia supernova multi-epoch high-spectral-resolution sample.

preprint2013arXiv

Nebular Excitation in z~2 Star-forming Galaxies from the SINS and LUCI Surveys: The Influence of Shocks and AGN

Based on high-resolution, spatially resolved data of 10 z~2 star-forming galaxies from the SINS/zC-SINF survey and LUCI data for 12 additional galaxies, we probe the excitation properties of high-z galaxies and the impact of active galactic nuclei (AGN), shocks and photoionization. We explore how these spatially-resolved line ratios can inform our interpretation of integrated emission line ratios obtained at high redshift. Many of our galaxies fall in the `composite' region of the z~0 [NII]/Ha versus [OIII]/Hb diagnostic (BPT) diagram, between star-forming galaxies and those with AGN. Based on our resolved measurements, we find that some of these galaxies likely host an AGN, while others appear to be affected by the presence of shocks possibly caused by an outflow or from enhanced ionization parameter as compared with HII regions in normal local star-forming galaxies. We find that the Mass-Excitation (MEx) diagnostic, which separates purely star-forming and AGN hosting local galaxies in the [OIII]/Hb versus stellar mass plane, does not properly separate z~2 galaxies classified according to the BPT diagram. However, if we shift the galaxies based on the offset between the local and z~2 mass-metallicity relation (i.e. to the mass they would have at z~0 with the same metallicity), we find better agreement between the MEx and BPT diagnostics. Finally, we find that metallicity calibrations based on [NII]/Ha are more biased by shocks and AGN at high-z than the [OIII]/Hb/[NII]/Ha calibration.

preprint2013arXiv

On the Source of the Dust Extinction in Type Ia Supernovae and the Discovery of Anomalously Strong Na I Absorption

High-dispersion observations of the Na I D 5890, 5896 and K I 7665, 7699 interstellar lines, and the diffuse interstellar band at 5780 Angstroms in the spectra of 32 Type Ia supernovae are used as an independent means of probing dust extinction. We show that the dust extinction of the objects where the diffuse interstellar band at 5780 Angstroms is detected is consistent with the visual extinction derived from the supernova colors. This strongly suggests that the dust producing the extinction is predominantly located in the interstellar medium of the host galaxies and not in circumstellar material associated with the progenitor system. One quarter of the supernovae display anomalously large Na I column densities in comparison to the amount of dust extinction derived from their colors. Remarkably, all of the cases of unusually strong Na I D absorption correspond to "Blueshifted" profiles in the classification scheme of Sternberg et al. (2011). This coincidence suggests that outflowing circumstellar gas is responsible for at least some of the cases of anomalously large Na I column densities. Two supernovae with unusually strong Na I D absorption showed essentially normal K I column densities for the dust extinction implied by their colors, but this does not appear to be a universal characteristic. Overall, we find the most accurate predictor of individual supernova extinction to be the equivalent width of the diffuse interstellar band at 5780 Angstroms, and provide an empirical relation for its use. Finally, we identify ways of producing significant enhancements of the Na abundance of circumstellar material in both the single-degenerate and double-degenerate scenarios for the progenitor system.

preprint2013arXiv

PHIBSS: molecular gas, extinction, star formation and kinematics in the z=1.5 star forming galaxy EGS13011166

We report matched resolution, imaging spectroscopy of the CO J=3-2 line (with the IRAM Plateau de Bure millimeter interferometer) and of the H-alpha line (with LUCI at the Large Binocular Telescope)in the massive z=1.53 main-sequence galaxy EGS 13011166, as part of the "Plateau de Bure high-z, blue sequence survey (PHIBSS). We combine these data with HST V-J-J-H-band maps to derive spatially resolved distributions of stellar surface density, star formation rate, molecular gas surface density, optical extinction and gas kinematics. The spatial distribution and kinematics of the ionized and molecular gas are remarkably similar and are well modeled by a turbulent, globally Toomre unstable rotating disk. The stellar surface density distribution is smoother than the clumpy rest-frame UV/optical light distribution, and peaks in an obscured, star forming massive bulge near the dynamical center. The molecular gas surface density and the effective optical screen extinction track each other and are well modeled by a 'mixed' extinction model. The inferred slope of the spatially resolved molecular gas to star formation rate relation depends strongly on the adopted extinction model and can vary from 0.8 to 1.7. For the preferred mixed dust-gas model we find a near linear slope.

preprint2013arXiv

The SINS/zC-SINF survey of z~2 galaxy kinematics: evidence for gravitational quenching

As part of the SINS/zC-SINF surveys of high-z galaxy kinematics, we derive the radial distributions of H-alpha surface brightness, stellar mass surface density, and dynamical mass at ~2 kpc resolution in 19 z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. From these data we infer the radial distribution of the Toomre Q-parameter for these main-sequence star forming galaxies (SFGs), covering almost two decades of stellar mass (10^9.6 to 10^11.5 solar masses). In more than half of our SFGs, the H-alpha distributions cannot be fit by a centrally peaked distribution, such as an exponential, but are better described by a ring, or the combination of a ring and an exponential. At the same time the kinematic data indicate the presence of a mass distribution more centrally concentrated than a single exponential distribution for 5 of the 19 galaxies. The resulting Q-distributions are centrally peaked for all, and significantly exceed unity there for three quarters of the SFGs. The occurrence of H-alpha rings and of large nuclear Q-values is strongly correlated, and is more common for the more massive SFGs. While our sample is small and there remain substantial uncertainties and caveats, our observations are consistent with a scenario in which cloud fragmentation and global star formation are secularly suppressed in gas rich high-z disks from the inside out, as the central stellar mass density of the disks grows.

preprint2013arXiv

Towards a resolved Kennicutt-Schmidt law at high redshift

Massive galaxies in the distant Universe form stars at much higher rates than today. Although direct resolution of the star forming regions of these galaxies is still a challenge, recent molecular gas observations at the IRAM Plateau de Bure interferometer enable us to study the star formation efficiency on subgalactic scales around redshift z = 1.2. We present a method for obtaining the gas and star formation rate (SFR) surface densities of ensembles of clumps composing galaxies at this redshift, even though the corresponding scales are not resolved. This method is based on identifying these structures in position-velocity diagrams corresponding to slices within the galaxies. We use unique IRAM observations of the CO(3-2) rotational line and DEEP2 spectra of four massive star forming distant galaxies - EGS13003805, EGS13004291, EGS12007881, and EGS13019128 in the AEGIS terminology - to determine the gas and SFR surface densities of the identifiable ensembles of clumps that constitute them. The integrated CO line luminosity is assumed to be directly proportional to the total gas mass, and the SFR is deduced from the [OII] line. We identify the ensembles of clumps with the angular resolution available in both CO and [OII] spectroscopy; i.e., 1-1.5". SFR and gas surface densities are averaged in areas of this size, which is also the thickness of the DEEP2 slits and of the extracted IRAM slices, and we derive a spatially resolved Kennicutt-Schmidt (KS) relation on a scale of ~8 kpc. The data generally indicates an average depletion time of 1.9 Gyr, but with significant variations from point to point within the galaxies.

preprint2012arXiv

Analysis of the Early-Time Optical Spectra of SN 2011fe in M101

The nearby Type Ia supernova SN 2011fe in M101 (cz=241 km s^-1) provides a unique opportunity to study the early evolution of a "normal" Type Ia supernova, its compositional structure, and its elusive progenitor system. We present 18 high signal-to-noise spectra of SN 2011fe during its first month beginning 1.2 days post-explosion and with an average cadence of 1.8 days. This gives a clear picture of how various line-forming species are distributed within the outer layers of the ejecta, including that of unburned material (C+O). We follow the evolution of C II absorption features until they diminish near maximum light, showing overlapping regions of burned and unburned material between ejection velocities of 10,000 and 16,000 km s^-1. This supports the notion that incomplete burning, in addition to progenitor scenarios, is a relevant source of spectroscopic diversity among SNe Ia. The observed evolution of the highly Doppler-shifted O I 7774 absorption features detected within five days post-explosion indicate the presence of O I with expansion velocities from 11,500 to 21,000 km s^-1. The fact that some O I is present above C II suggests that SN 2011fe may have had an appreciable amount of unburned oxygen within the outer layers of the ejecta.

preprint2012arXiv

Excited OH+, H2O+, and H3O+ in NGC 4418 and Arp 220

We report on Herschel/PACS observations of absorption lines of OH+, H2O+ and H3O+ in NGC 4418 and Arp 220. Excited lines of OH+ and H2O+ with E_lower of at least 285 and \sim200 K, respectively, are detected in both sources, indicating radiative pumping and location in the high radiation density environment of the nuclear regions. Abundance ratios OH+/H2O+ of 1-2.5 are estimated in the nuclei of both sources. The inferred OH+ column and abundance relative to H nuclei are (0.5-1)x10^{16} cm-2 and \sim2x10^{-8}, respectively. Additionally, in Arp 220, an extended low excitation component around the nuclear region is found to have OH+/H2O+\sim5-10. H3O+ is detected in both sources with N(H3O+)\sim(0.5-2)x10^{16} cm-2, and in Arp 220 the pure inversion, metastable lines indicate a high rotational temperature of ~500 K, indicative of formation pumping and/or hot gas. Simple chemical models favor an ionization sequence dominated by H+ - O+ - OH+ - H2O+ - H3O+, and we also argue that the H+ production is most likely dominated by X-ray/cosmic ray ionization. The full set of observations and models leads us to propose that the molecular ions arise in a relatively low density (\gtrsim10^4 cm-3) interclump medium, in which case the ionization rate per H nucleus (including secondary ionizations) is zeta>10^{-13} s-1, a lower limit that is severalx10^2 times the highest rate estimates for Galactic regions. In Arp 220, our lower limit for zeta is compatible with estimates for the cosmic ray energy density inferred previously from the supernova rate and synchrotron radio emission, and also with the expected ionization rate produced by X-rays. In NGC 4418, we argue that X-ray ionization due to an AGN is responsible for the molecular ion production.

preprint2012arXiv

Herschel-PACS Observations of Far-IR CO Line Emission in NGC 1068: Highly Excited Molecular Gas in the Circumnuclear Disk

We report the detection of far-IR CO rotational emission from the prototypical Seyfert 2 galaxy NGC 1068. Using Herschel-PACS, we have detected 11 transitions in the J_upper=14-30 (E_upper/k_B = 580-2565 K) range, all of which are consistent with arising from within the central 10" (700 pc). The detected transitions are modeled as arising from 2 different components: a moderate excitation (ME) component close to the galaxy systemic velocity, and a high excitation (HE) component that is blueshifted by ~80 km s^{-1}. We employ a large velocity gradient (LVG) model and derive n_H2~10^{5.6} cm^{-3}, T_kin~170 K, and M_H2~10^{6.7} M_sun for the ME component, and n_H2~10^{6.4} cm^{-3}, T_kin~570 K, and M_H2~10^{5.6} M_sun for the HE component, although for both components the uncertainties in the density and mass are plus/minus (0.6-0.9) dex. We compare the CO line profiles with those of other molecular tracers observed at higher spatial and spectral resolution, and find that the ME transitions are consistent with these lines arising in the ~200 pc diameter ring of material traced by H_2 1-0 S(1) observations. The blueshift of the HE lines may also be consistent with the bluest regions of this H_2 ring, but a better kinematic match is found with a clump of infalling gas ~40 pc north of the AGN. We discuss the prospects of placing the HE component near the AGN, and conclude that while the moderate thermal pressure precludes an association with the ~1 pc radius H_2O maser disk, the HE component could potentially be located only a few parsecs more distant from the AGN, and might then provide the N_H~10^{25} cm^{-2} column obscuring the nuclear hard X-rays. Finally, we also report sensitive upper limits extending up to J_upper=50, which place constraints on a previous model prediction for the CO emission from the X-ray obscuring torus. [Abridged]

preprint2012arXiv

Herschel/PACS spectroscopy of NGC 4418 and Arp 220: H2O, H2^{18}O, OH, ^{18}OH, O I, HCN and NH3

Herschel/PACS spectroscopy of the luminous infrared galaxies NGC4418 and Arp220 reveals high excitation in H2O, OH, HCN, and NH3. In NGC4418, absorption lines were detected with E_low>800 K (H2O), 600 K (OH), 1075 K (HCN), and 600 K (NH3), while in Arp220 the excitation is somewhat lower. While outflow signatures in moderate excitation lines are seen in Arp220 as reported in previous studies, in NGC4418 the lines tracing its outer regions are redshifted relative to the nucleus, suggesting an inflow with Mdot<~12 Msun yr^{-1}. Both galaxies have warm (Tdust>~100 K) nuclear continuum components, together with a more extended component that is much more prominent and massive in Arp220. A chemical dichotomy is found in both sources: on the one hand, the nuclear regions have high H2O abundances, ~10^{-5}, and high HCN/H2O and HCN/NH3 column density ratios of 0.1-0.4 and 2-5, respectively, indicating a chemistry typical of evolved hot cores where grain mantle evaporation has occurred. On the other hand, the high OH abundance, with OH/H2O ratios of ~0.5, indicates the effects of X-rays and/or cosmic rays. The nuclear media have surface brightnesses >~10^{13} Lsun/kpc^2 and are estimated to be thick (N_H>~10^{25} cm^{-2}). While NGC4418 shows weak absorption in H2^{18}O and ^{18}OH, with a ^{16}O-to-^{18}O ratio of >~250-500, the strong absorption of the rare isotopologues in Arp220 indicates ^{16}O-to-^{18}O of 70-130. Further away from the nuclear regions, the H2O abundance decreases to <~10^{-7} and the OH/H2O ratio is 2.5-10. Despite the different scales of NGC4418, Arp220, and Mrk231, preliminary evidence is found for an evolutionary sequence from infall, hot-core like chemistry, and solar oxygen isotope ratio to high velocity outflow, disruption of the hot core chemistry and cumulative high mass stellar processing of 18O.

preprint2012arXiv

PHIBSS: molecular gas content and scaling relations in z~1-3 normal star forming galaxies

We present PHIBSS, the IRAM Plateau de Bure high-z blue sequence CO 3-2 survey of the molecular gas properties in normal star forming galaxies (SFGs) near the cosmic star formation peak. PHIBSS provides 52 CO detections in two redshift slices at z~1.2 and 2.2, with log(M*(M_solar))>10.4 and log(SFR(M_solar/yr))>1.5. Including a correction for the incomplete coverage of the M*-SFR plane, we infer average gas fractions of ~0.33 at z~1.2 and ~0.47 at z~2.2. Gas fractions drop with stellar mass, in agreement with cosmological simulations including strong star formation feedback. Most of the z~1-3 SFGs are rotationally supported turbulent disks. The sizes of CO and UV/optical emission are comparable. The molecular gas - star formation relation for the z=1-3 SFGs is near-linear, with a ~0.7 Gyrs gas depletion timescale; changes in depletion time are only a secondary effect. Since this timescale is much less than the Hubble time in all SFGs between z~0 and 2, fresh gas must be supplied with a fairly high duty cycle over several billion years. At given z and M*, gas fractions correlate strongly with the specific star formation rate. The variation of specific star formation rate between z~0 and 3 is mainly controlled by the fraction of baryonic mass that resides in cold gas.

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.

preprint2012arXiv

PTF11kx: A Type-Ia Supernova with a Symbiotic Nova Progenitor

There is a consensus that Type-Ia supernovae (SNe Ia) arise from the thermonuclear explosion of white dwarf stars that accrete matter from a binary companion. However, direct observation of SN Ia progenitors is lacking, and the precise nature of the binary companion remains uncertain. A temporal series of high-resolution optical spectra of the SN Ia PTF 11kx reveals a complex circumstellar environment that provides an unprecedentedly detailed view of the progenitor system. Multiple shells of circumsteller are detected and the SN ejecta are seen to interact with circumstellar material (CSM) starting 59 days after the explosion. These features are best described by a symbiotic nova progenitor, similar to RS Ophiuchi.

preprint2012arXiv

Spectroscopic FIR mapping of the disk and galactic wind of M82 with Herschel-PACS

[Abridged] We present maps of the main cooling lines of the neutral atomic gas ([OI] at 63 and 145 micron and [CII] at 158 micron) and in the [OIII] 88 micron line of the starburst galaxy M82, carried out with the PACS spectrometer on board the Herschel satellite. By applying PDR modeling we derive maps of the main ISM physical parameters, including the [CII] optical depth, at unprecedented spatial resolution (~300 pc). We can clearly kinematically separate the disk from the outflow in all lines. The [CII] and [OI] distributions are consistent with PDR emission both in the disk and in the outflow. Surprisingly, in the outflow, the atomic and the ionized gas traced by the [OIII] line both have a deprojected velocity of ~75 km/s, very similar to the average velocity of the outflowing cold molecular gas (~ 100 km/s) and several times smaller than the outflowing material detected in Halpha (~ 600 km/s). This suggests that the cold molecular and neutral atomic gas and the ionized gas traced by the [OIII] 88 micron line are dynamically coupled to each other but decoupled from the Halpha emitting gas. We propose a scenario where cold clouds from the disk are entrained into the outflow by the winds where they likely evaporate, surviving as small, fairly dense cloudlets (n_H\sim 500-1000 cm^-3, G_0\sim 500- 1000, T_gas\sim300 K). We show that the UV photons provided by the starburst are sufficient to excite the PDR shells around the molecular cores. The mass of the neutral atomic gas in the outflow is \gtrsim 5-12x 10^7 M_sun to be compared with that of the molecular gas (3.3 x 10^8 M_sun) and of the Halpha emitting gas (5.8 x 10^6 M_sun). The mass loading factor, (dM/dt)/SFR, of the molecular plus neutral atomic gas in the outflow is ~ 2. Energy and momentum driven outflow models can explain the data equally well, if all the outflowing gas components are taken into account.

preprint2011arXiv

Constraints on the assembly and dynamics of galaxies. II. Properties of kiloparsec-scale clumps in rest-frame optical emission of z ~ 2 star-forming galaxies

We study the properties of luminous stellar clumps identified in deep, high resolution HST/NIC2 F160W imaging at 1.6um of six z~2 star-forming galaxies with existing near-IR integral field spectroscopy from SINFONI at the VLT. Individual clumps contribute ~0.5%-15% of the galaxy-integrated rest-frame ~5000A emission, with median of about 2%; the total contribution of clump light ranges from 10%-25%. The median intrinsic clump size and stellar mass are ~1kpc and log(Mstar[Msun])~9, in the ranges for clumps identified in rest-UV or line emission in other studies. The clump sizes and masses in the subset of disks are broadly consistent with expectations for clump formation via gravitational instabilities in gas-rich, turbulent disks given the host galaxies' global properties. By combining the NIC2 data with ACS/F814W imaging available for one source, and AO-assisted SINFONI Halpha data for another, we infer modest color, M/L, and stellar age variations within each galaxy. In these two objects, sets of clumps identified at different wavelengths do not fully overlap; NIC2-identified clumps tend to be redder/older than ACS- or Halpha-identified clumps without rest-frame optical counterparts. There is evidence for a systematic trend of older ages at smaller galactocentric radii among the clumps, consistent with scenarios where inward migration of clumps transports material towards the central regions. From constraints on a bulge-like component at radii <1-3kpc, none of the five disks in our sample appears to contain a compact massive stellar core, and we do not discern a trend of bulge stellar mass fraction with stellar age of the galaxy. Further observations are necessary to probe the build-up of stellar bulges and the role of clumps in this process.

preprint2011arXiv

Dense Molecular Gas around AGN: HCN/CO in NGC3227

There is now convincing evidence that the intensity of HCN molecular line emission is enhanced around active galactic nuclei. In this paper we examine the specific case of the Seyfert galaxy NGC3227, for which there are subarcsecond resolution data for the HCN (1-0) 88 GHz and CO (2-1) 230 GHz rotational lines, enabling us to spatially separate a circumnuclear ring at a radius of 140pc and an inner nuclear region within 40pc of the AGN. The HCN(1-0)/CO(2-1) flux ratio differs by more than an order of magnitude between these two regions. We carry out large velocity gradient (LVG) computations to determine the range of parameters (gas temperature and density, HCN/CO abundance ratio, column densities and velocity gradients) that yield physically plausible solutions for the observed flux ratio in the central 100pc. The observed HCN/CO intensity ratio in the nucleus is consistent with very optically thick thermalized emission in very dense (>=10^5cm^{-3}) gas, in which case the HCN/CO abundance ratio there is unconstrained. Alternatively, the HCN/CO intensity ratio could be due to optically thinner emissions but with very high (~10^{-2}) HCN/CO abundance ratios. This possibility is more consistent with the CO and HCN emissions seen in the nuclei of the Seyfert galaxies NGC1068 and NGC6951. It would imply the velocity gradients are large and the clouds may be gravitationally unbound. We estimate that the X-ray ionisation rate at radii less than 20pc in the centre of NGC3227 exceeds 10^{-13}s^{-1}. X-ray ionisation and heating may lead to high HCN/CO ratios in warm gas in a high-ionisation molecular phase near the AGN.

preprint2011arXiv

Far-Infrared Line Deficits in Galaxies with Extreme Lfir/MH2 Ratios

We report initial results from the far-infrared fine structure line observations of a sample of 44 local starbursts, Seyfert galaxies and infrared luminous galaxies obtained with the PACS spectrometer on board Herschel. We show that the ratio between the far-infrared luminosity and the molecular gas mass, Lfir/MH2, is a much better proxy for the relative brightness of the far-infrared lines than Lfir alone. Galaxies with high Lfir/MH2 ratios tend to have weaker fine structure lines relative to their far-infrared continuum than galaxies with Lfir/MH2 < 80 Lsun/Msun. A deficit of the [C II] 158 micron line relative to Lfir was previously found with the ISO satellite, but now we show for the first time that this is a general aspect of all far-infrared fine structure lines, regardless of their origin in the ionized or neutral phase of the interstellar medium. The Lfir/MH2 value where these line deficits start to manifest is similar to the limit that separates between the two modes of star formation recently found in galaxies on the basis of studies of their gas-star formation relations. Our finding that the properties of the interstellar medium are also significantly different in these regimes provides independent support for the different star forming relations in normal disk galaxies and major merger systems. We use the spectral synthesis code Cloudy to model the emission of the lines. The expected increase of the ionization parameter with Lfir/MH2 can simultaneously explain the line deficits in the [C II], [N II] and [O I] lines.

preprint2011arXiv

How well can we Measure the Intrinsic Velocity Dispersion of Distant Disk Galaxies?

The kinematics of distant galaxies, from z=0.1 to z>2, play a key role in our understanding of galaxy evolution from early times to the present. One of the important parameters is the intrinsic, or local, velocity dispersion of a galaxy, which allows one to quantify the degree of non-circular motions such as pressure support. However, this is difficult to measure because the observed dispersion includes the effects of (often severe) beam smearing on the velocity gradient. Here we investigate four methods of measuring the dispersion that have been used in the literature, to assess their effectiveness at recovering the intrinsic dispersion. We discuss the biasses inherent in each method, and apply them to model disk galaxies in order to determine which methods yield meaningful quantities, and under what conditions. All the mean weighted dispersion estimators are affected by (residual) beam smearing. In contrast, the dispersion recovered by fitting a spatially and spectrally convolved disk model to the data is unbiassed by the beam smearing it is trying to compensate. Because of this, and because the bias it does exhibit depends only on the signal-to-noise, it can be considered reliable. However, at very low signal-to-noise, all methods should be used with caution.

preprint2011arXiv

Massive molecular outflows and negative feedback in ULIRGs observed by Herschel-PACS

Mass outflows driven by stars and active galactic nuclei are a key element in many current models of galaxy evolution. They may produce the observed black hole-galaxy mass relation and regulate and quench both star formation in the host galaxy and black hole accretion. However, observational evidence of such feedback processes through outflows of the bulk of the star forming molecular gas is still scarce. Here we report the detection of massive molecular outflows, traced by the hydroxyl molecule (OH), in far-infrared spectra of ULIRGs obtained with Herschel-PACS as part of the SHINING key project. In some of these objects the (terminal) outflow velocities exceed 1000 km/s, and their outflow rates (up to ~1200 M_sol/yr) are several times larger than their star formation rates. We compare the outflow signatures in different types of ULIRGs and in starburst galaxies to address the issue of the energy source (AGN or starburst) of these outflows. We report preliminary evidence that ULIRGs with a higher AGN luminosity (and higher AGN contribution to L_IR) have higher terminal velocities and shorter gas depletion time scales. The outflows in the observed ULIRGs are able to expel the cold gas reservoirs from the centres of these objects within ~1E6-1E8 years.

preprint2011arXiv

The metallicity dependence of the CO {\rightarrow} H_2 conversion factor in z>1 star forming galaxies

We use the first systematic samples of CO millimeter emission in z>1 'main-sequence' star forming galaxies (SFGs) to study the metallicity dependence of the conversion factor αCO, from CO line luminosity to molecular gas mass. The molecular gas depletion rate inferred from the ratio of the star formation rate (SFR) to CO luminosity, is ~1 Gyr-1 for near-solar metallicity galaxies with stellar masses above M_S~1e11 M_sun. In this regime the depletion rate does not vary more than a factor of two to three as a function of molecular gas surface density, or redshift between z~0 and 2. Below M_S the depletion rate increases rapidly with decreasing metallicity. We argue that this trend is not caused by starburst events, by changes in the physical parameters of the molecular clouds, or by the impact of the fundamental metallicity-SFR-stellar mass relation. A more probable explanation is that the conversion factor is metallicity dependent and that star formation can occur in 'CO-dark' gas. The trend is also expected theoretically from the effect of enhanced photodissociation of CO by ultraviolet radiation at low metallicity. From the available z~0 and z~1-3 samples we constrain the slope of the log(αCO) -log (metallicity) relation to range between -1 and -2, fairly insensitive to the assumed slope of the gas-star formation rate relation. Because of the lower metallicities near the peak of the galaxy formation activity at z~1-2 compared to z~0, we suggest that molecular gas masses estimated from CO luminosities have to be substantially corrected upward for galaxies below M_S.

preprint2011arXiv

The SINS survey of z~2 galaxy kinematics: properties of the giant star forming clumps

We have studied the properties of giant star forming clumps in five z~2 star-forming disks with deep SINFONI AO spectroscopy at the ESO VLT. The clumps reside in disk regions where the Toomre Q-parameter is below unity, consistent with their being bound and having formed from gravitational instability. Broad Hα/[NII] line wings demonstrate that the clumps are launching sites of powerful outflows. The inferred outflow rates are comparable to or exceed the star formation rates, in one case by a factor of eight. Typical clumps may lose a fraction of their original gas by feedback in a few hundred million years, allowing them to migrate into the center. The most active clumps may lose much of their mass and disrupt in the disk. The clumps leave a modest imprint on the gas kinematics. Velocity gradients across the clumps are 10-40 km/s/kpc, similar to the galactic rotation gradients. Given beam smearing and clump sizes, these gradients may be consistent with significant rotational support in typical clumps. Extreme clumps may not be rotationally supported; either they are not virialized, or they are predominantly pressure supported. The velocity dispersion is spatially rather constant and increases only weakly with star formation surface density. The large velocity dispersions may be driven by the release of gravitational energy, either at the outer disk/accreting streams interface, and/or by the clump migration within the disk. Spatial variations in the inferred gas phase oxygen abundance are broadly consistent with inside-out growing disks, and/or with inward migration of the clumps.

preprint2011arXiv

The star formation history of the Milky Way's Nuclear Star Cluster

We present spatially resolved imaging and integral field spectroscopy data for 450 cool giant stars within 1\,pc from Sgr\,A*. We use the prominent CO bandheads to derive effective temperatures of individual giants. Additionally we present the deepest spectroscopic observation of the Galactic Center so far, probing the number of B9/A0 main sequence stars ($2.2-2.8\,M_\odot$) in two deep fields. From spectro-photometry we construct a Hertzsprung-Russell diagram of the red giant population and fit the observed diagram with model populations to derive the star formation history of the nuclear cluster. We find that (1) the average nuclear star-formation rate dropped from an initial maximum $\sim10$\,Gyrs ago to a deep minimum 1-2\,Gyrs ago and increased again during the last few hundred Myrs, and (2) that roughly 80% of the stellar mass formed more than 5\,Gyrs ago; (3) mass estimates within $\rm R\sim1\,pc$ from Sgr\,A* favor a dominant star formation mode with a 'normal' Chabrier/Kroupa initial mass function for the majority of the past star formation in the Galactic Center. The bulk stellar mass seems to have formed under conditions significantly different from the young stellar disks, perhaps because at the time of the formation of the nuclear cluster the massive black hole and its sphere of influence was much smaller than today.

preprint2010arXiv

A Study of the Gas-Star Formation Relation over Cosmic Time

We use the first systematic data sets of CO molecular line emission in z~1-3 normal star forming galaxies for a comparison of the dependence of galaxy-averaged star formation rates on molecular gas masses at low and high redshifts, and in different galactic environments. Although the current high-z samples are still small and biased toward the luminous and massive tail of the actively star-forming 'main-sequence', a fairly clear picture is emerging. Independent of whether galaxy integrated quantities or surface densities are considered, low- and high-z SFG galaxy populations appear to follow similar molecular gas-star formation relations with slopes 1.1 to 1.2. The gas-depletion time scale in these SFGs grows from 0.5 Gyrs at z~2 to 1.5 Gyrs at z~0. Because star formation depletion times are significantly smaller than the Hubble time at all redshifts sampled, star formation rates and gas fractions are set by the balance between gas accretion from the halo and stellar feedback. In contrast, very luminous gas rich major mergers at both low-z and high-z produce on average 4 to10 times more far-infrared luminosity per unit gas mass. Only some fraction of this difference can be explained by uncertainties in gas-mass or luminosity estimators; much of it must be intrinsic. The most likely interpretation is that the star formation relation is driven by global dynamical effects. For a given mass, the more compact merger systems produce stars more rapidly because their gas clouds are more compressed with shorter dynamical times, so that they churn more quickly through the available gas reservoir than the typical normal disk galaxies. When the dependence on galactic dynamical time scale is explicitly included, disk galaxies and mergers appear to follow similar gas to star-formation relations. The mergers may be forming stars at slightly higher efficiencies than the disks.

preprint2010arXiv

Gravity-driven Lyman-alpha blobs from cold streams into galaxies

We use high-resolution cosmological hydrodynamical AMR simulations to predict the characteristics of La emission from the cold gas streams that fed galaxies in massive haloes at high redshift. The La luminosity in our simulations is powered by the release of gravitational energy as gas flows from the intergalactic medium into the halo potential wells. The UV background contributes only <20% to the gas heating. The La emissivity is due primarily to electron-impact excitation cooling radiation in gas ~2x10^4K. We calculate the La emissivities assuming collisional ionisation equilibrium (CIE) at all gas temperatures. The simulated streams are self-shielded against the UV background, so photoionisation and recombination contribute negligibly to the La line formation. We produce theoretical maps of the La surface brightnesses, assuming that ~85% of the La photons are directly observable. We find that typical haloes of mass Mv~10^12-13 Msun at z~3 emit as La blobs (LABs) with luminosities 10^43-44 erg/s. Most of the La comes from the extended narrow, partly clumpy, inflowing, cold streams that feed the growing galaxies. The predicted LAB morphology is therefore irregular, with dense clumps and elongated extensions. The linewidth is expected to range from 10^2 to more than 10^3 km/s with a large variance. The typical La surface brightness profile is proportional to r^-1.2 where r is the distance from the halo centre. Our simulated LABs are similar in luminosity, morphology and extent to the observed LABs, with distinct kinematic features. The predicted La luminosity function is consistent with observations, and the predicted areas and linewidths roughly recover the observed scaling relations. This mechanism for producing LABs appears inevitable in many high-z galaxies. Some of the LABs may thus be regarded as direct detections of the cold streams that drove galaxy evolution at high z.

preprint2010arXiv

Herschel-PACS spectroscopy of IR-bright galaxies at high redshift

We present Herschel-PACS observations of rest-frame mid-infrared and far-infrared spectral line emissions from two lensed, ultra-luminous infrared galaxies at high redshift: MIPS J142824.0+352619 (MIPS J1428), a starburst-dominated system at z = 1.3, and IRAS F10214+4724 (F10214), a source at z = 2.3 hosting both star-formation and a luminous AGN. We have detected [OI]63 micron and [OIII]52 micron in MIPS J1428, and tentatively [OIII]52 micron in F10214. Together with the recent ZEUS-CSO [CII]158 micron detection in MIPS J1428 we can for the first time combine [OI], [CII] and far-IR (FIR) continuum measurements for photo-dissociation (PDR) modeling of an ultra-luminous (L_IR > 10^12 L_sun) star forming galaxy at the peak epoch of cosmic star formation. We find that MIPS J1428, contrary to average local ULIRGs, does not show a deficit in [OI] relative to FIR. The combination of far-UV flux G_0 and gas density n (derived from the PDR models), as well as the star formation efficiency (derived from CO and FIR) is similar to normal or starburst galaxies, despite the high infrared luminosity of this system. In contrast, F10214 has stringent upper limits on [OIV] and [SIII], and an [OIII]/FIR ratio at least an order of magnitude lower than local starbursts or AGN, similar to local ULIRGs.

preprint2010arXiv

High molecular gas fractions in normal massive star forming galaxies in the young Universe

Stars form from cold molecular interstellar gas. Since this is relatively rare in the local Universe, galaxies like the Milky Way form only a few new stars per year. Typical massive galaxies in the distant Universe formed stars an order of magnitude more rapidly. Unless star formation was significantly more efficient, this difference suggests that young galaxies were much more gas rich. Molecular gas observations in the distant Universe have so far been largely restricted to very luminous, rare objects, including mergers and quasars. Here we report the results of a systematic survey of molecular gas in samples of typical massive star forming galaxies at <z>~1.2 and 2.3, when the Universe was 40% and 24% of its current age. Our measurements provide empirical evidence that distant star forming galaxies indeed were gas rich, and that the star formation efficiency is not strongly dependent on cosmic epoch. The average fraction of cold gas relative to total galaxy baryonic mass at z= 2.3 and z=1.2 is ~44% and 34%, three to ten times higher than in today's massive spiral galaxies. The slow decrease between z~2 and 1 probably requires a mechanism of semi-continuous replenishment of fresh gas to the young galaxies.

preprint2010arXiv

High-Redshift Star-Forming Galaxies: Angular Momentum and Baryon Fraction, Turbulent Pressure Effects and the Origin of Turbulence

The structure of a sample of high-redshift (z=2), rotating galaxies with high star formation rates and turbulent gas velocities of sigma=40-80 km/s is investigated. Fitting the observed disk rotational velocities and radii with a Mo, Mao, White (1998) (MMW) model requires unusually large disk spin parameters lambda_d>0.1 and disk-to-dark halo mass fraction m_d=0.2, close to the cosmic baryon fraction. The galaxies segregate into dispersion-dominated systems with 1<vmax/sigma<3, maximum rotational velocities vmax<200 km/s and disk half-light radii rd=1-3 kpc and rotation-dominated systems with vmax>200 km/s, vmax/sigma>3 and rd=4-8 kpc. For the dispersion-dominated sample, radial pressure gradients partly compensate the gravitational force, reducing the rotational velocities. Including this pressure effect in the MMW model, dispersion-dominated galaxies can be fitted well with spin parameters lf lambda_d=0.03-0.05 for high disk mass fractions of m_d=0.2 and with lambda_d=0.01-0.03 for m_d=0.05. These values are in good agreement with cosmological expectations. For the rotation-dominated sample however pressure effects are small and better agreement with theoretically expected disk spin parameters can only be achieved if the dark halo mass contribution in the visible disk regime (2-3*rd) is smaller than predicted by the MMW model. We argue that these galaxies can still be embedded in standard cold dark matter halos if the halos did not contract adiabatically in response to disk formation. It is shown that the observed high turbulent gas motions of the galaxies are consistent with a Toomre instability parameter Q=1 which is equal to the critical value, expected for gravitational disk instability to be the major driver of turbulence. The dominant energy source of turbulence is then the potential energy of the gas in the disk.

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

An Extremely Top-Heavy IMF in the Galactic Center Stellar Disks

We present new observations of the nuclear star cluster in the central parsec of the Galaxy with the adaptive optics assisted, integral field spectrograph SINFONI on the ESO/VLT. Our work allows the spectroscopic detection of early and late type stars to m_K >= 16, more than 2 magnitudes deeper than our previous data sets. Our observations result in a total sample of 177 bona fide early-type stars. We find that most of these Wolf Rayet (WR), O- and B- stars reside in two strongly warped disks between 0.8" and 12" from SgrA*, as well as a central compact concentration (the S-star cluster) centered on SgrA*. The later type B stars (m_K>15) in the radial interval between 0.8" and 12" seem to be in a more isotropic distribution outside the disks. The observed dearth of late type stars in the central few arcseconds is puzzling, even when allowing for stellar collisions. The stellar mass function of the disk stars is extremely top heavy with a best fit power law of dN/dm ~ m^(-0.45+/-0.3). Since at least the WR/O-stars were formed in situ in a single star formation event ~6 Myrs ago, this mass function probably reflects the initial mass function (IMF). The mass functions of the S-stars inside 0.8" and of the early-type stars at distances beyond 12" are compatible with a standard Salpeter/Kroupa IMF (best fit power law of dN/dm ~ m^(-2.15+/-0.3)).

preprint2009arXiv

Evidence for Warped Disks of Young Stars in the Galactic Center

The central parsec around the super-massive black hole in the Galactic Center hosts more than 100 young and massive stars. Outside the central cusp (R~1") the majority of these O and Wolf-Rayet (WR) stars reside in a main clockwise system, plus a second, less prominent disk or streamer system at large angles with respect to the main system. Here we present the results from new observations of the Galactic Center with the AO-assisted near-infrared imager NACO and the integral field spectrograph SINFONI on the ESO/VLT. These include the detection of 27 new reliably measured WR/O stars in the central 12" and improved measurements of 63 previously detected stars, with proper motion uncertainties reduced by a factor of four compared to our earlier work. We develop a detailed statistical analysis of their orbital properties and orientations. Half of the WR/O stars are compatible with being members of a clockwise rotating system. The rotation axis of this system shows a strong transition as a function of the projected distance from SgrA*. The main clockwise system either is either a strongly warped single disk with a thickness of about 10 degrees, or consists of a series of streamers with significant radial variation in their orbital planes. 11 out of 61 clockwise moving stars have an angular separation of more than 30 degrees from the clockwise system. The mean eccentricity of the clockwise system is 0.36+/-0.06. The distribution of the counter-clockwise WR/O star is not isotropic at the 98% confidence level. It is compatible with a coherent structure such as stellar filaments, streams, small clusters or possibly a disk in a dissolving state. The observed disk warp and the steep surface density distribution favor in situ star formation in gaseous accretion disks as the origin of the young stars.

preprint2009arXiv

Spitzer Quasar and ULIRG Evolution Study (QUEST). IV. Comparison of 1-Jy Ultraluminous Infrared Galaxies with Palomar-Green Quasars

We report the results from a comprehensive study of 74 ultraluminous infrared galaxies (ULIRGs) and 34 Palomar-Green (PG) quasars within z ~ 0.3$ observed with the Spitzer Infrared Spectrograph (IRS). The contribution of nuclear activity to the bolometric luminosity in these systems is quantified using six independent methods that span a range in wavelength and give consistent results within ~ +/-10-15% on average. The average derived AGN contribution in ULIRGs is ~35-40%, ranging from ~15-35% among "cool" (f_25/f_60 =< 0.2) optically classified HII-like and LINER ULIRGs to ~50 and ~75% among warm Seyfert 2 and Seyfert 1 ULIRGs, respectively. This number exceeds ~80% in PG QSOs. ULIRGs fall in one of three distinct AGN classes: (1) objects with small extinctions and large PAH equivalent widths are highly starburst-dominated; (2) systems with large extinctions and modest PAH equivalent widths have larger AGN contributions, but still tend to be starburst-dominated; and (3) ULIRGs with both small extinctions and small PAH equivalent widths host AGN that are at least as powerful as the starbursts. The AGN contributions in class 2 ULIRGs are more uncertain than in the other objects, and we cannot formally rule out the possibility that these objects represent a physically distinct type of ULIRGs. A morphological trend is seen along the sequence (1)-(2)-(3), in general agreement with the standard ULIRG - QSO evolution scenario and suggestive of a broad peak in extinction during the intermediate stages of merger evolution. However, the scatter in this sequence, implies that black hole accretion, in addition to depending on the merger phase, also has a strong chaotic/random component, as in local AGN. (abridged)

preprint2009arXiv

The halo merger rate in the Millennium Simulation and implications for observed galaxy merger fractions

We developed a new method to extract halo merger rates from the Millennium Simulation. First, by removing superfluous mergers that are artifacts of the FOF halo finder, we find a lower merger rate compared to previous work. The reductions are more significant at lower redshifts, lower halo masses, and minor mergers. Our approach agrees better with predictions from the EPS model. Second, we find that the FOF halo finder overestimates the halo mass by up to 50% for halos that are about to merge, which leads to an additional ~20% overestimate of the merger rate. Therefore, we define halo masses by including only gravitationally bound particles. We provide new best-fitting parameters for a global formula to account for these improvements. In addition, we extract the merger rate per progenitor halo, as well as per descendant halo. The former is the quantity that is related to observed galaxy merger fractions when they are measured via pair counting. At low mass/redshift the merger rate increases with mass and redshift. At high mass/redshift (for halos with masses a few times the "knee" of the mass function) these trends break down, and the merger rate per progenitor halo decreases with mass and increases only moderately with redshift. Defining the merger rate per progenitor halo also allows us to quantify the rate at which halos are being accreted onto larger halos. We provide an analytic formula to convert given merger rates per descendant halo into merger rates per progenitor halo. Finally, we perform a comparison between observed merger fractions and the fraction of halos that have undergone a major merger during the recent dynamical friction time, and find a fair agreement, within the large uncertainties of the observations. Our new halo merger trees are available at http://www.mpe.mpg.de/ir/MillenniumMergerTrees/.

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.

preprint2009arXiv

The SINS survey: modeling the dynamics of z~2 galaxies and the high-z Tully-Fisher relation

We present the modeling of SINFONI integral field dynamics of 18 star forming galaxies at z ~ 2 from Halpha line emission. The galaxies are selected from the larger sample of the SINS survey, based on the prominence of ordered rotational motions with respect to more complex merger induced dynamics. The quality of the data allows us to carefully select systems with kinematics dominated by rotation, and to model the gas dynamics across the whole galaxy using suitable exponential disk models. We obtain a good correlation between the dynamical mass and the stellar mass, finding that large gas fractions Mgas~M*) are required to explain the difference between the two quantities. We use the derived stellar mass and maximum rotational velocity Vmax from the modeling to construct for the first time the stellar mass Tully-Fisher relation at z ~ 2.2. The relation obtained shows a slope similar to what is observed at lower redshift, but we detect an evolution of the zero point. We find that at z ~ 2.2 there is an offset in log(M*) for a given rotational velocity of 0.41+-0.11 with respect to the local Universe. This result is consistent with the predictions of the latest N-body/hydrodynamical simulations of disk formation and evolution, which invoke gas accretion onto the forming disk in filaments and cooling flows. This scenario is in agreement with other dynamical evidence from SINS, where gas accretion from the halo is required to reproduce the observed properties of a large fraction of the z ~ 2 galaxies.

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.

preprint2007arXiv

A PDR-Code Comparison Study

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

preprint2006arXiv

SINFONI Integral Field Spectroscopy of z~2 UV-selected Galaxies: Rotation Curves and Dynamical Evolution

We present 0.5" resolution near-IR integral field spectroscopy of the Ha line emission of 14 z~2 UV-selected BM/BX galaxies obtained with SINFONI at ESO/VLT. The mean Ha half-light radius r_1/2 is about 4kpc and line emission is detected over > ~20kpc in several sources. In 9 sources, we detect spatially-resolved velocity gradients, from 40 to 410 km/s over ~10kpc. The observed kinematics of the larger systems are consistent with orbital motions. Four galaxies are well described by rotating disks with clumpy morphologies and we extract rotation curves out to radii > ~10kpc. One or two galaxies exhibit signatures more consistent with mergers. Analyzing all 14 galaxies in the framework of rotating disks, we infer mean inclination- and beam-corrected maximum circular velocities v_c of 180+-90 km/s and dynamical masses of (0.5-25)x10^10 Msun within r_1/2. On average, the dynamical masses are consistent with photometric stellar masses assuming a Chabrier/Kroupa IMF but too small for a 0.1-100 Msun Salpeter IMF. The specific angular momenta of our BM/BX galaxies are similar to those of local late-type galaxies. The specific angular momenta of their baryons are comparable to those of their dark matter halos. Extrapolating from the average v_c at 10kpc, the virial mass of the typical halo of a galaxy in our sample is 10^(11.7+-0.5) Msun. Kinematic modeling of the 3 best cases implies a ratio of v_c to local velocity dispersion of order 2-4 and accordingly a large geometric thickness. We argue that this suggests a mass accretion (alternatively, gas exhaustion) timescale of ~500Myr. We also argue that if our BM/BX galaxies were initially gas rich, their clumpy disks will subsequently lose their angular momentum and form compact bulges on a timescale of ~1 Gyr. [ABRIDGED]

preprint2002arXiv

Atomic Hydrogen Gas in Dark-Matter Minihalos and the Compact High Velocity Clouds

We calculate the coupled hydrostatic and ionization structures of pressure-supported gas clouds that are confined by gravitationally dominant dark-matter (DM) mini-halos and by an external bounding pressure provided by a hot medium. We focus on clouds that are photoionized and heated by the present-day background metagalactic field and determine the conditions for the formation of warm (WNM), and multi-phased (CNM/WNM) neutral atomic hydrogen (HI) cores in the DM-dominated clouds. We consider LCDM dark-matter halos, and we compute models for a wide range of halo masses, total cloud gas masses, and external bounding pressures. We present models for the pressure-supported HI structures observed in the Local Group dwarf galaxies Leo A and Sag DIG. We then construct minihalo models for the multi-phased (and low-metallicity) compact high-velocity HI clouds (CHVCs). If the CHVCs are drawn from the same family of halos that successfully reproduce the dwarf galaxy observations, then the CHVCs must be "circumgalactic objects" with characteristic distances of 150 kpc. For such systems we find that multi-phased behavior occurs for peak WNM HI column densities between 2e19 and 1e20 cm^-2, consistent with observations. If the large population of CHVCs represent "missing low-mass satellites" of the Galaxy, then these clouds must be pressure-confined to keep the gas neutral. For an implied CHVC minihalo scale velocity of v_s=12 km s^-1, the confining pressure must exceed ~50 cm^-3 K. A hot (2e6 K) Galactic corona could provide the required pressure at 150 kpc.

preprint2000arXiv

Massive star formation and evolution in starburst galaxies: mid-infrared spectroscopy with ISO-SWS

We present new ISO-SWS data for a sample of 27 starburst galaxies, and with these data examine the issues of formation and evolution of the most massive stars in starburst galaxies. Using starburst models which incorporate time evolution, new stellar atmosphere models for massive stars, and a starburst model geometry derived from observations of the prototypical starburst M82, we model the integrated mid-infrared line ratio [NeIII](15.6 microns)/[NeII](12.8 microns). This line ratio is sensitive to the hardness of the stellar energy distribution and therefore to the most massive stars present. We conclude from our models, with consideration of recent determinations of the stellar census in local, high-mass star forming regions, that the [NeIII]/[NeII] ratios we measure are consistent with the formation of massive (~50-100 solar mass) stars in most starbursts. In this framework, the low nebular excitation inferred from the measured line ratios can be attributed to aging effects. By including estimates of the ratio of infrared-to-Lyman continuum luminosity for the galaxies in our sample, we further find that most starbursts are relatively short-lived (1-10 million years), only a few O-star lifetimes. We discuss a possible cause of such short events: the effectiveness of stellar winds and supernovae in destroying the starburst environment.

preprint1998arXiv

Infrared spectroscopy of NGC 4151: ISO observations and NLR line profiles

We present ISO-SWS and ISOPHOT-S spectroscopy of the Seyfert galaxy NGC 4151. We detect a total of 17 fine structure emission lines emitted by a wide range of low- and high-excitation ions, two rotational lines of molecular hydrogen, and the Br beta HI line. We find that the mid-IR fine structure line profiles display blue asymmetries which are very similar to those observed in the optical lines produced in the narrow line region. Because the mid-infrared lines are much less sensitive to extinction than are the optical lines this similarity places strong constraints on scenarios which have been invoked to explain the optical line asymmetries. For example, we are able to rule out the simplest radial-motion-plus-dust scenarios for the production of the line asymmetries. Our preferred model is that of a central, geometrically thin but optically thick, obscuring screen of sub-arcsecond extension, enclosing a total hydrogen gas mass of >~5E6 Msun. This mass may be molecular. The weakness of `PAH' emission features in the low resolution spectrum is evidence that star formation plays a minor role in the circumnuclear region of NGC 4151. In a companion paper (Alexander et al., accepted by ApJ), we use the rich set of mid-infrared lines to determine the obscured photoionizing continuum produced by the active galactic nucleus.

preprint1997arXiv

What Powers Ultra-luminous IRAS Galaxies?

We present an ISO SWS and ISOPHOT-S, mid-infrared spectroscopic survey of 15 ultra-luminous IRAS galaxies. We combine the survey results with a detailed case study, based on near-IR and mm imaging spectroscopy, of one of the sample galaxies (UGC 5101). We compare the near- and mid-IR characteristics of these ultra-luminous galaxies to ISO and literature data of thirty starburst and active galactic nuclei (AGN), template galaxies. We find that 1) 70-80% of the ultra-luminous IRAS galaxies in our sample are predominantly powered by recently formed massive stars. 20-30% are powered by a central AGN. These conclusions are based on a new infrared 'diagnostic diagram' involving the ratio of high to low excitation mid-IR emission lines on the one hand, and on the strength of the 7.7um PAH feature on the other hand. 2) at least half of the sources probably have simultaneously an active nucleus and starburst activity in a 1-2 kpc diameter circum-nuclear disk/ring. 3) the mid-infrared emitting regions are highly obscured. After correction for these extinctions, we estimate that the star forming regions in ULIRGs have ages between 10^7 and 10^8 years, similar to but somewhat larger than those found in lower luminosity starburst galaxies. 4) in the sample we have studied there is no obvious trend for the AGN component to dominate in the most compact, and thus most advanced mergers. Instead, at any given time during the merger evolution, the time dependent compression of the circum-nuclear interstellar gas, the accretion rate onto the central black hole and the associated radiation efficiency may determine whether star formation or AGN activity dominates the luminosity of the system.