Source author record

M. Spaans

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

ResearcherUnclaimed source record

astro-ph.GA astro-ph.CO astro-ph.SR astro-ph.EP astro-ph gr-qc hep-th hep-ex physics.chem-ph physics.gen-ph

Catalog footprint

What is connected

77works

10topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2016arXiv

ALMA observations of Lyman-alpha Blob 1: halo sub-structure illuminated from within

We present new Atacama Large Millimeter/Submillimeter Array (ALMA) 850um continuum observations of the original Lyman-alpha Blob (LAB) in the SSA22 field at z=3.1 (SSA22-LAB01). The ALMA map resolves the previously identified submillimeter source into three components with total flux density S_850 = 1.68+/-0.06 mJy, corresponding to a star formation rate of ~150 M_sun/yr. The submillimeter sources are associated with several faint (m~27 mag) rest-frame ultraviolet sources identified in Hubble Space Telescope Imaging Spectrograph (STIS) clear filter imaging (~5850A). One of these companions is spectroscopically confirmed with Keck MOSFIRE to lie within 20 projected kpc and 250 km/s of one of the ALMA components. We postulate that some of these STIS sources represent a population of low-mass star-forming satellites surrounding the central submillimeter sources, potentially contributing to their growth and activity through accretion. Using a high resolution cosmological zoom simulation of a 10^13 M_sun halo at z=3, including stellar, dust and Ly-alpha radiative transfer, we can model the ALMA+STIS observations and demonstrate that Ly-alpha photons escaping from the central submillimeter sources are expected to resonantly scatter in neutral hydrogen, the majority of which is predicted to be associated with halo substructure. We show how this process gives rise to extended Ly-alpha emission with similar surface brightness and morphology to observed giant LABs.

preprint2016arXiv

Constraining cloud parameters using high density gas tracers in galaxies

Far-infrared molecular emission is an important tool used to understand the excitation mechanisms of the gas in the inter-stellar medium of star-forming galaxies. In the present work, we model the emission from rotational transitions with critical densities n >~ 10^4 cm-3. We include 4-3 < J <= 15-14 transitions of CO and 13CO, in addition to J <= 7-6 transitions of HCN, HNC, and HCO+ on galactic scales. We do this by re-sampling high density gas in a hydrodynamic model of a gas-rich disk galaxy, assuming that the density field of the interstellar medium of the model galaxy follows the probability density function (PDF) inferred from the resolved low density scales. We find that in a narrow gas density PDF, with a mean density of ~10 cm-3 and a dispersion σ= 2.1 in the log of the density, most of the emission of molecular lines, emanates from the 10-1000 cm-3 part of the PDF. We construct synthetic emission maps for the central 2 kpc of the galaxy and fit the line ratios of CO and 13CO up to J = 15-14, as well as HCN, HNC, and HCO+ up to J = 7-6, using one photo-dissociation region (PDR) model. We attribute the goodness of the one component fits for our model galaxy to the fact that the distribution of the luminosity, as a function of density, is peaked at gas densities between 10 and 1000 cm-3. We explore the impact of different log-normal density PDFs on the distribution of the line-luminosity as a function of density, and we show that it is necessary to have a broad dispersion, corresponding to Mach numbers >~ 30 in order to obtain significant emission from n > 10^4 cm-3 gas. Such Mach numbers are expected in star-forming galaxies, LIRGS, and ULIRGS. By fitting line ratios of HCN(1-0), HNC(1-0), and HCO+(1-0) for a sample of LIRGS and ULIRGS using mechanically heated PDRs, we constrain the Mach number of these galaxies to 29 < M < 77.

preprint2016arXiv

HCO+ Detection of Dust-Depleted Gas in the Inner Hole of the LkCa 15 Pre-Transitional Disk

LkCa 15 is an extensively studied star in the Taurus region known for its pre-transitional disk with a large inner cavity in dust continuum and normal gas accretion rate. The most popular hypothesis to explain the LkCa 15 data invokes one or more planets to carve out the inner cavity, while gas continues to flow across the gap from the outer disk onto the central star. We present spatially unresolved HCO+ J=4-3 observations of the LkCa 15 disk from the JCMT and model the data with the ProDiMo code. We find that: (1) HCO+ line-wings are clearly detected, certifying the presence of gas in the cavity within <50 AU of the star. (2) Reproducing the observed line-wing flux requires both a significant suppression of cavity dust (by a factor >10^4 compared to the ISM) and a substantial increase in the gas scale-height within the cavity (H_0/R_0 ~ 0.6). An ISM dust-to-gas ratio (d:g=10^-2) yields too little line-wing flux regardless of the scale-height or cavity gas geometry, while a smaller scale-height also under predicts the flux even with a reduced d:g. (3) The cavity gas mass is consistent with the surface density profile of the outer disk extended inwards to the sublimation radius (corresponding to mass M_d ~ 0.03 M_sun), and masses lower by a factor >10 appear to be ruled out.

preprint2016arXiv

Modelling Mechanical Heating in Star-Forming Galaxies: CO and 13CO Line Ratios as Sensitive Probes

We apply photo-dissociation region (PDR) molecular line emission models, that have varying degrees of enhanced mechanical heating rates, to the gaseous component of simulations of star-forming galaxies taken from the literature. Snapshots of these simulations are used to produce line emission maps for the rotational transitions of the CO molecule and its 13CO isotope up to J = 4-3. We consider two galaxy models: a small disk galaxy of solar metallicity and a lighter dwarf galaxy with 0.2 \zsun metallicity. Elevated excitation temperatures for CO(1 - 0) correlate positively with mechanical feedback, that is enhanced towards the central region of both model galaxies. The emission maps of these model galaxies are used to compute line ratios of CO and 13CO transitions. These line ratios are used as diagnostics where we attempt to match them These line ratios are used as diagnostics where we attempt to match them to mechanically heated single component (i.e. uniform density, Far-UV flux, visual extinction and velocity gradient) equilibrium PDR models. We find that PDRs ignoring mechanical feedback in the heating budget over-estimate the gas density by a factor of 100 and the far-UV flux by factors of ~10 - 1000. In contrast, PDRs that take mechanical feedback into account are able to fit all the line ratios for the central < 2 kpc of the fiducial disk galaxy quite well. The mean mechanical heating rate per H atom that we recover from the line ratio fits of this region varies between $10^{-27}$ -- $10^{-26}$~erg s$^{-1}$.

preprint2016arXiv

Submillimeter H2O and H2O+ emission in lensed ultra- and hyper-luminous infrared galaxies at z ~ 2-4

(abridged) We report rest-frame submillimeter H2O emission line observations of 11 HyLIRGs/ULIRGs at z~2-4 selected among the brightest lensed galaxies discovered in the Herschel-ATLAS. Using the IRAM NOEMA, we have detected 14 new H2O emission lines. The apparent luminosities of the H2O emission lines are $μL_{\rm{H_2O}} \sim 6-21 \times 10^8 L_\odot$, with velocity-integrated line fluxes ranging from 4-15 Jy km s$^{-1}$. We have also observed CO emission lines using EMIR on the IRAM 30m telescope in seven sources. The velocity widths for CO and H2O lines are found to be similar. With almost comparable integrated flux densities to those of the high-J CO line, H2O is found to be among the strongest molecular emitters in high-z Hy/ULIRGs. We also confirm our previously found correlation between luminosity of H2O ($L_{\rm{H_2O}}$) and infrared ($L_{\rm{IR}}$) that $L_{\rm{H_2O}} \sim L_{\rm{IR}}^{1.1-1.2}$, with our new detections. This correlation could be explained by a dominant role of far-infrared (FIR) pumping in the H2O excitation. Modelling reveals the FIR radiation fields have warm dust temperature $T_\rm{warm}$~45-75 K, H2O column density per unit velocity interval $N_{\rm{H_2O}}/ΔV \gtrsim 0.3 \times 10^{15}$ cm$^{-2}$ km$^{-1}$ s and 100 $μ$m continuum opacity $τ_{100} > 1$ (optically thick), indicating that H2O is likely to trace highly obscured warm dense gas. However, further observations of $J\geq4$ H2O lines are needed to better constrain the continuum optical depth and other physical conditions of the molecular gas and dust. We have also detected H2O+ emission in three sources. A tight correlation between $L_{\rm{H_2O}}$ and $L_{\rm{H_2O^+}}$ has been found in galaxies from low to high redshift. The velocity-integrated flux density ratio between H2O+ and H2O suggests that cosmic rays generated by strong star formation are possibly driving the H2O+ formation.

preprint2016arXiv

Surface chemistry in photodissociation regions

The presence of dust can strongly affect the chemical composition of the interstellar medium. We model the chemistry in photodissociation regions (PDRs) using both gas-phase and dust-phase chemical reactions. Our aim is to determine the chemical compositions of the interstellar medium (gas/dust/ice) in regions with distinct (molecular) gas densities that are exposed to radiation fields with different intensities. We have significantly improved the Meijerink PDR code by including 3050 new gas-phase chemical reactions and also by implementing surface chemistry. In particular, we have included 117 chemical reactions occurring on grain surfaces covering different processes, such as adsorption, thermal desorption, chemical desorption, two-body reactions, photo processes, and cosmic-ray processes on dust grains. We obtain abundances for different gas and solid species as a function of visual extinction, depending on the density and radiation field. We also analyse the rates of the formation of CO2 and H2O ices in different environments. In addition, we study how chemistry is affected by the presence/absence of ice mantles (bare dust or icy dust) and the impact of considering different desorption probabilities. The type of substrate (bare dust or icy dust) and the probability of desorption can significantly alter the chemistry occurring on grain surfaces, leading to differences of several orders of magnitude in the abundances of gas-phase species, such as CO, H2CO, and CH3OH. The type of substrate, together with the density and intensity of the radiation field, also determine the threshold extinction to form ices of CO2 and H2O. We also conclude that H2CO and CH3OH are mainly released into the gas phase of low, far-ultraviolet illuminated PDRs through chemical desorption upon two-body surface reactions, rather than through photodesorption.

preprint2016arXiv

The average submillimetre properties of Lyman-alpha Blobs at z=3

