BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

P. A. M. van Hoof

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

ResearcherUnclaimed source record
astro-ph.SRastro-ph.GAastro-ph.COphysics.atom-phastro-ph.IMastro-phastro-ph.HEphysics.acc-phphysics.plasm-phphysics.space-ph

Catalog footprint

What is connected

35works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

35 published item(s)

preprint2020arXiv

Current and future development of the photoionization code Cloudy

The gas that is present in the interstellar medium is usually very far removed from (local) thermodynamic equilibrium, and in some cases may also not be in a steady-state equilibrium with its surroundings. The physics of this material is complex and one needs a sophisticated numerical code to study it. For this purpose the open-source photoionization code Cloudy was created. It models the physical state of the gas and predicts the spectrum that it emits. Cloudy is continually being developed to improve the treatment of the microphysical processes and the database of fundamental data that it uses. In this paper we will discuss how we are developing the code to improve our high-density predictions by implementing better collisional-radiative models for all ions. We will also briefly discuss the experimental mode in Cloudy to model gas that is not in steady-state equilibrium and present a preliminary model of recombining gas in a planetary nebula that is on the cooling track. We finish with a short discussion of how we are speeding up the code by using parallelization.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Improved Fe II emission line models for AGN using new atomic datasets

Understanding the Fe II emission from Active Galactic Nuclei (AGN) has been a grand challenge for many decades. The rewards from understanding the AGN spectra would be immense, involving both quasar classification schemes such as "Eigenvector 1" and tracing the chemical evolution of the cosmos. Recently, three large Fe II atomic datasets with radiative and electron collisional rates have become available. We have incorporated these into the spectral synthesis code Cloudy and examine predictions using a new generation of AGN Spectral Energy Distribution (SED), which indicates that the UV emission can be quite different depending on the dataset utilized. The Smyth et al dataset better reproduces the observed Fe II template of the I ZW 1 Seyfert galaxy in the UV and optical regions, and we adopt these data. We consider both thermal and microturbulent clouds and show that a microturbulence of $\approx$ 100 km/s reproduces the observed shape and strength of the so-called Fe II "UV bump". Comparing our predictions with the observed Fe II template, we derive a typical cloud density of $10^{11}$ cm$^{-3}$ and photon flux of $10^{20}$ cm$^{-2}$ s$^{-1}$, and show that these largely reproduce the observed Fe II emission in the UV and optical. We calculate the $I$(Fe II)/$I$(Mg II) emission-line intensity ratio using our best-fitting model and obtain log($I$(Fe II)/$I$(Mg II)) $\sim$ 0.7, suggesting many AGNs have a roughly solar Fe/Mg abundance ratio. Finally, we vary the Eddington ratio and SED shape as a step in understanding the Eigenvector 1 correlation.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

H, He-like recombination spectra I: $l$-changing collisions for hydrogen

