The solid accretion rate, necessary to grow gas giant planetary cores within the disk lifetime, has been a major constraint for theories of planet formation. We tested the solid accretion rate efficiency on planetary cores of different masses embedded in their birth disk, by means of 3D radiation-hydrodynamics, where we followed the evolution of a swarm of embedded solids of different sizes. We found that using a realistic equation of state and radiative cooling, the disk at 5 au is able to cool efficiently and reduce its aspect ratio. As a result, the pebble isolation mass is reached before the core grows to 10 Earth masses, stopping efficiently the pebble flux and creating a transition disk. Moreover, the reduced isolation mass halts the solid accretion before the core reaches the critical mass, leading to a barrier to giant planet formation, and it explains the large abundance of super-Earth planets in the observed population.

High resolution spectra of typical wind diagnostics ([OI] 6300 \r{A} and other forbidden emission lines) can often be decomposed into multiple components: high-velocity components with blueshifts up to several 100 km/s are usually attributed to fast jets, while narrow (NLVC) and broad (BLVC) low-velocity components are believed to trace slower disc winds. Under the assumption that the line-broadening is dominated by Keplerian rotation, several studies have found that the BLVCs should trace gas launched between 0.05 and 0.5 au and correlations between the properties of BLVCs and NLVCs have been interpreted as evidence for the emission tracing an extended MHD wind and not a photoevaporative wind. We calculated synthetic line profiles obtained from detailed photoionisation calculations of an X-ray photoevaporation model and a simple MHD wind model and analyzed the emission regions of different diagnostic lines and the resulting spectral profiles. The photoevaporation model reproduces most of the observed NLVCs but not the BLVCs or HVCs. The MHD model is able to reproduce all components but produces Keplerian double peaks at average inclinations that are rarely observed. The combination of MHD and photoevaporative winds could solve this problem. Our results suggest that the Gaussian decomposition does not allow for a clear distinction of flux from different wind regions and that the line broadening is often dominated by the velocity gradient in the outflow rather than by Keplerian rotation. We show that observed correlations between BLVC and NLVC do not necessarily imply a common origin in an extended MHD wind.

The dispersal phase of planet-forming disks via winds driven by irradiation from the central star and/or magnetic fields in the disk itself is likely to play an important role in the formation and evolution of planetary systems. Current theoretical models lack predictive power to adequately constrain observations. We present PRIZMO, a code for evolving thermochemistry in protoplanetary disks capable of being coupled with hydrodynamical and multi-frequency radiative transfer codes. We describe the main features of the code, including gas and surface chemistry, photochemistry, microphysics, and the main cooling and heating processes. The results of a suite of benchmarks, which include photon-dominated regions, slabs illuminated by radiation spectra that include X-ray, and well-established cooling functions evaluated at different temperatures show good agreement both in terms of chemical and thermal structures. The development of this code is an important step to perform quantitative spectroscopy of disk winds, and ultimately the calculation of line profiles, which is urgently needed to shed light on the nature of observed disk winds.

We calculated hydrogen recombination line luminosities (H-$\alpha$, Paschen-$\beta$ and Brackett-$\gamma$) from three dimensional thermo-hydrodynamical simulations of forming planets from 1 to 10 Jupiter-masses. We explored various opacities to estimate the line emissions with extinction in each cases. When realistic opacities are considered, only lines from planets $\ge$10 Jupiter-mass can be detected with current instrumentation, highlighting that from most planets one cannot expect detectable emission. This might explain the very low detection rate of H-$\alpha$ from forming planets from observations. While the line emission comes from both the forming planet and its circumplanetary disk, we found that only the disk component could be detected due to extinction. We examined the line variability as well, and found that it is higher for higher mass planets. Furthermore, we determine for the first time, the parametric relationship between the mass of the planet and the luminosity of the hydrogen recombination lines, as well as the equation between the accretion luminosity and hydrogen recombination line luminosities. These relationships for planets appear to be steeper than the analogous relation for young stars, suggesting that the accretion process around planets proceeds differently than around stars.

X-ray- and EUV- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may crucially impact disk evolution, affecting both gas and dust distributions. We investigate the dust entrainment in XEUV-driven photoevaporative winds and compare our results to existing MHD and EUV-only models. For an X-ray luminosity of $L_X = 2 \cdot 10^{30}\,\mathrm{erg/s}$ emitted by a $M_* = 0.7\,\mathrm{M}_\odot$ star, corresponding to a wind mass-loss rate of $\dot{M}_\mathrm{w} \simeq 2.6 \cdot 10^{-8} \,\mathrm{M_\odot/yr}$, we find dust entrainment for sizes $a_0 \lesssim 11\,\mu$m ($9\,\mu$m) from the inner $25\,$AU ($120\,$AU). This is an enhancement over dust entrainment in less vigorous EUV-driven winds with $\dot{M}_\mathrm{w} \simeq 10^{-10}\,\mathrm{M_\odot/yr}$. Our numerical model also shows deviations of dust grain trajectories from the gas streamlines even for $\mu$m-sized particles. In addition, we find a correlation between the size of the entrained grains and the maximum height they reach in the outflow.

We predict magnitudes for young planets embedded in transition discs, still affected by extinction due to material in the disc. We focus on Jupiter-size planets at a late stage of their formation, when the planet has carved a deep gap in the gas and dust distributions and the disc starts being transparent to the planet flux in the infrared (IR). Column densities are estimated by means of three-dimensional hydrodynamical models, performed for several planet masses. Expected magnitudes are obtained by using typical extinction properties of the disc material and evolutionary models of giant planets. For the simulated cases located at $5.2$ AU in a disc with local unperturbed surface density of $127$ $\mathrm{g} \cdot \mathrm{cm}^{-2}$, a $1$ $M_J$ planet is highly extincted in J-, H- and K-bands, with predicted absolute magnitudes $\ge 50$ mag. In L- and M-bands extinction decreases, with planet magnitudes between $25$ and $35$ mag. In the N-band, due to the silicate feature on the dust opacities, the expected magnitude increases to $40$ mag. For a $2$ $M_J$ planet, the magnitudes in J-, H- and K-bands are above $22$ mag, while for L-, M- and N-bands the planet magnitudes are between $15$ and $20$ mag. For the $5$ $M_J$ planet, extinction does not play a role in any IR band, due to its ability to open deep gaps. Contrast curves are derived for the transition discs in CQ Tau, PDS70, HL Tau, TW Hya and HD163296. Planet mass upper-limits are estimated for the known gaps in the last two systems.

We present new 890 $\mu m$ continuum ALMA observations of 5 brown dwarfs (BDs) with infrared excess in Lupus I and III -- which, in combination with 4 BDs previously observed, allowed us to study the mm properties of the full known BD disk population of one star-forming region. Emission is detected in 5 out of the 9 BD disks. Dust disk mass, brightness profiles and characteristic sizes of the BD population are inferred from continuum flux and modeling of the observations. Only one source is marginally resolved, allowing for the determination of its disk characteristic size. We conduct a demographic comparison between the properties of disks around BDs and stars in Lupus. Due to the small sample size, we cannot confirm or disprove if the disk mass over stellar mass ratio drops for BDs, as suggested for Ophiuchus. Nevertheless, we find that all detected BD disks have an estimated dust mass between 0.2 and 3.2 $M_{\bigoplus}$; these results suggest that the measured solid masses in BD disks can not explain the observed exoplanet population, analogous to earlier findings on disks around more massive stars. Combined with the low estimated accretion rates, and assuming that the mm-continuum emission is a reliable proxy for the total disk mass, we derive ratios of $\dot{M}_{\mathrm{acc}} / M_{\mathrm{disk}}$ significantly lower than in disks around more massive stars. If confirmed with more accurate measurements of disk gas masses, this result could imply a qualitatively different relationship between disk masses and inward gas transport in BD disks.

Photoevaporation of planet forming discs by high energy radiation from the central star is potentially a crucial mechanism for disc evolution and it may play an important role in the formation and evolution of planetary system. We present here a new generation of X-ray photoevaporation models for solar-type stars, based on a new set of hydrodynamical simulations, which account for stellar irradiation via a new, significantly improved, parameterisation of gas temperatures, based on detailed photoionisation and radiation transfer calculations. This is the first of a series of papers aiming at providing a library of models which cover the observed parameter space in stellar and disc mass, metallicity and stellar X-ray properties. We focus here on solar-type stars (0.7 M$_\odot$) with relatively low-mass discs (1% of the stellar mass) and explore the dependence of the wind mass loss rates on stellar X-ray luminosity. We model primordial discs as well as transition discs at various stages of evolution. Our 2D hydrodynamical models are then used to derive simple recipes for the mass loss rates that are suitable for inclusion in one-dimensional disc evolution and/or planet formation models typically employed for population synthesis studies. Line profiles from typical wind diagnostics ([OI] 6300 $\overset{\lower.5em\circ}{\mathrm{A}}$ and [NeII] 12.8 $\mu$m) are also calculated for our models and found to be roughly in agreement with previous studies. Finally, we perform a population study of transition discs by means of one-dimensional viscous evolution models including our new photoevaporation prescription and find that roughly a half of observed transition discs cavities and accretion rates could be reproduced by our models.

The so-called transition discs provide an important tool to probe various mechanisms that might influence the evolution of protoplanetary discs and therefore the formation of planetary systems. One of these mechanisms is photoevaporation due to energetic radiation from the central star, which can in principal explain the occurrence of discs with inner cavities like transition discs. Current models, however, fail to reproduce a subset of the observed transition discs, namely objects with large measured cavities and vigorous accretion. For these objects the presence of (multiple) giant planets is often invoked to explain the observations. In our work we explore the possibility of X-ray photoevaporation operating in discs with different gas-phase depletion of carbon and show that the influence of photoevaporation can be extended in such low-metallicity discs. As carbon is one of the main contributors to the X-ray opacity, its depletion leads to larger penetration depths of X-rays in the disc and results in higher gas temperatures and stronger photoevaporative winds. We present radiation-hydrodynamical models of discs irradiated by internal X-ray+EUV radiation assuming Carbon gas-phase depletions by factors of 3,10 and 100 and derive realistic mass-loss rates and profiles. Our analysis yields robust temperature prescriptions as well as photoevaporative mass-loss rates and profiles which may be able to explain a larger fraction of the observed diversity of transition discs.

The combination of high resolution and sensitivity offered by ALMA is revolutionizing our understanding of protoplanetary discs, as their bulk gas and dust distributions can be studied independently. In this paper we present resolved ALMA observations of the continuum emission ($\lambda=1.3$ mm) and CO isotopologues ($^{12}$CO, $^{13}$CO, C$^{18}$O $J=2-1$) integrated intensity from the disc around the nearby ($d = 162$ pc), intermediate mass ($M_{\star}=1.67\,M_{\odot}$) pre-main-sequence star CQ Tau. The data show an inner depression in continuum, and in both $^{13}$CO and C$^{18}$O emission. We employ a thermo-chemical model of the disc reproducing both continuum and gas radial intensity profiles, together with the disc SED. The models show that a gas inner cavity with size between 15 and 25 au is needed to reproduce the data with a density depletion factor between $\sim 10^{-1}$ and $\sim 10^{-3}$. The radial profile of the distinct cavity in the dust continuum is described by a Gaussian ring centered at $R_{\rm dust}=53\,$au and with a width of $\sigma=13\,$au. Three dimensional gas and dust numerical simulations of a disc with an embedded planet at a separation from the central star of $\sim20\,$au and with a mass of $\sim 6\textrm{-} 9\,M_{\rm Jup}$ reproduce qualitatively the gas and dust profiles of the CQ Tau disc. However, a one planet model appears not to be able to reproduce the dust Gaussian density profile predicted using the thermo-chemical modeling.

From molecular clouds to protoplanetary disks, non-ideal magnetic effects are important in many astrophysical environments. Indeed, in star and disk formation processes, it has become clear that these effects are critical to the evolution of the system. The efficacy of non-ideal effects are, however, determined by the complex interplay between magnetic fields, ionising radiation, cosmic rays, microphysics, and chemistry. In order to understand these key microphysical parameters, we present a one-dimensional non-ideal magnetohydrodynamics code and apply it to a model of a time-dependent, oblique, magnetic shock wave. By varying the microphysical ingredients of the model, we find that cosmic rays and dust play a major role, and that, despite the uncertainties, the inclusion of microphysics is essential to obtain a realistic outcome in magnetic astrophysical simulations.

High energy radiation from a planet host star can have strong influence on the final habitability of a system through several mechanisms. In this context we have constructed a catalogue containing the X-ray luminosities, as well as basic stellar and planetary properties of all known stars hosting giant planets (> 0.1 Mjup) that have been observed by the Chandra X-ray Observatory, XMM-Newton and/or ROSAT. Specifically in this paper we present a first application of this catalogue to search for a possible imprint of X-ray photoevaporation of planet-forming discs on the present-day orbital distribution of the observed giant planets. We found a suggestive void in the semi-major axis, a, versus X-ray luminosity, Lx, plane, roughly located between a ~ 0.05-1 au and Lx ~10^27-10^29 erg/s, which would be expected if photoevaporation played a dominant role in the migration history of these systems. However, due to the small observational sample size, the statistical significance of this feature cannot be proven at this point.

We present combined observations of the NH3 (J,K) = (1,1) and (2,2) inversion transitions towards OMC1 in Orion A obtained by the Karl G. Jansky Very Large Array (VLA) and the 100 m Robert C. Byrd Green Bank Telescope (GBT). With an angular resolution of 6" (0.01 pc), these observations reveal with unprecedented detail the complex filamentary structure extending north of the active Orion BN/KL region in a field covering 6' x 7'. We find a 0.012 pc wide filament within OMC1, with an aspect ratio of ~37:1, that was missed in previous studies. Its orientation is directly compared to the relative orientation of the magnetic field from the James Clerk Maxwell Telescope BISTRO survey in Orion A. We find a small deviation of ~11 deg between the mean orientation of the filament and the magnetic field, suggesting that they are almost parallel to one another. The filament's column density is estimated to be 2-3 orders of magnitude larger than the filaments studied with Herschel and is possibly self-gravitating given the low values of turbulence found. We further produce maps of the gas kinematics by forward modeling the hyperfine structure of the NH3 (J,K) = (1,1) and (2,2) lines. The resulting distribution of velocity dispersions peaks at ~0.5 km/s, close to the subsonic regime of the gas. This value is about 0.2 km/s smaller than previously measured in single-dish observations of the same region, suggesting that higher angular and spectral resolution observations will identify even lower velocity dispersions that might reach the subsonic turbulence regime in dense gas filaments.

We analyze the X-ray spectra of the $\sim$8000 sources detected in the Cygnus OB2 Chandra Legacy Survey (Drake et al., this issue), with the goals of characterizing the coronal plasma of the young low-mass stars in the region and estimating their intrinsic X-ray luminosities. We adopt two different strategies for X-ray sources for which more or less than 20 photons were detected. For the brighter sample we fit the spectra with absorbed isothermal models. In order to limit uncertainties, for most of the fainter Cygnus OB2 members in this sample, we constrain the spectral parameters to characteristic ranges defined from the brightest stars. For X-ray sources with $<$20 net photons we adopt a conversion factor from detected photon flux to intrinsic flux. This was defined, building on the results for the previous sample, as a function of the 20% quantile of the detected photon energy distributions, which we prove to also correlate well with extinction. We then use the X-ray extinction from the spectral fits to constrain the ratio between optical and X-ray extinction toward Cygnus OB2, finding it consistent with standard "Galactic"' values, when properly accounting for systematics. Finally we exploit the large number of sources to constrain the average coronal abundances of several elements, through two different ensemble analyses of the X-ray spectra of low-mass Cygnus OB2 members. We find the pattern of abundances to be largely consistent with that derived for the young stellar coronae in the Orion Nebula Cluster.

Gas giants' early ($\lesssim 5$ Myr) orbital evolution occurs in a disc losing mass in part to photoevaporation driven by high energy irradiance from the host star. This process may ultimately overcome viscous accretion to disperse the disc and halt migrating giants by starving their orbits of gas, imprinting on giant planet separations in evolved systems. Inversion of this distribution could then give insight into whether stellar FUV, EUV or X-ray flux dominates photoevaporation, constraining planet formation and disc evolution models. We use a 1D hydrodynamic code in population syntheses for gas giants undergoing Type II migration in a viscously evolving disc subject to either a primarily FUV, EUV or X-ray flux from a pre-solar T Tauri star. The photoevaporative mass loss profile's unique peak location and width in each energetic regime produces characteristic features in the distribution of giant separations: a severe dearth of $\lesssim$ 2 M$_{\rm J}$ planets interior to 5 AU in the FUV scenario, a sharp concentration of $\lesssim$ 3 M$_{\rm J}$ planets between $\approx 1.5 - 2$ AU in the EUV case, and a relative abundance of $\approx 2 - 3.5$ M$_{\rm J}$ giants interior to 0.5 AU in the X-ray model. These features do not resemble the observational sample of gas giants with mass constraints, though our results do show some weaker qualitative similarities. We thus assess how the differing photoevaporative profiles interact with migrating giants and address the effects of large model uncertainties as a step to better connect disc models with trends in the exoplanet population.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

