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

Cornelis P. Dullemond

Cornelis P. Dullemond 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.GAastro-phastro-ph.IM

Catalog footprint

What is connected

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

32 published item(s)

preprint2022arXiv

Dust entrainment in magnetically and thermally driven disk winds

Magnetically and thermally driven disk winds have gained popularity in the light of the current paradigm of low viscosities in protoplanetary disks that nevertheless present large accretion rates even in the presence of inner cavities. The possibility of dust entrainment in these winds may explain recent scattered light observations and constitutes a way of dust transport towards outer regions of the disk. We aim to study the dust dynamics in these winds and explore the differences between photoevaporation and magnetically driven disk winds in this regard. We quantify maximum entrainable grain sizes, the flow angle, and the general detectability. We used the FARGO3D code to perform global, 2.5D axisymmetric, nonideal MHD simulations including ohmic and ambipolar diffusion. Dust was treated as a pressureless fluid. Synthetic observations were created with the radiative transfer code RADMC-3D. We find a significant difference in the dust entrainment efficiency of warm, ionized winds such as photoevaporation and magnetic winds including X-ray and extreme ultraviolet (XEUV) heating compared to cold magnetic winds. The maximum entrainable grain size varies from $3\,μ\mathrm{m}$ to $6\,μ\mathrm{m}$ for ionized winds to $1\,μ\mathrm{m}$ for cold magnetic winds. Dust grains in cold magnetic winds tend to flow along a shallower angle compared to the warm winds. With increasing distance to the central star, the dust entrainment efficiency decreases. Larger values of the turbulent viscosity increase the maximum grain size radius of possible dust entrainment. Our simulations indicate that diminishing dust content in the outer regions of the wind can be mainly attributed to the dust settling in the disk. In the synthetic images, the dusty wind appears as a faint, conical emission region which is brighter for a cold magnetic wind.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A multi-frequency ALMA characterization of substructures in the GM Aur protoplanetary disk

The protoplanetary disk around the T Tauri star GM Aur was one of the first hypothesized to be in the midst of being cleared out by a forming planet. As a result, GM Aur has had an outsized influence on our understanding of disk structure and evolution. We present 1.1 and 2.1 mm ALMA continuum observations of the GM Aur disk at a resolution of ~50 mas (~8 au), as well as HCO$^+$ $J=3-2$ observations at a resolution of ~100 mas. The dust continuum shows at least three rings atop faint, extended emission. Unresolved emission is detected at the center of the disk cavity at both wavelengths, likely due to a combination of dust and free-free emission. Compared to the 1.1 mm image, the 2.1 mm image shows a more pronounced "shoulder" near R~40 au, highlighting the utility of longer-wavelength observations for characterizing disk substructures. The spectral index $α$ features strong radial variations, with minima near the emission peaks and maxima near the gaps. While low spectral indices have often been ascribed to grain growth and dust trapping, the optical depth of GM Aur's inner two emission rings renders their dust properties ambiguous. The gaps and outer disk ($R>100$ au) are optically thin at both wavelengths. Meanwhile, the HCO$^+$ emission indicates that the gas cavity is more compact than the dust cavity traced by the millimeter continuum, similar to other disks traditionally classified as "transitional."

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Global axisymmetric simulations of photoevaporation and magnetically driven protoplanetary disk winds

Photoevaporation and magnetically driven winds are two independent mechanisms to remove mass from protoplanetary disks. In addition to accretion, the effect of these two principles acting concurrently could be significant and the transition between those two has not been extensively studied and quantified in the literature yet. In order to contribute to the understanding of disk winds, we present the phenomena emerging in the framework of two-dimensional axisymmetric, non-ideal magnetohydrodynamic simulations including EUV-/ X-ray driven photoevaporation. Of particular interest are the examination of the transition region between photoevaporation and magnetically driven wind, the possibility of emerging magneto-centrifugal wind effects, as well as the morphology of the wind itself depending on the strength of the magnetic field. We use the PLUTO code in a 2.5D axisymmetric configuration with additional treatment of EUV-/ X-ray heating and dynamic ohmic diffusion based on a semi-analytical chemical model. We identify the transition between both outflow types to occur for values of the initial plasma beta $β\geq 10^7$, while magnetically driven winds generally outperform photoevaporation for stronger fields. In our simulations we observe irregular and asymmetric outflows for stronger magnetic fields. In the weak field regime the photoevaporation rates are slightly lowered by perturbations of the gas density in the inner regions of the disk. Overall, our results predict a wind with a lever arm smaller than 1.5, consistent with a hot magneto-thermal wind. Stronger accretion flows are present for values of $β< 10^7$.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Importance of radiative effects in gap opening by planets in protoplanetary disks

