The long-standing difference in chemical abundances determined from optical recombination lines and collisionally excited lines raises questions about our understanding of atomic physics, as well as the assumptions made when determining physical conditions and chemical abundances in astrophysical nebulae. Here, we study the recombination contribution of [O III] 4363 and the validity of the line ratio [O III] 4363/4959 as a temperature diagnostic in planetary nebulae with a high abundance discrepancy. We derive a fit for the recombination coefficient of [O III] 4363 that takes into account the radiative and dielectronic recombinations, for electron temperatures from 200 to 30,000 K. We estimate the recombination contribution of [O III] 4363 for the planetary nebulae Abell 46 and NGC 6778 by subtracting the collisional contribution from the total observed flux. We find that the spatial distribution for the estimated recombination contribution in [O III] 4363 follows that of the O II 4649 recombination line, both peaking in the central regions of the nebula, especially in the case of Abell 46 which has a much higher abundance discrepancy. The estimated recombination contribution reaches up to 70% and 40% of the total [O III] 4363 observed flux, for Abell 46 and NGC 6778, respectively.

We present a new catalog of HII regions based on the integral field spectroscopy (IFS) data of the extended CALIFA and PISCO samples. The selection of HII regions was based on two assumptions: a clumpy structure with high contrast of H$\alpha$ emission and an underlying stellar population comprising young stars. The catalog provides the spectroscopic information of 26,408 individual regions corresponding to 924 galaxies, including the flux intensities and equivalent widths of 51 emission lines covering the wavelength range between 3745-7200A. To our knowledge, this is the largest catalog of spectroscopic properties of HII regions. We explore a new approach to decontaminate the emission lines from diffuse ionized gas contribution. This diffuse gas correction was estimated to correct every emission line within the considered spectral range. With the catalog of HII regions corrected, new demarcation lines are proposed for the classical diagnostic diagrams. Finally, we study the properties of the underlying stellar populations of the HII regions. It was found that there is a direct relationship between the ionization conditions on the nebulae and the properties of stellar populations besides the physicals condition on the ionized regions.

We present a detailed analysis of deep VLT/X-Shooter observations of the planetary nebula Tc 1. We calculate gas temperature, density, extinction, and abundances for several species from the empirical analysis of the total line fluxes. In addition, a spatially resolved analysis of the most intense lines provides the distribution of such quantities across the nebula. The new data reveal that several lines exhibit a double peak spectral profile consistent with the blue- and red-shifted components of an expanding spherical shell. The study of such components allowed us to construct for the first time a three-dimensional morphological model, which reveals that Tc 1 is a slightly elongated spheroid with an equatorial density enhancement seen almost pole on. A few bright lines present extended wings (with velocities up to a few hundred km/s), but the mechanism producing them is not clear. We constructed photoionization models for the main shell of Tc 1. The models predict the central star temperature and luminosity, as well as the nebular density and abundances similar to previous studies. Our models indicate that Tc 1 is located at a distance of approximately 2 kpc. We report the first detection of the [Kr III] 6825 A emission line, from which we determine the Krypton abundance. Our model indicates that the main shell of Tc 1 is matter bounded; leaking H ionizing photons may explain the ionization of its faint AGB-remnant halo.

We explore the production and escape of ionizing photons in young galaxies by investigating the ultraviolet and optical emission-line properties of models of ionization-bounded and density-bounded HII regions, active-galactic-nucleus (AGN) narrow-line regions and radiative shocks computed all using the same physically-consistent description of element abundances and depletion on to dust grains down to very low metallicities. We compare these models with a reference sample of metal-poor star-forming galaxies and Lyman-continuum (LyC) leakers at various redshifts, which allows the simultaneous exploration of more spectral diagnostics than typically available at once for individual subsamples. We confirm that current single- and binary-star population synthesis models do not produce hard-enough radiation to account for the high-ionization emission of the most metal-poor galaxies. Introducing either an AGN or radiative-shock component brings models into agreement with observations. A published model including X-ray binaries is an attractive alternative to reproduce the observed rise in HeII4686/Hbeta ratio with decreasing oxygen abundance in metal-poor star-forming galaxies, but not the high observed HeII4686/Hbeta ratios of galaxies with large EW(Hbeta). A source of harder ionizing radiation appears to be required in these extreme objects, such as an AGN or radiative-shock component, perhaps linked to an initial-mass-function bias toward massive stars at low metallicity. This would also account for the surprisingly high [OI]/[OIII] ratios of confirmed LyC leakers relative to ionization-bounded models. We find no simple by-eye diagnostic of the nature of ionizing sources and the escape of LyC photon, which require proper simultaneous fits of several lines to be discriminated against.

SIGNALS, the Star formation, Ionized Gas, and Nebular Abundances Legacy Survey, is a large observing program designed to investigate massive star formation and HII regions in a sample of local extended galaxies. The program will use the imaging Fourier transform spectrograph SITELLE at the Canada-France-Hawaii Telescope. Over 355 hours (54.7 nights) have been allocated beginning in fall 2018 for eight consecutive semesters. Once completed, SIGNALS will provide a statistically reliable laboratory to investigate massive star formation, including over 50 000 resolved HII regions : the largest, most complete, and homogeneous database of spectroscopically and spatially resolved extragalactic HII regions ever assembled. For each field observed, three datacubes covering the spectral bands of the filters SN1 (363 -386 nm), SN2 (482 - 513 nm), and SN3 (647 - 685 nm) are gathered. The spectral resolution selected for each spectral band is 1000, 1000, and 5000, respectively. As defined, the project sample will facilitate the study of small-scale nebular physics and many other phenomena linked to star formation at a mean spatial resolution of 20 pc. This survey also has considerable legacy value for additional topics including planetary nebulae, diffuse ionized gas, andsupernova remnants. The purpose of this paper is to present a general outlook of the survey, notably the observing strategy, galaxy sample, and science requirements.

We present a new database of fully radiative shock models calculated with the shock and photoionization code MAPPINGS V. The database architecture is built to contain diverse shock grids comprising of multiple shock parameters. It can be easily accessible through the MySQL protocol. Intensities of spectral lines from infrared to X-rays are stored along with other useful outputs such as the ionic fractions/temperature, integrated densities, etc. A web page was created in other to explore interactively the database as it evolves with time. Examples of its usage is given using the Python language.

Planetary nebulae are ionized clouds of gas formed by the hydrogen-rich envelopes of low- and intermediate-mass stars ejected at late evolutionary stages. The strong UV flux from their central stars causes a highly stratified ionization structure, with species of higher ionization potential closer to the star. Here we report on the exceptional case of HuBi 1, a double-shell planetary nebula whose inner shell presents emission from low-ionization species close to the star and emission from high-ionization species farther away. Spectral analysis demonstrates that the inner shell of HuBi 1 is excited by shocks, whereas its outer shell is recombining. The anomalous excitation of these shells can be traced to its low-temperature [WC10] central star whose optical brightness has declined continuously by 10 magnitudes in a period of 46 years. Evolutionary models reveal that this star is the descendent of a low-mass star ($\simeq$1.1 $M_\odot$) that has experienced a born-again event whose ejecta shock-excite the inner shell. HuBi 1 represents the missing link in the formation of metal-rich central stars of planetary nebulae from low-mass progenitors, offering unique insight regarding the future evolution of the born-again Sakurai's object. Coming from a solar-mass progenitor, HuBi 1 represents a potential end-state for our Sun.

We investigated the infrared (IR) dust emission from PN IC 418, using a detailed model controlled by a previous determination of the stellar properties and the characteristics of the photoionized nebula, keeping as free parameters the dust types, amounts and distributions relative to the distance of the central star. The model includes the ionized region and the neutral region beyond the recombination front (Photodissociation region, or PDR), where the [OI] and [CII] IR lines are formed. We succeeded in reproducing the observed infrared emission from 2 to 200~\mm. The global energy budget is fitted by summing up contributions from big grains of amorphous carbon located in the neutral region and small graphite grains located in the ionized region (closer to the central star). Two emission features seen at 11.5 and 30~\mm are also reproduced by assuming them to be due to silicon carbide (SiC) and magnesium and iron sulfides (Mg$_x$Fe$_{1-x}$S), respectively. For this, we needed to consider ellipsoidal shapes for the grains to reproduce the wavelength distribution of the features. Some elements are depleted in the gaseous phase: Mg, Si, and S have sub-solar abundances (-0.5 dex below solar by mass), while the abundance of C+N+O+Ne by mass is close to solar. Adding the abundances of the elements present in the dusty and gaseous forms leads to values closer to but not higher than solar, confirming that the identification of the feature carriers is plausible. Iron is strongly depleted (3 dex below solar) and the small amount present in dust in our model is far from being enough to recover the solar value. A remaining feature is found as a residue of the fitting process, between 12 and 25~\mm, for which we do not have identification.

The morphology and optical spectrum of IPHASXJ210205+471015, a nebula classified as a possible planetary nebula, are however strikingly similar to those of ATCnc, a classical nova shell around a dwarf nova. To investigate its true nature, we have obtained high-resolution narrow-band [O III] and [N II] images and deep GTC OSIRIS optical spectra. The nebula shows an arc of [N II]-bright knots notably enriched in nitrogen, whilst an [O III]-bright bow-shock is progressing throughout the ISM. Diagnostic line ratios indicate that shocks are associated with the arc and bow-shock. The central star of this nebula has been identified by its photometric variability. Time-resolved photometric and spectroscopic data of this source reveal a period of 4.26 hours, which is attributed to a binary system. The optical spectrum is notably similar to that of RWSex, a cataclysmic variable star (CV) of the UXUMa nova-like (NL) type. Based on these results, we propose that IPHASX J210205+471015 is a classical nova shell observed around a CV-NL system in quiescence.

