P. J. Storey

Bio: P. J. Storey is an academic researcher from University College London. The author has contributed to research in topics: Planetary nebula & Doubly ionized oxygen. The author has an hindex of 43, co-authored 139 publications receiving 7775 citations. Previous affiliations of P. J. Storey include ASTRON.

Recombination line intensities for hydrogenic ions-IV. Total recombination coefficients and machine-readable tables for Z=1 to 8

Abstract: Line emissivities, effective recombination coefficients, opacity factors, departure coefficients and total recombination coefficients are calculated for hydrogenic ions with Z less-than-or-equal-to 8. Results are obtained for Cases A and B for n less-than-or-equal-to 50. Collisional transitions among individual n and l states are fully treated. Calculations were made for log N(e) = 2(1)14 for Case B and log N(e) = 2(1)10 for Case A. The electron temperature takes between nine and 12 values, lying within the range 500 to 100 000 K, depending on the ion. All results are available in the form of machine-readable files.Secondary files containing only effective emissivities for transitions for n less-than-or-equal-to 25 and total recombination coefficients are also available for use with an interactive data server. The server produces tables of relative intensities of any two specified transitions or emissivities for any transition at all temperatures and densities in the data set. Extensive facilities for two-dimensional interpolation of relative intensities, emissivities and total recombination coefficients are provided.

Recombination-line intensities for hydrogenic ions – I. Case B calculations for H I and He II

Abstract: Assuming Case B of Baker and Menzel (1938), the relative intensities of H I and He II recombination lines are determined for a larger range of temperature and density than previously considered and with the upper principle quantum number up to 50 and the lower number up to 29. Full collisional effects are included, and new collision strengths for the n = 1, 2, and 3 states of He(+) are presented. Case B theory is found to be invalidated for collisional excitation of the n = 3 levels from both n = 1 and n = 2 states under certain conditions which depend upon electron density and the Lyman-alpha escape probability, and the regimes of temperature and density for which Case B is valid are also considered.

Theoretical values for the [O iii] 5007/4959 line-intensity ratio and homologous cases

Abstract: The magnetic-dipole probabilities for the D-1(2)-P-3(2) and D-1(2)-P-3(1) transitions in carbon-like and oxygen-like ions are calculated up to atomic number Z = 12, including relativistic corrections to the magnetic dipole operator. The ratio of the probabilities for these two transitions is found to change by up to 5 per cent compared with previous theoretical work, none of which included these relativistic corrections, with the effect being largest for the near neutral ions. The transition probability ratio for the [O III] 5007 and 4959 Angstrom lines is found to be 3.01, implying an intensity ratio of 2.98, in significantly better agreement with the observed ratio than the earlier theoretical work.

Atomic data for opacity calculations. II. Computational methods

Abstract: For pt.I see ibid., vol.20, p.6363-78 (1987). A general description of the data requirements for opacity calculations has been given in paper I. The present paper gives a detailed description of the methods being used in a collaborative effort which is referred to as the Opacity Project. The close-coupling approximation of electron-atom collision theory is used to calculate energies and wavefunctions for bound states, oscillator strengths, photoionisation cross sections and parameters for line broadening by electron impact. The computations are made using the R-matrix method together with new codes for calculating outer-region solutions and dipole integrals. Use of these techniques provides an efficient means of calculating large amounts of accurate atomic data.

NGC 6153: a super-metal-rich planetary nebula?