Recent ALMA observations revealed concentric annular structures in several young class-II objects. In an attempt to produce the rings and gaps in some of these systems, they have been modeled numerically with a single embedded planet assuming a locally isothermal equation of state. This is often justified by observations targeting the irradiation-dominated outer regions of disks (approximately 100 au). We test this assumption by conducting hydrodynamics simulations of embedded planets in thin locally isothermal and radiative disks that mimic the systems HD 163296 and AS 209 in order to examine the effect of including the energy equation in a seemingly locally isothermal environment as far as planet-disk interaction is concerned. We find that modeling such disks with an ideal equation of state makes a difference in terms of the number of produced rings and the spiral arm contrast in the disk. Locally isothermal disks produce sharper annular or azimuthal features and overestimate a single planet's gap-opening capabilities by producing multiple gaps. In contrast, planets in radiative disks carve a single gap for typical disk parameters. Consequently, for accurate modeling of planets with semimajor axes up to about 100 au, radiative effects should be taken into account even in seemingly locally isothermal disks. In addition, for the case of AS 209, we find that the primary gap is significantly different between locally isothermal and radiative models. Our results suggest that multiple planets are required to explain the ring-rich structures in such systems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Disk Substructures at High Angular Resolution Project (DSHARP): II. Characteristics of Annular Substructures

The Disk Substructures at High Angular Resolution Project used ALMA to map the 1.25 millimeter continuum of protoplanetary disks at a spatial resolution of ~5 au. We present a systematic analysis of annular substructures in the 18 single-disk systems targeted in this survey. No dominant architecture emerges from this sample; instead, remarkably diverse morphologies are observed. Annular substructures can occur at virtually any radius where millimeter continuum emission is detected and range in widths from a few au to tens of au. Intensity ratios between gaps and adjacent rings range from near-unity to just a few percent. In a minority of cases, annular substructures co-exist with other types of substructures, including spiral arms (3/18) and crescent-like azimuthal asymmetries (2/18). No clear trend is observed between the positions of the substructures and stellar host properties. In particular, the absence of an obvious association with stellar host luminosity (and hence the disk thermal structure) suggests that substructures do not occur preferentially near major molecular snowlines. Annular substructures like those observed in DSHARP have long been hypothesized to be due to planet-disk interactions. A few disks exhibit characteristics particularly suggestive of this scenario, including substructures in possible mean-motion resonance and "double gap" features reminiscent of hydrodynamical simulations of multiple gaps opened by a planet in a low-viscosity disk.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The impact of planet wakes on the location and shape of the water iceline in a protoplanetary disk

Context: Planets in accretion disks can excite spiral shocks, and---if massive enough---open gaps in their vicinity. Both of these effects can influence the overall disk thermal structure. Aims: We model planets of different masses and semimajor axes in disks of various viscosities and accretion rates to examine their impact on disk thermodynamics and highlight the mutable, non-axisymmetric nature of icelines in systems with massive planets. Methods: We conduct a parameter study using numerical hydrodynamics simulations where we treat viscous heating, thermal cooling and stellar irradiation as additional source terms in the energy equation, with some runs including radiative diffusion. Our parameter space consists of a grid containing different combinations of planet and disk parameters. Results: Both gap opening and shock heating can displace the iceline, with the effects being amplified for massive planets in optically thick disks. The gap region can split an initially hot (T>170 K) disk into a hot inner disk and a hot ring just outside of the planet's location, while shock heating can reshape the originally axisymmetric iceline into water-poor islands along spirals. We also find that radiative diffusion does not alter the picture significantly in this context. Conclusions: Shock heating and gap opening by a planet can effectively heat up optically thick disks and in general move and/or reshape the water iceline. This can affect the gap structure and migration torques. It can also produce azimuthal features that follow the trajectory of spiral arms, creating hot zones, "islands" of vapor and ice around spirals which could affect the accretion or growth of icy aggregates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Spiral Density Waves in a Young Protoplanetary Disk

Gravitational forces are expected to excite spiral density waves in protoplanetary disks, disks of gas and dust orbiting young stars. However, previous observations that showed spiral structure were not able to probe disk midplanes, where most of the mass is concentrated and where planet formation takes place. Using the Atacama Large Millimeter/submillimeter Array we detected a pair of trailing symmetric spiral arms in the protoplanetary disk surrounding the young star Elias 2-27. The arms extend to the disk outer regions and can be traced down to the midplane. These millimeter-wave observations also reveal an emission gap closer to the star than the spiral arms. We argue that the observed spirals trace shocks of spiral density waves in the midplane of this young disk.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Submillimeter polarization observation of the protoplanetary disk around HD 142527

We present the polarization observations toward the circumstellar disk around HD 142527 by using Atacama Large Millimeter/submillimeter Array (ALMA) at the frequency of 343 GHz. The beam size is $0.51 " \times 0.44 "$, which corresponds to the spatial resolution of $\sim$ 71 $\times$ 62 AU. The polarized intensity displays a ring-like structure with a peak located on the east side with a polarization fraction of $P= 3.26 \pm 0.02$ %, which is different from the peak of the continuum emission from the northeast region. The polarized intensity is significantly weaker at the peak of the continuum where $P= 0.220 \pm 0.010$ %. The polarization vectors are in the radial direction in the main ring of the polarized intensity, while there are two regions outside at the northwest and northeast areas where the vectors are in the azimuthal direction. If the polarization vectors represent the magnetic field morphology, the polarization vectors indicate the toroidal magnetic field configuration on the main ring and the poloidal fields outside. On the other hand, the flip of the polarization vectors is predicted by the self-scattering of thermal dust emission due to the change of the direction of thermal radiation flux. Therefore, we conclude that self-scattering of thermal dust emission plays a major role in producing polarization at millimeter wavelengths in this protoplanetary disk. Also, this puts a constraint on the maximum grain size to be approximately 150 ${\rm μm}$ if we assume compact spherical dust grains.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A tunnel and a traffic jam: How transition disks maintain a detectable warm dust component despite the presence of a large planet-carved gap

