We present new collisional-radiative models (CRMs) for helium in the quiescent solar corona, and predict the emissivities of the He and He$^+$ lines to be observed by DKIST, Solar Orbiter, and Proba-3. We discuss in detail the rates we selected for these models, highlighting several shortcomings we have found in previous work. As no previous complete and self-consistent coronal CRM for helium existed, we have benchmarked our largest model at a density of 10$^{6}$ cm$^{-3}$ and temperature of 20,000 K against recent CRMs developed for photoionised nebulae. We then present results for the outer solar corona, using new dielectronic recombination rates we have calculated, which increase the abundance of neutral helium by about a factor of two. We also find that all the optical triplet He I lines, and in particular the well known He I 10830 and 5876 A lines are strongly affected by both photo-excitation and photo-ionisation from the disk radiation, and that extensive CRM models are required to obtain correct estimates. Close to the Sun, at an electron density of 10$^{8}$ cm$^{-3}$ and temperature of 1 MK, we predict the emissivity of the He I 10830 A to be comparable to that of the strong Fe XIII coronal line at 10798 A. However, we expect the He I emissivity to sharply fall in the outer corona, with respect to Fe XIII. We confirm that the He$^+$ Lyman $\alpha$ at 304 A is also significantly affected by photo-excitation and is expected to be detectable as a strong coronal line up to several solar radii.

We present the results of a calculation of recombination coefficients for O^{2+} + e^- using an intermediate coupling treatment that fully accounts for the dependence of the distribution of population among the ground levels of O^{2+} on electron density and temperature. The calculation is extended down to low electron temperatures where dielectronic recombination arising from Rydberg states converging on the O^{2+} ground levels is an important process. The data, which consist of emission coefficients for 8889 recombination lines and recombination coefficients for the ground and metastable states of O^+ are in Cases A, B and C, and are organised as a function of the electron temperature and number density, as well as wavelength. An interactive fortran 77 data server is also provided as an accessory for mining the line emission coefficients and obtaining Lagrange interpolated values for any choice of the two variables between the explicitly provided values for any set of wavelengths. Some illustrations of the application of the new data to nebular observations are also provided.

In this paper we compute the collision strengths and their thermally-averaged Maxwellian values for electron transitions between the fifteen lowest levels of doubly-ionised cobalt, Co^{2+}, which give rise to forbidden emission lines in the visible and infrared region of spectrum. The calculations also include transition probabilities and predicted relative line emissivities. The data are particularly useful for analysing the thermodynamic conditions of supernova ejecta.

We calculate collision strengths and their thermally-averaged Maxwellian values for electron excitation and de-excitation between the fifteen lowest levels of singly-ionised cobalt, Co+, which give rise to emission lines in the near- and mid-infrared. Transition probabilities are also calculated and relative line intensities predicted for conditions typical of supernova ejecta. The diagnostic potential of the 10.52, 15.46 and 14.74 micro-metre transition lines is briefly discussed.

We present effective collision strengths for electron excitation and de-excitation of the ten forbidden transitions between the five lowest energy levels of the astronomically abundant doubly-ionised oxygen ion, O^{2+}. The raw collision strength data were obtained from an R-matrix intermediate coupling calculation using the Breit-Pauli relativistic approximation published previously by the authors. The effective collision strengths were calculated with kappa-distributed electron energies and are tabulated as a function of the electron temperature and kappa.

We provide a data set of emission and recombination coefficients of hydrogen using a kappa-distribution of electron energies rather than the more traditional Maxwell-Boltzmann (MB) distribution. The data are mainly relevant to thin and relatively cold plasma found in planetary nebulae and H II regions. The data set extends the previous data sets provided by Storey and Hummer which were computed using a MB distribution. The data set, which is placed in the public domain, is structured as a function of electron number density, temperature and kappa. Interactive fortran 77 and C++ data servers are also provided as an accessory to probe the data and obtain Lagrange-interpolated values for any choice of all three variables between the explicitly computed values.

