BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

T. Birnstiel

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

ResearcherUnclaimed source record
astro-ph.EPastro-ph.SRastro-ph.GA

Catalog footprint

What is connected

31works
3topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

31 published item(s)

preprint2022arXiv

Dust entrainment in photoevaporative winds: Synthetic observations of transition disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Setting the Stage for Planet Formation: Measurements and Implications of the Fundamental Disk Properties

The field of planet formation is in an exciting era, where recent observations of disks around low- to intermediate-mass stars made with state of the art interferometers and high-contrast optical and IR facilities have revealed a diversity of substructures, some possibly planet-related. It is therefore important to understand the physical and chemical nature of the protoplanetary building blocks, as well as their spatial distribution, to better understand planet formation. Since PPVI, the field has seen tremendous improvements in observational capabilities, enabling both surveys of large samples of disks and high resolution imaging studies of a few bright disks. Improvements in data quality and sample size have, however, opened up many fundamental questions about properties such as the mass budget of disks, its spatial distribution, and its radial extent. Moreover, the vertical structure of disks has been studied in greater detail with spatially resolved observations, providing new insights on vertical layering and temperature stratification, yet also bringing rise to questions about other properties, such as material transport and viscosity. Each one of these properties - disk mass, surface density distribution, outer radius, vertical extent, temperature structure, and transport - is of fundamental interest as they collectively set the stage for disk evolution and corresponding planet formation theories. In this chapter, we will review our understanding of the fundamental properties of disks including the relevant observational techniques to probe their nature, modelling methods, and the respective caveats. Finally, we discuss the implications for theories of disk evolution and planet formation underlining what new questions have since arisen as our observational facilities have improved.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Late infall causing disk misalignment and dynamic structures in SU Aur

Gas-rich circumstellar disks are the cradles of planet formation. As such, their evolution will strongly influence the resulting planet population. In the ESO DESTINYS large program, we study these disks within the first 10 Myr of their development with near-infrared scattered light imaging. Here we present VLT/SPHERE polarimetric observations of the nearby class II system SU Aur in which we resolve the disk down to scales of ~7 au. In addition to the new SPHERE observations, we utilize VLT/NACO, HST/STIS and ALMA archival data. The new SPHERE data show the disk around SU Aur and extended dust structures in unprecedented detail. We resolve several dust tails connected to the Keplerian disk. By comparison with ALMA data, we show that these dust tails represent material falling onto the disk. The disk itself shows an intricate spiral structure and a shadow lane, cast by an inner, misaligned disk component. Our observations suggest that SU Aur is undergoing late infall of material, which can explain the observed disk structures. SU Aur is the clearest observational example of this mechanism at work and demonstrates that late accretion events can still occur in the class II phase, thereby significantly affecting the evolution of circumstellar disks. Constraining the frequency of such events with additional observations will help determine whether this process is responsible for the spin-orbit misalignment in evolved exoplanet systems.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Dust entrainment in photoevaporative winds: Densities and imaging

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

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

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

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Constraining turbulence mixing strength in transitional discs with planets using SPHERE and ALMA

We investigate the effect that the turbulent mixing strength parameter $α_{\rm{turb}}$ plays on near-infrared polarimetric and sub-millimetre interferometric imaging observations of transitional discs (TDs) with a gap carved by a planet. We generate synthetic observations of these objects with ALMA and VLT/SPHERE-ZIMPOL by combining hydrodynamical, dust evolution, radiative transfer and instrument models for values of $α_{\rm{turb}}=[10^{-4}, 10^{-3}, 10^{-2}]$. We find that, through a combination of effects on the viscosity of the gas, the turbulent mixing and dust evolution processes, $α_{\rm{turb}}$ strongly affects the morphology of the dust distribution that can be traced with these observations. We constrain the value of $α_{\rm{turb}}$ to be within an order of magnitude of $10^{-3}$ in TD sources that show cavities in sub-mm continuum images while featuring continuous distribution of dust or smaller cavities in NIR-polarimetric images.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Dust Evolution and the Formation of Planetesimals

The solid content of circumstellar disks is inherited from the interstellar medium: dust particles of at most a micrometer in size. Protoplanetary disks are the environment where these dust grains need to grow at least 13 orders of magnitude in size. Our understanding of this growth process is far from complete, with different physics seemingly posing obstacles to this growth at various stages. Yet, the ubiquity of planets in our galaxy suggests that planet formation is a robust mechanism. This chapter focuses on the earliest stages of planet formation, the growth of small dust grains towards the gravitationally bound "planetesimals", the building blocks of planets. We will introduce some of the key physics involved in the growth processes and discuss how they are expected to shape the global behavior of the solid content of disks. We will consider possible pathways towards the formation of larger bodies and conclude by reviewing some of the recent observational advances in the field.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Investigating dust trapping in transition disks with millimeter-wave polarization