We combined hydrodynamical simulations of planet-disk interactions with dust evolution models that include coagulation and fragmentation of dust grains over a large range of radii and derived observational properties using radiative transfer calculations. We studied the role of the snow line in the survival of the inner disk of transition disks. Inside the snow line, the lack of ice mantles in dust particles decreases the sticking efficiency between grains. As a consequence, particles fragment at lower collision velocities than in regions beyond the snow line. This effect allows small particles to be maintained for up to a few Myrs within the first astronomical unit. These particles are closely coupled to the gas and do not drift significantly with respect to the gas. For lower mass planets (1$M_{\rm{Jup}}$), the pre-transition appearance can be maintained even longer because dust still trickles through the gap created by the planet, moves invisibly and quickly in the form of relatively large grains through the gap, and becomes visible again as it fragments and gets slowed down inside of the snow line. The global study of dust evolution of a disk with an embedded planet, including the changes of the dust aerodynamics near the snow line, can explain the concentration of millimetre-sized particles in the outer disk and the survival of the dust in the inner disk if a large dust trap is present in the outer disk. This behaviour solves the conundrum of the combination of both near-infrared excess and ring-like millimetre emission observed in several transition disks.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Grain Growth in the Circumstellar Disks of the Young Stars CY Tau and DoAr 25

We present new results from the Disks@EVLA program for two young stars: CY Tau and DoAr 25. We trace continuum emission arising from their circusmtellar disks from spatially resolved observations, down to tens of AU scales, at λ = 0.9, 2.8, 8.0, and 9.8 mm for DoAr25 and at λ = 1.3, 2.8, and 7.1 mm for CY Tau. Additionally, we constrain the amount of emission whose origin is different from thermal dust emission from 5 cm observations. Directly from interferometric data, we find that observations at 7 mm and 1 cm trace emission from a compact disk while millimeter-wave observations trace an extended disk structure. From a physical disk model, where we characterize the disk structure of CY Tau and DoAr 25 at wavelengths shorter than 5 cm, we find that (1) dust continuum emission is optically thin at the observed wavelengths and over the spatial scales studied, (2) a constant value of the dust opacity is not warranted by our observations, and (3) a high-significance radial gradient of the dust opacity spectral index, β, is consistent with the observed dust emission in both disks, with low-β in the inner disk and high-β in the outer disk. Assuming that changes in dust properties arise solely due to changes in the maximum particle size (amax), we constrain radial variations of amax in both disks, from cm-sized particles in the inner disk (R < 40 AU) to millimeter sizes in the outer disk (R > 80 AU). These observational constraints agree with theoretical predictions of the radial-drift barrier, however, fragmentation of dust grains could explain our amax(R) constraints if these disks have lower turbulence and/or if dust can survive high-velocity collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Millimeter-wave polarization of protoplanetary disks due to dust scattering

We present a new method to constrain the grain size in protoplanetary disks with polarization observations at millimeter wavelengths. If dust grains are grown to the size comparable to the wavelengths, the dust grains are expected to have a large scattering opacity and thus the continuum emission is expected to be polarized due to self-scattering. We perform 3D radiative transfer calculations to estimate the polarization degree for the protoplanetary disks having radial Gaussian-like dust surface density distributions, which have been recently discovered. The maximum grain size is set to be $100 {\rm~μm}$ and the observing wavelength to be 870 ${\rm μm}$. We find that the polarization degree is as high as 2.5 % with a subarcsec spatial resolution, which is likely to be detected with near-future ALMA observations. The emission is polarized due to scattering of anisotropic continuum emission. The map of the polarization degree shows a double peaked distribution and the polarization vectors are in the radial direction in the inner ring and in the azimuthal direction in the outer ring. We also find the wavelength dependence of the polarization degree: the polarization degree is the highest if dust grains have a maximum size of $a_{\rm max}\simλ/2π$, where $λ$ is the observing wavelength. Hence, multi-wave and spatially resolved polarization observations toward protoplanetary disks enable us to put a constraint on the grain size. The constraint on the grain size from polarization observations is independent of or may be even stronger than that from the opacity index.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A critical analysis of shock models for chondrule formation

In recent years many models of chondrule formation have been proposed. One of those models is the processing of dust in shock waves in protoplanetary disks. In this model, the dust and the chondrule precursors are overrun by shock waves, which heat them up by frictional heating and thermal exchange with the gas. In this paper we reanalyze the nebular shock model of chondrule formation and focus on the downstream boundary condition. We show that for large-scale plane-parallel chondrule-melting shocks the postshock equilibrium temperature is too high to avoid volatile loss. Even if we include radiative cooling in lateral directions out of the disk plane into our model (thereby breaking strict plane-parallel geometry) we find that for a realistic vertical extent of the solar nebula disk the temperature decline is not fast enough. On the other hand, if we assume that the shock is entirely optically thin so that particles can radiate freely, the cooling rates are too high to produce the observed chondrules textures. Global nebular shocks are therefore problematic as the primary sources of chondrules.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Can dust coagulation trigger streaming instability?

