Source author record

B. Ercolano

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

49works
5topics
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

49 published item(s)

preprint2024arXiv

Small but mighty: High-resolution spectroscopy of ultra-hot Jupiter atmospheres with compact telescopes. KELT-9 b's transmission spectrum with Wendelstein's FOCES Spectrograph

When observing transmission spectra produced by atmospheres of ultra-hot Jupiters, large telescopes are typically the instrument of choice due to the very weak signal of the planet's atmosphere. This study aims to alleviate the desire for large telescopes by illustrating that the same science is possible with smaller telescope classes. We use the cross-correlation technique to showcase the potential of the high-resolution spectrograph FOCES at Wendelstein Observatory and demonstrate its potential to resolve the atmosphere of the ultra-hot Jupiter, KELT-9 b. A performance comparison is conducted between FOCES and HARPS-N spectrographs, considering both single transit and combined observations over three nights. With FOCES, we have detected seven species in KELT-9 b's atmosphere: Ti II, Fe I, Fe II, Na I, Mg I, Na II, Cr II, Sc II. Although HARPS-N surpasses FOCES in performance, our results reveal that smaller telescope classes are capable of resolving ultra-hot Jupiter atmospheres. This broadens the scope of potential studies, allowing for investigations into phenomena like temporal variations in atmospheric signals and the atmospheric loss characteristics of these close-in planets.

preprint2022arXiv

Dust entrainment in photoevaporative winds: Synthetic observations of transition disks

X-ray- and extreme-ultraviolet- (XEUV-) driven photoevaporative winds acting on protoplanetary disks around young T-Tauri stars may strongly impact disk evolution, affecting both gas and dust distributions. We compute dust densities for the wind regions of XEUV-irradiated transition disks with gap sizes of 20 and 30 AU, and determine whether they can be observed at wavelengths $0.7 \lesssim λ_\mathrm{obs} [μ\mathrm{m}] \lesssim 1.8$ in scattered and polarised light with current instrumentation. For an XEUV-driven outflow around a $M_* = 0.7 \mathrm{M}_\odot$ T-Tauri star with $L_X = 2 \cdot 10^{30} \mathrm{erg/s}$, we find dust mass-loss rates $\dot{M}_\mathrm{dust} \lesssim 2.0 \cdot 10^{-3} \dot{M}_\mathrm{gas}$, and if we invoke vertical settling, the outflow is quite collimated. The synthesised images exhibit a distinct chimney-like structure. The relative intensity of these chimneys is low, but under optimal conditions, their detection may still be feasible with current instrumentation such as JWST NIRCam and SPHERE IRDIS.

preprint2022arXiv

Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks

The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks. This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future.

preprint2021arXiv

Dust entrainment in photoevaporative winds: Densities and imaging

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 constrain the dust densities in a typical XEUV-driven outflow, and determine whether these winds can be observed at $μ\mathrm{m}$-wavelengths in scattered and polarised light. For an XEUV-driven outflow around a $M_* = 0.7\,\mathrm{M}_\odot$ T-Tauri star with $L_X = 2 \cdot 10^{30}\,\mathrm{erg/s}$, we find a dust mass-loss rate $\dot{M}_\mathrm{dust} \lesssim 4.1 \cdot 10^{-11}\,\mathrm{M_\odot / yr}$ for an optimistic estimate of dust densities in the wind (compared to $\dot{M}_\mathrm{gas} \approx 3.7 \cdot 10^{-8}\,\mathrm{M_\odot / yr}$). The synthesised scattered-light images suggest a distinct chimney structure emerging at intensities $I/I_{\max} < 10^{-4.5}$ ($10^{-3.5}$) at $λ_\mathrm{obs} = 1.6$ ($0.4$) $μ\mathrm{m}$, while the features in the polarised-light images are even fainter. Observations synthesised from our model do not exhibit clear features for SPHERE IRDIS, but show a faint wind signature for JWST NIRCam under optimal conditions. In conclusion, unambiguous detections of photoevaporative XEUV winds launched from primordial disks are at least challenging with current instrumentation; this provides a possible explanation as to why disk winds are not routinely detected in scattered or polarised light.

preprint2020arXiv

Detectability of embedded protoplanets from hydrodynamical simulations

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.

preprint2020arXiv