Spatially resolved polarized (sub-)mm emission has been observed for example in the protoplanetary disk around HL Tau. Magnetically aligned grains are commonly interpreted as the source of polarization. However, self-scattering by large dust grains with a high enough albedo is another polarization mechanism, becoming a compelling method independent of the spectral index to constrain the dust grain size in protoplanetary disks. We study the dust polarization at mm wavelength in the dust trapping scenario proposed for transition disks, when a giant planet opens a gap in the disk. We investigate the characteristic polarization patterns and their dependence on disk inclination, dust size evolution, planet position, and observing wavelength. We combine two-dimensional hydrodynamical simulations of planet-disk interactions with self-consistent dust growth models. These size-dependent dust density distributions are used for follow-up three-dimensional radiative transfer calculations to predict the polarization degree at ALMA bands due to scattered thermal emission. We find that the polarization pattern of a disk with a planetary gap after 1 Myr of dust evolution shows a distinctive three ring structure. Two narrow inner rings are located at the planet gap edges. For increasing observing wavelengths all three rings slightly change their position, where the innermost and outermost rings move inward. This distance is detectable comparing the results at ALMA bands 3, 6 and 7. Within the highest polarized intensity regions the polarization vectors are oriented in the azimuthal direction. For an inclined disk there is an interplay between polarization originating from a flux gradient and inclination-induced quadrupole polarization. For intermediate inclined transition disks the polarization degree is as high as ~ 2% at band 3, which is well above the detection limit of future ALMA observations.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The shadow of the Flying Saucer: A very low temperature for large dust grains

Dust determines the temperature structure of protoplanetary disks. However, dust temperature determinations almost invariably rely on a complex modeling of the Spectral Energy Distribution. We attempt a direct determination of the temperature of large grains emitting at mm wavelengths.} We observe the edge-on dust disk of the Flying Saucer, which appears in silhouette against the CO J=2-1 emission from a background molecular cloud in $ρ$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity allows us to directly measure the %absorbing dust temperature. The dust opacity can then be derived from the emitted continuum radiation. The dust disk absorbs the radiation from the CO clouds at several velocities. We derive very low dust temperatures, 5 to 7 K at radii around 100 au, which is much lower than most model predictions. The dust optical depth is $> 0.2$ at 230 GHz, and the scale height at 100 au is at least 8 au (best fit 13 au). However, the dust disk is very flat (flaring index -0.35), which is indicative of dust settling in the outer parts.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Asymmetric features in the protoplanetary disk MWC758

The study of dynamical processes in protoplanetary disks is essential to understand planet formation. In this context, transition disks are prime targets because they are at an advanced stage of disk clearing and may harbor direct signatures of disk evolution. In this paper, we aim to derive new constraints on the structure of the transition disk MWC 758, to detect non-axisymmetric features and understand their origin. We obtained infrared polarized intensity observations of the protoplanetary disk MWC 758 with SPHERE/VLT at 1.04 microns to resolve scattered light at a smaller inner working angle (0.093") and a higher angular resolution (0.027") than previously achieved. We observe polarized scattered light within 0.53" (148 au) down to the inner working angle (26 au) and detect distinct non-axisymmetric features but no fully depleted cavity. The two small-scale spiral features that were previously detected with HiCIAO are resolved more clearly, and new features are identified, including two that are located at previously inaccessible radii close to the star. We present a model based on the spiral density wave theory with two planetary companions in circular orbits. The best model requires a high disk aspect ratio (H/r~0.20 at the planet locations) to account for the large pitch angles which implies a very warm disk. Our observations reveal the complex morphology of the disk MWC758. To understand the origin of the detected features, the combination of high-resolution observations in the submillimeter with ALMA and detailed modeling is needed.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sequential planet formation in the HD 100546 protoplanetary disk?