We report deep spectroscopy of ten planetary nebulae (PNe) in the Andromeda Galaxy (M31) using the 10.4m GTC. Our targets reside in different regions of M31, including halo streams and dwarf satellite M32, and kinematically deviate from the extended disk. The temperature-sensitive [O III] 4363 auroral line is observed in all targets. For four PNe, the GTC spectra extend beyond 1 micron, enabling explicit detection of the [S III] 6312 and 9069,9531 lines and thus determination of the [S III] temperature. Abundance ratios are derived and generally consistent with AGB model predictions. Our PNe probably all evolved from low-mass (<2 M_sun) stars, as analyzed with the most up-to-date post-AGB evolutionary models, and their main-sequence ages are mostly ~2-5 Gyr. Compared to the underlying, smooth, metal-poor halo of M31, our targets are uniformly metal-rich ([O/H]> -0.4), and seem to resemble the younger population in the stream. We thus speculate that our halo PNe formed in the Giant Stream's progenitor through extended star formation. Alternatively, they might have formed from the same metal-rich gas as did the outer-disk PNe, but was displaced into their present locations as a result of galactic interactions. These interpretations are, although speculative, qualitatively in line with the current picture, as inferred from previous wide-field photometric surveys, that M31's halo is the result of complex interactions and merger processes. The behavior of N/O of the combined sample of the outer-disk and our halo/substructure PNe signifies that hot bottom burning might actually occur at <3 M_sun, but careful assessment is needed.

We identify [Se III] 1.0994 micron in the planetary nebula (PN) NGC 5315 and [Kr VI] 1.2330 micron in three PNe, from spectra obtained with the FIRE spectrometer on the 6.5-m Baade Telescope. Se and Kr are the two most widely-detected neutron-capture elements in astrophysical nebulae, and can be enriched by s-process nucleosynthesis in PN progenitor stars. The detection of [Se III] 1.0994 micron is particularly valuable when paired with observations of [Se IV] 2.2858 micron, as it can be used to improve the accuracy of nebular Se abundance determinations, and allows Se ionization correction factor (ICF) schemes to be empirically tested for the first time. We present new effective collision strength calculations for Se^{2+} and Kr^{5+}, which we use to compute ionic abundances. In NGC 5315, we find that the Se abundance computed from Se^{3+}/H^+ is lower than that determined with ICFs that incorporate Se^{2+}/H^+. We compute new Kr ICFs that take Kr^{5+}/H^+ into account, by fitting correlations found in grids of Cloudy models between Kr ionic fractions and those of more abundant elements, and use these to derive Kr abundances in four PNe. Observations of [Se III] and [Kr VI] in a larger sample of PNe, with a range of excitation levels, are needed to rigorously test the ICF prescriptions for Se and our new Kr ICFs.

We present the first part of an investigation on the planetary nebula (PN) PC 22 which focuses on the use of deep imaging and high resolution echelle spectroscopy to perform a detailed morpho-kinematical analysis. PC 22 is revealed to be a multipolar PN emitting predominantly in [O III] and displaying multiple non-symmetric outflows. Its central region is found to be also particularly inhomogeneous with a series of low ionization structures (knots) located on the path of the outflows. The morpho-kinematical model obtained with SHAPE indicates that i) the de-projected velocities of the outflows are rather large, > 100 km/s, while the central region has expansion velocities in the range ~25 to ~45 km/s following the "Wilson effect", ii) the majority of the measured structures share similar inclination, ~100 degrees, i.e. they are coplanar, and iii) all outflows and lobes are coeval (within the uncertainties). All these results make us to suggest that PC 22 is an evolved starfish PN. We propose that the mechanism responsible for the morphology of PC 22 consists of a wind-shell interaction, where the fast post-AGB wind flows through a filamentary AGB shell with some large voids.

We study the galactic wind in the edge-on spiral galaxy UGC 10043 with the combination of the CALIFA integral field spectroscopy data, scanning Fabry-Perot interferometry (FPI), and multiband photometry. We detect ionized gas in the extraplanar regions reaching a relatively high distance, up to ~ 4 kpc above the galactic disk. The ionized gas line ratios ([N ii]/Ha, [S ii]/Ha and [O i]/Ha) present an enhancement along the semi minor axis, in contrast with the values found at the disk, where they are compatible with ionization due to H ii-regions. These differences, together with the biconic symmetry of the extra-planar ionized structure, makes UGC 10043 a clear candidate for a galaxy with gas outflows ionizated by shocks. From the comparison of shock models with the observed line ratios, and the kinematics observed from the FPI data, we constrain the physical properties of the observed outflow. The data are compatible with a velocity increase of the gas along the extraplanar distances up to < 400 km/s and the preshock density decreasing in the same direction. We also observe a discrepancy in the SFR estimated based on Ha (0.36 Msun/yr ) and the estimated with the CIGALE code, being the latter 5 times larger. Nevertheless, this SFR is still not enough to drive the observed galactic wind if we do not take into account the filling factor. We stress that the combination of the three techniques of observation with models is a powerful tool to explore galactic winds in the Local Universe.

We present a quick-start guide to BOND, a statistical method to derive oxygen and nitrogen abundances in H II regions. BOND compares a set of carefully selected strong and semistrong emission lines to a grid photoionization models. The first novelty, in comparison to other statistical methods, is that BOND relies on the [Ar III]/[Ne III] emission line ratio to break the oxygen abundance bimodality. In doing so, we can measure oxygen and nitrogen abundances without assuming any a priori relation between N/O and O/H. The second novelty is that BOND takes into account changes in the hardness of the ionizing radiation field, which can come about due to the ageing of H II regions or the stochastically sampling of the IMF. We use the emission line ratio He I/Hb, in addition to commonly used strong lines, to constrain the hardness of the ionizing radiation field. Finally, we also stress the pragmatic considerations behind our Bayesian inference.

The understanding of astronomical nebulae is based on observational data (images, spectra, 3D data-cubes) and theoretical models. In this review, I present my very biased view on photoionization modeling of planetary nebulae, focusing on 1D multi-component models, on 3D models and on big database of models.

I present recent and forthcoming works to model the CALIFA HII region using photoionization models. The first results are obtained with ad-hoc models (combining parameter determination by model fitting and strong line methods) while the next ones will use a Genetic Algorithm to fit the observations in a multi-dimensional space.

Photoionization models of HII regions require as input a description of the ionizing SED and of the gas distribution, in terms of ionization parameter U and chemical abundances (e.g. O/H and N/O). A strong degeneracy exists between the hardness of the SED and U, which in turn leads to high uncertainties in the determination of the other parameters, including abundances. One way to resolve the degeneracy is to fix one of the parameters using additional information. For each of the ~ 20000 sources of the CALIFA HII regions catalog, a grid of photoionization models is computed assuming the ionizing SED being described by the underlying stellar population obtained from spectral synthesis modeling. The ionizing SED is then defined as the sum of various stellar bursts of different ages and metallicities. This solves the degeneracy between the shape of the ionizing SED and U. The nebular metallicity (associated to O/H) is defined using the classical strong line method O3N2 (which gives to our models the status of "hybrids"). The remaining free parameters are the abundance ratio N/O and the ionization parameter U, which are determined by looking for the model fitting [NII]/Ha and [OIII]/Hb. The models are also selected to fit [OII]/Hb. This process leads to a set of ~ 3200 models that reproduce simultaneously the three observations. We find that the regions associated to young stellar bursts suffer leaking of the ionizing photons, the proportion of escaping photons having a median of 80\%. The set of photoionization models satisfactorily reproduces the electron temperature derived from the [OIII]4363/5007 line ratio. We determine new relations between the ionization parameter U and the [OII]/[OIII] or [SII]/[SIII] line ratios. New relations between N/O and O/H and between U and O/H are also determined. All the models are publicly available on the 3MdB database.

We present deep spectra of ring nebulae associated with Wolf-Rayet (WR) and O-type stars: NGC 6888, G2.4+1.4, RCW 58, S 308, NGC 7635 and RCW 52. The data have been taken with the 10m Gran Telescopio Canarias and the 6.5m Clay Telescope. We extract spectra of several apertures in some of the objects. We derive C$^{++}$ and O$^{++}$ abundances from faint recombination lines in NGC 6888 and NGC 7635, permitting to derive their C/H and C/O ratios and estimate the abundance discrepancy factor (ADF) of O$^{++}$. The ADFs are larger than the typical ones of normal HII regions but similar to those found in the ionised gas of star-forming dwarf galaxies. We find that chemical abundances are rather homogeneous in the nebulae where we have spectra of several apertures: NGC 6888, NGC 7635 and G2.4+1.4. We obtain very high values of electron temperature in a peripheral zone of NGC 6888, finding that shock excitation can reproduce its spectral properties. We find that all the objects associated with WR stars show N enrichment. Some of them also show He enrichment and O deficiency as well as a lower Ne/O than expected, this may indicate the strong action of the ON and NeNa cycles. We have compared the chemical composition of NGC 6888, G2.4+1.4, RCW 58 and S 308 with the nucleosynthesis predicted by stellar evolution models of massive stars. We find that non-rotational models of stars of initial masses between 25 and 40 solar masses seem to reproduce the observed abundance ratios of most of the nebulae.