Dust entrainment in photoevaporative winds: The impact of X-rays

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\,μ$m ($9\,μ$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 $μ$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.

preprint2020arXiv

Modelling thermochemical processes in protoplanetary disks I: numerical methods

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.

preprint2019arXiv

Demographics of disks around young very low-mass stars and brown dwarfs in Lupus

We present new 890 $μ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.

preprint2016arXiv

Photoevaporation and close encounters: how the environment around Cygnus OB2 affects the evolution of protoplanetary disks

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.

preprint2015arXiv

A multi-wavelength analysis for interferometric (sub-)mm observations of protoplanetary disks: radial constraints on the dust properties and the disk structure

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.

preprint2015arXiv

Controversial Age Spreads from the Main Sequence Turn-Off and Red Clump in Intermediate-Age Clusters in the LMC

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.

preprint2015arXiv

Early evolution of embedded clusters

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.

preprint2015arXiv

Magnetospheres of hot Jupiters: hydrodynamic models & ultraviolet absorption

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 ($τ\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.

preprint2015arXiv

The Pillars of Creation revisited with MUSE: gas kinematics and high-mass stellar feedback traced by optical spectroscopy

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.

preprint2014arXiv

Before the first supernova: combined effects of HII regions and winds on molecular clouds

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.

preprint2014arXiv

Observations and three-dimensional photoionization modelling of the Wolf-Rayet planetary nebula Abell 48

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.

preprint2014arXiv

On the gas content of transitional disks: a VLT/X-Shooter study of accretion and winds

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.

preprint2014arXiv

The timing and location of dust formation in the remnant of SN 1987A

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\,$μ$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.

preprint2014arXiv

Turbulence in Giant Molecular Clouds: The effect of photoionisation feedback

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.

preprint2013arXiv

Ionization--induced star formation V: Triggering in partially unbound clusters

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.

preprint2013arXiv

Massive stars in massive clusters IV: Disruption of clouds by momentum-driven winds

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.

preprint2013arXiv

Optical IFU observations of gas pillars surrounding the super star cluster NGC 3603

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.

preprint2013arXiv

The Mdot - M* relation of pre-main sequence stars: a consequence of X-ray driven disc evolution

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.

preprint2013arXiv

X-ray ionization rates in protoplanetary discs

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

preprint2012arXiv

A simple model for the evolution of the dust population in protoplanetary disks

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.

preprint2012arXiv

Can grain growth explain transition disks?

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.

preprint2012arXiv

Disc Clearing of Young Stellar Objects: Evidence for Fast Inside-out Dispersal

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.

preprint2012arXiv

Ionization--induced star formation IV: Triggering in bound clusters

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.

preprint2012arXiv

Massive stars in massive clusters II: Disruption of bound clusters by photoionization

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.

preprint2012arXiv

On the origin of [Ne II] emission in young stars: mid-infrared and optical observations with the Very Large Telescope

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

preprint2012arXiv

Optical and Infrared Analysis of Type II SN 2006BC

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$α$ line profile to be red-shifted with respect to the host galaxy H$α$ 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$α$ 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.

preprint2012arXiv

The dynamically disrupted gap in HD 142527

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.

preprint2011arXiv

Evidence for Pre-Existing Dust in the Bright Type IIn SN 2010jl

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 η Carinae.

preprint2011arXiv

High resolution spectroscopy of Ne II emission from young stellar objects

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

preprint2011arXiv

Photometric and Spectroscopic Evolution of the Type IIP SN 2007it to Day 944

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.

preprint2011arXiv

Photometric determination of the mass accretion rates of pre-main sequence stars. III. Results in the Large Magellanic Cloud

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.

preprint2011arXiv

Planet Formation Around M-dwarfs: From Young Disks to Planets

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.

preprint2011arXiv

The effects of dust on the optical and infrared evolution of SN 2004et

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.

preprint2011arXiv

The Photoevaporative Wind from the Disk of TW Hya

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

preprint2011arXiv

The ultracool dwarf DENIS-P J104814.7-395606. Chromospheres and coronae at the low-mass end of the main-sequence

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.

preprint2011arXiv

When Shape Matters: correcting the ICFs to derive the chemical abundances of bipolar and elliptical PNe

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.

preprint2010arXiv

SN 2007od: A Type IIP SN with Circumstellar Interaction

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.

preprint2010arXiv

The evolution of stellar structures in dwarf galaxies

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.

preprint2010arXiv

Three-dimensional chemically homogeneous and bi-abundance photoionization models of the "super-metal-rich" planetary nebula NGC 6153

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.

preprint2009arXiv

Effects of X-ray irradiation and disk flaring on the [NeII] 12.8 micron emission from young stellar objects

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.

preprint2009arXiv

Radiation-Hydrodynamic Models of X-Ray & EUV Photoevaporating Protoplanetary Discs

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

preprint2009arXiv

The destruction and survival of dust in the shell around SN 2008S

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.

preprint2008arXiv

Evolution of stellar structure in the Small Magellanic Cloud

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.

preprint2008arXiv

Gas Phase diagnostics of Protoplanetary disk extension

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.