Source author record

Andreas Efstathiou

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

19 published item(s)

preprint2022arXiv

A hyperluminous obscured quasar at a redshift of z ~ 4.3

In this work we report the discovery of the hyperluminous galaxy HELP_J100156.75+022344.7 at the photometric redshift of z ~ 4.3. The galaxy was discovered in the Cosmological Evolution Survey (COSMOS) field, one of the fields studied by the Herschel Extragalactic Legacy Project (HELP). We present the spectral energy distribution (SED) of the galaxy and fit it with the CYprus models for Galaxies and their NUclear Spectra (CYGNUS) multi-component radiative transfer models. We find that its emission is dominated by an obscured quasar with a predicted total 1-1000um luminosity of $3.91^{+1.69}_{-0.55} \times 10^{13} L_\odot$ and an active galactic nucleus (AGN) fraction of ~89%. We also fit HELP_J100156.75+022344.7 with the Code Investigating GALaxy Emission (CIGALE) code and find a similar result. This is only the second z > 4 hyperluminous obscured quasar discovered to date. The discovery of HELP_J100156.75+022344.7 in the ~ 2deg^2 COSMOS field implies that a large number of obscured hyperluminous quasars may lie in the HELP fields which cover ~ 1300deg^2. If this is confirmed, tension between supermassive black hole evolution models and observations will be alleviated. We estimate the space density of objects like HELP_J100156.75+022344.7 at z ~ 4.5 to be $\sim 1.8 \times 10^{-8}$Mpc$^{-3}$. This is slightly higher than the space density of coeval hyperluminous optically selected quasars suggesting that the obscuring torus in z > 4 quasars may have a covering factor $\gtrsim 50\%$.

preprint2022arXiv

Detecting and analysing the topology of the cosmic web with spatial clustering algorithms I: Methods

In this paper we explore the use of spatial clustering algorithms as a new computational approach for modeling the cosmic web. We demonstrate that such algorithms are efficient in terms of computing time needed. We explore three distinct spatial methods which we suitably adjust for (i) detecting the topology of the cosmic web and (ii) categorizing various cosmic structures as voids, walls, clusters and superclusters based on a variety of topological and physical criteria such as the physical distance between objects, their masses and local densities. The methods explored are (1) a new spatial method called Gravity Lattice ; (2) a modified version of another spatial clustering algorithm, the ABACUS; and (3) the well known spatial clustering algorithm HDBSCAN. We utilize HDBSCAN in order to detect cosmic structures and categorize them using their overdensity. We demonstrate that the ABACUS method can be combined with the classic DTFE method to obtain similar results in terms of the achieved accuracy with about an order of magnitude less computation time. To further solidify our claims, we draw insights from the computer science domain and compare the quality of the results with and without the application of our method. Finally, we further extend our experiments and verify their effectiveness by showing their ability to scale well with different cosmic web structures that formed at different redshifts.

preprint2022arXiv

Radio Spectra of Luminous, Heavily Obscured WISE-NVSS Selected Quasars

We present radio spectra spanning $0.1 - 10$ GHz for the sample of heavily obscured luminous quasars with extremely red mid-infrared-optical colors and compact radio emission. The spectra are constructed from targeted 10 GHz observations and archival radio survey data, which together yield $6-11$ flux density measurements for each object. Our suite of Python tools for modeling the radio spectra is publicly available on Github. Our primary result is that most (61%) of the sample have peaked or curved radio spectra and many (36%) could be classified as Gigahertz Peaked Spectrum (GPS) sources. This indicates compact emission regions likely arising from recently triggered radio jets. Assuming synchrotron self-absorption (SSA) generates the peaks, we infer compact source sizes ($3 - 100$ pc) with strong magnetic fields ($6 - 100$ mG) and young ages ($30 - 10^4$ years). Conversely, free-free absorption (FFA) could also create peaks due to the high column densities associated with the deeply embedded nature of the sample. However, we find no correlations between the existence or frequency of the peaks and any parameters of the MIR emission. The high-frequency spectral indices are steep ($α\approx -1$) and correlate, weakly, with the ratio of MIR photon energy density to magnetic energy density, suggesting that the spectral steepening could arise from inverse Compton scattering off the intense MIR photon field. This study provides a foundation for combining multi-frequency and mixed-resolution radio survey data for understanding the impact of young radio jets on the ISM and star formation rates of their host galaxies.

