Spectral Energy Distribution Analysis Research Articles

Spectrophotometric analysis of gamma-ray burst afterglow extinction curves with X-Shooter

In this work we use gamma-ray burst (GRB) afterglow spectra observed with the VLT/X-shooter spectrograph to measure rest-frame extinction in GRB lines-of-sight by modeling the broadband near-infrared (NIR) to X-ray afterglow spectral energy distributions (SEDs). Our sample consists of nine Swift GRBs, eight of them belonging to the long-duration and one to the short-duration class. Dust is modeled using the average extinction curves of the Milky Way and the two Magellanic Clouds. We derive the rest-frame extinction of the entire sample, which fall in the range $0 \lesssim {\it A}_{\rm V} \lesssim 1.2$. Moreover, the SMC extinction curve is the preferred extinction curve template for the majority of our sample, a result which is in agreement with those commonly observed in GRB lines-of-sights. In one analysed case (GRB 120119A), the common extinction curve templates fail to reproduce the observed extinction. To illustrate the advantage of using the high-quality X-shooter afterglow SEDs over the photometric SEDs, we repeat the modeling using the broadband SEDs with the NIR-to-UV photometric measurements instead of the spectra. The main result is that the spectroscopic data, thanks to a combination of excellent resolution and coverage of the blue part of the SED, are more successful in constraining the extinction curves and therefore the dust properties in GRB hosts with respect to photometric measurements. In all cases but one the extinction curve of one template is preferred over the others. We show that the modeled values of the extinction and the spectral slope, obtained through spectroscopic and photometric SED analysis, can differ significantly for individual events. Finally we stress that, regardless of the resolution of the optical-to-NIR data, the SED modeling gives reliable results only when the fit is performed on a SED covering a broader spectral region.

A CONNECTION BETWEEN OBSCURATION AND STAR FORMATION IN LUMINOUS QUASARS

We present a measurement of the star formation properties of a uniform sample of mid-IR selected, unobscured and obscured quasars (QSO1s and QSO2s) in the Bo\"otes survey region. We use an spectral energy distribution (SED) analysis for photometric data spanning optical to far-IR wavelengths to decompose AGN and host galaxy components. We find that when compared to a matched sample of QSO1s, the QSO2s have higher far-IR detection fractions, far-IR fluxes and infrared star formation luminosities ($L_{\rm IR}^{\rm SF}$) by a factor of $\sim2$. Correspondingly, we show that the AGN obscured fraction rises from 0.3 to 0.7 between $4-40\times10^{11}L_\odot$. We also find evidence associating the absorption in the X-ray emission with the presence of far-IR emitting dust. Overall, these results are consistent with galaxy evolution models in which quasar obscurations can be associated with a dust-enshrouded starburst galaxies.

Star formation inz> 1 3CR host galaxies as seen byHerschel

We present Herschel (PACS and SPIRE) far-infrared (FIR) photometry of a complete sample of z>1 3CR sources, from the Herschel GT project The Herschel Legacy of distant radio-loud AGN (PI: Barthel). Combining these with existing Spitzer photometric data, we perform an infrared (IR) spectral energy distribution (SED) analysis of these landmark objects in extragalactic research to study the star formation in the hosts of some of the brightest active galactic nuclei (AGN) known at any epoch. Accounting for the contribution from an AGN-powered warm dust component to the IR SED, about 40% of our objects undergo episodes of prodigious, ULIRG-strength star formation, with rates of hundreds of solar masses per year, coeval with the growth of the central supermassive black hole. Median SEDs imply that the quasar and radio galaxy hosts have similar FIR properties, in agreement with the orientation-based unification for radio-loud AGN. The star-forming properties of the AGN hosts are similar to those of the general population of equally massive non-AGN galaxies at comparable redshifts, thus there is no strong evidence of universal quenching of star formation (negative feedback) within this sample. Massive galaxies at high redshift may be forming stars prodigiously, regardless of whether their supermassive black holes are accreting or not.

Study of morphology and stellar content of the Galactic H ii region IRAS 16148−5011