Streaming instability can be a very efficient way of overcoming growth and drift barriers to planetesimal formation. However, it was shown that strong clumping, which leads to planetesimal formation, requires a considerable number of large grains. State-of-the-art streaming instability models do not take into account realistic size distributions resulting from the collisional evolution of dust. We investigate whether a sufficient quantity of large aggregates can be produced by sticking and what the interplay of dust coagulation and planetesimal formation is. We develop a semi-analytical prescription of planetesimal formation by streaming instability and implement it in our dust coagulation code based on the Monte Carlo algorithm with the representative particles approach. We find that planetesimal formation by streaming instability may preferentially work outside the snow line, where sticky icy aggregates are present. The efficiency of the process depends strongly on local dust abundance and radial pressure gradient, and requires a super-solar metallicity. If planetesimal formation is possible, the dust coagulation and settling typically need ~100 orbits to produce sufficiently large and settled grains and planetesimal formation lasts another ~1000 orbits. We present a simple analytical model that computes the amount of dust that can be turned into planetesimals given the parameters of the disk model.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Modeling dust growth in protoplanetary disks: The breakthrough case

Simple toy models are often not sufficient to cover the complexity of the dust coagulation process, and a number of numerical approaches are therefore used, among which integration of the Smoluchowski equation and various versions of Monte Carlo algorithm are the most popular. In this paper, we directly compare the Smoluchowski and Monte Carlo approaches and we find a general agreement for most of the coagulation problems. However, for the sweep-up growth driven by the "lucky" breakthrough mechanism, the methods exhibit very different resolution dependencies. With too few mass bins, the Smoluchowski algorithm tends to overestimate the growth rate and the probability of breakthrough. The Monte Carlo method is less resolution dependent in the growth timescale aspect but it tends to underestimate the breakthrough chance due to its limited dynamic mass range. We discuss the features and drawbacks of both the approaches, which may limit their astrophysical applications.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Resolved Multifrequency Radio Observations of GG Tau

We present sub-arcsecond resolution observations of continuum emission associated with the GG Tau quadruple star system at wavelengths of 1.3, 2.8, 7.3, and 50 mm. These data confirm that the GG Tau A binary is encircled by a circumbinary ring at a radius of 235 AU with a FWHM width of ~60 AU. We find no clear evidence for a radial gradient in the spectral shape of the ring, suggesting that the particle size distribution is spatially homogeneous on angular scales of ~0.1". A central point source, likely associated with the primary component (GG Tau Aa), exhibits a composite spectrum from dust and free-free emission. Faint emission at 7.3 mm is observed toward the low-mass star GG Tau Ba, although its origin remains uncertain. Using these measurements of the resolved, multifrequency emission structure of the GG Tau A system, models of the far-infrared to radio spectrum are developed to place constraints on the grain size distribution and dust mass in the circumbinary ring. The non-negligible curvature present in the ring spectrum implies a maximum particle size of 1-10 mm, although we are unable to place strong constraints on the distribution shape. The corresponding dust mass is 30-300 M_earth, at a temperature of 20-30 K. We discuss how this significant concentration of relatively large particles in a narrow ring at a large radius might be produced in a local region of higher gas pressures (i.e., a particle "trap") located near the inner edge of the circumbinary disk.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

A Major Asymmetric Dust Trap in a Transition Disk

The statistics of discovered exoplanets suggest that planets form efficiently. However, there are fundamental unsolved problems, such as excessive inward drift of particles in protoplanetary disks during planet formation. Recent theories invoke dust traps to overcome this problem. We report the detection of a dust trap in the disk around the star Oph IRS 48 using observations from the Atacama Large Millimeter/submillimeter Array (ALMA). The 0.44-millimeter-wavelength continuum map shows high-contrast crescent-shaped emission on one side of the star originating from millimeter-sized grains, whereas both the mid-infrared image (micrometer-sized dust) and the gas traced by the carbon monoxide 6-5 rotational line suggest rings centered on the star. The difference in distribution of big grains versus small grains/gas can be modeled with a vortex-shaped dust trap triggered by a companion.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Breaking through: The effects of a velocity distribution on barriers to dust growth

It is unknown how far dust growth can proceed by coagulation. Obstacles to collisional growth are the fragmentation and bouncing barriers. However, in all previous simulations of the dust-size evolution in protoplanetary disks, only the mean collision velocity has been considered, neglecting that a small but possibly important fraction of the collisions will occur at both much lower and higher velocities. We study the effect of the probability distribution of impact velocities on the collisional dust growth barriers. Assuming a Maxwellian velocity distribution for colliding particles to determine the fraction of sticking, bouncing, and fragmentation, we implement this in a dust-size evolution code. We also calculate the probability of growing through the barriers and the growth timescale in these regimes. We find that the collisional growth barriers are not as sharp as previously thought. With the existence of low-velocity collisions, a small fraction of the particles manage to grow to masses orders of magnitude above the main population. A particle velocity distribution softens the fragmentation barrier and removes the bouncing barrier. It broadens the size distribution in a natural way, allowing the largest particles to become the first seeds that initiate sweep-up growth towards planetesimal sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Constraints on the Radial Variation of Grain Growth in the AS 209 Circumstellar Disk

