Source author record

Cornelis Dullemond

Cornelis Dullemond 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

11works
4topics
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

11 published item(s)

preprint2020arXiv

The Dynamic Proto-atmospheres around Low-Mass Planets with Eccentric Orbits

Protoplanets are able to accrete primordial atmospheres when embedded in the gaseous protoplanetary disk. The formation and structure of the proto-atmosphere are subject to the planet--disk environment and orbital effects. Especially, when planets are on eccentric orbits, their relative velocities to the gas can exceed the sound speed. The planets generate atmosphere-stripping bow shocks. We investigate the proto-atmospheres on low-mass planets with eccentric orbits with radiation-hydrodynamics simulations. A 2D radiative model of the proto-atmosphere is established with tabulated opacities for the gas and dust. The solutions reveal large-scale gas recycling inside a bow shock structure. The atmospheres on eccentric planets are typically three to four orders of magnitude less massive than those of planets with circular orbits. Overall, however, a supersonic environment is favorable for planets to keep an early stable atmosphere, rather than harmful, due to the steady gas supply through the recycling flow. We also quantitatively explore how such atmospheres are affected by the relative velocity of the planet to the gas, the planet mass, and the background gas density. Our time-dependent simulations track the orbital evolution of the proto-atmosphere with the planet--disk parameters changing throughout the orbit. Atmospheric properties show oscillatory patterns as the planet travels on an eccentric orbit, with a lag in phase. To sum up, low-mass eccentric planets can retain small proto-atmospheres despite the stripping effects of bow shocks. The atmospheres are always connected to and interacting with the disk gas. These findings provide important insights into the impacts of migration and scattering on planetary proto-atmospheres.

preprint2020arXiv

The efficiency of dust trapping in ringed proto-planetary discs

When imaged at high-resolution, many proto-planetary discs show gaps and rings in their dust sub-mm continuum emission profile. These structures are widely considered to originate from local maxima in the gas pressure profile. The properties of the underlying gas structures are however unknown. In this paper we present a method to measure the dust-gas coupling $α/St$ and the width of the gas pressure bumps affecting the dust distribution, applying high-precision techniques to extract the gas rotation curve from emission lines data-cubes. As a proof-of-concept, we then apply the method to two discs with prominent sub-structure, HD163296 and AS 209. We find that in all cases the gas structures are larger than in the dust, confirming that the rings are pressure traps. Although the grains are sufficiently decoupled from the gas to be radially concentrated, we find that the degree of coupling of the dust is relatively good ($α/St \sim 0.1$). We can therefore reject scenarios in which the disc turbulence is very low and the dust has grown significantly. If we further assume that the dust grain sizes are set by turbulent fragmentation, we find high values of the $α$ turbulent parameter ($α\sim 10^{-2}$). Alternatively, solutions with smaller turbulence are still compatible with our analysis if another process is limiting grain growth. For HD163296, recent measurements of the disc mass suggest that this is the case if the grain size is 1mm. Future constraints on the dust spectral indices will help to discriminate between the two alternatives.

preprint2019arXiv

Midplane temperature and outer edge of the protoplanetary disk around HD 163296

Knowledge of the midplane temperature of protoplanetary disks is one of the key ingredients in theories of dust growth and planet formation. However, direct measurement of this quantity is complicated, and often depends on the fitting of complex models to the data. In this paper we demonstrate a method to directly measure the midplane gas temperature from an optically thick molecular line, if the disk is moderately inclined. The only model assumption that enters is that the line is very optically thick, also in the midplane region where we wish to measure the temperature. Freeze-out of the molecule onto dust grains could thwart this. However, in regions that are expected to be warm enough to avoid freeze-out, this method should work. We apply the method to the CO 2-1 line channel maps of the disk around HD 163296. We find that the midplane temperature between 100 and 400 au drops only mildly from 25 K down to 18 K. While we see no direct evidence of the midplane being optically thin due to strong CO depletion by freeze-out, we cannot rule it out either. The fact that the inferred temperatures are close to the expected CO freeze-out temperature could be an indication of this. Incidently, for the disk around HD 163296 we also find dynamic evidence for a rather abrupt outer edge of the disk, suggestive of outside-in photoevaporation or truncation by an unseen companion.