We present BOND, a Bayesian code to simultaneously derive oxygen and nitrogen abundances in giant H II regions. It compares observed emission lines to a grid of photoionization models without assuming any relation between O/H and N/O. Our grid spans a wide range in O/H, N/O and ionization parameter U, and covers different starburst ages and nebular geometries. Varying starburst ages accounts for variations in the ionizing radiation field hardness, which arise due to the ageing of H II regions or the stochastic sampling of the initial mass function. All previous approaches assume a strict relation between the ionizing field and metallicity. The other novelty is extracting information on the nebular physics from semi-strong emission lines. While strong lines ratios alone ([O III]/Hbeta, [O II]/Hbeta and [N II]/Hbeta) lead to multiple O/H solutions, the simultaneous use of [Ar III]/[Ne III] allows one to decide whether an H II region is of high or low metallicity. Adding He I/Hbeta pins down the hardness of the radiation field. We apply our method to H II regions and blue compact dwarf galaxies, and find that the resulting N/O vs O/H relation is as scattered as the one obtained from the temperature-based method. As in previous strong-line methods calibrated on photoionization models, the BOND O/H values are generally higher than temperature-based ones, which might indicate the presence of temperature fluctuations or kappa distributions in real nebulae, or a too soft ionizing radiation field in the models.

The $^3$He isotope is important to many fields of astrophysics, including stellar evolution, chemical evolution, and cosmology. The isotope is produced in low-mass stars which evolve through the planetary nebula (PN) phase. $^3$He abundances in PNe can help test models of the chemical evolution of the Galaxy. We present the detection of the $^3$He$^+$ emission line using the single dish Deep Space Station 63, towards the PN IC$\,$418. We derived a $^3$He/H abundance in the range 1.74$\pm$0.8$\times$10$^{-3}$ to 5.8$\pm$1.7$\times$10$^{-3}$, depending on whether part of the line arises in an outer ionized halo. The lower value for $^3$He/H ratio approaches values predicted by stellar models which include thermohaline mixing, but requires that large amounts of $^3$He are produced inside low-mass stars which enrich the interstellar medium (ISM). However, this over-predicts the $^3$He abundance in HII regions, the ISM, and proto-solar grains, which is known to be of the order of 10$^{-5}$. This discrepancy questions our understanding of the evolution of the $^3$He, from circumstellar environments to the ISM.

The chemical content of the planetary nebula NGC 3918 is investigated through deep, high-resolution (R~40000) UVES at VLT spectrophotometric data. We identify and measure more than 750 emission lines, making ours one of the deepest spectra ever taken for a planetary nebula. Among these lines we detect very faint lines of several neutron-capture elements (Se, Kr, Rb, and Xe), which enable us to compute their chemical abundances with unprecedented accuracy, thus constraining the efficiency of the s-process and convective dredge-up in the progenitor star of NGC 3918.

We apply an axially symmetric pseudo-3D photoionization model, pyCloudy, to derive the structures of 6 bipolar nebulae and 2 suggested post-bipolars in a quest to constrain the bipolar planetary nebulae evolution. HST images and VLT/UVES spectroscopy are used for the modelling. The targets are located in the direction of the Galactic bulge. A 3D model structure is used as input to the photoionization code, so as to fit the HST images. Line profiles of different ions constrain the velocity field. The model and associated velocity fields allow us to derive masses, velocities, and ages. The 3D models find much lower ionized masses than required in 1D models: ionized masses are reduced by factors of 2-7. The selected bi-lobed planetary nebulae show a narrow range of ages: the averaged radii and velocities result in values between 1300 and 2000 yr. The lobes are fitted well with velocities linearly increasing with radius. These Hubble-type flows have been found before, and suggest that the lobes form at a defined point in time. The lobes appear to be slightly younger than the main (host) nebulae, by ~500 yr, they seem to form at an early phase of PN evolution, and fade after 1-2 kyr. We find that 30-35% of bulge PNe pass through a bipolar phase.

We analyse whether a stellar atmosphere model computed with the code CMFGEN provides an optimal description of the stellar observations of WR 136 and simultaneously reproduces the nebular observations of NGC 6888, such as the ionization degree, which is modelled with the pyCloudy code. All the observational material available (far and near UV and optical spectra) were used to constrain such models. We found that even when the stellar luminosity and the mass-loss rate were well constrained, the stellar temperature T_* at tau = 20, can be in a range between 70 000 and 110 000 K. When using the nebula as an additional restriction we found that the stellar models with T_* \sim 70 000 K represent the best solution for both, the star and the nebula. Results from the photoionization model show that if we consider a chemically homogeneous nebula, the observed N^+/O^+ ratios found in different nebular zones can be reproduced, therefore it is not necessary to assume a chemical inhomogeneous nebula. Our work shows the importance of calculating coherent models including stellar and nebular constraints. This allowed us to determine, in a consistent way, all the physical parameters of both the star and its associated nebula. The chemical abundances derived are 12 + log(N/H) = 9.95, 12 + log(C/H) = 7.84 and 12 + log(O/H) = 8.76 for the star and 12 + log(N/H) = 8.40, 12 + log(C/H) = 8.86 and 12 + log(O/H) = 8.20. Thus the star and the nebula are largely N- and C- enriched and O-depleted.

(Abridged) The chemical content of the planetary nebula NGC 3918 is investigated through deep, high-resolution UVES at VLT spectrophotometric data. We identify and measure more than 750 emission lines, making ours one of the deepest spectra ever taken for a planetary nebula. Among these lines we detect very faint lines of several neutron-capture elements (Se, Kr, Rb, and Xe), which enable us to compute their chemical abundances with unprecedented accuracy, thus constraining the efficiency of the s-process and convective dredge-up in the progenitor star of NGC 3918. We find that Kr is strongly enriched in NGC 3918 and that Se is less enriched than Kr, in agreement with the results of previous papers and with predicted s-process nucleosynthesis. We also find that Xe is not as enriched by the s-process in NGC 3918 as is Kr and, therefore, that neutron exposure is typical of modestly sub-solar metallicity AGB stars. A clear correlation is found when representing [Kr/O] vs. log(C/O) for NGC 3918 and other objects with detection of multiple ions of Kr in optical data, confirming that carbon is brought to the surface of AGB stars along with s-processed material during third dredge-up episodes, as predicted by nucleosynthesis models. We also detect numerous refractory element lines (Ca, K, Cr, Mn, Fe, Co, Ni, and Cu). We compute physical conditions from a large number of diagnostics. Thanks to the high ionization of NGC 3918 we detect a large number of recombination lines of multiple ionization stages of C, N, O and Ne. The abundances obtained for these elements by using recently-determined state-of-the-art ICF schemes or simply adding ionic abundances are in very good agreement, demonstrating the quality of the recent ICF scheme for high ionization planetary nebulae.

The possibility that star-forming galaxies may leak ionizing photons is at the heart of many present-day studies that investigate the reionization of the Universe. We test this hypothesis on local blue compact dwarf galaxies of very high excitation. We assembled a sample of such galaxies by examining the spectra from Data Releases 7 and 10 of the Sloan Digital Sky Survey. We argue that reliable conclusions cannot be based on strong lines alone, and adopt a strategy that includes important weak lines such as [OI] and the high-excitation HeII and [ArIV] lines. Our analysis is based on purely observational diagrams and on a comparison of photoionization models with well-chosen emission-line ratio diagrams. We show that spectral energy distributions from current stellar population synthesis models cannot account for all the observational constraints, which led us to mimick several scenarios that could explain the data. These include the additional presence of hard X-rays or of shocks. We find that only ionization-bounded models (or models with an escape fraction of ionizing photons lower than 10%) are able to simultaneously explain all the observational constraints.

We study the relation between the chemical composition and the type of dust present in a group of 20 Galactic planetary nebulae (PNe) that have high quality optical and infrared spectra. The optical spectra are used, together with the best available ionization correction factors, to calculate the abundances of Ar, C, Cl, He, N, Ne, and O relative to H. The infrared spectra are used to classify the PNe in two groups depending on whether the observed dust features are representative of oxygen-rich or carbon-rich environments. The sample contains one object from the halo, eight from the bulge, and eleven from the local disc. We compare their chemical abundances with nucleosynthesis model predictions and with the ones obtained in seven Galactic H II regions of the solar neighbourhood. We find evidence of O enrichment (by $\sim$ 0.3 dex) in all but one of the PNe with carbon-rich dust (CRD). Our analysis shows that Ar, and especially Cl, are the best metallicity indicators of the progenitors of PNe. There is a tight correlation between the abundances of Ar and Cl in all the objects, in agreement with a lockstep evolution of both elements. The range of metallicities implied by the Cl abundances covers one order of magnitude and we find significant differences in the initial masses and metallicities of the PNe with CRD and oxygen-rich dust (ORD). The PNe with CRD tend to have intermediate masses and low metallicities, whereas most of the PNe with ORD show higher enrichments in N and He, suggesting that they had high-mass progenitors.