The continuum spectrum of the planetary nebula Hf 2-2 close to the Balmer discontinuity is modeled in the context of the long standing problem of the abundance and temperature discrepancy found when analyzing optical recombination lines and collisionally excited forbidden lines in nebulae. Models are constructed using single and double Maxwell-Boltzmann distributions as well as kappa-distributions for the energies of the free electrons. New results for the necessary continuum and line emission coefficients are presented calculated with kappa-distributed energies. The best fit to the observed continuum spectrum is found to be a model comprising two components with dramatically different temperatures and with a Maxwell-Boltzmann distribution of electron energies. On the basis of a chi-squared analysis, this model is strongly favored over a model with kappa-distributed electron energies.

We present electron collision strengths and their thermally averaged values for the forbidden lines of the astronomically abundant doubly-ionized oxygen ion, O^{2+}, in an intermediate coupling scheme using the Breit-Pauli relativistic terms as implemented in an R-matrix atomic scattering code. We use several atomic targets for the R-matrix scattering calculations including one with 72 atomic terms. We also compare with new results obtained using the intermediate coupling frame transformation method. We find spectroscopically significant differences against a recent Breit-Pauli calculation for the excitation of the [O III] lambda 4363 transition but confirm the results of earlier calculations.

Errors were discovered in our recent calculations of HeI emissivities (Porter et al. 2012; arXiv:1206.4115). This corrigendum discusses the errors and reports corrected emissivities.

The current paper presents atomic data generated to investigate the recombination lines of C II in the spectra of planetary nebulae. These data include energies of bound and autoionizing states, oscillator strengths and radiative transition probabilities, autoionization probabilities, and recombination coefficients. The R-matrix method of electron scattering theory was used to describe the C2+ plus electron system.

A recently generated theoretical line list of C II dielectronic recombination lines together with observational data gathered from the literature is used to investigate the electron temperature in a range of astronomical objects, mainly planetary nebulae. The electron temperature is obtained by a least-squares optimisation using all the reliable observed lines in each object. In addition, the subset of lines arising directly from autoionising states is used to directly determine the free-electron energy distribution which is then compared with various theoretical possibilities. The method described here can potentially determine whether there are departures from Maxwell-Boltzmann distributions in some nebulae, as has been recently proposed. Using published observations of the three planetary nebulae where the relevant lines are recorded, we find that the data are best matched by Maxwell-Boltzmann distributions but that the uncertainties are sufficiently large at present that kappa-distributions or two-component nebular models are not excluded.

We carry out plasma diagnostic analyses for 123 planetary nebulae (PNe) and 42 H II regions using the N II and O II optical recombination lines (ORLs). New effective recombination coefficients for the N II and O II optical recombination spectra are used. These data were calculated under the intermediate coupling scheme for a number of electron temperature (Te) and density (Ne) cases. We used a new method to determine the Te and Ne for the nebular sample, combining the ORLs with the most reliable measurements for each ion and the predicted intensities that are based on the new atomic data. Uncertainties of the derived Te and Ne are estimated for each object. The diagnostic results from heavy element ORLs show reasonable agreement with previous calculations in the literature. We compare the electron temperatures derived from the N II and O II ORLs, Te(ORLs), and those from the collisionally excited lines (CELs), Te(CELs), as well as the hydrogen Balmer jump, Te(H I BJ), especially for the PNe with large abundance discrepancies. Temperatures from He I recombination lines, Te(He I), are also used for comparison if available. For all the objects included in our sample, Te(ORLs) are lower than Te(H I BJ), which are in turn systematically lower than Te(CELs). Nebulae with Te(He I) available show the relation Te(ORLs) < Te(He I) < Te(H I BJ) < Te(CELs), which is consistent with predictions from the bi-abundance nebular model postulated by Liu et al. (2000).

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.

