Source author record

Barbara Ercolano

Barbara Ercolano appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

46works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

46 published item(s)

preprint2026arXiv

Can photoevaporation open gaps in protoplanetary discs?

We investigate whether photoevaporation alone can open and sustain gaps in protoplanetary discs by coupling the evolving disc structure with the photoevaporative flow in two dimensional radiation hydrodynamical simulations. Our results show that once a density depression forms, the local mass-loss rate decreases sharply, suppressing further gap deepening. Viscous inflow and radial mass transport along the disc surface act to partially refill the depleted region, preventing complete clearing. The resulting configuration is a persistent, partially depleted zone whose evolution is largely insensitive to the initial disc morphology. This behaviour challenges the standard paradigm that photoevaporation efficiently carves clean inner cavities and directly produces transition discs. However, the pressure maximum at the outer edge of the depression may still trap dust grains, giving rise to transition disc like observational signatures. We also present a first-order prescription to approximate this behaviour in one dimensional disc evolution models, suitable for use in planet formation and population synthesis studies. Although the prescription improves upon static mass-loss treatments, it remains approximate, underscoring the need for further multidimensional simulations and parameter-space exploration to derive robust recipes for global disc and planet population models.

preprint2022arXiv

Observations of PAHs in the atmospheres of discs and exoplanets

Polycyclic aromatic hydrocarbons (PAHs) play a key role in the chemical and hydrodynamical evolution of the atmospheres of exoplanets and planet-forming discs. If they can survive the planet formation process, PAHs are likely to be involved in pre-biotic chemical reactions eventually leading to more complex molecules such as amino acids and nucleotides, which form the basis for life as we know it. However, the abundance and specific role of PAHs in these environments is largely unknown due to limitations in sensitivity and range of wavelength of current and previous space-borne facilities. Upcoming infrared space spectroscopy missions, such as Twinkle and Ariel, present a unique opportunity to detect PAHs in the atmospheres of exoplanets and planet-forming discs. In this work we present synthetic observations based on conservative numerical modeling of typical planet-forming discs and a transiting hot Saturnian planet around solar type star. Our models show that Twinkle and Ariel might both be able to detect the 3.3 micron PAH feature within reasonable observing time in discs and transiting planets, assuming that PAHs are present with an abundance of at least one tenth of the interstellar medium value.

preprint2022arXiv

The importance of X-ray frequency in driving photoevaporative winds

Photoevaporative winds are a promising mechanism for dispersing protoplanetary discs, but so far theoretical models have been unable to agree on the relative roles that the X-ray, Extreme Ultraviolet or Far Ultraviolet play in driving the winds. This has been attributed to a variety of methodological differences between studies, including their approach to radiative transfer and thermal balance, the choice of irradiating spectrum employed, and the processes available to cool the gas. We use the \textsc{mocassin} radiative transfer code to simulate wind heating for a variety of spectra on a static density grid taken from simulations of an EUV-driven wind. We explore the impact of choosing a single representative X-ray frequency on their ability to drive a wind by measuring the maximum heated column as a function of photon energy. We demonstrate that for reasonable luminosities and spectra, the most effective energies are at a few $100~\mathrm{eV}$, firmly in the softer regions of the X-ray spectrum, while X-rays with energies $\sim1000~\mathrm{eV}$ interact too weakly with disc gas to provide sufficient heating to drive a wind. We develop a simple model to explain these findings. We argue that further increases in the cooling above our models - for example due to molecular rovibrational lines - may further restrict the heating to the softer energies but are unlikely to prevent X-ray heated winds from launching entirely; increasing the X-ray luminosity has the opposite effect. The various results of photoevaporative wind models should therefore be understood in terms of the choice of irradiating spectrum.

preprint2022arXiv

Twinkle -- a small satellite spectroscopy mission for the next phase of exoplanet science

With a focus on off-the-shelf components, Twinkle is the first in a series of cost competitive small satellites managed and financed by Blue Skies Space Ltd. The satellite is based on a high-heritage Airbus platform that will carry a 0.45 m telescope and a spectrometer which will provide simultaneous wavelength coverage from 0.5-4.5 $\rm{μm}$. The spacecraft prime is Airbus Stevenage while the telescope is being developed by Airbus Toulouse and the spectrometer by ABB Canada. Scheduled to begin scientific operations in 2025, Twinkle will sit in a thermally-stable, sun-synchronous, low-Earth orbit. The mission has a designed operation lifetime of at least seven years and, during the first three years of operation, will conduct two large-scale survey programmes: one focused on Solar System objects and the other dedicated to extrasolar targets. Here we present an overview of the architecture of the mission, refinements in the design approach, and some of the key science themes of the extrasolar survey.

preprint2021arXiv

Giant planet migration during the disc dispersal phase