Circumstellar discs with large dust depleted cavities and vigorous accretion onto the central star are often considered signposts for (multiple) giant planet formation. In this letter we show that X-ray photoevaporation operating in discs with modest (factors 3-10) gas-phase depletion of Carbon and Oxygen at large radii (> 15 AU) yield the inner radius and accretion rates for most of the observed discs, without the need to invoke giant planet formation. We present one-dimensional viscous evolution models of discs affected by X-ray photoevaporation assuming moderate gas-phase depletion of Carbon and Oxygen, well within the range reported by recent observations. Our models use a simplified prescription for scaling the X-ray photoevaporation rates and profiles at different metallicity, and our quantitative result depends on this scaling. While more rigorous hydrodynamical modelling of mass loss profiles at low metallicities is required to constrain the observational parameter space that can be explained by our models, the general conclusion that metal sequestering at large radii may be responsible for the observed diversity of transition discs is shown to be robust. Gap opening by giant planet formation may still be responsible for a number of observed transition discs with large cavities and very high accretion rate.

The streaming instability is often invoked as solution to the fragmentation and drift barriers in planetesimal formation, catalyzing the aggregation of dust on kyr timescales to grow km-sized cores. However there remains a lack of consensus on the physical mechanism(s) responsible for initiating it. One potential avenue is disc photoevaporation, wherein the preferential removal of relatively dust-free gas increases the disc metallicity. Late in the disc lifetime, photoevaporation dominates viscous accretion, creating a gradient in the depleted gas surface density near the location of the gap. This induces a local pressure maximum that collects drifting dust particles, which may then become susceptible to the streaming instability. Using a one-dimensional viscous evolution model of a disc subject to internal X-ray photoevaporation, we explore the efficacy of this process to build planetestimals. Over a range of parameters we find that the amount of dust mass converted into planetesimals is often < 1 M_Earth and at most a few M_Earth spread across tens of AU. We conclude that photoevaporation may at best be relevant for the formation of debris discs, rather than a common mechanism for the formation of planetary cores. Our results are in contrast to a recent, similar investigation that considered an FUV-driven photoevaporation model and reported the formation of tens of M_Earth at large (> 100 AU) disc radii. The discrepancies are primarily a consequence of the different photoevaporation profiles assumed. Until observations more tightly constrain photoevaporation models, the relevance of this process to the formation of planets remains uncertain.

(Abridged) We explore the role of X-ray photoevaporation in the evolution and dispersal of viscously evolving T-Tauri discs. We show that the X-ray photoevaporation wind rates scale linearly with X-ray luminosity, such that the observed range of X-ray luminosities for solar-type T-Tauri stars (10e28-10e31 erg\s) gives rise to vigorous disc winds with rates of order 10e-10-10e-7 M_sun/yr. We use the wind solutions from radiation-hydrodynamic models, coupled to a viscous evolution model to construct a population synthesis model so that we may study the physical properties of evolving discs and so-called `transition discs'. Current observations of disc lifetimes and accretion rates can be matched by our model assuming a viscosity parameter alpha = 2.5e-3. Our models confirm that X-rays play a dominant role in the evolution and dispersal of protoplanetary discs giving rise to the observed diverse population of inner hole `transition' sources which include those with massive outer discs, those with gas in their inner holes and those with detectable accretion signatures. To help understand the nature of observed transition discs we present a diagnostic diagram based on accretion rates versus inner hole sizes that demonstrate that, contrary to recent claims, many of the observed accreting and non accreting transition discs can easily be explained by X-ray photoevaporation. Finally, we confirm the conjecture of Drake et al. (2009), that accretion is suppressed by the X-rays through `photoevaporation starved accretion' and predict this effect can give rise to a negative correlation between X-ray luminosity and accretion rate, as reported in the Orion data.

Discs of gas and dust around Myr-old stars are a by-product of the star formation process and provide the raw material to form planets. Hence, their evolution and dispersal directly impact what type of planets can form and affect the final architecture of planetary systems. Here, we review empirical constraints on disc evolution and dispersal with special emphasis on transition discs, a subset of discs that appear to be caught in the act of clearing out planet-forming material. Along with observations, we summarize theoretical models that build our physical understanding of how discs evolve and disperse and discuss their significance in the context of the formation and evolution of planetary systems. By confronting theoretical predictions with observations, we also identify the most promising areas for future progress.

The dispersal of the circumstellar discs of dust and gas surrounding young low- mass stars has important implications for the formation of planetary systems. Photo- evaporation from energetic radiation from the central object is thought to drive the dispersal in the majority of discs, by creating a gap which disconnects the outer from the inner regions of the disc and then disperses the outer disc from the inside-out, while the inner disc keeps draining viscously onto the star. In this Letter we show that the disc around TW Hya, the closest protoplanetary disc to Earth, may be the first object where a photoevaporative gap has been imaged around the time at which it is being created. Indeed the detected gap in the ALMA images is consistent with the expectations of X-ray photoevaporation models, thus not requiring the presence of a planet. The photoevaporation model is also consistent with a broad range of properties of the TW Hya system, e.g. accretion rate and the location of the gap at the onset of dispersal. We show that the central, unresolved 870 {\mu}m continuum source might be produced by free free emission from the gas and/or residual dust inside the gap.

Using broadband photometry from the Hubble Space Telescope in combination with Very Large Telescope narrowband Halpha observations of the nearby spiral galaxy NGC 300, we explore a method for estimating the escape fractions of hydrogen-ionizing photons from HII regions within this galaxy. Our goal in this concept study is to evaluate the spectral types of the most massive stars using the broadband data and estimating their ionizing photon output with the help of stellar atmosphere models. A comparison with the Halpha flux that gives the amount of ionized gas in the HII region provides a measure of the escape fraction of ionizing photons from that region. We performed some tests with a number of synthetic young clusters with varying parameters to assess the reliability of the method. However, we found that the derived stellar spectral types and consequently the expected ionizing photon luminosity of a region is highly uncertain. The tests also show that on one hand we tended to overestimate the integrated photon output of a region for young ages and low numbers of stars, and on the other hand we mostly underestimated the combined ionizing luminosity for a large stellar number and older cluster ages. We conclude that the proposed method of using stellar broadband photometry to infer the leakage of ionizing photons from HII regions is highly uncertain and dominated by the errors of the resulting stellar spectral types. Therefore this method is not suitable. Stellar spectra are needed to reliably determine the stellar types and escape fractions. Studies to this end have been carried out for the Magellanic Clouds.

Photoevaporation of protoplanetary discs by high energy radiation from the central young stellar object is currently the favourite model to explain the sudden dispersal of discs from the inside out. While several theoretical works have provided a detailed pictured of this process, the direct observational validation is still lacking. Emission lines produced in these slow moving protoplanetary disc winds may bear the imprint of the wind structure and thus provide a potential diagnostic of the underlying dispersal process. In this paper we primarily focus on the collisionally excited neutral oxygen line at 6300A. We compare our models predictions to observational data and demonstrate a thermal origin for the observed blueshifted low-velocity component of this line from protoplanetary discs. Furthermore our models show that while this line is a clear tell-tale-sign of a warm, quasi-neutral disc wind, typical of X-ray photoevaporation, its strong temperature dependence makes it unsuitable to measure detailed wind quantities like mass-loss-rate.

In our Galaxy, star formation occurs in a variety of environments, with a large fraction of stars formed in clusters hosting massive stars. OB stars have an important feedback on the evolution of protoplanetary disks around nearby young stars and likely on the process of planet formation occurring in them. The nearby massive association Cygnus OB2 is an outstanding laboratory to study this feedback. It is the closest massive association to our Sun, and hosts hundreds of massive stars and thousands of low mass members. In this paper, we analyze the spatial variation of the disk fraction in Cygnus OB2 and we study its correlation with the local values of Far and Extreme ultraviolet radiation fields and the local stellar surface density. We present definitive evidence that disks are more rapidly dissipated in the regions of the association characterized by intense local UV field and large stellar density. In particular, the FUV radiation dominates disks dissipation timescales in the proximity (i.e. within 0.5 pc) of the O stars. In the rest of the association, EUV photons potentially induce a significant mass loss from the irradiated disks across the entire association, but the efficiency of this process is reduced at increasing distances from the massive stars due to absorption by the intervening intracluster material. We find that disk dissipation due to close stellar encounters is negligible in Cygnus OB2, and likely to have affected 1% or fewer of the stellar population. Disk dissipation is instead dominated by photoevaporation. We also compare our results to what has been found in other young clusters with different massive populations, concluding that massive associations like Cygnus OB2 are potentially hostile to protoplanetary disks, but that the environments where disks can safely evolve in planetary systems are likely quite common in our Galaxy.

Theoretical models of grain growth predict dust properties to change as a function of protoplanetary disk radius, mass, age and other physical conditions. We lay down the methodology for a multi-wavelength analysis of (sub-)mm and cm continuum interferometric observations to constrain self-consistently the disk structure and the radial variation of the dust properties. The computational architecture is massively parallel and highly modular. The analysis is based on the simultaneous fit in the uv-plane of observations at several wavelengths with a model for the disk thermal emission and for the dust opacity. The observed flux density at the different wavelengths is fitted by posing constraints on the disk structure and on the radial variation of the grain size distribution. We apply the analysis to observations of three protoplanetary disks (AS 209, FT Tau, DR Tau) for which a combination of spatially resolved observations in the range ~0.88mm to ~10mm is available (from SMA, CARMA, and VLA), finding evidence of a decreasing maximum dust grain size (a_max) with radius. We derive large a_max values up to 1 cm in the inner disk between 15 and 30 AU and smaller grains with a_max~1 mm in the outer disk (R > 80AU). In this paper we develop a multi-wavelength analysis that will allow this missing quantity to be constrained for statistically relevant samples of disks and to investigate possible correlations with disk or stellar parameters.

We present hydrodynamic simulations of stellar wind-magnetosphere interactions in hot Jupiters such as WASP-12b. For fiducial stellar wind rates we find that a planetary magnetic field of a few G produces a large magnetospheric cavity, which is typically 6-9 planetary radii in size. A bow shock invariably forms ahead of the magnetosphere, but the pre-shock gas is only mildly supersonic (with typical Mach numbers of $\simeq$1.6-1.8) so the shock is weak. This results in a characteristic signature in the ultraviolet light curve: a broad absorption feature that leads the optical transit by 10-20% in orbital phase. The shapes of our synthetic light-curves are consistent with existing observations of WASP-12b, but the required near-UV optical depth ($\tau \sim 0.1$) can only be achieved if the shocked gas cools rapidly. We further show that radiative cooling is inefficient, so we deem it unlikely that a magnetospheric bow shock is responsible for the observed near-UV absorption. Finally, we apply our model to two other well-studied hot Jupiters (WASP-18b and HD209458b), and suggest that UV observations of more massive short-period planets (such as WASP-18b) will provide a straightforward test to distinguish between different models of circumplanetary absorption.

Ionising feedback from massive stars dramatically affects the interstellar medium local to star forming regions. Numerical simulations are now starting to include enough complexity to produce morphologies and gas properties that are not too dissimilar from observations. The comparison between the density fields produced by hydrodynamical simulations and observations at given wavelengths relies however on photoionisation/chemistry and radiative transfer calculations. We present here an implementation of Monte Carlo radiation transport through a Voronoi tessellation in the photoionisation and dust radiative transfer code MOCASSIN. We show for the first time a synthetic spectrum and synthetic emission line maps of an hydrodynamical simulation of a molecular cloud affected by massive stellar feedback. We show that the approach on which previous work is based, which remapped hydrodynamical density fields onto Cartesian grids before performing radiative transfer/photoionisation calculations, results in significant errors in the temperature and ionisation structure of the region. Furthermore, we describe the mathematical process of tracing photon energy packets through a Voronoi tessellation, including optimisations, treating problematic cases and boundary conditions. We perform various benchmarks using both the original version of MOCASSIN and the modified version using the Voronoi tessellation. We show that for uniform grids, or equivalently a cubic lattice of cell generating points, the new Voronoi version gives the same results as the original Cartesian-grid version of MOCASSIN for all benchmarks. For non-uniform initial conditions, such as using snapshots from Smoothed Particle Hydrodynamics simulations, we show that the Voronoi version performs better than the Cartesian grid version, resulting in much better resolution in dense regions.

Most star clusters at an intermediate age (1-2 Gyr) in the Large and Small Magellanic Clouds show a puzzling feature in their color-magnitude diagrams (CMD) that is not in agreement with a simple stellar population. The main sequence turn-off of these clusters is much broader than would be expected from photometric uncertainties. One interpretation of this feature is that age spreads of the order 200-500 Myr exist within individual clusters, although this interpretation is highly debated. Such large age spreads should affect other parts of the CMD, which are sensitive to age, as well. In this study, we analyze the CMDs of a sample of 12 intermediate-age clusters in the Large Magellanic Cloud that all show an extended turn-off using archival optical data taken with the Hubble Space Telescope. We fit the star formation history of the turn-off region and the red clump region independently with two different theoretical isochrone models. We find that in most of the cases, the age spreads inferred from the red clumps are smaller than the ones resulting from the turn-off region. However, the age ranges resulting from the red clump region are broader than would be expected for a single age. Only two out of 12 clusters in our sample show a red clump which seems to be consistent with a single age. As our results are not unambiguous, we can not ultimately tell if the extended main sequence turn-off feature is due to an age spread, or not, by fitting the star formation histories to the red clump regions. However, we find that the width of the extended main sequence turn-off feature is correlated with the age of the clusters in a way which would be unexplained in the "age spread" interpretation, but which may be expected if stellar rotation is the cause of the spread at the turn-off.

Photo-evaporation and planet formation have both been proposed as mechanisms responsible for the creation of a transition disc. We have studied their combined effect through a suite of 2d simulations of protoplanetary discs undergoing X-ray photoevaporation with an embedded giant planet. In a previous work we explored how the formation of a giant planet triggers the dispersal of the inner disc by photo-evaporation at earlier times than what would have happened otherwise. This is particularly relevant for the observed transition discs with large holes and high mass accretion rates that cannot be explained by photo-evaporation alone. In this work we significantly expand the parameter space investigated by previous simulations. In addition, the updated model includes thermal sweeping, needed for studying the complete dispersal of the disc. After the removal of the inner disc the disc is a non accreting transition disc, an object that is rarely seen in observations. We assess the relative length of this phase, to understand if it is long lived enough to be found observationally. Depending on the parameters, especially on the X-ray luminosity of the star, we find that the fraction of time spent as a non-accretor greatly varies. We build a population synthesis model to compare with observations and find that in general thermal sweeping is not effective enough to destroy the outer disc, leaving many transition discs in a relatively long lived phase with a gas free hole, at odds with observations. We discuss the implications for transition disc evolution. In particular, we highlight the current lack of explanation for the missing non-accreting transition discs with large holes, which is a serious issue in the planet hypothesis.

By means of radiative transfer simulation we study the evolution of the far-infrared colours of protoplanetary discs undergoing inside-out dispersal, often referred to as transition discs. We show that a brightening of the mid and far-infrared emission from these objects is a natural consequence of the removal of the inner disc. Our results can fully explain recent observations of transition discs in the Chamaleon and Lupus star forming regions from the Herschel Gould Belt Survey, which show a higher median for the 70?um (Herschel PACS 1) band of known transition objects compared with primordial discs. Our theoretical results hence support the suggestion that the 70?um band may be a powerful diagnostic for the identification of transition discs from photometry data, provided that the inner hole is larger than tens of AU, depending on spectral type. Furthermore we show that a comparison of photometry in the K , 12?um and 7u0?m bands to model tracks can provide a rough, but quick estimate of the inner hole size of these objects, provided their inclination is below ?85 degrees and the inner hole size is again larger than tens of AU.