Ly-alpha blobs (LABs) offer insight into the complex interface between galaxies and their circumgalactic medium. Whilst some LABs have been found to contain luminous star-forming galaxies and active galactic nuclei that could potentially power the Ly-alpha emission, others appear not to be associated with obvious luminous galaxy counterparts. It has been speculated that LABs may be powered by cold gas streaming on to a central galaxy, providing an opportunity to directly observe the `cold accretion' mode of galaxy growth. Star-forming galaxies in LABs could be dust obscured and therefore detectable only at longer wavelengths. We stack deep SCUBA-2 observations of the SSA22 field to determine the average 850um flux density of 34 LABs. We measure S_850 = 0.6 +/- 0.2mJy for all LABs, but stacking the LABs by size indicates that only the largest third (area > 1794 kpc^2) have a mean detection, at 4.5 sigma, with S_850 = 1.4 +/- 0.3mJy. Only two LABs (1 and 18) have individual SCUBA-2 > 3.5 sigma detections at a depth of 1.1mJy/beam. We consider two possible mechanisms for powering the LABs and find that central star formation is likely to dominate the emission of Ly-alpha, with cold accretion playing a secondary role.

preprint2016arXiv

The mid-IR water and silicate relation in protoplanetary disks

Mid-IR water lines from protoplanetary disks around T Tauri stars have a detection rate of 50\%. Models have identified multiple physical properties of disks such as dust-to-gas mass ratio, dust size power law distribution, disk gas mass, disk inner radius, and disk scale height as potential explanation for the current detection rate. We search for a connection between mid-IR water line fluxes and the strength of the 10~$μ$m silicate feature. We analyse observed water line fluxes from three blends and compute the 10~$μ$m silicate feature strength from Spitzer spectra. We use a series of published models, exploring disk dust and gas properties, and the effects of different stars. The models also show that the increasing stellar luminosity enhance simultaneously the strength of this dust feature and the water lines fluxes. No correlation is found between the observed mid-IR water lines and the 10~$μ$m silicate. Our sample shows the same difference in the peak strength between amorphous and crystalline silicates that was noted in earlier studies, but our models do not support this intrinsic difference in silicate peak strength. Individual properties of our model series are not able to reproduce the most extreme observations, suggesting that more complex dust properties are required. A parametrized settling prescription is able to boost the peak strength by a factor 2 for the standard model. Water line fluxes are unrelated to the composition of the dust. The pronounced regular trends seen in the model results are washed out in the data due to the larger diversity in stellar and disk properties compared to our model series. The disks with with weaker mid-IR water line fluxes are depleted in gas or enhanced in dust in the inner 10~au. In the case of gas depleted disks, settling produces very strong 10~$μ$m silicate features, with strong peak strength.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts

We present a catalogue of nearly 3,000 submillimetre sources detected at 850um over ~5 square degrees surveyed as part of the James Clerk Maxwell Telescope (JCMT) SCUBA-2 Cosmology Legacy Survey (S2CLS). This is the largest survey of its kind at 850um, probing a meaningful cosmic volume at the peak of star formation activity and increasing the sample size of submillimetre galaxies selected at 850um by an order of magnitude. We describe the wide 850um survey component of S2CLS, which covers the key extragalactic survey fields: UKIDSS-UDS, COSMOS, Akari-NEP, Extended Groth Strip, Lockman Hole North, SSA22 and GOODS-North. The average 1-sigma depth of S2CLS is 1.2 mJy/beam, approaching the SCUBA-2 850um confusion limit, which we determine to be ~0.8 mJy/beam. We measure the single dish 850um number counts to unprecedented accuracy, reducing the Poisson errors on the differential counts to approximately 4% at S_850~3mJy. With several independent fields, we investigate field-to-field variance, finding that the number counts on 0.5-1 degree scales are generally within 50% of the S2CLS mean for S_850>3mJy, with scatter consistent with the Poisson and estimated cosmic variance uncertainties, although there is a marginal (2-sigma) density enhancement in the GOODS-North field. The observed number counts are in reasonable agreement with recent phenomenological and semi-analytic models. Finally, the large solid angle of S2CLS allows us to measure the bright-end counts: at S_850>10mJy there are approximately ten sources per square degree, and we detect the distinctive up-turn in the number counts indicative of the detection of local sources of 850um emission and strongly lensed high-redshift galaxies. Here we describe the data collection and reduction procedures and present calibrated maps and a catalogue of sources; these are made publicly available.

preprint2016arXiv

The SCUBA-2 Cosmology Legacy Survey: The EGS deep field I - Deep number counts and the redshift distribution of the recovered Cosmic Infrared Background at 450 and 850 um

We present deep observations at 450 um and 850 um in the Extended Groth Strip field taken with the SCUBA-2 camera mounted on the James Clerk Maxwell Telescope as part of the deep SCUBA-2 Cosmology Legacy Survey (S2CLS), achieving a central instrumental depth of $σ_{450}=1.2$ mJy/beam and $σ_{850}=0.2$ mJy/beam. We detect 57 sources at 450 um and 90 at 850 um with S/N > 3.5 over ~70 sq. arcmin. From these detections we derive the number counts at flux densities $S_{450}>4.0$ mJy and $S_{850}>0.9$ mJy, which represent the deepest number counts at these wavelengths derived using directly extracted sources from only blank-field observations with a single-dish telescope. Our measurements smoothly connect the gap between previous shallower blank-field single-dish observations and deep interferometric ALMA results. We estimate the contribution of our SCUBA-2 detected galaxies to the cosmic infrared background (CIB), as well as the contribution of 24 um-selected galaxies through a stacking technique, which add a total of $0.26\pm0.03$ and $0.07\pm0.01$ MJy/sr, at 450 um and 850 um, respectively. These surface brightnesses correspond to $60\pm20$ and $50\pm20$ per cent of the total CIB measurements, where the errors are dominated by those of the total CIB. Using the photometric redshifts of the 24 um-selected sample and the redshift distributions of the submillimetre galaxies, we find that the redshift distribution of the recovered CIB is different at each wavelength, with a peak at $z\sim1$ for 450 um and at $z\sim2$ for 850um, consistent with previous observations and theoretical models.

preprint2015arXiv

Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs

We propose a set of standard assumptions for the modelling of Class II and III protoplanetary disks, which includes detailed continuum radiative transfer, thermo-chemical modelling of gas and ice, and line radiative transfer from optical to cm wavelengths. We propose new standard dust opacities for disk models, we present a simplified treatment of PAHs sufficient to reproduce the PAH emission features, and we suggest using a simple treatment of dust settling. We roughly adjust parameters to obtain a model that predicts typical Class II T Tauri star continuum and line observations. We systematically study the impact of each model parameter (disk mass, disk extension and shape, dust settling, dust size and opacity, gas/dust ratio, etc.) on all continuum and line observables, in particular on the SED, mm-slope, continuum visibilities, and emission lines including [OI] 63um, high-J CO lines, (sub-)mm CO isotopologue lines, and CO fundamental ro-vibrational lines. We find that evolved dust properties (large grains) often needed to fit the SED, have important consequences for disk chemistry and heating/cooling balance, leading to stronger emission lines in general. Strong dust settling and missing disk flaring have similar effects on continuum observations, but opposite effects on far-IR gas emission lines. PAH molecules can shield the gas from stellar UV radiation because of their strong absorption and negligible scattering opacities. The observable millimetre-slope of the SED can become significantly more gentle in the case of cold disk midplanes, which we find regularly in our T Tauri models. We propose to use line observations of robust chemical tracers of the gas, such as O, CO, and H2, as additional constraints to determine some key properties of the disks, such as disk shape and mass, opacities, and the dust/gas ratio, by simultaneously fitting continuum and line observations.

preprint2015arXiv

Detection of a large fraction of atomic gas not associated with star-forming material in M17 SW

We probe the column densities and masses traced by the ionized and neutral atomic carbon with spectrally resolved maps, and compare them to the diffuse and dense molecular gas traced by [C I] and low-$J$ CO lines toward the star-forming region M17SW. We mapped a 4.1pc x 4.7pc region in the [C I] 609 m$μ$ line using the APEX telescope, as well as the CO isotopologues with the IRAM 30m telescope. We analyze the data based on velocity channel maps that are 1 km/s wide. We correlate their spatial distribution with that of the [C II] map obtained with SOFIA/GREAT. Optically thin approximations were used to estimate the column densities of [C I] and [C II] in each velocity channel. The spatial distribution of the [C I] and all CO isotopologues emission was found to be associated with that of [C II] in about 20%-80% of the mapped region, with the high correlation found in the central (15-23 km/s ) velocity channels. The excitation temperature of [C I] ranges between 40 K and 100 K in the inner molecular region of M17 SW. Column densities in 1 km/s channels between ~10$^{15}$ and ~10$^{17}$ cm$^{-2}$ were found for [C I]. Just ~20% of the velocity range (~40 km/s) that the [C II] line spans is associated with the star-forming material traced by [C I] and CO. The total gas mass estimated from the [C II] emission gives a lower limit of ~4.4x10$^3$ $M_{\odot}$. At least 64% of this mass is not associated with the star-forming material in M17SW. We also found that about 36%, 17%, and 47% of the [C II] emission is associated with the HII, HI, and H_2 regimes, respectively. Comparisons with the H41$α$ line shows an ionization region mixed with the neutral and part of the molecular gas, in agreement with the clumped structure and dynamical processes at play in M17SW. These results are also relevant to extra-galactic studies in which [C II] is often used as a tracer of star-forming material.

preprint2015arXiv

Disentangling the excitation conditions of the dense gas in M17 SW

We probe the chemical and energetic conditions in dense gas created by radiative feedback through observations of multiple CO, HCN and HCO$^+$ transitions toward the dense core of M17 SW. We used the dual band receiver GREAT on board the SOFIA airborne telescope to obtain maps of the $J=16-15$, $J=12-11$, and $J=11-10$ transitions of $^{12}$CO. We compare these maps with corresponding APEX and IRAM 30m telescope data for low- and mid-$J$ CO, HCN and HCO$^+$ emission lines, including maps of the HCN $J=8-7$ and HCO$^+$ $J=9-8$ transitions. The excitation conditions of $^{12}$CO, HCO$^+$ and HCN are estimated with a two-phase non-LTE radiative transfer model of the line spectral energy distributions (LSEDs) at four selected positions. The energy balance at these positions is also studied. We obtained extensive LSEDs for the CO, HCN and HCO$^+$ molecules toward M17 SW. The LSED shape, particularly the high-$J$ tail of the CO lines observed with SOFIA/GREAT, is distinctive for the underlying excitation conditions. The critical magnetic field criterion implies that the cold cloudlets at two positions are partially controlled by processes that create and dissipate internal motions. Supersonic but sub-Alfvénic velocities in the cold component at most selected positions indicates that internal motions are likely MHD waves. Magnetic pressure dominates thermal pressure in both gas components at all selected positions, assuming random orientation of the magnetic field. The magnetic pressure of a constant magnetic field throughout all the gas phases can support the total internal pressure of the cold components, but it cannot support the internal pressure of the warm components. If the magnetic field scales as $B \propto n^{2/3}$, then the evolution of the cold cloudlets at two selected positions, and the warm cloudlets at all selected positions, will be determined by ambipolar diffusion.

preprint2015arXiv

Exploring the molecular chemistry and excitation in obscured luminous infrared galaxies: An ALMA mm-wave spectral scan of NGC 4418

We obtained an ALMA Cycle 0 spectral scan of the dusty LIRG NGC 4418, spanning a total of 70.7 GHz in bands 3, 6, and 7. We use a combined local thermal equilibrium (LTE) and non-LTE (NLTE) fit of the spectrum in order to identify the molecular species and derive column densities and excitation temperatures. We derive molecular abundances and compare them with other Galactic and extragalactic sources by means of a principal component analysis. We detect 317 emission lines from a total of 45 molecular species, including 15 isotopic substitutions and six vibrationally excited variants. Our LTE/NLTE fit find kinetic temperatures from 20 to 350 K, and densities between 10$^5$ and 10$^7$ cm$^{-3}$. The spectrum is dominated by vibrationally excited HC$_3$N, HCN, and HNC, with vibrational temperatures from 300 to 450 K. We find high abundances of HC$_3$N, SiO, H$_2$S, and c-HCCCH and a low CH$_3$OH abundance. A principal component analysis shows that NGC 4418 and Arp 220 share very similar molecular abundances and excitation, which clearly set them apart from other Galactic and extragalactic environments. The similar molecular abundances observed towards NCG 4418 and Arp 220 are consistent with a hot gas-phase chemistry, with the relative abundances of SiO and CH$_3$OH being regulated by shocks and X-ray driven dissociation. The bright emission from vibrationally excited species confirms the presence of a compact IR source, with an effective diameter $<$5 pc and brightness temperatures $>$350 K. The molecular abundances and the vibrationally excited spectrum are consistent with a young AGN/starburst system. We suggest that NGC 4418 may be a template for a new kind of chemistry and excitation, typical of compact obscured nuclei (CON). Because of the narrow line widths and bright molecular emission, NGC 4418 is the ideal target for further studies of the chemistry in CONs.

preprint2015arXiv

How chemistry influences cloud structure, star formation, and the IMF

In the earliest phases of star-forming clouds, stable molecular species, such as CO, are important coolants in the gas phase. Depletion of these molecules on dust surfaces affects the thermal balance of molecular clouds and with that their whole evolution. For the first time, we study the effect of grain surface chemistry (GSC) on star formation and its impact on the initial mass function (IMF). We follow a contracting translucent cloud in which we treat the gas-grain chemical interplay in detail, including the process of freeze-out. We perform 3d hydrodynamical simulations under three different conditions, a pure gas-phase model, a freeze-out model, and a complete chemistry model. The models display different thermal evolution during cloud collapse. The equation of state (EOS) of the gas becomes softer with CO freeze-out and the results show that at the onset of star formation, the cloud retains its evolution history such that the number of formed stars differ (by 7%) between the three models. While the stellar mass distribution results in a different IMF when we consider pure freeze-out, with the complete treatment of the GSC, the divergence from a pure gas-phase model is minimal. We find that the impact of freeze-out is balanced by the non-thermal processes; chemical and photodesorption. We also find an average filament width of 0.12 pc ($\pm$0.03 pc), and speculate that this may be a result from the changes in the EOS caused by the gas-dust thermal coupling. We conclude that GSC plays a big role in the chemical composition of molecular clouds and that surface processes are needed to accurately interpret observations, however, that GSC does not have a significant impact as far as star formation and the IMF is concerned.

preprint2015arXiv

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

We present high resolution (0."4) IRAM PdBI and ALMA mm and submm observations of the (ultra) luminous infrared galaxies ((U)LIRGs) IRAS17208-0014, Arp220, IC860 and Zw049.057 that reveal intense line emission from vibrationally excited ($ν_2$=1) J=3-2 and 4-3 HCN. The emission is emerging from buried, compact (r<17-70 pc) nuclei that have very high implied mid-infrared surface brightness $>$$5\times 10^{13}$ L$_{\odot}$ kpc$^{-2}$. These nuclei are likely powered by accreting supermassive black holes (SMBHs) and/or hot (>200 K) extreme starbursts. Vibrational, $ν_2$=1, lines of HCN are excited by intense 14 micron mid-infrared emission and are excellent probes of the dynamics, masses, and physical conditions of (U)LIRG nuclei when H$_2$ column densities exceed $10^{24}$ cm$^{-2}$. It is clear that these lines open up a new interesting avenue to gain access to the most obscured AGNs and starbursts. Vibrationally excited HCN acts as a proxy for the absorbed mid-infrared emission from the embedded nuclei, which allows for reconstruction of the intrinsic, hotter dust SED. In contrast, we show strong evidence that the ground vibrational state ($ν$=0), J=3-2 and 4-3 rotational lines of HCN and HCO$^+$ fail to probe the highly enshrouded, compact nuclear regions owing to strong self- and continuum absorption. The HCN and HCO$^+$ line profiles are double-peaked because of the absorption and show evidence of non-circular motions - possibly in the form of in- or outflows. Detections of vibrationally excited HCN in external galaxies are so far limited to ULIRGs and early-type spiral LIRGs, and we discuss possible causes for this. We tentatively suggest that the peak of vibrationally excited HCN emission is connected to a rapid stage of nuclear growth, before the phase of strong feedback.

preprint2015arXiv

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

(Ultra) Luminous Infrared Galaxies ((U)LIRGs) are objects characterized by their extreme infrared (8-1000 $μ$m) luminosities ($L_{LIRG}>10^{11} $L$_\odot$ and $L_{ULIRG}>10^{12}$ L$_\odot$). The Herschel Comprehensive ULIRG Emission Survey (HerCULES; PI van der Werf) presents a representative flux-limited sample of 29 (U)LIRGs that spans the full luminosity range of these objects (10$^{11}\leq L_\odot \geq10^{13}$). With the \emph{Herschel Space Observatory}, we observe [CII] 157 $μ$m, [OI] 63 $μ$m, and [OI] 145 $μ$m line emission with PACS, CO J=4-3 through J=13-12, [CI] 370 $μ$m, and [CI] 609 $μ$m with SPIRE, and low-J CO transitions with ground-based telescopes. The CO ladders of the sample are separated into three classes based on their excitation level. In 13 of the galaxies, the [OI] 63 $μ$m emission line is self absorbed. Comparing the CO excitation to the IRAS 60/100 $μ$m ratio and to far infrared luminosity, we find that the CO excitation is more correlated to the far infrared colors. We present cooling budgets for the galaxies and find fine-structure line flux deficits in the [CII], [SiII], [OI], and [CI] lines in the objects with the highest far IR fluxes, but do not observe this for CO $4\leq J_{upp}\leq13$. In order to study the heating of the molecular gas, we present a combination of three diagnostic quantities to help determine the dominant heating source. Using the CO excitation, the CO J=1-0 linewidth, and the AGN contribution, we conclude that galaxies with large CO linewidths always have high-excitation CO ladders, and often low AGN contributions, suggesting that mechanical heating is important.

preprint2015arXiv

The JCMT Plane Survey: early results from the l = 30 degree field

We present early results from the JCMT Plane Survey (JPS), which has surveyed the northern inner Galactic plane between longitudes l=7 and l=63 degrees in the 850-μm continuum with SCUBA-2, as part of the James Clerk Maxwell Telescope Legacy Survey programme. Data from the l=30 degree survey region, which contains the massive star-forming regions W43 and G29.96, are analysed after approximately 40% of the observations had been completed. The pixel-to-pixel noise is found to be 19 mJy/beam, after a smooth over the beam area, and the projected equivalent noise levels in the final survey are expected to be around 10 mJy/beam. An initial extraction of compact sources was performed using the FellWalker method resulting in the detection of 1029 sources above a 5-σ surface-brightness threshold. The completeness limits in these data are estimated to be around 0.2 Jy/beam (peak flux density) and 0.8 Jy (integrated flux density) and are therefore probably already dominated by source confusion in this relatively crowded section of the survey. The flux densities of extracted compact sources are consistent with those of matching detections in the shallower ATLASGAL survey. We analyse the virial and evolutionary state of the detected clumps in the W43 star-forming complex and find that they appear younger than the Galactic-plane average.

preprint2014arXiv

[OI] disk emission in the Taurus star forming region

The structure of protoplanetary disks is thought to be linked to the temperature and chemistry of their dust and gas. Whether the disk is flat or flaring depends on the amount of radiation that it absorbs at a given radius, and on the efficiency with which this is converted into thermal energy. The understanding of these heating and cooling processes is crucial to provide a reliable disk structure for the interpretation of dust continuum emission and gas line fluxes. Especially in the upper layers of the disk, where gas and dust are thermally decoupled, the infrared line emission is strictly related to the gas heating/cooling processes. We aim to study the thermal properties of the disk in the oxygen line emission region, and to investigate the relative importance of X-ray (1-120 Angstrom) and far-UV radiation (FUV, 912-2070 Angstrom) for the heating balance there. We use [OI] 63 micron line fluxes observed in a sample of protoplanetary disks of the Taurus/Auriga star forming region and compare it to the model predictions presented in our previous work. The data were obtained with the PACS instrument on board the Herschel Space Observatory as part of the Herschel Open Time Key Program GASPS (GAS in Protoplanetary diskS), published in Howard et al. (2013). Our theoretical grid of disk models can reproduce the [OI] absolute fluxes and predict a correlation between [OI] and the sum Lx+Lfuv. The data show no correlation between the [OI] line flux and the X-ray luminosity, the FUV luminosity or their sum. The data show that the FUV or X-ray radiation has no notable impact on the region where the [OI] line is formed. This is in contrast with what is predicted from our models. Possible explanations are that the disks in Taurus are less flaring than the hydrostatic models predict, and/or that other disk structure aspects that were left unchanged in our models are important. ..abridged..

preprint2014arXiv

A submillimeter galaxy illuminating its circumgalactic medium: Ly-alpha scattering in a cold, clumpy outflow

We report the detection at 850um of the central source in SSA22-LAB1, the archetypal Lyman-alpha Blob (LAB), a 100kpc-scale radio-quiet emission-line nebula at z=3.1. The flux density of the source, $S_{850}=4.6\pm1.1$mJy implies the presence of a galaxy, or group of galaxies, with a total luminosity of $L_{\rm IR}\approx10^{12}L_\odot$. The position of an active source at the center of a ~50kpc-radius ring of linearly polarized Ly-alpha emission detected by Hayes et al. (2011) suggests that the central source is leaking Ly-alpha photons preferentially in the plane of the sky, which undergo scattering in HI clouds at large galactocentric radius. The Ly-alpha morphology around the submillimeter detection is reminiscent of biconical outflow, and the average Ly-alpha line profiles of the two `lobes' are dominated by a red peak, expected for a resonant line emerging from a medium with a bulk velocity gradient that is outflowing relative to the line center. Taken together, these observations provide compelling evidence that the central active galaxy (or galaxies) is responsible for a large fraction of the extended Ly-alpha emission and morphology. Less clear is the history of the cold gas in the circumgalactic medium being traced by Ly-alpha: is it mainly pristine material accreting into the halo that has not yet been processed through an interstellar medium (ISM), now being blown back as it encounters an outflow, or does it mainly comprise gas that has been swept-up within the ISM and expelled from the galaxy?

