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

Spectral variation in the X-ray pulsar GX 1+4 during a low-flux episode

The X-ray pulsar GX 1+4 was observed with the RXTE satellite for a total of 51ks between 1996 July 19 - 21. During this period the flux decreased smoothly from an initial mean level of ~ 6 X 10^36 erg/s to a minimum of ~ 4 X 10^35 erg/s (2-60 keV, assuming a source distance of 10 kpc) before partially recovering towards the initial level at the end of the observation. BATSE pulse timing measurements indicate that a torque reversal took place approximately 10 d after this observation. Both the mean pulse profile and the photon spectrum varied significantly. The observed variation in the source may provide important clues as to the mechanism of torque reversals. The single best-fitting spectral model was based on a component originating from thermal photons with kT ~ 1 keV Comptonised by a plasma of temperature kT \~ 7 keV. Both the flux modulation with phase during the brightest interval and the evolution of the mean spectra over the course of the observation are consistent with variations in this model component; with, in addition, a doubling of the column density nH contributing to the mean spectral change. A strong flare of duration 50 s was observed during the interval of minimum flux, with the peak flux ~ 20 times the mean level. Although beaming effects are likely to mask the true variation in Mdot thought to give rise to the flare, the timing of a modest increase in flux prior to the flare is consistent with dual episodes of accretion resulting from successive orbits of a locally dense patch of matter in the accretion disc.Comment: 8 pages, 3 figures, submitted to MNRA

Pulse Profiles, Accretion Column Dips and a Flare in GX 1+4 During a Faint State

The Rossi X-ray Timing Explorer (RXTE) spacecraft observed the X-ray pulsar GX 1+4 for a period of 34 hours on July 19/20 1996. The source faded from an intensity of ~20 mCrab to a minimum of <~0.7 mCrab and then partially recovered towards the end of the observation. This extended minimum lasted ~40,000 seconds. Phase folded light curves at a barycentric rotation period of 124.36568 +/- 0.00020 seconds show that near the center of the extended minimum the source stopped pulsing in the traditional sense but retained a weak dip feature at the rotation period. Away from the extended minimum the dips are progressively narrower at higher energies and may be interpreted as obscurations or eclipses of the hot spot by the accretion column. The pulse profile changed from leading-edge bright before the extended minimum to trailing-edge bright after it. Data from the Burst and Transient Source Experiment (BATSE) show that a torque reversal occurred <10 days after our observation. Our data indicate that the observed rotation departs from a constant period with a Pdot/P value of ~-1.5% per year at a 4.5 sigma significance. We infer that we may have serendipitously obtained data, with high sensitivity and temporal resolution about the time of an accretion disk spin reversal. We also observed a rapid flare which had some precursor activity, close to the center of the extended minimum.Comment: 19 pages, 6 figures, accepted for publication in Astrophysical Journal (tentatively scheduled for vol. 529 #1, 20 Jan 2000

Overcoming the boundary layer turbulence at Dome C: ground-layer adaptive optics versus tower

The unique atmospheric conditions present at sites such as Dome C on the Antarctic plateau are very favorable for high spatial resolution astronomy. At Dome C, the majority of the optical turbulence is confined to a 30 to 40 m thick stable boundary layer that results from the strong temperature inversion created by the heat exchange between the air and the ice-covered ground. To fully realize the potential of the exceptionally calm free atmosphere, this boundary layer must be overcome. In this article we compare the performance of two methods proposed to beat the boundary layer: mounting a telescope on a tower that physically puts it above the turbulent layer, and installing a telescope at ground level with a ground-layer adaptive optics system. A case is also made to combine these two methods to further improve the image quality

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