An investigation of the IRAS 16148−5011 region – a cluster at a distance of 3.6 kpc – is presented here, carried out using multiwavelength data in near-infrared (NIR) from the 1.4 m Infrared Survey Facility telescope, mid-infrared (MIR) from the archival Spitzer GLIMPSE (Galactic Legacy Infrared Midplane Survey Extraordinaire) survey, far-infrared (FIR) from the Herschel archive, and low-frequency radio continuum observations at 1280 and 843 MHz from the Giant Metrewave Radio Telescope and Molonglo Survey archive, respectively. A combination of NIR and MIR data is used to identify 7 Class I and 133 Class II sources in the region. Spectral energy distribution (SED) analysis of selected sources reveals a 9.6 M⊙ high-mass source embedded in nebulosity. However, Lyman continuum luminosity calculation using radio emission – which shows a compact H ii region – indicates the spectral type of the ionizing source to be earlier than B0-O9.5. Free–free emission SED modelling yields the electron density as 138 cm−3, and thus the mass of the ionized hydrogen as ∼16.4 M⊙. Thermal dust emission modelling, using the FIR data from Herschel and performing modified blackbody fits, helped us construct the temperature and column density maps of the region, which show peak values of 30 K and 3.3 × 1022 cm−2, respectively. The column density maps reveal an AV > 20 mag extinction associated with the nebular emission, and weak filamentary structures connecting dense clumps. The clump associated with this IRAS object is found to have dimensions of ∼ 1.1 pc × 0.8 pc, and a mass of 1023 M⊙.

Modelling and interpretation of SEDs

Circumstellar disks are mostly detected by larger continuum fluxes in the in- frared to mm spectral regions as compared to naked stars (a flux excess). The analysis of the Spectral Energy Distribution (SED) of that flux excess was crucial for the develop- ment of the first theories about protoplanetary disks, and even nowadays, it is still one of the major tools to physically characterise the disks in terms of their mass, inner holes and gaps, vertical extension & shape, dust properties, and evolutionary state. In this chapter, we will review some of the early simple theories, show some examples, discuss the influ- ence of typical disk shape and dust size parameters in modern SED analysis, and discuss how degenerate the results can be.

Iris: An extensible application for building and analyzing spectral energy distributions

Iris is an extensible application that provides astronomers with a user-friendly interface capable of ingesting broad-band data from many different sources in order to build, explore, and model spectral energy distributions (SEDs). Iris takes advantage of the standards defined by the International Virtual Observatory Alliance, but hides the technicalities of such standards by implementing different layers of abstraction on top of them. Such intermediate layers provide hooks that users and developers can exploit in order to extend the capabilities provided by Iris. For instance, custom Python models can be combined in arbitrary ways with the Iris built-in models or with other custom functions. As such, Iris offers a platform for the development and integration of SED data, services, and applications, either from the user's system or from the web. In this paper we describe the built-in features provided by Iris for building and analyzing SEDs. We also explore in some detail the Iris framework and software development kit, showing how astronomers and software developers can plug their code into an integrated SED analysis environment.

HERSCHELPHOTOMETRY OF DISKS AROUND LOW-MASS STARS IN THE R CrA CLOUD

We report photometric results from a subset of a Herschel-PACS program to observe cool dust in disks around low-mass stars as a complement to our earlier program to measure far-infrared emission from brown dwarfs. In this latest study we observed five low-mass objects in the nearby R Corona Australis region and detected at least three at 70 mu m. Using a Monte Carlo radiative transfer code we have investigated the disk masses and geometry based on detailed spectral energy distribution (SED) modeling, and we compare these new results to those from our earlier larger sample of brown dwarfs. In particular, our SED analysis for these five objects shows again that disk geometries of brown dwarfs or low-mass stars are generally similar to their higher mass counterparts like T Tauri disks, but the range of disk mass extends to well below the value found in T Tauri stars.

The dark nature of GRB 130528A and its host galaxy