preprint2014arXiv

A UV flux constraint on the formation of direct collapse black holes

The ability of metal free gas to cool by molecular hydrogen in primordial halos is strongly associated with the strength of ultraviolet (UV) flux produced by the stellar populations in the first galaxies. Depending on the stellar spectrum, these UV photons can either dissociate $\rm H_{2}$ molecules directly or indirectly by photo-detachment of $\rm H^{-}$ as the latter provides the main pathway for $\rm H_{2}$ formation in the early universe. In this study, we aim to determine the critical strength of the UV flux above which the formation of molecular hydrogen remains suppressed for a sample of five distinct halos at $z>10$ by employing a higher order chemical solver and a Jeans resolution of 32 cells. We presume that such flux is emitted by PopII stars implying atmospheric temperatures of $\rm 10^{4}$~K. We performed three-dimensional cosmological simulations and varied the strength of the UV flux below the Lyman limit in units of $\rm J_{21}$. Our findings show that the value of $\rm J_{21}^{crit}$ varies from halo to halo and is sensitive to the local thermal conditions of the gas. For the simulated halos it varies from 400-700 with the exception of one halo where $\rm J_{21}^{crit} \geq 1500$. This has important implications for the formation of direct collapse black holes and their estimated population at z > 6. It reduces the number density of direct collapse black holes by almost three orders of magnitude compared to the previous estimates.

preprint2014arXiv

Diagnostics of the molecular component of PDRs with mechanical heating. II: line intensities and ratios

CO observations in active galactic nuclei and star-bursts reveal high kinetic temperatures. Those environments are thought to be very turbulent due to dynamic phenomena such as outflows and high supernova rates. We investigate the effect of mechanical heating (MH) on atomic fine-structure and molecular lines, and their ratios. We use those ratios as a diagnostic to constrain the amount of MH in an object and also study its significance on estimating the H2 mass. Equilibrium PDRs models were used to compute the thermal and chemical balance for the clouds. The equilibria were solved for numerically using the optimized version of the Leiden PDR-XDR code. Large velocity gradient calculations were done as post-processing on the output of the PDR models using RADEX. High-J CO line ratios are very sensitive to MH. Emission becomes at least one order of magnitude brighter in clouds with n~10^5~cm^-3 and a star formation rate of 1 Solar Mass per year (corresponding to a MH rate of 2 * 10^-19 erg cm^-3 s^-1). Emission of low-J CO lines is not as sensitive to MH, but they do become brighter in response to MH. Generally, for all of the lines we considered, MH increases excitation temperatures and decreases the optical depth at the line centre. Hence line ratios are also affected, strongly in some cases. Ratios involving HCN are a good diagnostic for MH, such as HCN(1-0)/CO(1-0) and HCN(1-0)/HCO^+(1-0). Both ratios increase by a factor 3 or more for a MH equivalent to > 5 percent of the surface heating, as opposed to pure PDRs. The first major conclusion is that low-J to high-J intensity ratios will yield a good estimate of the MH rate (as opposed to only low-J ratios). The second one is that the MH rate should be taken into account when determining A_V or equivalently N_H, and consequently the cloud mass. Ignoring MH will also lead to large errors in density and radiation field estimates.

preprint2014arXiv

Effects of turbulence and rotation on protostar formation as a precursor to seed black holes

Context. The seeds of the first supermassive black holes may have resulted from the direct collapse of hot primordial gas in $\gtrsim 10^4$ K haloes, forming a supermassive or quasistar as an intermediate stage. Aims. We explore the formation of a protostar resulting from the collapse of primordial gas in the presence of a strong Lyman-Werner radiation background. Particularly, we investigate the impact of turbulence and rotation on the fragmentation behaviour of the gas cloud. We accomplish this goal by varying the initial turbulent and rotational velocities. Methods. We performed 3D adaptive mesh refinement simulations with a resolution of 64 cells per Jeans length using the ENZO code, simulating the formation of a protostar up to unprecedentedly high central densities of $10^{21}$ cm$^{-3}$, and spatial scales of a few solar radii. To achieve this goal, we employed the KROME package to improve modelling of the chemical and thermal processes. Results. We find that the physical properties of the simulated gas clouds become similar on small scales, irrespective of the initial amount of turbulence and rotation. After the highest level of refinement was reached, the simulations have been evolved for an additional ~5 freefall times. A single bound clump with a radius of $2 \times 10^{-2}$ AU and a mass of ~$7 \times 10^{-2}$ M$_{\odot}$ is formed at the end of each simulation, marking the onset of protostar formation. No strong fragmentation is observed by the end of the simulations, regardless of the initial amount of turbulence or rotation, and high accretion rates of a few solar masses per year are found. Conclusions. Given such high accretion rates, a quasistar of $10^5$ M$_{\odot}$ is expected to form within $10^5$ years.

preprint2014arXiv

Future mmVLBI Research with ALMA: A European vision

Very long baseline interferometry at millimetre/submillimetre wavelengths (mmVLBI) offers the highest achievable spatial resolution at any wavelength in astronomy. The anticipated inclusion of ALMA as a phased array into a global VLBI network will bring unprecedented sensitivity and a transformational leap in capabilities for mmVLBI. Building on years of pioneering efforts in the US and Europe the ongoing ALMA Phasing Project (APP), a US-led international collaboration with MPIfR-led European contributions, is expected to deliver a beamformer and VLBI capability to ALMA by the end of 2014 (APP: Fish et al. 2013, arXiv:1309.3519). This report focuses on the future use of mmVLBI by the international users community from a European viewpoint. Firstly, it highlights the intense science interest in Europe in future mmVLBI observations as compiled from the responses to a general call to the European community for future research projects. A wide range of research is presented that includes, amongst others: - Imaging the event horizon of the black hole at the centre of the Galaxy - Testing the theory of General Relativity an/or searching for alternative theories - Studying the origin of AGN jets and jet formation - Cosmological evolution of galaxies and BHs, AGN feedback - Masers in the Milky Way (in stars and star-forming regions) - Extragalactic emission lines and astro-chemistry - Redshifted absorption lines in distant galaxies and study of the ISM and circumnuclear gas - Pulsars, neutron stars, X-ray binaries - Testing cosmology - Testing fundamental physical constants

preprint2014arXiv

High-lying OH absorption, [C II] deficits, and extreme $L_{\mathrm{FIR}}/M_{\mathrm{H2}}$ ratios in galaxies

Herschel/PACS observations of 29 local (Ultra-)Luminous Infrared Galaxies, including both starburst and AGN-dominated sources as diagnosed in the mid-infrared/optical, show that the equivalent width of the absorbing OH 65 um Pi_{3/2} J=9/2-7/2 line (W_{eq}(OH65)) with lower level energy E_{low}~300 K, is anticorrelated with the [C ii]158 um line to far-infrared luminosity ratio, and correlated with the far-infrared luminosity per unit gas mass and with the 60-to-100 um far-infrared color. While all sources are in the active L_{IR}/M_{H2}>50 Lsun/Msun mode as derived from previous CO line studies, the OH65 absorption shows a bimodal distribution with a discontinuity at L_{FIR}/M_{H2}~100 Lsun/Msun. In the most buried sources, OH65 probes material partially responsible for the silicate 9.7 um absorption. Combined with observations of the OH 71 um Pi_{1/2} J=7/2-5/2 doublet (E_{low}~415 K), radiative transfer models characterized by the equivalent dust temperature, Tdust, and the continuum optical depth at 100 um, tau_{100}, indicate that strong [C ii]158 um deficits are associated with far-IR thick (tau_{100}>~0.7, N_{H}>~10^{24} cm^{-2}), warm (Tdust>~60 K) structures where the OH 65 um absorption is produced, most likely in circumnuclear disks/tori/cocoons. With their high L_{FIR}/M_{H2} ratios and columns, the presence of these structures is expected to give rise to strong [C ii] deficits. W_{eq}(OH65) probes the fraction of infrared luminosity arising from these compact/warm environments, which is >~30-50% in sources with high W_{eq}({OH65}). Sources with high W_{eq}({OH65}) have surface densities of both L_{IR} and M_{H2} higher than inferred from the half-light (CO or UV/optical) radius, tracing coherent structures that represent the most buried/active stage of (circum)nuclear starburst-AGN co-evolution.

preprint2014arXiv

How realistic UV spectra and X-rays suppress the abundance of direct collapse black holes

Observations of high redshift quasars at $z>6$ indicate that they harbor supermassive black holes (SMBHs) of a billion solar masses. The direct collapse scenario has emerged as the most plausible way to assemble SMBHs. The nurseries for the direct collapse black holes are massive primordial halos illuminated with an intense UV flux emitted by population II (Pop II) stars. In this study, we compute the critical value of such a flux ($J_{21}^{\rm crit}$) for realistic spectra of Pop II stars through three-dimensional cosmological simulations. We derive the dependence of $J_{21}^{\rm crit}$ on the radiation spectra, on variations from halo to halo, and on the impact of X-ray ionization. Our findings show that the value of $J_{21}^{\rm crit}$ is a few times $\rm 10^4$ and only weakly depends on the adopted radiation spectra in the range between $T_{\rm rad}=2 \times 10^4-10^5$ K. For three simulated halos of a few times $\rm 10^{7}$~M$_{\odot}$, $J_{21}^{\rm crit}$ varies from $\rm 2 \times 10^4 - 5 \times 10^4$. The impact of X-ray ionization is almost negligible and within the expected scatter of $J_{21}^{\rm crit}$ for background fluxes of $J_{\rm X,21} \leq 0.1$. The computed estimates of $J_{21}^{\rm crit}$ have profound implications for the quasar abundance at $z=10$ as it lowers the number density of black holes forming through an isothermal direct collapse by a few orders of magnitude below the observed black holes density. However, the sites with moderate amounts of $\rm H_2$ cooling may still form massive objects sufficient to be compatible with observations.

preprint2014arXiv

Molecular line emission in NGC1068 imaged with ALMA: II. The chemistry of the dense molecular gas

We present a detailed analysis of ALMA Bands 7 and 9 data of CO, HCO+, HCN and CS, augmented with Plateau de Bure Interferometer (PdBI) data of the ~ 200 pc circumnuclear disk (CND) and the ~ 1.3 kpc starburst ring (SB ring) of NGC~1068, a nearby (D = 14 Mpc) Seyfert 2 barred galaxy. We aim at determining the physical characteristics of the dense gas present in the CND and whether the different line intensity ratios we find within the CND as well as between the CND and the SB ring are due to excitation effects (gas density and temperature differences) or to a different chemistry. We estimate the column densities of each species in Local Thermodynamic Equilibrium (LTE). We then compute large one-dimensional non-LTE radiative transfer grids (using RADEX) by using first only the CO transitions, and then all the available molecules in order to constrain the densities, temperatures and column densities within the CND. We finally present a preliminary set of chemical models to determine the origin of the gas. We find that in general the gas in the CND is very dense (> 10^5 cm^-3) and hot (T> 150K), with differences especially in the temperature across the CND. The AGN position has the lowest CO/HCO+, CO/HCN and CO/CS column density ratios. RADEX analyses seem to indicate that there is chemical differentiation across the CND. We also find differences between the chemistry of the SB ring and some regions of the CND; the SB ring is also much colder and less dense than the CND. Chemical modelling does not succeed in reproducing all the molecular ratios with one model per region, suggesting the presence of multi-gas phase components. The LTE, RADEX and chemical analyses all indicate that more than one gas-phase component is necessary to uniquely fit all the available molecular ratios within the CND.

preprint2014arXiv

Molecular line emission in NGC1068 imaged with ALMA. I An AGN-driven outflow in the dense molecular gas