We examine the combined effects of winds and photoionizing radiation from O--type stars on embedded stellar clusters formed in model turbulent molecular clouds covering a range of masses and radii. We find that feedback is able to increase the quantities of dense gas present, but decreases the rate and efficiency of the conversion of gas to stars relative to control simulations in which feedback is absent. Star formation in these calculations often proceeds at a rate substantially slower than the freefall rate in the dense gas. This decoupling is due to the weakening of, and expulsion of gas from, the deepest parts of the clouds' potential wells where most of the star formation occurs in the control simulations. This results in large fractions of the stellar populations in the feedback simulation becoming dissociated from dense gas. However, where star formation \emph{does} occur in both control and feedback simulations, it does so in dense gas, so the correlation between star formation activity and dense gas is preserved. The overall dynamical effects of feedback on the \emph{clusters} are minimal, with only small fraction of stars becoming unbound, despite large quantities of gas being expelled from some clouds. This owes to the settling of the stars into virialised and stellar--dominated configurations before the onset of feedback. By contrast, the effects of feedback on the observable properties of the clusters -- their U--, B-- and V--band magnitudes -- are strong and sudden. The timescales on which the clusters become visible and unobscured are short compared with the timescales which the clouds are actually destroyed.

We investigate the influence of photoevaporation of protoplanetary discs on the final distribution of exoplanets semi-major axis distances. We model giant planet migration in viscous discs affected by photoevaporation driven by either pure EUV or soft X-ray radiation (XEUV). We show that the final exoplanet distributions are strongly dependant on the choice of the photoevaporation model. In particular, we find that XEUV is more efficient than pure EUV radiation at parking planets at approximately 1-2 AU distance from their central star, hence roughly reproducing the observed peak in the exoplanets semi-major axis distributions. We note however that a more quantitative comparison with the observations is hindered by the oversimplified treatment of planetary accretion, which severely affects migration rates. For this reason, caution should be used when using these models to constrain details of disc clearing and/or migration from the observations. Nevertheless our results indicate that disc dispersal by photoevaporation may be the main driver of the features in the exoplanets semi-major axis distribution observed by recent surveys.

Integral field unit (IFU) data of the iconic Pillars of Creation in M16 are presented. The ionisation structure of the pillars was studied in great detail over almost the entire visible wavelength range, and maps of the relevant physical parameters, e.g. extinction, electron density, electron temperature, line-of-sight velocity of the ionised and neutral gas are shown. In agreement with previous authors, we find that the pillar tips are being ionised and photo-evaporated by the massive members of the nearby cluster NGC 6611. They display a stratified ionisation structure where the emission lines peak in a descending order according to their ionisation energies. The IFU data allowed us to analyse the kinematics of the photo-evaporative flow in terms of the stratified ionisation structure, and we find that, in agreement with simulations, the photo-evaporative flow is traced by a blueshift in the position-velocity profile. The gas kinematics and ionisation structure have allowed us to produce a sketch of the 3D geometry of the Pillars, positioning the pillars with respect to the ionising cluster stars. We use a novel method to detect a previously unknown bipolar outflow at the tip of the middle pillar and suggest that it has an embedded protostar as its driving source. Furthermore we identify a candidate outflow in the leftmost pillar. With the derived physical parameters and ionic abundances, we estimate a mass loss rate due to the photo-evaporative flow of 70 M$_{\odot}$ Myr$^{-1}$ which yields an expected lifetime of approximately 3 Myr.

The planetary nebula SuWt 2 (PN G311.0+02.4), is an unusual object with a prominent, inclined central emission ellipse and faint bipolar extensions. It has two A-type stars in a proven binary system at the centre. However, the radiation from these two central stars is too soft to ionize the surrounding material leading to a so far fruitless search for the responsible ionizing source. Such a source is clearly required and has already been inferred to exist via an observed temporal variation of the centre-of-mass velocity of the A-type stars. Moreover, the ejected nebula is nitrogen-rich which raises question about the mass-loss process from a likely intermediate-mass progenitor. We use optical integral-field spectroscopy to study the emission lines of the inner nebula ring. This has enabled us to perform an empirical analysis of the optical collisionally excited lines, together with a fully three-dimensional photoionization modelling. Our empirical results are used to constrain the photoionization models, which determine the evolutionary stage of the responsible ionizing source and its likely progenitor. The time-scale for the evolutionary track of a hydrogen-rich model atmosphere is inconsistent with the dynamical age obtained for the ring. This suggests that the central star has undergone a very late thermal pulse. We conclude that the ionizing star could be hydrogen-deficient and compatible with what is known as a PG 1159-type star. The evolutionary tracks for the very late thermal pulse models imply a central star mass of ~ 0.64M${}_{\odot}$, which originated from a ~ 3M${}_{\odot}$ progenitor. The evolutionary time-scales suggest that the central star left the asymptotic giant branch about 25,000 years ago, which is consistent with the nebula's age.

Giant Molecular Clouds (GMCs) are observed to be turbulent, but theory shows that without a driving mechanism turbulence should quickly decay. The question arises by which mechanisms turbulence is driven or sustained. It has been shown that photoionising feedback from massive stars has an impact on the surrounding GMC and can for example create vast HII bubbles. We therefore address the question of whether turbulence is a consequence of this effect of feedback on the cloud. To investigate this, we analyse the velocity field of simulations of high mass star forming regions by studying velocity structure functions and power spectra. We find that clouds whose morphology is strongly affected by photoionising feedback also show evidence of driving of turbulence by preserving or recovering a Kolmogorov-type velocity field. On the contrary, control run simulations without photoionising feedback have a velocity distribution that bears the signature of gravitational collapse and of the dissipation of energy, where the initial Kolmogorov-type structure function is erased.

We take the end result of smoothed particle hydrodynamics (SPH) simulations of star formation which include feedback from photoionisation and stellar winds and evolve them for a further 10Myr using $N$-body simulations. We compare the evolution of each simulation to a control run without feedback, and to a run with photoionisation feedback only. In common with previous work, we find that the presence of feedback prevents the runaway growth of massive stars, and the resulting star-forming regions are less dense, and preserve their initial substructure for longer. The addition of stellar winds to the feedback produces only marginal differences compared to the simulations with just photoionisation feedback. We search for mass segregation at different stages in the simulations; before feedback is switched on in the SPH runs, at the end of the SPH runs (before $N$-body integration) and during the $N$-body evolution. Whether a simulation is primordially mass segregated (i.e. before dynamical evolution) depends extensively on how mass segregation is defined, and different methods for measuring mass segregation give apparently contradictory results. Primordial mass segregation is also less common in the simulations when star formation occurs under the influence of feedback. Further dynamical mass segregation can also take place during the subsequent (gas-free) dynamical evolution. Taken together, our results suggest that extreme caution should be exercised when interpreting the spatial distribution of massive stars relative to low-mass stars in simulations.

The discovery with the {\it Herschel Space Observatory} of bright far infrared and submm emission from the ejecta of the core collapse supernova SN\,1987A has been interpreted as indicating the presence of some 0.4--0.7\,M$_\odot$ of dust. We have constructed radiative transfer models of the ejecta to fit optical to far-infrared observations from the literature at epochs between 615 days and 24 years after the explosion, to determine when and where this unexpectedly large amount of dust formed. We find that the observations by day 1153 are consistent with the presence of 3$\times$10$^{-3}$M$_\odot$ of dust. Although this is a larger amount than has previously been considered possible at this epoch, it is still very small compared to the amount present in the remnant after 24 years, and significantly higher dust masses at the earlier epochs are firmly ruled out by the observations, indicating that the majority of the dust must have formed at very late times. By 8515-9200 days after the explosion, 0.6--0.8\,M$_\odot$ of dust is present, and dust grains with radii greater than 2\,$\mu$m are required to obtain a fit to the observed SED. This suggests that the dust mass increase at late times was caused by accretion onto and coagulation of the dust grains formed at earlier epochs. These findings provide further confirmation that core collapse supernovae can create large quantities of dust, and indicate that the reason for small dust masses being estimated in many cases is that the vast majority of the dust forms long after most supernovae have been detectable at mid-infrared wavelengths.

Transitional disks (TDs) are thought to be a late evolutionary stage of protoplanetary disks with dust depleted inner regions. The mechanism responsible for this depletion is still under debate. To constrain the models it is mandatory to have a good understanding of the properties of the gas content of the inner disk. Using X-Shooter broad band -UV to NIR- medium resolution spectroscopy we derive the stellar, accretion, and wind properties of a sample of 22 TDs. The analysis of these properties allows us to put strong constraints on the gas content in a region very close to the star (<0.2 AU) which is not accessible with any other observational technique. We fit the spectra with a self-consistent procedure to derive simultaneously SpT,Av,and mass accretion rates (Macc) of the targets. From forbidden emission lines we derive the wind properties of the targets. Comparing our findings to values for cTTs, we find that Macc and wind properties of 80% of the TDs in our sample, which is strongly biased towards strongly accreting objects, are comparable to those of cTTs. Thus, there are (at least) some TDs with Macc compatible with those of cTTs, irrespective of the size of the dust inner hole.Only in 2 cases Macc are much lower, while the wind properties are similar. We do not see any strong trend of Macc with the size of the dust depleted cavity, nor with the presence of a dusty optically thick disk close to the star. In the TDs in our sample there is a gas rich inner disk with density similar to that of cTTs disks. At least for some TDs, the process responsible of the inner disk clearing should allow for a transfer of gas from the outer disk to the inner region. This should proceed at a rate that does not depend on the physical mechanism producing the gap seen in the dust emission and results in a gas density in the inner disk similar to that of unperturbed disks around stars of similar mass.

We model the combined effects of photoionization and momentum--driven winds from O--stars on molecular clouds spanning a parameter space of initial conditions. The dynamical effects of the winds are very modest. However, in the lower--mass clouds, they influence the morphologies of the HII regions by creating 10pc--scale central cavities.\\ The inhomogeneous structures of the model GMCs make them highly permeable to photons, ionized gas and supernova ejecta, and the leaking of ionized gas in particular strongly affects their evolution, reducing the effectiveness of feedback. Nevertheless, feedback is able to expel large fractions of the mass of the lower escape--velocity clouds. Its impact on star formation is more modest, decreasing final star formation efficiencies by 10--20$\%$, and the rate of change of the star formation efficiency per freefall time by about one third. However, the clouds still form stars substantially faster than observed GMCs.

Most stars form in a clustered environment. Therefore, it is important to assess how this environment influences the evolution of protoplanetary discs around young stars. In turn, this affects their ability to produce planets and ultimately life. We present here for the first time 3D SPH/N-body simulations that include both the hydrodynamical evolution of the discs around their natal stars, as well as the dynamics of the stars themselves. The discs are viscously evolving, accreting mass onto the central star and spreading. We find penetrating encounters to be very destructive for the discs as in previous studies, although the frequency of such encounters is low. We also find, however, that encounter influence the disc radii more strongly than other disc properties such as the disc mass. The disc sizes are set by the competition between viscous spreading and the disruptive effect of encounters. As discs spread, encounters become more and more important. In the regime of rapid spreading encounters simply truncate the discs, stripping the outer portions. In the opposite regime, we find that the effect of many distant encounters is able to limit the disc size. Finally, we predict from our simulations that disc sizes are limited by encounters at stellar densities exceeding $\sim 2-3 \times 10^3 \ \mathrm{pc}^{-2}$.

Recent observations reveal that the central star of the planetary nebula Abell 48 exhibits spectral features similar to massive nitrogen-sequence Wolf-Rayet stars. This raises a pertinent question, whether it is still a planetary nebula or rather a ring nebula of a massive star. In this study, we have constructed a three-dimensional photoionization model of Abell 48, constrained by our new optical integral field spectroscopy. An analysis of the spatially resolved velocity distributions allowed us to constrain the geometry of Abell 48. We used the collisionally excited lines to obtain the nebular physical conditions and ionic abundances of nitrogen, oxygen, neon, sulphur and argon, relative to hydrogen. We also determined helium temperatures and ionic abundances of helium and carbon from the optical recombination lines. We obtained a good fit to the observations for most of the emission-line fluxes in our photoionization model. The ionic abundances deduced from our model are in decent agreement with those derived by the empirical analysis. However, we notice obvious discrepancies between helium temperatures derived from the model and the empirical analysis, as overestimated by our model. This could be due to the presence of a small fraction of cold metal-rich structures, which were not included in our model. It is found that the observed nebular line fluxes were best reproduced by using a hydrogen-deficient expanding model atmosphere as the ionizing source with an effective temperature of 70 kK and a stellar luminosity of 5500 L_sun, which corresponds to a relatively low-mass progenitor star (~ 3 M_sun) rather than a massive Pop I star.

We analyse current measurements of accretion rates onto pre-main sequence stars as a function of stellar mass, and conclude that the steep dependance of accretion rates on stellar mass is real and not driven by selection/detection threshold, as has been previously feared. These conclusions are reached by means of statistical tests including a survival analysis which can account for upper limits. The power-law slope of the Mdot-M* relation is found to be in the range of 1.6-1.9 for young stars with masses lower than 1 Msun. The measured slopes and distributions can be easily reproduced by means of a simple disc model which includes viscous accretion and X-ray photoevaporation. We conclude that the Mdot-M* relation in pre-main sequence stars bears the signature of disc dispersal by X-ray photoevaporation, suggesting that the relation is a straight- forward consequence of disc physics rather than an imprint of initial conditions.

We examine the effect of momentum-driven OB-star stellar winds on a parameter space of simulated turbulent Giant Molecular Clouds using SPH hydrodynamical simulations. By comparison with identical simulations in which ionizing radiation was included instead of winds, we show that momentum-driven winds are considerably less effective in disrupting their host clouds than are HII regions. The wind bubbles produced are smaller and generally smoother than the corresponding ionization-driven bubbles. Winds are roughly as effective in destroying the very dense gas in which the O-stars are embedded, and thus shutting down the main regions of star-forming activity in the model clouds. However, their influence falls off rapidly with distance from the sources, so they are not as good at sweeping up dense gas and triggering star formation further out in the clouds. As a result, their effect on the star formation rate and efficiency is generally more negative than that of ionization, if they exert any effect at all.

Low-mass young stellar objects are powerful emitters of X-rays that can ionize and heat the disks and the young planets they harbour. The X-rays produce molecular ions that affect the chemistry of the disk atmospheres and their spectroscopic signatures. Deeper down, X-rays are the main ionization source and influence the operation of the magnetorotational instability, believed to be the main driver for the angular momentum redistribution crucial for the accretion and formation of these pre main-sequence stars. X-ray ionization also affects the character of the dead zones around the disk midplane where terrestrial planets are likely to form. To obtain the physical and chemical effects of the stellar X-rays, their propagation through the disk has to be calculated taking into account both absorption and scattering. To date the only calculation of this type was done almost 15 years ago, and here we present new three-dimensional radiative transfer calculations of X-ray ionization rates in protoplanetary discs. Our study confirms the results from previous work for the same physical parameters. It also updates them by including a more detailed treatment of the radiative transfer and by using ionizing spectra and elemental abundance more appropriate for what is currently known about protoplanetary disks and their host stars. The new calculations for a typical ionizing spectrum yield respectively lower and higher ionisation rates at high and low column densities at a given radius in a disc. The differences can be up to an order of magnitude near 1 AU, depending on the abundances used

We consider continuum free-free radio emission from the upper atmosphere of protoplanetary discs as a probe of the ionized luminosity impinging upon the disc. Making use of previously computed hydrodynamic models of disc photoevaporation within the framework of EUV and X-ray irradiation, we use radiative transfer post-processing techniques to predict the expected free-free emission from protoplanetary discs. In general, the free-free luminosity scales roughly linearly with ionizing luminosity in both EUV and X-ray driven scenarios, where the emission dominates over the dust tail of the disc and is partial optically thin at cm wavelengths. We perform a test observation of GM Aur at 14-18 Ghz and detect an excess of radio emission above the dust tail to a very high level of confidence. The observed flux density and spectral index are consistent with free-free emission from the ionized disc in either the EUV or X-ray driven scenario. Finally, we suggest a possible route to testing the EUV and X-ray driven dispersal model of protoplanetary discs, by combining observed free-free flux densities with measurements of mass-accretion rates. On the point of disc dispersal one would expect to find a M_dot^2 scaling with free-free flux in the case of EUV driven disc dispersal or a M_dot scaling in the case of X-ray driven disc dispersal.

We present optical integral field unit (IFU) observations of two gas pillars surrounding the Galactic young massive star cluster NGC 3603. The high S/N and spectral resolution of these data have allowed us to accurately quantify the H-alpha, [NII] and [SII] emission line shapes, and we find a mixture of broad (FWHM~70-100 km/s) and narrow (<50 km/s) components. The broad components are found close to the edges of both pillars, suggesting that they originate in turbulent mixing layers (TMLs) driven by the effect of the star cluster wind. Both pillars exhibit surprisingly high ionized gas densities of >10000 cm-3. In one pillar we found that these high densities are only found in the narrow component, implying it must originate from deeper within the pillar than the broad component. From this, together with our kinematical data, we conclude that the narrow component traces a photoevaporation flow, and that the TML forms at the interface with the hot wind. On the pillar surfaces we find a consistent offset in radial velocity between the narrow (brighter) components of H-alpha and [NII] of ~5-8 km/s, for which we were unable to find a satisfactory explanation. We urge the theoretical community to simulate mechanical and radiative cloud interactions in more detail to address the many unanswered questions raised by this study.