preprint2015arXiv

Spiral arms in scattered light images of protoplanetary discs: Are they the signposts of planets?

One of the striking discoveries of protoplanetary disc research in recent years are the spiral arms seen in several transitional discs in polarised scattered light. An interesting interpretation of the observed spiral features is that they are density waves launched by one or more embedded (proto-)planets in the disc. In this paper we investigate whether planets can be held responsible for the excitation mechanism of the observed spirals. We use locally isothermal hydrodynamic simulations as well as analytic formulae to model the spiral waves launched by planets. Then H-band scattered light images are calculated using a 3D continuum radiative transfer code to study the effect of surface density and pressure scale height perturbation on the detectability of the spirals. We find that a relative change of about 3.5 in the surface density is required for the spirals to be detected with current telescopes in the near-infrared for sources at the distance of typical star-forming regions (140pc). This value is a factor of eight higher than what is seen in hydrodynamic simulations. We also find that a relative change of only 0.2 in pressure scale height is sufficient to create detectable signatures under the same conditions. Therefore, we suggest that the spiral arms observed to date in protoplanetary discs are the results of changes in the vertical structure of the disc (e.g. pressure scale height perturbation) instead of surface density perturbations.

preprint2015arXiv

The Impact of Dust Evolution and Photoevaporation on Disk Dispersal

Protoplanetary disks are dispersed by viscous evolution and photoevaporation in a few million years; in the interim small, sub-micron sized dust grains must grow and form planets. The time-varying abundance of small grains in an evolving disk directly affects gas heating by far-ultraviolet photons, while dust evolution affects photoevaporation by changing the disk opacity and resulting penetration of FUV photons in the disk. Photoevaporative flows, in turn, selectively carry small dust grains leaving the larger particles---which decouple from the gas---behind in the disk. We study these effects by investigating the evolution of a disk subject to viscosity, photoevaporation by EUV, FUV and X-rays, dust evolution, and radial drift using a 1-D multi-fluid approach (gas + different dust grain sizes) to solve for the evolving surface density distributions. The 1-D evolution is augmented by 1+1D models constructed at each epoch to obtain the instantaneous disk structure and determine photoevaporation rates. The implementation of a dust coagulation/fragmentation model results in a marginal decrease in disk lifetimes when compared to models with no dust evolution; the disk lifetime is thus found to be relatively insensitive to the evolving dust opacity. We find that photoevaporation can cause significant reductions in the gas/dust mass ratio in the planet-forming regions of the disk as it evolves, and may result in a corresponding increase in heavy element abundances relative to hydrogen. We discuss implications for theories of planetesimal formation and giant planet formation, including the formation of gas-poor giants. After gas disk dispersal, $\sim 3\times 10^{-4}$ \ms\ of mass in solids typically remain, comparable to the solids inventory of our solar system.

preprint2012arXiv

Optical spectroscopy of EX Lupi during quiescence and outburst: Infall, wind, and dynamics in the accretion flow