Context. The disk around the Herbig Ae star, HD 100546, shows structures that suggest the presence of two companions in the disk at $\sim10$ and $\sim70$ AU. The outer companion seems to be in the act of formation. Aims. Our aims are to provide constraints on the age of the planets in HD 100546 and to explore the potential evidence for sequential planet formation in transition disks such as HD 100546. Methods. We compare the recent resolved continuum observations of the disk around HD 100546 with the results of dust evolution simulations using an analytical prescription for the shapes of gaps carved by massive planets. Results. An inner pressure bump must have been present since early in the disk lifetime to have good agreement between the dust evolution models and the continuum observations of HD 100546. This pressure bump may have resulted from the presence of a very massive planet ($\sim20 M_{\rm{Jup}}$), which formed early in the inner disk ($r\sim$10 AU). If only this single planet exists, the disk is likely to be old, comparable to the stellar age ($\sim$5-10 Myr). Another possible explanation is an additional massive planet in the outer disk ($r\sim70$ AU): either a low-mass outer planet ($\lesssim5 M_{\rm{Jup}}$) injected at early times, or a higher mass outer planet ($\gtrsim15 M_{\rm{Jup}}$) formed very recently, traps the right amount of dust in pressure bumps to reproduce the observations. In the latter case, the disk could be much younger ($\sim3.0$ Myr). Conclusions. In the case in which two massive companions are embedded in the disk around HD 100546, as suggested in the literature, the outer companion could be at least $\gtrsim$2.5 Myr younger than the inner companion.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Testing particle trapping in transition disks with ALMA

We present new Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations at 336GHz of two transition disks, SR21 and HD135344B. In combination with previous ALMA observations from Cycle 0 at 689GHz, we compare the visibility profiles at the two frequencies and calculate the spectral index ($α_{\rm{mm}}$). The observations of SR21 show a clear shift in the visibility nulls, indicating radial variations of the inner edge of the cavity at the two wavelengths. Notable radial variations of the spectral index are also detected for SR21 with values of $α_{\rm{mm}}{\sim}3.8-4.2$ in the inner region ($r<35$ AU) and $α_{\rm{mm}}{\sim}2.6-3.0$ outside. An axisymmetric ring (which we call the ring model) or a ring with the addition of an azimuthal Gaussian profile, for mimicking a vortex structure (which we call the vortex model), is assumed for fitting the disk morphology. For SR21, the ring model better fits the emission at 336GHz, conversely the vortex model better fits the 689GHz emission. For HD135344B, neither a significant shift in the null of the visibilities nor radial variations of $α_{\rm{mm}}$ are detected. Furthermore, for HD135344B, the vortex model fits both frequencies better than the ring model. However, the azimuthal extent of the vortex increases with wavelength, contrary to model predictions for particle trapping by anticyclonic vortices. For both disks, the azimuthal variations of $α_{\rm{mm}}$ remain uncertain to confirm azimuthal trapping. The comparison of the current data with a generic model of dust evolution that includes planet-disk interaction suggests that particles in the outer disk of SR21 have grown to millimetre sizes and have accumulated in a radial pressure bump, whereas with the current resolution there is not clear evidence of radial trapping in HD135344B, although it cannot be excluded either.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Variability and dust filtration in the transition disk J160421.7-213028 observed in optical scattered light

Context. Some of transition disks show asymmetric structures in thermal sub-millimetre emission and optical scattered light. These structures can be the result of planet(s) or companions embedded in the disk. Aims. We aim to detect and analyse the scattered light of the transition disk J160421.7-213028, identify disk structures, and compare the results with previous observations of this disk at other wavelengths. Methods. We obtained and analysed new polarised intensity observations of the transition disk J160421.7-213028 with VLT/SPHERE using the visible light instrument ZIMPOL at $R'$ band (0.626$μ$m). We probe the disk gap down to a radius of confidence of 0.1'' (${\sim}15$ AU at 145 pc). We interpret the results in the context of dust evolution when planets interact with the parental disk. Results. We observe a gap from 0.1 to 0.3'' (${\sim}15$ to 40AU) and a bright annulus as previously detected by HiCIAO $H$ band observations at $1.65μ$m. The radial width of the annulus is around $40$AU, and its centre is at ${\sim}61$AU from the central star. The peak of the reflected light at 0.626$μ$m is located 20 AU inward of the cavity detected in the submillimetre. In addition, we detect a dip at a position angle of ${\sim}46.2 \pm 5.4^\circ$. A dip was also detected with HiCIAO, but located at ${\sim}85^\circ$. If the dip observed with HiCIAO is the same, this suggests an average dip rotation of ${\sim}12^\circ/$year, which is inconsistent with the local Keplerian angular velocity of $\sim$0.8$^\circ$/yr at $\sim$61AU. Conclusions. The spatial discrepancy in the radial emission in J160421.7-213028 at different wavelengths is consistent with dust filtration at the outer edge of a gap carved by a massive planet. The dip rotation can be interpreted by fast variability of the inner disk and/or the presence of a warp or circumplanetary material of a planet at ${\sim}9.6$ AU.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

ALMA hints at the presence of two companions in the disk around HD 100546