Emission line studies from Hii regions in galaxies require tools for the inversion of line ratios into desired physical properties. These tools generally come in form of diagnostic ratios/diagrams that are based on grids of photoionisation models. An important input to the photoionisation models is the stellar atmosphere spectrum of the ionising sources. The current omission of shocks in the calculation of the former set of models could threaten the accuracy of the physical interpretation of emission line ratios from Hii regions. Current stellar atmosphere models that are crucial inputs to the grid of photoionisation models used to generate nebular emission line diagnostic diagrams might produce significant biases due to the omission of shocks. We investigate whether a new generation of photoionisation model grids, taking shocks into account, is required to compensate for the biases. We make use of the WM-Basic stellar atmosphere code, which can account for the extra energetic emission in the stellar spectral energy distribution produced by shocks, in conjunction with the photoionisation code MOCASSIN to determine whether shocks produce significant biases in the determination of the physical parameters of the interstellar medium and/or ionising stellar parameters. We conclude that these effects are only important for stellar sources with effective temperatures lower than 30kK and in this case they yield artificially high stellar temperatures, electron temperatures and nebular ionisation parameters. The magnitude of the effect is also obviously dependent on the strength of the shock and is likely to be unimportant for the majority of stellar sources. Nevertheless, we find our 20kK and 30kK shock models to strongly enhance the He ii 4686 nebular emission line. This result is however not strong enough to explain previously observed He ii 4686 line emission in the spectra of Hii galaxies.

Gaseous giant planet formation is thought to occur in the first few million years following stellar birth. Models predict that giant planet formation carves a deep gap in the dust component (shallower in the gas). Infrared observations of the disk around the young star HD142527, at ~140pc, found an inner disk ~10AU in radius, surrounded by a particularly large gap, with a disrupted outer disk beyond 140AU, indicative of a perturbing planetary-mass body at ~90 AU. From radio observations, the bulk mass is molecular and lies in the outer disk, whose continuum emission has a horseshoe morphology. The vigorous stellar accretion rate would deplete the inner disk in less than a year, so in order to sustain the observed accretion, matter must flow from the outer-disk into the cavity and cross the gap. In dynamical models, the putative protoplanets channel outer-disk material into gap-crossing bridges that feed stellar accretion through the inner disk. Here we report observations with the Atacama Large Millimetre Array (ALMA) that reveal diffuse CO gas inside the gap, with denser HCO+ gas along gap-crossing filaments, and that confirm the horseshoe morphology of the outer disk. The estimated flow rate of the gas is in the range 7E-9 to 2E-7 Msun/yr, which is sufficient to maintain accretion onto the star at the present rate.

We compare emission line intensities from photoionization models of smooth and fractal shell geometries for low density H II regions, with particular focus on the low-ionization diagnostic diagram [N II]/H-alpha vs H-alpha. Building on previously published models and observations of Barnard's Loop, we show that the observed range of intensities and variations in the line intensity ratios may be reproduced with a three dimensional shell geometry. Our models adopt solar abundances throughout the model nebula, in contrast with previous one dimensional modeling which suggested the variations in line intensity ratios could only be reproduced if the heavy element abundances were increased by a factor of 1.4. For spatially resolved H II regions, the multiple sightlines that pierce and sample different ionization and temperature conditions within smooth and fractal shells produce a range of line intensities that are easily overlooked if only the total integrated intensities from the entire nebula model are computed. Our conclusion is that inference of H II region properties, such as elemental abundances, via photoionization models of one dimensional geometries must be treated with caution and further tested through three dimensional modeling.

We present the fourth in a series of papers detailing our SPH study of the effects of ionizing feedback from O--type stars on turbulent star forming clouds. Here, we study the effects of photoionization on a series of initially partially unbound clouds with masses ranging from $10^{4}$--$10^{6}$M$_{\odot}$ and initial sizes from 2.5-45pc. We find that ionizing feedback profoundly affects the structure of the gas in most of our model clouds, creating large and often well-cleared bubble structures and pillars. However, changes in the structures of the embedded clusters produced are much weaker and not well correlated to the evolution of the gas. We find that in all cases, star formation efficiencies and rates are reduced by feedback and numbers of objects increased, relative to control simulations. We find that local triggered star formation does occur and that there is a good correlation between triggered objects and pillars or bubble walls, but that triggered objects are often spatially-mixed with those formed spontaneously. Some triggered objects acquire large enough masses to become ionizing sources themselves, lending support to the concept of propagating star formation. We find scant evidence for spatial age gradients in most simulations, and where we do see them, they are not a good indicator of triggering, as they apply equally to spontaneously-formed objects as triggered ones. Overall, we conclude that inferring the global or local effects of feedback on stellar populations from observing a system at a single epoch is very problematic.

We assess the potential of planet formation instigating the early formation of a photoevaporation driven gap, up to radii larger than typical for photoevaporation alone. For our investigation we make use of hydrodynamics models of photoevaporating discs with a giant planet embedded. We find that, by reducing the mass accretion flow onto the star, discs that form giant planets will be dispersed at earlier times than discs without planets by X-ray photoevaporation. By clearing the portion of the disc inner of the planet orbital radius, planet formation induced photoevaporation (PIPE) is able to produce transition disc that for a given mass accretion rate have larger holes when compared to standard X-ray photoevaporation. This constitutes a possible route for the formation of the observed class of accreting transition discs with large holes, which are otherwise difficult to explain by planet formation or photoevaporation alone. Moreover, assuming that a planet is able to filter dust completely, PIPE produces a transition disc with a large hole and may provide a mechanism to quickly shut down accretion. This process appears to be too slow however to explain the observed desert in the population of transition disc with large holes and low mass accretion rates.

We investigate the viability of the magnetorotational instability (MRI) in accretion disks around both solar-type stars and very low mass stars. In particular, we determine the disk regions where the MRI can be shut off either by Ohmic resistivity (the so-called Dead and Undead Zones) or by ampipolar diffusion (a region we term the Zombie Zone). We consider 2 stellar masses: Mstar = 0.7 and 0.1 Msun. In each case, we assume that: the disk surface density profile is that of a scaled Minimum Mass Solar Nebula, with Mdisk/Mstar ~ 0.01 as currently estimated; disk ionisation is driven primarily by stellar X-rays, complemented by cosmic rays and radionuclides; and the stellar X-ray luminosity scales with bolometric luminosity as Lx/Lstar ~ 10^-3.5, as observed. Ionization rates are calculated with the MOCCASIN code, and ionisation balance determined using a simplified chemical network, including well-mixed 0.1 um grains at various levels of depletion. We find that (1) ambipolar diffusion is the primary factor controlling MRI activity in disks around both solar-type and very low mass stars. Assuming that the MRI yields the maximum possible field strength at each radius, we further find that: (2) the MRI-active layer constitutes only ~ 5-10% of the total disk mass; (3) the accretion rate (Mdot) varies radially in both magnitude and sign (inward or outward), implying time-variable accretion as well as the creation of disk gaps and overdensities, with consequences for planet formation and migration; (4) achieving the empirical accretion rates in solar-type and very low mass stars requires a depletion of well-mixed small grains by a factor of 10-1000 relative to the standard dust-to-gas mass ratio of 10^-2; and (5) the current non-detection of polarized emission from field-aligned grains in the outer disk regions is consistent with active MRI at those radii.

We extend our previous SPH parameter study of the effects of photoionization from O-stars on star-forming clouds to include initially unbound clouds. We generate a set of model clouds in the mass range $10^{4}-10^{6}$M$_{\odot}$ with initial virial ratios $E_{\rm kin}/E_{\rm pot}$=2.3, allow them to form stars, and study the impact of the photoionizing radiation produced by the massive stars. We find that, on the 3Myr timescale before supernovae are expected to begin detonating, the fractions of mass expelled by ionizing feedback is a very strong function of the cloud escape velocities. High-mass clouds are largely unaffected dynamically, while lower-mass clouds have large fractions of their gas reserves expelled on this timescale. However, the fractions of stellar mass unbound are modest and significant portions of the unbound stars are so only because the clouds themselves are initially partially unbound. We find that ionization is much more able to create well-cleared bubbles in the unbound clouds, owing to their intrinsic expansion, but that the presence of such bubbles does not necessarily indicate that a given cloud has been strongly influenced by feedback. We also find, in common with the bound clouds from our earlier work, that many of the systems simulated here are highly porous to photons and supernova ejecta, and that most of them will likely survive their first supernova explosions.

The time-scale over which and the modality by which young stellar objects (YSOs) disperse their circumstellar discs dramatically influences the eventual formation and evolution of planetary systems. By means of extensive radiative transfer (RT) modelling, we have developed a new set of diagnostic diagrams in the infrared colour-colour plane (K-[24] vs. K-[8]), to aid with the classification of the evolutionary stage of YSOs from photometric observations. Our diagrams allow the differentiation of sources with unevolved (primordial) discs from those evolving according to different clearing scenarios (e.g. homologous depletion vs. inside-out dispersal), as well as from sources that have already lost their disc. Classification of over 1500 sources in 15 nearby star-forming regions reveals that approximately 39 % of the sources lie in the primordial disc region, whereas between 31 % and 32 % disperse from the inside-out and up to 22 % of the sources have already lost their disc. Less than 2 % of the objects in our sample lie in the homogeneous draining regime. Time-scales for the transition phase are estimated to be typically a few 10^5 years independent of stellar mass. Therefore, regardless of spectral type, we conclude that currently available infrared photometric surveys point to fast (of order 10 % of the global disc lifetime) inside-out clearing as the preferred mode of disc dispersal.

We investigate the effect of ionising radiation from Massive Young Stellar Objects impinging on their emerging spectral energy distribution. By means of detailed radiative transfer calculations including both the gaseous and dust phase of their surrounding protoplanetary discs we highlight the importance of modelling both phases simultaneously when interpreting observations from such objects. In particular we find that models that only include dust may lead to incorrect conclusions about the inner disc evolution. Furthermore the omission of gas from models overproduces far-infrared and sub-millimiter fluxes with the result that derived dust masses may be underestimated by a factor of two in some cases. Finally free-free emission from the ionised component of gaseous discs causes the slope of the dust emission in the sub-mm and mm regime to appear flatter, resulting in incorrectly modelling the dust properties, with consequences on the derived disc masses, power law index of the surface density profile and other disc properties.

We perform a comparison between two radiative transfer algorithms commonly employed in hydrodynamical calculations of star formation: grey flux limited diffusion and the hybrid scheme, in addition we compare these algorithms to results from the Monte-Carlo radiative transfer code MOCASSIN. In disc like density structures the hybrid scheme performs significantly better than the FLD method in the optically thin regions, with comparable results in optically thick regions. In the case of a forming high mass star we find the FLD method significantly underestimates the radiation pressure by a factor of ~100.

The time-scale over which and modality by which young stellar objects (YSOs) disperse their circumstellar discs dramatically influences the eventual formation and evolution of planetary systems. By means of extensive radiative transfer (RT) modelling, we have developed a new set of diagnostic diagrams in the infrared colour-colour plane (K-[24] vs. K-[8]), to aid with the classiffication of the evolutionary stage of YSOs from photometric observations. Our diagrams allow the differentiation of sources with un-evolved (primordial) discs from those evolving according to different clearing scenarios (e.g. homologous depletion vs. inside-out dispersal), as well as from sources that have already lost their disc. Classification of over 1500 sources in 15 nearby star-forming regions reveals that approximately 39% of the sources lie in the primordial disc region, whereas between 31% and 32% disperse from the inside-out and up to 22% of the sources have already lost their disc. Less than 2% of the objects in our sample lie in the homogeneous draining regime. Time-scales for the transition phase are estimated to be typically a few 10^5 years independent of stellar mass. Therefore, regardless of spectral type, we conclude that currently available infrared photometric surveys point to fast (of order 10% of the global disc lifetime) inside-out clearing as the preferred mode of disc dispersal.

We present a detailed study of star formation occurring in bound star--forming clouds under the influence of internal ionizing feedback from massive stars across a spectrum of cloud properties. We infer which objects are triggered by comparing our feedback simulations with control simulations in which no feedback was present. We find feedback always results in a lower star--formation efficiency and usually but not always results in a larger number of stars or clusters. Cluster mass functions are not strongly affected by feedback, but stellar mass functions are biased towards lower masses. Ionization also affects the geometrical distribution of stars in ways that are robust against projection effects, but may make the stellar associations more or less subclustered depending on the background cloud environment. We observe a prominent pillar in one simulation which is the remains of an accretion flow feeding the central ionizing cluster of its host cloud and suggest that this may be a general formation mechanism for pillars such as those observed in M16. We find that the association of stars with structures in the gas such as shells or pillars is a good but by no means foolproof indication that those stars have been triggered and we conclude overall that it is very difficult to deduce which objects have been induced to form and which formed spontaneously simply from observing the system at a single time.

We discuss three different ways in which stellar feedback may alter the outcome of star cluster formation: triggering or suppressing star formation, and redistributing the stellar population in space. We use detailed Smoothed Particle Hydrodynamics (SPH) simulations of HII regions in turbulent molecular clouds to show that all three of these may happen in the same system, making inferences about the effects of feedback problematic.

The vestiges of planet formation have been observed in debris disks harboring young and massive gaseous giants. The process of giant planet formation is terminated by the dissipation of gas in the protoplanetary disk. The gas-rich disk around HD142527 features a small inner disk, a large gap from \sim10 to \sim140AU, and a massive outer disk extending out to \sim300AU. The gap could have been carved-out by a giant planet. We have imaged the outer regions of this gap using the adaptive-optics camera NICI on Gemini South. Our images reveal that the disk is dynamically perturbed. The outer boundary of the roughly elliptical gap appears to be composed of several segments of spiral arms. The stellar position is offset by 0.17+-0.02" from the centroid of the cavity, consistent with earlier imaging at coarser resolutions. These transient morphological features are expected in the context of disk evolution in the presence of a perturbing body located inside the cavity. We perform hydro-dynamical simulations of the dynamical clearing of a gap in a disk. A 10Mjup body in a circular orbit at r = 90AU, perturbs the whole disks, even after thousands of orbits. By then the model disk has an eccentric and irregular cavity, flanked by tightly wound spiral arms, but it is still evolving far from steady state. A particular transient configuration that is a qualitative match to HD142527 is seen at 1.7Myr.

Aims: Grain growth has been suggested as one possible explanation for the diminished dust optical depths in the inner regions of protoplanetary "transition" disks. In this work, we directly test this hypothesis in the context of current models of grain growth and transport. Methods: A set of dust evolution models with different disk shapes, masses, turbulence parameters, and drift efficiencies is combined with radiative transfer calculations in order to derive theoretical spectral energy distributions (SEDs) and images. Results: We find that grain growth and transport effects can indeed produce dips in the infrared SED, as typically found in observations of transition disks. Our models achieve the necessary reduction of mass in small dust by producing larger grains, yet not large enough to be fragmenting efficiently. However, this population of large grains is still detectable at millimeter wavelengths. Even if perfect sticking is assumed and radial drift is neglected, a large population of dust grains is left behind because the time scales on which they are swept up by the larger grains are too long. This mechanism thus fails to reproduce the large emission cavities observed in recent millimeter-wave interferometric images of accreting transition disks.

We present nebular phase optical imaging and spectroscopy and near/mid-IR imaging of the Type II SN 2006bc. Observations reveal the central wavelength of the symmetric H$\alpha$ line profile to be red-shifted with respect to the host galaxy H$\alpha$ emission by day 325. Such an phenomenon has been argued to result from an asymmetric explosion in the iron-peak elements resulting in a larger mass of $^{56}$Ni and higher excitation of hydrogen on the far side of the SN explosion. We also observe a gradual blue-shifting of this H$\alpha$ peak which is indicative of dust formation in the ejecta. Although showing a normal peak brightness, V $\sim$ -17.2, for a core-collapse SN, 2006bc fades by $\sim$6 mag during the first 400 days suggesting either a relatively low $^{56}$Ni yield, an increase in extinction due to new dust, or both. A short duration flattening of the light curve is observed from day 416 to day 541 suggesting an optical light echo. Based on the narrow time window of this echo, we discuss implications on the location and geometry of the reflecting ISM. With our radiative transfer models, we find an upper limit of 2 x 10$^{-3}$ M$_{\odot}$ of dust around SN 2006bc. In the event that all of this dust were formed during the SN explosion, this quantity of dust is still several orders of magnitude lower than that needed to explain the large quantities of dust observed in the early universe.