We explore the accretion mechanisms in EX Lupi, prototype of EXor variables, during its quiescence and outburst phases. We analyse high-resolution optical spectra taken before, during, and after its 2008 outburst. In quiescence and outburst, the star presents many permitted emission lines, including typical CTTS lines and numerous neutral and ionized metallic lines. During the outburst, the number of emission lines increases to over a thousand, with narrow plus broad component structure (NC+BC). The BC profile is highly variable on short timescales (24-72h). An active chromosphere can explain the metallic lines in quiescence and the outburst NC. The dynamics of the BC line profiles suggest an origin in a hot, dense, non-axisymmetric, and non-uniform accretion column that suffers velocity variations along the line-of-sight on timescales of days. Assuming Keplerian rotation, the emitting region would be located at ~0.1-0.2 AU, consistent with the inner disk rim, but the velocity profiles of the lines reveal a combination of rotation and infall. Line ratios of ions and neutrals can be reproduced with a temperature of T~6500 K for electron densities of a few times 10$^{12}$cm$^{-3}$ in the line-emitting region. The data confirm that the 2008 outburst was an episode of increased accretion, albeit much stronger than previous EX Lupi and typical EXors outbursts. The line profiles are consistent with the infall/rotation of a non-axisymmetric structure that could be produced by clumpy accretion during the outburst phase. A strong inner disk wind appears in the epochs of higher accretion. The rapid recovery of the system after the outburst and the similarity between the pre-outburst and post-outburst states suggest that the accretion channels are similar during the whole period, and only the accretion rate varies, providing a superb environment for studying the accretion processes.

preprint2011arXiv

Hot HCN around young massive stars at 0.1" resolution

Massive stars form deeply embedded in dense molecular gas, which they stir and heat up and ionize. During an early phase, the ionization is confined to hypercompact HII regions, and the stellar radiation is entirely absorbed by dust, giving rise to a hot molecular core. To investigate the innermost structure of such high-mass star-forming regions, we observed vibrationally excited HCN (via the direct $\ell$-type transition of v2=1, $ΔJ$=0, J=13, which lies 1400 K above ground) toward the massive hot molecular cores G10.47+0.03, SgrB2-N, and SgrB2-M with the Very Large Array (VLA) at 7 mm, reaching a resolution of about 1000 AU (0.1"). We detect the line both in emission and in absorption against HII regions. The latter allows to derive lower limits on the column densities of hot HCN, which are several times $10^{19}$ cm$^{-2}$. We see indication of expansion motions in G10.47+0.03 and detect velocity components in SgrB2-M at 50, 60, and 70 km/s relative to the Local Standard of Rest. The emission originates in regions of less than 0.1 pc diameter around the hypercompact HII regions G10.47+0.03 B1 and SgrB2-N K2, and reaches brightness temperatures of more than 200 K. Using the three-dimensional radiative transfer code RADMC-3D, we model the sources as dense dust cores heated by stars in the HII regions, and derive masses of hot (>300 K) molecular gas of more than 100 solar masses (for an HCN fractional abundance of 10$^{-5}$), challenging current simulations of massive star formation. Heating only by the stars in the HII regions is sufficient to produce such large quantities of hot molecular gas, provided that dust is optically thick to its own radiation, leading to high temperatures through diffusion of radiation.

preprint2011arXiv

The 2008 outburst of EX Lup - silicate crystals in motion

EX Lup is the prototype of the EXor class of eruptive young stars. These objects show optical outbursts which are thought to be related to runaway accretion onto the star. In a previous study we observed in-situ crystal formation in the disk of EX Lup during its latest outburst in 2008, making the object an ideal laboratory to investigate circumstellar crystal formation and transport. This outburst was monitored by a campaign of ground-based and Spitzer Space Telescope observations. Here we modeled the spectral energy distribution of EX Lup in the outburst from optical to millimeter wavelengths with a 2D radiative transfer code. Our results showed that the shape of the SED at optical wavelengths was more consistent with a single temperature blackbody than a temperature distribution. We also found that this single temperature component emitted 80-100 % of the total accretion luminosity. We concluded that a thermal instability, the most widely accepted model of EXor outbursts, was likely not the triggering mechanism of the 2008 outburst of EX Lup. Our mid-infrared Spitzer spectra revealed that the strength of all crystalline bands between 8 and 30 um increased right after the end of the outburst. Six months later, however, the crystallinity in the 10 um silicate feature complex decreased. Our modeling of the mid-infrared spectral evolution of EXLup showed that, although vertical mixing should be stronger during the outburst than in the quiescent phase, fast radial transport of crystals (e.g., by stellar/disk wind) was required to reproduce the observed mid-infrared spectra.

preprint2010arXiv

Fast and accurate frequency-dependent radiation transport for hydrodynamics simulations in massive star formation

Context: Radiative feedback plays a crucial role in the formation of massive stars. The implementation of a fast and accurate description of the proceeding thermodynamics in pre-stellar cores and evolving accretion disks is therefore a main effort in current hydrodynamics simulations. Aims: We introduce our newly implemented three-dimensional frequency dependent radiation transport algorithm for hydrodynamics simulations of spatial configurations with a dominant central source. Methods: The module combines the advantage of the speed of an approximate Flux Limited Diffusion (FLD) solver with the high accuracy of a frequency dependent first order ray-tracing routine. Results: We prove the viability of the scheme in a standard radiation benchmark test compared to a full frequency dependent Monte-Carlo based radiative transfer code. The setup includes a central star, a circumstellar flared disk, as well as an envelope. The test is performed for different optical depths. Considering the frequency dependence of the stellar irradiation, the temperature distributions can be described precisely in the optically thin, thick, and irradiated transition regions. Resulting radiative forces onto dust grains are reproduced with high accuracy. The achievable parallel speedup of the method imposes no restriction on further radiative (magneto-) hydrodynamics simulations. Conclusions: The proposed approximate radiation transport method enables frequency dependent radiation hydrodynamics studies of the evolution of pre-stellar cores and circumstellar accretion disks around an evolving massive star in a highly efficient and accurate manner.

preprint2010arXiv

The first stages of planet formation in binary systems: How far can dust coagulation proceed?

We examine the first phase of the core accretion model, namely the dust growth/fragmentation in binary systems. In our model, a gas and dust disk is present around the primary star and is perturbed by the secondary. We study the effects of a secondary with/without eccentricity on the dust population to determine what sizes the aggregates can reach and how that compares to the dust population in disks around single stars. We find that the secondary star has two effects on the dust population. 1.) The disk is truncated due to the presence of the secondary star and the maximum mass of the particles is decreased in the lowered gas densities. This effect is dominant in the outer disk. 2.) The perturbation of the secondary pumps up the eccentricity of the gas disk, which in turn increases the relative velocity between the dust and the gas. Therefore the maximum particle sizes are further decreased. The second effect of the secondary influences the entire disk. Coagulation is efficiently reduced even at the very inner parts of the disk. The average mass of the particles is reduced by four orders of magnitude (as a consequence, the stopping time is reduced by one order of magnitude) in disks around binary systems compared to dust in disks around single stars.