Here we report new ${\it ab initio}$ calculations of the effective recombination coefficients for the \ion{N}{ii} recombination spectrum. We have taken into account the density dependence of the coefficients arising from the relative populations of the fine-structure levels of the ground state of the recombining ion, an elaboration that has not been attempted before for this ion, and it opens up the possibility of electron density determination via recombination line analysis. Photoionization cross-sections, bound state energies, and the oscillator strengths of \ion{N}{ii} with $n \leq 11$ and $l \leq 4$ have been obtained using the close-coupling R-matrix method in the intermediate coupling scheme. Photoionization data were computed that accurately map out the near-threshold resonances and were used to derive recombination coefficients, including radiative and dielectronic recombination. Also new is including the effects of dielectronic recombination via high-$n$ resonances lying between the $^2$P$^{\rm o}$\,$_{1/2}$ and $^2$P$^{\rm o}$\,$_{3/2}$ levels. The new calculations are valid for temperatures down to an unprecedentedly low level (approximately 100 K). The newly calculated effective recombination coefficients allow us to construct plasma diagnostics based on the measured strengths of the \ion{N}{ii} optical recombination lines (ORLs). The derived effective recombination coefficients are fitted with analytic formulae as a function of electron temperature for different electron densities. The dependence of the emissivities of the strongest transitions of \ion{N}{ii} on electron density and temperature is illustrated. Potential applications of the current data to electron density and temperature diagnostics for photoionized gaseous nebulae are discussed. We also present a method of determining electron temperature and density simultaneously.

Effective collision strengths for forbidden transitions among the 5 energetically lowest finestructure levels of O II are calculated in the Breit-Pauli approximation using the R-matrix method. Results are presented for the electron temperature range 100 to 100 000 K. The accuracy of the calculations is evaluated via the use of different types of radial orbital sets and a different configuration expansion basis for the target wavefunctions. A detailed assessment of previous available data is given, and erroneous results are highlighted. Our results reconfirm the validity of the original Seaton and Osterbrock scaling for the optical O II ratio, a matter of some recent controversy. Finally we present plasma diagnostic diagrams using the best collision strengths and transition probabilities.

A number of recent works in astronomy and cosmology have relied upon theoretical He I emissivities, but we know of no effort to quantify the uncertainties in the atomic data. We analyze and assign uncertainties to all relevant atomic data, perform Monte Carlo analyses, and report standard deviations in the line emissivities. We consider two sets of errors, which we call "optimistic" and "pessimistic." We also consider three different conditions, corresponding to prototypical Galactic and extragalactic H II regions and the epoch of cosmological recombination. In the extragalactic H II case, the errors we obtain are comparable to or larger than the errors in some recent $Y_p$ calculations, including those derived from CMB observations. We demonstrate a systematic effect on primordial abundance calculations; this effect cannot be reduced by observing a large number of objects. In the cosmological recombination case, the errors are comparable to many of the effects considered in recent calculations.

We have analysed deep optical spectra of the `born-again' planetary nebula Abell 58 and its hydrogen-deficient knot, surrounding V605 Aql, which underwent a nova-like eruption in 1919. The electron temperature we derive for the central knot varies widely depending on the diagnostic used. The [O III] nebular-to-auroral transition ratio gives a temperature of 20800 K, while the ratio of the [N II] nebular and auroral lines give Te=15200 K. The helium line ratios 5876/4471 and 6678/4471 imply temperatures of 350 K and 550 K respectively. Weakly temperature-sensitive O II recombination line ratios imply similarly low electron temperatures. Abundances derived from recombination lines are vastly higher than those found from collisionally excited lines, with the abundance discrepancy factor (adf) for O2+ reaching 89 -- the second highest known value after that found for the hydrogen deficient knots in Abell 30. The observed temperature diagnostics and abundances support the idea that, like Abell 30, the knot of Abell 58 contains some very cold ionised material. Although the central star is carbon-rich (C/O>1), the knot is found to be oxygen-rich, a situation not predicted by the single-star `born again' theory of its formation. We argue that the abundances in the ejecta observed in A 30 and A 58 have more in common with neon novae than with Sakurai's Object, which is believed to have undergone a final helium flash. In particular, the C/O ratio of less than unity and presence of substantial quantities of neon in the ejecta of both Abell 30 and Abell 58 are not predicted by very late thermal pulse models.

We present coefficients for the calculation of the continuous emission spectra of HI, HeI and HeII due to electron-ion recombination. Coefficients are given for photon energies from the first ionization threshold for each ion to the n=20 threshold of hydrogen (36.5um), and for temperatures 100 K <= Te <=10^5 K. The emission coefficients for HeI are derived from accurate ab initio photoionization data. The coefficients are scaled in such a way that they may be interpolated by a simple scheme with uncertainties less than 1% in the whole temperature and wavelength domain. The data are suitable for incorporation into photoionisation/plasma codes and should aid with the interpretation of spectra from the very cold ionised gas phase inferred to exist in a number of gaseous clouds.