{Abridged version for ArXiv}. We provide direct constraints on the origin of the [Ne II] emission in 15 young stars using high-spatial and spectral resolution observations with VISIR at the VLT that allow us to study the kinematics of the emitting gas. In addition we compare the [Ne II] line with optical forbidden lines observed for three stars with UVES. The [Ne II] line was detected in 7 stars, among them the first confirmed detection of [Ne II] in a Herbig Be star, V892 Tau. In four cases, the large blueshifted lines indicate an origin in a jet. In two stars, the small shifts and asymmetric profiles indicate an origin in a photo-evaporative wind. CoKu Tau 1, seen close to edge-on, shows a spatially unresolved line centered at the stellar rest velocity, although cross-dispersion centroids move within 10 AU from one side of the star to the other as a function of wavelength. The line profile is symmetric with wings extending up to about +-80 km/s. The origin of the [Ne II] line could either be due to the bipolar jet or to the disk. For the stars with VLT-UVES observations, in several cases, the optical forbidden line profiles and shifts are very similar to the profile of the [Ne II] line, suggesting that the lines are emitted in the same region. A general trend observed with VISIR is a lower line flux when compared with the fluxes obtained with Spitzer. We found no correlation between the line full-width at half maximum and the line peak velocity. The [Ne II] line remains undetected in a large part of the sample, an indication that the emission detected with Spitzer in those stars is likely extended.

We present an SPH parameter study of the dynamical effect of photoionization from O--type stars on star--forming clouds of a range of masses and sizes during the time window before supernovae explode. Our model clouds all have the same degree of turbulent support initially, the ratio of turbulent kinetic energy to gravitational potential energy being set to $E_{\rm kin}/|E_{\rm pot}|$=0.7. We allow the clouds to form stars and study the dynamical effects of the ionizing radiation from the massive stars or clusters born within them. We find that dense filamentary structures and accretion flows limit the quantities of gas that can be ionized, particularly in the higher density clusters. More importantly, the higher escape velocities in our more massive (10$^{6}$M$_{\odot}$) clouds prevent the HII regions from sweeping up and expelling significant quantities of gas, so that the most massive clouds are largely dynamically unaffected by ionizing feedback. However, feedback has a profound effect on the lower--density 10$^{4}$ and 10$^{5}$M$_{\odot}$ clouds in our study, creating vast evacuated bubbles and expelling tens of percent of the neutral gas in the 3Myr timescale before the first supernovae are expected to detonate, resulting in clouds highly porous to both photons and supernova ejecta.

We report optical and mid-infrared photometry of SN 1980K between 2004 and 2010, which show slow monotonic fading consistent with previous spectroscopic and photometric observations made 8 to 17 years after outburst. The slow rate-of-change over two decades suggests that this evolution may result from scattered and thermal light echoes off of extended circumstellar material. We present a semi- analytic dust radiative-transfer model that uses an empirically corrected effective optical depth to provide a fast and robust alternative to full Monte-Carlo radiative transfer modeling for homogenous dust at low to intermediate optical depths. We find that unresolved echoes from a thin circumstellar shell 14-15 lt-yr from the progenitor, and containing about 0.02 Msun of carbon-rich dust, can explain the broadband spectral and temporal evolution. The size, mass and dust composition are in good agreement with the contact discontinuity observed in scattered echoes around SN 1987A. The origin of slowly-changing high-velocity [O I] and Halpha lines is also considered. We propose an origin in shocked high-velocity metal-rich clumps of ejecta, rather than arising in the impact of ejecta on slowly-moving circumstellar material, as is the case with hot spots in SN 1987A.

Context: The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary disks and in the formation of larger bodies, such as planetesimals. Aims: We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the grain size distribution and which can readily be used for further modeling or for interpreting of observational data. Methods: We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmentation, and transport in viscously evolving gas disks. Results: We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the radial flux by solid material "rain out", which is crucial for triggering any gravity assisted formation of planetesimals. We show that fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of -1.5, while the outer regions at later times can become drift-dominated, yielding a dust surface density exponent of -0.75. Our results show that radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation due to the turbulent state of the disk.

Constraining the spatial and thermal structure of the gaseous component of circumstellar disks is crucial to understand star and planet formation. Models predict that the [Ne II] line at 12.81 {\mu}m detected in young stellar objects with Spitzer traces disk gas and its response to high energy radiation, but such [Ne II] emission may also originate in shocks within powerful outflows. To distinguish between these potential origins for mid-infrared [Ne II] emission and to constrain disk models, we observed 32 young stellar objects using the high resolution (R~30000) mid-infrared spectrograph VISIR at the VLT. We detected the 12.81 {\mu}m [Ne II] line in 12 objects, tripling the number of detections of this line in young stellar objects with high spatial and spectral resolution spectrographs. We obtain the following main results: a) In Class I objects the [Ne II] emission observed from Spitzer is mainly due to gas at a distance of more than 20-40 AU from the star, where neon is, most likely, ionized by shocks due to protostellar outflows. b) In transition and pre-transition disks, most of the emission is confined to the inner disk, within 20-40 AU from the central star. c) Detailed analysis of line profiles indicates that, in transition and pre-transition disks, the line is slightly blue-shifted (2-12 km s{^-1}) with respect to the stellar velocity, and the line width is directly correlated with the disk inclination, as expected if the emission is due to a disk wind. d) Models of EUV/X-ray irradiated disks reproduce well the observed relation between the line width and the disk inclination, but underestimate the blue-shift of the line.

We discuss a hydrodynamical model for the dispersal of protoplanetary discs around young, low mass (<1.5 M_sun) stars by photoevaporation from the central object's energetic radiation, which considers the far-ultraviolet as well as the X-ray component of the radiation field. We present analytical scaling relations and derive estimates for the total mass-loss rates, as well as discussing the existence of similarity solutions for flows from primordial discs and discs with inner holes. Furthermore, we perform numerical calculations, which span a wide range of parameter space and allow us to provide accurate scalings of the mass-loss rates with the physical parameters of the systems (X-ray and FUV luminosity, stellar mass, disc mass, disc temperature and inner hole radius). The model suggest that the X-ray component dominates the photoevaporative mass-loss rates from the inner disc. The mass-loss rates have values in the range from 10e-11 to 10e-7 M_sun/yr and scale linearly with X-ray luminosity, with only a weak dependence on the other parameters explored. However, in the case of high FUV to X-ray (L_FUV/L_X>100) luminosity ratios, the FUV constricts the X-ray flow and may dominate the mass-loss. Simulations of low mass discs with inner holes demonstrate a further stage of disc clearing, which we call `thermal sweeping'. This process occurs when the mid-plane pressure drops to sufficiently low values. At this stage a bound, warm, X-ray heated region becomes sufficiently large and unstable, such that the remaining disc material is cleared on approximately dynamical time-scales. This process significantly reduces the time taken to clear the outer regions of the disc, resulting in an expected transition disc population that will be dominated by accreting objects, as indicated by recent observations.

We have obtained an XMM-Newton observation and a broad-band spectrum from the ultraviolet to the near infrared with X-Shooter for one of the nearest M9 dwarfs, DENIS-P J1048-3956 (4pc). We integrate these data by a compilation of activity parameters for ultracool dwarfs from the literature with the aim to advance our understanding of these objects by comparing them to early-M type dwarf stars and the Sun. Our deep XMM-Newton observation has led to the first X-ray detection of DENIS-P J1048-3956 (log Lx = 25.1) as well as the first measurement of its V band brightness (V = 17.35mag). Flux-flux relations between X-ray and chromospheric activity indicators are here for the first time extended into the regime of the ultracool dwarfs. The approximate agreement of DENIS-P J1048-3956 and other ultracool dwarfs with flux-flux relations for early-M dwarfs suggests that the same heating mechanisms work in the atmospheres of ultracool dwarfs, albeit weaker as judged from their lower fluxes. The observed Balmer decrements of DENIS-P J1048-3956 are compatible with optically thick plasma in LTE at low, nearly photospheric temperature or optically thin LTE plasma at 20000K. Describing the decrements with CaseB recombination requires different emitting regions for Halpha and the higher Balmer lines. The high observed Halpha/Hbeta flux ratio is also poorly fitted by the optically thin models. We derive a similarly high value for the Halpha/Hbeta ratio of vB10 and LHS2065 and conclude that this may be a characteristic of ultracool dwarfs. We add DENIS-P J1048-3956 to the list of ultracool dwarfs detected in both the radio and the X-ray band. The Benz-Guedel relation between radio and X-ray luminosity of late-type stars is well-known to be violated by ultracool dwarfs. We speculate on the presence of two types of ultracool dwarfs with distinct radio and X-ray behavior.

We present BVRIJHK band photometry of 6 core-collapse supernovae, SNe 1999bw, 2002hh, 2003gd, 2004et, 2005cs, and 2006bc measured at late epochs (>2 yrs) based on Hubble Space Telescope (HST), Gemini north, and WIYN telescopes. We also show the JHK lightcurves of a supernova impostor SN 2008S up to day 575. Of our 43 HST observations in total, 36 observations are successful in detecting the light from the SNe alone and measuring magnitudes of all the targets. HST observations show a resolved scattered light echo around SN 2003gd at day 1520 and around SN 2002hh at day 1717. Our Gemini and WIYN observations detected SNe 2002hh and 2004et, as well. Combining our data with previously published data, we show VRIJHK-band lightcurves and estimate decline magnitude rates at each band in 4 different phases. Our prior work on these lightcurves and other data indicate that dust is forming in our targets from day ~300-400, supporting SN dust formation theory. In this paper we focus on other physical properties derived from the late time light curves. We estimate 56Ni masses for our targets (0.5-14 x 10^{-2} Msun) from the bolometric lightcurve of each for days ~150-300 using SN 1987A as a standard (7.5 x 10^{-2} Msun). The flattening or sometimes increasing fluxes in the late time light curves of SNe 2002hh, 2003gd, 2004et and 2006bc indicate the presence of light echos. We estimate the circumstellar hydrogen density of the material causing the light echo and find that SN 2002hh is surrounded by relatively dense materials (n(H) >400 cm^{-3}) and SNe 2003gd and 2004et have densities more typical of the interstellar medium (~1 cm^{-3}). The 56Ni mass appears well correlated with progenitor mass with a slope of 0.31 x 10^{-2}, supporting the previous work by Maeda et al. (2010), who focus on more massive Type II SNe. The dust mass does not appear to be correlated with progenitor mass.

We present a multi-wavelength study of three star forming regions, spanning the age range 1-14 Myr, located between the 30 Doradus complex and supernova SN1987A in the Large Magellanic Cloud (LMC). We reliably identify about 1000 pre-main sequence (PMS) star candidates actively undergoing mass accretion and estimate their stellar properties and mass accretion rate (Macc). Our measurements represent the largest Macc dataset of low-metallicity stars presented so far. As such, they offer a unique opportunity to study on a statistical basis the mass accretion process in the LMC and, more in general, the evolution of the mass accretion process around low-metallicity stars. We find that the typical \dot{M} of PMS stars in the LMC is higher than for galactic PMS stars of the same mass, independently of their age. Taking into account the caveats of isochronal age and \dot{M} estimates, the difference in Macc between the LMC and our Galaxy appears to be about an order of magnitude. We review the main mechanisms of disk dispersal and find indications that typically higher Macc are to be expected in low-metallicity environments. However, many issues of this scenario need to be clarified by future observations and modeling. We also find that, in the mass range 1-2 M_sun, the Macc of PMS stars in the LMC increases with stellar mass as \dot{M}acc\proptoM^b, with b \approx1, i.e. slower than the second power low found for galactic PMS stars in the same mass regime.

We present synthetic observations of star-forming interstellar medium structures obtained by hydrodynamic calculations of a turbulent box under the influence of an ionising radiation field. The morphological appearance of the pillar-like structures in optical emission lines is found to be very similar to observations of nearby star forming regions. We calculate line profiles as a function of position along the pillars for collisionally excited [OIII]5007, [NII]6584 and [SII]6717, which show typical FWHM of 2--4 km/s. Spatially resolved emission line diagnostic diagrams are also presented which show values in general agreement with observations of similar regions. The diagrams, however, also highlight significant spatial variations in the line ratios, including values that would be classically interpreted as shocked regions based on one-dimensional photoionisation calculations. These values tend to be instead the result of lines of sight intersecting which intersect for large portions of their lenghts the ionised-to-neutral transition regions in the gas. We caution therefore against a straightforward application of classical diagnostic diagrams and one--dimensional photoionisation calculations to spatially resolved observations of complex three-dimensional star forming regions.

We present an analysis of multi-epoch observations of the Type II-P supernova SN 2004et. New and archival optical spectra of SN 2004et are used to study the evolution of the Halpha and [O I] 6300A line profiles between days 259 and 646. Mid-infrared imaging was carried out between 2004 to 2010. We include Spitzer `warm' mission photometry at 3.6 and 4.5um obtained on days 1779, 1931 and 2151, along with ground-based and HST optical and near-infrared observations obtained between days 79 and 1803. Multi-wavelength light curves are presented, as well as optical-infrared spectral energy distributions (SEDs) for multiple epochs. Starting from about day 300, the optical light curves provide evidence for an increasing amount of circumstellar extinction attributable to newly formed dust, with the additional extinction reaching 0.8-1.5 magnitudes in the V-band by day 690. The overall SEDs were fitted with multiple blackbody components, in order to investigate the luminosity evolution of the supernova, and then with Monte Carlo radiative transfer models using smooth or clumpy dust distributions, in order to estimate how much new dust condensed in the ejecta. The luminosity evolution was consistent with the decay of 56Co in the ejecta up until about day 690, after which an additional emission source is required, in agreement with the findings of Kotak et al. (2009). Clumped dust density distributions consisting of 20% amorphous carbons and 80% silicates by mass were able to match the observed optical and infrared SEDs, with dust masses that increased from 8x10^{-5} Msun on day 300 to 1.5x10^{-3} Msun on day 690, still significantly lower than the values needed for core collapse supernovae to make a significant contribution to the dust enrichment of galaxies.

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.

We present a 3D photoionization model of the PN NGC 6302, one of the most complex objects of its kind. Our Mocassin model is composed of an extremely dense circumstellar disk and a large pair of diffuse bipolar lobes, a combination necessary to reproduce the observed spectrum. The masses of these components gives a total nebular mass of 4.7Mo. Discrepancies between our model fit and the observations are attributed to complex density inhomogeneities in the nebula. The potential to resolve such discrepancies with more complex models is confirmed by a range of models introducing small-scale structures. Compared to solar abundances He is enhanced by 50%, C is slightly subsolar, O is solar, and N is enhanced by a factor of 6. These imply a significant 3rd dredge-up coupled with hot-bottom burning CN-cycle conversion of dredged-up C to N. The central star is partly obscured by the edge-on circumstellar disk and its properties are not well constrained. Emission from a number of high-ionization `coronal' lines provides constraints on the form of the high-energy ionizing flux. Using a solar abundance stellar atmosphere we are unable to fit all of the observed line fluxes, but a substantially better fit was obtained using a 220,000K H-deficient stellar atmosphere with L*=14,300 Lo. The H-deficient nature of the central star suggests it has undergone a late thermal pulse, and fits to evolutionary tracks imply a central star mass of 0.73-0.82Mo. Timescales for these tracks suggest the object left the top of the AGB ~2100 years ago, in agreement with studies of the recent mass-loss event that formed the bipolar lobes. Based on the modelled nebular and central star masses we estimate the initial mass of the central star to be 5.5Mo, in agreement with that derived from evolutionary tracks. (Abstract truncated)

SN 2010jl was an extremely bright, Type IIn SNe which showed a significant IR excess no later than 90 days after explosion. We have obtained Spitzer 3.6 and 4.5 \mum and JHK observations of SN 2010jl \sim90 days post explosion. Little to no reddening in the host galaxy indicated that the circumstellar material lost from the progenitor must lie in a torus inclined out of the plane of the sky. The likely cause of the high mid-IR flux is the reprocessing of the initial flash of the SN by pre-existing circumstellar dust. Using a 3D Monte Carlo Radiative Transfer code, we have estimated that between 0.03-0.35 Msun of dust exists in a circumstellar torus around the SN located 6 \times 10 ^17 cm away from the SN and inclined between 60-80\cdot to the plane of the sky. On day 90, we are only seeing the illumination of approximately 5% of this torus, and expect to see an elevated IR flux from this material up until day \sim 450. It is likely this dust was created in an LBV-like mass loss event of more than 3 Msun, which is large but consistent with other LBV progenitors such as {\eta} Carinae.