We present dust continuum observations of the protoplanetary disk surrounding the pre-main sequence star AS 209, spanning more than an order of magnitude in wavelength from 0.88 to 9.8 mm. The disk was observed with sub-arcsecond angular resolution (0.2"-0.5") to investigate radial variations in its dust properties. At longer wavelengths, the disk emission structure is notably more compact, providing model-independent evidence for changes in the grain properties across the disk. We find that physical models which reproduce the disk emission require a radial dependence of the dust opacity κ_ν. Assuming that the observed wavelength-dependent structure can be attributed to radial variations in the dust opacity spectral index (β), we find that β(R) increases from β<0.5 at \sim20 AU to β>1.5 for R>80 AU, inconsistent with a constant value of β across the disk (at the 10σ level). Furthermore, if radial variations of κ_ν are caused by particle growth, we find that the maximum size of the particle-size distribution (a_{max}) increases from sub-millimeter-sized grains in the outer disk (R>70 AU) to millimeter and centimeter-sized grains in the inner disk regions (R< 70 AU). We compare our observational constraint on a_{max}(R) with predictions from physical models of dust evolution in proto-planetary disks. For the dust composition and particle-size distribution investigated here, our observational constraints on a_{max}(R) are consistent with models where the maximum grain size is limited by radial drift.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Planetesimal formation by sweep-up: How the bouncing barrier can be beneficial to growth

The formation of planetesimals is often accredited to collisional sticking of dust grains. The exact process is unknown, as collisions between larger aggregates tend to lead to fragmentation or bouncing rather than sticking. Recent laboratory experiments have however made great progress in the understanding and mapping of the complex physics involved in dust collisions. We want to study the possibility of planetesimal formation using the results from the latest laboratory experiments, particularly by including the fragmentation with mass transfer effect, which might lead to growth even at high impact velocities. We present a new experimentally and physically motivated dust collision model capable of predicting the outcome of a collision between two particles of arbitrary masses and velocities. It is used together with a continuum dust-size evolution code that is both fast in terms of execution time and able to resolve the dust well at all sizes, allowing for all types of interactions to be studied without biases. We find that for the general dust population, bouncing collisions prevent the growth above millimeter-sizes. However, if a small number of cm-sized particles are introduced, for example due to vertical mixing or radial drift, they can act as a catalyst and start to sweep up the smaller particles. At a distance of 3 AU, 100-meter-sized bodies are formed on a timescale of 1 Myr. We conclude that direct growth of planetesimals might be a possibility thanks to a combination of the existence of a bouncing barrier and the fragmentation with mass transfer effect. The bouncing barrier is here even beneficial, as it prevents the growth of too many large particles that would otherwise only fragment among each other, and creates a reservoir of small particles that can be swept up by larger bodies. However, for this process to work, a few seeds of cm in size or larger have to be introduced.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing Interstellar Dust with Infrared Echoes from the Cas A Supernova

We present the analysis of an IRS 5-38 μm spectrum and MIPS photometric measurements of an infrared echo near the Cassiopeia A supernova remnant observed with the Spitzer Space Telescope. We have modeled the recorded echo accounting for PAHs, quantum-heated carbon and silicate grains, as well as thermal carbon and silicate particles. Using the fact that optical light echo spectroscopy has established that Cas A originated from a type IIb supernova explosion showing an optical spectrum remarkably similar to the prototypical type IIb SN 1993J, we use the latter to construct template data input for our simulations. We are then able to reproduce the recorded infrared echo spectrum by combining the emission of dust heated by the UV burst produced at the shock breakout after the core-collapse and dust heated by optical light emitted near the visual maximum of the supernova light curve, where the UV burst and optical light curve characteristics are based on SN 1993J. We find a mean density of \sim680 H cm^{-3} for the echo region, with a size of a few light years across. We also find evidence of dust processing in the form of a lack of small PAHs with less than \sim300 carbon atoms, consistent with a scenario of PAHs destruction by the UV burst via photodissociation at the estimated distance of the echo region from Cas A. Furthermore, our simulations suggest that the weak 11 μm features of our recorded infrared echo spectrum are consistent with a strong dehydrogenated state of the PAHs. This exploratory study highlights the potential of investigating dust processing in the interstellar medium through infrared echoes.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Synthetic observations of first hydrostatic cores in collapsing low-mass dense cores. I. Spectral energy distributions and evolutionary sequence

