Ralf Siebenmorgen

Massive star formation in galaxies: radiative transfer models of the UV to millimetre emission of starburst galaxies

We present illustrative models for the UV to millimeter emission of starburst galaxies which are treated as an ensemble of optically thick giant molecular clouds (GMCs) centrally illuminated by recently formed stars. The models follow the evolution of the GMCs due to the ionization-induced expansion of the HII regions and the evolution of the stellar population within the GMC according to the Bruzual&Charlot stellar population synthesis models. The effect of transiently heated dust grains/PAHs to the radiative transfer, as well as multiple scattering, is taken into account. The expansion of the HII regions and the formation of a narrow neutral shell naturally explains why the emission from PAHs dominates over that from hot dust in the near to mid-IR, an emerging characteristic of the infrared spectra of starburst galaxies. The models allow us to relate the observed properties of a galaxy to its age and star formation history. We find that exponentially decaying 10^7-10^8 yrs old bursts can explain the IRAS colours of starburst galaxies. The models are also shown to account satisfactorily for the multiwavelength data on the prototypical starburst galaxy M82 and NGC6090, a starburst galaxy recently observed by ISO. In M82 we find evidence for two bursts separated by 10^7yrs. In NGC6090 we find that at least part of the far-IR excess may be due to the age of the burst (6.4 x 10^7yrs). We also make predictions about the evolution of the luminosity of starbursts at different wavelengths which indicate that far-IR surveys may preferentially detect older starbursts than mid-IR surveys.

CRIRES: A High Resolution Infrared Spectrograph for ESO’s VLT

CRIRES is a cryogenic, pre-dispersed, infrared echelle spectrograph designed to provide a resolving power lambda/(Delta lambda) of 105 between 1 and 5mu m at the Nasmyth focus B of the 8m VLT unit telescope #1 (Antu). A curvature sensing adaptive optics system feed is used to minimize slit losses and to provide diffraction limited spatial resolution along the slit. A mosaic of 4 Aladdin~III InSb-arrays packaged on custom-fabricated ceramics boards has been developed. This provides for an effective 4096x512 pixel focal plane array, to maximize the free spectral range covered in each exposure. Insertion of gas cells to measure high precision radial velocities is foreseen. For measurement of circular polarization a Fresnel rhomb in combination with a Wollaston prism for magnetic Doppler imaging is foreseen. The implementation of full spectropolarimetry is under study. This is one result of a scientific workshop held at ESO in late 2003 to refine the science-case of CRIRES. Installation at the VLT is scheduled during the first half of 2005. Here we briefly recall the major design features of CRIRES and describe its current development status including a report of laboratory testing.

Dust in starburst nuclei and ULIRGs: SED models for observers

Aims. We provide a library of some 7000 SEDs for the nuclei of starburst and ultra luminous galaxies. Its purpose is to quickly ob tain estimates of the basic parameters, such as luminosity, size and dust or gas mass and to predict the flux at yet unobserved wavelengths . The procedure is simple and consists of finding an element in the library that m atches the observations. The objects may be in the local universe or at high z. Methods. We calculate the radiative transfer in spherical symmetry for a stellar cluster permeated by an interstellar medium with standard (Milky Way) dust properties. The cluster contains two stellar populations: old bulge stars and OB stars. Because the latter are young, a certain fraction of them will be embedded in compact clouds which constitute hot spots that determine the MIR fluxes. Results. We present SEDs for a broad range of luminosities, sizes and obscurations. We argue that the assumption of spherical symmetry and the neglect of clumpiness of the medium are not severe shortcomings for computing the dust emission. The validity of the approach is demonstrated by matching the SED of seven of the best studied galaxies, including M82 and Arp220, by library elements. In all cases, one finds an element which fits the observed SED very well, and the p arameters defining the element are in full accord with what is known about the galaxy from detailed studies. We also compare our method of computing SEDs with other techniques described in the literature.

Mid–infrared emission of galactic nuclei - TIMMI2 versus ISO observations and models