preprint2022arXiv

Stellar and black hole assembly in z<0.3 infrared-luminous mergers: intermittent starbursts vs. super-Eddington accretion

We study stellar and black hole mass assembly in a sample of 42 infrared-luminous galaxy mergers at z<0.3 by combining results from radiative transfer modelling with archival measures of molecular gas and black hole mass. The ratios of stellar mass, molecular gas mass, and black hole mass to each other are consistent with those of massive gas-rich galaxies at z<0.3. The advanced mergers may show increased black hole mass to stellar mass ratios, consistent with the transition from AGN to ellipticals and implying substantial black hole mass growth over the course of the merger. Star formation rates are enhanced relative to the local main sequence, by factors of ~100 in the starburst and ~1.8 in the host. The starburst star formation rates appear distinct to star formation in the main sequence at all redshifts up to at least z~5. Starbursts may prefer late-stage mergers, but are observed at any merger stage. We do not find evidence that the starbursts in these low-redshift systems substantially increase the total stellar mass, with a soft upper limit on the stellar mass increase from starburst activity of about a factor of two. In contrast, 12 objects show evidence for super-Eddington accretion, associated with late-stage mergers, suggesting that many AGN in infrared-luminous mergers go through a super-Eddington phase. The super-Eddington phase may increase black hole mass by up to an order of magnitude at an accretion efficiency of 42+/-33% over a period of 44+/-22Myr. Our results imply that super-Eddington accretion is an important black hole growth channel in infrared-luminous galaxies at all redshifts.

preprint2021arXiv

The viewing angle in AGN SED models, a data-driven analysis

The validity of the unified active galactic nuclei (AGN) model has been challenged in the last decade, especially when different types of AGNs are considered to only differ in the viewing angle to the torus. We aim to assess the importance of the viewing angle in classifying different types of Seyfert galaxies in spectral energy distribution (SED) modelling. We retrieve photometric data from publicly available astronomical databases: CDS and NED, to model SEDs with X-CIGALE in a sample of 13 173 Seyfert galaxies located at redshift range from $z=0$ to $z=3.5$, with a median redshift of $z\approx0.2$. We assess whether the estimated viewing angle from the SED models reflects different Seyfert classifications. Two AGN models with either a smooth or clumpy torus structure are adopted in this paper. We find that the viewing angle in Type-1 AGNs is better constrained than in Type-2 AGNs. Limiting the viewing angles representing these two types of AGNs do not affect the physical parameter estimates such as star-formation rate (SFR) or AGN fractional contribution ($f_{\rm{AGN}}$). In addition, the viewing angle is not the most discriminating physical parameter to differentiate Seyfert types. We suggest that the observed and intrinsic AGN disc luminosity can: i) be used in $z<0.5$ studies to distinguish between Type-1 and Type-2 AGNs, and ii) explain the probable evolutionary path between these AGN types. Finally, we propose the use of X-CIGALE for AGN galaxy classification tasks. All data from the 13 173 SED fits are available at https://doi.org/10.5281/zenodo.5221764

preprint2020arXiv

High-resolution VLA Imaging of Obscured Quasars: Young Radio Jets Caught in a Dense ISM