The assumption of spherical symmetry is not justified for the vast majority of PNe. The interpretation of spatially-resolved observations cannot rely solely on the application of 1D codes, which may yield incorrect abundances determinations resulting in misleading conclusions. The 3D photoionisation code MOCASSIN (Monte CAarlo SimulationS of ionised Nebulae) is designed to remedy these shortcomings. The 3D transfer of both primary and secondary radiation is treated self-consistently without the need of approximations. The code was benchmarked and has been applied to the study of several PNe. The current version includes a fully self-consistent radiative transfer treatment for dust grains mixed within the gas, taking into account the microphysics of dust-gas interactions within the geometry-independent Monte Carlo transfer. The new code provides an excellent tool for the self-consistent analysis of dusty ionised regions showing asymmetries and/or density and chemical inhomogeneities. Work is currently in progress to incorporate the processes that dominate the thermal balance of photo-dissociation regions (PDRs), as well as the formation and destruction processes for all the main molecular species.

We have acquired high spectral resolution observations (R=150,000) of the planetary nebulae NGC 7009 and NGC 6153, using bHROS on Gemini South. Observations of this type may provide a key to understanding why optical recombination lines (ORLs) yield systematically higher heavy element abundances for photoionized nebulae than do the classical forbidden collisionally excited lines (CELs) emitted by the same ions; NGC 7009 and NGC 6153 have notably high ORL/CEL abundance discrepancy factors (ADFs) of 5 and 10, respectively. Due to the opposite temperature dependences of ORLs and CELs, ORLs should be preferentially emitted by colder plasma. Our bHROS observations of NGC 7009 reveal that the [O III] 4363A CEL has a FWHM linewidth that is 1.5 times larger than that shown by O II ORLs in the same spectrum, despite the fact that all of these lines are emitted by the O2+ ion. The bHROS spectra of NGC 6153 also show that its O II ORLs have significantly narrower linewidths than do the [O III] 4363A and 5007A lines but, in addition, the [O III] 4363A and 5007A lines show very different velocity profiles, implying the presence of large temperature variations in the nebula.

We present high quality optical spectroscopic observations of the planetary nebula (PN) Hf 2-2. The spectrum exhibits many prominent optical recombination lines (ORLs) from heavy element ions. Analysis of the H {\sc i} and He {\sc i} recombination spectrum yields an electron temperature of $\sim 900$ K, a factor of ten lower than given by the collisionally excited [O {\sc iii}] forbidden lines. The ionic abundances of heavy elements relative to hydrogen derived from ORLs are about a factor of 70 higher than those deduced from collisionally excited lines (CELs) from the same ions, the largest abundance discrepancy factor (adf) ever measured for a PN. By comparing the observed O {\sc ii} $\lambda$4089/$\lambda$4649 ORL ratio to theoretical value as a function of electron temperature, we show that the O {\sc ii} ORLs arise from ionized regions with an electron temperature of only $\sim 630$ K. The current observations thus provide the strongest evidence that the nebula contains another previously unknown component of cold, high metallicity gas, which is too cool to excite any significant optical or UV CELs and is thus invisible via such lines. The existence of such a plasma component in PNe provides a natural solution to the long-standing dichotomy between nebular plasma diagnostics and abundance determinations using CELs on the one hand and ORLs on the other.

We present the first 3D Monte Carlo (MC) photoionisation code to include a fully self-consistent treatment of dust radiative transfer (RT) within a photoionised region. This is the latest development (Version 2.0) of the gas-only photoionisation code MOCASSIN (Ercolano et al., 2003a), and employs a stochastic approach to the transport of radiation, allowing both the primary and secondary components of the radiation field to be treated self-consistently, whilst accounting for the scattering of radiation by dust grains mixed with the gas, as well as the absorption and emission of radiation by both the gas and the dust components. A set of rigorous benchmark tests have been carried out for dust-only spherically symmetric geometries and 2D disk configurations. MOCASSIN's results are found to be in agreement with those obtained by well established dust-only RT codes that employ various approaches to the solution of the RT problem. A model of the dust and of the photoionised gas components of the planetary nebula (PN) NGC 3918 is also presented as a means of testing the correct functioning of the RT procedures in a case where both gas and dust opacities are present. The two components are coupled via the heating of dust grains by the absorption of both UV continuum photons and resonance line photons emitted by the gas. The MOCASSIN results show agreement with those of a 1D dust and gas model of this nebula published previously, showing the reliability of the new code, which can be applied to a variety of astrophysical environments.