Photoionization models obtained with numerical codes are widely used to study the physics of the interstellar medium (Planetary Nebulae, H II regions, etc). Grid of models are performed to understand what are the effects of the different parameters used to describe the regions on the observables (mainly emission line intensities). Most of the time, only a small part of the computed results of such grids are published, and they are sometimes hard to obtain in a user-friendly format. We present here the Mexican Million Models dataBase (3MdB), an effort of resolving both of these issues in the form of a database of photoionization models, easily accessible throught the MySQL protocol, and containing a lot of usefull outputs from the models, such as the intensities of 178 emission lines, the ionic fractions of all the ions, etc. Some examples of the use of the 3MdB are also presented.

We present deep, high-resolution (R~40000) UVES at VLT spectrophotometric data of the planetary nebula NGC 3918. This is one of the deepest spectra ever taken of a planetary nebula. We have identified and measured more than 700 emission lines and, in particular, we have detected very faint lines of several neutron-capture elements (s-process elements: Kr, Xe and Rb) that enable us to compute their chemical abundances with unprecedented accuracy, thus constraining the efficiency of the s-process and convective dredge-up.

Analysis of emission lines in gaseous nebulae yields direct measures of physical conditions and chemical abundances and is the cornerstone of nebular astrophysics. Although the physical problem is conceptually simple, its practical complexity can be overwhelming since the amount of data to be analyzed steadily increases; furthermore, results depend crucially on the input atomic data, whose determination also improves each year. To address these challenges we created PyNeb, an innovative code for analyzing emission lines. PyNeb computes physical conditions and ionic and elemental abundances, and produces both theoretical and observational diagnostic plots. It is designed to be portable, modular, and largely customizable in aspects such as the atomic data used, the format of the observational data to be analyzed, and the graphical output. It gives full access to the intermediate quantities of the calculation, making it possible to write scripts tailored to the specific type of analysis one wants to carry out. In the case of collisionally excited lines, PyNeb works by solving the equilibrium equations for an n-level atom; in the case of recombination lines, it works by interpolation in emissivity tables. The code offers a choice of extinction laws and ionization correction factors, which can be complemented by user-provided recipes. It is entirely written in the python programming language and uses standard python libraries. It is fully vectorized, making it apt for analyzing huge amounts of data. The code is stable and has been benchmarked against IRAF/NEBULAR. It is public, fully documented, and has already been satisfactorily used in a number of published papers.

This paper describes the Second Public Data Release (DR2) of the Calar Alto Legacy Integral Field Area (CALIFA) survey. The data for 200 objects are made public, including the 100 galaxies of the First Public Data Release (DR1). Data were obtained with the integral-field spectrograph PMAS/PPak mounted on the 3.5 m telescope at the Calar Alto observatory. Two different spectral setups are available for each galaxy, (i) a low-resolution V500 setup covering the wavelength range 3745-7500 \AA\ with a spectral resolution of 6.0 \AA\ (FWHM), and (ii) a medium-resolution V1200 setup covering the wavelength range 3650-4840 \AA\ with a spectral resolution of 2.3 \AA\ (FWHM). The sample covers a redshift range between 0.005 and 0.03, with a wide range of properties in the Color-Magnitude diagram, stellar mass, ionization conditions, and morphological types. All released cubes were reduced with the latest pipeline, including improved spectrophotometric calibration, spatial registration and spatial resolution. The spectrophotometric calibration is better than 6% and the median spatial resolution is 2.5". Altogether the second data release contains over 1.5 million spectra. It is available at http://califa.caha.es/DR2.

We use optical integral field spectroscopy and 8 and 24 micron mid-IR observations of the giant HII region NGC 588 in the disc of M33 as input and constraints for two-dimensional tailor-made photoionisation models. Two different geometrical approaches are followed for the modelling structure: i) Each spatial element of the emitting gas is studied individually using models which assume that the ionisation structure is complete in each element to look for azimuthal variations across gas and dust. ii) A single model is considered, and the two-dimensional structure of the gas and the dust are assumed to be due to the projection of an emitting sphere onto the sky. The models in both assumptions reproduce the radial profiles of Hbeta surface brightness, the observed number of ionising photons, and the strong optical emission-line relative intensities. The first approach produces a constant-density matter-bounded thin shell of variable thickness and dust-to-gas ratio, while the second gives place to a radiation-bounded thick shell sphere of decreasing particle density. However, the radial profile of the 8/24 microns IR ratio, depending on the gas and dust geometry, only fits well when the thick-shell model is used. The resulting dust-to-gas mass ratio, which was obtained empirically from the derived dust mass using data from Spitzer, also has a better fit using the thick-shell solution. In both approaches, models support the chemical homogeneity, and the ionisation-parameter radial decrease, These results must be taken with care in view of the very low extinction values that are derived from the IR, as compared to that derived from the Balmer decrement. Besides, the IR can be possibly contaminated with the emission from a cloud of diffuse gas and dust above the plane of the galaxy detected at 250 micron Herschel image.

We present new results on the chemical composition of the Galactic ring nebula NGC6888 surrounding the WN6(h) star WR136. The data are based on deep spectroscopical observations taken with the High Dispersion Spectrograph at the 8.2m Subaru Telescope. The spectra cover the optical range from 3700 to 7400 A. The effect of the CNO cycle is well identified in the abundances of He, N, and O, while elements not involved in the synthesis such as Ar, S, and Fe present values consistent with the solar vicinity and the ambient gas. The major achievement of this work is the first detection of the faint CII 4267 recombination line in a Wolf-Rayet nebula. This allows to estimate the C abundance in NGC6888 and therefore investigate for the first time the trace of the CNO cycle in a ring nebula around a Wolf-Rayet star. Although the detection of the CII line has a low signal-to-noise ratio, the C abundance seems to be higher than the predictions of recent stellar evolution models of massive stars. The Ne abundance also show a puzzling pattern with an abundance of about 0.5 dex lower than the solar vicinity, which may be related to the action of the NeNa cycle. Attending to the constraints imposed by the dynamical timescale and the He/H and N/O ratios of the nebula, the comparison with stellar evolution models indicates that the initial mass of the stellar progenitor of NGC6888 is between 25 Msun and 40 Msun.

We compute a large grid of photoionization models that covers a wide range of physical parameters and is representative of most of the observed PNe. Using this grid, we derive new formulae for the ionization correction factors (ICFs) of He, O, N, Ne, S, Ar, Cl, and C. Analytical expressions to estimate the uncertainties arising from our ICFs are also provided. This should be useful since these uncertainties are usually not considered when estimating the error bars in element abundances. Our ICFs are valid over a variety of assumptions such as the input metallicities, the spectral energy distribution of the ionizing source, the gas distribution, or the presence of dust grains. Besides, the ICFs are adequate both for large aperture observations and for pencil-beam observations in the central zones of the nebulae. We test our ICFs on a large sample of observed PNe that extends as far as possible in ionization, central star temperature, and metallicity, by checking that the Ne/O, S/O, Ar/O, and Cl/O ratios show no trend with the degree of ionization. Our ICFs lead to significant differences in the derived abundance ratios as compared with previous determinations, especially for N/O, Ne/O, and Ar/O.

(Abridged) We present the abundance analysis of 12 PNe ionized by [WC]-type stars and wels obtained from high-resolution spectrophotometric data. Our main aims are to determine the chemical composition of the PNe and to study the behaviour of the abundance discrepancy problem (ADF) in this type of planetary nebulae. The detection of a large number of optical recombination lines (ORLs) and collisionally excited lines (CELs) from different ions were presented previously. Most of the ORLs were reported for the first time in these PNe. Ionic abundances were derived from the available CELs and ORLs, using previously determined physical conditions. Based on these two sets of ionic abundances, we derived the total chemical abundances in the nebulae using suitable ICFs (when available). In spite of the [WC] nature of the central stars, moderate ADF(O^++), in the range from 1.2 to 4, were found for all the objects. We found that when the quality of the spectra is high enough the ORLs O^++/H^+ abundance ratios obtained from different multiplets excited mainly by recombination are very similar. Possible dependence of ADFs with some nebular characteristics were analysed, finding no correlation. Abundances derived from CELs were corrected by determining the t^2 parameter. O abundances for PNe, derived from ORLs, are in general larger than the solar abundance. We derived the C/O ratio from ORLs and N/O and alpha-element/O ratios from CELs and found that these PNe are, in average, N-and C-richer than the average of large PN samples. About half of our sample is C-rich (C/O>1). The alpha-elements grow in lockstep with O abundance. Comparing the N/O and C /O ratios with those derived from stellar evolution models, we estimate that about half of our PNe have progenitors with initial masses > 4 M_sun. No correlation was found between the stellar [WC]-type and the nebular abundances.