We investigate the fueling and the feedback of star formation and nuclear activity in NGC1068, a nearby (D=14Mpc) Seyfert 2 barred galaxy, by analyzing the distribution and kinematics of the molecular gas in the disk. We have used ALMA to map the emission of a set of dense molecular gas tracers (CO(3-2), CO(6-5), HCN(4-3), HCO+(4-3) and CS(7-6)) and their underlying continuum emission in the central r ~ 2kpc of NGC1068 with spatial resolutions ~ 0.3"-0.5" (~ 20-35pc). Molecular line and dust continuum emissions are detected from a r ~ 200pc off-centered circumnuclear disk (CND), from the 2.6kpc-diameter bar region, and from the r ~ 1.3kpc starburst (SB) ring. Most of the emission in HCO+, HCN and CS stems from the CND. Molecular line ratios show dramatic order-of-magnitude changes inside the CND that are correlated with the UV/X-ray illumination by the AGN, betraying ongoing feedback. The gas kinematics from r ~ 50pc out to r ~ 400pc reveal a massive (M_mol ~ 2.7 (+0.9, -1.2) x 10^7 Msun) outflow in all molecular tracers. The tight correlation between the ionized gas outflow, the radio jet and the occurrence of outward motions in the disk suggests that the outflow is AGN-driven. The outflow rate estimated in the CND, dM/dt ~ 63 (+21, -37) Msun yr^-1, is an order of magnitude higher than the star formation rate at these radii, confirming that the outflow is AGN-driven. The power of the AGN is able to account for the estimated momentum and kinetic luminosity of the outflow. The CND mass load rate of the CND outflow implies a very short gas depletion time scale of <=1 Myr.

preprint2014arXiv

Probing the gaseous disk of T Tau N with CN 5-4 lines

We present spectrally resolved Herschel/HIFI observations of the young multiple system T Tau in atomic and molecular lines. While CO, H2O, [C II], and SO lines trace the envelope and the outflowing gas up to velocities of 33 km/s with respect to systemic, the CN 5-4 hyperfine structure lines at 566.7, 566.9 GHz show a narrow double-peaked profile centered at systemic velocity, consistent with an origin in the outer region of the compact disk of T Tau N. Disk modeling of the T Tau N disk with the thermo-chemical code ProDiMo produces CN line fluxes and profiles consistent with the observed ones and constrain the size of the gaseous disk (R_out = 110 (+10, -20) AU) and its inclination (i = 25 \pm 5 degree). The model indicates that the CN lines originate in a disk upper layer at 40-110 AU from the star, which is irradiated by the stellar UV field and heated up to temperatures of 50-700 K. With respect to previously observed CN 2-1 millimeter lines, the CN 5-4 lines appear to be less affected by envelope emission, due to their larger critical density and excitation temperature. Hence, high-J CN lines are a unique confusion-free tracer of embedded disks, such as the disk of T Tau N.

preprint2014arXiv

Star formation relations and CO SLEDs across the J-ladder and redshift