The core of planetary nebula NGC 6302 is filled with high-excitation photoionised gas at low expansion velocities. It represents a unique astrophysical situation in which to search for hyperfine structure (HFS) in coronal emission lines from highly ionised species. HFS is otherwise blended by thermal or velocity broadening. Spectra containing [AlVI]~3.66um 3P2-3P1, obtained with Phoenix on Gemini-South at resolving powers of up to 75000, resolve the line into five hyperfine components separated by 20 to 60km/s due to the coupling of the I=5/2 nuclear spin of 27Al with the total electronic angular momentum J. 26Al has a different nuclear spin of I=5, and a different HFS, which allows us to place a 3sigma upper limit on the 26Al/27Al abundance ratio of 1/33. We measure the HFS magnetic-dipole coupling constants for [AlVI], and provide the first estimates of the electric-quadrupole HFS coupling constants obtained through astronomical observations of an atomic transition.

We present the first 3D Monte Carlo photoionisation code to include a fully self-consistent treatment of dust radiative transfer (RT) within the photoionised region. This is the latest development of the recently published pure photoionisation code MOCASSIN (Ercolano et al., 2003a) and it is currently undergoing several benchmarking tests. The preliminary results of these tests are presented in these conference proceedings. The new code provides the ideal tool for the analysis of dusty ionised regions showing asymmetries and/or density and chemical inhomogeneities

[Abridged] Deep optical observations of the spectra of 12 Galactic planetary nebulae (PNe) and 3 Magellanic Cloud PNe were presented in Paper I by Tsamis et al. (2003b), who carried out an abundance analysis using the collisionally excited forbidden lines. Here, the relative intensities of faint optical recombination lines (ORLs) from ions of carbon, nitrogen and oxygen are analysed in order to derive the abundances of these ions relative to hydrogen. We define an abundance discrepancy factor (ADF) as the ratio of the abundance derived for a heavy element ion from its recombination lines to that derived for the same ion from its ultraviolet, optical or infrared collisionally excited lines (CELs). All of the PNe in our sample are found to have ADF's that exceed unity. There is no dependence of the magnitude of the ADF upon the excitation energy of the UV, optical or IR CEL transition used, indicating that classical nebular temperature fluctuations--i.e. in a chemically homogeneous medium--are not the cause of the observed abundance discrepancies. Instead, we conclude that the main cause of the discrepancy is enhanced ORL emission from cold ionized gas located in hydrogen-deficient clumps inside the main body of the nebulae. We have developed a new electron temperature diagnostic, based upon the relative intensities of the OII 4f-3d 4089A and 3p-3s 4649A recombination transitions. For six out of eight PNe for which both transitions are detected, we derive O2+ ORL electron temperatures of <300 K, very much less than the O2+ forbidden-line and Balmer jump temperatures derived for the same nebulae. These results provide direct observational evidence for the presence of H-deficient, cold plasma regions within the nebulae, consistent with gas cooled largely by infrared fine structure and recombination transitions.

A method is presented to derive electron temperatures and densities of planetary nebulae (PNe) simultaneously, using the observed hydrogen recombination spectrum, which includes continuum and line emission. By matching theoretical spectra to observed spectra around the Balmer jump at about 3646A, we determine electron temperatures and densities for 48 Galactic PNe. The electron temperatures based on this method are found to be systematically lower than those derived from [O III] 4959/4363 and [O III] (88\mu+52\mu)/4959 ratios. And the electron densities based on this method are found to be systematically higher than those derived from forbidden line ratios. These results suggest that temperature and density fluctuations are generally present within nebulae. In addition, He/H abundances of PNe are found to be positively correlated with the difference between T_e([O III]) and T_e(Bal), suggesting that He/H abundances might have been overestimated generally because of the possible existence of H-deficient knots.