We address the question of whether protoplanetary discs around low mass stars (e.g. M-dwarfs) may be longer lived than their solar-type counterparts. This question is particularly relevant to assess the planet-making potential of these lower mass discs. Given the uncertainties inherent to age-dating young stars, we propose an alternative approach that is to analyse the spatial distribution of disc-bearing low-mass stars and compare it to that of disc-bearing solar-type stars in the same cluster. A significant age difference between the two populations should be reflected in their average nearest neighbour distance (normalised to the number of sources), where the older population should appear more spread out. To this aim, we perform a Minimum Spanning Tree (MST) analysis on the spatial distribution of disc-bearing young stellar objects (YSOs) in six nearby low mass star forming regions. We find no evidence for significant age differences between the disc-bearing low-mass (later than M2) and 'solar-type' (earlier than M2) stars in these regions. We model our results by constructing and analysing synthetic fractal distributions that we evolve for typical values of the velocity dispersions. A comparison of simple models to our MST analysis suggests that the lifetime of discs around M-stars is similar to that of discs around solar-type stars. Furthermore, a model-independent spatial analysis of the observations robustly shows that any age differences between the two samples must be smaller than the average age difference between disc-bearing classical T-Tauri stars and disc-less Weak-Lined T-Tauri stars.

Photoevaporation driven by the central star is expected to be a ubiquitous and important mechanism to disperse the circumstellar dust and gas from which planets form. Here, we present a detailed study of the circumstellar disk surrounding the nearby star TW Hya and provide observational constraints to its photoevaporative wind. Our new high-resolution (R ~ 30,000) mid-infrared spectroscopy in the [Ne II] 12.81 {\mu}m line confirms that this gas diagnostic traces the unbound wind component within 10AU from the star. From the blueshift and asymmetry in the line profile, we estimate that most (>80%) of the [Ne II] emission arises from disk radii where the midplane is optically thick to the redshifted outflowing gas, meaning beyond the 1 or 4AU dust rim inferred from other observations. We re-analyze high-resolution (R ~ 48, 000) archival optical spectra searching for additional transitions that may trace the photoevaporative flow. Unlike the [Ne II] line, optical forbidden lines from OI, SII, and MgI are centered at the stellar velocity and have symmetric profiles. The only way these lines could trace the photoevaporative flow is if they arise from a disk region physically distinct from that traced by the [Ne II] line, specifically from within the optically thin dust gap. However, the small (~10 km/s) FWHM of these lines suggest that most of the emitting gas traced at optical wavelengths is bound to the system rather than unbound. We discuss the implications of our results for a planet-induced versus a photoevaporation-induced gap.

SN 2007it is a bright, Type IIP supernova which shows indications of both pre-existing and newly formed dust. The visible photometry shows a bright late-time luminosity, powered by the 0.09 M$_{\sun}$ of $^{56}$Ni present in the ejecta. There is also a sudden drop in optical brightness after day 339, and a corresponding brightening in the IR due to new dust forming in the ejecta. CO and SiO emission, generally thought to be precursors to dust formation, may have been detected in the mid-IR photometry of SN 2007it. The optical spectra show stronger than average [O I] emission lines and weaker than average [Ca II] lines, which may indicate a 16 - 27 M$_{\sun}$ progenitor, on the higher end of expected Type IIP masses. Multi-component [O I] lines are also seen in the optical spectra, most likely caused by an asymmetric blob or a torus of oxygen core material being ejected during the SN explosion. Interaction with circumstellar material prior to day 540 may have created a cool dense shell between the forward and reverse shocks where new dust is condensing. At late times there is also a flattening of the visible lightcurve as the ejecta luminosity fades and a surrounding light echo becomes visible. Radiative transfer models of SN 2007it SEDs indicate that up to 10$^{-4}$ M$_{\sun}$ of new dust has formed in the ejecta, which is consistent with the amount of dust formed in other core collapse supernovae.

Cool M dwarfs outnumber sun-like G stars by ten to one in the solar neighborhood. Due to their proximity, small size, and low mass, M-dwarf stars are becoming attractive targets for exoplanet searches via almost all current search methods. But what planetary systems can form around M dwarfs? Following up on the Cool Stars~16 Splinter Session "Planet Formation Around M Dwarfs", we summarize here our knowledge of protoplanetary disks around cool stars, how they disperse, what planetary systems might form and can be detected with current and future instruments.

We compare the three-dimensional gas temperature distributions obtained by a dedicated radiative transfer and photoionisation code, MOCASSIN, against those obtained by the recently-developed Smooth Particle Hydrodynamics (SPH) plus ionisation code iVINE for snapshots of an hydrodynamical simulation of a turbulent interstellar medium (ISM) irradiated by a nearby O star. Our tests demonstrate that the global ionisation properties of the region are correctly reproduced by iVINE, hence validating further application of this code to the study of feedback in star forming regions. However we highlight potentially important discrepancies in the detailed temperature distribution. In particular we show that in the case of highly inhomogenous density distributions the commonly employed on-the-spot (OTS) approximation yields unrealistically sharp shadow regions which can affect the dynamical evolution of the system. We implement a simple strategy to include the effects of the diffuse field in future calculations, which makes use of physically motivated temperature calibrations of the diffuse-field dominated regions and can be readily applied to similar codes. We find that while the global qualitative behaviour of the system is captured by simulations with the OTS approximation, the inclusion of the diffuse field in iVINE (called DiVINE) results in a stronger confinement of the cold gas, leading to denser and less coherent structures. This in turn leads to earlier triggering of star formation. We confirm that turbulence is being driven in simulations that include the diffuse field, but the efficiency is slightly lower than in simulations that use the OTS approximation.

Feedback from photoionisation may dominate on parsec scales in massive star-forming regions. Such feedback may inhibit or enhance the star formation efficiency and sustain or even drive turbulence in the parent molecular cloud. Photoionisation feedback may also provide a mechanism for the rapid expulsion of gas from young clusters' potentials, often invoked as the main cause of 'infant mortality'. There is currently no agreement, however, with regards to the efficiency of this process and how environment may affect the direction (positive or negative) in which it proceeds. The study of the photoionisation process as part of hydrodynamical simulations is key to understanding these issues, however, due to the computational demand of the problem, crude approximations for the radiation transfer are often employed. We will briefly review some of the most commonly used approximations and discuss their major drawbacks. We will then present the results of detailed tests carried out using the detailed photoionisation code MOCASSIN and the SPH+ionisation code iVINE code, aimed at understanding the error introduced by the simplified photoionisation algorithms. This is particularly relevant as a number of new codes have recently been developed along those lines. We will finally propose a new approach that should allow to efficiently and self-consistently treat the photoionisation problem for complex radiation and density fields.

We present a study of the variation of spatial structure of stellar populations within dwarf galaxies as a function of the population age. We use deep Hubble Space Telescope/Advanced Camera for Surveys imaging of nearby dwarf galaxies in order to resolve individual stars and create composite colour-magnitude diagrams (CMDs) for each galaxy. Using the obtained CMDs, we select Blue Helium Burning stars (BHeBs), which can be unambiguously age-dated by comparing the absolute magnitude of individual stars with stellar isochrones. Additionally, we select a very young (<10 Myr) population of OB stars for a subset of the galaxies based on the tip of the young main-sequence. By selecting stars in different age ranges we can then study how the spatial distribution of these stars evolves with time. We find, in agreement with previous studies, that stars are born within galaxies with a high degree of substructure which is made up of a continuous distribution of clusters, groups and associations from parsec to hundreds of parsec scales. These structures disperse on timescales of tens to hundreds of Myr, which we quantify using the two-point correlation function and the Q-parameter developed by Cartwright & Whitworth (2004). On galactic scales, we can place lower limits on the time it takes to remove the original structure (i.e., structure survives for at least this long), tevo, which varies between ~100~Myr (NGC~2366) and ~350 Myr (DDO~165). This is similar to what we have found previously for the SMC (~80~Myr) and the LMC (~175 Myr). We do not find any strong correlations between tevo and the luminosity of the host galaxy.

We have performed hydrodynamic and radiative transfer calculations of a photoevaporating disc around a Herbig Ae/Be star to determine the evolution and observational impact of dust entrained in the wind. We find that the wind selectively entrains grains of different sizes at different radii resulting in a dust population that varies spatially and increases with height above the disc at radii > 10 AU. This variable grain population results in a 'wingnut' morphology to the dust density distribution. We calculate images of this dust distribution at NIR wavelengths that also show a wingnut morphology at all wavelengths considered. We have also considered the contribution that crystalline dust grains will have in the wind and show that a photoevaporative wind can result in a significant crystallinity fraction at all radii, when the disc is edge-on. However, when the disc's photosphere is unobscured, a photoevaporative wind makes no contribution to the observable crystallinity fraction in the disc. Finally, we conclude that the analysis of extended emission around edge-on discs could provide a new and independent method of testing photoevaporation models.

Deep spectroscopy of the planetary nebula (PN) NGC\,6153 shows that its heavy element abundances derived from optical recombination lines (ORLs) are ten times higher than those derived from collisionally excited lines (CELs), and points to the existence of H-deficient inclusions embedded in the diffuse nebula. In this study, we have constructed chemically homogeneous and bi-abundance three-dimensional photoionization models, using the Monte Carlo photoionization code {\sc mocassin}. We attempt to reproduce the multi-waveband spectroscopic and imaging observations of NGC\,6153, and investigate the nature and origin of the postulated H-deficient inclusions, as well as their impacts on the empirical nebular analyses assuming a uniform chemical composition. Our results show that chemically homogeneous models yield small electron temperature fluctuations and fail to reproduce the strengths of ORLs from C, N, O and Ne ions. In contrast, bi-abundance models incorporating a small amount of metal-rich inclusions ($\sim 1.3$ per cent of the total nebular mass) are able to match all the observations within the measurement uncertainties. The metal-rich clumps, cooled down to a very low temperature ($\sim 800$~K) by ionic infrared fine-structure lines, dominate the emission of heavy element ORLs, but contribute almost nil to the emission of most CELs. We find that the abundances of C, N, O and Ne derived empirically from CELs, assuming a uniform chemical composition, are about 30 per cent lower than the corresponding average values of the whole nebula, including the contribution from the H-deficient inclusions. Ironically, in the presence of H-deficient inclusions, the traditional standard analysis of the optical helium recombination lines, assuming a chemically homogeneous nebula, overestimates the helium abundance by 40 per cent.

We have undertaken SED modeling of discs around low mass T Tauri stars in order to delineate regions of the infrared two colour plane (K - [8] versus K - [24]) that correspond to discs in different evolutionary stages. This provides a ready tool for classifying the nature of star-disc systems based on infrared photometry. In particular we demonstrate the distinct loci followed by discs that undergo `uniform draining' (reduction in surface density by a spatially uniform factor) from those that clear from the inside out. We draw attention to the absence of objects on this `draining locus' in those star forming regions where the 24um sensitivity would permit their detection, as compared with the ~20 objects in these regions with colours suggestive of inner holes. We thus conclude that discs predominantly clear from the inside out. We also apply our classification of the infrared two colour plane to stars of spectral type M3-M5 in the IC 348 cluster and show that some of its members (dubbed `weak excess' sources by Muzerolle et al 2010) that had previously been assumed to be in a state of clearing are instead likely to be optically thick discs in which the dust is well settled towards the mid-plane. Nevertheless, there are many discs in a state of partial clearing in IC 348, with their abundance (relative to the total population of disc bearing stars) being around four times higher than for comparable stars in Taurus. However, the number of partially cleared discs relative to the total number of late type stars is similarly low in both regions (10 and 20 % respectively). We therefore conclude that IC 348 represents a more evolved version of the Taurus population (with more of its discs being highly settled or partially cleared) but that the timescale for clearing is similarly short (a few times 10^5 years) in both cases.

We present a calculation of the atomic and low-ionisation emission line spectra of photoevaporating protoplanetary discs. Line luminosities and profiles are obtained from detailed photoionisation calculations of the disc and wind structures surrounding young active solar-type stars. The disc and wind density and velocity fields were obtained from the recently developed radiation-hydrodynamic models of Owen et al., that include stellar X-ray and EUV irradiation of protoplanetary discs at various stages of clearing, from primordial sources to inner hole sources of various hole sizes. Our models compare favourably with currently available observations, lending support to an X-ray driven photoevaporation model for disc dispersal. In particular, we find that X-rays drive a warm, predominantly neutral flow where the OI 6300A line can be produced by neutral hydrogen collisional excitation. Our models can, for the first time, provide a very good match to both luminosities and profiles of the low-velocity component of the OI 6300A line and other forbidden lines observed by Hartigan et al., which covered a large sample of T-Tauri stars. We find that the OI 6300A and the NeII 12.8um lines are predominantly produced in the X-ray-driven wind and thus appear blue-shifted by a few km/s for some of the systems when observed at non-edge-on inclinations. We note however that blue-shifts are only produced under certain conditions: X-ray luminosity, spectral shape and inner hole size all affect the location of the emitting region and the physical conditions in the wind. We caution therefore that while a blueshifted line is a tell-tale sign of an outflow, the lack of a blueshift should not be necessarily interpreted as a lack of outflow.

SN 2007od exhibits characteristics that have rarely been seen in a Type IIP supernova (SN). Optical V band photometry reveals a very steep brightness decline between the plateau and nebular phases of ~4.5 mag, likely due to SN 2007od containing a low mass of 56Ni. The optical spectra show an evolution from normal Type IIP with broad Halpha emission, to a complex, four component Halpha emission profile exhibiting asymmetries caused by dust extinction after day 232. This is similar to the spectral evolution of the Type IIn SN 1998S, although no early-time narrow (~200 km s-1) Halpha component was present in SN 2007od. In both SNe, the intermediate-width Halpha emission components are thought to arise in the interaction between the ejecta and its circumstellar medium (CSM). SN 2007od also shows a mid-IR excess due to new dust. The evolution of the Halpha profile and the presence of the mid-IR excess provide strong evidence that SN 2007od formed new dust before day 232. Late-time observations reveal a flattening of the visible lightcurve. This flattening is a strong indication of the presence of a light echo, which likely accounts for much of the broad, underlying Halpha component seen at late-times. We believe the multi-peaked Halpha emission is consistent with the interaction of the ejecta with a circumstellar ring or torus (for the inner components at \pm1500 km s-1), and a single blob or cloud of circumstellar material out of the plane of the CSM ring (for the outer component at -5000 km s-1). The most probable location for the formation of new dust is in the cool dense shell created by the interaction between the expanding ejecta and its CSM. Monte Carlo radiative transfer modeling of the dust emission from SN 2007od implies that up to 4x 10-4Msun of new dust has formed. This is similar to the amounts of dust formed in other CCSNe such as SNe 1999em, 2004et, and 2006jc.

SN 2008S erupted in early 2008 in the grand design spiral galaxy NGC 6946. The progenitor was detected by Prieto et al. in Spitzer Space Telescope images taken over the four years prior to the explosion, but was not detected in deep optical images, from which they inferred a self-obscured object with a mass of about 10 Msun. We obtained Spitzer observations of SN 2008S five days after its discovery, as well as coordinated Gemini and Spitzer optical and infrared observations six months after its outburst. We have constructed radiative transfer dust models for the object before and after the outburst, using the same r^-2 density distribution of pre-existing amorphous carbon grains for all epochs and taking light-travel time effects into account for the early post-outburst epoch. We rule out silicate grains as a significant component of the dust around SN 2008S. The inner radius of the dust shell moved outwards from its pre-outburst value of 85 AU to a post-outburst value of 1250 AU, attributable to grain vaporisation by the light flash from SN 2008S. Although this caused the circumstellar extinction to decrease from Av = 15 before the outburst to 0.8 after the outburst, we estimate that less than 2% of the overall circumstellar dust mass was destroyed. The total mass-loss rate from the progenitor star is estimated to have been (0.5-1.0)x10^-4 Msun yr^-1. The derived dust mass-loss rate of 5x10^-7 Msun yr^-1 implies a total dust injection into the ISM of up to 0.01 Msun over the suggested duration of the self-obscured phase. We consider the potential contribution of objects like SN 2008S to the dust enrichment of galaxies.