We present FIR-CO luminosity relations ($\log L_{\rm FIR} = α\log L'_{\rm CO} + β$) for the full CO rotational ladder from J=1-0 to J=13-12 for 62 local (z < 0.1) (Ultra) Luminous Infrared Galaxies (LIRGs) using data from Herschel SPIRE-FTS and ground-based telescopes. We extend our sample to high redshifts (z > 1) by including 35 (sub)-millimeter selected dusty star forming galaxies from the literature with robust CO observations. The addition of luminous starbursts at high redshifts enlarge the range of the FIR-CO luminosity relations towards the high-IR-luminosity end while also significantly increasing the small amount of mid-/high-J CO line data available prior to Herschel. This new data-set (both in terms of IR luminosity and J-ladder) reveals linear FIR-CO luminosity relations ($α\sim 1$) for J=1-0 up to J=5-4, with a nearly constant normalisation ($β\sim 2$). This is expected from the (also) linear FIR-(molecular line) relations found for the dense gas tracer lines (HCN and CS), as long as the dense gas mass fraction does not vary strongly within our (merger/starburst)-dominated sample. However from J=6-5 and up to J=13-12 we find an increasingly sub-linear slope and higher normalization constant with increasing J. We argue that these are caused by a warm (~100K) and dense ($>10^4{\rm cm^{-3}}$) gas component whose thermal state is unlikely to be maintained by star formation powered far-UV radiation fields (and thus is no longer directly tied to the star formation rate). We suggest that mechanical heating (e.g., supernova driven turbulence and shocks), and not cosmic rays, is the more likely source of energy for this component. The global CO spectral line energy distributions (SLEDs), which remain highly excited from J=6-5 up to J=13-12, are found to be a generic feature of the (U)LIRGs in our sample, and further support the presence of this gas component.

preprint2014arXiv

The formation of massive primordial stars in the presence of moderate UV backgrounds

Radiative feedback from populations II stars played a vital role in early structure formation. Particularly, photons below the Lyman limit can escape the star forming regions and produce a background ultraviolet (UV) flux which consequently may influence the pristine halos far away from the radiation sources. These photons can quench the formation of molecular hydrogen by photo-detachment of $\rm H^{-}$. In this study, we explore the impact of such UV radiation on fragmentation in massive primordial halos of a few times $\rm 10^{7}$~M${_\odot}$. To accomplish this goal, we perform high resolution cosmological simulations for two distinct halos and vary the strength of the impinging background UV field in units of $\rm J_{21}$. We further make use of sink particles to follow the evolution for 10,000 years after reaching the maximum refinement level. No vigorous fragmentation is observed in UV illuminated halos while the accretion rate changes according to the thermal properties. Our findings show that a few 100-10, 000 solar mass protostars are formed when halos are irradiated by $\rm J_{21}=10-500$ at $\rm z>10$ and suggest a strong relation between the strength of UV flux and mass of a protostar. This mode of star formation is quite different from minihalos, as higher accretion rates of about $\rm 0.01-0.1$ M$_{\odot}$/yr are observed by the end of our simulations. The resulting massive stars are the potential cradles for the formation of intermediate mass black holes at earlier cosmic times and contribute to the formation of a global X-ray background.

preprint2014arXiv

The nature of the ISM in galaxies during the star-formation activity peak of the Universe

We combine a semi-analytic model of galaxy formation, tracking atomic and molecular phases of cold gas, with a three-dimensional radiative-transfer and line tracing code to study the sub-mm emission from atomic and molecular species (CO, HCN, [CI], [CII], [OI]) in galaxies. We compare the physics that drives the formation of stars at the epoch of peak star formation (SF) in the Universe (z = 2.0) with that in local galaxies. We find that normal star-forming galaxies at high redshift have much higher CO-excitation peaks than their local counterparts and that CO cooling takes place at higher excitation levels. CO line ratios increase with redshift as a function of galaxy star-formation rate, but are well correlated with H2 surface density independent of redshift. We find an increase in the [OI]/[CII] line ratio in typical star-forming galaxies at z = 1.2 and z = 2.0 with respect to counterparts at z = 0. Our model results suggest that typical star-forming galaxies at high redshift consist of much denser and warmer star-forming clouds than their local counterparts. Galaxies belonging to the tail of the SF activity peak at z = 1.2 are already less dense and cooler than counterparts during the actual peak of SF activity (z = 2.0). We use our results to discuss how future ALMA surveys can best confront our predictions and constrain models of galaxy formation.

preprint2014arXiv

The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the rest-frame far-infrared emission of sub-millimeter galaxies

We present high-resolution (0.3'') ALMA 870um imaging of 52 sub-millimeter galaxies (SMGs) in the Ultra Deep Survey (UDS) field and investigate the size and morphology of the sub-millimeter (sub-mm) emission on 2-10kpc scales. We derive a median intrinsic angular size of FWHM=0.30$\pm$0.04'' for the 23 SMGs in the sample detected at a signal-to-noise ratio (SNR) >10. Using the photometric redshifts of the SMGs we show that this corresponds to a median physical half-light diameter of 2.4$\pm$0.2kpc. A stacking analysis of the SMGs detected at an SNR <10 shows they have sizes consistent with the 870um-bright SMGs in the sample. We compare our results to the sizes of SMGs derived from other multi-wavelength studies, and show that the rest-frame ~250um sizes of SMGs are consistent with studies of resolved 12CO (J=3-2 to 7-6) emission lines, but that sizes derived from 1.4GHz imaging appear to be approximately two times larger on average, which we attribute to cosmic ray diffusion. The rest-frame optical sizes of SMGs are around four times larger than the sub-millimeter sizes, indicating that the star formation in these galaxies is compact relative to the pre-existing stellar distribution. The size of the starburst region in SMGs is consistent with the majority of the star formation occurring in a central region, a few kpc in extent, with a median star formation rate surface density of 90$\pm$30Msol/yr/kpc$^2$, which may suggest that we are witnessing an intense period of bulge growth in these galaxies.

preprint2013arXiv

On Particle Mass Changes and GR: Space-time Topology Causes LHC Leakage

It is argued that a change in particle rest mass must involve a multiply connected space-time topology. The LHC can probe these topological effects through the particle leakage (~18%) that it experiences in particle mass changing interactions. With a probability of ~82%, this 4-space leakage causes a downward shift in the observable Higgs mass, from its physical value of 131.6 GeV to 125.2 GeV.

preprint2013arXiv

The impact of freeze-out on collapsing molecular clouds

Atoms and molecules, and in particular CO, are important coolants during the evolution of interstellar star-forming gas clouds. The presence of dust grains, which allow many chemical reactions to occur on their surfaces, strongly impacts the chemical composition of a cloud. At low temperatures, dust grains can lock-up species from the gas phase which freeze out and form ices. In this sense, dust can deplete important coolants. Our aim is to understand the effects of freeze-out on the thermal balance and the evolution of a gravitationally bound molecular cloud. For this purpose, we perform 3D hydrodynamical simulations with the adaptive mesh code FLASH. We simulate a gravitationally unstable cloud under two different conditions, with and without grain surface chemistry. We let the cloud evolve until one free-fall time is reached and track the thermal evolution and the abundances of species during this time. We see that at a number density of 10$^4$ cm$^{-3}$ most of the CO molecules are frozen on dust grains in the run with grain surface chemistry, thereby depriving the most important coolant. As a consequence, we find that the temperature of the gas rises up to $\sim$25 K. The temperature drops once again due to gas-grain collisional cooling when the density reaches a few$\times$10$^4$ cm$^{-3}$. We conclude that grain surface chemistry not only affects the chemical abundances in the gas phase, but also leaves a distinct imprint in the thermal evolution that impacts the fragmentation of a star-forming cloud. As a final step, we present the equation of state of a collapsing molecular cloud that has grain surface chemistry included.

preprint2013arXiv

The influence of magnetic fields, turbulence, and UV radiation on the formation of supermassive black holes

Context. The seeds of the supermassive black holes (SMBHs) with masses of ~10^9 M_Sun observed already at z ~ 6 may have formed through the direct collapse of primordial gas in T_vir >~ 10^4 K halos, whereby the gas must stay hot (~10^4 K) in order to avoid fragmentation. Aims. The interplay between magnetic fields, turbulence, and a UV radiation background during the gravitational collapse of primordial gas in a halo is explored; in particular, the possibilities for avoiding fragmentation are examined. Methods. Using an analytical one-zone model, the evolution of a cloud of primordial gas is followed from its initial cosmic expansion through turnaround, virialization, and collapse up to a density of 10^7 cm-3. Results. It was found that in halos with no significant turbulence, the critical UV background intensity (J_21^crit) for keeping the gas hot is lower by a factor ~10 for an initial comoving magnetic field B_0 ~ 2 nG than for the zero-field case, and even lower for stronger fields. In turbulent halos, J_21^crit is found to be a factor ~10 lower than for the zero-field-zero-turbulence case, and the stronger the turbulence (more massive halo and/or stronger turbulent heating) the lower J_21^crit. Conclusions. The reduction in J_21^crit is particularly important, since it exponentially increases the number of halos exposed to a supercritical radiation background.

preprint2013arXiv

The Response of Metal Rich Gas to X-Ray Irradiation from a Massive Black Hole at High Redshift: Proof of Concept

Observational studies show that there is a strong link between the formation and evolution of galaxies and the growth of supermassive black holes (SMBH) at their centers. However, the underlying physics of this observed relation is poorly understood. In order to study the effects of X-ray radiation on the surroundings of the black hole, we implement X-ray Dominated Region (XDR) physics into Enzo and use the radiation transport module Moray to calculate the radiative transfer for a polychromatic spectrum. In this work, we investigate the effects of X-ray irradiation, produced by a central massive black hole (MBH) with a mass of M = 5x10^4 M_(solar), on ambient gas with solar and zero metallicity. We find that in the solar metallicity case, due to high opacity of the metals, the energy deposition rate in the central region (< 20 pc) is high and hence the temperatures in the center are on the order of 10^(5-7) K. Moreover, due to the cooling ability and high intrinsic opacity of solar metallicity gas, column densities of 10^(24) cm^(-2) are reached at a radius of 20 pc from the MBH. These column densities are about 3 orders of magnitudes higher than in the zero metallicity case. Furthermore, in the zero metallicity case an X-ray induced H II region is formed already after 5.8 Myr. This causes a significant outflow of gas (~8x10^6 M_(solar) from the central region, with the gas reaching outflow velocities up to ~100 km s^(-1). At later times, ~23 Myr after we insert the MBH, we find that the solar metallicity case also develops an X-ray induced H II region, but delayed by ~17 Myr.

preprint2013arXiv

The SCUBA-2 Cosmology Legacy Survey: demographics of the 450 μm-population

We investigate the multi-wavelength properties of a sample of 450-μm selected sources from the SCUBA-2 Cosmology Legacy Survey (S2CLS). A total of 69 sources were identified above 4σ in deep SCUBA-2 450-μm observations overlapping the UDS and COSMOS fields and covering 210 sq. arcmin to a typical depth of σ450=1.5 mJy. Reliable cross identification are found for 58 sources (84 per cent) in Spitzer and Hubble Space Telescope WFC3/IR data. The photometric redshift distribution (dN/dz) of 450μm-selected sources is presented, showing a broad peak in the redshift range 1<z<3, and a median of z=1.4. Combining the SCUBA-2 photometry with Herschel SPIRE data from HerMES, the submm spectral energy distribution (SED) is examined via the use of modified blackbody fits, yielding aggregate values for the IR luminosity, dust temperature and emissivity of <LIR>=10^12 +/- 0.8 L_sol, <T_D>=42 +/- 11 K and <β_D>=1.6 +/- 0.5, respectively. The relationship between these SED parameters and the physical properties of galaxies is investigated, revealing correlations between T_D and LIR and between β_D and both stellar mass and effective radius. The connection between star formation rate and stellar mass is explored, with 24 per cent of 450 μm sources found to be ``star-bursts'', i.e. displaying anomalously high specific SFRs. However, both the number density and observed properties of these ``star-burst'' galaxies are found consistent with the population of normal star-forming galaxies.

preprint2013arXiv

Uncertainties in water chemistry in disks: An application to TW Hya

Context. This paper discusses the sensitivity of water lines to chemical processes and radiative transfer for the protoplanetary disk around TW Hya. The study focuses on the Herschel spectral range in the context of new line detections with the PACS instrument from the Gas in Protoplanetary Systems project (GASPS). Aims. The paper presents an overview of the chemistry in the main water reservoirs in the disk around TW Hya. It discusses the limitations in the interpretation of observed water line fluxes. Methods. ... (abbreviated) Results. We report new line detections of p-H2O (3_22-2_11) at 89.99 micron and CO J=18-17 at 144.78 micron for the disk around TW Hya. Disk modeling shows that the far-IR fine structure lines ([OI], [CII]) and molecular submm lines are very robust to uncertainties in the chemistry, while the water line fluxes can change by factors of a few. The water lines are optically thick, sub-thermally excited and can couple to the background continuum radiation field. The low-excitation water lines are also sensitive to uncertainties in the collision rates, e.g. with neutral hydrogen. The gas temperature plays an important role for the [OI] fine structure line fluxes, the water line fluxes originating from the inner disk as well as the high excitation CO, CH+ and OH lines. Conclusions. Due to their sensitivity on chemical input data and radiative transfer, water lines have to be used cautiously for understanding details of the disk structure. Water lines covering a wide range of excitation energies provide access to the various gas phase water reservoirs (inside and outside the snow line) in protoplanetary disks and thus provide important information on where gas-phase water is potentially located. Experimental and/or theoretical collision rates for H2O with atomic hydrogen are needed to diminish uncertainties from water line radiative transfer.

preprint2013arXiv

Water vapor in the protoplanetary disk of DG Tau

Water is key in the evolution of protoplanetary disks and the formation of comets and icy/water planets. While high excitation water lines originating in the hot inner disk have been detected in several T Tauri stars (TTSs), water vapor from the outer disk, where most of water ice reservoir is stored, was only reported in the closeby TTS TW Hya. We present spectrally resolved Herschel/HIFI observations of the young TTS DG Tau in the ortho- and para- water ground-state transitions at 557, 1113 GHz. The lines show a narrow double-peaked profile, consistent with an origin in the outer disk, and are ~19-26 times brighter than in TW Hya. In contrast, CO and [C II] lines are dominated by emission from the envelope/outflow, which makes H2O lines a unique tracer of the disk of DG Tau. Disk modeling with the thermo-chemical code ProDiMo indicates that the strong UV field, due to the young age and strong accretion of DG Tau, irradiates a disk upper layer at 10-90 AU from the star, heating it up to temperatures of 600 K and producing the observed bright water lines. The models suggest a disk mass of 0.015-0.1 Msun, consistent with the estimated minimum mass of the solar nebula before planet formation, and a water reservoir of ~1e2-1e3 Earth oceans in vapour, and ~100 times larger in the form of ice. Hence, this detection supports the scenario of ocean delivery on terrestrial planets by impact of icy bodies forming in the outer disk.

preprint2012arXiv

Diagnostics of the Molecular Component of PDRs with Mechanical Heating

Context. Multitransition CO observations of galaxy centers have revealed that significant fractions of the dense circumnuclear gas have high kinetic temperatures, which are hard to explain by pure photon excitation, but may be caused by dissipation of turbulent energy. Aims. We aim to determine to what extent mechanical heating should be taken into account while modelling PDRs. To this end, the effect of dissipated turbulence on the thermal and chemical properties of PDRs is explored. Methods. Clouds are modelled as 1D semi-infinite slabs whose thermal and chemical equilibrium is solved for using the Leiden PDR-XDR code. Results. In a steady-state treatment, mechanical heating seems to play an important role in determining the kinetic temperature of the gas in molecular clouds. Particularly in high-energy environments such as starburst galaxies and galaxy centers, model gas temperatures are underestimated by at least a factor of two if mechanical heating is ignored. The models also show that CO, HCN and H2 O column densities increase as a function of mechanical heating. The HNC/HCN integrated column density ratio shows a decrease by a factor of at least two in high density regions with n \sim 105 cm-3, whereas that of HCN/HCO+ shows a strong dependence on mechanical heating for this same density range, with boosts of up to three orders of magnitude. Conclusions. The effects of mechanical heating cannot be ignored in studies of the molecular gas excitation whenever the ratio of the star formation rate to the gas density is close to, or exceeds, 7 \times 10-6 M yr-1 cm4.5 . If mechanical heating is not included, predicted column densities are underestimated, sometimes even by a few orders of magnitude. As a lower bound to its importance, we determined that it has non-negligible effects already when mechanical heating is as little as 1% of the UV heating in a PDR.

preprint2012arXiv

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

We present Herschel SPIRE FTS spectroscopy of the nearby luminous infrared galaxy NGC 6240. In total 20 lines are detected, including CO J=4-3 through J=13-12, 6 H2O rotational lines, and [CI] and [NII] fine-structure lines. The CO to continuum luminosity ratio is 10 times higher in NGC 6240 than Mrk 231. Although the CO ladders of NGC 6240 and Mrk 231 are very similar, UV and/or X-ray irradiation are unlikely to be responsible for the excitation of the gas in NGC 6240. We applied both C and J shock models to the H2 v=1-0 S(1) and v=2-1 S(1) lines and the CO rotational ladder. The CO ladder is best reproduced by a model with shock velocity v_s=10 km s^-1 and a pre-shock density n_H=5 * 10^4 cm^-3. We find that the solution best fitting the H2 lines is degenerate: The shock velocities and number densities range between v_s = 17 - 47 km s^-1 and n_H=10^7 - 5 * 10^4 cm^-3, respectively. The H2 lines thus need a much more powerful shock than the CO lines. We deduce that most of the gas is currently moderately stirred up by slow (10 km s^-1) shocks while only a small fraction (< 1 percent) of the ISM is exposed to the high velocity shocks. This implies that the gas is rapidly loosing its highly turbulent motions. We argue that a high CO line-to-continuum ratio is a key diagnostic for the presence of shocks.

preprint2012arXiv

Extended warm and dense gas towards W49A: starburst conditions in our Galaxy?

The star formation rates in starburst galaxies are orders of magnitude higher than in local star-forming regions, and the origin of this difference is not well understood. We use sub-mm spectral line maps to characterize the physical conditions of the molecular gas in the luminous Galactic star-forming region W49A and compare them with the conditions in starburst galaxies. We probe the temperature and density structure of W49A using H_2CO and HCN line ratios over a 2'x2' (6.6x6.6 pc) field with an angular resolution of 15" (~0.8 pc) provided by the JCMT Spectral Legacy Survey. We analyze the rotation diagrams of lines with multiple transitions with corrections for optical depth and beam dilution, and estimate excitation temperatures and column densities. Comparing the observed line intensity ratios with non-LTE radiative transfer models, our results reveal an extended region (about 1'x1', equivalent to ~3x3 pc at the distance of W49A) of warm (> 100 K) and dense (>10^5 cm^-3) molecular gas, with a mass of 2x10^4 - 2x10^5 M_Sun (by applying abundances derived for other regions of massive star-formation). These temperatures and densities in W49A are comparable to those found in clouds near the center of the Milky Way and in starburst galaxies. The highly excited gas is likely to be heated via shocks from the stellar winds of embedded, O-type stars or alternatively due to UV irradiation, or possibly a combination of these two processes. Cosmic rays, X-ray irradiation and gas-grain collisional heating are less likely to be the source of the heating in the case of W49A.

preprint2012arXiv

FUV and X-ray irradiated protoplanetary disks: a grid of models I. The disk structure

Context. Planets are thought to eventually form from the mostly gaseous (~99% of the mass) disks around young stars. The density structure and chemical composition of protoplanetary disks are affected by the incident radiation field at optical, FUV, and X-ray wavelengths, as well as by the dust properties. Aims. The effect of FUV and X-rays on the disk structure and the gas chemical composition are investigated. This work forms the basis of a second paper, which discusses the impact on diagnostic lines of, e.g., C+, O, H2O, and Ne+ observed with facilities such as Spitzer and Herschel. Methods. A grid of 240 models is computed in which the X-ray and FUV luminosity, minimum grain size, dust size distribution, and surface density distribution are varied in a systematic way. The hydrostatic structure and the thermo-chemical structure are calculated using ProDiMo. Results. The abundance structure of neutral oxygen is stable to changes in the X-ray and FUV luminosity, and the emission lines will thus be useful tracers of the disk mass and temperature. The C+ abundance distribution is sensitive to both X-rays and FUV. The radial column density profile shows two peaks, one at the inner rim and a second one at a radius r=5-10 AU. Ne+ and other heavy elements have a very strong response to X-rays, and the column density in the inner disk increases by two orders of magnitude from the lowest (LX = 1e29 erg/s) to the highest considered X-ray flux (LX = 1e32 erg/s). FUV confines the Ne+ ionized region to areas closer to the star at low X-ray luminosities (LX = 1e29 erg/s). H2O abundances are enhanced by X-rays due to higher temperatures in the inner disk and higher ionization fractions in the outer disk. The line fluxes and profiles are affected by the effects on these species, thus providing diagnostic value in the study of FUV and X-ray irradiated disks around T Tauri stars. (abridged)

preprint2012arXiv

Hydrogenation of PAH cations: a first step towards H2 formation

Molecular hydrogen is the most abundant molecule in the universe. A large fraction of H2 forms by association of hydrogen atoms adsorbed on polycyclic aromatic hydrocarbons (PAHs), where formation rates depend crucially on the H sticking probability. We have experimentally studied PAH hydrogenation by exposing coronene cations, confined in a radiofrequency ion trap, to gas phase atomic hydrogen. A systematic increase of the number of H atoms adsorbed on the coronene with the time of exposure is observed. Odd coronene hydrogenation states dominate the mass spectrum up to 11 H atoms attached. This indicates the presence of a barrier preventing H attachment to these molecular systems. For the second and fourth hydrogenation, barrier heights of 72 +- 6 meV and 40 +- 10 meV, respectively are found which is in good agreement with theoretical predictions for the hydrogenation of neutral PAHs. Our experiments however prove that the barrier does not vanish for higher hydrogenation states. These results imply that PAH cations, as their neutral counterparts, exist in highly hydrogenated forms in the interstellar medium. Due to this catalytic activity, PAH cations and neutrals seem to contribute similarly to the formation of H2.

preprint2012arXiv

Testing the Proposed Connection between Dark Energy and Black Holes

In 1997, an extension of general relativity was proposed that predicts the dark energy density Λ to vary linearly with the total number of macroscopic black holes in the universe. We explore this prediction and find that Λ must be roughly constant after the bulk of the stellar mass black holes are in place, so for a redshift z smaller than unity. Conversely, the highest black hole formation rate corresponds to the peak in the cosmic star formation history, earlier than z=1. This yields a fast declining Λ, by a factor of about 5, from redshift 1 to 3. At even earlier times, before many stars were formed, the value of Λ should be much smaller than its current value. These predicted effects are all consistent with current data, and near future observations can definitively confirm or disproof the link between the dark energy density and the total number of black holes in the universe.

preprint2012arXiv

The impact of magnetic fields on the IMF in star-forming clouds near a supermassive black hole

Star formation in the centers of galaxies is thought to yield massive stars with a possibly top-heavy stellar mass distribution. It is likely that magnetic fields play a crucial role in the distribution of stellar masses inside star-forming molecular clouds. In this context, we explore the effects of magnetic fields, with a typical field strength of 38 microG, such as in RCW 38, and a field strength of 135 microG, similar to NGC 2024 and the infrared dark cloud G28.34+0.06, on the initial mass function (IMF) near (< 10 pc) a 10^7 solar mass black hole. Using these conditions, we perform a series of numerical simulations with the hydrodynamical code FLASH to elucidate the impact of magnetic fields on the IMF and the star-formation efficiency (SFE) emerging from an 800 solar mass cloud. We find that the collapse of a gravitationally unstable molecular cloud is slowed down with increasing magnetic field strength and that stars form along the field lines. The total number of stars formed during the simulations increases by a factor of 1.5-2 with magnetic fields. The main component of the IMF has a lognormal shape, with its peak shifted to sub-solar (< 0.3 M_sun) masses in the presence of magnetic fields, due to a decrease in the accretion rates from the gas reservoir. In addition, we see a top-heavy, nearly flat IMF above ~2 solar masses, from regions that were supported by magnetic pressure until high masses are reached. We also consider the effects of X-ray irradiation if the central black hole is active. X-ray feedback inhibits the formation of sub-solar masses and decreases the SFEs even further. Thus, the second contribution is no longer visible. We conclude that magnetic fields potentially change the SFE and the IMF both in active and inactive galaxies, and need to be taken into account in such calculations.

preprint2012arXiv

The implications of dust for high-redshift protogalaxies and the formation of binary disks

Numerical simulations suggest that the first galaxies are formed in protogalactic halos with virial temperatures >= 10^4 K. It is likely that such halos are polluted with trace amounts of metals produced by the first generation of stars. The presence of dust can significantly change the chemistry and dynamics of early galaxies. In this article, we aim to assess the role of dust on the thermal and dynamical evolution of the first galaxies in the presence of a background UV flux, and its implications for the observability of Lyman alpha emitters and sub-mm sources. We have performed high resolution cosmological simulations using the adaptive mesh refinement code FLASH to accomplish this goal. We have developed a chemical network appropriate for these conditions and coupled it with the FLASH code. The main ingredients of our chemical model include the formation of molecules, a multi-level treatment of atomic hydrogen, line trapping of Lyman alpha photons and, photoionization and photodissociation processes in a UV background. We found that the formation of molecules (H_{2} and HD) is significantly enhanced in the presence of dust grains. The presence of a background UV flux strongly influences the formation of molecules by photodissociating them. We explore the evolution after a major merger, leading to the formation of a binary disk. These disks have gas masses of ~10^{7} M_sun at a redshift of 5.4. Each disk lies in a separate subhalo as a result of the merger event. The disks are supported by turbulent pressure due to the highly supersonic turbulence present in the halo. The presence of dust does not significantly reduce the Lyman alpha emission. The emission of Lyman alpha is extended and originates from the envelope of the halo due to line trapping effects. We also find that dust masses of 10^8 M_sun are required to observe the dust continuum emission from z 5 galaxies with ALMA.

preprint2012arXiv

The SCUBA-2 Cosmology Legacy Survey: blank-field number counts of 450um-selected galaxies and their contribution to the cosmic infrared background

The first deep blank-field 450um map (1-sigma~1.3mJy) from the SCUBA-2 Cosmology Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is presented. Our map covers 140 arcmin^2 of the COSMOS field, in the footprint of the HST CANDELS area. Using 60 submillimetre galaxies (SMGs) detected at >3.75-sigma, we evaluate the number counts of 450um-selected galaxies with flux densities S_450>5mJy. The 8-arcsec JCMT beam and high sensitivity of SCUBA-2 now make it possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking at ~200um) into the individual galaxies responsible for its emission than has previously been possible at this wavelength. At S_450>5mJy we resolve (7.4[+/-]0.7)x10^-2 MJy/sr of the CIB at 450um (equivalent to 16[+/-]7% of the absolute brightness measured by COBE at this wavelength) into point sources. A further ~40% of the CIB can be recovered through a statistical stack of 24um emitters in this field, indicating that the majority (~60%) of the CIB at 450um is emitted by galaxies with S_450>2mJy. The average redshift of 450um emitters identified with an optical/near-infrared counterpart is estimated to be <z>=1.3, implying that the galaxies in the sample are in the ultraluminous class (L_IR~1.1x10^12 L_sun). If the galaxies contributing to the statistical stack lie at similar redshifts, then the majority of the CIB at 450um is emitted by galaxies in the LIRG class with L_IR>3.6x10^11 L_sun.

preprint2011arXiv

Diffuse interstellar bands in Upper Scorpius: Probing variations in the DIB spectrum due to changing environmental conditions

We study the effects of local environmental conditions affecting the diffuse interstellar band (DIB) carriers within the Upper Scorpius subgroup of the Sco OB2 association. We obtained optical spectra with FEROS, and measured the equivalent widths of five DIBs (at 5780, 5797, 6196, 6379, and 6613 Angstrom) as well as those of absorption lines of di-atomic molecules (CH, CH+, CN) and atoms (KI, CaI) towards 89 targets in the direction of Upper Scorpius. We construct a simple radiative transfer and chemical network model of the diffuse interstellar medium (ISM) sheet in front of Upp Sco to infer the effective radiation field. Both the interstellar radiation field strength, I_UV, and the molecular hydrogen fraction, f_H2, have been derived for 55 sightlines probing the Upp Sco ISM. We discuss the relations between DIB strengths, CH and CH+ line strengths, E(B-V), I_UV, and f_H2. The ratio between the 5780 and 5797 DIBs reveals a (spatial) dependence on the local environment in terms of cloud density and exposure to the interstellar radiation field, reflecting the molecular nature of these DIB carriers.

preprint2011arXiv

Enhanced H2O formation through dust grain chemistry in X-ray exposed environments

The ULIRG Mrk 231 exhibits very strong water rotational lines between λ= 200-670μm, comparable to the strength of the CO rotational lines. High redshift quasars also show similar CO and H2O line properties, while starburst galaxies, such as M82, lack these very strong H2O lines in the same wavelength range, but do show strong CO lines. We explore the possibility of enhancing the gas phase H2O abundance in X-ray exposed environments, using bare interstellar carbonaceous dust grains as a catalyst. Cloud-cloud collisions cause C and J shocks, and strip the grains of their ice layers. The internal UV field created by X-rays from the accreting black hole does not allow to reform the ice. We determine formation rates of both OH and H2O on dust grains, having temperature T_dust=10-60 K, using both Monte Carlo as well as rate equation method simulations. The acquired formation rates are added to our X-ray chemistry code, that allows us to calculate the thermal and chemical structure of the interstellar medium near an active galactic nucleus. We derive analytic expressions for the formation of OH and H2O on bare dust grains as a catalyst. Oxygen atoms arriving on the dust are released into the gas phase under the form of OH and H2O. The efficiencies of this conversion due to the chemistry occurring on dust are of order 30 percent for oxygen converted into OH and 60 percent for oxygen converted into H_2O between T_dust=15-40 K. At higher temperatures, the efficiencies rapidly decline. When the gas is mostly atomic, molecule formation on dust is dominant over the gas-phase route, which is then quenched by the low H2 abundance. Here, it is possible to enhance the warm (T> 200 K) water abundance by an order of magnitude in X-ray exposed environments. This helps to explain the observed bright water lines in nearby and high-redshift ULIRGs and Quasars.

preprint2011arXiv

Gas modelling in the disc of HD 163296

We present detailed model fits to observations of the disc around the Herbig Ae star HD 163296. This well-studied object has an age of ~ 4 Myr, with evidence of a circumstellar disc extending out to ~ 540AU. We use the radiation thermo-chemical disc code ProDiMo to model the gas and dust in the circumstellar disc of HD 163296, and attempt to determine the disc properties by fitting to observational line and continuum data. These include new Herschel/PACS observations obtained as part of the open-time key program GASPS (Gas in Protoplanetary Systems), consisting of a detection of the [OI]63mic line and upper limits for several other far infrared lines. We complement this with continuum data and ground-based observations of the 12CO 3-2, 2-1 and 13CO J=1-0 line transitions, as well as the H2 S(1) transition. We explore the effects of stellar ultraviolet variability and dust settling on the line emission, and on the derived disc properties. Our fitting efforts lead to derived gas/dust ratios in the range 9-100, depending on the assumptions made. We note that the line fluxes are sensitive in general to the degree of dust settling in the disc, with an increase in line flux for settled models. This is most pronounced in lines which are formed in the warm gas in the inner disc, but the low excitation molecular lines are also affected. This has serious implications for attempts to derive the disc gas mass from line observations. We derive fractional PAH abundances between 0.007 and 0.04 relative to ISM levels. Using a stellar and UV excess input spectrum based on a detailed analysis of observations, we find that the all observations are consistent with the previously assumed disc geometry.

preprint2011arXiv

Interstellar ices as witnesses of star formation: selective deuteration of water and organic molecules unveiled

Observations of star forming environments revealed that the abundances of some deuterated interstellar molecules are markedly larger than the cosmic D/H ratio of 10-5. Possible reasons for this pointed to grain surface chemistry. How- ever, organic molecules and water, which are both ice constituents, do not enjoy the same deuteration. For example, deuterated formaldehyde is very abundant in comets and star forming regions, while deuterated water rarely is. In this article, we explain this selective deuteration by following the formation of ices (using the rate equation method) in translucent clouds, as well as their evolu- tion as the cloud collapses to form a star. Ices start with the deposition of gas phase CO and O onto dust grains. While reaction of oxygen with atoms (H or D) or molecules (H2) yields H2O (HDO), CO only reacts with atoms (H and D) to form H2CO (HDCO, D2CO). As a result, the deuteration of formaldehyde is sensitive to the gas D/H ratio as the cloud undergoes gravitational collapse, while the deuteration of water strongly depends on the dust temperature at the time of ice formation. These results reproduce well the deuterium fractionation of formaldehyde observed in comets and star forming regions and can explain the wide spread of deuterium fractionation of water observed in these environments.

preprint2011arXiv

Lyman alpha emission from the first galaxies: Implications of UV backgrounds and the formation of molecules

The Lyman alpha line is a robust tracer of high redshift galaxies. We present estimates of Lyman alpha emission from a protogalactic halo illuminated by UV background radiation fields with various intensities. For this purpose, we performed cosmological hydrodynamics simulations with the adaptive mesh refinement code FLASH, including a detailed network for primordial chemistry,comprising the formation of primordial molecules, a multi-level model for the hydrogen atom as well as the photo-ionization and photo-dissociation processes in a UV background. We find that the presence of a background radiation field J_21 excites the emission of Lyman alpha photons, increasing the Lyman alpha luminosity up to two orders of magnitude. For a halo of \sim 10^10 M_sun, we find that a maximum flux of 5 \times 10^-15 erg cm^-2 s^-1 is obtained for J21 \times f_esc = 0.1, where f_esc is the escape fraction of the ionizing radiation. Depending on the environmental conditions, the flux may vary by three orders of magnitude. For J_21 \times f_esc > 0.1 the Lyman alpha luminosity decreases as the atomic hydrogen abundance becomes rather small. The fluxes derived here can be probed using Subaru and the upcoming James Webb Space Telescope. The emission of Lyman alpha photons is extended and comes from the envelope of the halo rather than its core. In the center of the halo, line trapping becomes effective above columns of 10^22 cm^-2 and suppresses the emission of Lyman alpha. In addition, cooling by primordial molecules may decrease the gas temperature in the central region, which further reduces Lyman alpha emission. In the central core, H_2 is photo-dissociated for a background flux of J_21 \geq 1000. For weaker radiation fields, i.e. J_21 < 0.1, H_2 and HD cooling are particularly strong in the center of the halo, leading to gas temperatures as low as \sim 100 K.

preprint2011arXiv

Lyman alpha emission from the first galaxies: Signatures of accretion and infall in the presence of line trapping

The formation of the first galaxies is accompanied by large accretion flows and virialization shocks, during which the gas is shock-heated to temperatures of $\sim10^4$ K, leading to potentially strong fluxes in the Lyman alpha line. Indeed, a number of Lyman alpha blobs has been detected at high redshift. In this letter, we explore the origin of such Lyman alpha emission using cosmological hydrodynamical simulations that include a detailed model of atomic hydrogen as a multi-level atom and the effects of line trapping with the adaptive mesh refinement code FLASH. We see that baryons fall into the center of a halo through cold streams of gas, giving rise to a Lyman alpha luminosity of at least $\rm 10^{44} erg s^{-1}$ at $\rm z=4.7$, similar to observed Lyman alpha blobs. We find that a Lyman alpha flux of $\rm 5.0\times 10^{-17} erg cm^{-2} s^{-1}$ emerges from the envelope of the halo rather than its center, where the photons are efficiently trapped. Such emission can be probed in detail with the upcoming James Webb Space Telescope (JWST) and will constitute an important probe of gas infall and accretion.

preprint2011arXiv

Star formation near an obscured AGN: Variations in the initial mass function

The conditions that affect the formation of stars in radiatively and mechanically active environments are quite different than the conditions that apply to our local interstellar neighborhood. In such galactic environments, a variety of feedback processes can play a significant role in shaping the initial mass function (IMF). Here, we present a numerical study on the effects of an accreting black hole and the influence of nearby massive stars to a collapsing, 800 M_sun, molecular cloud at 10 pc distance from the black hole. We parametrize and study radiative feedback effects of hard X-rays emanating from the black hole broad line region, increased cosmic ray rates due to supernovae in starbursts, and strong UV radiation produced by nearby massive stars. We also investigate the importance of shear from the supermassive, 10^6-10^8 M_sun, black hole as the star-forming cloud orbits around it. We find that thermal pressure from X-rays compresses the cloud, which induces a high star formation rate early on, but reduces the overall star formation efficiency to about 7% due to gas depletion by evaporation. We see that the turn-over mass of the IMF increases up to a factor of 2.3, M_turn = 1-1.5 M_sun, for the model with the highest X-ray flux (160 erg s^-1 cm^-2), while the high-mass slope of the IMF becomes Gamma > -1. This results in more high mass stars and a non-Salpeter IMF. Cosmic rays penetrate deeply into the cloud and increase the gas temperature (50-200 K), which leads to a reduced formation efficiency of low mass stars. High cosmic ray rates increase the average mass of stars, thereby shifting the turn-over mass to higher values, i.e., up to several solar masses. Due to this process, the onset of star formation is also delayed. We conclude that the IMF inside active galaxies is different than the one obtained from local environments.

preprint2011arXiv

The deeply obscured AGN of NGC4945 I. Spitzer-IRS maps of [Ne V], [Ne II], H2 0-0 S(1), S(2), and other tracers

The nearly edge-on galaxy NGC4945 is one of the closest galaxies where an AGN and starburst coexist, and is one of the brightest sources at 100 keV. Near and mid-infrared spectroscopy have shown very strong obscuration of its central region, rivaled only in strength by some of the most deeply obscured ULIRGs. We aim to determine the spatial distribution of ISM features in the central 426x426 pc^2 of NGC4945. We map the central region of NGC4945 in three Spitzer-IRS modules (SH, SL and LL). We produce maps of the flux distribution of the starburst tracers [Ne II], [Ne III], [S III] and [S IV] at 12.81, 15.56, 18.71 and 10.51 mum, respectively, and a map of the AGN narrow-line region tracer [Ne V] at 14.32 mum. We also mapped the S(1), S(2) and S(3) pure rotational lines of H2, which trace the distribution of warm molecular hydrogen. We obtained an extinction map (A_V) based on the apparent strength of the 9.7 mum silicate absorption feature. Our A_V map traces the contours of the starburst ring but the highest extinction (A_V(9.85 mum)~60) is found at the nucleus. Within the uncertainty of the astrometry all emission lines are found to peak on the nucleus, except for the warm molecular hydrogen emission which shows a maximum 60-100 pc NW of the nucleus. We favour a scenario in which the lower H2 0-0 S(1) and S(2) rotational lines originate mainly from an unobscured extra-nuclear component associated with the super-wind cone observed in the HST NICMOS map of the H2 1-0 S(1) vibrational line. For the [Ne V] emission we infer an attenuation of a factor 12-160 (A_V=55-112) based on a comparison of the ratio of our [Ne V] flux and the absorption-corrected 14-195 keV Swift-BAT flux to the average [Ne V]/BAT ratio for Seyfert 1 nuclei. The high attenuation indicates that [Ne V] and [O IV] cannot be used as extinction-free tracers of AGN power in galaxies with deeply buried nuclei.

preprint2011arXiv

The JCMT Spectral Legacy Survey: physical structure of the molecular envelope of the high-mass protostar AFGL2591

The understanding of the formation process of massive stars (>8 Msun) is limited, due to theoretical complications and observational challenges. We investigate the physical structure of the large-scale (~10^4-10^5 AU) molecular envelope of the high-mass protostar AFGL2591 using spectral imaging in the 330-373 GHz regime from the JCMT Spectral Legacy Survey. Out of ~160 spectral features, this paper uses the 35 that are spatially resolved. The observed spatial distributions of a selection of six species are compared with radiative transfer models based on a static spherically symmetric structure, a dynamic spherical structure, and a static flattened structure. The maps of CO and its isotopic variations exhibit elongated geometries on scales of ~100", and smaller scale substructure is found in maps of N2H+, o-H2CO, CS, SO2, CCH, and methanol lines. A velocity gradient is apparent in maps of all molecular lines presented here, except SO, SO2, and H2CO. We find two emission peaks in warm (Eup~200K) methanol separated by 12", indicative of a secondary heating source in the envelope. The spherical models are able to explain the distribution of emission for the optically thin H13CO+ and C34S, but not for the optically thick HCN, HCO+, and CS, nor for the optically thin C17O. The introduction of velocity structure mitigates the optical depth effects, but does not fully explain the observations, especially in the spectral dimension. A static flattened envelope viewed at a small inclination angle does slightly better. We conclude that a geometry of the envelope other than an isotropic static sphere is needed to circumvent line optical depth effects. We propose that this could be achieved in envelope models with an outflow cavity and/or inhomogeneous structure at scales smaller than ~10^4 AU. The picture of inhomogeneity is supported by observed substructure in at least six species.

preprint2011arXiv

When sticking influences H2 formation

Aims. Interstellar dust grains, because of their catalytic properties, are crucial to the formation of H2, the most abundant molecule in the Universe. The formation of molecular hydrogen strongly depends on the ability of H atoms to stick on dust grains. In this study we determine the sticking coefficient of H atoms chemisorbed on graphitic surfaces, and estimate its impact on the forma- tion of H2. Methods. The sticking probability of H atoms chemisorbed onto graphitic surfaces is obtained using a mixed classical-quantum dynamics method. In this, the H atom is treated quantum- mechanically and the vibrational modes of the surface are treated classically. The implications of sticking for the formation of H2 are addressed by using Kinetic Monte Carlo simulations that follow how atoms stick, move and associate with each other on dust surfaces of different temper- ature. Results. In our model, molecular hydrogen forms very efficiently for dust temperatures lower than 15 K through the involvement of physisorbed H atoms. At dust temperatures higher than 15 K and gas temperatures lower than 2000 K, H2 formation differs strongly if the H atoms com- ing from the gas phase have to cross a square barrier (usually considered in previous studies) or a barrier obtained by DFT calculations to become chemisorbed. The product of sticking times efficiency can be increased by many orders of magnitude when realistic barriers are considered. If graphite phonons are taken into account in the dynamics calculations, then H atoms stick better on the surface at high energies, but the overall H2 formation efficiency is only slightly affected. Our results suggest that H2 formation can proceed efficiently in photon dominated regions, X-ray dominated regions, hot cores and in the early Universe when the first dust is available.

preprint2010arXiv

A new condition for the transition from runaway to oligarchic growth

Accretion among macroscopic bodies of ~km size or larger is enhanced significantly due to gravitational focusing. Two regimes can be distinguished. Initially, the system experiences runaway growth, in which the gravitational focusing factors increase, and bodies at the high-mass tail of the distribution grow fastest. However, at some point the runaway body dynamically heats its environment, gravitational focusing factors decrease, and runaway growth passes into oligarchic growth. Based on the results of recent simulations, we reconsider the runaway growth-oligarchy transition. In contrast to oligarchy, we find that runaway growth cannot be approximated with a two component model (of small and large bodies) and that the criterion of Ida & Makino (1993), which is frequently adopted as the start of oligarchy, is not a sufficient condition to signify the transition. Instead, we propose a new criterion based on timescale arguments. We then find a larger value for the runaway growth-oligarchy transition: from several hundreds of km in the inner disk regions up to ~10^3 km for the outer disk. These findings are consistent with the view that runaway growth has been responsible for the size distribution of the present day Kuiper belt objects. Our finding furthermore outlines the proper initial conditions at the start of the oligarchy stage.

preprint2010arXiv

Accretion among preplanetary bodies: the many faces of runaway growth

(abridged) When preplanetary bodies reach proportions of ~1 km or larger in size, their accretion rate is enhanced due to gravitational focusing (GF). We have developed a new numerical model to calculate the collisional evolution of the gravitationally-enhanced growth stage. We validate our approach against existing N-body and statistical codes. Using the numerical model, we explore the characteristics of the runaway growth and the oligarchic growth accretion phases starting from an initial population of single planetesimal radius R_0. In models where the initial random velocity dispersion (as derived from their eccentricity) starts out below the escape speed of the planetesimal bodies, the system experiences runaway growth. We find that during the runaway growth phase the size distribution remains continuous but evolves into a power-law at the high mass end, consistent with previous studies. Furthermore, we find that the largest body accretes from all mass bins; a simple two component approximation is inapplicable during this stage. However, with growth the runaway body stirs up the random motions of the planetesimal population from which it is accreting. Ultimately, this feedback stops the fast growth and the system passes into oligarchy, where competitor bodies from neighboring zones catch up in terms of mass. Compared to previous estimates, we find that the system leaves the runaway growth phase at a somewhat larger radius. Furthermore, we assess the relevance of small, single-size fragments on the growth process. In classical models, where the initial velocity dispersion of bodies is small, these do not play a critical role during the runaway growth; however, in models that are characterized by large initial relative velocities due to external stirring of their random motions, a situation can emerge where fragments dominate the accretion.

preprint2010arXiv

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

We present a full high resolution SPIRE FTS spectrum of the nearby ultraluminous infrared galaxy Mrk231. In total 25 lines are detected, including CO J=5-4 through J=13-12, 7 rotational lines of H2O, 3 of OH+ and one line each of H2O+, CH+, and HF. We find that the excitation of the CO rotational levels up to J=8 can be accounted for by UV radiation from star formation. However, the approximately flat luminosity distribution of the CO lines over the rotational ladder above J=8 requires the presence of a separate source of excitation for the highest CO lines. We explore X-ray heating by the accreting supermassive black hole in Mrk231 as a source of excitation for these lines, and find that it can reproduce the observed luminosities. We also consider a model with dense gas in a strong UV radiation field to produce the highest CO lines, but find that this model strongly overpredicts the hot dust mass in Mrk231. Our favoured model consists of a star forming disk of radius 560 pc, containing clumps of dense gas exposed to strong UV radiation, dominating the emission of CO lines up to J=8. X-rays from the accreting supermassive black hole in Mrk231 dominate the excitation and chemistry of the inner disk out to a radius of 160 pc, consistent with the X-ray power of the AGN in Mrk231. The extraordinary luminosity of the OH+ and H2O+ lines reveals the signature of X-ray driven excitation and chemistry in this region.

preprint2010arXiv

Evolution of the ISM in Luminous IR Galaxies

Molecules that trace the high-density regions of the interstellar medium may be used to evaluate the changing physical and chemical environment during the ongoing nuclear activity in (Ultra-)Luminous Infrared Galaxies. The changing ratios of the HCN(1-0), HNC(1-0), HCO+(1-0), CN(1-0) and CN(2-1), and CS(3-2) transitions were compared with the HCN(1-0)/CO(1-0) ratio, which is proposed to represent the progression time scale of the starburst. These diagnostic diagrams were interpreted using the results of theoretical modeling using a large physical and chemical network to describe the state of the nuclear ISM in the evolving galaxies. Systematic changes are seen in the line ratios as the sources evolve from early stage for the nuclear starburst (ULIRGs) to later stages. These changes result from changing environmental conditions and particularly from the lowering of the average density of the medium. A temperature rise due to mechanical heating of the medium by feedback explains the lowering of the ratios at later evolutionary stages. Infrared pumping may affect the CN and HNC line ratios during early evolutionary stages. Molecular transitions display a behavior that relates to changes of the environment during an evolving nuclear starburst. Molecular properties may be used to designate the evolutionary stage of the nuclear starburst. The HCN(1-0)/CO(1-0) and HCO+(1-0)/HCN(1-0) ratios serve as indicators of the time evolution of the outburst.

preprint2010arXiv

Gravity and Information: Putting a Bit of Quantum into GR

It is shown that quantum aspects of the energy-momentum tensor reveal dark matter and dark energy behavior for general relativity. Special meaning is assigned to the operator X^2d^2/dX^2 +2Xd/dX, for any field X that is part of the energy-momentum tensor, and to the distinction between information on gravitational and non-gravitational effects. Changes in such information cause the Sun-Earth distance to increase by 11 cm/yr.

preprint2010arXiv

Herschel observations of water vapour in Markarian 231

The Ultra Luminous InfraRed Galaxy Mrk 231 reveals up to seven rotational lines of water (H2O) in emission, including a very high-lying (E_{upper}=640 K) line detected at a 4sigma level, within the Herschel/SPIRE wavelength range, whereas PACS observations show one H2O line at 78 microns in absorption, as found for other H2O lines previously detected by ISO. The absorption/emission dichotomy is caused by the pumping of the rotational levels by far-infrared radiation emitted by dust, and subsequent relaxation through lines at longer wavelengths, which allows us to estimate both the column density of H2O and the general characteristics of the underlying far-infrared continuum source. Radiative transfer models including excitation through both absorption of far-infrared radiation emitted by dust and collisions are used to calculate the equilibrium level populations of H2O and the corresponding line fluxes. The highest-lying H2O lines detected in emission, with levels at 300-640 K above the ground state, indicate that the source of far-infrared radiation responsible for the pumping is compact (radius=110-180 pc) and warm (T_{dust}=85-95 K), accounting for at least 45% of the bolometric luminosity. The high column density, N(H2O)~5x10^{17} cm^{-2}, found in this nuclear component, is most probably the consequence of shocks/cosmic rays, an XDR chemistry, and/or an "undepleted chemistry" where grain mantles are evaporated. A more extended region, presumably the inner region of the 1-kpc disk observed in other molecular species, could contribute to the flux observed in low-lying H2O lines through dense hot cores, and/or shocks. The H2O 78 micron line observed with PACS shows hints of a blue-shifted wing seen in absorption, possibly indicating the occurrence of H2O in the prominent outflow detected in OH (Fischer et al., this volume).

preprint2010arXiv

Star Formation in Extreme Environments: The Effects of Cosmic Rays and Mechanical Heating

Context: Molecular data of extreme environments, such as Arp 220, but also NGC 253, show evidence for extremely high cosmic ray (CR) rates (10^3-10^4 * Milky Way) and mechanical heating from supernova driven turbulence. Aims: The consequences of high CR rates and mechanical heating on the chemistry in clouds are explored. Methods: PDR model predictions are made for low, n=10^3, and high, n=10^5.5 cm^-3, density clouds using well-tested chemistry and radiation transfer codes. Column densities of relevant species are discussed, and special attention is given to water related species. Fluxes are shown for fine-structure lines of O, C+, C, and N+, and molecular lines of CO, HCN, HNC, and HCO+. A comparison is made to an X-ray dominated region model. Results: Fine-structure lines of [CII], [CI], and [OI] are remarkably similar for different mechanical heating and CR rates, when already exposed to large amounts of UV. HCN and H2O abundances are boosted for very high mechanical heating rates, while ionized species are relatively unaffected. OH+ and H2O+ are enhanced for very high CR rates zeta > 5 * 10^-14 s^-1. A combination of OH+, OH, H2O+, H2O, and H3O+ trace the CR rates, and are able to distinguish between enhanced cosmic rays and X-rays.

preprint2010arXiv

The Impact of Lyman alpha Trapping on the Formation of Primordial Objects

Numerous cosmological simulations have been performed to study the formation of the first objects. We present the results of high resolution 3-D cosmological simulations of primordial objects formation using the adaptive mesh refinement code FLASH by including in an approximate manner the radiative transfer effects of Lyman alpha photons. We compare the results of a Lyman alpha trapping case inside gas clouds with atomic and molecular hydrogen cooling cases.The principal objective of this research is to follow the collapse of a zero metallicity halo with an effective equation of state (that accounts for the trapping) and to explore the fate of a halo in each of the three cases, specifically, the impact of thermodynamics on fragmentation of halos.Our results show that in the case of Lyman alpha trapping, fragmentation is halted and a massive object is formed at the center of a halo. The temperature of the gas remains well above $10^{4}$ K and the halo is not able to fragment to stellar masses. In the atomic cooling case, gas collapses into one or two massive clumps in contrast to the Lyman alpha trapping case. For the molecular hydrogen cooling case, gas cools efficiently and fragments.The formation of massive primordial objects is thus strongly dependent on the thermodynamics of the gas. A salient feature of our results is that for the formation of massive objects, e.g. intermediate mass black holes, feedback effects are not required to suppress $H_{2}$ cooling, as molecular hydrogen is collisionally dissociated at temperatures higher than $10^{4}$ K as a consequence of Lyman alpha trapping.

preprint2010arXiv

The impact of X-rays on molecular cloud fragmentation and the IMF

Star formation is regulated through a variety of feedback processes. In this study, we treat feedback by X-rays and discuss its implications. Our aim is to investigate whether star formation is significantly affected when a star forming cloud resides in the vicinity of a strong X-ray source. We perform an Eulerian grid simulation with embedded Lagrangian sink particles of a collapsing molecular cloud near a massive, 10^7 M_o black hole. The chemical and thermal changes caused by radiation are incorporated into the FLASH code. When there is strong X-ray feedback the star forming cloud fragments into larger clumps whereby fewer but more massive protostellar cores are formed. Competitive accretion has a strong impact on the mass function and a near-flat, non-Salpeter IMF results.

preprint2010arXiv

X-ray impact on the protoplanetary disks around T Tauri stars

Context: T Tauri stars have X-ray luminosities ranging from L_X = 10^28-10^32 erg/s. These luminosities are similar to UV luminosities (L_UV 10^30-10^31 erg/s) and therefore X-rays are expected to affect the physics and chemistry of the upper layers of their surrounding protoplanetary disks. Aim: The effects and importance of X-rays on the chemical and hydrostatic structure of protoplanetary disks are investigated, species tracing X-ray irradiation (for L_X >= 10^29 erg/s) are identified and predictions for [OI], [CII] and [NII] fine structure line fluxes are provided. Methods: We have implemented X-ray physics and chemistry into the chemo-physical disk code ProDiMo. We include Coulomb heating and H2 ionization as heating processes and primary and secondary ionization due to X-rays in the chemistry. Results: X-rays heat up the gas causing it to expand in the optically thin surface layers. Neutral molecular species are not much affected in their abundance and spatial distribution, but charged species such as N+, OH+, H2O+ and H3O+ show enhanced abundances in the disk surface. Conclusions: Coulomb heating by X-rays changes the vertical structure of the disk, yielding temperatures of ~ 8000 K out to distances of 50 AU. The chemical structure is altered by the high electron abundance in the gas in the disk surface, causing an efficient ion-molecule chemistry. The products of this, OH+, H2O+ and H3O+, are of great interest for observations of low-mass young stellar objects with the Herschel Space Observatory. [OI] (at 63 and 145 mic) and [CII] (at 158 mic) fine structure emission are only affected for L_X > 10^30 erg/s.

preprint2009arXiv

CHAMP+ observations of warm gas in M17 SW

Since the main cooling lines of the gas phase are important tracers of the interstellar medium in Galactic and extragalactic sources, proper and detailed understanding of their emission, and the ambient conditions of the emitting gas, is necessary for a robust interpretation of the observations. With high resolution (7"-9") maps (~3x3 pc^2) of mid-J molecular lines we aim to probe the physical conditions and spatial distribution of the warm (50 to few hundred K) and dense gas (n(H_2)>10^5 cm^-3) across the interface region of M17 SW nebula. We have used the dual color multiple pixel receiver CHAMP+ on APEX telescope to obtain a 5'.3x4'.7 map of the J=6-5 and J=7-6 transitions of 12CO, the 13CO J=6-5 line, and the {^3P_2}-{^3P_1} 370 um fine-structure transition of [C I] in M17 SW. LTE and non-LTE radiative transfer models are used to constrain the ambient conditions. The warm gas extends up to a distance of ~2.2 pc from the M17 SW ridge. The 13CO J=6-5 and [C I] 370 um lines have a narrower spatial extent of about 1.3 pc along a strip line at P.A=63 deg. The structure and distribution of the [C I] {^3P_2}-{^3P_1} 370 um map indicate that its emission arises from the interclump medium with densities of the order of 10^3 cm^-3. The warmest gas is located along the ridge of the cloud, close to the ionization front. An LTE approximation indicates that the excitation temperature of the embedded clumps goes up to ~120 K. The non-LTE model suggests that the kinetic temperature at four selected positions cannot exceed 230 K in clumps of density n(H_2)~5x10^5 cm^-3, and that the warm T_k>100 K and dense (n(H_2)>10^4 cm^-3) gas traced by the mid-J 12CO lines represent just about 2% of the bulk of the molecular gas. The clump volume filling factor ranges between 0.04 and 0.11 at these positions.

preprint2009arXiv

Cooling lines as probes of the formation and buildup of galaxies and black holes

We discuss the use of SPICA to study the cosmic history of star formation and accretion by supermassive black holes. The cooling lines, in particular the high-J rotational lines of CO, provide a clear-cut and unique diagnostic for separating the contributions of star formation and AGN accretion to the total infrared luminosity of active, gas-rich galaxies. We briefly review existing efforts for studying high-J CO emission from galaxies at low and high redshift. We finally comment on the detectability of cooling radiation from primordial (very low metallicity) galaxies containing an accreting supermassive black hole with SPICA/SAFARI.

preprint2009arXiv

Probing X-ray irradiation in the nucleus of NGC 1068 with observations of high-J lines of dense gas tracers

With the incorporation of high-J molecular lines, we aim to constrain the physical conditions of the dense gas in the central region of the Seyfert 2 galaxy NGC 1068 and to determine signatures of the AGN or the starburst contribution. We used the James Clerk Maxwell Telescope to observe the J=4-3 transition of HCN, HNC, and HCO+, as well as the CN N_J=2_{5/2}-1_{3/2} and N_J=3_{5/2}-2_{5/2}, in NGC 1068. We estimate the excitation conditions of HCN, HNC, and CN, based on the line intensity ratios and radiative transfer models. We find that the bulk emission of HCN, HNC, CN, and the high-J HCO+ emerge from dense gas n(H_2)>=10^5 cm^-3). However, the low-J HCO+ lines (dominating the HCO+ column density) trace less dense (n(H_2)<10^5 cm^-3) and colder (T_K<=20 K) gas, whereas the high-J HCO+ emerges from warmer (>30 K) gas than the other molecules. The HCO+ J=4-3 line intensity, compared with the lower transition lines and with the HCN J=4-3 line, support the influence of a local XDR environment. The estimated N(CN)/N(HCN)~1-4 column density ratios are indicative of an XDR/AGN environment with a possible contribution of grain-surface chemistry induced by X-rays or shocks.