Hydrogen and helium emission lines in nebulae form by radiative recombination. This is a simple process which, in principle, can be described to very high precision. Ratios of He I and H I emission lines can be used to measure the He$^+$/H$^+$ abundance ratio to the same precision as the recombination rate coefficients. This paper investigates the controversy over the correct theory to describe dipole $l$-changing collisions ($nl\rightarrow nl'=l\pm 1$) between energy-degenerate states within an $n$-shell. The work of Pengelly & Seaton (1964) has, for half-a-century, been considered the definitive study which "solved" the problem. Recent work by Vrinceanu et al.(2012) recommended the use of rate coefficients from a semi-classical approximation which are nearly an order of magnitude smaller than those of Pengelly & Seaton (1964), with the result that significantly higher densities are needed for the $nl$ populations to come into local thermodynamic equilibrium. Here, we compare predicted H~I emissivities from the two works and find widespread differences, of up to $\approx 10$%. This far exceeds the 1% precision required to obtain the primordial He/H abundance ratio from observations so as to constrain Big Bang cosmologies. We recommend using the rate coefficients of Pengelly & Seaton (1964) for $l$-changing collisions, to describe the H recombination spectrum, based-on their quantum mechanical representation of the long-range dipole interaction.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

H, He-like recombination spectra II: $l$-changing collisions for He Rydberg states

Cosmological models can be constrained by determining primordial abundances. Accurate predictions of the He I spectrum are needed to determine the primordial helium abundance to a precision of $< 1$% in order to constrain Big Bang Nucleosynthesis models. Theoretical line emissivities at least this accurate are needed if this precision is to be achieved. In the first paper of this series, which focused on H I, we showed that differences in $l$-changing collisional rate coefficients predicted by three different theories can translate into 10% changes in predictions for H I spectra. Here we consider the more complicated case of He atoms, where low-$l$ subshells are not energy degenerate. A criterion for deciding when the energy separation between $l$ subshells is small enough to apply energy-degenerate collisional theories is given. Moreover, for certain conditions, the Bethe approximation originally proposed by Pengelly & Seaton (1964) is not sufficiently accurate. We introduce a simple modification of this theory which leads to rate coefficients which agree well with those obtained from pure quantal calculations using the approach of Vrinceanu et al. (2012). We show that the $l$-changing rate coefficients from the different theoretical approaches lead to differences of $\sim 10$% in He I emissivities in simulations of H II regions using spectral code Cloudy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Accurate determination of the free-free Gaunt factor. II - relativistic Gaunt factors

When modelling an ionised plasma, all spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum. Until now no data set exists that would meet these needs completely. We have therefore produced a table of relativistic Gaunt factors over a much wider range of parameter space than has ever been produced before. We present tables of the thermally averaged Gaunt factor covering the range log10(gamma^2) = -6 to 10 and log10(u) = -16 to 13 for all atomic numbers Z = 1 through 36. The data were calculated using the relativistic Bethe-Heitler-Elwert (BHE) approximation and were subsequently merged with accurate non-relativistic results in those parts of the parameter space where the BHE approximation is not valid. These data will be incorporated in the next major release of the spectral synthesis code Cloudy. We also produced tables of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to 10 for all values of Z between 1 and 36. All the data presented in this paper are available online.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Collisional excitation of [C II], [O I] and CO in Massive Galaxies

Many massive galaxies at the centres of relaxed galaxy clusters and groups have vast reservoirs of cool (~10,000 K) and cold (~100 K) gas. In many low redshift brightest group and cluster galaxies this gas is lifted into the hot ISM in filamentary structures, which are long lived and are typically not forming stars. Two important questions are how far do these reservoirs cool and if cold gas is abundant what is the cause of the low star formation efficiency? Heating and excitation of the filaments from collisions and mixing of hot particles in the surrounding X-ray gas describes well the optical and near infra-red line ratios observed in the filaments. In this paper we examine the theoretical properties of dense, cold clouds emitting in the far infra-red and submillimeter through the bright lines of [C II]157 μm , [O I]63 μm and CO, exposed to these energetic ionising particles. While some emission lines may be optically thick we find this is not sufficient to model the emission line ratios. Models where the filaments are supported by thermal pressure support alone also cannot account for the cold gas line ratios but a very modest additional pressure support, either from turbulence or magnetic fields can fit the observed [O I]/[C II] line ratios by decreasing the density of the gas. This may also help stabilise the filaments against collapse leading to the low rates of star formation. Finally we make predictions for the line ratios expected from cold gas under these conditions and present diagnostic diagrams for comparison with further observations. We provide our code as an Appendix.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Stout: Cloudy's Atomic and Molecular Database

We describe a new atomic and molecular database we developed for use in the spectral synthesis code Cloudy. The design of Stout is driven by the data needs of Cloudy, which simulates molecular, atomic, and ionized gas with kinetic temperatures 2.8 K < T < 1e10 K and densities spanning the low to high-density limits. The radiation field between photon energies $10^{-8}$ Ry and 100 MeV is considered, along with all atoms and ions of the lightest 30 elements, and ~100 molecules. For ease of maintenance, the data are stored in a format as close as possible to the original data sources. Few data sources include the full range of data we need. We describe how we fill in the gaps in the data or extrapolate rates beyond their tabulated range. We tabulate data sources both for the atomic spectroscopic parameters and for collision data for the next release of Cloudy. This is not intended as a review of the current status of atomic data, but rather a description of the features of the database which we will build upon.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A Far-Infrared Spectral Sequence of Galaxies: Trends and Models

We present a framework for the interpretation of the far-infrared spectra of galaxies in which we have expanded the model parameters compared with previous work by varying the ionization parameter $U$, column density $N$(H), and gas density at the cloud face $n$(H$^{+}$) for a central starburst or AGN. We compare these models carried out with the $Cloudy$ spectral synthesis code to trends in line-to-total far-infrared luminosity ratios, far-infrared fine-structure line ratios, IRAS colors, and OH and H$_{2}$O column densities with trends found in the well-studied sample of ten nearby galaxies from the IRAS Bright Galaxy Sample with infrared luminosities greater than 10$^{10}$ L$_{\odot}$ and IRAS 60 micron fluxes equal to or greater than that of the nearby ULIRG Arp 220. We find that the spectral sequence extending from normal starburst-type emission line spectra to ULIRG-type absorption-dominated spectra with significant absorption from excited levels, can be best explained by simultaneously increasing the hydrogen column density, from as low as 10$^{21}$ cm$^{-2}$ to as high as 10$^{24.8}$ cm$^{-2}$ or greater, and the ionization parameter, from as low as 10$^{-4}$ to as high as 1. The starburst models best reproduce most of the sequence, while AGN models are somewhat better able to produce the high OH and H$_{2}$O column densities in Arp 220. Our results suggest that the molecular ISM in ULIRG-like, absorption-dominated systems is located close to and at least partially obscures the source of power throughout much of the far-infrared, which must be taken into account in order to properly interpret diagnostics of both their sources of power and of feedback.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A stubbornly large mass of cold dust in the ejecta of Supernova 1987A

We present new Herschel photometric and spectroscopic observations of Supernova 1987A, carried out in 2012. Our dedicated photometric measurements provide new 70 micron data and improved imaging quality at 100 and 160 micron compared to previous observations in 2010. Our Herschel spectra show only weak CO line emission, and provide an upper limit for the 63 micron [O I] line flux, eliminating the possibility that line contaminations distort the previously estimated dust mass. The far-infrared spectral energy distribution (SED) is well fitted by thermal emission from cold dust. The newly measured 70 micron flux constrains the dust temperature, limiting it to nearly a single temperature. The far-infrared emission can be fitted by 0.5+-0.1 Msun of amorphous carbon, about a factor of two larger than the current nucleosynthetic mass prediction for carbon. The observation of SiO molecules at early and late phases suggests that silicates may also have formed and we could fit the SED with a combination of 0.3 Msun of amorphous carbon and 0.5 Msun of silicates, totalling 0.8 Msun of dust. Our analysis thus supports the presence of a large dust reservoir in the ejecta of SN 1987A. The inferred dust mass suggests that supernovae can be an important source of dust in the interstellar medium, from local to high-redshift galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Accurate determination of the free-free Gaunt factor; I - non-relativistic Gaunt factors

Modern spectral synthesis codes need the thermally averaged free-free Gaunt factor defined over a very wide range of parameter space in order to produce an accurate prediction for the spectrum emitted by an ionized plasma. Until now no set of data exists that would meet this need in a fully satisfactory way. We have therefore undertaken to produce a table of very accurate non-relativistic Gaunt factors over a much wider range of parameters than has ever been produced before. We first produced a table of non-averaged Gaunt factors, covering the parameter space log10(epsilon_i) = -20 to +10 and log10(w) = -30 to +25. We then continued to produce a table of thermally averaged Gaunt factors covering the parameter space log10(gamma^2) = -6 to +10 and log10(u) = -16 to +13. Finally we produced a table of the frequency integrated Gaunt factor covering the parameter space log10(gamma^2) = -6 to +10. All the data presented in this paper are available online.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Herschel imaging of the dust in the Helix Nebula (NGC 7293)

In our series of papers presenting the Herschel imaging of evolved planetary nebulae, we present images of the dust distribution in the Helix nebula (NGC 7293). Images at 70, 160, 250, 350, and 500 micron were obtained with the PACS and SPIRE instruments on board the Herschel satellite. The broadband maps show the dust distribution over the main Helix nebula to be clumpy and predominantly present in the barrel wall. We determined the spectral energy distribution of the main nebula in a consistent way using Herschel, IRAS, and Planck flux values. The emissivity index of 0.99 +/- 0.09, in combination with the carbon rich molecular chemistry of the nebula, indicates that the dust consists mainly of amorphous carbon. The dust excess emission from the central star disk is detected at 70 micron and the flux measurement agree with previous measurement. We present the temperature and dust column density maps. The total dust mass across the Helix nebula (without its halo) is determined to be 0.0035 solar mass at a distance of 216 pc. The temperature map shows dust temperatures between 22 and 42 K, which is similar to the kinetic temperature of the molecular gas, strengthening the fact that the dust and gas co-exist in high density clumps. Archived images are used to compare the location of the dust emission in the far infrared (Herschel) with the ionized (GALEX, Hbeta) and molecular hydrogen component. The different emission components are consistent with the Helix consisting of a thick walled barrel-like structure inclined to the line of sight. The radiation field decreases rapidly through the barrel wall.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Herschel Planetary Nebula Survey (HerPlaNS) - First Detection of OH+ in Planetary Nebulae

We report the first detections of OH$^+$ emission in planetary nebulae (PNe). As part of an imaging and spectroscopy survey of 11 PNe in the far-IR using the PACS and SPIRE instruments aboard the Herschel Space Observatory, we performed a line survey in these PNe over the entire spectral range between 51 and 672$μ$m to look for new detections. OH$^+$ rotational emission lines at 152.99, 290.20, 308.48, and 329.77$μ$m were detected in the spectra of three planetary nebulae: NGC 6445, NGC 6720, and NGC 6781. Excitation temperatures and column densities derived from these lines are in the range of 27 to 47 K and 2$\times$10$^{10}$ to 4 $\times$10$^{11}$ cm$^{-2}$, respectively. In PNe, the OH+ rotational line emission appears to be produced in the photodissociation region (PDR) in these objects. The emission of OH+ is observed only in PNe with hot central stars (T$_{eff}$ > 100000 K), suggesting that high-energy photons may play a role in the OH+ formation and its line excitation in these objects, as it seems to be the case for ultraluminous galaxies.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Herschel spectral-mapping of the Helix Nebula (NGC 7293): Extended CO photodissociation and OH+ emission

The Helix Nebula (NGC 7293) is the closest planetary nebulae. Therefore, it is an ideal template for photochemical studies at small spatial scales in planetary nebulae. We aim to study the spatial distribution of the atomic and the molecular gas, and the structure of the photodissociation region along the western rims of the Helix Nebula as seen in the submillimeter range with Herschel. We use 5 SPIRE FTS pointing observations to make atomic and molecular spectral maps. We analyze the molecular gas by modeling the CO rotational lines using a non-local thermodynamic equilibrium (non-LTE) radiative transfer model. For the first time, we have detected extended OH+ emission in a planetary nebula. The spectra towards the Helix Nebula also show CO emission lines (from J= 4 to 8), [NII] at 1461 GHz from ionized gas, and [CI] (2-1), which together with the OH+ lines, trace extended CO photodissociation regions along the rims. The estimated OH+ column density is (1-10)x1e12 cm-2. The CH+ (1-0) line was not detected at the sensitivity of our observations. Non-LTE models of the CO excitation were used to constrain the average gas density (n(H2)=(1-5)x1e5 cm-3) and the gas temperature (Tk= 20-40 K). The SPIRE spectral-maps suggest that CO arises from dense and shielded clumps in the western rims of the Helix Nebula whereas OH+ and [CI] lines trace the diffuse gas and the UV and X-ray illuminated clumps surface where molecules reform after CO photodissociation. [NII] traces a more diffuse ionized gas component in the interclump medium.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Implications of Coronal Line Emission in NGC 4696

We announce a new facility in the spectral code CLOUDY that enables tracking the evolution of a cooling parcel of gas with time. For gas cooling from temperatures relevant to galaxy clusters, earlier calculations estimated the [Fe XIV] λ5303 / [Fe X] λ6375 luminosity ratio, a critical diagnostic of a cooling plasma, to slightly less than unity. By contrast, our calculations predict a ratio ~3. We revisit recent optical coronal line observations along the X-ray cool arc around NGC 4696 by Canning et al. (2011), which detected [Fe X] λ6375, but not [Fe XIV] λ5303. We show that these observations are not consistent with predictions of cooling flow models. Differential extinction could in principle account for the observations, but it requires extinction levels (A_V > 3.625) incompatible with previous observations. The non-detection of [Fe XIV] implies a temperature ceiling of 2.1 million K. Assuming cylindrical geometry and transonic turbulent pressure support, we estimate the gas mass at ~1 million solar masses. The coronal gas is cooling isochorically. We propose that the coronal gas has not condensed out of the intracluster medium, but instead is the conductive or mixing interface between the X-ray plume and the optical filaments. We present a number of emission lines that may be pursued to test this hypothesis and constrain the amount of intermediate temperature gas in the system.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radiative Cooling II: Effects of Density and Metallicity

This work follows Lykins et al. discussion of classic plasma cooling function at low density and solar metallicity. Here we focus on how the cooling function changes over a wide range of density (n_H<10^12 cm^(-3)) and metallicity (Z<30Z _sun ). We find that high densities enhance the ionization of elements such as hydrogen and helium until they reach local thermodynamic equilibrium. By charge transfer, the metallicity changes the ionization of hydrogen when it is partially ionized. We describe the total cooling function as a sum of four parts: those due to H&He, the heavy elements, electron-electron bremsstrahlung and grains. For the first 3 parts, we provide a low-density limit cooling function, a density dependence function, and a metallicity dependence function. These functions are given with numerical tables and analytical fit functions. For grain cooling, we only discuss in ISM case. We then obtain a total cooling function that depends on density, metallicity and temperature. As expected, collisional de-excitation suppresses the heavy elements cooling. Finally, we provide a function giving the electron fraction, which can be used to convert the cooling function into a cooling rate.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The evolving spectrum of the planetary nebula Hen 2-260

We analysed the planetary nebula Hen 2-260 using optical spectroscopy and photometry. We compared our observations with the data from literature to search for evolutionary changes. The nebular line fluxes were modelled with the Cloudy photoionization code to derive the stellar and nebular parameters. The planetary nebula shows a complex structure and possibly a bipolar outflow. The nebula is relatively dense and young. The central star is just starting $\rm O^+$ ionization ($\rm T_{eff} \approx 30,000 \, K$). Comparison of our observations with literature data indicates a 50% increase of the [OIII] 5007 Å line flux between 2001 and 2012. We interpret it as the result of the progression of the ionization of $\rm O^{+}$. The central star evolves to higher temperatures at a rate of $\rm 45 \pm 7\,K\, yr^{-1}$. The heating rate is consistent with a final mass of $\rm 0.626 ^{+0.003}_{-0.005} \, M_{\odot}$ or $\rm 0.645 ^{+0.008}_{-0.008} \, M_{\odot}$ for two different sets of post-AGB evolutionary tracks from literature. The photometric monitoring of Hen 2-260 revealed variations on a timescale of hours or days. The variability may be caused by pulsations of the star. The temperature evolution of the central star can be traced using spectroscopic observations of the surrounding planetary nebula spanning a timescale of roughly a decade. This allows us to precisely determine the stellar mass, since the pace of the temperature evolution depends critically on the core mass. The kinematical age of the nebula is consistent with the age obtained from the evolutionary track. The final mass of the central star is close to the mass distribution peak for central stars of planetary nebulae found in other studies. The object belongs to a group of young central stars of planetary nebulae showing photometric variability.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The Very Fast Evolution of the VLTP Object V4334 Sgr

V4334 Sgr (Sakurai's object) is an enigmatic evolved star that underwent a very late thermal pulse a few years before its discovery in 1996. It ejected a new, hydrogen-deficient nebula in the process. Emission lines from the newly ejected gas were first discovered in 1998 (He I 1083 nm) and 2001 (optical). We have monitored the optical emission spectrum since. From 2001 through 2007 the optical spectrum showed an exponential decline in flux, consistent with a shock that occurred around 1998 and started cooling soon after that. In this paper we show that since 2008 the line fluxes have been continuously rising again. Our preliminary interpretation is that this emission comes from a region close to the central star, and is excited by a second shock. This shock may have been induced by an increase in the stellar mass loss and wind velocity associated with a rise in the stellar temperature.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Herschel study of NGC 650

As part of the Herschel Guaranteed Time Key Project MESS (Mass loss of Evolved StarS) we have imaged a sample of planetary nebulae. In this paper we present the PACS and SPIRE images of the classical bipolar planetary nebula NGC 650. We used these images to derive a temperature map of the dust. We also constructed a photoionization and dust radiative transfer model using the spectral synthesis code Cloudy. To constrain this model, we used the PACS and SPIRE fluxes and combined these with hitherto unpublished IUE and Spitzer IRS spectra as well as various other data from the literature. The temperature map combined with the photoionization model were used to study various aspects of the central star, the nebula, and in particular the dust grains in the nebula. The central star parameters are determined to be T_eff = 208 kK and L = 261 L_sol assuming a distance of 1200 pc. The stellar temperature is much higher than previously published values. We confirm that the nebula is carbon-rich with a C/O ratio of 2.1. The nebular abundances are typical for a type IIa planetary nebula. With the photoionization model we determined that the grains in the ionized nebula are large (assuming single-sized grains, they would have a radius of 0.15 micron. Most likely these large grains were inherited from the asymptotic giant branch phase. The PACS 70/160 micron temperature map shows evidence for two radiation components heating the grains. The first component is direct emission from the central star, while the second component is diffuse emission from the ionized gas (mainly Ly alpha). We show that previous suggestions that there is a photo-dissociation region surrounding the ionized region are incorrect. The neutral material resides in dense clumps inside the ionized region. These may also harbor stochastically heated very small grains in addition to the large grains.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Carbon Monoxide in the Cold Debris of Supernova 1987A

We report spectroscopic and imaging observations of rotational transitions of cold CO and SiO in the ejecta of SN1987A, the first such emission detected in a supernova remnant. In addition to line luminosities for the CO J=1-0, 2-1, 6-5, and 7-6 transitions, we present upper limits for all other transitions up to J=13-12, collectively measured from the Atacama Large Millimeter Array (ALMA), the Atacama Pathfinder EXperiment (APEX), and the Herschel Spectral and Photometric Imaging REceiver (SPIRE). Simple models show the lines are emitted from at least 0.01 solar masses of CO at a temperature > 14 K, confined within at most 35% of a spherical volume expanding at ~ 2000 km/s. Moreover, we locate the emission within 1'' of the central debris. These observations, along with a partial observation of SiO, confirm the presence of cold molecular gas within supernova remnants and provide insight into the physical conditions and chemical processes in the ejecta. Furthermore, we demonstrate the powerful new window into supernova ejecta offered by submillimeter observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

CKVul: evolving nebula and three curious background stars

We analyse the remnants of CK Vul (Nova Vul 1670) using optical imaging and spectroscopy. The imaging, obtained between 1991 and 2010, spans 5.6% of the life-time of the nebula. The flux of the nebula decreased during the last 2 decades. The central source still maintains the ionization of the innermost part of the nebula, but recombination proceeds in more distant parts of the nebula. Surprisingly, we discovered two stars located within 10 arcsec of the expansion centre of the radio emission that are characterized by pronounced long term variations and one star with high proper motion. The high proper motion star is a foreground object, and the two variable stars are background objects. The photometric variations of two variables are induced by a dusty cloud ejected by CK Vul and passing through the line of sight to those stars. The cloud leaves strong lithium absorption in the spectra of the stars. We discuss the nature of the object in terms of recent observations.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Expanded Iron UTA spectra -- probing the thermal stability limits in AGN clouds

The Fe unresolved transition array (UTAs) produce prominent features in the 15-17?A wavelength range in the spectra of Active Galactic Nuclei (AGN). Here we present new calculations of the energies and oscillator strengths of inner- shell lines from Fe XIV, Fe XV, and Fe XVI. These are crucial ions since they are dominant at inflection points in the gas thermal stability curve, and UTA excitation followed by autoionization is an important ionization mechanism for these species. We incorporate these, and data reported in previous papers, into the plasma simulation code Cloudy. This updated physics is subsequently employed to reconsider the thermally stable phases in absorbing media in Active Galactic Nuclei. We show how the absorption profile of the Fe XIV UTA depends on density, due to the changing populations of levels within the ground configuration.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Herschel SPIRE and PACS observations of the red supergiant VY CMa: analysis of the molecular line spectra

We present an analysis of the far-infrared and submillimetre molecular emission line spectrum of the luminous M-supergiant VY CMa, observed with the SPIRE and PACS spectrometers aboard the Herschel Space Observatory. Over 260 emission lines were detected in the 190-650-micron SPIRE FTS spectra, with one-third of the observed lines being attributable to H2O. Other detected species include CO, 13CO, H2^18O, SiO, HCN, SO, SO2, CS, H2S, and NH3. Our model fits to the observed 12CO and 13CO line intensities yield a 12C/13C ratio of 5.6+-1.8, consistent with measurements of this ratio for other M supergiants, but significantly lower than previously estimated for VY CMa from observations of lower-J lines. The spectral line energy distribution for twenty SiO rotational lines shows two temperature components: a hot component at 1000 K, which we attribute to the stellar atmosphere and inner wind, plus a cooler ~200 K component, which we attribute to an origin in the outer circumstellar envelope. We fit the line fluxes of 12CO, 13CO, H2O and SiO, using the SMMOL non-LTE line transfer code, with a mass-loss rate of 1.85x10^-4 Msun yr^-1 between 9 R* and 350 R*. To fit the observed line fluxes of 12CO, 13CO, H2O and SiO with SMMOL non-LTE line radiative transfer code, along with a mass-loss rate of 1.85x10^-4 Msun yr^-1. To fit the high rotational lines of CO and H2O, the model required a rather flat temperature distribution inside the dust condensation radius, attributed to the high H2O opacity. Beyond the dust condensation radius the gas temperature is fitted best by an r^-0.5 radial dependence, consistent with the coolant lines becoming optically thin. Our H2O emission line fits are consistent with an ortho:para ratio of 3 in the outflow.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Near-infrared photometry of Galactic planetary nebulae with the VVV Survey

Planetary nebulae (PNe) are powerful tracers of evolved stellar populations. Among the 3000 known PNe in the Galaxy, about 600 are located within the 520 square-degree area covered by the VVV survey. The VVV photometric catalogue provides an important new dataset for the study of PNe, with high-resolution imaging in five near-infrared bands. Aperture photometry of known PNe in the VVV area was retrieved from source catalogues. Care was taken to minimise any confusion with field stars. The colours of the PNe we are determined for H-Ks, J-H, Z-Y, and Y-J, and compared to stars and to other types of emission line objects. Cloudy photo-ionisation models were used to predict colours for typical PNe. We present near-infrared photometry for 353 known PNe. The best separation from other objects is obtained in the H-Ks vs. J-H diagram. We calculated the emission-line contribution to the in-band flux based on a model for NGC 6720: we find that this is highest in the Z and Y bands at over 50%, lower in the J band at 40%, and lowest in the H and Ks bands at 20%. A new view of PNe in the wavelength domain of the Z and Y bands is shown. Photo-ionisation models are used to explore the observed colours in these bands. The Y band is shown to be dominated by HeI 1.083 mu and HeII 1.012 mu, and colours involving this band are very sensitive to the temperature of the ionizing star. The VVV survey represents a unique dataset for studing crowded and obscured regions in the Galactic plane. The diagnostic diagrams presented here allow one to study the properties of known PNe and to uncover objects not previously classified.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The 2013 Release of Cloudy

This is a summary of the 2013 release of the plasma simulation code Cloudy. Cloudy models the ionization, chemical, and thermal state of material that may be exposed to an external radiation field or other source of heating, and predicts observables such as emission and absorption spectra. It works in terms of elementary processes, so is not limited to any particular temperature or density regime. This paper summarizes advances made since the last major review in 1998. Much of the recent development has emphasized dusty molecular environments, improvements to the ionization / chemistry solvers, and how atomic and molecular data are used. We present two types of simulations to demonstrate the capability of the code. We consider a molecular cloud irradiated by an X-ray source such as an Active Nucleus and show how treating EUV recombination lines and the full SED affects the observed spectrum. A second example illustrates the very wide range of particle and radiation density that can be considered.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Hydrogen Two-Photon Continuum Emission from the Horseshoe Filament in NGC 1275

Far ultraviolet emission has been detected from a knot of Halpha emission in the Horseshoe filament, far out in the NGC 1275 nebula. The flux detected relative to the brightness of the Halpha line in the same spatial region is very close to that expected from Hydrogen two-photon continuum emission in the particle heating model of Ferland et al. (2009) if reddening internal to the filaments is taken into account. We find no need to invoke other sources of far ultraviolet emission such as hot stars or emission lines from CIV in intermediate temperature gas to explain these data.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Pumping up the [N I] nebular lines

The optical [N I] doublet near 5200 Å is anomalously strong in a variety of emission-line objects. We compute a detailed photoionization model and use it to show that pumping by far-ultraviolet (FUV) stellar radiation previously posited as a general explanation applies to the Orion Nebula (M42) and its companion M43; but, it is unlikely to explain planetary nebulae and supernova remnants. Our models establish that the observed nearly constant equivalent width of [N I] with respect to the dust-scattered stellar continuum depends primarily on three factors: the FUV to visual-band flux ratio of the stellar population; the optical properties of the dust; and the line broadening where the pumping occurs. In contrast, the intensity ratio [N I]/Hβ depends primarily on the FUV to extreme-ultraviolet ratio, which varies strongly with the spectral type of the exciting star. This is consistent with the observed difference of a factor of five between M42 and M43, which are excited by an O7 and B0.5 star respectively. We derive a non-thermal broadening of order 5 km/s for the [N I] pumping zone and show that the broadening mechanism must be different from the large-scale turbulent motions that have been suggested to explain the line-widths in this H II region. A mechanism is required that operates at scales of a few astronomical units, which may be driven by thermal instabilities of neutral gas in the range 1000 to 3000 K. In an appendix, we describe how collisional and radiative processes are treated in the detailed model N I atom now included in the Cloudy plasma code.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Physical Conditions in Barnard's Loop, Components of the Orion-Eridanus Bubble, and Implications for the WIM Component of the ISM

We have supplemented existing spectra of Barnard's Loop with high accuracy spectrophotometry of one new position. Cloudy photoionization models were calculated for a variety of ionization parameters and stellar temperatures and compared with the observations. After testing the procedure with recent observations of M43, we establish that Barnard's Loop is photoionized by four candidate ionizing stars, but agreement between the models and observations is only possible if Barnard's Loop is enhanced in heavy elements by about a factor of 1.4. Barnard's Loop is very similar in properties to the brightest components of the Orion-Eridanus Bubble and the Warm Ionized Medium (WIM). We are able to establish models that bound the range populated in low-ionization color-color diagrams (I([SII])/I(Hα) versus I([NII])/I(Hα)) using only a limited range of ionization parameters and stellar temperatures. Previously established variations in the relative abundance of heavy elements render uncertain the most common method of determining electron temperatures for components of the Orion-Eridanus Bubble and the WIM based on only the I([NII])/I(Hα) ratio, although we confirm that the lowest surface brightness components of the WIM are on average of higher electron temperature. The electron temperatures for a few high surface brightness WIM components determined by direct methods are comparable to those of classical bright H II regions. In contrast, the low surface brightness HII regions studied by the Wisconsin Hα Mapper are of lower temperatures than the classical bright HII regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Herschel images of NGC 6720: H2 formation on dust grains

Herschel PACS and SPIRE images have been obtained of NGC 6720 (the Ring Nebula). This is an evolved planetary nebula with a central star that is currently on the cooling track, due to which the outer parts of the nebula are recombining. From the PACS and SPIRE images we conclude that there is a striking resemblance between the dust distribution and the H2 emission, which appears to be observational evidence that H2 forms on grain surfaces. We have developed a photoionization model of the nebula with the Cloudy code which we used to determine the physical conditions of the dust and investigate possible formation scenarios for the H2. We conclude that the most plausible scenario is that the H2 resides in high density knots which were formed after the recombination of the gas started when the central star entered the cooling track. Hydrodynamical instabilities due to the unusually low temperature of the recombining gas are proposed as a mechanism for forming the knots. H2 formation in the knots is expected to be substantial after the central star underwent a strong drop in luminosity about one to two thousand years ago, and may still be ongoing at this moment, depending on the density of the knots and the properties of the grains in the knots.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Dust-Bounded ULIRGs? Model Predictions for Infrared Spectroscopic Surveys

The observed faintness of infrared fine-structure line emission along with the warm far-infrared (FIR) colors of ultraluminous infrared galaxies (ULIRGs) is a long-standing problem. In this work, we calculate the line and continuum properties of a cloud exposed to an Active Galactic Nucleus (AGN) and starburst spectral energy distribution (SED). We use an integrated modeling approach, predicting the spectrum of ionized, atomic, and molecular environments in pressure equilibrium. We find that the effects of high ratios of impinging ionizing radiation density to particle density (i.e. high ionization parameters, or U) can reproduce many ULIRG observational characteristics. Physically, as U increases, the fraction of UV photons absorbed by dust increases, corresponding to fewer photons available to photoionize and heat the gas, producing what is known as a "dust-bounded" nebula. We show that high U effects can explain the "[C II] deficit", the ~1 dex drop in the [C II] 158 micron /FIR ratio seen in ULIRGs when compared to starburst or normal galaxies. Additionally, by increasing U through increasing the ionizing photon flux, warmer dust and thus higher IRAS F(60)/F(100) ratios result. High U effects also predict an increase in [O I]63 micron /[C II] 158 micron and a gradual decline in [O III] 88 micron /FIR, similar to the magnitude of the trends observed, and yield a reasonable fit to [Ne V]14 micron /FIR ratio AGN observations.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

A PDR-Code Comparison Study

We present a comparison between independent computer codes, modeling the physics and chemistry of interstellar photon dominated regions (PDRs). Our goal was to understand the mutual differences in the PDR codes and their effects on the physical and chemical structure of the model clouds, and to converge the output of different codes to a common solution. A number of benchmark models have been created, covering low and high gas densities and far ultraviolet intensities. The benchmark models were computed in two ways: one set assuming constant temperatures, thus testing the consistency of the chemical network and photo-processes, and a second set determining the temperature selfconsistently. We investigated the impact of PDR geometry and agreed on the comparison of results from spherical and plane-parallel PDR models. We identified a number of key processes governing the chemical network which have been treated differently in the various codes, and defined a proper common treatment. We established a comprehensive set of reference models for ongoing and future PDR model bench-marking and were able to increase the agreement in model predictions for all benchmark models significantly.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The central stars of the planetary nebulae NGC 7027 and NGC 6543

Infrared spectra of NGC 7027 and NGC 6543 ranging from 2.4 to 45 micron were obtained with the Short Wavelength Spectrometer on board the Infrared Space Observatory. A first analysis of these spectra, with the aid of photo-ionization models, is presented. We report the first detection of the Ar VI 4.53 micron and Ne VI 7.65 micron lines in the spectrum of NGC 7027. When compared with older observations it is clear that the Ar VI line and possibly also other lines have increased in strength since 1981. We argue that a likely explanation for this variability is a change in the spectral energy distribution of the central star, possibly an increase in effective temperature. However, this result needs to be confirmed by further observations. We also report a non-detection of the O IV 25.9 micron line and the first detection of the Na III 7.32 micron line in the spectrum of NGC 6543. The non-detection is not expected based on a blackbody approximation for the spectrum of the central star. The ionization threshold for O(3+) is just beyond the He II limit, and the absence of this line shows that the stellar flux drops at least by a factor 350 at the He II limit. Modeling the O IV line may prove to be a valuable test for atmosphere models.

0 citations0 reviewsNo code linkedOpen access
preprint1999arXiv

The spectral evolution of post-AGB stars

A parameter study of the spectral evolution of a typical post-AGB star, with particular emphasis on the evolution of the infrared colours, is presented. The models are based on the latest evolutionary tracks for hydrogen burning post-AGB stars. For such tracks the evolutionary rate is very dependent on the assumed mass loss rate as a function of time. We investigate this effect by modifying the mass loss prescription. The newly calculated evolutionary rates and density distributions are used to model the spectral evolution of a post-AGB star with the photo-ionization code Cloudy, including dust in the radiative transfer. Different assumptions for the dust properties and dust formation are considered. It is shown that by varying these parameters in a reasonable way, entirely different paths are followed in the IRAS colour-colour diagram. First of all, the effects of the evolution of the central star on the expanding dust shell can not be neglected. Also the dust properties and the definition of the end of the AGB phase have an important effect. The model tracks show that objects occupying the same location in the IRAS colour-colour diagram can have a different evolutionary past, and therefore the position in the IRAS colour-colour diagram alone can not a priori give a unique determination of the evolutionary status of an object. An alternative colour-colour diagram, the K-[12] vs. [12]-[25] diagram, is presented. The tracks in this diagram seem less affected by particulars of the grain emission. This diagram may be a valuable additional tool for studying post-AGB evolution.

0 citations0 reviewsNo code linkedOpen access