Abstract: We have obtained deep optical spectra of the planetary nebula NGC 6153, both along its minor axis and by uniformly scanning a long slit across the whole nebula. The scanned spectra, when combined with the nebular total H beta flux, yield integrated fluxes for all the lines (similar to 400) in our spectra, which are rich in strong recombination lines from C, N, O and Ne ions. A weak O vi lambda 3811 emission line from the central star has been detected, suggesting that the nucleus of NGC 6153 has a hydrogen-deficient surface. The optical data, together with the ISO LWS 43-197 mu m spectrum and the archival IUE and IRAS LRS spectra, are used to study the thermal and density structure and to derive the heavy-element abundances from lines produced by different excitation mechanisms. In all cases, the C2+/H+, N2+/H+, O2+/H+ and Ne2+/H+ abundances derived from multiple optical recombination lines (ORLs) are consistently higher, by about a factor of 10, than the corresponding values deduced from optical, UV or infrared (IR) collisionally excited lines (CELs), regardless of the excitation energies or critical densities of the latter. The agreement between the temperature-sensitive optical forbidden lines and the temperature-insensitive IR fine-structure lines rules out temperature fluctuations as the cause of the large difference between the ORL and CEL abundances.We present the results of a new calculation of recombination coefficients for [O II] which lead to good agreement between the observed and predicted [O II] lambda lambda 7320, 7330 forbidden line intensities if these lines are solely excited by recombination at the Balmer jump temperature. Recombination excitation is also found to be important in exciting the [N ii] lambda 5754 line, which, if unaccounted for, would lead to an overestimated [N ii] temperature from the observed (lambda 6548+lambda 6584)/lambda 5754 ratio. Analysis of a number of C ii lines arising from levels as high as 7g in the recombination ladder reveals excellent agreement between their reddening-corrected relative intensities and those predicted by recombination theory. Spatial analysis of the long-slit spectra taken along the nebular minor axis yields a varying [O iii] temperature, whereas the hydrogen Balmer jump temperature of 6000 K is approximately constant across the nebula, and is 2000-3000 K lower than the [O iii] temperature. The observed high-n Balmer line decrement indicates that the hydrogen lines arise from material having an electron density of 2000(-1000)(+2000) cm(-3), consistent with the optical and IR forbidden-line density diagnostics, which yield average line-of-sight electron densities along the minor axis varying between 2000 and 4000 cm(-3).While the He/H ratio mapped by He I and He ii recombination lines is constant within 5 per cent across the nebula, the C2+/H+ and O2+/H+ recombination-line abundances decrease by a factor of 2-3 over a radius of 15 arcsec from the centre, pointing to the presence of abundance gradients. We consider a variety of hypotheses to account for the observed behaviour of the various thermal, density and abundance diagnostics. Empirical nebular models containing two components with differing densities and temperatures are able to account for many of the observed patterns, but only if one of the components is significantly hydrogen-deficient. One such model, which gives a good fit to the observed line intensities and patterns, has 500-K H-depleted material, presumed to be evaporating from dense neutral inclusions, embedded in 9500-K material with 'normal' abundances. An alternative model, which appears more physically plausible on a number of grounds, has high-density (2x10(6) cm(-3)), fully ionized, H-deficient knots embedded in the 'normal' component, although this model fails to account adequately for the observed low (6000 K) hydrogen Balmer jump temperature. However, the observed fact that the ORLs and CELs yield heavy-element abundance ratios that are identical within the uncertainties finds no obvious explanation in the context of H-deficient knot models.

CLOUDY 90: Numerical Simulation of Plasmas and Their Spectra

Abstract: CLOUDY is a large‐scale spectral synthesis code designed to simulate fully physical conditions within an astronomical plasma and then predict the emitted spectrum. Here we describe version 90 (C90) of the code, paying particular attention to changes in the atomic database and numerical methods that have affected predictions since the last publicly available version, C84. The computational methods and uncertainties are outlined together with the direction future development will take. The code is freely available and is widely used in the analysis and interpretation of emission‐line spectra. Web access to the Fortran source for CLOUDY, its documentation Hazy, and an independent electronic form of the atomic database is also described.

A grid of MARCS model atmospheres for late-type stars. I. Methods and general properties

Abstract: Context. In analyses of stellar spectra and colours, and for the analysis of integrated light from galaxies, a homogeneous grid of model atmospheres of late-type stars and corresponding flux spectra is needed. Aims. We construct an extensive grid of spherically-symmetric models (supplemented with plane-parallel ones for the highest surface gravities), built on up-to-date atomic and molecular data, and make it available for public use. Methods. The most recent version of the MARCS program is used. Results. We present a grid of about 104 model atmospheres for stars with 2500K <= T-eff <= 8000 K, -1 <= log g = log (GM/R-2) <= 5 (cgs) with various masses and radii, -5 <= [Me/H] <= + 1, with [alpha/Fe] = 0.0 and 0.4 and different choices of C and N abundances. This includes "CN-cycled" models with C/N=4.07 (solar), 1.5 and 0.5, C/O ranging from 0.09 to (normally) 5.0 to also represent stars of spectral types R, S and N, and with 1.0 <= xi(t) = 5km s(-1). We also list thermodynamic quantities (T, P-g, P-e, rho, partial pressures of molecules, etc.) and provide them on the World Wide Web, as well as calculated fluxes in approximately 108 000 wavelength points. Underlying assumptions in addition to 1D stratification (spherical or plane-parallel) include hydrostatic equilibrium, mixing-length convection and local thermodynamic equilibrium. We discuss a number of general properties of the models, in particular in relation to the effects of changing abundances, of blanketing, and of sphericity. We illustrate positive and negative feedbacks between sphericity and molecular blanketing. We compare the models with those of other available grids and find excellent agreement with planeparallel models of Castelli & Kurucz (if convection is treated consistently) within the overlapping parameter range. Although there are considerable departures from the spherically-symmetric NextGen models, the agreement with more recent PHOENIX models is gratifying. Conclusions. The models of the grid show considerable regularities, but some interesting departures from general patterns occur for the coolest models due to the molecular opacities. We have tested a number of approximate "rules of thumb" concerning effects of blanketing and sphericity and often found them to be astonishingly accurate. Some interesting new phenomena have been discovered and explored, such as the intricate coupling between blanketing and sphericity, and the strong effects of carbon enhancement on metal-poor models. We give further details of line absorption data for molecules, as well as details of models and comparisons with observations in subsequent papers.