The formation of planets within a disc must operate within the time frame of disc dispersal, it is thus crucial to establish what is the dominant process that disperses the gaseous component of discs around young stars. Planet formation itself as well as photoevaporation by energetic radiation from the central young stellar object have been proposed as plausible dispersal mechanisms. [abridged]. In this paper we use the different metallicity dependance of X-ray photoevaporation and planet formation to discriminate between these two processes. We study the effects of metallicity, Z, on the dispersal timescale, t_phot, in the context of a photoevaporation model, by means of detailed thermal calculations of a disc in hydrostatic equilibrium irradiated by EUV and X-ray radiation from the central source. Our models show t_phot \propto Z^0.52 for a pure photoevaporation model. By means of analytical estimates we derive instead a much stronger negative power dependance on metallicity of the disc lifetime for a dispersal model based on planet formation. A census of disc fractions in lower metallicity regions should therefore be able to distinguish between the two models. A recent study by Yasui et al. in low metallicity clusters of the extreme outer Galaxy ([O/H] ~- 0.7dex and dust to gas ratio of ~0.001) provides preliminary observational evidence for shorter disc lifetimes at lower metallicities, in agreement with the predictions of a pure photoevaporation model. [abridged] We finally develop an analytical framework to study the effects of metallicity dependent photoevaporation on the formation of gas giants in the core accretion scenario. We show that accounting for this effect strengthens the conclusion that planet formation is favoured at higher metallicity. [abridged]

The [Ne II] fine-structure emission line at 12.8 micron has been detected in several young stellar objects (YSO) spectra. This line is thought to be produced by X-ray irradiation of the warm protoplanetary disk atmospheres, however the observational correlation between [Ne II] luminosities and measured X-ray luminosities shows a large scatter. Such spread limits the utility of this line as a probe of the gaseous phase of disks, as several authors have suggested pollution by outflows as a probable cause of the observed scatter. In this work we explore the possibility that the large variations in the observed [Ne II] luminosity may be caused instead by different star-disk parameters. In particular we study the effects that the hardness of the irradiating source and the structure (flaring) of the disk have on the luminosity and spectral profile of the [Ne II] 12.8 micron line. We find that varying these parameter can indeed cause up to an order of magnitude variation in the emission luminosities which may explain the scatter observed, although our models predict somewhat smaller luminosities than those recently reported by other authors who observed the line with the Spitzer Space Telescope. Our models also show that the hardness of the spectrum has only a limited (undetectable) effect on the line profiles, while changes in the flaring power of the disk significantly affect the size of the [Ne II] emission region and, as a consequence, its line profile. In particular we suggest that broad line profiles centred on the stellar radial velocity may be indicative of flat disks seen at large inclination angles.

(Abridged) We present the first radiation-hydrodynamic model of a protoplanetary disc irradiated with an X-EUV spectrum. In a model where the total ionizing luminosity is divided equally between X-ray and EUV luminosity, we find a photoevaporation rate of 1.4e-8 M_sun/yr, which is two orders of magnitude greater than the case of EUV photoevaporation alone. Thus it is clear that the X-rays are the dominant driving mechanism for photoevaporation. This can be understood inasmuch as X-rays are capable of penetrating much larger columns (~1e22 cm^-2) and can thus effect heating in denser regions and at larger radius than the EUV can. The radial extent of the launching region of the X-ray heated wind is 1-70AU compared with the pure EUV case where the launch region is concentrated around a few AU. When we couple our wind mass-loss rates with models for the disc's viscous evolution, we find that, as in the pure EUV case, there is a photoevaporative switch, such that an inner hole develops at ~ 1 AU at the point that the accretion rate in the disc drops below the wind mass loss rate. At this point, the remaining disc material is quickly removed in the final 15-20% of the disc's lifetime. This is consistent with the 1e5 yr transitional timescale estimated from observations of T-Tauri stars. We also caution that although our mass-loss rates are high compared to some accretion rates observed in young stars, our model has a rather large X-ray luminosity of 2e30 erg/s; further modeling is required in order to investigate the evolutionary implications of the large observed spread of X-ray luminosities in T-Tauri stars.

We present a theoretical investigation of the effect of multiple ionisation sources in HII regions on the total elemental abundances derived from the analysis of collisionally excited emission lines. We focus on empirical methods based on direct temperature measurements that are commonly employed in cases when the temperature of the nebular gas can be determined from the ratio of nebular to auroral lines of (e.g.) doubly ionised oxygen. We find that direct temperature methods that employ a two-temperature zone approach (DT2T methods) are very robust against the spatial distribution of sources. Errors smaller than 0.15 dex are estimated for regions where the metallicity is twice solar and errors below 0.05 dex for solar metallicities and below. The biases introduced by the spatial distribution of the ionisation sources are thus much smaller for DT2T methods than for strong line methods, previously investigated by Ercolano, Bastian & Stasinska. Our findings are in agreement with the recent study of HII regions in NGC 300 by Bresolin et al.

We re-examine the well-known discrepancy between ionic abundances determined via the analysis of recombination lines (RLs) and collisionally excited lines (CELs). We show that abundance variations can be mimicked in a {\it chemically homogeneous} medium by the presence of dense X-ray irradiated regions which present different ionisation and temperature structures from those of the more diffuse medium they are embedded in, which is predominantly ionised by extreme-ultraviolet radiation. The presence of X-ray ionised dense clumps or filaments also naturally explains the lower temperatures often measured from O {\sc ii} recombination lines and from the Balmer jump when compared to temperatures determined by CELs. We discuss the implications for abundances determined via the analysis of CELs and RLs and provide a simple analytical procedure to obtain upwards corrections for CEL-determined abundance. While we show that the abundance discrepancy factor (ADF) and the Balmer Jump temperature determined from observations of the Orion Nebula can simultaneously be reproduced by this model (implying upward corrections for CELs by a factor of 1.15), we find that the required X-ray fluxes exceed the known Orion's stellar and diffuse X-ray budget, if we assume that the clumps are located at the edge of the blister. We propose, however, that spatially resolved observations may be used to empirically test the model, and we outline how the framework developed in this letter may be applied in the future to objects with better constrained geometries (e.g. planetary nebulae).

X-ray luminosities of accreting T Tauri stars are observed to be systematically lower than those of non-accretors. There is as yet no widely accepted physical explanation for this effect, though it has been suggested that accretion somehow suppresses, disrupts or obscures coronal X-ray activity. Here, we suggest that the opposite might be the case: coronal X-rays modulate the accretion flow. We re-examine the X-ray luminosities of T Tauri stars in the Orion Nebula Cluster and find that not only are accreting stars systematically fainter, but that there is a correlation between mass accretion rate and stellar X-ray luminosity. We use the X-ray heated accretion disk models of Ercolano et al. to show that protoplanetary disk photoevaporative mass loss rates are strongly dependent on stellar X-ray luminosity and sufficiently high to be competitive with accretion rates. X-ray disk heating appears to offer a viable mechanism for modulating the gas accretion flow and could be at least partially responsible for the observed correlation between accretion rates and X-ray luminosities of T Tauri stars.

We present new models for the X-ray photoevaporation of circumstellar discs which suggest that the resulting mass loss (occurring mainly over the radial range 10-40 AU) may be the dominant dispersal mechanism for gas around low mass pre-main sequence stars, contrary to the conclusions of previous workers. Our models combine use of the MOCASSIN Monte Carlo radiative transfer code and a self-consistent solution of the hydrostatic structure of the irradiated disc. We estimate the resulting photoevaporation rates assuming sonic outflow at the surface where the gas temperature equals the local escape temperature and derive mass loss rates of ~10^{-9} M_sun/yr, typically a factor 2-10 times lower than the corresponding rates in our previous work (Ercolano et al., 2008) where we did not adjust the density structure of the irradiated disc. The somewhat lower rates, and the fact that mass loss is concentrated towards slightly smaller radii, result from the puffing up of the heated disc at a few AU which partially screens the disc at tens of AU. (.....) We highlight the fact that X-ray photoevaporation has two generic advantages for disc dispersal compared with photoevaporation by extreme ultraviolet (EUV) photons that are only modestly beyond the Lyman limit: the demonstrably large X-ray fluxes of young stars even after they have lost their discs and the fact that X-rays are effective at penetrating much larger columns of material close to the star (abridged).

We investigate the effect of line of sight temperature variations and noise on two commonly used methods to determine dust properties from dust continuum observations of dense cores. One method employs a direct fit to a modified blackbody SED; the other involves a comparison of flux ratios to an analytical prediction. Fitting fluxes near the SED peak produces inaccurate temperature and dust spectral index estimates due to the line of sight temperature (and density) variations. Longer wavelength fluxes in the Rayleigh-Jeans part of the spectrum (>~ 600 micron for typical cores) may more accurately recover the spectral index, but both methods are very sensitive to noise. The temperature estimate approaches the density weighted temperature, or "column temperature," of the source as short wavelength fluxes are excluded. An inverse temperature - spectral index correlation naturally results from SED fitting, due to the inaccurate isothermal assumption, as well as noise uncertainties. We show that above some "threshold" temperature, the temperatures estimated through the flux ratio method can be highly inaccurate. In general, observations with widely separated wavelengths, and including shorter wavelengths, result in higher threshold temperatures; such observations thus allow for more accurate temperature estimates of sources with temperatures less than the threshold temperature. When only three fluxes are available, a constrained fit, where the spectral index is fixed, produces less scatter in the temperature estimate when compared to the estimate from the flux ratio method.

We re-examine the recent suggestion of a high fraction of transition discs (i.e. those with a cleared inner hole) in M stars, motivated by the fact that we expect that, for M stars, even discs without inner holes should exhibit very weak excess shortward of around 10um. Our analysis of spectral energy distribution models suggest that this indeed means that M stars where a detectable excess begins at around 6um may be mis-classified as transition discs when in fact they have optically thick dust extending in to the dust sublimation radius. Consequently, we estimate that the transition disc fraction among M stars in the Coronet cluster is ~15 +/-10 % (rather than the recently claimed value of 50%). This revised figure would imply that the transition disc fraction is not after all markedly higher in later type stars. We suggest that for M stars, transition discs can only be readily identified if they have emission that is close to photospheric out to > 10um.

We investigate the potential of using ratios of fine structure and near-infrared forbidden line transitions of atomic carbon to diagnose protoplanetary disk extension. Using results from 2D photoionisation and radiative transfer modeling of a realistic protoplanetary disk structure irradiated by X-rays from a T Tauri star, we obtain theoretical emission maps from which we construct radial distributions of the strongest emission lines produced in the disk. We show that ratios of fine structure to near-infrared forbidden line emission of atomic carbon are especially promising to constrain the minimum size of gaseous protoplanetary disks. While theoretically viable, the method presents a number of observational difficulties that are also discussed here.

