Source author record

R. Siebenmorgen

R. Siebenmorgen 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
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

19 published item(s)

preprint2022arXiv

Broadened profiles of Diffuse Interstellar Bands

The paper describes profile broadening and peak wavelength variation of diffuse interstellar bands (DIBs) measured for 46 lines of sight, probably caused by physical properties of intervening clouds. Full width at half maximum of four studied diffuse bands (5780, 5797, 6196 and 6614 Å) demonstrate strong variability sometimes doubling the features' width. Despite the high magnitude of the effect, our current analysis is restricted to the strongest diffuse bands because the weaker ones require a much higher S/N ratio. The profile broadening in the studied DIBs moves the profile's centers towards longer wavelengths, probably due to the excitation of higher levels of the P branch of the unknown molecular carrier. Moreover, diffuse bands are broader in clouds with abundantly populated vibrationally excited states of the hydrogen molecules, i.e. DIB's broadening correlates with the rotational temperature estimated on H$_2$ $ν$=2 vibrational level. However, objects demonstrating extremely broadened profiles of DIBs are scarce. The extreme peculiarity of DIBs' profiles was detected in Herschel~36. Here we show the gradual growths of the widths of diffuse bands, confirmed in spectra from different instruments.

preprint2021arXiv

FORS-Up: Making the most versatile instrument in Paranal ready for 15 more years of operations

The FORS Upgrade project (FORS-Up) aims at bringing a new life to the highly demanded workhorse instrument attached to ESO's Very Large Telescope (VLT). FORS2 is a multimode optical instrument, which started regular science operations in 2000 and since then, together with its twin, FORS1, has been one of the most demanded and most productive instruments of the VLT. In order to ensure that a FORS shall remain operational for at least another 15 years, an upgrade has been planned. This is required as FORS2 is using technology and software that is now obsolete and cannot be put and maintained to the standards in use at the Observatory. The project - carried out as a collaboration between ESO and INAF-Astronomical Observatory of Trieste - aims at bringing to the telescope in 2023/2024 a refurbished instrument with a new scientific detector, an upgrade of the instrument control software and electronics, a new calibration unit, as well as additional filters and grisms. The new FORS will also serve as a test bench for the Extremely Large Telescope (ELT) standard technologies (among them the use of programmable logic controllers and of the features of the ELT Control Software). The project aims at minimising the downtime of the instrument by performing the upgrade on the currently decommissioned instrument FORS1 and retrofitting the Mask Exchange Unit and polarisation optics from FORS2 to FORS1.

preprint2021arXiv

Mid-infrared photometry of the T Tauri triple system with kernel phase interferometry

T Tauri has long been the prototypical young pre-main-sequence star. However, it has now been decomposed into a triple system with a complex disk and outflow geometry. We aim to measure the brightness of all three components of the T Tauri system (T Tau N, T Tau Sa, T Tau Sb) in the mid-infrared in order to obtain photometry around the $\sim 9.7~μm$ silicate feature. This allows us to study their variability and to investigate the distribution of dust and the geometry of circumstellar and circumbinary disks in this complex system. We observe T Tauri with the VLT/VISIR-NEAR instrument. With kernel phase interferometry post-processing of the data, and using the astrometric positions of all three components from VLT/SPHERE, we measure the three components' individual brightnesses (including the southern binary at an angular separation down to $\sim 0.2~λ/D$) and obtain their photometry. In order to validate our methods, we simulate and recover mock data of the T Tauri system using the observed reference point-spread function of HD 27639. We find that T Tau N is rather stable and shows weak silicate emission, while T Tau Sa is highly variable and shows prominent silicate absorption. T Tau Sb became significantly fainter compared to data from 2004 and 2006, suggesting increased extinction by dust. The precision of our photometry is limited by systematic errors, which is consistent with previous studies using kernel phase interferometry. Our results confirm the complex scenario of misaligned disks in the T Tauri system that had been observed previously, and they are in agreement with the recently observed dimming of T Tau Sb in the near-infrared. Our mid-infrared photometry supports the interpretation that T Tau Sb has moved behind the dense region of the Sa-Sb circumbinary disk on its tight orbit around Sa, therefore suffering increased extinction.