preprint2009arXiv

The Thermodynamics of Molecular Cloud Fragmentation: Star Formation Under Non-Milky Way Conditions

Properties of candidate stars, forming out of molecular clouds, depend on the ambient conditions of the parent cloud. We present a series of 2D and 3D simulations of fragmentation of molecular clouds in starburst regions as well as clouds under conditions in dwarf galaxies, leading to the formation of protostellar cores. We explore in particular the metallicity dependence of molecular cloud fragmentation and the possible variations in the dense core mass function, as the expression of a multi-phase ISM, due to dynamic and thermodynamic effects in starburst and metal-poor environments. To study the level of fragmentation during the collapse, the adaptive mesh refinement code FLASH is used. With this code, including self-gravity, thermal balance, turbulence and shocks, collapse is simulated with four different metallicity dependent cooling functions. Turbulent and rotational energies are considered as well. During the simulations, number densities of 10^8-10^9 cm^-3 are reached. The influences of dust and cosmic ray heating are investigated and a comparison to isothermal cases is made. The presented results indicate that fragmentation increases with metallicity, while cosmic ray and gas-grain collisional heating counteract this. We also find that modest rotation and turbulence can impact the cloud evolution as far as fragmentation is concerned. In this light, we conclude that radiative feedback, in starburst regions, will inhibit fragmentation, while low-metallicity dwarfs also should enjoy a star formation mode in which fragmentation is suppressed.