We present new sub-arcsecond-resolution Karl G. Jansky Very Large Array (VLA) imaging at 10 GHz of 155 ultra-luminous ($L_{\rm bol}\sim10^{11.7-14.2} L_\odot$) and heavily obscured quasars with redshifts $z \sim0.4-3$. The sample was selected to have extremely red mid-infrared (MIR)-optical color ratios based on data from Wide-Field Infrared Survey Explorer (WISE) along with a detection of bright, unresolved radio emission from the NRAO VLA Sky Survey (NVSS) or Faint Images of the Radio Sky at Twenty-Centimeters (FIRST) Survey. Our high-resolution VLA observations have revealed that the majority of the sources in our sample (93 out of 155) are compact on angular scales $<0.2^{\prime \prime}$ ($\leq 1.7$ kpc at $z \sim2$). The radio luminosities, linear extents, and lobe pressures of our sources are similar to young radio active galactic nuclei (AGN; e.g., Gigahertz Peaked Spectrum, GPS, and Compact Steep Spectrum, CSS, sources), but their space density is considerably lower. Application of a simple adiabatic lobe expansion model suggests relatively young dynamical ages ($\sim10^{4-7}$ years), relatively high ambient ISM densities ($\sim1-10^4$ cm$^{-3}$), and modest lobe expansion speeds ($\sim30-10,000$ km s$^{-1}$). Thus, we find our sources to be consistent with a population of newly triggered, young jets caught in a unique evolutionary stage in which they still reside within the dense gas reservoirs of their hosts. Based on their radio luminosity function and dynamical ages, we estimate only $\sim20\%$ of classical large scale FRI/II radio galaxies could have evolved directly from these objects. We speculate that the WISE-NVSS sources might first become GPS or CSS sources, of which some might ultimately evolve into larger radio galaxies.

preprint2016arXiv

HERUS: A CO Atlas from SPIRE Spectroscopy of local ULIRGs

We present the Herschel SPIRE Fourier Transform Spectroscopy (FTS) atlas for a complete flux limited sample of local Ultra-Luminous Infra-Red Galaxies as part of the HERschel ULIRG Survey (HERUS). The data reduction is described in detail and was optimized for faint FTS sources with particular care being taken with the subtraction of the background which dominates the continuum shape of the spectra. Special treatment in the data reduction has been given to any observation suffering from artefacts in the data caused by anomalous instrumental effects to improve the final spectra. Complete spectra are shown covering $200 - 671μ$m with photometry in the SPIRE bands at 250$μ$m, 350$μ$m and 500$μ$m. The spectra include near complete CO ladders for over half of our sample, as well as fine structure lines from [CI] 370 $μ$m, [CI] 609 $μ$m, and [NII] 205 $μ$m. We also detect H$_{2}$O lines in several objects. We construct CO Spectral Line Energy Distributions (SLEDs) for the sample, and compare their slopes with the far-infrared colours and luminosities. We show that the CO SLEDs of ULIRGs can be broadly grouped into three classes based on their excitation. We find that the mid-J (5$<$J$<$8) lines are better correlated with the total far-infrared luminosity, suggesting that the warm gas component is closely linked to recent star-formation. The higher J transitions do not linearly correlate with the far-infrared luminosity, consistent with them originating in hotter, denser gas unconnected to the current star-formation. {\bf We conclude that in most cases more than one temperature components are required to model the CO SLEDs.}

preprint2016arXiv

Star Formation in Luminous Quasars at 2<z<3