The projected distribution of stars in the Small Magellanic Cloud (SMC) from the Magellanic Clouds Photometric Survey is analysed. Stars of different ages are selected via criteria based on V magnitude and V-I colour, and the degree of `grouping' as a function of age is studied. We quantify the degree of structure using the two-point correlation function and a method based on the Minimum Spanning Tree and find that the overall structure of the SMC is evolving from a high degree of sub-structure at young ages (~10 Myr) to a smooth radial density profile. This transition is gradual and at ~75 Myr the distribution is statistically indistinguishable from the background SMC distribution. This time-scale corresponds to approximately the dynamical crossing time of stars in the SMC. The spatial positions of the star clusters in the SMC show a similar evolution of spatial distribution with age. Our analysis suggests that stars form with a high degree of (fractal) sub-structure, probably imprinted by the turbulent nature of the gas from which they form, which is erased by random motions in the galactic potential on a time-scale of a galactic crossing time.

We present an analysis of the spatial distribution of various stellar populations within the Large Magellanic Cloud. We combine mid-infrared selected young stellar objects, optically selected samples with mean ages between ~9 and ~1000 Myr, and existing stellar cluster catalogues to investigate how stellar structures form and evolve within the LMC. For the analysis we use Fractured Minimum Spanning Trees, the statistical Q parameter, and the two-point correlation function. Restricting our analysis to young massive (OB) stars we confirm our results obtained for M33, namely that the luminosity function of the groups is well described by a power-law with index -2, and that there is no characteristic length-scale of star-forming regions. We find that stars in the LMC are born with a large amount of substructure, consistent with a 2D fractal distribution with dimension ~1.8 and evolve towards a uniform distribution on a timescale of ~175 Myr. This is comparable to the crossing time of the galaxy and we suggest that stellar structure, regardless of spatial scale, will be eliminated in a crossing time. This may explain the smooth distribution of stars in massive/dense young clusters in the Galaxy, while other, less massive, clusters still display large amounts of structure at similar ages. By comparing the stellar and star cluster distributions and evolving timescales, we show that infant mortality of clusters (or 'popping clusters') have a negligible influence on galactic structure. Finally, we quantify the influence of the elongation, differential extinction, and contamination of a population on the measured Q value.

We present an analysis of the spatial distribution of various stellar populations within the Large and Small Magellanic Clouds. We use optically selected stellar samples with mean ages between ~9 and ~1000 Myr, and existing stellar cluster catalogues to investigate how stellar structures form and evolve within the LMC/SMC. We use two statistical techniques to study the evolution of structure within these galaxies, the $Q$-parameter and the two-point correlation function (TPCF). In both galaxies we find the stars are born with a high degree of substructure (i.e. are highly fractal) and that the stellar distribution approaches that of the 'background' population on timescales similar to the crossing times of the galaxy (~80/150 Myr for the SMC/LMC respectively). By comparing our observations to simple models of structural evolution we find that 'popping star clusters' do not significantly influence structural evolution in these galaxies. Instead we argue that general galactic dynamics are the main drivers, and that substructure will be erased in approximately the crossing time, regardless of spatial scale, from small clusters to whole galaxies. This can explain why many young Galactic clusters have high degrees of substructure, while others are smooth and centrally concentrated. We conclude with a general discussion on cluster 'infant mortality', in an attempt to clarify the time/spatial scales involved.

The observation by the Swift X-ray Telescope of the Fe K alpha_1, alpha_2 doublet during a large flare on the RS CVn binary system II Peg represents one of only two firm detections to date of photospheric Fe K alpha from a star other than our Sun. We present models of the Fe K alpha equivalent widths reported in the literature for the II Peg observations and show that they are most probably due to fluorescence following inner shell photoionisation of quasi-neutral Fe by the flare X-rays. Our models constrain the maximum height of flare the to 0.15 R_* assuming solar abundances for the photospheric material, and 0.1 R_* and 0.06 R_* assuming depleted photospheric abundances ([M/H]=-0.2 and [M/H]=-0.4, respectively). Accounting for an extended loop geometry has the effect of increasing the estimated flare heights by a factor of ~3. These predictions are consistent with those derived using results of flaring loop models, which are also used to estimate the flaring loop properties and energetics. From loop models we estimate a flare loop height of 0.13 R_*, plasma density of ~4 * 10^12 cm^-3 and emitting volume of ~6 * 10^30 cm^3. Our estimates for the flare dimensions and density allow us to estimate the conductive energy losses to E_cond <= 2 * 10^36 erg, consistent with upper limits previously obtained in the literature. Finally, we estimate the average energy output of this large flare to be ~10^33 erg sec^-1, or 1/10th of the stellar bolometric luminosity.

We present MOCASSIN 2D photoionisation and dust radiative transfer models of a prototypical T Tauri disk irradiated by X-rays from the young pre-main sequence star. The calculations demonstrate a layer of hot gas reaching temperatures of ~10^6 K at small radii and ~10^4 K at a distance of 1 AU. The gas temperatures decrease sharply with depth, but appear to be completely decoupled from dust temperatures down to a column depth of ~5*10^21 cm^-2. We predict that several fine-structure and forbidden lines of heavy elements, as well as recombination lines of hydrogen and helium, should be observable with current and future instrumentation, although optical lines may be smothered by the stellar spectrum. Predicted line luminosities are given for the the brightest collisionally excited lines (down to ~10^-8L_sun, and for recombination transitions from several levels of HI and HeI. The mass loss rate due to X-ray photoevaporation estimated from our models is of the order of 10^-8 M_sun yr^-1, implying a dispersal timescale of a few Myr for a disk of mass 0.027 M_sun, which is the mass of the disk structure model we employed. We discuss the limitations of our model and highlight the need for further calculations that should include the simultaneous solution of the 2D radiative transfer problem and the 1D hydrostatic equilibrium in the polar direction.

Monte Carlo calculations of the O Kalpha line fluoresced by coronal X-rays and emitted just above the temperature minimum region of the solar atmosphere have been employed to investigate the use of this feature as an abundance diagnostic. While quite weak, we estimate line equivalent widths in the range 0.02-0.2 AA, depending on the X-ray plasma temperature. The line remains essentially uncontaminated by blends for coronal temperatures T =< 3e6 K and should be quite observable, with a flux >~ 2 ph/s/arcmin^2. Model calculations for solar chemical mixtures with an O abundance adjusted up and down by a factor of 2 indicate 35-60% changes in O Kalpha line equivalent width, providing a potentially useful O abundance diagnostic. Sensitivity of equivalent width to differences between recently recommended chemical compositions with ``high'' and ``low'' complements of the CNO trio important for interpreting helioseismological observations is less accute, amounting to 20-26% at coronal temperatures T ~< 2e6 K. While still feasible for discriminating between these two mixtures, uncertainties in measured line equivalent widths and in the models used for interpretation would need to be significantly less than 20%. Provided a sensitive X-ray spectrometer with resolving power >= 1000 and suitably well-behaved instrumental profile can be built, X-ray fluorescence presents a viable means for resolving the solar ``oxygen crisis''.

We present evidence of Fe fluorescent emission in the Chandra HETGS spectrum of the single G-type giant HR 9024 during a large flare. In analogy to solar X-ray observations, we interpret the observed Fe K$\alpha$ line as being produced by illumination of the photosphere by ionizing coronal X-rays, in which case, for a given Fe photospheric abundance, its intensity depends on the height of the X-ray source. The HETGS observations, together with 3D Monte Carlo calculations to model the fluorescence emission, are used to obtain a direct geometric constraint on the scale height of the flaring coronal plasma. We compute the Fe fluorescent emission induced by the emission of a single flaring coronal loop which well reproduces the observed X-ray temporal and spectral properties according to a detailed hydrodynamic modeling. The predicted Fe fluorescent emission is in good agreement with the observed value within observational uncertainties, pointing to a scale height $\lesssim 0.3$\rstar. Comparison of the HR 9024 flare with that recently observed on II Peg by Swift indicates the latter is consistent with excitation by X-ray photoionization.

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

We present a new version of the fully 3D photoionization and dust radiative transfer code, MOCASSIN, that uses a Monte Carlo approach for the transfer of radiation. The X-ray enabled MOCASSIN allows a fully geometry independent description of low-density gaseous environments strongly photoionized by a radiation field extending from radio to gamma rays. The code has been thoroughly benchmarked against other established codes routinely used in the literature, using simple plane parallel models designed to test performance under standard conditions. We show the results of our benchmarking exercise and discuss applicability and limitations of the new code, which should be of guidance for future astrophysical studies with MOCASSIN.

X-ray spectra from stellar coronae are reprocessed by the underlying photosphere through scattering and photoionization events. While reprocessed X-ray spectra reaching a distant observer are at a flux level of only a few percent of that of the corona itself, characteristic lines formed by inner shell photoionization of some abundant elements can be significantly stronger. The emergent photospheric spectra are sensitive to the distance and location of the fluorescing radiation and can provide diagnostics of coronal geometry and abundance. Here we present Monte Carlo simulations of the photospheric Kalpha doublet arising from quasi-neutral Fe irradiated by a coronal X-ray source. Fluorescent line strengths have been computed as a function of the height of the radiation source, the temperature of the ionising X-ray spectrum, and the viewing angle. We also illustrate how the fluorescence efficiencies scale with the photospheric metallicity and the Fe abundance. Based on the results we make three comments: (1) fluorescent Fe lines seen from pre-main sequence stars mostly suggest flared disk geometries and/or super-solar disk Fe abundances; (2) the extreme ~1400 mA line observed from a flare on V1486 Ori can be explained entirely by X-ray fluorescence if the flare itself were partially eclipsed by the limb of the star; and (3) the fluorescent Fe line detected by Swift during a large flare on II Peg is consistent with X-ray excitation and does not require a collisional ionisation contribution. There is no convincing evidence supporting the energetically challenging explanation of electron impact excitation for observed stellar Fe Kalpha lines.

We present a three-dimensional photoionisation and dust radiative transfer model of NGC 6302, an extreme, high-excitation planetary nebula. We use the 3D photoionisation code Mocassin} to model the emission from the gas and dust. We have produced a good fit to the optical emission-line spectrum, from which we derived a density distribution for the nebula. A fit to the infrared coronal lines places strong constraints on the properties of the unseen ionising source. We find the best fit comes from using a 220,000 K hydrogen-deficient central star model atmosphere, indicating that the central star of this PN may have undergone a late thermal pulse. We have also fitted the overall shape of the ISO spectrum of NGC 6302 using a dust model with a shallow power-law size distribution and grains up to 1.0 micron in size. To obtain a good fit to the infrared SED the dust must be sufficiently recessed within the circumstellar disk to prevent large amounts of hot dust at short wavelengths, a region where the ISO spectrum is particularly lacking. These and other discoveries are helping to unveil many properties of this extreme object and trace it's evolutionary history.

Monte Carlo calculations of the Ne K$\alpha$ line fluoresced by coronal x-rays and emitted near the temperature minimum region of the solar atmosphere have been employed to investigate the use of this feature to measure directly the solar photospheric Ne abundance. Though very weak, comparison with spectral line databases indicates that at plasma temperatures typical of the quiet Sun and cool active regions ($\leq 2\times 10^6$ K) the line is isolated and unblended. A canonical solar chemical composition yields an equivalent width of $\sim 6$ m\AA (0.3 eV) when observed at heliocentric angles $\sim 0$. For a 1 arcmin field of view, photon fluxes at Earth are of order 0.2 ph s$^{-1}$ for the quiet Sun, rendering the Ne K$\alpha$ fluorescent line a quite feasible means for determining the solar photospheric Ne content.

Star-formation within galaxies appears on multiple scales, from spiral structure, to OB associations, to individual star clusters, and often sub-structure within these clusters. This multitude of scales calls for objective methods to find and classify star-forming regions, regardless of spatial size. To this end, we present an analysis of star-forming groups in the local group spiral galaxy M33, based on a new implementation of the Minimum Spanning Tree (MST) method. Unlike previous studies which limited themselves to a single spatial scale, we study star-forming structures from the effective resolution limit (~20pc) to kpc scales. We find evidence for a continuum of star-forming group sizes, from pc to kpc scales. We do not find a characteristic scale for OB associations, unlike that found in previous studies, and we suggest that the appearance of such a scale was caused by spatial resolution and selection effects. The luminosity function of the groups is found to be well represented by a power-law with an index, -2, similar to that found for clusters and GMCs. Additionally, the groups follow a similar mass-radius relation as GMCs. The size distribution of the groups is best described by a log-normal distribution and we show that within a hierarchical distribution, if a scale is selected to find structure, the resulting size distribution will have a log-normal distribution. We find an abrupt drop of the number of groups outside a galactic radius of ~4kpc, suggesting a change in the structure of the star-forming ISM, possibly reflected in the lack of GMCs beyond this radius. (abridged)

We present a new fast algorithm which allows the simulation of ionising radiation emitted from point sources to be included in high-resolution three-dimensional smoothed particle hydrodynamics simulations of star cluster formation. We employ a Str\"omgren volume technique in which we use the densities of particles near the line-of-sight between the source and a given target particle to locate the ionisation front in the direction of the target. Along with one--dimensional tests, we present fully three--dimensional comparisons of our code with the three--dimensional Monte-Carlo radiative transfer code, MOCASSIN, and show that we achieve good agreement, even in the case of highly complex density fields.

Spatially resolved studies of star forming regions show that the assumption of spherical geometry is not realistic in most cases, with a major complication posed by the gas being ionised by multiple non-centrally located stars or star clusters. We try to isolate the effects of multiple non-centrally located stars on the temperature and ionisation structure of HII regions, via the construction of 3D photoionisation models using the 3D Monte Carlo photoionisation code MOCASSIN. We find that the true temperature fluctuations due to the stellar distribution (as opposed to the large-scale temperature gradients due to other gas properties) are small in all cases and not a significant cause of error in metallicity studies. Strong emission lines from HII regions are often used to study the metallicity of star-forming regions. We compare integrated emission line spectra from our models and quantify any systematic errors caused by the simplifying assumption of a single, central location for all ionising sources. We find that the dependence of the metallicity indicators on the ionisation parameter causes a clear bias, due to the fact that models with a fully distributed configuration of stars always display lower ionisation parameters than their fully concentrated counterparts. The errors found imply that the geometrical distribution of ionisation sources may partly account for the large scatter in metallicities derived using model-calibrated empirical methods.

Aims: We investigate temperature fluctuations in H II regions in terms of a two-phase model, which assumes that the nebular gas consists of a hot and a cold phase. Methods: We derive general formulae for T([O III), the [O III] forbidden line temperature, and T(H I), the hydrogen Balmer jump temperature, in terms of the temperatures of the hot and cold phases, T_h and T_c. Results: For large temperature differences, the values of t2 required to account for the observed difference between T([O III]) and T(H I) are much lower than those deduced using the classical formulae that assume random and small amplitude temperature fluctuations. One should therefore be cautious when using a two-phase model to account for empirically derived $t^2$ values. We present a correction of a recent work by Giammanco & Beckman, who use a two-phase model to estimate the ionization rate of H II regions by cosmicrays. We show that a very small amount of cold gas is sufficient to account for t2 values typically inferred for H II regions.

We present a study of the effects of clumping on the emergent spectral energy distribution (SED) from dusty supernova (SN) shells illuminated by a diffuse radiation source distributed throughout the medium. (...) The fully 3D radiation transport problem is solved using a Monte Carlo code, MOCASSIN, and we present a set of models aimed at investigating the sensitivity of the SEDs to various clumping parameters. We find that, contrary to the predictions of analytical prescriptions, the combination of an optical and IR observational data set is sufficient to constrain dust masses even in the case where optically thick clumps are present. Using both smoothly varying and clumped grain density distributions, we obtain new estimates for the mass of dust condensed by the Type II SN 1987A by fitting the optical and infrared spectrophotometric data of Wooden et al. (1993) at two epochs (day 615 and day 775). (...) From our numerical models we derive dust masses for SN 1987A that are comparable to previous analytic clumped graphite grain mass estimates, and at least two orders of magnitude below the 0.1-0.3 Msol that have been predicted to condense as dust grains in primordial core collapse supernova ejecta. This low condensation efficiency for SN 1987A is in contrast to the case of SN 2003gd, for which a dust condensation efficiency as large as 0.12 has recently been estimated. (Abridged)

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.

We present late-time optical and mid-infrared observations of the Type-II supernova 2003gd in NGC 628. Mid-infrared excesses consistent with cooling dust in the ejecta are observed 499-678 days after outburst, and are accompanied by increasing optical extinction and growing asymmetries in the emission-line profiles. Radiative-transfer models show that up to 0.02 solar masses of dust has formed within the ejecta, beginning as early as 250 days after outburst. These observations show that dust formation in supernova ejecta can be efficient and that massive-star supernovae can be major dust producers throughout the history of the Universe.

Emission line spectra from HII regions are often used to study properties of the gas in star-forming regions, as well as temperatures and luminosities of the ionising sources. Empirical diagnostics for the interpretation of observational data must often be calibrated with the aid of photoionisation models. Most studies so far have been carried out by assuming spherical or plane-parallel geometries, with major limitations on allowed gas and dust density distributions and with the spatial distribution of multiple, non-centrally-located ionising sources not being accounted for. We present the first results of our theoretical study of geometric effects, via the construction of a number of 3D photoionisation models using the MOCASSIN code for a variety of spatial configurations and ionisation sources. We compare integrated emission line spectra from such configurations and show evidence of systematic errors caused by the simplifying assumption of a single, central location for all ionising sources.

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.

The study of photoionised gas in planetary nebulae (PNe) has played a major role achieving, over the years, a better understanding of a number of physical processes pertinent to a broader range of fields than just PNe studies, ranging from atomic physics to stellar evolution. Whilst empirical techniques are routinely employed for the analysis of the emission line spectra of such objects, the accurate interpretation of the observational data often requires the solution of the radiative transfer (RT) problem in the nebula, via the application of a photoionisation code.The development of the codes has been driven by the observational constraints available, but also compromised by the available computer power. Modern codes are faster and more flexible, with the ultimate goal being the achievement of a description of the observations relying on the smallest number of parameters possible. In this light, recent developments have been focused on the inclusion of new atomic data, the inclusion of a realistic treatment for dust grains mixed in the ionised and photon dominated regions (PDRs) and the expansion of some codes to PDRs with the inclusion of chemical reaction networks. Furthermore the last few years have seen the development of fully 3D photoionisation codes based on the Monte Carlo method. A review of the photoionisation codes currently in use is given here, with emphasis on recent developments, including the expansion of the models to the 3D domain, the identification of new observational constraints and how these can be used to extract useful information from realistic models.[Abridged]

Planetary nebulae (PNe) and their central stars (CSs) are ideal tools to test evolutionary theory: photospheric properties of their exciting stars give stringent constraints for theoretical predictions of stellar evolution. The nebular abundances display the star's photosphere at the time of the nebula's ejection which allows to look back into the history of stellar evolution - but, more importantly, they even provide a possibility to investigate on the chemical evolution of our Galaxy because most of the nuclear processed material goes back into the interstellar medium via PNe. The recent developments in observation techniques and a new three-dimensional photoionization code MOCASSIN enable us to analyze PNe properties precisely by the construction of consistent models of PNe and CSs. In addition to PNe imaging and spectroscopy, detailed information about the velocity field within the PNe is a pre-requisite to employ de-projection techniques in modeling the physical structureof the PNe.

The study of photoionised gas in planetary nebulae (PNe) has played a major role in the achievement, over the years, of a better understanding of a number of physical processes, pertinent to a broader range of fields than that of PNe studies, spanning from atomic physics to stellar evolution theories. Whilst empirical techniques are routinely employed for the analysis of the emission line spectra of these objects, the accurate interpretation of the observational data often requires the solution of a set of coupled equations, via the application of a photoionisation/plasma code. A number of large-scale codes have been developed since the late sixties, using various analytical or statistical techniques for the transfer of continuum radiation, mainly under the assumption of spherical symmetry and a few in 3D. These codes have been proved to be powerful and in many cases essential tools, but a clear idea of the underlying physical processes and assumptions is necessary in order to avoid reaching misleading conclusions. A brief review of the field of photoionisation today is given here, with emphasis on the recent developments, including the expansion of the models to the 3D domain. Attention is given to the identification of new available observational constraints and how these can used to extract useful information from realistic models. (abridged)

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.

Dust emission from the Type II supernova SN 2002hh in NGC 6946 has been detected at mid-infrared wavelengths by the Spitzer Space Telescope from 590 to 758 days after outburst and confirmed by higher angular resolution Gemini-N mid-IR observations. The day-600 5.8-24 micron emission can be fit by a 290-K blackbody having a luminosity of 1.6E+7 L_Sol. The minimum emitting radius of 1.1E+17 cm is too large for the emitting dust to have been formed in the supernova ejecta. Using radiative transfer models and realistic dust grain parameters, fits to the observed flux distribution could be obtained with an optically thick dust shell having a mass of 0.10-0.15 M_Sol, corresponding to a total dust+gas mass in excess of 10 M_Sol, suggesting a massive M supergiant or luminous blue variable precursor to this self-obscured object.

Deep optical spectra of the high excitation planetary nebula NGC 1501 and its W04 central star are presented. A recombination line abundance analysis of the central star's emission-line spectrum yields He:C:O mass fractions of 0.36:0.48:0.16, similar to those of PG1159 stars. A detailed empirical analysis of the nebular collisionally excited line (CEL) and optical recombination line (ORL) spectrum is presented, together with fully three-dimensional photoionisation modelling of the nebula. We found very large ORL-CEL abundance discrepancy factors (ADFs) for O^2+ (32) and Ne^2+ (33). The mean value of 5100K for the T_e derived from HeI recombination lines ratios is 6000K lower than the value of 11100K implied by the [OIII] line ratio. A three-dimensional photoionisation model of NGC 1501 was constructed using the photoionisation code MOCASSIN, based on our new spectroscopic data and using the three-dimensional electron density distribution determined from long-slit echellograms of the nebula by Ragazzoni et al. (2001). The initial models showed higher degrees of ionisation of heavy elements than indicated by observations. We investigated the importance of the missing low-temperature dielectronic recombination rates for third-row elements and have estimated upper limits to their rate coefficients. Our model heavily under-predicts the optical recombination line emission. We conclude that the presence of a hydrogen-deficient, metal-rich component is necessary to explain the observed ORL spectrum of this object. The existence of such knots could also provide a softening of the radiation field, via the removal of ionising photons by absorption in the knots, thereby helping to alleviate the over-ionisation of the heavy elements in our models.

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

We have developed a three dimensional Monte Carlo photoionization code tailored for the study of Galactic H II regions and the percolation of ionizing photons in diffuse ionized gas. We describe the code, our calculation of photoionization, heating & cooling, and the approximations we have employed for the low density H II regions we wish to study. Our code gives results in agreement with the Lexington H II region benchmarks. We show an example of a 2D shadowed region and point out the very significant effect that diffuse radiation produced by recombinations of helium has on the temperature within the shadow.

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.

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

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.

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.