We have constructed a photoionization model, using the 3-D Monte Carlo code MOCASSIN for one of the hydrogen-deficient knots (J3) of the born-again planetary nebula Abell 30. The model consists of spherical knots, comprising a cold, dense, hydrogen-deficient core with very high metal abundances. The inner core is surrounded by a less dense hydrogen-deficient and metal-enriched gas envelope, with less extreme abundances. The envelope of the knot might have been formed by the mixing of the knot material with the surrounding nebular gas. This bi-chemistry, bi-density model did not produce enough heating to match the fluxes of the collisionally excited emission lines (CELs) and of the optical recombination lines (ORLs) observed in the spectrum of the knot. We therefore included heating by photoelectric emission from dust grains in the thermal equilibrium calculations, and found that dust-to-gas ratios of 0.077 and 0.107 by mass for the central core and the envelope of the knot, respectively, are sufficient to fit the spectrum. Surprisingly, photoelectric emission from grains is the dominant source of heating in the hot envelope of the knot, while heating by photoionization of helium and heavy elements dominates in the cold core. The chemical abundances inferred from our modelling indicate, in agreement with the empirical analysis of Wesson et al. (2003), that the C/O ratio in the two regions of the knot is less than unity, contrary to theoretical predictions for born-again nebulae.

We present spectrophotometry of 12 Galactic and 3 Magellanic Cloud planetary nebulae (PNe). Nine of the Galactic PNe were observed by scanning the slit across the PN. We use the fluxes of collisionally excited lines (CELs) to derive electron densities (D's) and temperatures (T's), and ionic abundances. We find that the D's derived from optical CEL ratios are systematically higher than those derived from the ratios of the IR fine-structure (FS) lines of [OIII], indicating the presence of significant density variations within the PNe. We also compare T's obtained from the ratio of optical nebular to auroral [OIII] lines with those obtained from the ratio of [OIII] optical to IR FS lines. We find that when the latter are derived using D's based on the [OIII] 52um/88um ratio, they yield values that are significantly higher than the optical [OIII] T's. Contrasting this, [OIII] optical/IR T's derived using the higher D's obtained from [ClIII] 5517A/5537A ratios show much closer agreement with optical [OIII] T's, implying that the observed [OIII] optical/IR ratios are significantly weighted by D's in excess of the critical densities of both [OIII] FS lines. Consistent with this, ionic abundances derived from [OIII] and [NIII] FS lines using D's from optical CELs show much better agreement with abundances derived for the same ions from optical and UV CELs than do abundances derived from the FS lines using the lower D's obtained from the 52um/88um ratios. The behaviour of T's obtained making use of the T-insensitive IR FS lines provides no support for significant T-fluctuations within the PNe that could be responsible for derived Balmer jump T's being lower than those obtained from the much more T-sensitive [OIII] optical lines.

The study of photoionized environments is fundamental to many astrophysical problems. Up to the present most photoionization codes have numerically solved the equations of radiative transfer by making the extreme simplifying assumption of spherical symmetry. Unfortunately very few real astronomical nebulae satisfy this requirement. To remedy these shortcomings, a self-consistent, three-dimensional radiative transfer code has been developed using Monte Carlo techniques. The code, Mocassin, is designed to build realistic models of photoionized nebulae having arbitraries geometry and density distributions with both the stellar and diffuse radiation fields treated self-consistently. In addition, the code is capable of tretating on or more exciting stars located at non-central locations. The gaseous region is approximated by a cuboidal Cartesian grid composed of numerous cells. The physical conditions within each grid cell are determined by solving the thermal equilibrium and ionization balance equations This requires a knowledge of the local primary and secondary radiation fields, which are calculated self-consistently by locally simulating the individual processes of ionization and recombination. The main structure and computational methods used in the Mocassin code are described in this paper. Mocassin has been benchmarked against established one-dimensional spherically symmetric codes for a number of standard cases, as defined by the Lexington/Meudon photoionization workshops (Pequignot et al., 1986; Ferland et al., 1995; Pequignot et al., 2001)\citep{pequignot86,ferland95, pequignot01}. The results obtained for the benchmark cases are satisfactory and are presented in this paper. A performance analysis has also been carried out and is discussed here.