We investigate the mid-infrared radiation of galaxies that are powered by a starburst or by an AGN. For this end, we compare the spectra obtained at different spatial scales in a sample of infrared bright galaxies. ISO observations which include emission of the nucleus as well as most of the host galaxy are compared with TIMMI2 spectra of the nuclear region. We find that ISO spectra are generally dominated by strong PAH bands. However, this is no longer true when inspecting the mid-infrared emission of the pure nucleus. Here PAH emission is detected in starbursts whereas it is significantly reduced or completely absent in AGNs. A physical explanation of these new observational results is presented by examining the temperature fluctuation of a PAH after interaction with a photon. It turns out that the hardness of the radiation field is a key parameter for quantifying the photo-destruction of small grains. Our theoretical study predicts PAH evaporation in soft X-ray environments. Radiative transfer calculations of clumpy starbursts and AGN corroborate the observational fact that PAH emission is connected to starburst activity whereas PAHs are destroyed near an AGN. The radiative transfer models predict for starbursts a much larger mid-infrared size than for AGN. This is confirmed by our TIMMI2 acquisition images: We find that the mid-infrared emission of Seyferts is dominated by a compact core while most of the starbursts are spatially resolved.

ISOCAM survey and dust models of 3CR radio galaxies and quasars

We present a survey of all 3CR sources imaged with ISOCAM onboard the Infrared Space Observatory (ISO). For the source, we present spatially integrated mid--infrared (MIR, 5 - 18mic.) fluxes measured from newly calibrated ISOCAM images. In total, we detected 68 AGN of the 3CR catalogue, at redshifts z<2.5. The one with the highest redshift is 4C+72.26 at z = 3.53. ISOCAM data are combined with other photometric measurements to construct the spectral energy distribution (SED) from optical to radio wavelengths. In order to describe dust emission we apply new radiative transfer models. By varying three parameters, luminosity, effective size and extinction, we obtain a fit to the SED for our objects. In the MIR the hot dust component is mainly due to small grains and PAHs. In the models, a type 1 AGN is represented by a compact dust distribution, the dust is therefore very warm and emission of PAHs is weak because of photo--destruction. In AGNs of type 2, the dust is relatively colder but PAH bands are strong.

Subarcsecond mid-infrared imaging of warm dust in the narrow-line region of NGC 1068

Subarcsecond 8 and 10 μm and diffraction-limited 19 μm imaging of the inner few hundred parsecs of the Seyfert nucleus in NGC 1068 shows the emission to be extended over a region of ∼70×140 pc. In particular, 10.3 μm images with spatial resolutions of 0″.5 or better reveal that the warm dust is associated with the narrow-line clouds and is probably partially mixed with the photoionized gas. Extinction considerations, however, imply that the bulk of the warm dust is located deeper in neighboring molecular clouds, the exposed surfaces of which form the narrow-line clouds

Mid-infrared spectral evidence for a luminous dust enshrouded source in Arp 220

We have re-analyzed the 6-12 µm ISO spectrum of the ultra-luminous infrared galaxy Arp 220 with the conclusion that it is not consistent with that of a scaled up version of a typical starburst. Instead, both template fitting with spectra of the galaxies NGC 4418 and M 83 and with dust models suggest that it is best represented by combinations of a typical starburst component, exhibiting PAH emission features, and a heavily absorbed dust continuum which contributes ∼40% of the 6-12 µm flux and likely dominates the luminosity. Of particular significance relative to previous studies of Arp 220 is the fact that the emission feature at 7.7 µm comprises both PAH emission and a broader component resulting from ice and silicate absorption against a heavily absorbed continuum. Extinction to the PAH emitting source, however, appears to be relatively low. We ten- tatively associate the PAH emitting and heavily dust/ice absorbed components with the diffuse emission region and the two compact nuclei respectively identified by Soifer et al. (2002) in their higher spatial resolution 10 µm study. Both the similarity of the absorbed continuum with that of the embedded Galactic protostars and results of the dust models imply that the embedded source(s) in Arp 220 could be powered by, albeit extremely dense, starburst activity. Due to the high extinction, it is not possible with the available data to exclude that AGN(s) also contribute some or all of the observed luminosity. In this case, however, the upper limit measured for its hard X-ray emission would require Arp 220 to be the most highly obscured AGN known.