The low-mass star formation evolutionary sequence is relatively well-defined both from observations and theoretical considerations. The first hydrostatic core is the first protostellar equilibrium object that is formed during the star formation process. Using state-of-the-art radiation-magneto-hydrodynamic 3D adaptive mesh refinement calculations, we aim to provide predictions for the dust continuum emission from first hydrostatic cores. We investigate the collapse and the fragmentation of magnetized one solar mass prestellar dense cores and the formation and evolution of first hydrostatic cores using the RAMSES code. We use three different magnetization levels for the initial conditions, which cover a large variety of early evolutionary morphology, e.g., the formation of a disk or a pseudo-disk, outflow launching, and fragmentation. We post-process the dynamical calculations using the 3D radiative transfer code RADMC-3D. We compute spectral energy distributions and usual evolutionary stage indicators such as bolometric luminosity and temperature. We find that the first hydrostatic core lifetimes depend strongly on the initial magnetization level of the parent dense core. We derive, for the first time, spectral energy distribution evolutionary sequences from high-resolution radiation-magneto-hydrodynamic calculations. We show that under certain conditions, first hydrostatic cores can be identified from dust continuum emission at 24 microns and 70 microns. We also show that single spectral energy distributions cannot help to distinguish between the formation scenarios of the first hydrostatic core, i.e., between the magnetized and non-magnetized models. Spectral energy distributions are a first useful and direct way to target first hydrostatic core candidates but high-resolution interferometry is definitively needed to determine the evolutionary stage of the observed sources.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Understanding hydrogen recombination line observations with ALMA and EVLA

Hydrogen recombination lines are one of the major diagnostics of H II region physical properties and kinematics. In the near future, the Expanded Very Large Array (EVLA) and the Atacama Large Millimeter Array (ALMA) will allow observers to study recombination lines in the radio and sub-mm regime in unprecedented detail. In this paper, we study the properties of recombination lines, in particular at ALMA wavelengths. We find that such lines will lie in almost every wideband ALMA setup and that the line emission will be equally detectable in all bands. Furthermore, we present our implementation of hydrogen recombination lines in the adaptive-mesh radiative transfer code RADMC-3D. We particularly emphasize the importance of non-LTE (local thermodynamical equilibrium) modeling since non-LTE effects can drastically affect the line shapes and produce asymmetric line profiles from radially symmetric H II regions. We demonstrate how these non-LTE effects can be used as a probe of systematic motions (infall & outflow) in the gas. We use RADMC-3D to produce synthetic observations of model H II regions and study the necessary conditions for observing such asymmetric line profiles with ALMA and EVLA.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dust properties and disk structure of evolved protoplanetary disks in Cep OB2: Grain growth, settling, gas and dust mass, and inside-out evolution

We present Spitzer/IRS spectra of 31 TTS and IRAM/1.3mm observations for 34 low- and intermediate-mass stars in the Cep OB2 region. Including our previously published data, we analyze 56 TTS and the 3 intermediate-mass stars with silicate features in Tr 37 (~4 Myr) and NGC 7160 (~12 Myr). The silicate emission features are well reproduced with a mixture of amorphous (with olivine, forsterite, and silica stoichiometry) and crystalline grains (forsterite, enstatite). We explore grain size and disk structure using radiative transfer disk models, finding that most objects have suffered substantial evolution (grain growth, settling). About half of the disks show inside-out evolution, with either dust-cleared inner holes or a radially-dependent dust distribution, typically with larger grains and more settling in the innermost disk. The typical strong silicate features require nevertheless the presence of small dust grains, and could be explained by differential settling according to grain size, anomalous dust distributions, and/or optically thin dust populations within disk gaps. M-type stars tend to have weaker silicate emission and steeper SEDs than K-type objects. The inferred low dust masses are in a strong contrast with the relatively high gas accretion rates, suggesting global grain growth and/or an anomalous gas to dust ratio. Transition disks (TD) in the Cep OB2 region display strongly processed grains, suggesting that they are dominated by dust evolution and settling. Finally, the presence of rare but remarkable disks with strong accretion at old ages reveals that some very massive disks may still survive to grain growth, gravitational instabilities, and planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Modeling CO Emission: II. The Physical Characteristics that Determine the X factor in Galactic Molecular Clouds

We investigate how the X factor, the ratio of H_2 column density (NH2) to velocity-integrated CO intensity (W), is determined by the physical properties of gas in model molecular clouds (MCs). We perform radiative transfer calculations on chemical-MHD models to compute X. Using integrated NH2 and W reproduces the limited range in X found in observations, resulting in a mean value X=2\times10^20 s/cm^2/K^1/km^1 from the Galactic MC model. However, in limited velocity intervals, X can take on a much larger range due to CO line saturation. Thus, X strongly depends on both the range in gas velocities and volume densities. The temperature (T) variations within individual MCs do not strongly affect X, as dense gas contributes most to setting X. For fixed velocity and density structure, gas with higher T has higher W, yielding X ~ T^-1/2 for T~20-100 K. We demonstrate that the linewidth-size scaling relation does not influence the X factor - only the range in velocities is important. Clouds with larger linewidths, regardless of the linewidth-size relation, have a higher W, corresponding to a lower value of X, scaling roughly as X ~ sigma^-1/2. The "mist" model, consisting of optically thick cloudlets with well-separated velocities, does not accurately reflect the conditions in a turbulent MC. We propose that the observed cloud-average values of X ~ XGal is simply a result of the limited range in NH2, temperatures, and velocities found in Galactic MCs - a ~constant value of X therefore does not require any linewidth-size relation, or that MCs are virialized objects. Since gas properties likely differ (slightly) between clouds, masses derived through a standard X should only be considered as a rough first estimate. For temperatures T~10-20 K, velocity dispersions ~1-6 km/s, and NH2~2-20\times10^21 cm^-2, we find cloud-averaged X ~ 2-4\times10^20 s/cm^2/K^1/km^1 for Solar-metallicity models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Equilibrium Initialization and Stability of Three-Dimensional Gas Disks