We investigate the relation between star formation rates ($\dot{M}_{s}$) and AGN properties in optically selected type 1 quasars at $2<z<3$ using data from Herschel and the SDSS. We find that $\dot{\rm{M}}_s$ remains approximately constant with redshift, at $300\pm100~\rm{M}_{\odot}$yr$^{-1}$. Conversely, $\dot{\rm{M}}_s$ increases with AGN luminosity, up to a maximum of $\sim600~\rm{M}_{\odot}$yr$^{-1}$, and with CIV FWHM. In context with previous results, this is consistent with a relation between $\dot{\rm{M}}_s$ and black hole accretion rate ($\dot{\rm{M}}_{bh}$) existing in only parts of the $z-\dot{\rm{M}}_{s}-\dot{\rm{M}}_{bh}$ plane, dependent on the free gas fraction, the trigger for activity, and the processes that may quench star formation. The relations between $\dot{\rm{M}}_s$ and both AGN luminosity and CIV FWHM are consistent with star formation rates in quasars scaling with black hole mass, though we cannot rule out a separate relation with black hole accretion rate. Star formation rates are observed to decline with increasing CIV equivalent width. This decline can be partially explained via the Baldwin effect, but may have an additional contribution from one or more of three factors; $M_i$ is not a linear tracer of L$_{2500}$, the Baldwin effect changes form at high AGN luminosities, and high CIV EW values signpost a change in the relation between $\dot{\rm{M}}_s$ and $\dot{\rm{M}}_{bh}$. Finally, there is no strong relation between $\dot{\rm{M}}_s$ and Eddington ratio, or the asymmetry of the CIV line. The former suggests that star formation rates do not scale with how efficiently the black hole is accreting, while the latter is consistent with CIV asymmetries arising from orientation effects.

preprint2015arXiv

Radio Jet Feedback and Star Formation in Heavily Obscured Quasars at Redshifts ~0.3-3, I: ALMA Observations

We present ALMA 870 micron (345 GHz) data for 49 high redshift (0.47<z<2.85), luminous (11.7 < log L(bol) (Lsun) < 14.2) radio-powerful AGN, obtained to constrain cool dust emission from starbursts concurrent with highly obscured radiative-mode black hole (BH) accretion in massive galaxies which possess a small radio jet. The sample was selected from WISE with extremely steep (red) mid-infrared (MIR) colors and with compact radio emission from NVSS/FIRST. Twenty-six sources are detected at 870 microns, and we find that the sample has large mid- to far-infrared luminosity ratios consistent with a dominant and highly obscured quasar. The rest-frame 3 GHz radio powers are 24.7 < log P3.0 GHz (W/Hz) < 27.3, and all sources are radio-intermediate or radio-loud. BH mass estimates are 7.7 < log M(BH) (Msun) < 10.2. The rest frame 1-5 um SEDs are very similar to the "Hot DOGs" (Hot Dust Obscured Galaxies), and steeper (redder) than almost any other known extragalactic sources. ISM masses estimated for the ALMA detected sources are 9.9 < log M(ISM) (Msun) < 11.75 assuming a dust temperature of 30K. The cool dust emission is consistent with star formation rates (SFRs) reaching several thousand Msun/yr, depending on the assumed dust temperature, however we cannot rule out the alternative that the AGN powers all the emission in some cases. Our best constrained source has radiative transfer solutions with ~ equal contributions from an obscured AGN and a young (10-15 Myr) compact starburst.

preprint2015arXiv

Self-consistent 2-phase AGN torus models: SED library for observers

We assume that dust near active galactic nuclei (AGN) is distributed in a torus-like geometry, which may be described by a clumpy medium or a homogeneous disk or as a combination of the two (i.e. a 2-phase medium). The dust particles considered are fluffy and have higher submillimeter emissivities than grains in the diffuse ISM. The dust-photon interaction is treated in a fully self-consistent three dimensional radiative transfer code. We provide an AGN library of spectral energy distributions (SEDs). Its purpose is to quickly obtain estimates of the basic parameters of the AGN, such as the intrinsic luminosity of the central source, the viewing angle, the inner radius, the volume filling factor and optical depth of the clouds, and the optical depth of the disk midplane, and to predict the flux at yet unobserved wavelengths. The procedure is simple and consists of finding an element in the library that matches the observations. We discuss the general properties of the models and in particular the 10mic. silicate band. The AGN library accounts well for the observed scatter of the feature strengths and wavelengths of the peak emission. AGN extinction curves are discussed and we find that there is no direct one-to-one link between the observed extinction and the wavelength dependence of the dust cross sections. We show that objects of the library cover the observed range of mid IR colors of known AGN. The validity of the approach is demonstrated by matching the SEDs of a number of representative objects: Four Seyferts and two quasars for which we present new Herschel photometry, two radio galaxies, and one hyperluminous infrared galaxy. Strikingly, for the five luminous objects we find pure AGN models fit the SED without a need to postulate starburst activity.