preprint2020arXiv

{\sc precision}: A fast python pipeline for high-contrast imaging -- application to SPHERE observations of the red supergiant VX Sagitariae

The search for extrasolar planets has driven rapid advances in instrumentation, resulting in cameras such as SPHERE at the VLT, GPI at Gemini South and SCExAO at Subaru, capable of achieving very high contrast ($\sim10^{6}$) around bright stars with small inner working angles ($\sim 0\farcs{1}$). The optimal exploitation of data from these instruments depends on the availability of easy-to-use software to process and analyse their data products. We present a pure-python pipeline, {\sc precision}, which provides fast, memory-efficient reduction of data from the SPHERE/IRDIS near-infrared imager, and can be readily extended to other instruments. We apply {\sc precision} to observations of the extreme red supergiant VX~Sgr, the inner outflow of which is revealed to host complex, asymmetric structure in the near-IR. In addition, optical polarimetric imaging reveals clear extended polarised emission on $\sim0.5^{\prime\prime}$ scales which varies significantly with azimuth, confirming the asymmetry. While not conclusive, this could suggest that the ejecta are confined to a disc or torus, which we are viewing nearly face on, although other non-spherical or clumpy configurations remain possible. VX~Sgr has no known companions, making such a geometry difficult to explain, as there is no obvious source of angular momentum in the system.

preprint2020arXiv

Dark Dust and single-cloud sightlines in the ISM

The precise characteristics of clouds and the nature of dust in the diffuse interstellar medium can only be extracted by inspecting the rare cases of single-cloud sightlines. In our nomenclature such objects are identified by interstellar lines, such as K I, that show at a resolving power of $λ/Δλ\sim 75,000$ one dominating Doppler component that accounts for more than half of the observed column density. We searched for such sightlines using high-resolution spectroscopy towards reddened OB stars for which far-UV extinction curves are known. We compiled a sample of 186 spectra, 100 of which were obtained specifically for this project with UVES. In our sample we identified 65 single-cloud sightlines, about half of which were previously unknown. We used the CH/CH$^+$ line ratio of our targets to establish whether the sightlines are dominated by warm or cold clouds. We found that CN is detected in all cold (CH/CH$^+ >1$) clouds, but is frequently absent in warm clouds. We inspected the WISE ($3-22\, μ$m) observed emission morphology around our sightlines and excluded a circumstellar nature for the observed dust extinction. We found that most sightlines are dominated by cold clouds that are located far away from the heating source. For 132 stars, we derived the spectral type and the associated spectral type-luminosity distance. We also applied the interstellar Ca II distance scale, and compared these two distance estimates with GAIA parallaxes. These distance estimates scatter by 40%. By comparing spectral type-luminosity distances with those of GAIA, we detected a hidden dust component that amounts to a few mag of extinction for eight sightlines. This Dark Dust is populated by $\ge 1 μ$m large grains and predominately appears in the field of the cold interstellar medium.

preprint2019arXiv

RefPlanets: Search for reflected light from extra-solar planets with SPHERE/ZIMPOL

RefPlanets is a guaranteed time observation (GTO) programme that uses the Zurich IMaging POLarimeter (ZIMPOL) of SPHERE/VLT for a blind search for exoplanets in wavelengths from 600-900 nm. The goals of this study are the characterization of the unprecedented high polarimetic contrast and polarimetric precision capabilities of ZIMPOL for bright targets, the search for polarized reflected light around some of the closest bright stars to the Sun and potentially the direct detection of an evolved cold exoplanet for the first time. For our observations of Alpha Cen A and B, Sirius A, Altair, Eps Eri and Tau Ceti we used the polarimetric differential imaging (PDI) mode of ZIMPOL which removes the speckle noise down to the photon noise limit for angular separations >0.6". We describe some of the instrumental effects that dominate the noise for smaller separations and explain how to remove these additional noise effects in post-processing. We then combine PDI with angular differential imaging (ADI) as a final layer of post-processing to further improve the contrast limits of our data at these separations. For good observing conditions we achieve polarimetric contrast limits of 15.0-16.3 mag at the effective inner working angle of about 0.13", 16.3-18.3 mag at 0.5" and 18.8-20.4 mag at 1.5". The contrast limits closer in (<0.6") depend significantly on the observing conditions, while in the photon noise dominated regime (>0.6"), the limits mainly depend on the brightness of the star and the total integration time. We compare our results with contrast limits from other surveys and review the exoplanet detection limits obtained with different detection methods. For all our targets we achieve unprecedented contrast limits. Despite the high polarimetric contrasts we are not able to find any additional companions or extended polarized light sources in the data that has been taken so far.