While collisionally excited lines in HII regions allow one to easily probe the chemical composition of the interstellar medium in galaxies, the possible presence of important temperature fluctuations casts some doubt on the derived abundances. To provide new insights into this question, we have carried out a detailed study of a giant HII region, H 1013, located in the galaxy M101, for which many observational data exist and which has been claimed to harbour temperature fluctuations at a level of t^2 = 0.03-0.06. We have first complemented the already available optical observational datasets with a mid-infrared spectrum obtained with the Spitzer Space Telescope. Combined with optical data, this spectrum provides unprecedented information on the temperature structure of this giant HII region. A preliminary analysis based on empirical temperature diagnostics suggests that temperature fluctuations should be quite weak. We have then performed a detailed modelling using the pyCloudy package based on the photoionization code Cloudy. We have been able to produce photoionization models constrained by the observed Hb surface brightness distribution and by the known properties of the ionizing stellar population than can account for most of the line ratios within their uncertainties. Since the observational constraints are both strong and numerous, this argues against the presence of significant temperature fluctuations in H 1013. The oxygen abundance of our best model is 12 + log O/H = 8.57, as opposed to the values of 8.73 and 8.93 advocated by Esteban et al. (2009) and Bresolin (2007), respectively, based on the significant temperature fluctuations they derived. However, our model is not able to reproduce the intensities of the oxygen recombination lines . This cannot be attributed to observational uncertainties and requires an explanation other than temperature fluctuations.

An analysis of the OIII and SIII temperatures measurements compiled by Perez-Montero et al. of emission line objects consisting of HII galaxies, giant extragalactic HII regions, Galactic HII regions and HII regions from the Magellanic Clouds, reveals that the OIII temperatures are higher than the corresponding values from SIII in most objects with gas metallicities in excess of 0.2 solar. We explore the possibility of temperature inhomogeneities. We also explored metallicity inhomogeneities by combining two models of widely different metallicity. We calculate models that consider a non-Mawell-Boltzmann distributions for the electron energies (kappa parametrization). We also consider shock heating within the photoionized nebula. Varying the various input parameters in the pure photoionization case does not reproduce the observed nebular temperatures and neither does having local temperature inhomogeneities. We find that (1) metallicity inhomogeneities of the nebular gas, (2) shock waves of velocities < 60 km/s propagating in a photoionized plasma, and (3) an electron energy distribution given by a kappa-distribution, are successful in reproducing the observed excess in the [OIII] temperatures. Shock models, however, would require proper 3D hydrodynamical simulations to become a fully developed alternative while models with metallicity inhomogeneities appear to fail in metal poor nebulae, since they result in T_O++(rec)> T_OIII.

We construct a detailed model for the IR dust emission from the PN IC 418. We succeed to reproduce the emission from 2 to 200$\mu$m. We can determine the amount of emitting dust as well as its composition, and compare to the depletion of elements determined for the photoionized region.

We compare calculated intensities of lines of CII, NI, NII, OI and OII with a published deep spectroscopic survey of planetary nebula IC 418. Our calculations use a self-consistent nebular model and a synthetic spectrum of the central star atmosphere to take into account line excitation by continuum fluorescence and electron recombination. We found that the NII spectrum of s, p and most d states is excited by fluorescence due to the low excitation conditions of the nebula. Many CII and OII lines have significant amounts of excitation by fluorescence. Recombination excites all the lines from f and g states and most OII lines. In the neutral-ionized boundary the NI quartet and OI triplet dipole allowed lines are excited by fluorescence, while the quintet OI lines are excited by recombination. Electron excitation produces the forbidden optical lines of OI, and continuum fluorescence enhances the NI forbidden line intensities. Lines excited by fluorescence of light below the Lyman limit thus suggest a new diagnostic to explore the inner boundary of the photodissociation region of the nebula.

(Abridged) Planetary nebulae (PNe) around Wolf-Rayet [WR] central stars ([WR]PNe) constitute a particular photoionized nebula class that represents about 10% of the PNe with classified central stars. We analyse deep high-resolution spectrophotometric data of 12 [WR]PNe. This sample represents the most extensive analysed so far, at such high spectral resolution. We aim to select the optimal physical conditions in the PNe to be used in ionic abundance calculations that will be presented in a forthcoming paper. We acquired spectra at LCO with the 6.5-m telescope and the MIKE spectrograph, covering a wavelength range from 3350 to 9400 A. The spectra were exposed deep enough to detect, with signal-to-noise ratio higher than three, the weak ORLs of OII, CII, and other species. We detect and identify about 2980 emission lines, which, to date, is the most complete set of spectrophotometric data published for this type of objects. From our deep data, numerous diagnostic line ratios for Te and Ne are determined from CELs, ORLs, and continuum measurements. Densities are closely described by the average of all determined values for objects with ne<10^4 cm-3, and by ne([Cl III]) for the densest objects. For Te, we adopt a three-zone ionization scheme, where the low ionization zone is characterized by Te([NII]), the medium ionization zone by Te([OIII]) and the highest ionization one by Te([ArIV]) when available. We compute Te from the HI Paschen discontinuity and from HeI lines. We do not find evidences of the presence of low-temperature, high-density clumps in our [WR]PNe from the analysis of faint O II and N II plasma diagnostics. The behaviour of Te([OIII])/Te([NII]) being very low for high ionization degrees can be reproduced by a set of combined matter-bounded and radiation-bounded models, however, for the lowest temperature ratios, too high metallicity seem required.

We revisit the question of the ionization of the diffuse medium in late type galaxies, by studying NGC 891. The most important challenge for the models considered so far was the observed increase of [OIII]/Hb, [OII]/Hb, and [NII]/Ha with increasing distance to the galactic plane. We propose a scenario based on the expected population of massive OB stars and hot low-mass evolved stars (HOLMES) in this galaxy to explain this observational fact. In the framework of this scenario we construct a finely meshed grid of photoionization models. For each value of the galactic latitude z we look for the models which simultaneously fit the observed values of the [OIII]/Hb, [OII]/Hb, and [NII]/Ha ratios. For each value of z we find a range of solutions which depends on the value of the oxygen abundance. The models which fit the observations indicate a systematic decrease of the electron density with increasing z. They become dominated by the HOLMES with increasing z only when restricting to solar oxygen abundance models, which argues that the metallicity above the galactic plane should be close to solar. They also indicate that N/O increases with increasing z.

The extraction of chemical abundances of ionised nebulae from a limited spectral range is usually hampered by the lack of emission lines corresponding to certain ionic stages. So far, the missing emission lines have been accounted for by the ionisation correction factors (ICFs), constructed under simplistic assumptions like spherical geometry by using 1-D photoionisation modelling. In this contribution we discuss the results (Goncalves et al. 2011, in prep.) of our ongoing project to find a new set of ICFs to determine total abundances of N, O, Ne, Ar, and S, with optical spectra, in the case of non-spherical PNe. These results are based on a grid of 3-D photoionisation modelling of round, elliptical and bipolar shaped PNe, spanning the typical PN luminosities, effective temperatures and densities. We show that the additional corrections --to the widely used Kingsburgh and Barlow (1994) ICFs-- are always higher for bipolars than for ellipticals. Moreover, these additional corrections are, for bipolars, up to: 17% for oxygen, 33% for nitrogen, 40% for neon, 28% for argon and 50% for sulphur. Finally, on top of the fact that corrections change greatly with shape, they vary also greatly with the central star temperature, while the luminosity is a less important parameter.