preprint2009arXiv

Water and ammonia abundances in S140 with the Odin satellite

We have used the Odin satellite to obtain strip maps of the ground-state rotational transitions of ortho-water and ortho-ammonia, as well as CO(5-4) and 13CO(5-4) across the PDR, and H218O in the central position. A physi-chemical inhomogeneous PDR model was used to compute the temperature and abundance distributions for water, ammonia and CO. A multi-zone escape probability method then calculated the level populations and intensity distributions. These results are compared to a homogeneous model computed with an enhanced version of the RADEX code. H2O, NH3 and 13CO show emission from an extended PDR with a narrow line width of ~3 kms. Like CO, the water line profile is dominated by outflow emission, however, mainly in the red wing. The PDR model suggests that the water emission mainly arises from the surfaces of optically thick, high density clumps with n(H2)>10^6 cm^-3 and a clump water abundance, with respect to H2, of 5x10^-8. The mean water abundance in the PDR is 5x10^-9, and between ~2x10^-8 -- 2x10^-7 in the outflow derived from a simple two-level approximation. Ammonia is also observed in the extended clumpy PDR, likely from the same high density and warm clumps as water. The average ammonia abundance is about the same as for water: 4x10^-9 and 8x10^-9 given by the PDR model and RADEX, respectively. The similarity of water and ammonia PDR emission is also seen in the almost identical line profiles observed close to the bright rim. Around the central position, ammonia also shows some outflow emission although weaker than water in the red wing. Predictions of the H2O(110-101) and (111-000) antenna temperatures across the PDR are estimated with our PDR model for the forthcoming observations with the Herschel Space Observatory.

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

