Abstract
The composition of gas and solids in protoplanetary discs sets the composition of planets that form out of them. Recent chemical models have shown that the composition of gas and dust in ...discs evolves on Myr time-scales, with volatile species disappearing from the gas phase. However, discs evolve due to gas accretion and radial drift of dust on time-scales similar to these chemical time-scales. Here, we present the first model coupling the chemical evolution in the disc mid-planes with the evolution of discs due to accretion and radial drift of dust. Our models show that transport will always overcome the depletion of CO2 from the gas phase, and can also overcome the depletion of CO and CH4 unless both transport is slow (viscous α ≲ 10−3) and the ionization rate is high (ζ ≈ 10−17). Including radial drift further enhances the abundances of volatile species because they are carried in on the surface of grains before evaporating left at their ice lines. Due to large differences in the abundances within 10 au for models with and without efficient radial drift, we argue that composition can be used to constrain models of planet formation via pebble accretion.
Herbig Ae/Be stars (HAeBes) span a key mass range that links low- and high-mass stars, and thus provide an ideal window from which to explore their formation. This paper presents Very Large ...Telescope/X-shooter spectra of 91 HAeBes, the largest spectroscopic study of HAeBe accretion to date. A homogeneous approach to determining stellar parameters is undertaken for the majority of the sample. Measurements of the ultraviolet are modelled within the context of magnetospheric accretion, allowing a direct determination of mass accretion rates. Multiple correlations are observed across the sample between accretion and stellar properties: the youngest and often most massive stars are the strongest accretors, and there is an almost 1:1 relationship between the accretion luminosity and stellar luminosity. Despite these overall trends of increased accretion rates in HAeBes when compared to classical T Tauri stars, we also find noticeable differences in correlations when considering the Herbig Ae and Herbig Be subsets. This, combined with the difficulty in applying a magnetospheric accretion model to some of the Herbig Be stars, could suggest that another form of accretion may be occurring within Herbig Be mass range.
Abstract
In this paper, we model the chemical evolution of a 0.25 M⊙ protoplanetary disc surrounding a 1 M⊙ star that undergoes fragmentation due to self-gravity. We use smoothed particle ...hydrodynamics including a radiative transfer scheme, along with a time-dependent chemical evolution code to follow the composition of the disc and resulting fragments over approximately 4000 yr. Initially, four quasi-stable fragments are formed, of which two are eventually disrupted by tidal torques in the disc. From the results of our chemical modelling, we identify species that are abundant in the fragments (e.g. H2O, H2S, HNO, N2, NH3, OCS, SO), species that are abundant in the spiral shocks within the disc (e.g. CO, CH4, CN, CS, H2CO) and species that are abundant in the circumfragmentary material (e.g. HCO+). Our models suggest that in some fragments it is plausible for grains to sediment to the core before releasing their volatiles into the planetary envelope, leading to changes in, e.g., the C/O ratio of the gas and ice components. We would therefore predict that the atmospheric composition of planets generated by gravitational instability should not necessarily follow the bulk chemical composition of the local disc material.
We present high-resolution (∼300 au) Atacama Large Millimeter/submillimeter Array observations of the massive young stellar object G11.92-0.61 MM 1. We resolve the immediate circumstellar environment ...of MM 1 in 1.3 mm continuum emission and CH3CN emission for the first time. The object divides into two main sources-MM 1a, which is the source of a bipolar molecular outflow, and MM 1b, located 0 57 (1920 au) to the southeast. The main component of MM 1a is an elongated continuum structure, perpendicular to the bipolar outflow, with a size of 0 141 × 0 050 (480 × 170 au). The gas kinematics toward MM 1a probed via CH3CN trace a variety of scales. The lower energy J = 12-11 K = 3 line traces extended, rotating gas within the outflow cavity, while the v8 = 1 line shows a clearly resolved Keplerian rotation signature. Analysis of the gas kinematics and dust emission shows that the total enclosed mass in MM 1a is 40 5 M (where between 2.2 and 5.8 M is attributed to the disk), while MM 1b is <0.6 M . The extreme mass ratio and orbital properties of MM 1a and MM 1b suggest that MM 1b is one of the first observed examples of the formation of a binary star via disk fragmentation around a massive young (proto)star.
We present high-resolution millimeter continuum imaging of the disk surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disk system. The system has ...extended mm emission on which are superposed three prominent annular gaps at radii ∼13, 39, and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including the hot Jupiter, the system contains four gas giant planets at an age of only 2 Myr. Two of the new planets are similarly located to those inferred in the famous HL Tau protoplanetary disk; in CI Tau, additional observational data enables a more complete analysis of the system properties than was possible for HL Tau. Our dust and gas dynamical modeling satisfies every available observational constraint and points to the most massive ensemble of exoplanets ever detected at this age, with its four planets spanning a factor 1000 in orbital radius. Our results show that the association between hot Jupiters and gas giants on wider orbits, observed in older stars, is apparently in place at an early evolutionary stage.
Context.
Chemical modelling serves two purposes in dynamical models: accounting for the effect of microphysics on the dynamics and providing observable signatures. Ideally, the former must be done as ...part of the hydrodynamic simulation but this comes with a prohibitive computational cost that leads to many simplifications being used in practice.
Aims.
We aim to produce a statistical emulator that replicates a full chemical model capable of solving the temperature and abundances of a gas through time. This emulator should suffer only a minor loss of accuracy when compared to a full chemical solver and would have a fraction of the computational cost allowing it to be included in a dynamical model.
Methods.
The gas-grain chemical code UCLCHEM was updated to include heating and cooling processes, and a large dataset of model outputs from possible starting conditions was produced. A neural network was then trained to map directly from inputs to outputs.
Results.
Chemulator replicates the outputs of UCLCHEM with an overall mean squared error (MSE) of 1.7 × 10
−4
for a single time step of 1000 yr, and it is shown to be stable over 1000 iterations with an MSE of 3 × 10
−3
on the log-scaled temperature after one timzze step and 6 × 10
−3
after 1000 time steps. Chemulator was found to be approximately 50 000 times faster than the time-dependent model it emulates but can introduce a significant error to some models.
Abstract
The European FP7 project DIANA has performed a coherent analysis of a large set of observational data of protoplanetary disks by means of thermo-chemical disk models. The collected data ...include extinction-corrected stellar UV and X-ray input spectra (as seen by the disk), photometric fluxes, low and high resolution spectra, interferometric data, emission line fluxes, line velocity profiles and line maps, which probe the dust, polycyclic aromatic hydrocarbons (PAHs) and the gas in these objects. We define and apply a standardized modeling procedure to fit these data by state-of-the-art modeling codes (
ProDiMo
,
MCFOST
,
MCMax
), solving continuum and line radiative transfer (RT), disk chemistry, and the heating and cooling balance for both the gas and the dust. 3D diagnostic RT tools (e.g., FLiTs) are eventually used to predict all available observations from the same disk model, the DIANA-standard model. Our aim is to determine the physical parameters of the disks, such as total gas and dust masses, the dust properties, the disk shape, and the chemical structure in these disks. We allow for up to two radial disk zones to obtain our best-fitting models that have about 20 free parameters. This approach is novel and unique in its completeness and level of consistency. It allows us to break some of the degeneracies arising from pure Spectral Energy Distribution (SED) modeling. In this paper, we present the results from pure SED fitting for 27 objects and from the all inclusive DIANA-standard models for 14 objects. Our analysis shows a number of Herbig Ae and T Tauri stars with very cold and massive outer disks which are situated at least partly in the shadow of a tall and gas-rich inner disk. The disk masses derived are often in excess to previously published values, since these disks are partially optically thick even at millimeter wavelength and so cold that they emit less than in the Rayleigh–Jeans limit. We fit most infrared to millimeter emission line fluxes within a factor better than 3, simultaneously with SED, PAH features and radial brightness profiles extracted from images at various wavelengths. However, some line fluxes may deviate by a larger factor, and sometimes we find puzzling data which the models cannot reproduce. Some of these issues are probably caused by foreground cloud absorption or object variability. Our data collection, the fitted physical disk parameters as well as the full model output are available to the community through an online database (
http://www.univie.ac.at/diana
).
There have been several recent detections of candidate Keplerian discs around massive young protostars. Given the relatively large disc-to-star mass ratios in these systems, and their young ages, it ...is worth investigating their propensity to becoming self-gravitating. To this end, we compute self-consistent, semi-analytic models of putative self-gravitating discs for five candidate disc systems. Our aim is not to fit exactly the observations, but to demonstrate that the expected dust continuum emission from marginally unstable self-gravitating discs can be quite weak, due to high optical depth at the mid-plane even at millimetre wavelengths. In the best cases, the models produce ‘observable’ disc masses within a factor of <1.5 of those observed, with mid-plane dust temperatures comparable to measured temperatures from molecular line emission. We find in two cases that a self-gravitating disc model compares well with observations. If these discs are self-gravitating, they satisfy the conditions for disc fragmentation in their outer regions. These systems may hence have as-yet-unresolved low-mass stellar companions, and are thus promising targets for future high angular resolution observations.
Abstract
To date, there is no explanation as to why disc-tracing CO first overtone (or ‘bandhead’) emission is not a ubiquitous feature in low- to medium-resolution spectra of massive young stellar ...objects (MYSOs), but instead is only detected towards approximately 25 per cent of their spectra. In this paper, we investigate the hypothesis that only certain mass accretion rates result in detectable bandhead emission in the near-infrared spectra of MYSOs. Using an analytic disc model combined with an LTE model of the CO emission, we find that high accretion rates (≳10−4 M⊙ yr−1) result in large dust sublimation radii, a larger contribution to the K-band continuum from hot dust at the dust sublimation radius, and therefore correspondingly lower CO emission with respect to the continuum. On the other hand, low accretion rates (≲10−6 M⊙ yr−1) result in smaller dust sublimation radii, a correspondingly smaller emitting area of CO, and thus also lower CO emission with respect to the continuum. In general, moderate accretion rates produce the most prominent, and therefore detectable, CO first overtone emission. We compare our findings to a recent near-infrared spectroscopic survey of MYSOs, finding results consistent with our hypothesis. We conclude that the detection rate of CO bandhead emission in the spectra of MYSOs could be the result of MYSOs exhibiting a range of mass accretion rates, perhaps due to the variable accretion suggested by recent multi-epoch observations of these objects.