We study the dark nature of GRB 130528A through multi-wavelength observations and conclude that the main reason for the optical darkness is local extinction inside of the host galaxy. Automatic observations were performed at BOOTES-4/MET robotic telescope. We also triggered target of opportunity (ToO) observation at the OSN, IRAM PdBI and the GTC+OSIRIS. The host galaxy photometric observations in optical to near-infrared (nIR) wavelengths were achieved through large ground-based aperture telescopes, such as the 10.4m GTC, the 4.2m WHT, 6m BTA, and the 2m LT. Based on these observations, spectral energy distributions (SED) for the host galaxy and afterglow were constructed. Thanks to mm observations at PdBI, we confirm the presence of a mm source within the XRT error circle that faded over the course of our observations and identify the host galaxy. However, we do not find any credible optical source within early observations with BOOTES-4/MET and 1.5m OSN telescopes. Spectroscopic observation of this galaxy by GTC showed a single faint emission line that likely corresponds to [OII] 3727\{AA} at a redshift of 1.250+/-0.001 implying a SFR(M_sun/yr) > 6.18 M_sun/yr without correcting for dust extinction. The probable extinction was revealed through analysis of the afterglow SED, resulting in a value of AV >= ~ 0.9 at the rest frame, this is comparable to extinction levels found among other dark GRBs. The SED of the host galaxy is explained well (chi2/d.o.f.=0.564) by a luminous (MB=-21.16), low-extinction (AV =0, rest frame), and aged (2.6 Gyr) stellar population. We can explain this apparent contradiction in global and line-of-sight extinction if the GRB birth place happened to lie in a local dense environment. In light of having relatively small specific SFR (SSFR) ~ 5.3 M_sun/yr (L/L_star)-1, this also could explain the age of the old stellar population of host galaxy.

A new method for classifying galaxy SEDs from multiwavelength photometry

We present a new method to classify the broad band optical-NIR spectral energy distributions (SEDs) of galaxies using three shape parameters (super-colours) based on a Principal Component Analysis of model SEDs. As well as providing a compact representation of the wide variety of SED shapes, the method allows for easy visualisation of information loss and biases caused by the incomplete sampling of the rest-frame SED as a function of redshift. We apply the method to galaxies in the UKIDSS Ultra Deep Survey with 0.9<z<1.2, and confirm our classifications by stacking rest-frame optical spectra for a fraction of objects in each class. As well as cleanly separating a tight red-sequence from star-forming galaxies, three unusual populations are identifiable by their unique colours: very dusty star-forming galaxies with high metallicity and old mean stellar age; post-starburst galaxies which have formed greater than around 10% of their mass in a recent unsustained starburst event; and metal-poor quiescent dwarf galaxies. We find that quiescent galaxies account for 45% of galaxies with log(M*/Msol)>11, declining steadily to 13% at log(M*/Msol)=10. The properties and mass-function of the post-starburst galaxies are consistent with a scenario in which gas-rich mergers contribute to the growth of the low and intermediate mass range of the red sequence.

NuSTAROBSERVATIONS OF HEAVILY OBSCURED QUASARS ATz∼ 0.5

We present NuSTAR hard X-ray (3-79 keV) observations of three Type 2 quasars at z ~ 0.4-0.5, optically selected from the Sloan Digital Sky Survey (SDSS). Although the quasars show evidence for being heavily obscured Compton-thick systems on the basis of the 2-10 keV to [OIII] luminosity ratio and multiwavelength diagnostics, their X-ray absorbing column densities (N_H) are poorly known. In this analysis: (1) we study X-ray emission at >10 keV, where X-rays from the central black hole are relatively unabsorbed, in order to better constrain N_H; (2) we further characterize the physical properties of the sources through broad-band near-UV to mid-IR spectral energy distribution (SED) analyses. One of the quasars is detected with NuSTAR at >8 keV with a no-source probability of <0.1%, and its X-ray band ratio suggests near Compton-thick absorption with N_H \gtrsim 5 x 10^23 cm^-2. The other two quasars are undetected, and have low X-ray to mid-IR luminosity ratios in both the low energy (2-10 keV) and high energy (10-40 keV) X-ray regimes that are consistent with extreme, Compton-thick absorption (N_H \gtrsim 10^24 cm^-2). We find that for quasars at z ~ 0.5, NuSTAR provides a significant improvement compared to lower energy (<10 keV) Chandra and XMM-Newton observations alone, as higher column densities can now be directly constrained.

Studying the molecular ambient towards the young stellar object EGO G35.04−0.47