A simple model for the absorption of starlight by dust in galaxies

Abstract: We present a new model for computing the effects of dust on the integrated spectral properties of galaxies, based on an idealized description of the main features of the interstellar medium (ISM). The model includes the ionization of H II regions in the interiors of the dense clouds in which stars form and the influence of the finite lifetime of these clouds on the absorption of radiation. We compute the production of emission lines and the absorption of continuum radiation in the H II regions and the subsequent transfer of line and continuum radiation in the surrounding H I regions and the ambient ISM. This enables us to interpret simultaneously all the observations of an ultraviolet-selected sample of nearby starburst galaxies, including the ratio of far-infrared to ultraviolet luminosities, the ratio of Hα to Hβ luminosities, the Hα equivalent width, and the ultraviolet spectral slope. We show that the finite lifetime of stellar birth clouds is a key ingredient for resolving an apparent discrepancy between the attenuation of line and continuum photons in starburst galaxies. In addition, we find that an effective absorption curve proportional to λ-0.7 reproduces the observed relation between the ratio of far-infrared to ultraviolet luminosities and the ultraviolet spectral slope. We interpret this relation most simply as a sequence in the overall dust content of the galaxies. The shallow wavelength dependence of the effective absorption curve is compatible with the steepness of known extinction curves if the dust has a patchy distribution. In particular, we find that a random distribution of discrete clouds with optical depths similar to those in the Milky Way provides a consistent interpretation of all the observations. A noteworthy outcome of our detailed analysis is that the observed mean relations for starburst galaxies can be closely approximated by the following simple recipe: use an effective absorption curve proportional to λ-0.7 to attenuate the line and continuum radiation from each stellar generation, and lower the normalization of the curve, typically by a factor of 3 after 107 yr, to account for the dispersal of the birth clouds. This recipe or our full model for absorption can be incorporated easily into any population synthesis model.

Interstellar abundances from absorption-line observations with the hubble space telescope

Abstract: ▪ Abstract The Goddard High-Resolution Spectrograph (GHRS) aboard the Hubble Space Telescope (HST) has yielded precision abundance results for a range of interstellar environments, including gas in the local medium, in the warm neutral medium, in cold diffuse clouds, and in distant halo clouds. Through GHRS studies, investigators have determined the abundances of elements such as C, N, O, Mg, Si, S, and Fe in individual interstellar clouds. These studies have provided new information about the composition of interstellar dust grains, the origin of the Galactic high-velocity cloud system, and the processes that transport gas between the disk and the halo. Precision measurements of the interstellar D to H ratio and of the abundances of r- and s-process elements have also provided fiducial reference values for cosmological and stellar evolutionary observations and theoretical models.

THE SINS SURVEY: SINFONI INTEGRAL FIELD SPECTROSCOPY OF z ∼ 2 STAR-FORMING GALAXIES*

Abstract: We present the Spectroscopic Imaging survey in the near-infrared (near-IR) with SINFONI (SINS) of high-redshift galaxies. With 80 objects observed and 63 detected in at least one rest-frame optical nebular emission line, mainly Hα, SINS represents the largest survey of spatially resolved gas kinematics, morphologies, and physical properties of star-forming galaxies at z ~ 1-3. We describe the selection of the targets, the observations, and the data reduction. We then focus on the "SINS Hα sample," consisting of 62 rest-UV/optically selected sources at 1.3 < z < 2.6 for which we targeted primarily the Hα and [N II] emission lines. Only ≈30% of this sample had previous near-IR spectroscopic observations. The galaxies were drawn from various imaging surveys with different photometric criteria; as a whole, the SINS Hα sample covers a reasonable representation of massive M_* ≳ 10^(10) M_☉ star-forming galaxies at z ≈ 1.5-2.5, with some bias toward bluer systems compared to pure K-selected samples due to the requirement of secure optical redshift. The sample spans 2 orders of magnitude in stellar mass and in absolute and specific star formation rates, with median values ≈3 × 10^(10) M_☉, ≈70 M_☉ yr^(–1), and ≈3 Gyr^(–1). The ionized gas distribution and kinematics are spatially resolved on scales ranging from ≈1.5 kpc for adaptive optics assisted observations to typically ≈4-5 kpc for seeing-limited data. The Hα morphologies tend to be irregular and/or clumpy. About one-third of the SINS Hα sample galaxies are rotation-dominated yet turbulent disks, another one-third comprises compact and velocity dispersion-dominated objects, and the remaining galaxies are clear interacting/merging systems; the fraction of rotation-dominated systems increases among the more massive part of the sample. The Hα luminosities and equivalent widths suggest on average roughly twice higher dust attenuation toward the H II regions relative to the bulk of the stars, and comparable current and past-averaged star formation rates.