Search

Theoretical models indicate that photoevaporative and magnetothermal winds play a crucial role in the evolution and dispersal of protoplanetary disks and affect the formation of planetary systems. However, it is still unclear what wind-driving mechanism is dominant or if both are at work, perhaps at different stages of disk evolution. Recent spatially resolved observations by Fang et al. (2023) of the [OI] 6300 Angstrom spectral line, a common disk wind tracer, in TW Hya revealed that about 80% of the emission is confined to the inner few au of the disk. In this work, we show that state-of-the-art X-ray driven photoevaporation models can reproduce the compact emission and the line profile of the [OI] 6300 Angstrom line. Furthermore, we show that the models also simultaneously reproduce the observed line luminosities and detailed spectral profiles of both the [OI] 6300 Angstrom and the [NeII] 12.8 micron lines. While MHD wind models can also reproduce the compact radial emission of the [OI] 6300 Angstrom line, they fail to match the observed spectral profile of the [OI] 6300 Angstrom line and underestimate the luminosity of the [NeII] 12.8 micron line by a factor of three. We conclude that, while we cannot exclude the presence of an MHD wind component, the bulk of the wind structure of TW Hya is predominantly shaped by a photoevaporative flow., Comment: 7 pages, 4 figures, accepted for publication in Astrophysical Journal Letters

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.

Winds in protoplanetary disks play an important role in their evolution and dispersal. However, what physical process is driving the winds is still unclear (i.e. magnetically vs thermally driven), and can only be understood by directly confronting theoretical models with observational data. We use hydrodynamic photoevaporative disk-wind models and post-process them with a thermo-chemical model to produce synthetic observables for the o-H$_2$ at 2.12 micron and [OI] at 0.63 micron spectral lines and directly compare the results to a sample of observations. Our photoevaporative disk-wind model is consistent with the observed signatures of the blueshifted narrow low-velocity component (NLVC), which is usually associated with slow disk winds, for both tracers. Only for one out of seven targets that show blueshifted NLVCs does the photoevaporative model fail to explain the observed line kinematics. Our results also indicate that interpreting spectral line profiles by simple methods, such as the thin-disk approximation, to determine the line emitting region can yield misleading conclusions. The photoevaporative disk-wind models are largely consistent with the studied observational data set, but it is not possible to clearly discriminate between different wind-driving mechanisms. Further improvements to the models, such as consistent modelling of the dynamics and chemistry and detailed modelling of individual targets would be beneficial. Furthermore, a direct comparison of magnetically driven disk-wind models to the observational data set is necessary in order to determine whether or not spatially unresolved observations of multiple wind tracers are sufficient to discriminate between theoretical models., Comment: 14 pages, 13 figures, accepted for publication in A&A