The three-dimensional Monte Carlo photoionization code Mocassin has been applied to construct a realistic model of the planetary nebula NGC 3918. Three different geometric models were tried. The effects of the interaction of the diffuse fields coming from two adjacent regions of different densities were investigated. These are found to be non-negligible, even for the relatively uncomplicated case of a biconical geometry. We found that the ionization structure of low ionization species near the boundaries is particularly affected. It is found that all three models provided acceptable matches to the integrated nebular optical and ultraviolet spectrum. Large discrepancies were found between all of the model predictions of infrared fine-structure line fluxes and ISO SWS measurements. This was found to be largely due to an offset of ~14 arcsec from the centre of the nebula that affected all of the ISO observations of NGC 3918. For each model, we also produced projected emission-line maps and position-velocity diagrams from synthetic long-slit spectra, which could be compared to recent HST images and ground-based long-slit echelle spectra. This comparison showed that spindle-like model B provided the best match to the observations. We have therefore shown that although the integrated emission line spectrum of NGC 3918 can be reproduced by all three of the three-dimensional models investigated in this work, the capability of creating projected emission-line maps and position-velocity diagrams from synthetic long-slit spectra was crucial in allowing us to constrain the structure of this object.

We have obtained deep optical, long-slit spectrophotometry of the Galactic HII regions M 17, NGC 3576 and of the Magellanic Cloud HII regions 30 Doradus, LMC N11B and SMC N66, recording the optical recombination lines (ORLs) of CII, NII and OII. Temperature-insensitive ORL C2+/O2+ and N2+/O2 ratios are obtained for all nebulae except SMC N66. The ORL C2+/O2+ ratios show remarkable agreement within each galactic system, while also being in agreement with the corresponding CEL ratios. For all five nebulae, the O2+/H+ abundance derived from multiple OII ORLs is found to be higher than the corresponding value derived from the strong [OIII] 4959, 5007A CELs, by factors of 1.8--2.7 for four of the nebulae. The LMC N11B nebula exhibits a more extreme discrepancy factor for the O2+ ion, ~5. Thus these HII regions exhibit ORL/CEL abundance discrepancy factors that are similar to those previously encountered amongst planetary nebulae. Our optical CEL O2+/H+ abundances agree to within 20-30 per cent with published O2+/H+ abundances that were obtained from observations of infrared fine-structure lines. Since the low excitation energies of the latter make them insensitive to variations about typical nebular temperatures, fluctuations in temperature are ruled out as the cause of the observed ORL/CEL O2+ abundance discrepancies. We present evidence that the observed OII ORLs from these HII regions originate from gas of very similar density (<3500 cm-3) to that emitting the observed heavy-element optical and infrared CELs, ruling out models that employ high-density ionized inclusions in order to explain the abundance discrepancy. We consider a scenario whereby much of the heavy-element ORL emission originates from cold (<=500 K) metal-rich ionized regions.

We present new HST/WFPC2 imagery and STIS long-slit spectroscopy of the planetary nebula NGC 7009. The primary goal was to obtain high spatial resolution of the intrinsic line ratio [O III] 4364/5008 and thereby evaluate the electron temperature (T_e) and the fractional mean-square T_e variation (t_A^2) across the nebula. The WFPC2 T_e map is rather uniform; almost all values are between 9000 - 11,000 K, with the higher T_e's closely coinciding with the inner He^++-zone. The results indicate very small values - <~ 0.01 - for t_A^2 throughout. Our STIS data allow an even more direct determination of T_e and t_A^2, albeit for a much smaller area than with WFPC2. We present results from binning the data along the slit into tiles that are 0.5'' square (matching the slit width). The average [O III] temperature using 45 tiles (excluding the central star and STIS fiducial bars) is 10,139 K; t_A^2 is 0.0035. The measurements of T_e reported here are an average along each line of sight. Therefore, despite finding remarkably low t_A^2, we cannot completely rule out temperature fluctuations along the line of sight as the cause of the large abundance discrepancy between heavy element abundances inferred from collisionally excited emission lines compared to those derived from recombination lines.