preprint2014arXiv

A multi-wavelength view of the central kiloparsec region in the Luminous Infrared Galaxy NGC1614

The Luminous Infrared Galaxy NGC1614 hosts a prominent circumnuclear ring of star formation. However, the nature of the dominant emitting mechanism in its central ~100 pc is still under debate. We present sub-arcsecond angular resolution radio, mid-infrared, Pa-alpha, optical, and X-ray observations of NGC1614, aimed at studying in detail both the circumnuclear ring and the nuclear region. The 8.4 GHz continuum emission traced by the Very Large Array (VLA) and the Gemini/T-ReCS 8.7 micron emission, as well as the Pa-alpha line emission, show remarkable morphological similarities within the star-forming ring, suggesting that the underlying emission mechanisms are tightly related. We used an HST/NICMOS Pa-alpha map of similar resolution to our radio maps to disentangle the thermal free-free and non-thermal synchrotron radio emission, from which we obtained the intrinsic synchrotron power-law for each individual region within the central kpc of NGC1614. The radio ring surrounds a relatively faint, steep-spectrum source at the very center of the galaxy, suggesting that the central source is not powered by an AGN, but rather by a compact (r < 90 pc) starburst. Chandra X-ray data also show that the central kpc region is dominated by starburst activity, without requiring the existence of an AGN. We also used publicly available infrared data to model-fit the spectral energy distribution of both the starburst ring and a putative AGN in NGC1614. In summary, we conclude that there is no need to invoke an AGN to explain the observed bolometric properties of the galaxy.

preprint2013arXiv

Far-infrared Fine-Structure Line Diagnostics of Ultraluminous Infrared Galaxies

We present Herschel observations of six fine-structure lines in 25 Ultraluminous Infrared Galaxies at z<0.27. The lines, [O III]52, [N III]57, [O I]63, [N II]122, [O I]145, and [C II]158, are mostly single gaussians with widths <600 km s-1 and luminosities of 10^7 - 10^9 Solar. There are deficits in the [O I]63/L_IR, [N II]/L_IR, [O I]145/L_IR, and [C II]/L_IR ratios compared to lower luminosity systems. The majority of the line deficits are consistent with dustier H II regions, but part of the [C II] deficit may arise from an additional mechanism, plausibly charged dust grains. This is consistent with some of the [C II] originating from PDRs or the ISM. We derive relations between far-IR line luminosities and both IR luminosity and star formation rate. We find that [N II] and both [O I] lines are good tracers of IR luminosity and star formation rate. In contrast, [C II] is a poor tracer of IR luminosity and star formation rate, and does not improve as a tracer of either quantity if the [C II] deficit is accounted for. The continuum luminosity densities also correlate with IR luminosity and star formation rate. We derive ranges for the gas density and ultraviolet radiation intensity of 10^1 < n < 10^2.5 and 10^2.2 < G_0 < 10^3.6, respectively. These ranges depend on optical type, the importance of star formation, and merger stage. We do not find relationships between far-IR line properties and several other parameters; AGN activity, merger stage, mid-IR excitation, and SMBH mass. We conclude that these far-IR lines arise from gas heated by starlight, and that they are not strongly influenced by AGN activity.

preprint2012arXiv

A new model for the infrared emission of IRAS F10214+4724