Irradiated ISM: Discriminating between Cosmic Rays and X-rays

The ISM of active galaxy centers is exposed to a combination of cosmic ray, FUV and X-ray radiation. We apply PDR models to this ISM with both `normal' and highly elevated (5\times 10^{-15}s^-1) cosmic-ray rates and compare the results to those obtained for XDRs. Our existing PDR-XDR code is used to construct models over a 10^3-10^5 cm^-3 density range and for 0.16-160 erg s^-1 cm^-2 impingent fluxes. We obtain larger high J (J>10) CO ratios in PDRs when we use the highly elevated cosmic ray rate, but these are always exceeded by the corresponding XDR ratios. The [CI] 609 mum/13CO(2-1) line ratio is boosted by a factor of a few in PDRs with n~10^3 cm^-3 exposed to a high cosmic ray rate. At higher densities ratios become identical irrespective of cosmic ray flux, while XDRs always show elevated [CI] emission per CO column. The HCN/CO and HCN/HCO+ line ratios, combined with high J CO emission lines, are good diagnostics to distinguish between PDRs under either low or high cosmic ray irradiation conditions, and XDRs. Hence, the HIFI instrument on Herschel, which can detect these CO lines, will be crucial in the study of active galaxies.

preprint2004arXiv

Molecular hydrogen formation on dust grains in the high redshift universe

We study the formation of molecular hydrogen on dust grain surfaces and apply our results to the high redshift universe. We find that a range of physical parameters, in particular dust temperature and gas temperature, but not so much dust surface composition, influence the formation rate of H$_2$. The H$_2$ formation rate is found to be suppressed above gas kinetic temperatures of a few hundred K and for dust temperatures above 200-300 K and below 10 K. We highlight the differences between our treatment of the H$_2$ formation process and other descriptions in the literature. We also study the relative importance of H$_2$ formation on dust grains with respect to molecular hydrogen formation in the gas phase, through the H$^-$ route. The ratio of formation rates of these two routes depends to a large part on the dust abundance, on the electron abundance, and also on the relative strength of the FUV (extra-)galactic radiation field. We find that for a cosmological evolution of the star formation rate and dust density consistent with the Madau plot, a positive feedback effect on the abundance of H$_2$ due to the presence of dust grains can occur at redshifts $z\ge 3$. This effect occurs for a dust-to-gas mass ratio as small as 10$^{-3}$ of the galactic value.

M. Spaans

What is connected

Connect this record

See the researcher in context

Building this map preview

77 published item(s)

ALMA observations of Lyman-alpha Blob 1: halo sub-structure illuminated from within

Constraining cloud parameters using high density gas tracers in galaxies

HCO+ Detection of Dust-Depleted Gas in the Inner Hole of the LkCa 15 Pre-Transitional Disk

Modelling Mechanical Heating in Star-Forming Galaxies: CO and 13CO Line Ratios as Sensitive Probes

Submillimeter H2O and H2O+ emission in lensed ultra- and hyper-luminous infrared galaxies at z ~ 2-4

Surface chemistry in photodissociation regions

The average submillimetre properties of Lyman-alpha Blobs at z=3

The mid-IR water and silicate relation in protoplanetary disks

The SCUBA-2 Cosmology Legacy Survey: 850um maps, catalogues and number counts

The SCUBA-2 Cosmology Legacy Survey: The EGS deep field I - Deep number counts and the redshift distribution of the recovered Cosmic Infrared Background at 450 and 850 um

Consistent dust and gas models for protoplanetary disks. I. Disk shape, dust settling, opacities, and PAHs

Detection of a large fraction of atomic gas not associated with star-forming material in M17 SW

Disentangling the excitation conditions of the dense gas in M17 SW

Exploring the molecular chemistry and excitation in obscured luminous infrared galaxies: An ALMA mm-wave spectral scan of NGC 4418

How chemistry influences cloud structure, star formation, and the IMF

Probing highly obscured, self-absorbed galaxy nuclei with vibrationally excited HCN

The Herschel Comprehensive (U)LIRG Emission Survey (HerCULES): CO Ladders, fine structure lines, and neutral gas cooling

The JCMT Plane Survey: early results from the l = 30 degree field

[OI] disk emission in the Taurus star forming region

A submillimeter galaxy illuminating its circumgalactic medium: Ly-alpha scattering in a cold, clumpy outflow

A UV flux constraint on the formation of direct collapse black holes

Diagnostics of the molecular component of PDRs with mechanical heating. II: line intensities and ratios

Effects of turbulence and rotation on protostar formation as a precursor to seed black holes

Future mmVLBI Research with ALMA: A European vision

High-lying OH absorption, [C II] deficits, and extreme $L_{\mathrm{FIR}}/M_{\mathrm{H2}}$ ratios in galaxies

How realistic UV spectra and X-rays suppress the abundance of direct collapse black holes

Molecular line emission in NGC1068 imaged with ALMA: II. The chemistry of the dense molecular gas

Molecular line emission in NGC1068 imaged with ALMA. I An AGN-driven outflow in the dense molecular gas

Probing the gaseous disk of T Tau N with CN 5-4 lines

Star formation relations and CO SLEDs across the J-ladder and redshift

The formation of massive primordial stars in the presence of moderate UV backgrounds

The nature of the ISM in galaxies during the star-formation activity peak of the Universe

The SCUBA-2 Cosmology Legacy Survey: ALMA resolves the rest-frame far-infrared emission of sub-millimeter galaxies

On Particle Mass Changes and GR: Space-time Topology Causes LHC Leakage

The impact of freeze-out on collapsing molecular clouds

The influence of magnetic fields, turbulence, and UV radiation on the formation of supermassive black holes

The Response of Metal Rich Gas to X-Ray Irradiation from a Massive Black Hole at High Redshift: Proof of Concept

The SCUBA-2 Cosmology Legacy Survey: demographics of the 450 μm-population

Uncertainties in water chemistry in disks: An application to TW Hya

Water vapor in the protoplanetary disk of DG Tau

Diagnostics of the Molecular Component of PDRs with Mechanical Heating

Evidence for CO shock excitation in NGC 6240 from Herschel SPIRE spectroscopy

Extended warm and dense gas towards W49A: starburst conditions in our Galaxy?

FUV and X-ray irradiated protoplanetary disks: a grid of models I. The disk structure

Hydrogenation of PAH cations: a first step towards H2 formation

Testing the Proposed Connection between Dark Energy and Black Holes

The impact of magnetic fields on the IMF in star-forming clouds near a supermassive black hole

The implications of dust for high-redshift protogalaxies and the formation of binary disks

The SCUBA-2 Cosmology Legacy Survey: blank-field number counts of 450um-selected galaxies and their contribution to the cosmic infrared background

Diffuse interstellar bands in Upper Scorpius: Probing variations in the DIB spectrum due to changing environmental conditions

Enhanced H2O formation through dust grain chemistry in X-ray exposed environments

Gas modelling in the disc of HD 163296

Interstellar ices as witnesses of star formation: selective deuteration of water and organic molecules unveiled

Lyman alpha emission from the first galaxies: Implications of UV backgrounds and the formation of molecules

Lyman alpha emission from the first galaxies: Signatures of accretion and infall in the presence of line trapping

Star formation near an obscured AGN: Variations in the initial mass function

The deeply obscured AGN of NGC4945 I. Spitzer-IRS maps of [Ne V], [Ne II], H2 0-0 S(1), S(2), and other tracers

The JCMT Spectral Legacy Survey: physical structure of the molecular envelope of the high-mass protostar AFGL2591

When sticking influences H2 formation

A new condition for the transition from runaway to oligarchic growth

Accretion among preplanetary bodies: the many faces of runaway growth

Black hole accretion and star formation as drivers of gas excitation and chemistry in Mrk231

Evolution of the ISM in Luminous IR Galaxies

Gravity and Information: Putting a Bit of Quantum into GR

Herschel observations of water vapour in Markarian 231

Star Formation in Extreme Environments: The Effects of Cosmic Rays and Mechanical Heating

The Impact of Lyman alpha Trapping on the Formation of Primordial Objects

The impact of X-rays on molecular cloud fragmentation and the IMF

X-ray impact on the protoplanetary disks around T Tauri stars

CHAMP+ observations of warm gas in M17 SW

Cooling lines as probes of the formation and buildup of galaxies and black holes

Probing X-ray irradiation in the nucleus of NGC 1068 with observations of high-J lines of dense gas tracers

The Thermodynamics of Molecular Cloud Fragmentation: Star Formation Under Non-Milky Way Conditions

Water and ammonia abundances in S140 with the Odin satellite