preprint2016arXiv

Appearance of Dusty Filaments at Different Viewing Angles

Context: In the last years, there have been many studies on the omnipresence and structures of filaments in star-forming regions, as well as their role in the process of star formation. Those filaments are normally identified as elongated fibres across the plane of the sky. But how would we detect filaments that are inclined? Aims: We aim to learn more about whether, and how, total column density or dust temperature change with respect to the line of sight. Such variations would enable observers to use dust observations to identify and study filaments at any inclination and gain more insight on the distribution and orientations of filaments within the Galactic plane. Methods: As a first step, we perform numerical calculations on simple cylindrical models to evaluate the influence of filament geometry on the average flux density. After that, we apply our three-dimensional Monte Carlo dust radiative transfer code on two models of star-forming regions and derive maps of effective total column density and dust temperature at different viewing angles. Results: We see only slight changes of average flux density for all cylinders we study. For our more complex models, we find that the effective dust temperature is not sensitive to viewing angle, while the total column density is strongly influenced, with differences exceeding an order of magnitude. The variations are not injective with the viewing angle and depend on the structure of the object. Conclusions: We conclude that there is no single quantity in our analysis that can uniquely trace the inclination and three-dimensional structure of a filament based on dust observations alone. However, observing wide variations in total column density at a given effective dust temperature is indicative of inclined filaments.

preprint2016arXiv

Infrared polarimetry of Mrk 231: Scattering off hot dust grains in the central core