HD 100546 is a well-studied Herbig Be star-disk system that likely hosts a close-in companion with compelling observational evidence for an embedded protoplanet at 68 AU. We present ALMA observations of the HD 100546 disk which resolve the gas and dust structure at (sub)mm wavelengths. The CO emission (at 345.795 GHz) originates from an extensive molecular disk (390+/-20 AU in radius) whereas the continuum emission is more compact (230+/-20 AU in radius) suggesting radial drift of the mm-sized grains. The CO emission is similar in extent to scattered light images indicating well-mixed gas and um-sized grains in the disk atmosphere. Assuming azimuthal symmetry, a single-component power-law model cannot reproduce the continuum visibilities. The visibilities and images are better reproduced by a double-component model: a compact ring with a width of 21 AU centered at 26 AU and an outer ring with a width of 75+/-3 AU centered at 190+/-3 AU. The influence of a companion and protoplanet on the dust evolution is investigated. The companion at 10 AU facilitates the accumulation of mm-sized grains within a compact ring, ~ 20 - 30 AU, by ~ 10 Myr. The injection of a protoplanet at 1 Myr hastens the ring formation (~ 1.2 Myr) and also triggers the development of an outer ring (~ 100 - 200 AU). These observations provide additional evidence for the presence of a close-in companion and hint at dynamical clearing by a protoplanet in the outer disk.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Gas and dust structures in protoplanetary disks hosting multiple planets

Transition disks have dust-depleted inner regions and may represent an intermediate step of an on-going disk dispersal process, where planet formation is probably in progress. Recent millimetre observations of transition disks reveal radially and azimuthally asymmetric structures, where micron- and millimetre-sized dust particles may not spatially coexist. These properties can be the result of particle trapping and grain growth in pressure bumps originating from the disk interaction with a planetary companion. The multiple features observed in some transition disks, such as SR 21, suggest the presence of more than one planet. We study the gas and dust distributions of a disk hosting two massive planets as function of different disk and dust parameters. Observational signatures, such as the spectral energy distribution, sub-millimetre, and polarised images are simulated for the various parameters. We confirm that planets can lead to particle trapping, although for a disk with high viscosity ($α_{\rm{turb}}=10^{-2}$), the planet should be more massive than $5 M_{\rm{Jup}}$ and dust fragmentation should occur with low efficiency ($v_{f}\sim30\rm{m s}^{-1}$). This will lead to a ring-like feature as observed in transition disks in the millimetre. When trapping occurs, we find that a smooth distribution of micron sized grains throughout the disk, sometimes observed in scattered light, can only happen if the combination of planet mass and turbulence is such that small grains are not fully filtered out. A high disk viscosity ($α_{\rm{turb}}=10^{-2}$) ensures a replenishment of the cavity in micron-sized dust, while for lower viscosity ($α_{\rm{turb}}=10^{-3}$), the planet mass is constrained to be less than $5 M_{\rm{Jup}}$. In these cases, the gas distribution is likely to show low-amplitude azimuthal asymmetries caused by disk eccentricity rather than by long-lived vortices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Millimetre spectral indices of transition disks and their relation to the cavity radius

Transition disks are protoplanetary disks with inner depleted dust cavities and excellent candidates to investigate the dust evolution under the existence of a pressure bump. A pressure bump at the outer edge of the cavity allows dust grains from the outer regions to stop their rapid inward migration towards the star and efficiently grow to millimetre sizes. Dynamical interactions with planet(s) have been one of the most exciting theories to explain the clearing of the inner disk. We look for evidence of the presence of millimetre dust particles in transition disks by measuring their spectral index with new and available photometric data. We investigate the influence of the size of the dust depleted cavity on the disk integrated millimetre spectral index. We present the 3mm photometric observations carried out with PdBI of four transition disks: LkHa330, UXTauA, LRLL31, and LRLL67. We use available values of their fluxes at 345GHz to calculate their spectral index, as well as the spectral index for a sample of twenty transition disks. We compare the observations with two kind of models. In the first set of models, we consider coagulation and fragmentation of dust in a disk in which a cavity is formed by a massive planet located at different positions. The second set of models assumes disks with truncated inner parts at different radius and with power-law dust size distributions, where the maximum size of grains is calculated considering turbulence as the source of destructive collisions. We show that the integrated spectral index is higher for transition disks than for regular protoplanetary disks. For transition disks, the probability that the measured spectral index is positively correlated with the cavity radius is 95%. High angular resolution imaging of transition disks is needed to distinguish between the dust trapping scenario and the truncated disk case.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Explaining millimeter-sized particles in brown dwarf disks