Disappearance of stellar debris disks around main-sequence stars after 400 million years

Almost 5 billion years ago, the Sun formed in a local contraction of a cloud of molecular gas. A rotating disk of gas and dust is believed to have fed material onto the proto-Sun for the first few million years of its life, and to have formed the planets, comets and other Solar System objects. Similar disks, but with less mass, have been observed around a few main-sequence stars such as Vega. The dust particles orbiting stars like Vega will be removed on timescales of the order of 1 Myr (Vega is about 350 Myr old), and therefore must be resupplied, at least for a time. But earlier surveys lacked the sensitivity to determine how many nearby stars have dust disks, and to investigate how long such disks survive. Here we report infrared observations indicating that most stars younger than 300 Myr have dust disks, while most older than 400 Myr do not: ninety per cent of the disks disappear when the star is between 300 and 400 Myr old. Several events that are related to the ‘clean up’ of debris in the early history of our Solar System have a similar timescale.

Synthesising, using, and correcting for telluric features in high-resolution astronomical spectra - A near-infrared case study using CRIRES

We present a technique to synthesise telluric absorption and emission features both for in-situ wavelength calibration and for their removal from astronomical spectra. While the presented technique is applicable for a wide variety of optical and infrared spectra, we concentrate in this paper on selected high-resolution near-infrared spectra obtained with the CRIRES spectrograph to demonstrate its performance and limitation. We find that synthetic spectra reproduce telluric absorption features to about 2%, even close to saturated line cores. Thus, synthetic telluric spectra could be used to replace the observation of telluric standard stars, saving valuable observing time. This technique also provides a precise in-situ wavelength calibration, especially useful for high-resolution near-infrared spectra in the absence of other calibration sources.

First light of the VLT planet finder SPHERE I. Detection and characterization of the substellar companion GJ 758 B

GJ 758 B is a brown dwarf companion to a nearby (15.76%) solar-type, metal-rich (M/H = +0.2 dex) main-sequence star (G9V) that was discovered with Subaru/HiCIAO in 2009. From previous studies, it has drawn attention as being the coldest (similar to 600 K) companion ever directly imaged around a neighboring star. We present new high-contrast data obtained during the commissioning of the SPHERE instrument at the Very Large Telescope (VLT). The data was obtained in Y-, J-, H-, and K-s-bands with the dual-band imaging (DBI) mode of IRDIS, thus providing a broad coverage of the full near-infrared (near-IR) range at higher contrast and better spectral sampling than previously reported. In this new set of high-quality data, we report the re-detection of the companion, as well as the first detection of a new candidate closer-in to the star. We use the new eight photometric points for an extended comparison of GJ 758 B with empirical objects and four families of atmospheric models. From comparison to empirical object, we estimate a T8 spectral type, but none of the comparison objects can accurately represent the observed near-IR fluxes of GJ 758 B. From comparison to atmospheric models, we attribute a T-eff = 600 +/- 100 K, but we find that no atmospheric model can adequately fit all the fluxes of GJ 758 B. The lack of exploration of metal enrichment in model grids appears as a major limitation that prevents an accurate estimation of the companion physical parameters. The photometry of the new candidate companion is broadly consistent with L-type objects, but a second epoch with improved photometry is necessary to clarify its status. The new astrometry of GJ 758 B shows a significant proper motion since the last epoch. We use this result to improve the determination of the orbital characteristics using two fitting approaches: Least-Squares Monte Carlo and Markov chain Monte Carlo. We confirm the high-eccentricity of the orbit (peak at 0.5), and find a most likely semi-major axis of 46.05 AU. We also use our imaging data, as well as archival radial velocity data, to reject the possibility that this is a false positive effect created by an unseen, closer-in, companion. Finally, we analyze the sensitivity of our data to additional closer-in companions and reject the possibility of other massive brown dwarf companions down to 4-5 AU.