The understanding of the accretion process has a central role in the understanding of star and planet formation. We aim to test how accretion variability influences previous correlation analyses of the relation between X-ray activity and accretion rates, which is important for understanding the evolution of circumstellar disks and disk photoevaporation. We monitored accreting stars in the Orion Nebula Cluster from November 24, 2014, until February 17, 2019, for 42 epochs with the Wendelstein Wide Field Imager in the Sloan Digital Sky Survey u'g'r' filters on the 2 m Fraunhofer Telescope on Mount Wendelstein. Mass accretion rates were determined from the measured ultraviolet excess. The influence of the mass accretion rate variability on the relation between X-ray luminosities and mass accretion rates was analyzed statistically. We find a typical interquartile range of ~ 0.3 dex for the mass accretion rate variability on timescales from weeks to ~ 2 years. The variability has likely no significant influence on a correlation analysis of the X-ray luminosity and the mass accretion rate observed at different times when the sample size is large enough. The observed anticorrelation between the X-ray luminosity and the mass accretion rate predicted by models of photoevaporation-starved accretion is likely not due to a bias introduced by different observing times.

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., Comment: 32 pages, 7 figures, to appear in Protostars and Planets VII, eds: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura

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 \lambda_\mathrm{obs} [\mu\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., Comment: Accepted for publication in A&A. 20+10 pages, 17+12 figures (References updated, minor adjustments.)

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 $\mu\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 $\lambda_\mathrm{obs} = 1.6$ ($0.4$) $\mu\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., Comment: Accepted for publication in A&A. 15 pages, 11 figures. (Affiliations updated, one reference corrected.)

In many astrophysical applications, the cost of solving a chemical network represented by a system of ordinary differential equations (ODEs) grows significantly with the size of the network, and can often represent a significant computational bottleneck, particularly in coupled chemo-dynamical models. Although standard numerical techniques and complex solutions tailored to thermochemistry can somewhat reduce the cost, more recently, machine learning algorithms have begun to attack this challenge via data-driven dimensional reduction techniques. In this work, we present a new class of methods that take advantage of machine learning techniques to reduce complex data sets (autoencoders), the optimization of multi-parameter systems (standard backpropagation), and the robustness of well-established ODE solvers to to explicitly incorporate time-dependence. This new method allows us to find a compressed and simplified version of a large chemical network in a semi-automated fashion that can be solved with a standard ODE solver, while also enabling interpretability of the compressed, latent network. As a proof of concept, we tested the method on an astrophysically-relevant chemical network with 29 species and 224 reactions, obtaining a reduced but representative network with only 5 species and 12 reactions, and a x65 speed-up., Comment: accepted A&A, code available at https://bitbucket.org/tgrassi/latent_ode_paper/

The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msun) young (age<10 Myr) stars at UV wavelengths. Here we present the PENELLOPE Large Program that is being carried out at the ESO Very Large Telescope (VLT) to acquire, contemporaneous to HST, optical ESPRESSO/UVES high-resolution spectra to investigate the kinematics of the emitting gas, and UV-to-NIR X-Shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time, and the fully reduced spectra are made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of thirteen targets located in the Orion OB1 association and in the sigma-Orionis cluster, observed in Nov-Dec 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days of <3. The comparison of the fits to the continuum excess emission obtained with a slab model on the X-Shooter spectra and the HST/STIS spectra shows a shortcoming in the X-Shooter estimates of <10%, well within the assumed uncertainty. Its origin can be either a wrong UV extinction curve or due to the simplicity of this modelling, and will be investigated in the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key for a better understanding of the accretion/ejection mechanisms in young stars., Comment: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

We perform a comprehensive demographic study of the CO extent relative to dust of the disk population in the Lupus clouds, in order to find indications of dust evolution and possible correlations with other properties. We increase up to 42 the number of disks of the region with measured CO and dust sizes ($R_{\mathrm{CO}}$, $R_{\mathrm{dust}}$) from observations with the Atacama Large Millimeter/submillimeter Array (ALMA). The sizes are obtained from modeling the ${^{12}}$CO $J = 2-1$ line emission and continuum emission at $\sim 0.89$ mm with an empirical function (Nuker profile or Gaussian function). The CO emission is more extended than the dust continuum, with a $R_{68\%}^{\mathrm{CO}}$/$R_{68\%}^{\mathrm{dust}}$ median value of 2.5, for the entire population and for a sub-sample with high completeness. 6 disks, around $15\%$ of the Lupus disk population have a size ratio above 4. Based on thermo-chemical modeling, this value can only be explained if the disk has undergone grain growth and radial drift. These disks do not have unusual properties in terms of stellar mass ($M_{\star}$), disk mass ($M_{\mathrm{disk}}$), CO and dust sizes ($R_{\mathrm{CO}}$, $R_{\mathrm{dust}}$), and mass accretion. We search for correlations between the size ratio and $M_{\star}$, $M_{\mathrm{disk}}$, $R_{\mathrm{CO}}$ and $R_{\mathrm{dust}}$: only a weak monotonic anti-correlation with the $R_{\mathrm{dust}}$ is found. The lack of strong correlations is remarkable and suggests that the bulk of the population may be in a similar evolutionary stage, independent of the stellar and disk properties. These results should be further investigated, since the optical depth difference between CO and dust continuum may play a role in the inferred size ratios. Lastly, the CO emission for the majority of the disks is consistent with optically thick emission and an average CO temperature of around 30 K., Comment: Accepted for publication on A&A, 14 pages of main text with 5 figures, and 11 pages of appendices A, B, C, D, E and F with 13 figures