We present a new model for the infrared emission of the high redshift hyperluminous infrared galaxy IRAS F10214+4724 which takes into account recent photometric data from Spitzer and Herschel that sample the peak of its spectral energy distribution. We first demonstrate that the combination of the AGN tapered disc and starburst models of Efstathiou and coworkers, while able to give an excellent fit to the average spectrum of type 2 AGN measured by Spitzer, fails to match the spectral energy distribution of IRAS F10214+4724. This is mainly due to the fact that the nuSnu distribution of the galaxy falls very steeply with increasing frequency (a characteristic of heavy absorption by dust) but shows a silicate feature in emission. We propose a model that assumes two components of emission: clouds that are associated with the narrow-line region and a highly obscured starburst. The emission from the clouds must suffer significantly stronger gravitational lensing compared to the emission from the torus to explain the observed spectral energy distribution.

preprint2012arXiv

Principal Component Analysis and Radiative Transfer modelling of Spitzer IRS Spectra of Ultra Luminous Infrared Galaxies

The mid-infrared spectra of ultraluminous infrared galaxies (ULIRGs) contain a variety of spectral features that can be used as diagnostics to characterise the spectra. However, such diagnostics are biased by our prior prejudices on the origin of the features. Moreover, by using only part of the spectrum they do not utilise the full information content of the spectra. Blind statistical techniques such as principal component analysis (PCA) consider the whole spectrum, find correlated features and separate them out into distinct components. We further investigate the principal components (PCs) of ULIRGs derived in Wang et al.(2011). We quantitatively show that five PCs is optimal for describing the IRS spectra. These five components (PC1-PC5) and the mean spectrum provide a template basis set that reproduces spectra of all z<0.35 ULIRGs within the noise. For comparison, the spectra are also modelled with a combination of radiative transfer models of both starbursts and the dusty torus surrounding active galactic nuclei. The five PCs typically provide better fits than the models. We argue that the radiative transfer models require a colder dust component and have difficulty in modelling strong PAH features. Aided by the models we also interpret the physical processes that the principal components represent. The third principal component is shown to indicate the nature of the dominant power source, while PC1 is related to the inclination of the AGN torus. Finally, we use the 5 PCs to define a new classification scheme using 5D Gaussian mixtures modelling and trained on widely used optical classifications. The five PCs, average spectra for the four classifications and the code to classify objects are made available at: http://www.phys.susx.ac.uk/~pdh21/PCA/

preprint2011arXiv

Direct Evidence for Termination of Obscured Star Formation by Radiatively Driven Outflows in FeLoBAL QSOs