preprint2009arXiv

Do we really know the dust? Systematics and uncertainties of the mid-infrared spectral analysis methods

The spectral region around 10 micrometer, showing prominent dust emission bands, is commonly used to derive the chemical composition of protoplanetary dust. Different analysis methods have been proposed for this purpose, but so far, no comparative test has been performed to test their validity. We calculated model spectra of disk models with different geometries and central sources, using a 2D radiative transfer code. These spectra were then fitted in a blind test using four different spectral decomposition methods. We studied the effect of disk structure (flared vs. flat), inclination angle, size of the inner disk hole and stellar luminosity on the fitted chemical composition. Our results show that the derived dust compositions by all methods deviate systematically from the real chemical composition. Out of the four tested spectral decomposition methods, our new two-layer temperature distribution method, differs the least from the input dust composition and the results show the weakest systematic effects. The reason for the deviations of the results given by other methods lies in their simplifying assumptions (e.g. single average grain temperature or one component continuum). We also tested the influence of different noise levels on the results of the spectral decomposition methods. We find that, for ground-based observations (8-13 micrometer) the expected uncertainty in the value of the crystallinity is about 11% for a signal-to-noise ratio of 100, while for space-based observations (7-17 micrometer) the uncertainty is about 5%. On the basis of our results, we propose a recipe for the analysis and interpretation of dust spectroscopy data in the mid-infrared which should be especially valuable for analysing Spitzer data and ground-based infrared spectroscopy data in the 10 micrometer window.