The evaporation of molecules from dust grains is crucial to understand some key aspects of the star- and the planet-formation processes. During the warm-up phase the presence of young protostellar objects induces molecules to evaporate from the dust surface into the gas phase, enhancing its chemical complexity. Similarly, in circumstellar disks, the position of the so-called snow-lines is determined by evaporation, with important consequences for the formation of planets. The amount of molecules that are desorbed depends on the interaction between the species and the grain surface, which is controlled by the binding energy. Recent theoretical and experimental works point towards a distribution of values for this parameter instead of the single value often employed in astrochemical models.We present here a new "multi-binding energy" framework, to assess the effects that a distribution of binding energies has on the amount of species bound to the grains. We find that the efficiency of the surface chemistry is significantly influenced by this process with crucial consequences on the theoretical estimates of the desorbed species., Comment: Accepted A&A

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., Comment: MNRAS accepted

We calculated hydrogen recombination line luminosities (H-$\alpha$, Paschen-$\beta$ and Brackett-$\gamma$) from three dimensional thermo-hydrodynamical simulations of forming planets from 1 to 10 Jupiter-masses. We explored various opacities to estimate the line emissions with extinction in each cases assuming boundary layer accretion. When realistic opacities are considered, only lines from planets $\ge$10 Jupiter-mass can be detected with current instrumentation, highlighting that from most planets one cannot expect detectable emission. This might explain the very low detection rate of H-$\alpha$ from forming planets from observations. While the line emission comes from both the forming planet and its circumplanetary disk, we found that only the disk component could be detected due to extinction. We examined the line variability as well, and found that it is higher for higher mass planets. Furthermore, we determine for the first time, the parametric relationship between the mass of the planet and the luminosity of the hydrogen recombination lines, as well as the equation between the accretion luminosity and hydrogen recombination line luminosities., Comment: ApJ accepted. 13 pages, 5 figures, 4 tables. This version: typo corrections

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

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., Comment: Accepted for publication on January 8, 2020 in MNRAS. 15 pages of main text with 14 figures, and 5 pages of appendices A and B with 4 figures

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

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

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

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., Comment: Accepted for publication in ApJS as part of the special issue on the Chandra Cygnus OB2 Legacy Project

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., Comment: 19 pages, 15 figures, accepted for publication in A&A

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

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

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., Comment: 17 pages, 18 figures, accepted for publication in A&A

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., Comment: 19 pages, 17 figures (many degraded to get below submission size limit), accepted by MNRAS

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., Comment: 25 pages, 21 figures, appendix 10 pages and 10 figures, MNRAS in press

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., Comment: accepted by MNRAS, 13 pages, 12 figures

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

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., Comment: Accepted on Astronomy & Astrophysics. Abstract shortened to fit arXiv constraints

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., Comment: 21 pages, 19 figures, accepted for publication in MNRAS

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., Comment: 12 pages, 7 figures, 12 tables, accepted by MNRAS

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., Comment: 9 pages, 6 Figures, 2 Tables, accepted by MNRAS

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., Comment: Accepted by MNRAS, 18 pages, 18 figures

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, Comment: 6 pages, 6 figures, accepted by MNRAS

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., Comment: 17 pages, 19 figures, accepted to MNRAS

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., Comment: 17 pages, 11 figures (mostly degraded to get under the submission size limit), accepted by MNRAS

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

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., Comment: 31 pages, 21 figures, 6 tables, accepted for publication in MNRAS

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., Comment: 15 pages, 12 figures, to appear in MNRAS

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., Comment: ApJL, accepted

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., Comment: 11 pages, 5 figures, accepted to A&A

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

{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., Comment: Accepted for publication in Astronomy & Astrophysics; revised version: corrected minor typos, corrected center values (col 3) for CoKuTau1 in Table 5

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., Comment: 18 pages, 17 figures (mostly degraded to get under the filesize limit), accepted by MNRAS

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., Comment: 12 pages, 10 figures, accepted to A&A

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

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., Comment: accepted for publication in Astronomy & Astrophysics