We use SDSS spectra and optical to far-infrared photometry for a sample of 31 FeLoBAL QSOs to study the relationship between the AGN-driven outflows, and obscured star formation in their host galaxies. We find that FeLoBAL QSOs invariably have IR luminosities exceeding 10^12 Solar luminosities. The AGN supplies 75% of the total IR emission, on average, but with a range from 20% to 100%. We find a clear anticorrelation between the strength of the AGN-driven outflows and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN act to curtail obscured star formation in the host galaxies of reddened QSOs to less than ~25% of the total IR luminosity. This is the most direct evidence yet obtained for `quasar mode' AGN feedback.

preprint2011arXiv

Direct Evidence for Termination of Obscured Star Formation by Radiatively Driven Outflows in Reddened QSOs

We present optical to far-infrared photometry of 31 reddened QSOs that show evidence for radiatively driven outflows originating from AGN in their rest-frame UV spectra. We use these data to study the relationships between the AGN-driven outflows, and the AGN and starburst infrared luminosities. We find that FeLoBAL QSOs are invariably IR-luminous, with IR luminosities exceeding 10^{12} Solar luminosities in all cases. The AGN supplies 76% of the total IR emission, on average, but with a range from 20% to 100%. We find no evidence that the absolute luminosity of obscured star formation is affected by the AGN-driven outflows. Conversely, we find an anticorrelation between the strength of AGN-driven outflows, as measured from the range of outflow velocities over which absorption exceeds a minimal threshold, and the contribution from star formation to the total IR luminosity, with a much higher chance of seeing a starburst contribution in excess of 25% in systems with weak outflows than in systems with strong outflows. Moreover, we find no convincing evidence that this effect is driven by the IR luminosity of the AGN. We conclude that radiatively driven outflows from AGN can have a dramatic, negative impact on luminous star formation in their host galaxies. We find that such outflows act to curtail star formation such that star formation contributes less than ~25% of the total IR luminosity. We also propose that the degree to which termination of star formation takes place is not deducible from the IR luminosity of the AGN.

preprint2009arXiv

A far-infrared survey at the North Galactic Pole I: Nearby star-forming galaxies and the effect of confused sources on source counts

We present follow-up observations of the far-infrared (FIR) sources at 90, 150 and 180 mu detected as part of the ISOPHOT EBL project, which has recently measured the absolute surface brightness of the cosmic infrared background radiation (CIRB) for the first time independently from COBE data. We have observed the fields at the North Galactic Pole region in the optical and near-IR, and complement these data with SDSS photometry, and spectroscopy where available, and present identifications of the 25 FIR sources which reach down to ~150 mJy in all three ISOPHOT bands. Identifications are done by means of full spectral energy density fitting to all sources in the FIR error circle areas. Approximately 80 per cent are identified as star-forming or star-bursting galaxies at z<0.3. We also find that more than half of the counterparts have disturbed morphologies, with signs of past or present interactions. However, only 20 per cent of all the sources are uniquely matched with a single galaxy -- 40 per cent are blends of two or more of these nearby star-forming galaxies, while another 20 per cent are likely blends of nearby and fainter galaxies. The final 20 per cent are likely to be more luminous IR galaxies at higher redshifts. The blended sources have an effect on the FIR source counts. In particular, taking into account realistic confusion or blending of sources, the differential FIR counts move down by a factor of ~1.5 and steepen in the 100 to 400 mJy range.

preprint2009arXiv

Modelling the Spoon IRS diagnostic diagram

We explore whether our models for starbursts, quiescent star-forming galaxies and for AGN dust tori are able to model the full range of IRS spectra measured with Spitzer. The diagnostic plot of 9.7 mu silicate optical depth versus 6.2 mu PAH equivalent width, introduced by Spoon and coworkers in 2007, gives a good indication of the age and optical depth of a starburst, and of the contribution of an AGN dust torus. However there is aliasing between age and optical depth at later times in the evolution of a starburst, and between age and the presence of an AGN dust torus. Modeling the full IRS spectra and using broad-band 25-850 mu fluxes can help to resolve these aliases. The observed spectral energy distributions require starbursts of a range of ages with initial dust optical depth ranging from 50-200, optically thin dust emission ('cirrus') illuminated by a range of surface brightnesses of the interstellar radiation field, and AGN dust tori with a range of viewing angles.

preprint2009arXiv

Starburst and cirrus models for submillimeter galaxies

We present radiative transfer models for submillimeter galaxies with spectroscopic redshifts and mid-infrared spectroscopy from Spitzer/IRS and analyze available Spitzer/MIPS 24, 70 and 160mu data. We use two types of starburst models, a cirrus model and a model for the emission of an AGN torus in order to investigate the nature of these objects. We find that for three of the objects (25%) cirrus emission alone can account for the mid-infrared spectrum and the MIPS and submillimeter data. For the remaining objects we find that we need a combination of starburst and cirrus in order to fit simultaneously the multi--wavelength data. We find that the typical submillimeter galaxy has comparable luminosity in the starburst (median L=10^12.5 Lo) and cirrus (median L=10^12.4 Lo) components. This could arise if the galaxies have been forming stars continuously for the last 250Myr with the star formation occurring in the last 5Myr being shrouded by high-optical-depth molecular cloud dust, whereas the rest of the starlight is attenuated by diffuse dust or cirrus with an Av of about 1mag.