Planets have been detected around a variety of stars, including low-mass objects, such as brown dwarfs. However, such extreme cases are challenging for planet formation models. Recent sub-millimeter observations of disks around brown dwarf measured low spectral indices of the continuum emission that suggest that dust grains grow to mm-sizes even in these very low mass environments. To understand the first steps of planet formation in scaled-down versions of T-Tauri disks, we investigate the physical conditions that can theoretically explain the growth from interstellar dust to millimeter-sized grains in disks around brown dwarf. We modeled the evolution of dust particles under conditions of low-mass disks around brown dwarfs. We used coagulation, fragmentation and disk-structure models to simulate the evolution of dust, with zero and non-zero radial drift. For the non-zero radial drift, we considered strong inhomogeneities in the gas surface density profile that mimic long-lived pressure bumps in the disk. We studied different scenarios that could lead to an agreement between theoretical models and the spectral slope found by millimeter observations. We find that fragmentation is less likely and rapid inward drift is more significant for particles in brown dwarf disks than in T-Tauri disks. We present different scenarios that can nevertheless explain millimeter-sized grains. As an example, a model that combines the following parameters can fit the millimeter fluxes measured for brown dwarf disks: strong pressure inhomogeneities of $\sim$ 40% of amplitude, a small radial extent $\sim$ 15 AU, a moderate turbulence strength $α_{\mathrm{turb}}= 10^{-3}$, and average fragmentation velocities for ices $v_f = 10 m s^{-1}$.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Imaging diagnostics for transitional discs

Transitional discs are a special type of protoplanetary discs where planet formation is thought to be taking place. These objects feature characteristic inner cavities and/or gaps of a few tens of AUs in the sub-millimitre images of the disc. This signature suggests a localised depletion of matter in the disc that could be caused by planet formation processes. However, recent observations have revealed differences in the structures imaged at different wavelengths in some of these discs. In this paper, we aim to explain these observational differences using self-consistent physical 2-D hydrodynamical and dust evolution models of such objects, assuming their morphology is indeed generated by the presence of a planet. We use these models to derive the distribution of gas and dust in a theoretical planet-hosting disc, for various planet masses and orbital separations. We then simulate observations of the emitted and scattered light from these models with VLT/SPHERE ZIMPOL, Subaru/HiCIAO, VLT/VISIR and ALMA. We do this by first computing the full resolution images of the models at different wavelengths, and then simulating the observations accounting for the characteristics of each particular instrument. The presence of the planet generates pressure bumps in the gas distribution of the disc whose characteristics strongly depend on the planet mass and position. These bumps cause large grains to accumulate while small grains are allowed into inner regions. This spatial differentiation of the grain sizes explains the differences in the observations since different wavelengths and observing techniques trace different parts of the dust size distribution. Based on this effect, we conclude that the combination of visible/near-infrared polarimetric and sub-mm images is the best strategy to constrain the properties of the unseen planet responsible for the disc structure.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Lopsided dust rings in transition disks

Context. Particle trapping in local or global pressure maxima in protoplanetary disks is one of the new paradigms in the theory of the first stages of planet formation. However, finding observational evidence for this effect is not easy. Recent work suggests that the large ring-shaped outer disks observed in transition disk sources may in fact be lopsided and constitute large banana-shaped vortices. Aims. We wish to investigate how effective dust can accumulate along the azimuthal direction. We also want to find out if the size- sorting resulting from this can produce a detectable signatures at millimeter wavelengths. Methods. To keep the numerical cost under control we develop a 1+1D method in which the azimuthal variations are treated sepa- rately from the radial ones. The azimuthal structure is calculated analytically for a steady-state between mixing and azimuthal drift. We derive equilibration time scales and compare the analytical solutions to time-dependent numerical simulations. Results. We find that weak, but long-lived azimuthal density gradients in the gas can induce very strong azimuthal accumulations of dust. The strength of the accumulations depends on the Péclet number, which is the relative importance of advection and diffusion. We apply our model to transition disks and our simulated observations show that this effect would be easily observable with ALMA and in principle allows to put constraints on the strength of turbulence and the local gas density.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Protoplanetary Disk Structure With Grain Evolution: the ANDES Model