We present high-angular (0.17$-$0.35 arcsec) resolution imaging polarimetric observations of Mrk 231 in the 3.1 $μ$m filter using MMT-Pol on the 6.5-m MMT, and in the 8.7 $μ$m, 10.3 $μ$m, and 11.6 $μ$m filters using CanariCam on the 10.4-m Gran Telescopio CANARIAS. In combination with already published observations, we compile the 1$-$12 $μ$m total and polarized nuclear spectral energy distribution (SED). The total flux SED in the central 400 pc is explained as the combination of 1) a hot (731 $\pm$ 4 K) dusty structure, directly irradiated by the central engine, which is at 1.6 $\pm$ 0.1 pc away and attributed to be in the pc-scale polar region, 2) an optically-thick, smooth and disk-like dusty structure (`torus') with an inclination of 48 $\pm$ 23$^{\circ}$ surrounding the central engine, and 3) an extinguished (A$_{\mbox{V}} =$ 36 $\pm$ 5 mag) starburst component. The polarized SED decreases from 0.77 $\pm$ 0.14 per cent at 1.2 $μ$m to 0.31 $\pm$ 0.15 per cent at 11.6 $μ$m and follows a power-law function, $λ^{\sim0.57}$. The polarization angle remains constant ($\sim$108$^{\circ}$) in the 1$-$12 $μ$m wavelength range. The dominant polarization mechanism is explained as scattering off hot dust grains in the pc-scale polar regions.

preprint2015arXiv

Extinction and dust properties in a clumpy medium

(abridged) The dust content of the universe is primarily explored via its interaction with stellar photons, producing interstellar extinction. However, owing to the physical extension of the observing beam, observations may detect scattered photons, resulting in a change in the observed (or effective) extinction, depending on the spatial distribution of the dust and the resolution of the instrument. We investigate the influence of clumpy dust distributions on effective extinction toward embedded sources and those in the diffuse ISM. We use Monte Carlo radiative transfer to examine effective extinction for various geometries. By varying the number, optical depth and volume-filling factor of clumps in models of spherical shells and the diffuse ISM, we explore the evolution of extinction. Depending on the number of scattering events in the beam, the extinction curve steepens in homogeneous media and flattens in clumpy media. As a result, clumpy dust distributions can to reproduce extinction curves with arbitrary R_v, the effective ratio of total-to-selective extinction. The flattening `washes out' the 2175 Åbump and shifts the peak to shorter wavelengths. The mean R_v of a shell correlates with the optical depth of an individual clump and the wavelength at which a clump becomes optically thick. Similar behaviour is seen for edge-on discs or tori. However, at grazing inclinations the combination of extinction and strong forward scattering results in chaotic behaviour. Caution is therefore advised when measuring extinction in, for example, AGN tori or toward SNIa or GRB afterglows. In face-on discs, the shape of the scattered continuum changes with clumpiness, however, unlike absorption features, individual features in the scattering cross-sections are preserved. Finally, we show that diffuse interstellar extinction is not modified by scattering on scales of a few kpc.

preprint2015arXiv

Large dust grains in the wind of VY Canis Majoris

Massive stars live short lives, losing large amounts of mass through their stellar wind. Their mass is a key factor determining how and when they explode as supernovae, enriching the interstellar medium with heavy elements and dust. During the red supergiant phase, mass-loss rates increase prodigiously, but the driving mechanism has proven elusive. Here we present high-contrast optical polarimetric-imaging observations of the extreme red supergiant VY Canis Majoris and its clumpy, dusty, mass-loss envelope, using the new extreme-adaptive-optics instrument SPHERE at the VLT. These observations allow us to make the first direct and unambiguous detection of submicron dust grains in the ejecta; we derive an average grain radius $\sim$ 0.5 $μ$m, 50 times larger than in the diffuse ISM, large enough to receive significant radiation pressure by photon scattering. We find evidence for varying grain sizes throughout the ejecta, highlighting the dynamical nature of the envelope. Grains with 0.5 $μ$m sizes are likely to reach a safe distance from the eventual explosion of VY Canis Majoris; hence it may inject upwards of 10$^{-2}$ M$_\odot$ of dust into the ISM.

preprint2015arXiv

Star formation in z>1 3CR host galaxies as seen by Herschel

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.

preprint2014arXiv

Nuclear 11.3$μ$m PAH emission in local active galactic nuclei

We present Gran Telescopio CANARIAS CanariCam 8.7$μ$m imaging and 7.5-13$μ$m spectroscopy of six local systems known to host an active galactic nucleus (AGN) and have nuclear star formation. Our main goal is to investigate whether the molecules responsible for the 11.3$μ$m polyclyclic aromatic hydrocarbon (PAH) feature are destroyed in the close vicinity of an AGN. We detect 11.3$μ$m PAH feature emission in the nuclear regions of the galaxies as well as extended PAH emission over a few hundred parsecs. The equivalent width (EW) of the feature shows a minimum at the nucleus but increases with increasing radial distances, reaching typical star-forming values a few hundred parsecs away from the nucleus. The reduced nuclear EW are interpreted as due to increased dilution from the AGN continuum rather than destruction of the PAH molecules. We conclude that at least those molecules responsible for the 11.3$μ$m PAH feature survive in the nuclear environments as close as 10pc from the AGN and for Seyfert-like AGN luminosities. We propose that material in the dusty tori, nuclear gas disks, and/or host galaxies of AGN is likely to provide the column densities necessary to protect the PAH molecules from the AGN radiation field.

preprint2013arXiv

Active galactic nucleus torus models and the puzzling infrared spectrum of IRAS F10214+4724

We present a revised model for the infrared emission of the 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 present and discuss a grid of smooth active galactic nucleus (AGN) torus models computed with the method of Efstathiou & Rowan-Robinson and demonstrate that the combination of these models and the starburst models of Efstathiou and coworkers, while able to give an excellent fit to the average spectrum of Seyfert 2s and spectra of individual type 2 quasars 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 that is usually indicative of heavy absorption by dust) but shows a silicate feature in emission. Such emission features are not expected in sources with optical/near-infrared type 2 AGN spectral signatures. The Herschel data show that there is more power emitted in the rest-frame 20-50um wavelength range compared with the model presented by Efstathiou which assumes three components of emission: an edge-on torus, clouds (at a temperature of 610 and 200K) that are associated with the narrow-line region (NLR) and a highly obscured starburst that dominates in the submillimetre. We present a revised version of that model that assumes an additional component of emission which we associate with NLR clouds at a temperature of 100K. The 100K dust component could also be explained by a highly obscured hot starburst. The model suggests that the NLR of IRAS F10214+4724 has an unusually high covering factor (>= 17%) or more likely the magnification of the emission from the NLR clouds is significantly higher than that of the emission from the torus.