Transition discs are expected to be a natural outcome of the interplay between photoevaporation (PE) and giant planet formation. Massive planets reduce the inflow of material from the outer to the inner disc, therefore triggering an earlier onset of disc dispersal due to PE through a process known as Planet-Induced PhotoEvaporation (PIPE). In this case, a cavity is formed as material inside the planetary orbit is removed by PE, leaving only the outer disc to drive the migration of the giant planet. We investigate the impact of PE on giant planet migration and focus specifically on the case of transition discs with an evacuated cavity inside the planet location. This is important for determining under what circumstances PE is efficient at halting the migration of giant planets, thus affecting the final orbital distribution of a population of planets. For this purpose, we use 2D FARGO simulations to model the migration of giant planets in a range of primordial and transition discs subject to PE. The results are then compared to the standard prescriptions used to calculate the migration tracks of planets in 1D planet population synthesis models. The FARGO simulations show that once the disc inside the planet location is depleted of gas, planet migration ceases. This contradicts the results obtained by the impulse approximation, which predicts the accelerated inward migration of planets in discs that have been cleared inside the planetary orbit. These results suggest that the impulse approximation may not be suitable for planets embedded in transition discs. A better approximation that could be used in 1D models would involve halting planet migration once the material inside the planetary orbit is depleted of gas and the surface density at the 3:2 mean motion resonance location in the outer disc reaches a threshold value of $0.01\,\mathrm{g\,cm^{-2}}$.

preprint2020arXiv

Solid accretion onto planetary cores in radiative disks

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.

preprint2020arXiv

The interpretation of protoplanetary disc wind diagnostic lines from X-ray photoevaporation and analytical MHD models

High resolution spectra of typical wind diagnostics ([OI] 6300 Å 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.

preprint2016arXiv

Blueshifted [OI] lines from protoplanetary discs: the smoking gun of X-ray photoevaporation

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.

preprint2016arXiv

Constraining the escape fraction of ionizing photons from HII regions within NGC 300: A concept paper

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.

preprint2015arXiv

Far-infrared signatures and inner hole sizes of protoplanetary discs undergoing inside-out dust dispersal

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.

preprint2015arXiv

Observations and three-dimensional ionization structure of the planetary nebula SuWt 2

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.

preprint2015arXiv

The link between disc dispersal by photoevaporation and the semi-major axis distribution of exoplanets

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.

preprint2015arXiv

The long-term evolution of photoevaporating transition discs with giant planets

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.

preprint2014arXiv

Primordial mass segregation in simulations of star formation?

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.

preprint2014arXiv

Protoplanetary disc evolution affected by star-disc interactions in young stellar clusters

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}$.

preprint2013arXiv

Effects of shocks in stellar atmosphere models on the emission line spectrum of surrounding Hii regions

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.

preprint2013arXiv

Observations of gas flows inside a protoplanetary gap

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.

preprint2013arXiv

Testing protoplanetary disc dispersal with radio emission

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.

preprint2013arXiv

The interplay between X-ray photoevaporation and planet formation

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.

preprint2013arXiv

Three dimensional geometries and the analysis of H II regions

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.

preprint2012arXiv

Dead, Undead and Zombie Zones in Protostellar Disks as a Function of Stellar Mass

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.

preprint2012arXiv

Radiative Transfer in Star Formation: Testing FLD and Hybrid Methods

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.

preprint2012arXiv

The lifetime of protoplanetary discs: Observations and Theory

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.

preprint2012arXiv

The spectral energy distribution of protoplanetary diss around Massive Young Stellar Objects

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.

preprint2012arXiv

Thirty years of SN 1980K: Evidence for light echoes

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.

preprint2012arXiv

Triggering, suppressing and redistributing star formation

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.

preprint2011arXiv

A 3D Photoionization Model of the Extreme Planetary Nebula NGC 6302

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)

preprint2011arXiv

Late-time Light Curves of Type II Supernovae: Physical Properties of SNe and Their Environment

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.

preprint2011arXiv

On the lifetime of discs around late type stars

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.

preprint2011arXiv

On the theory of disc photoevaporation

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.

preprint2011arXiv

Synthetic observations of simulated pillars of creation

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.

preprint2010arXiv

Ionisation Feedback in Star and Cluster Formation Simulations

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.

preprint2010arXiv

Ionisation feedback in star formation simulations: The role of diffuse fields

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.

preprint2010arXiv

The clearing of discs around late type T Tauri stars: constraints from the infrared two colour plane

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.

preprint2010arXiv

The imprint of photoevaporation on edge-on discs

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.

preprint2010arXiv

Theoretical spectra of photoevaporating protoplanetary discs: An atlas of atomic and low-ionisation emission lines

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.

preprint2009arXiv

Can X-rays provide a solution to the abundance discrepancy problem in photoionised nebulae?

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).

preprint2009arXiv

Metallicity, planet formation, and disc lifetimes

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]

preprint2009arXiv

Multiple ionisation sources in HII regions and their effect on derived nebular abundances

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.

preprint2009arXiv

The Effect of Line of Sight Temperature Variation and Noise on Dust Continuum Observations

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.

preprint2009arXiv

X-ray irradiated protoplanetary disk atmospheres II: Predictions from models in hydrostatic equilibrium

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).

preprint2009arXiv

X-ray Photoevaporation-starved T Tauri Accretion

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.

preprint2007arXiv

A new algorithm for modelling photoionising radiation in smoothed particle hydrodynamics

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ömgren 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.

preprint2006arXiv

Massive-Star Supernovae as Major Dust Factories

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.

preprint2005arXiv

On the Velocity Field and the 3D Structure of the Galactic Soccer Ball Abell 43

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.