We present a new systematic way of setting up galactic gas disks based on the assumption of detailed hydrodynamic equilibrium. To do this, we need to specify the density distribution and the velocity field which supports the disk. We first show that the required circular velocity has no dependence on the height above or below the midplane so long as the gas pressure is a function of density only. The assumption of disks being very thin enables us to decouple the vertical structure from the radial direction. Based on that, the equation of hydrostatic equilibrium together with the reduced Poisson equation leads to two sets of second-order non-linear differential equation, which are easily integrated to set-up a stable disk. We call one approach `density method' and the other one `potential method'. Gas disks in detailed balance are especially suitable for investigating the onset of the gravitational instability. We revisit the question of global, axisymmetric instability using fully three-dimensional disk simulations. The impact of disk thickness on the disk instability and the formation of spontaneously induced spirals is studied systematically with or without the presence of the stellar potential. In our models, the numerical results show that the threshold value for disk instability is shifted from unity to 0.69 for self-gravitating thick disks and to 0.75 for combined stellar and gas thick disks. The simulations also show that self-induced spirals occur in the correct regions and with the right numbers as predicted by the analytic theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Modeling CO Emission: I. CO as a Column Density Tracer and the X-Factor in Molecular Clouds

Theoretical and observational investigations have indicated that the abundance of carbon monoxide (CO) is very sensitive to intrinsic properties of the gaseous medium, such as density, metallicity, and the background UV field. In order to accurately interpret CO observations, it is thus important to understand how well CO traces the gas, which in molecular clouds (MCs) is predominantly molecular hydrogen (H2). Recent hydrodynamic simulations by Glover & Mac Low have explicitly followed the formation and destruction of molecules in model MCs under varying conditions, confirming that CO formation strongly depends on the cloud properties. Conversely, the H2 formation is primarily determined by the age of the MC. We apply radiative transfer calculations to these MC models in order to investigate the properties of CO line emission. We focus on integrated CO (J=1-0) intensities emerging from individual clouds, including its relationship to the total, H2, and CO column densities, as well as the "X factor," the ratio of H2 column density to CO intensity. Models with high CO abundances have a threshold CO intensity ~65 K km/s at sufficiently large extinctions. Clouds with low CO abundances show no such intensity thresholds. The distribution of H2 column densities are well described as log-normal functions, though the distributions of CO intensities and column densities are usually not log-normal. In general, the PDFs of the integrated intensity do not follow the distribution functions of CO column densities. In the model with Milky Way-like conditions, the X factor is in agreement with the near constant value determined from observations. In clouds with lower CO abundances the X factor can vary appreciably - sometimes by > 4 orders of magnitude. In models with high densities, the CO line is fully saturated, so that the X factor is directly proportional to the molecular column density.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Stellar-Mass-Dependent Disk Structure in Coeval Planet-Forming Disks

Previous studies suggest that the planet-forming disks around very-low-mass stars/brown dwarfs may be flatter than those around more massive stars, in contrast to model predictions of larger scale heights for gas-disks around lower-mass stars. We conducted a statistically robust study to determine whether there is evidence for stellar-mass-dependent disk structure in planet-forming disks. We find a statistically significant difference in the Spitzer/IRAC color distributions of disks around very-low-mass and low-mass stars all belonging to the same star-forming region, the Chamaeleon I star-forming region. We show that self consistently calculated disk models cannot fit the median spectral energy distributions (SEDs) of the two groups. These SEDs can be only explained by flatter disk models, consistent with the effect of dust settling in disks. We find that relative to the disk structure predicted for flared disks the required reduction in disk scale height is anti-correlated with the stellar mass, i.e. disks around lower-mass stars are flatter. Our results show that the initial and boundary conditions of planet formation are stellar-mass-dependent, an important finding that must be considered in planet formation models.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Understanding Spatial and Spectral Morphologies of Ultracompact H II Regions

The spatial morphology, spectral characteristics, and time variability of ultracompact H II regions provide strong constraints on the process of massive star formation. We have performed simulations of the gravitational collapse of rotating molecular cloud cores, including treatments of the propagation of ionizing and non-ionizing radiation. We here present synthetic radio continuum observations of H II regions from our collapse simulations, to investigate how well they agree with observation, and what we can learn about how massive star formation proceeds. We find that intermittent shielding by dense filaments in the gravitationally unstable accretion flow around the massive star leads to highly variable H II regions that do not grow monotonically, but rather flicker, growing and shrinking repeatedly. This behavior appears able to resolve the well-known lifetime problem. We find that multiple ionizing sources generally form, resulting in groups of ultracompact H II regions, consistent with observations. We confirm that our model reproduces the qualitative H II region morphologies found in surveys, with generally consistent relative frequencies. We also find that simulated spectral energy distributions (SEDs) from our model are consistent with the range of observed H II region SEDs, including both regions showing a normal transition from optically thick to optically thin emission, and those with intermediate spectral slopes. In our models, anomalous slopes are solely produced by inhomogeneities in the H II region, with no contribution from dust emission at millimeter or submillimeter wavelengths. We conclude that many observed characteristics of ultracompact H II regions appear consistent with massive star formation in fast, gravitationally unstable, accretion flows.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