Abstract We are performing a systematic study of the interstellar medium around extended green objects (EGOs), likely massive young stellar objects driving outflows. EGO G35.04−0.47 is located towards a dark cloud at the north-western edge of an H ii region. Recently, H2 jets were discovered towards this source, mainly towards its south-west, where the H2 1–0 S(1) emission peaks. Therefore, the source was catalogued as the molecular hydrogen emission-line object MHO 2429. In order to study the molecular ambient towards this star-forming site, we observed a region around the aforementioned EGO using the Atacama Submillimeter Telescope Experiment in the 12CO J = 3–2, 13CO J = 3–2, HCO+ J = 4–3, and CS J = 7–6 lines with an angular and spectral resolution of 22′′ and 0.11 km s−1, respectively. The observations revealed a molecular clump where the EGO is embedded at vLSR ∼ 51 km s−1, in coincidence with the velocity of a Class I 95 GHz methanol maser previously detected. Analyzing the 12CO line we discovered high velocity molecular gas in the range from 34 to 47 km s−1, most likely a blueshifted outflow driven by the EGO. The alignment and shape of this molecular structure coincide with those of the south-west lobe of MHO 2429, mainly between 46 and 47 km s−1, confirming that we are mapping its CO counterpart. Performing a spectral energy distribution analysis of EGO G35.04−0.47, we found that its central object should be an intermediate-mass young stellar object accreting mass at a rate similar to those found in some massive YSOs. We suggest that this source can become a massive YSO.

VOSA: SED building and analysis of thousands of stars in the framework of Gaia