We present a self-consistent model of a protoplanetary disk: 'ANDES' ('AccretioN disk with Dust Evolution and Sedimentation'). ANDES is based on a flexible and extendable modular structure that includes 1) a 1+1D frequency-dependent continuum radiative transfer module, 2) a module to calculate the chemical evolution using an extended gas-grain network with UV/X-ray-driven processes surface reactions, 3) a module to calculate the gas thermal energy balance, and 4) a 1+1D module that simulates dust grain evolution. For the first time, grain evolution and time-dependent molecular chemistry are included in a protoplanetary disk model. We find that grain growth and sedimentation of large grains to the disk midplane lead to a dust-depleted atmosphere. Consequently, dust and gas temperatures become higher in the inner disk (R < 50 AU) and lower in the outer disk (R > 50 AU), in comparison with the disk model with pristine dust. The response of disk chemical structure to the dust growth and sedimentation is twofold. First, due to higher transparency a partly UV-shielded molecular layer is shifted closer to the dense midplane. Second, the presence of big grains in the disk midplane delays the freeze-out of volatile gas-phase species such as CO there, while in adjacent upper layers the depletion is still effective. Molecular concentrations and thus column densities of many species are enhanced in the disk model with dust evolution, e.g., CO2, NH2CN, HNO, H2O, HCOOH, HCN, CO. We also show that time-dependent chemistry is important for a proper description of gas thermal balance.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

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

Context: The global size and spatial distribution of dust is an important ingredient in the structure and evolution of protoplanetary disks and in the formation of larger bodies, such as planetesimals. Aims: We aim to derive simple equations that explain the global evolution of the dust surface density profile and the upper limit of the grain size distribution and which can readily be used for further modeling or for interpreting of observational data. Methods: We have developed a simple model that follows the upper end of the dust size distribution and the evolution of the dust surface density profile. This model is calibrated with state-of-the-art simulations of dust evolution, which treat dust growth, fragmentation, and transport in viscously evolving gas disks. Results: We find very good agreement between the full dust-evolution code and the toy model presented in this paper. We derive analytical profiles that describe the dust-to-gas ratios and the dust surface density profiles well in protoplanetary disks, as well as the radial flux by solid material "rain out", which is crucial for triggering any gravity assisted formation of planetesimals. We show that fragmentation is the dominating effect in the inner regions of the disk leading to a dust surface density exponent of -1.5, while the outer regions at later times can become drift-dominated, yielding a dust surface density exponent of -0.75. Our results show that radial drift is not efficient in fragmenting dust grains. This supports the theory that small dust grains are resupplied by fragmentation due to the turbulent state of the disk.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Can grain growth explain transition disks?

Aims: Grain growth has been suggested as one possible explanation for the diminished dust optical depths in the inner regions of protoplanetary "transition" disks. In this work, we directly test this hypothesis in the context of current models of grain growth and transport. Methods: A set of dust evolution models with different disk shapes, masses, turbulence parameters, and drift efficiencies is combined with radiative transfer calculations in order to derive theoretical spectral energy distributions (SEDs) and images. Results: We find that grain growth and transport effects can indeed produce dips in the infrared SED, as typically found in observations of transition disks. Our models achieve the necessary reduction of mass in small dust by producing larger grains, yet not large enough to be fragmenting efficiently. However, this population of large grains is still detectable at millimeter wavelengths. Even if perfect sticking is assumed and radial drift is neglected, a large population of dust grains is left behind because the time scales on which they are swept up by the larger grains are too long. This mechanism thus fails to reproduce the large emission cavities observed in recent millimeter-wave interferometric images of accreting transition disks.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Ring shaped dust accumulation in transition disks

Context.Transition disks are believed to be the final stages of protoplanetary disks, during which a forming planetary system or photoevaporation processes open a gap in the inner disk, drastically changing the disk structure. From theoretical arguments it is expected that dust growth, fragmentation and radial drift are strongly influenced by gas disk structure, and pressure bumps in disks have been suggested as key features that may allow grains to converge and grow efficiently. Aims. We want to study how the presence of a large planet in a disk influences the growth and radial distribution of dust grains, and how observable properties are linked to the mass of the planet. Methods. We combine two-dimensional hydrodynamical disk simulations of disk-planet interactions with state-of-the-art coagulation/fragmentation models to simulate the evolution of dust in a disk which has a gap created by a massive planet. We compute images at different wavelengths and illustrate our results using the example of the transition disk LkCa15. Results. The gap opened by a planet and the long-range interaction between the planet and the outer disk create a single large pressure bump outside the planetary orbit. Millimeter-sized particles form and accumulate at the pressure maximum and naturally produce ring-shaped sub-millimeter emission that is long-lived because radial drift no longer depletes the large grain population of the disk. For large planet masses around 9 $M_{\mathrm{Jup}}$, the pressure maximum and, therefore, the ring of millimeter particles is located at distances that can be more than twice the star-planet separation, creating a large spatial separation between the gas inner edge of the outer disk and the peak millimeter emission. Smaller grains do get closer to the gap and we predict how the surface brightness varies at different wavelengths.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Evolution of Gas and Dust in Protoplanetary Accretion Disks