A new raytracer for modeling AU-scale imaging of lines from protoplanetary disks

The material that formed the present-day Solar System originated in feeding zones in the inner Solar Nebula located at distances within ~20 AU from the Sun, known as the planet-forming zone. Meteoritic and cometary material contain abundant evidence for the presence of a rich and active chemistry in the planet-forming zone during the gas-rich phase of Solar System formation. It is a natural conjecture that analogs can be found amoung the zoo of protoplanetary disks around nearby young stars. The study of the chemistry and dynamics of planet formation requires: 1) tracers of dense gas at 100-1000 K and 2) imaging capabilities of such tracers with 5-100 (0.5-20 AU) milli-arcsec resolution, corresponding to the planet-forming zone at the distance of the closest star-forming regions. Recognizing that the rich infrared (2-200 micron) molecular spectrum recently discovered to be common in protoplanetary disks represents such a tracer, we present a new general raytracing code, RADLite, that is optimized for producing infrared line spectra and images from axisymmetric structures. RADLite can consistently deal with a wide range of velocity gradients, such as those typical for the inner regions of protoplanetary disks. The code is intended as a backend for chemical and excitation codes, and can rapidly produce spectra of thousands of lines for grids of models for comparison with observations. Such radiative transfer tools will be crucial for constraining both the structure and chemistry of planet-forming regions, including data from current infrared imaging spectrometers and extending to the Atacama Large Millimeter Array and the next generation of Extremely Large Telescopes, the James Webb Space Telescope and beyond.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

Evolutionary Signatures in the Formation of Low-Mass Protostars. II. Towards Reconciling Models and Observations

A long-standing problem in low-mass star formation is the "luminosity problem," whereby protostars are underluminous compared to the accretion luminosity expected both from theoretical collapse calculations and arguments based on the minimum accretion rate necessary to form a star within the embedded phase duration. Motivated by this luminosity problem, we present a set of evolutionary models describing the collapse of low-mass, dense cores into protostars, using the Young & Evans (2005) model as our starting point. We calculate the radiative transfer of the collapsing cores throughout the full duration of the collapse in two dimensions. From the resulting spectral energy distributions, we calculate standard observational signatures to directly compare to observations. We incorporate several modifications and additions to the original Young & Evans model in an effort to better match observations with model predictions. We find that scattering, 2-D geometry, mass-loss, and outflow cavities all affect the model predictions, as expected, but none resolve the luminosity problem. A cycle of episodic mass accretion, however, can resolve this problem and bring the model predictions into better agreement with observations. Standard assumptions about the interplay between mass accretion and mass loss in our model give star formation efficiencies consistent with recent observations that compare the core mass function (CMF) and stellar initial mass function (IMF). The combination of outflow cavities and episodic mass accretion reduce the connection between observational Class and physical Stage to the point where neither of the two common observational signatures (bolometric temperature and ratio of bolometric to submillimeter luminosity) can be considered reliable indicators of physical Stage.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The outcome of protoplanetary dust growth: pebbles, boulders, or planetesimals? I. Mapping the zoo of laboratory collision experiments

The growth processes from protoplanetary dust to planetesimals are not fully understood. Laboratory experiments and theoretical models have shown that collisions among the dust aggregates can lead to sticking, bouncing, and fragmentation. However, no systematic study on the collisional outcome of protoplanetary dust has been performed so far so that a physical model of the dust evolution in protoplanetary disks is still missing. We intend to map the parameter space for the collisional interaction of arbitrarily porous dust aggregates. This parameter space encompasses the dust-aggregate masses, their porosities and the collision velocity. With such a complete mapping of the collisional outcomes of protoplanetary dust aggregates, it will be possible to follow the collisional evolution of dust in a protoplanetary disk environment. We use literature data, perform own laboratory experiments, and apply simple physical models to get a complete picture of the collisional interaction of protoplanetary dust aggregates. In our study, we found four different types of sticking, two types of bouncing, and three types of fragmentation as possible outcomes in collisions among protoplanetary dust aggregates. We distinguish between eight combinations of porosity and mass ratio. For each of these cases, we present a complete collision model for dust-aggregate masses between 10^-12 and 10^2 g and collision velocities in the range 10^-4 to 10^4 cm/s for arbitrary porosities. This model comprises the collisional outcome, the mass(es) of the resulting aggregate(s) and their porosities. We present the first complete collision model for protoplanetary dust. This collision model can be used for the determination of the dust-growth rate in protoplanetary disks.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Chemistry and line emission of outer protoplanetary disks

The structure and chemistry of protoplanetary disks depends strongly on the nature of the central star around which it has formed. The dust temperature is mainly set by the stellar luminosity, while the chemistry of the upper disk layers depends on the amount of intercepted UV and X-ray flux. We will study the differences in chemistry, thermal sturcture and line emission around Herbig Ae/Be, T Tauri stars and low mass M dwarfs. Predictions will be made for future observations with SOFIA and Herschel.

0 citations0 reviewsNo code linkedOpen access