Three-Dimensional Ionisation, Dust RT and Chemical Modelling of Planetary Nebulae

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.Comment: 3 pages, to appear in Proc. IAU Symp. 234, Planetary Nebulae in Our Galaxy and Beyond (3-7 Apr 2006), eds. M.J. Barlow & R.H. Mendez (Cambridge Univ. Press

Mocassin: A fully three-dimensional Monte Carlo photoionization code

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.Comment: 17 pages, 4 figures, 1 appendix Changes: appendix adde

Theoretical calculations of the HI, HeI and HeII free-bound continuou emission spectra

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

A change in the geodynamics of continental growth 3 billion years ago

Continental Growth Spurts The appearance and persistence of continents through geologic time has influenced most processes on Earth, from the evolution of new species to the climate. The relative proportion of newly formed crust compared to reworked, or destroyed, older crust reveals which processes controlled continental growth. Based on the combined analyses of Hf-Pb and O isotopes in zircon minerals, Dhuime et al. (p. 1334 ) measured continuous but variable rates of new crustal production throughout Earth's history. Increased rates of crustal destruction starting around 3 billion years ago coincide with the onset of subduction-drive plate tectonics, slowing down the overall rate of crustal growth. </jats:p

The dusty MOCASSIN: fully self-consistent 3D photoionisation and dust radiative transfer models

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.Comment: 9 pages, 3 figure

The 21cm Signature of the First Stars

We predict the 21-cm signature of the first metal-free stars. The soft X-rays emitted by these stars penetrate the atomic medium around their host halos, generating Lyman alpha photons that couple the spin and kinetic temperatures. These creates a region we call the Lyman alpha sphere, visible in 21-cm against the CMB, which is much larger than the HII region produced by the same star. The spin and kinetic temperatures are strongly coupled before the X-rays can substantially heat the medium, implying that a strong 21-cm absorption signal from the adiabatically cooled gas in Hubble expansion around the star is expected when the medium has not been heated previously. A central region of emission from the gas heated by the soft X-rays is also present although with a weaker signal than the absorption. The Lyman alpha sphere is a universal signature that should be observed around any first star illuminating its vicinity for the first time. The 21-cm radial profile of the Lyman alpha sphere can be calculated as a function of the luminosity, spectrum and age of the star. For a star of a few hundred solar masses and zero metallicity (as expected for the first stars), the physical radius of the Lyman alpha sphere can reach tens of kiloparsecs. The first metal-free stars should be strongly clustered because of high cosmic biasing; this implies that the regions producing a 21-cm absorption signal may contain more than one star and will generally be irregular and not spherical, because of the complex distribution of the gas. We discuss the feasiblity of detecting these Lyman alpha spheres, which would be present at redshifts $z\sim 30$ in the Cold Dark Matter model. Their observation would represent a direct proof of the detection of a first star.Comment: replaced with ApJ accepted version. Many minor revisions and additional references, major results unchange

A New Look At Carbon Abundances In Planetary Nebulae. III. DDDM1, IC 3568, IC4593, NGC 6210, NGC 6720, NGC 6826, & NGC 7009

This paper is the third in a series reporting on a study of carbon abundances in a carefully chosen sample of planetary nebulae representing a large range in progenitor mass and metallicity. We make use of the IUE Final Archive database containing consistently-reduced spectra to measure line strengths of C III] 1909 along with numerous other UV lines for the planetary nebulae DDDM1, IC 3568, IC 4593, NGC 6210, NGC 6720, NGC 6826, & NGC 7009. We combine the IUE data with line strengths from optical spectra obtained specifically to match the IUE slit positions as closely as possible, to determine values for the abundance ratios He/H, O/H, C/O, N/O, and Ne/O. The ratio of C III] 1909/C II 4267 is found to be effective for merging UV and optical spectra when He II 1640/4686 is unavailable. Our abundance determination method includes a 5-level program whose results are fine-tuned by corrections derived from detailed photoionization models constrained by the same set of emission lines. All objects appear to have subsolar levels of O/H, and all but one show N/O levels above solar. In addition, the seven planetary nebulae span a broad range in C/O values. We infer that many of our objects are matter bounded, and thus the standard ionization correction factor for N/O may be inappropriate for these PNe. Finally, we estimate C/O using both collisionally-excited and recombination lines associated with C+2 and find the well established result that abundances from recombination lines usually exceed those from collisionally-excited lines by several times.Comment: 36 pages, 7 tables, 2 figures, latex. Tables and figures supplied as two separate postscript files. Accepted for publication in Ap

Momentum transfer for momentum transfer-free which-path experiments

We analyze the origin of interference disappearance in which-path double aperture experiments. We show that we can unambiguously define an observable momentum transfer between the quantum particle and the path detector and we prove in particular that the so called ``momentum transfer free'' experiments can be in fact logically interpreted in term of momentum transfer.Comment: to appear in Phys. Rev . A (2006). (7 pages, 2 figures