Dust constitutes only about one percent of the mass of circumstellar disks, yet it is of crucial importance for the modeling of planet formation, disk chemistry, radiative transfer and observations. The initial growth of dust from sub-micron sized grains to planetesimals and also the radial transport of dust in disks around young stars is the topic of this thesis. Circumstellar dust is subject to radial drift, vertical settling, turbulent mixing, collisional growth, fragmentation and erosion. We approach this subject from three directions: analytical calculations, numerical simulations, and comparison to observations. We describe the physical and numerical concepts that go into a model which is able to simulate the radial and size evolution of dust in a gas disk which is viscously evolving over several million years. The resulting dust size distributions are compared to our analytical predictions and a simple recipe for obtaining steady-state dust size distributions is derived. With the numerical model at hand, we show that grain fragmentation can explain the fact that circumstellar disks are observed to be dust-rich for several million years. Finally, we investigate the challenges that observations present to the theory of grain evolution, namely that grains of millimeter sizes are observed at large distances from the star. We have found that under the assumption that radial drift is ineffective, we can reproduce some of the observed spectral indices and fluxes. Fainter objects point towards a reduced dust-to-gas ratio or lower dust opacities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The mm-colors of a young binary disk system in the Orion Nebula Cluster

We present new EVLA continuum observations at 7 mm of the 253-1536 binary disk system in the Orion Nebula Cluster. The measured fluxes were combined with data in the sub-mm to derive the millimeter spectral index of each individual disk component. We show how these observations can be used to test the models of dust evolution and early growth of solids in protoplanetary disks. Our analysis indicates that the disk with lower density and higher temperature hosts larger grains than the companion disk. This result is the opposite of what predicted by the dust evolution models. The models and observational results can be reconciled if the viscosity $α$-parameter differs by more than a factor of ten in the two disks, or if the distribution of solids in the disks is strongly affected by radial motions. This analysis can be applied to future high-angular resolution observations of young disks with EVLA and ALMA to provide even stronger observational constraints to the models of dust evolution in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The TW Hya Disk at 870 microns: Comparison of CO and Dust Radial Structures