A preliminary VLT-UVES spectrum of NGC 6302 (Casassus et al. 2002, MN), which hosts one of the hottest PN nuclei known (Teff ~ 220000 K; Wright et al. 2011, MN), has been recently analysed by means of X-SSN, a spectrum synthesis code for nebulae (Morisset and P\'equignot). Permitted recombination lines from highly-ionized species are detected/identified for the first time in a PN, and some of them probably for the first time in (Astro)Physics. The need for a homogeneous, high signal-to-noise UVES spectrum for NGC 6302 is advocated.

We predict intensities of lines of CII, NI, NII, OI and OII and compare them with a deep spectroscopic survey of IC 418 to test the effect of excitation of nebular emission lines by continuum fluorescence of starlight. Our calculations use a nebular model and a synthetic spectrum of its central star to take into account excitation of the lines by continuum fluorescence and recombination. The NII spectrum is mostly produced by fluorescence due to the low excitation conditions of the nebula, but many CII and OII lines have more excitation by fluorescence than recombination. In the neutral envelope, the NI permitted lines are excited by fluorescence, and almost all the OI lines are excited by recombination. Electron excitation produces the forbidden optical lines of OI, but continuum fluorescence excites most of the NI forbidden line intensities. Lines excited by fluorescence of light below the Lyman limit thus suggest a new diagnostic to explore the photodissociation region of a nebula.

We present preliminary results obtained from the analysis of very deep echelle spectra of a dozen planetary nebulae with [WC] or weak emission lines (wels) central stars. The computed abundance discrepancy factors (ADFs) are moderate, with values lower than 4. In principle, no evidence of the H-poor metal enriched inclusions proposed by Liu et al. (2000) have been found. However, a detailed analysis of the data is in progress.

We have selected the Galactic HII region M43, a close-by apparently spherical nebula ionized by a single star (HD37061) to investigate several topics of recent interest in the field of HII regions and massive stars. In a series of two papers we perform a combined, comprehensive study of the nebula and its ionizing star by using as many observational constraints as possible. We collected for this study a set of high-quality observations, including the optical spectrum of HD37061, along with nebular optical imaging and long-slit spatially resolved spectroscopy. The first part of our study comprises a quantitative spectroscopic analysis of the ionizing star, and the empirical analysis of the nebular images and spectroscopy. We determine the stellar parameters of HD37061 and the total number of ionizing photons emitted by the star. We find observational evidence of the presence of scattered light from the Huygens region (brightest part of the Orion nebula) in the M43 region. We show the importance of an adequate correction of this scattered light from the imagery and spectroscopic observations of M43 for a proper determination of the total nebular H_alpha luminosity, the nebular physical conditions and chemical abundances. We perform a detailed nebular empirical analysis of 9 apertures extracted from a long-slit located to the west of HD37061, obtaining the spatial distribution of the physical conditions and ionic abundances. For three of the analyzed elements (O, S, and N) we could determine total abundances directly from observable ions (no ionization correction factors were needed). The comparison of these abundances with those derived from the spectrum of the Orion nebula indicates the importance of the atomic data and, specially in the case of M\,42, the considered ionization correction factors.

We revisit the question of the ionization of the diffuse medium in late type galaxies, by studying NGC 891, the prototype of edge-on spiral galaxies. The most important challenge for the models considered so far was the observed increase of [OIII]/Hbeta, [OII]/Hbeta, and [NII]/Halpha with increasing distance to the galactic plane. We propose a scenario based on the expected population of massive OB stars and hot low-mass evolved stars (HOLMES) in this galaxy to explain this observational fact. In the framework of this scenario we construct a finely meshed grid of photoionization models. For each value of the galactic latitude z we look for the models which simultaneously fit the observed values of the [OIII]/Hbeta, [OII]/Hbeta, and [NII]/Halpha ratios. For each value of z we find a range of solutions which depends on the value of the oxygen abundance. The models which fit the observations indicate a systematic decrease of the electron density with increasing z. They become dominated by the HOLMES with increasing z only when restricting to solar oxygen abundance models, which argues that the metallicity above the galactic plane should be close to solar. They also indicate that N/O increases with increasing z.

We present a flexible interactive 3D morpho-kinematical modeling application for astrophysics. Compared to other systems, our application reduces the restrictions on the physical assumptions, data type and amount that is required for a reconstruction of an object's morphology. It is one of the first publicly available tools to apply interactive graphics to astrophysical modeling. The tool allows astrophysicists to provide a-priori knowledge about the object by interactively defining 3D structural elements. By direct comparison of model prediction with observational data, model parameters can then be automatically optimized to fit the observation. The tool has already been successfully used in a number of astrophysical research projects.

We present a detailed investigation of SBS1150+599A, a close binary star hosted by the planetary nebula PN G135.9+55.9 (TS01, Stasinska et al, 2009). The nebula, located in the Galactic halo, is the most oxygen-poor one known to date and is the only one known to harbor a double degenerate core. We present XMM-Newton observations of this object, which allowed the detection of the previously invisible component of the binary core, whose existence was inferred so far only from radial velocity and photometric variations. The parameters of the binary system were deduced from a wealth of information via three independent routes using the spectral energy distribution (from the infrared to X-rays), the light and radial velocity curves, and a detailed model atmosphere fitting of the stellar absorption features of the optical/UV component. We find that the cool component must have a mass of 0.54+/-0.2 Msun, an average effective temperature, Teff, of 58000+/-3000 K, a mean radius of 0.43+/-0.3 Rsun, a gravity log g=5.0+/-0.3, and that it nearly fills its Roche lobe. Its surface elemental abundances are found to be: 12 + log He/H = 10.95+/-0.04 dex, 12 + log C/H = 7.20+/-0.3 dex, 12 + log N/H < 6.92 and 12 + log O/H < 6.80, in overall agreement with the chemical composition of the planetary nebula. The hot component has Teff = 160-180 kK, a luminosity of about ~10e4 Lsun and a radius slightly larger than that of a white dwarf. It is probably bloated and heated as a result of intense accretion and nuclear burning on its surface in the past. The total mass of the binary system is very close to Chandrasekhar limit. This makes TS01 one of the best type Ia supernova progenitor candidates. We propose two possible scenarios for the evolution of the system up to its present stage.

The planetary nebula TS 01 (also called PN G 135.9+55.9 or SBS 1150+599A), with its record-holding low oxygen abundance and its double degenerate close binary core (period 3.9 h), is an exceptional object located in the Galactic halo. We have secured observational data in a complete wavelength range in order to pin down the abundances of half a dozen elements in the nebula. The abundances are obtained via detailed photoionization modelling taking into account all the observational constraints (including geometry and aperture effects) using the pseudo-3D photoionization code Cloudy_3D. The spectral energy distribution of the ionizing radiation is taken from appropriate model atmospheres. Both stellar components contribute to the ionization: the ``cool'' one provides the bulk of hydrogen ionization, and the ``hot'' one is responsible for the presence of the most highly charged ions, which explains why previous attempts to model the nebula experienced difficulties. The nebular abundances of C, N, O, and Ne are found to be respectively, 1/3.5, 1/4.2, 1/70, and 1/11 of the Solar value, with uncertainties of a factor 2. Thus the extreme O deficiency of this object is confirmed. The abundances of S and Ar are less than 1/30 of Solar. Standard models of stellar evolution and nucleosynthesis cannot explain the abundance pattern observed in the nebula. To obtain an extreme oxygen deficiency in a star whose progenitor has an initial mass of about 1 msun requires an additional mixing process, which can be induced by stellar rotation and/or by the presence of the close companion. We have computed a stellar model with initial mass of 1 msun, appropriate metallicity, and initial rotation of 100 kms, and find that rotation greatly improves the agreement between the predicted and observed abundances.

Studies of the Diffuse Ionized Gas (DIG) have progressed without providing so far any strict criterion to distinguish DIGs from HII regions. In this work, we compile the emission line measurements of 29 galaxies that are available in the scientific literature, thereby setting up the first DIG database (DIGEDA). Making use of this database, we proceed to analyze the global properties of the DIG using the [NII]/Ha, [OI]/Ha, [OIII]/Hb and [SII]/Ha lines ratios, including the H alpha emission measure. This analysis leads us to conclude that the [NII]/Ha ratio provides an objective criterion for distinguishing whether an emission region is a DIG or an HII region, while the EM(Ha) is a useful quantity only when the galaxies are considered individually. Finally, we find that the emission regions of Irr galaxies classified as DIG in the literature appear in fact to be much more similar to HII regions than to the DIGs of spiral galaxies.

We present a coherent stellar and nebular model reproducing the observations of the Planetary Nebula IC418. We want to test whether a stellar model obtained by fitting the stellar observations is able to satisfactory ionize the nebula and reproduce the nebular observations, which is by no mean evident. This allows us to determine all the physical parameters of both the star and the nebula, including the abundances and the distance. We used all the observational material available (FUSE, IUE, STIS and optical spectra) to constrain the stellar atmosphere model performed using the CMFGEN code. The photoionization model is done with Cloudy_3D, and is based on CTIO, Lick, SPM, IUE and ISO spectra as well as HST images. More than 140 nebular emission lines are compared to the observed intensities. We reproduce all the observations for the star and the nebula. The 3D morphology of the gas distribution is determined. The effective temperature of the star is 36.7kK. Its luminosity is 7700 solar luminosity. We describe an original method to determine the distance of the nebula using evolutionary tracks. No clumping factor is need to reproduce the age-luminosity relation. The distance of 1.25 kpc is found in very good agreement with recent determination using parallax method. The chemical composition of both the star and the nebula are determined. Both are Carbon-rich. The nebula presents evidence of depletion of elements Mg, Si, S, Cl (0.5 dex lower than solar) and Fe (2.9 dex lower than solar). This is the first self-consistent stellar and nebular model for a Planetary Nebula that reproduces all the available observations ranging from IR to UV, showing that the combined approach for the modeling process leads to more restrictive constraints and, in principle, more trustworthy results.

The 3MdB (Mexican Million Models database) is a large database of photoionization models for H II regions. The number of free parameters for the models is close to 15, including the description of the ionizing Spectral Energy Distribution (effective temperature, luminosity, surface gravity, for different type of stellar atmosphere models) and the description of the ionized gas (distance to the ionizing source, density, abundances of the most common elements, dust). The outputs of the models are more than 70 emission line intensities, the ionic fractions and temperatures. All the parameters and outputs are included in the MySQL database, giving the possibility to the user to search into the database for example for all the models that reproduce a given set of observations.

In this paper, we present radio continuum observations of the planetary nebula IC 418 obtained at two epochs separated by more than 20 years. These data allow us to show that the angular expansion rate of the ionization front in IC 418 is 5.8 $\pm$ 1.5 mas yr$^{-1}$. If the expansion velocity of the ionization front is equal to the expansion velocity of the gas along the line of sight as measured by optical spectroscopy, then the distance to IC 418 must be 1.1 $\pm$ 0.3 kpc. Recent theoretical predictions appropriate for the case of IC 418, however, suggest that the ionization front may be expanding about 20% faster than the material. Under this assumption, the distance to IC 418 would increase to 1.3 $\pm$ 0.4 kpc.

It is generally believed that O stars, confined near the galactic midplane, are somehow able to photoionize a significant fraction of what is termed the "diffuse ionized gas" (DIG) of spiral galaxies, which can extend up to 1-2 kpc above the galactic midplane. The heating of the DIG remains poorly understood, however, as simple photoionization models do not reproduce the observed line ratio correlations well or the DIG temperature. We present turbulent mixing layer models in which warm photoionized condensations are immersed in a hot supersonic wind. Turbulent dissipation and mixing generate an intermediate region where the gas is accelerated, heated and mixed. The emission spectrum of such layers are compared with observations of Rand (ApJ 462, 712) of the DIG in the edge-on spiral NGC2363. We generate two sequence of models that fit the line ratio correlations between [SII]/H-alpha, [OI]/H-alpha, [NII]/[SII] and [OIII]/H-beta reasonably well. In one sequence of models the hot wind velocity increases while in the other the ionization parameter and layer opacity increases. Despite the success of the mixing layer models, the overall efficiency in reprocessing the stellar UV is much too low, much less than 1%, which compels us to reject the TML model in its present form.

We analyzed the regions of star formation in the spiral galaxies NGC 3184 and NGC 3938 from archive images at multiple wavelengths (NUV from GALEX, H-alpha from JKT and KPNO, 8 and 24 um from Spitzer, and CO from BIMA). We used the Clump Find Algorithm to extract the properties of the star forming tracers shown as emission regions at each wavelength. We obtained a power-law relation between the luminosity and the emission region volume that scales as expected, $L proportional to V, for the H-alpha and NUV emission; but much faster for the dust (8 and 24 um) and molecular gas emitting regions in CO. This is interpreted as a change on the emissivity with the size of the cloud, either by an augmentation of the overall density or due to the presence of high density clumps, with high local emissivity coefficients. Although the clumpy nature of molecular gas may not be surprising, the clumpy nature of mid-infrared emission regions, that could be explained by newly formed high to intermediate mass stars embedded into the dust and heating it, is clearly revealed in both galaxies.

The kinematics, shaping, density distribution, expansion distance, and ionized mass of the nebula Hen 2-104, and the nature of its symbiotic Mira are investigated. A combination of multi-epoch HST images and VLT integral field high-resolution spectroscopy is used to study the nebular dynamics both along the line of sight and in the plane of the sky. These observations allow us to construct a 3-D spatio-kinematical model of the nebula, which together with the measurement of its apparent expansion in the plane of the sky over a period of 4 years, provides the expansion parallax for the nebula. The integral field data featuring the [S{\sc ii}] $\lambda\lambda$671.7,673.1 emission line doublet provide us with a density map of the inner lobes of the nebula, which together with the distance estimation allow us to estimate its ionized mass. We find densities ranging from n$_e$=500 to 1000 cm$^{-3}$ in the inner lobes and from 300 to 500cm$^{-3}$ in the outer lobes. We determine an expansion parallax distance of 3.3$\pm$0.9 kpc to Hen 2-104, implying an unexpectedly large ionized mass for the nebula of the order of one tenth of a solar mass.

We have constructed a grid of photoionization models of spherical, elliptical and bipolar planetary nebulae. Assuming different velocity fields, we have computed line profiles corresponding to different orientations, slit sizes and positions. The atlas is meant both for didactic purposes and for the interpretation of data on real nebulae. As an application, we have shown that line profiles are often degenerate, and that recovering the geometry and velocity field from observations requires lines from ions with different masses and different ionization potentials. We have also shown that the empirical way to measure mass-weighted expansion velocities from observed line widths is reasonably accurate if considering the HWHM. For distant nebulae, entirely covered by the slit, the unknown geometry and orientation do not alter the measured velocities statistically. The atlas is freely accessible from internet. The Cloudy_3D suite and the associated VISNEB tool are available on request.

We present the first results of a project aimed at the combined study of massive stars and their surrounding nebulae by means of a detailed study of Galactic HII regions ionized by only one massive star. With this, we intend to check the validity of the new generation of massive star model atmosphere codes in terms of ionizing flux distribution. We take into account the effect of the nebular density distribution in our analyses. Various types of stellar and nebular observations have been collected for this purpose.

We developed a new quick pseudo-3D photoionization code based on Cloudy (G. Ferland) and IDL (RSI) tools. The code is running the 1D photoionization code Cloudy various times, changing at each run the input parameters (e.g. inner radius, density law) according to an angular law describing the morphology of the object. Then a cube is generated by interpolating the outputs of Cloudy. In each cell of the cube, the physical conditions (electron temperature and density, ionic fractions) and the emissivities of lines are determined. Associated tools (VISNEB and VELNEB_3D) are used to rotate the nebula and to compute surface brightness maps and emission line profiles, given a velocity law and taking into account the effect of the thermal broadening and eventually the turbulence. Integrated emission line profiles are computed, given aperture shapes and positions (seeing and instrumental width effects are included). The main advantage of this tool is the short time needed to compute a model (a few tens minutes).

Using the Cloudy_3D code (Morisset, this conference), we run a set of models to build a catalog of emission line profiles of PNe. The goal is to cover the main morphologies (spherical, ellipsoidal, bipolar, thin or thick shells) and to compute profiles changing the orientation of the nebula, the expansion velocity law (including or not turbulence), the position/size of the aperture. The results is a huge set of profiles and PV-diagrams that can be compared to observations as a help to derive morpho-kinematical properties of real objects.

We present a tool, VELNEB_3D, which can be applied to the results of 3D photoionization codes to generate emission line profiles, position-velocity maps and 3D maps in any emission line by assuming an arbitrary velocity field. We give a few examples, based on our pseudo-3D photoionization code NEBU_3D (Morisset, Stasinska and Pena, 2005) which show the potentiality and usefulness of our tool. One example shows how complex line profiles can be obtained even with a simple expansion law if the nebula is bipolar and the slit slightly off-center. Another example shows different ways to produce line profiles that could be attributed to a turbulent velocity field while there is no turbulence in the model. A third example shows how, in certain circumstances, it is possible to discriminate between two very different geometrical structures -- here a face-on blister and its ``spherical impostor'' -- when using appropriate high resolution spectra. Finally, we show how our tool is able to generate 3D maps, similar to the ones that can be obtained by observing extended nebulae with integral field units.

Expansion velocities (Vexp) of different ions and line widths at the base of the lines are measured and analyzed for 24 PNe with [WC]-type nuclei (WRPNe), 9 PNe ionized by WELS (WLPNe) and 14 ordinary PNe. A comparative study of the kinematical behavior of the sample clearly demonstrates that WRPNe have in average 40-45% larger Vexp, and possibly more turbulence than WLPNe and ordinary PNe. WLPNe have velocity fields very much alike the ones of ordinary PNe, rather than the ones of WRPNe. All the samples (WRPNe, WLPNe and ordinary PNe) show expansion velocities increasing with age indicators, for example <Vexp> is larger for low-density nebulae and also it is larger for nebulae around high-temperature stars. This age effect is much stronger for evolved WRPNe, suggesting that the [WC] winds have been accelerating the nebulae for a long time, while for non-WRPNe the acceleration seems to stop at some point when the star reaches a temperature of about 90,000-100,000 K. Non-WR nebulae reach a maximum Vexp < 30 km/s while evolved WRPNe reach maximum Vexp of about 40 km/s. For all kind of objects (WRPNe and non-WRPNe) it is found that in average Vexp(N+) is slightly larger than Vexp(O++), indicating that the nebulae present acceleration of the external shells.

The ultraviolet energy distribution of quasars shows a sharp steepening of the continuum shortward of 1000 A (rest-frame). We describe how we came to consider the possibility that this continuum break might be the result of absorption by carbon crystallite dust grains.

We describe a pseudo-3D photoionization code, NEBU_3D and its associated visualization tool, VIS_NEB3D, which are able to easily and rapidly treat a wide variety of nebular geometries, by combining models obtained with a 1D photoionization code. We also present a tool, VELNEB_3D, which can be applied to the results of 1D or 3D photoionization codes to generate emission line profiles, position-velocity maps and 3D maps in any emission line by assuming an arbitrary velocity field. As examples of the capabilities of these new tools, we consider three very different theoretical cases. The first one is a blister HII region, for which we have also constructed a spherical model (the spherical impostor) which has exactly the same Hbeta surface brightness distribution as the blister model and the same ionizing star. The second example shows how complex line profiles can be obtained even with a simple expansion law if the nebula is bipolar and the slit slightly off-center. The third example shows different ways to produce line profiles that could be attributed to a turbulent velocity field while there is no turbulence in the model.

The predicted emission spectrum of NII is compared with observations of permitted lines in the Orion nebula. Conventional nebular models show that the intensities of the more intense lines can be explained by fluorescence of starlight absorption with a N abundance that is consistent with forbidden lines. Lines excited mostly by recombination on the other hand predict high N abundances. The effects of stellar and nebular parameters and of the atomic data on the predicted intensities are examined.

We present high spatial resolution spectroscopic observations of the proplyd 167-317 (LV2) near the Trapezium cluster in the Orion nebula, obtained during the System Verification run of the Gemini Multi Object Spectrograph (GMOS) Integral Field Unit (IFU) at the Gemini South Observatory. We have detected 38 forbidden and permitted emission lines associated with the proplyd and its redshifted jet. We have been able to detect three velocity components in the profiles of some of these lines: a peak with a 28-33 km/s systemic velocity that is associated with the photoevaporated proplyd flow, a highly redshifted component associated with a previously reported jet (which has receding velocities of about 80-120 km/s with respect to the systemic velocity and is spatially distributed to the southeast of the proplyd) and a less obvious, approaching structure, which may possibly be associated with a faint counter-jet with systemic velocity of (-75 +/- 15) km/s. We find evidences that the redshifted jet has a variable velocity, with slow fluctuations as a function of the distance from the proplyd. We present several background subtracted, spatially distributed emission line maps and we use this information to obtain the dynamical characteristics over the observed field. Using a simple model and with the extinction corrected Halpha fluxes, we estimate the mass loss rate for both the proplyd photoevaporated flow and the redshifted microjet, obtaining (6.2 +/- 0.6) x 10^{-7} M_sun/year and (2.0 +/- 0.7) x 10^{-8} M_sun/year, respectively.

We explore the possibility that the steepening observed shortward of 1000A in the energy distribution of quasars may result from absorption by dust, being either intrinsic to the quasar environment or intergalactic. We find that a dust extinction curve consisting of nanodiamonds, composed of terrestrial cubic diamonds or with surface impurities as found in carbonaceous chondrite meteorites, such as Allende, is successful in reproducing the sharp break observed. The intergalactic dust model is partially successful in explaining the shape of the composite energy distribution, but must be discarded in the end, as the amount of crystalline dust required is unreasonable and would imply an improbable fine tuning among the dust formation processes. The alternative intrinsic dust model requires a mixture of both cubic diamonds and Allende nanodiamonds and provide a better fit of the UV break. The gas column densities implied are of the order 10^{20} cm^{-2} assuming solar metallicity for carbon and full depletion of carbon into dust. The absorption only occurs in the ultraviolet and is totally negligible in the visible. The minimum dust mass required is of the order ~ 0.003 r_{pc}^{2}M_o, where r_{pc} is the distance in parsec between the dust screen and the continuum source. The intrinsic dust model reproduces the flux {\it rise} observed around 660A in key quasar spectra quite well. We present indirect evidence of a shallow continuum break near 670A (18.5 eV), which would be intrinsic to the quasar continuum.

We describe a pseudo-3D photoionization code, NEBU_3D and its associated visualization tool, VIS_NEB3D, which are able to easily and rapidly treat a wide variety of nebular geometries, by combining models obtained with a 1D photoionization code. The only requirement for the code to work is that the ionization source is uniqu e and not extended. It is applicable as long as the diffuse ionizing radiation f ield is not dominant and strongly inhomogeneous. As examples of the capabilities of these new tools, we consider two very differ ent theoretical cases. One is that of a high excitation planetary nebula that ha s an ellipsoidal shape with two polar density knots. The other one is that of a blister HII region, for which we have also constructed a spherical model (the sp herical impostor) which has exactly the same Hbeta surface brightness distrib ution as the blister model and the same ionizing star. These two examples warn against preconceived ideas when interpreting spectroscop ic and imaging data of HII regions and planetary nebulae. The tools NEBU_3D and VIS_NEB3D, which will be made publicly available in the future, should facilitat e the performance of numerical experiments, to yield a better understanding of t he physics of aspherical ionized nebulae.

Extensive photoionization model grids are computed for single star HII regions using stellar atmosphere models from the WM-basic code. Mid-IR emission line intensities are predicted and diagnostic diagrams of [NeIII]/[NeII] and [SIV]/[SIII] excitation ratio are build, taking into account the metallicities of both the star and the HII region. The diagrams are used in conjunction with galactic HII region observations obtained with the ISO Observatory to determine the effective temperature Teff of the exciting O stars and the mean ionization parameter U. Teff and U are found to increase and decrease, respectively, with the metallicity of the HII region represented by the [Ne/Ne_sol] ratio. No evidence is found for gradients of Teff or U with galactocentric distance Rgal. The observed excitation sequence with Rgal is mainly due to the effect of the metallicity gradient on the spectral ionizing shape, upon which the effect of an increase in Teff with Z is superimposed. We show that not taking properly into account the effect of metallicity on the ionizing shape of the stellar atmosphere would lead to an apparent decrease of Teff with Z and an increase of Teff with Rgal.

Extensive photoionization model grids for single star HII regions using a variety of recent state-of-the-art stellar atmosphere models have been computed with the main aim of constraining/testing their predicted ionizing spectra against recent ISO mid-IR observations of Galactic HII regions, which probe the ionizing spectra between ~ 24 and 41 eV thanks to Ne, Ar, and S fine structure lines. Particular care has been paid to examining in detail the dependences of the nebular properties on the numerous nebular parameters (mean ionization parameter U, abundances, dust etc.) which are generally unconstrained for the objects considered here. Finally we have examined which parameters are chiefly responsible for the observed mid-IR excitation sequences. The galactic gradient of metallicity changing the shape of the stellar emission is found to be one of the drivers for the excitation sequence of Galactic HII regions, the actual contribution of this effect being finally atmosphere model dependent. We find that the dispersion of Teff between different HII regions due to statistical sampling of the IMF plus additional scatter in the ionization parameter are probably the dominant driver for the observed excitation scatter.

Most of the enveloppes of Planetary Nebulae (and other objects like novae) are far from beeing homogeneous: clumps, knots and tails are often observed. We present here the first attempt to build a 3D-photoionization model of a knot and the corresponding tail, ionized by diffused radiation issuing from surrounding material.

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.

Two independent groups (Giveon et al. 2002; Martin-Hernandez et al. 2002) have recently investigated the Galactic metallicity gradient as probed by ISO observations of mid-infrared emission lines from HII regions. We show that the different gradients inferred by the two groups are due to differing source selection and differing extinction corrections. We show that both data sets in fact provide consistent results if identical assumptions are made in the analysis. We present a consistent set of gradients in which we account for extinction and variation in electron temperature across the disk.

We present a detailed photoionization model of G29.96-0.02 (hereafter G29.96), one of the brightest Galactic Ultra Compact HII (UCHII) regions in the Galaxy. This source has been observed extensively at radio and infrared wavelengths. The most recent data include a complete ISO (SWS and LWS) spectrum, which displays a remarkable richness in atomic fine-structure lines. The number of observables is twice as much as the number available in previous studies. In addition, most atomic species are now observed in two ionization stages. The radio and infrared data on G29.96 are best reproduced using a nebular model with two density components: a diffuse (n_e~680cm-3) extended (~1 pc) component surrounding a compact (~0.1 pc) dense (n_e~57000cm-3) core. The properties of the ionizing star were derived using state-of-the-art stellar atmosphere models. CoStar models yield an effective temperature of \~30^{+2}_{-1} kK whereas more recent non-LTE line blanketed atmospheres with stellar winds indicate somewhat higher values, Teff~32--38 kK. This range in Teff is compatible with all observational constraints, including near-infrared photometry and bolometric luminosity. The range 33-36 kK is also compatible with the spectral type O5-O8 determined by Watson and Hanson (97) when recent downward revisions of the effective temperature scale of O stars are taken into account. The age of the ionizing star of G29.96 is found to be a few 10^6 yr, much older than the expected lifetime of UCHII regions. Accurate gas phase abundances are derived with the most robust results being Ne/S=7.5 and N/O=0.43 (1.3 and 3.5 times the solar values, respectively).

Based on the ISO spectral catalogue of compact HII regions by Peeters et al. (2001), we present a first analysis of the hydrogen recombination and atomic fine-structure lines originated in the ionized gas. The sample consists of 34 HII regions located at galactocentric distances between Rgal = 0 and 15 kpc. The SWS HI recombination lines between 2 and 8 mum are used to estimate the extinction law at these wavelengths for 14 HII regions. An extinction in the K band between 0 and $\sim$ 3 mag. has been derived. The fine-structure lines of N, O, Ne, S and Ar are detected in most of the sources. Most of these elements are observed in two different ionization stages probing a range in ionization potential up to 41 eV. The ISO data, by itself or combined with radio data taken from the literature, is used to derive the elemental abundances relative to hydrogen. The present data thus allow us to describe for each source its elemental abundance, its state of ionization and to constrain the properties of the ionizing star(s).

We use a 3D photoionization modeling tool to study the morpho-kinematic properties of the Planetary Nebula (PN) NGC 3132. We show that it is possible to reproduce the low resolution observations (spectra and images) with an ellipsoidal shell. However, high resolution observations, as those showing a density variation along the nebula and the [O III]5007 velocity profiles, definitively rule out this description. We show that a bipolar Diabolo shape with a 40 degres rotation of the symmetry axis relative to the line of sight successfully reproduces the observed images, as well as the high resolution observations.

We present a new modeling tool for planetary nebulae, based on 3D photoionization calculations. Our goal is to show that all the information provided by observations, regarding kinematics and morphology, have to be consistently accounted for, in order to get a real insight of the object. Only 3D simulations offer this possibility. From models for two theoretical PNe, we show that the enhancement in the equatorial zone observed in several PNe is not necessarily due to a density gradient, as usually interpreted. It is also shown that asymmetric velocity profiles often observed (e.g., Gesicki et al. 1998) can be easily reproduced. Observations providing a better insight on the morphology of the PN are discussed.

We present a new modeling tool for planetary nebulae (PNe), based on 3D photoionization calculations. From models for two theoretical PNe, we show that the enhancement in the equatorial zone observed in several PNe is not necessarily due to a density gradient, as usually interpreted. We study the morpho-kinematic properties of the PN NGC 3132 and show that a bipolar Diabolo shape successfully reproduces the observed images, as well as the high resolution observations.