preprint2013arXiv

Dust in the diffuse interstellar medium: Extinction, emission, linear and circular polarisation

We present a model for the diffuse interstellar dust that explains the observed wavelength-dependence of extinction, emission, linear and circular polarisation of light. The model is set-up with a small number of parameters. It consists of a mixture of amorphous carbon and silicate grains with sizes from the molecular domain of 0.5 up to about 500nm. Dust grains with radii larger than 6nm are spheroids. Spheroidal dust particles have a factor 1.5 - 3 larger absorption cross section in the far IR than spherical grains of the same volume. Mass estimates derived from submillimeter observations that ignore this effect are overestimated by the same amount. In the presence of a magnetic field, spheroids may be partly aligned and polarise light. We find that polarisation spectra help to determine the upper particle radius of the otherwise rather unconstrained dust size distribution. Stochastically heated small grains of graphite, silicates and polycyclic aromatic hydrocarbons (PAHs) are included. We tabulate parameters for PAH emission bands in various environments. They show a trend with the hardness of the radiation field that can be explained by the ionisation state or hydrogenation coverage of the molecules. For each dust component its relative weight is specified, so that absolute element abundances are not direct input parameters. The model is confronted with the average properties of the Milky Way, which seems to represent dust in the solar neighbourhood. It is then applied to four specific sight lines including the reflection nebula NGC2023. For these sight lines, we present linear and circular spectro-polarimetric observations obtained with FORS/VLT. Using prolate rather than oblate grains gives a better fit to observed spectra; the axial ratio of the spheroids is typically two and aligned silicates are the dominant contributor to the polarisation.

preprint2010arXiv

Mid-infrared interferometry of the massive young stellar object NGC3603 - IRS 9A

We present observations and models for one of these MYSO candidates, NGC3603 IRS 9A. Our goal is to investigate with infrared interferometry the structure of IRS 9A on scales as small as 200AU, exploiting the fact that a cluster of O and B stars has blown away much of the obscuring foreground dust and gas. Observations in the N-band were carried out with the MIDI beam combiner attached to the VLTI. Additional interferometric observations which probe the structure of IRS 9A on larger scales were performed with an aperture mask installed in the T-ReCS instrument of Gemini South. The spectral energy distribution (SED) is constrained by the MIDI N-band spectrum and by data from the Spitzer Space Telescope. Our efforts to model the structure and SED of IRS 9A range from simple geometrical models of the brightness distribution to one- and two-dimensional radiative transfer computations. The target is resolved by T-ReCS, with an equivalent (elliptical) Gaussian width of 330mas by 280mas (2300 AU by 2000 AU). Despite this fact, a warm compact unresolved component was detected by MIDI which is possibly associated with the inner regions of a flattened dust distribution. Based on our interferometric data, no sign of multiplicity was found on scales between about 200AU and 700AU projected separation. A geometric model consisting of a warm (1000 K) ring (400 AU diameter) and a cool (140 K) large envelope provides a good fit to the data. No single model fitting all visibility and photometric data could be found, with disk models performing better than spherical models. While the data are clearly inconsistent with a spherical dust distribution they are insufficient to prove the existence of a disk but rather hint at a more complex dust distribution.