VOSA (http://svo2.cab.inta-csic.es/theory/vosa/), is a web-based tool designed to combine private photometric measurements with data available in VO services distributed worldwide to build the observational spectral energy distributions (SEDs) of hundreds of objects. VOSA also accesses various collections of models to simulate the equivalent theoretical SEDs, allows the user to decide the range of physical parameters to explore, performs the SED comparison, provides the best fitting models to the user following two different approaches (chi square and Bayesian fitting), and, for stellar sources, compares these parameters with isochrones and evolutionary tracks to estimate masses and ages. In particular, VOSA offers the advantage of deriving physical parameters using all the available photometric information instead of a restricted subset of colors. VOSA was firstly released in 2008 and its functionalities are described in Bayo et al. (2008). At the time of writing there are more than 300 active users in VOSA who have published more than 60 refereed papers. In the framework of the GENIUS (https://gaia.am.ub.es/Twiki/bin/view/GENIUS) project we are upgrading VOSA to, on one hand, provide a seamless access to Gaia data and, on the other hand, handle thousands of objects at a time.In this poster, the main functionalities to be implemented in the Gaia context will be described. The poster can be found at: http://svo.cab.inta-csic.es/files/svo//Public/SVOPapers/posters/vosa-poster3.pdf.

Satmc: Spectral energy distribution Analysis Through Markov Chains

We present the general purpose spectral energy distribution (SED) fitting tool SED Analysis Through Markov Chains (SATMC). Utilizing Monte Carlo Markov Chain (MCMC) algorithms, SATMC fits an observed SED to SED templates or models of the user's choice to infer intrinsic parameters, generate confidence levels and produce the posterior parameter distribution. Here we describe the key features of SATMC from the underlying MCMC engine to specific features for handling SED fitting. We detail several test cases of SATMC, comparing results obtained to traditional least-squares methods, which highlight its accuracy, robustness and wide range of possible applications. We also present a sample of submillimetre galaxies that have been fitted using the SED synthesis routine GRASIL as input. In general, these SMGs are shown to occupy a large volume of parameter space, particularly in regards to their star formation rates which range from ~30-3000 M_sun yr^-1 and stellar masses which range from ~10^10-10^12 M_sun. Taking advantage of the Bayesian formalism inherent to SATMC, we also show how the fitting results may change under different parametrizations (i.e., different initial mass functions) and through additional or improved photometry, the latter being crucial to the study of high-redshift galaxies.

Dusty Universe viewed by AKARI far infrared detector

We present the results of the analysis of multiwavelength Spectral Energy Distributions (SEDs) of far-infrared galaxies detected in the AKARI Deep Field-South (ADF--S) Survey. The analysis uses a carefully selected sample of 186 sources detected at the 90 $\mu$m AKARI band, identified as galaxies with cross-identification in public catalogues. For sources without known spectroscopic redshifts, we estimate photometric redshifts after a test of two independent methods: one based on using mainly the optical -- mid infrared range, and one based on the whole range of ultraviolet -- far infrared data. We observe a vast improvement in the estimation of photometric redshifts when far infrared data are included, compared with an approach based mainly on the optical -- mid infrared range. We discuss the physical properties of our far-infrared-selected sample. We conclude that this sample consists mostly of rich in dust and young stars nearby galaxies, and, furthermore, that almost 25% of these sources are (Ultra)Luminous Infrared Galaxies. Average SEDs normalized at 90 $\mu$m for normal galaxies (138 sources), LIRGs (30 sources), and ULIRGs (18 galaxies) a the significant shift in the peak wavelength of the dust emission, and an increasing ratio between their bolometric and dust luminosities which varies from 0.39 to 0.73.

Star formation in the luminous YSO IRAS 18345-0641

Aims: We aim to understand the star formation associated with the luminous young stellar object (YSO) IRAS 18345-0641 and to address the complications arising from unresolved multiplicity in interpreting the observations of massive star-forming regions. Methods: New infrared imaging data at sub-arcsec spatial resolution are obtained for IRAS 18345-0641. The new data are used along with mid- and far-IR imaging data, and CO (J=3-2) spectral line maps downloaded from archives to identify the YSO and study the properties of the outflow. Available radiative-transfer models are used to analyze the spectral energy distribution (SED) of the YSO. Results: Previous tentative detection of an outflow in the H_2 (1-0) S1 line (2.122 micron) is confirmed through new and deeper observations. The outflow appears to be associated with a YSO discovered at infrared wavelengths. At high angular resolution, we see that the YSO is probably a binary. The CO (3--2) lines also reveal a well defined outflow. Nevertheless, the direction of the outflow deduced from the H_2 image does not agree with that mapped in CO. In addition, the age of the YSO obtained from the SED analysis is far lower than the dynamical time of the outflow. We conclude that this is probably caused by the contributions from a companion. High-angular-resolution observations at mid-IR through mm wavelengths are required to properly understand the complex picture of the star formation happening in this system, and generally in massive star forming regions, which are located at large distances from us.

Sequential clustering of star formations in IC 1396

We present a comprehensive study of the H II region IC 1396 and its star forming activity, in which multi-wavelength data ranging from the optical to the near- and far-infrared were employed. The surface density distribution of all the 2MASS sources with a certain detection toward IC 1396 indicates the existence of a compact cluster spatially consistent with the position of the exciting source of the H II region, HD 206267. The spatial distribution of the sources with excessive infrared emission, selected based on archived 2MASS data, reveals the existence of four sub-clusters in this region. One is associated with the open cluster Trumpler 37. The other three are found to be spatially coincident with the bright rims of the H II region. All the sources with excessive emission in the near infrared are cross-identified with AKARI IRC data. An analysis of the spectral energy distributions (SEDs) of the resultant sample leads to the identification of eight CLASS I, 15 CLASS II and 15 CLASS III sources in IC 1396. Optical identification of the sample sources with R magnitudes brighter than 17 mag corroborates the results from the SED analysis. Based on the spatial distribution of the infrared young stellar objects at different evolutionary stages, the surrounding sub-clusters located in the bright rims are believed to be younger than the central one. This is consistent with a scenario of sequential star formation in this region. Imaging data of a dark patch in IC 1396 by Herschel SPIRE, on the other hand, indicate the presence of two far-infrared cores in LDN 1111, which are likely to be a new generation of protostellar objects in formation. So we infer that the star formation process in this H II region was not continuous but rather episodic.

GOODS-Herschel: radio-excess signature of hidden AGN activity in distant star-forming galaxies

We present here a new spectral energy distribution (SED) fitting approach that we adopt to select radio-excess sources amongst distant star-forming galaxies in the GOODS-Herschel (North) field and to reveal the presence of hidden, highly obscured AGN. Through extensive SED analysis of 458 galaxies with radio 1.4 GHz and mid-IR 24 um detections using some of the deepest Chandra X-ray, Spitzer and Herschel infrared, and VLA radio data available to date, we have robustly identified a sample of 51 radio-excess AGN (~1300 deg^-2) out to redshift z~3. These radio-excess AGN have a significantly lower far-IR/radio ratio (q<1.68) than the typical relation observed for star-forming galaxies (q~2.2). We find that ~45% of these radio-excess sources have a dominant AGN component in the mid-IR band, while for the remainders the excess radio emission is the only indicator of AGN activity. The fraction of radio-excess AGN increases with X-ray luminosity reaching ~60% at Lx~10^44-10^45 erg/s, making these sources an important part of the total AGN population. However, almost half (24/51) of these radio-excess AGN are not detected in the deep Chandra X-ray data, suggesting that some of these sources might be heavily obscured. We also find that the specific star formation rates (sSFRs) of the radio-excess AGN are on average lower that those observed for X-ray selected AGN hosts, indicating that our sources are forming stars more slowly than typical AGN hosts, and possibly their star formation is progressively quenching.

AHUBBLE SPACE TELESCOPETREASURY STUDY OF STAR-FORMING REGIONS IN THE LOCAL GROUP. II. YOUNG STELLAR POPULATIONS IN M31

We studied the young stellar populations of 22 star-forming regions in the Andromeda galaxy (M31), with Hubble Space Telescope (HST) multi-band imaging from far-UV to I. The regions were selected from Galaxy Evolution Explorer (GALEX) wide-field far-UV imaging; they sample different environments and galactocentric distances from 6 to 22 kpc. They were imaged with 30 HST fields (360 distinct images, in six bandpasses), with a pixel scale of 0.38 pc projected on the sky, at the distance of M31. This study is part of HST treasury survey program HST-GO-11079, which includes star-forming regions in eight Local Group galaxies. We provide a merged catalog of six-band stellar photometry in the 30 M31 fields, containing 118,036 sources brighter than V and B ~ 23 mag. Each HST field covers about 0.3 kpc2 in M31, and contains up to ~7000 stars, of which the number varies by a factor of >7 among the target regions; a large fraction of the sample are hot massive stars, due to our choice of filters and exposures. We derived stellar physical parameters and interstellar extinction for individual sources by spectral energy distribution analysis with model-atmosphere colors, and used the results to infer ages, massive stars content, and extinction of the star-forming regions. Reddening is up to E(B – V) 0.6 mag in some OB associations, and lowest in the outermost regions (average of 0.12 mag in OB184 at 21.9 kpc). We examined the spatial distribution (clustering) of the hot massive stars, and defined OB associations on various spatial scales from compact to wider, more spread out ones. A hierarchical structuring is observed, with small compact groups arranged within large complexes. Their areas vary from less than 10 to 105 pc2, and masses are up to ≈105 M ☉, in the scales sampled by our analysis. Their cumulative mass distribution follows a power law, at least in part of the sampled regime. Hot-star counts in the young regions compare very well with integrated measurements of UV flux from GALEX.

A Positive Influence

Astronomy It is thought that most, if not all, large galaxies host massive black holes at their centers. These black holes become luminous when they accrete matter from their immediate surroundings, and during these phases of active accretion, they are thought to influence the evolution of their galaxies' stellar populations (see Volonteri, Special Issue Review, 3 Aug. 2012, p. [544][1]). Whether this influence is positive or negative (i.e., whether black hole accretion suppresses or induces star formation) is a matter of current debate. By combining observations obtained with the European Space Agency's Herschel Space Observatory with data from NASA's Spitzer Space Telescope, Barthel et al. inferred the star formation rates for three extremely massive galaxies, whose central massive black holes are being fed at a high rate. Analysis of these galaxies' spectral energy distributions implies that they are all forming stars at exceptionally high rates, showing that powerful black hole accretion and vigorous star formation can occur simultaneously, at least in the most massive galaxies. Astrophys. J. 757 , L26 (2012). [1]: /lookup/doi/10.1126/science.1220843

Bayesian analysis of galaxy spectral energy distributions with BayeSED

AbstractIn Han &amp; Han (2012), we have preliminarily built BayeSED and applied it to a sample of hyperluminous infrared galaxies. The physically reasonable results obtained from Bayesian model comparison and parameter estimation show that BayeSED could be a useful tool for understanding the nature of complex systems, such as dust obscured starburst-AGN composite galaxies, from decoding their complex SEDs. In this contribution, we present a more rigorous test of BayeSED by making a mock catalog from model SEDs with the value of all parameters to be known in advance.