We present high resolution (0.3" = 16 AU), high signal-to-noise ratio Submillimeter Array observations of the 870 microns (345 GHz) continuum and CO J=3--2 line emission from the protoplanetary disk around TW Hya. Using continuum and line radiative transfer calculations, those data and the multiwavelength spectral energy distribution are analyzed together in the context of simple two-dimensional parametric disk structure models. Under the assumptions of a radially invariant dust population and (vertically integrated) gas-to-dust mass ratio, we are unable to simultaneously reproduce the CO and dust observations with model structures that employ either a single, distinct outer boundary or a smooth (exponential) taper at large radii. Instead, we find that the distribution of millimeter-sized dust grains in the TW Hya disk has a relatively sharp edge near 60 AU, contrary to the CO emission (and optical/infrared scattered light) that extends to a much larger radius of at least 215 AU. We discuss some possible explanations for the observed radial distribution of millimeter-sized dust grains and the apparent CO-dust size discrepancy, and suggest that they may be hallmarks of substructure in the dust disk or natural signatures of the growth and radial drift of solids that might be expected for disks around older pre-main sequence stars like TW Hya.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Trapping dust particles in the outer regions of protoplanetary disks

In order to explain grain growth to mm sized particles and their retention in outer regions of protoplanetary disks, as it is observed at sub-mm and mm wavelengths, we investigate if strong inhomogeneities in the gas density profiles can slow down excessive radial drift and can help dust particles to grow. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million years time scales. In addition, we introduced a comparison between the pressure inhomogeneities considered in this work and the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study if there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. We present the favorable conditions to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of the Taurus Auriga or Ophiucus star forming regions.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Dust size distributions in coagulation/fragmentation equilibrium: Numerical solutions and analytical fits

Context. Grains in circumstellar disks are believed to grow by mutual collisions and subsequent sticking due to surface forces. Results of many fields of research involving circumstellar disks, such as radiative transfer calculations, disk chemistry, magneto-hydrodynamic simulations largely depend on the unknown grain size distribution. Aims. As detailed calculations of grain growth and fragmentation are both numerically challenging and computationally expensive, we aim to find simple recipes and analytical solutions for the grain size distribution in circumstellar disks for a scenario in which grain growth is limited by fragmentation and radial drift can be neglected. Methods. We generalize previous analytical work on self-similar steady-state grain distributions. Numerical simulations are carried out to identify under which conditions the grain size distributions can be understood in terms of a combination of power-law distributions. A physically motivated fitting formula for grain size distributions is derived using our analytical predictions and numerical simulations. Results. We find good agreement between analytical results and numerical solutions of the Smoluchowski equation for simple shapes of the kernel function. The results for more complicated and realistic cases can be fitted with a physically motivated "black box" recipe presented in this paper. Our results show that the shape of the dust distribution is mostly dominated by the gas surface density (not the dust-to-gas ratio), the turbulence strength and the temperature and does not obey an MRN type distribution.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Gas- and dust evolution in protoplanetary disks

Context. Current models of the size- and radial evolution of dust in protoplanetary disks generally oversimplify either the radial evolution of the disk (by focussing at one single radius or by using steady state disk models) or they assume particle growth to proceed monodispersely or without fragmentation. Further studies of protoplanetary disks - such as observations, disk chemistry and structure calculations or planet population synthesis models - depend on the distribution of dust as a function of grain size and radial position in the disk. Aims. We attempt to improve upon current models to be able to investigate how the initial conditions, the build-up phase, and the evolution of the protoplanetary disk influence growth and transport of dust. Methods. We introduce a new version of the model of Brauer et al. (2008) in which we now include the time-dependent viscous evolution of the gas disk, and in which more advanced input physics and numerical integration methods are implemented. Results. We show that grain properties, the gas pressure gradient, and the amount of turbulence are much more influencing the evolution of dust than the initial conditions or the build-up phase of the protoplanetary disk. We quantify which conditions or environments are favorable for growth beyond the meter size barrier. High gas surface densities or zonal flows may help to overcome the problem of radial drift, however already a small amount of turbulence poses a much stronger obstacle for grain growth.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Impact of grain evolution on the chemical structure of protoplanetary disks

We study the impact of dust evolution in a protoplanetary disk around a T Tauri star on the disk chemical composition. For the first time we utilize a comprehensive model of dust evolution which includes growth, fragmentation and sedimentation. Specific attention is paid to the influence of grain evolution on the penetration of the UV field in the disk. A chemical model that includes a comprehensive set of gas phase and grain surface chemical reactions is used to simulate the chemical structure of the disk. The main effect of the grain evolution on the disk chemical composition comes from sedimentation, and, to a lesser degree, from the reduction of the total grain surface area. The net effect of grain growth is suppressed by the fragmentation process which maintains a population of small grains, dominating the total grain surface area. We consider three models of dust properties. In model GS both growth and sedimentation are taken into account. In models A5 and A4 all grains are assumed to have the same size (10(-5) cm and 10(-4) cm, respectively) with constant gas-to-dust mass ratio of 100. Like in previous studies, the "three-layer" pattern (midplane, molecular layer, hot atmosphere) in the disk chemical structure is preserved in all models, but shifted closer to the midplane in models with increased grain size (GS and A4). Unlike other similar studies, we find that in models GS and A4 column densities of most gas-phase species are enhanced by 1-3 orders of magnitude relative to those in a model with pristine dust (A5), while column densities of their surface counterparts are decreased. We show that column densities of certain species, like C2H, HC(2n+1)N (n=0-3), H2O and some other molecules, as well as the C2H2/HCN abundance ratio which are accessible with Herschel and ALMA can be used as observational tracers of early stages of the grain evolution process in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Testing the theory of grain growth and fragmentation by millimeter observations of protoplanetary disks

Context. Observations at sub-millimeter and mm wavelengths will in the near future be able to resolve the radial dependence of the mm spectral slope in circumstellar disks with a resolution of around a few AU at the distance of the closest star-forming regions. Aims. We aim to constrain physical models of grain growth and fragmentation by a large sample of (sub-)mm observations of disks around pre-main sequence stars in the Taurus-Auriga and Ophiuchus star-forming regions. Methods. State-of-the-art coagulation/fragmentation and disk-structure codes are coupled to produce steady-state grain size distributions and to predict the spectral slopes at (sub-)mm wavelengths. Results. This work presents the first calculations predicting the mm spectral slope based on a physical model of grain growth. Our models can quite naturally reproduce the observed mm-slopes, but a simultaneous match to the observed range of flux levels can only be reached by a reduction of the dust mass by a factor of a few up to about 30 while keeping the gas mass of the disk the same. This dust reduction can either be due to radial drift at a reduced rate or during an earlier evolutionary time (otherwise the predicted fluxes would become too low) or due to efficient conversion of dust into larger, unseen bodies.

0 citations0 reviewsNo code linkedOpen access