preprint2010arXiv

Mid-infrared spectroscopy of high-redshift 3CRR sources

Using the Spitzer Space Telescope, we have obtained rest frame 9-16mu spectra of 11 quasars and 9 radio galaxies from the 3CRR catalog at redshifts 1.0<z<1.4. This complete flux-limited 178MHz-selected sample is unbiased with respect to orientation and therefore suited to study orientation-dependent effects in the most powerful active galactic nuclei (AGN). The mean radio galaxy spectrum shows a clear silicate absorption feature (tau_9.7mu = 1.1) whereas the mean quasar spectrum shows silicates in emission. The mean radio galaxy spectrum matches a dust-absorbed mean quasar spectrum in both shape and overall flux level. The data for individual objects conform to these results. The trend of the silicate depth to increase with decreasing core fraction of the radio source further supports that for this sample, orientation is the main driver for the difference between radio galaxies and quasars, as predicted by AGN unification. However, comparing our high-z sample with lower redshift 3CRR objects reveals that the absorption of the high-z radio galaxy MIR continuum is lower than expected from a scaled up version of lower luminosity sources, and we discuss some effects that may explain these trends.

preprint2010arXiv

Synthesising, using, and correcting for telluric features in high-resolution astronomical spectra

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.

preprint2009arXiv

VLT diffraction-limited imaging at 11 and 18 micron of the nearest active galactic nuclei

Mid-infrared imaging at resolutions of 300 mas of the central kpc region of 13 nearby, well-known active galaxies is presented. The bulk of the mid-IR emission is concentrated on an unresolved central source within a size of less than 5 to 130 pc, depending on the object distance. Further resolved emission is detected in 70% of the sample in the form of circumnuclear star-forming rings or diffuse nuclear extended emission. In the three cases with circumnuclear star formation, the stellar contribution is at least as important as that of the AGN. In those with extended nuclear emission -- a third of the sample -- this emission represents a few per cent of the total measured; however, this contribution may be underestimated because of the chopped nature of these observations. This extended emission is generally collimated in a preferential direction often coinciding with that of the extended ionized gas or the jet. In all cases, the nuclear fluxes measured at 11.8 and 18.7 micron represent a minor contribution of the flux levels measured by large aperture IRAS data at the nearest energy bands of 12 and 25 micron. This contribution ranges from 30% to less than 10%. In only three cases do the AGN fluxes agree with IRAS to within a factor of 2. In the AGNs with strong circumnuclear star formation, this component can well account for most of the IRAS flux measured in these objects. But in all other cases, either a low surface brightness component extending over galactic scales or strong extra-nuclear IR sources -- e.g. HII regions in spiral arms -- have to be the main source of the IRAS emission. In either case, the contribution of these components dwarfs that of the AGN at mid-IR wavelengths.

preprint1994arXiv

BD+30 3639: The Infrared Spectrum During Post-AGB Stellar Evolution

We present a radiative-transfer calculation which reproduces the infrared spectrum of the planetary nebula BD~+30$^{\circ}$3639. We calculate the transfer process through absorption and scattering in a spherical-symmetric multi-grain dust shell. The emission of transiently heated particles is taken into account, as well as polycyclic aromatic hydrocarbons. We obtain an acceptable fit to most of the spectrum, including the PAH infrared bands. At submillimetre wavelengths the observed emission is larger than the model predicts, indicating that large dust conglomerates (``f{}luffy grains'') may be needed as an additional constituent. The fit favours a distance of $\ge 2 \,$kpc, which implies that BD~+30$^\circ$3639 has evolved from a massive progenitor of several solar masses. A low dust-to-gas mass ratio is found in the ionised region. The calculations yield an original mass-loss rate of $2\times10^{-4} \msolar \peryr$ on the Asymptotic Giant Branch. Using this mass-loss rate, we calculate how the infrared spectrum has evolved during the post-AGB evolution. We show in particular the evolution of the IRAS colours during the preceding post-AGB evolution.