Source author record

R. L. Porter

R. L. Porter 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

16works
10topics
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

16 published item(s)

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.

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.

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.

preprint2015arXiv

The Abundances of Light Neutron-Capture Elements in Planetary Nebulae III. The Impact of New Atomic Data on Nebular Selenium and Krypton Abundance Determinations

The detection of neutron(n)-capture elements in several planetary nebulae (PNe) has provided a new means of investigating s-process nucleosynthesis in low-mass stars. However, a lack of atomic data has inhibited accurate trans-iron element abundance determinations in astrophysical nebulae. Recently, photoionization and recombination data were determined for Se and Kr, the two most widely detected n-capture elements in nebular spectra. We have incorporated these new data into the photoionization code Cloudy. To test the atomic data, numerical models were computed for 15 PNe that exhibit emission lines from multiple Kr ions. We found systematic discrepancies between the predicted and observed emission lines that are most likely caused by inaccurate photoionization and recombination data. These discrepancies were removed by adjusting the Kr$^+$--Kr$^{3+}$ photoionization cross sections within their cited uncertainties and the dielectronic recombination rate coefficients by slightly larger amounts. From grids of models spanning the physical conditions encountered in PNe, we derive new, broadly applicable ionization correction factor (ICF) formulae for calculating Se and Kr elemental abundances. The ICFs were applied to our previous survey of near-infrared [Kr III] and [Se IV] emission lines in 120 PNe. The revised Se and Kr abundances are 0.1-0.3 dex lower than former estimates, with average values of [Se/(O, Ar)]=0.12$\pm$0.27 and [Kr/(O, Ar)]=0.82$\pm$0.29, but correlations previously found between their abundances and other nebular and stellar properties are unaffected. We also find a tendency for high-velocity PNe that can be associated with the Galactic thick disk to exhibit larger s-process enrichments than low-velocity PNe belonging to the thin disk population.

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.

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.

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.

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.

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.

preprint2013arXiv

The primordial helium abundance from updated emissivities

Observations of metal-poor extragalactic H II regions allow the determination of the primordial helium abundance, Y_p. The He I emissivities are the foundation of the model of the H II region's emission. Porter, Ferland, Storey, & Detisch (2012) have recently published updated He I emissivities based on improved photoionization cross-sections. We incorporate these new atomic data and update our recent Markov Chain Monte Carlo analysis of the dataset published by Izotov, Thuan, & Stasinska (2007). As before, cuts are made to promote quality and reliability, and only solutions which fit the data within 95% confidence level are used to determine the primordial He abundance. The previously qualifying dataset is almost entirely retained and with strong concordance between the physical parameters. Overall, an upward bias from the new emissivities leads to a decrease in Y_p. In addition, we find a general trend to larger uncertainties in individual objects (due to changes in the emissivities) and an increased variance (due to additional objects included). From a regression to zero metallicity, we determine Y_p = 0.2465 +/- 0.0097, in good agreement with the Planck result of Y_p = 0.2485 +/- 0.0002. In the future, a better understanding of why a large fraction of spectra are not well fit by the model will be crucial to achieving an increase in the precision of the primordial helium abundance determination.

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.

preprint2012arXiv

Improved He I Emissivities in the Case B Approximation

We update our prior work on the case B collisional-recombination spectrum of He I to incorporate \textit{ab initio} photoionisation cross-sections. This large set of accurate, self-consistent cross-sections represents a significant improvement in He I emissivity calculations because it largely obviates the piecemeal nature that has marked all modern works. A second, more recent set of \textit{ab initio} cross-sections is also available, but we show that those are less consistent with bound-bound transition probabilities than our adopted set. We compare our new effective recombination coefficients with our prior work and our new emissivities with those by other researchers, and we conclude with brief remarks on the effects of the present work on the He I error budget. Our calculations cover temperatures $5000 \le T_e \le 25000$ K and densities $10^1 \le n_e \le 10^{14}$ cm$^{-3}$. Full results are available online.

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.

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.

preprint2008arXiv

Fluorescent excitation of Balmer lines in gaseous nebulae: case D

(abridged) Non-ionizing stellar continua are a source of photons for continuum pumping in the hydrogen Lyman transitions. In the environments where these transitions are optically thick, deexcitation occurs through higher series lines, so that the flux in these lines has a fluorescent contribution in addition to recombination; in particular, Balmer emissivities are systematically enhanced above case B. The effectiveness of such mechanism in HII regions and the adequacy of photoionization models as a tool to study it are the two main focuses of this work. We find that photoionization models of H II regions illuminated by low-resolution population synthesis models significantly overpredict the fluorescent contribution to the Balmer lines. Conversely, photoionization models in which the non-ionizing part of the continuum is omitted or is not transferred underpredict the fluorescent contribution to the Balmer lines, producing a bias of similar amplitude in the opposite direction. In this paper, we carry out realistic estimations of the fluorescent Balmer intensity and discuss the variations to be expected as the simulated observational setup and the stellar population's parameters are varied. In all the cases explored, we find that fluorescent excitation provides a significant contribution. We also show that differential fluorescent enhancement may produce line-of-sight differences in the Balmer decrement, mimicking interstellar extinction. Fluorescent excitation emerges from our study as a small but important mechanism for the enhancement of Balmer lines, which should be taken into account in the abundance analysis of photoionized regions, particularly in the case of high-precision applications such as the determination of primordial helium.