Source author record

Nienke van der Marel

Nienke van der Marel appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

35works
4topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

35 published item(s)

preprint2025arXiv

WIde Separation Planets In Time (WISPIT): Two directly imaged exoplanets around the Sun-like stellar binary WISPIT 1

Wide separation gas giant planets present a challenge to current planet formation theories, and the detection and characterisation of these systems enables us to constrain their formation pathways. The WIde Separation Planets In Time (WISPIT) survey aims to detect and characterise wide separation planetary-mass companions over a range of ages from <5 to 20 Myr around solar-type host stars at distances of 75-500 (median: 140) parsecs. The WISPIT survey carries out two 5 minute H-band exposures with the VLT/SPHERE instrument and IRDIS camera, separated by at least six months to identify co-moving companions via proper motion analysis. These two H-band observations in combination with a follow-up Ks-band observation were used to determine the colour-magnitude of the co-moving companions and to derive their masses by comparing to AMES-COND and AMES-DUSTY evolutionary tracks. We report the discovery of WISPIT 1b and WISPIT 1c, two gas giant exoplanets that are co-moving with the stellar binary WISPIT 1, which itself consists of a K4 star and M5.5 star in a multi-decadal orbit. The planets are at projected separations of 338 au and 840 au and have masses of 10 Mj and 4 Mj respectively. We identified two common proper motion planetary companions to a (previously unknown) stellar binary with a Sun-like primary. These targets are ideal for follow up characterisation with both ground and space-based telescopes. Monitoring of the orbit with the GRAVITY interferometer will place constraints on their eccentricity, and spectroscopic characterisation will identify the composition and metallicity, providing information on their formation pathways.

preprint2022arXiv

A major asymmetric ice trap in a planet-forming disk: III. First detection of dimethyl ether

The complex organic molecules (COMs) detected in star-forming regions are the precursors of the prebiotic molecules that can lead to the emergence of life. By studying COMs in more evolved protoplanetary disks we can gain a better understanding of how they are incorporated into planets. This paper presents ALMA band 7 observations of the dust and ice trap in the protoplanetary disk around Oph IRS 48. We report the first detection of dimethyl ether (CH3OCH3) in a planet-forming disk and a tentative detection of methyl formate (CH3OCHO). We determined column densities for the detected molecules and upper limits on non-detected species using the CASSIS spectral analysis tool. The inferred column densities of CH3OCH3 and CH3OCHO with respect to methanol (CH3OH) are of order unity, indicating unusually high abundances of these species compared to other environments. Alternatively, the 12CH3OH emission is optically thick and beam diluted, implying a higher CH3OH column density and a smaller emitting area than originally thought. The presence of these complex molecules can be explained by thermal ice sublimation, where the dust cavity edge is heated by irradiation and the full volatile ice content is observable in the gas phase. This work confirms the presence of oxygen-bearing molecules more complex than CH3OH in protoplanetary disks for the first time. It also shows that it is indeed possible to trace the full interstellar journey of COMs across the different evolutionary stages of star, disk, and planet formation.

preprint2022arXiv

A spatially-resolved large cavity of the J0337 protoplanetary disk in Perseus

We present Keck/NIRC2 $K_{\rm p}L_{\rm p}$ high-contrast imaging observations of a J0337 protoplanetary disk. The data discover the spatially-resolved large cavity, which is the second report among protoplanetary disks in the Perseus star forming region after the LkH$α$~330 system. Our data and forward modeling using RADMC-3D suggests $\sim80$~au for the cavity radius. There is discrepancy between J0337's SED and the modeled SED at $\sim10\micron$ and this suggests an unseen inner disk. We also searched for companions around J0337 but did not detect any companion candidates at separations between $0\farcs1$ and $2\farcs5$. The $L_{\rm p}$-band detection limit corresponds to $\sim20 M_{\rm Jup}$ at 60~au, $\sim9-10 M_{\rm Jup}$ at 90~au, and $\sim3 M_{\rm Jup}$ at $>120$~au. Compared with other young systems with large cavities such as PDS~70 and RX~J1604, multiple Jovian planets, a single eccentric Jovian planet, or a massive brown-dwarf at an inner separation could exist within the cavity.

preprint2022arXiv

Gap Opening and Inner Disk Structure in the Strongly Accreting Transition Disk of DM Tau

Large inner dust gaps in transition disks are frequently posited as evidence of giant planets sculpting gas and dust in the disk, or the opening of a gap by photoevaporative winds. Although the former hypothesis is strongly supported by the observations of planets and deep depletions in gas within the gap some disks, many T Tauri stars hosting transition disks accrete at rates typical for an undepleted disk, raising the question of how gap opening occurs in these objects. We thus present an analysis of the structure of the transition disk around the T Tauri star DM Tau, which is strongly accreting ($\sim 10^{-8.3}~\mathrm{M}_\odot~ \mathrm{yr}^{-1}$) and turbulent ($α=0.078 \pm 0.02$). Using the DALI thermochemical code, we fit disk models to simultaneously reproduce the accretion rate, high level of turbulence, the gas traced by ALMA band 6 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2--1 lines, and the observed dust emission from the mm continuum and spectral energy distribution. We find a shallow depletion in gas surface density of $\sim 10$ relative to the outer disk and a gas rich inner disk is consistent with the observations. The planet mass of $<1$ M$_\mathrm{Jup}$ implied by the gap depth is in tension with predictions for dust trapping in a highly viscous disk, which requires a more massive planet of of $\sim10$M$_\mathrm{Jup}$. Photoevaporative models including a dead zone can qualitatively reproduce some features of the DM Tau disk, but still struggle to explain the high accretion rates and the observed mm continuum flux.

preprint2022arXiv

High-resolution ALMA observations of transition disk candidates in Lupus

Transition disks with small inner dust cavities are interesting targets for the study of disk clearing mechanisms. Such disks have been identified through a deficit in the infrared part of their SED, but spatially resolved millimeter imaging is required to confirm the presence of an inner dust cavity. We use high-resolution ALMA observations of 30 mas resolution in Band 6 continuum and $^{12}$CO 2--1 emission of 10 transition disk candidates in the Lupus star forming region, in order to confirm the presence of inner dust cavities and infer the responsible mechanism. The continuum data are analyzed using visibility modeling and the SEDs are compared with radiative transfer models. Out of the six transition disk candidates selected from their SED, only one disk revealed an inner dust cavity of 4 au in radius. Three of the other disks are highly inclined, which limits the detectability of an inner dust cavity but it is also demonstrated to be the possible cause for the infrared deficit in their SED. The two remaining SED-selected disks are very compact, with dust radii of only $\sim$3 au. From the four candidates selected from low-resolution images, three new transition disks with large inner cavities $>$20 au are identified, bringing the total number of transition disks with large cavities in Lupus to 13. SED-selected transition disks with small cavities are biased towards highly inclined and compact disks, which casts doubt on the use of their occurrence rates in estimating dispersal timescales of photoevaporation. Using newly derived disk dust masses and radii, we re-evaluate the size-luminosity and $M_{\rm dust}-M_{\rm star}$ relations. These relations can be understood if the bright disks are dominated by disks with substructure whereas faint disks are dominated by drift-dominated disks. (Abridged)

preprint2022arXiv

Images of Embedded Jovian Planet Formation At A Wide Separation Around AB Aurigae

Direct images of protoplanets embedded in disks around infant stars provide the key to understanding the formation of gas giant planets like Jupiter. Using the Subaru Telescope and Hubble Space Telescope, we find evidence for a jovian protoplanet around AB Aurigae orbiting at a wide projected separation (93 au), likely responsible for multiple planet-induced features in the disk. Its emission is reproducible as reprocessed radiation from an embedded protoplanet. We also identify two structures located at 430-580 au that are candidate sites of planet formation. These data reveal planet formation in the embedded phase and a protoplanet discovery at wide, > 50 au separations characteristic of most imaged exoplanets. With at least one clump-like protoplanet and multiple spiral arms, the AB Aur system may also provide the evidence for a long-considered alternative to the canonical model for Jupiter's formation: disk (gravitational) instability.

preprint2022arXiv

The Origin of the Doppler-flip in HD 100546: a large scale spiral arm generated by an inner binary companion

Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner spiral and the previously reported CO emission Doppler-flip, which has been interpreted as the signature of an embedded protoplanet. Comparisons with hydrodynamical models indicate that this Doppler-flip is instead the kinematic counterpart of the spiral, which is likely generated by an inner companion inside the disk cavity.

preprint2022arXiv

Two Rings and a Marginally Resolved, 5 AU, Disk Around LkCa 15 Identified Via Near Infrared Sparse Aperture Masking Interferometry

Sparse aperture masking interferometry (SAM) is a high resolution observing technique that allows for imaging at and beyond a telescope's diffraction limit. The technique is ideal for searching for stellar companions at small separations from their host star; however, previous analysis of SAM observations of young stars surrounded by dusty disks have had difficulties disentangling planet and extended disk emission. We analyse VLT/SPHERE-IRDIS SAM observations of the transition disk LkCa\,15, model the extended disk emission, probe for planets at small separations, and improve contrast limits for planets. We fit geometrical models directly to the interferometric observables and recover previously observed extended disk emission. We use dynamic nested sampling to estimate uncertainties on our model parameters and to calculate evidences to perform model comparison. We compare our extended disk emission models against point source models to robustly conclude that the system is dominated by extended emission within 50 au. We report detections of two previously observed asymmetric rings at $\sim$17 au and $\sim$45 au. The peak brightness location of each model ring is consistent with the previous observations. We also, for the first time with imaging, robustly recover an elliptical Gaussian inner disk, previously inferred via SED fitting. This inner disk has a FWHM of ~5 au and a similar inclination and orientation as the outer rings. Finally, we recover no clear evidence for candidate planets. By modelling the extended disk emission, we are able to place a lower limit on the near infrared companion contrast of at least 1000.

preprint2021arXiv

Dust Traps and the Formation of Cavities in Transition Discs: A millimetre to sub-millimetre comparison survey

The origin of the inner dust cavities observed in transition discs remains unknown. The segregation of dust and size of the cavity is expected to vary depending on which clearing mechanism dominates grain evolution. We present the results from the Discs Down Under program, an 8.8 mm continuum Australia Telescope Compact Array (ATCA) survey targeting 15 transition discs with large (> 20 au) cavities, and compare the resulting dust emission to Atacama Large millimetre/sub-millimetre Array (ALMA) observations. Our ATCA observations resolve the inner cavity for 8 of the 14 detected discs. We fit the visibilities and reconstruct 1D radial brightness models for 10 sources with a S/N > 5sigma. We find that, for sources with a resolved cavity in both wavebands, the 8.8 mm and sub-mm brightness distributions peak at the same radius from the star. We suggest that a similar cavity size for 8.8 mm and sub-mm dust grains is due to a dust trap induced by the presence of a companion.

preprint2020arXiv

ALMA Observations of the Inner Cavity in the Protoplanetary Disk around Sz 84

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of a protoplanetary disk around the T Tauri star Sz~84 and analyses of the structures of the inner cavity in the central region of the dust disk. Sz~84's spectral energy distribution (SED) has been known to exhibit negligible infrared excess at $λ\lesssim$10~$μ$m due to the disk's cavity structure. Analyses of the observed visibilities of dust continuum at 1.3~mm and the SED indicate that the size of the cavity in the disk of large (millimeter size) dust grains is 8~au in radius and that in the disk of small (sub-micron size) dust grains is 60~au in radius. Furthermore, from the SED analyses, we estimate that the upper limit mass of small dust grains at $r<$60~au is less than $\sim$10$^{-3}$~$M_{\rm earth}$, which is $\lesssim$0.01~\% of the total (small~$+$~large) dust mass at $r<$60~au. These results suggest that large dust grains are dominant at $r<$60~au, implying that dust grains efficiently grow with less efficient fragmentation in this region, potentially due to weak turbulence and/or stickier dust grains. The balance of grain growth and dust fragmentation is an important factor for determining the size of large dust grains in protoplanetary disks, and thus Sz~84 could serve as a good testbed for investigations of grain growth in such disks.

preprint2020arXiv

Are the spiral arms in the MWC 758 protoplanetary disc driven by a companion inside the cavity?

Spiral arms in protoplanetary discs are thought to be linked to the presence of companions. We test the hypothesis that the double spiral arm morphology observed in the transition disc MWC 758 can be generated by an $\approx 10$ M$_{\rm Jup}$ companion on an eccentric orbit internal to the spiral arms. Previous studies on MWC 758 have assumed an external companion. We compare simulated observations from three dimensional hydrodynamics simulations of disc-companion interaction to scattered light, infrared and CO molecular line observations, taking into account observational biases. The inner companion hypothesis is found to explain the double spiral arms, as well as several additional features seen in MWC 758 -- the arc in the northwest, substructures inside the spiral arms, the cavity in CO isotopologues, and the twist in the kinematics. Testable predictions include detection of fainter spiral structure, detection of a point source south-southeast of the primary, and proper motion of the spiral arms.

preprint2020arXiv

Dust depleted inner disks in a large sample of transition disks through long-baseline ALMA observations

Transition disks with large inner dust cavities are thought to host massive companions. However, the disk structure inside the companion orbit and how material flows toward an actively accreting star remain unclear. We present a high resolution continuum study of inner disks in the cavities of 38 transition disks. Measurements of the dust mass from archival Atacama Large Millimeter/Submillimeter Array observations are combined with stellar properties and spectral energy distributions to assemble a detailed picture of the inner disk. An inner dust disk is detected in 18 of 38 disks in our sample. Of the 14 resolved disks, 9 are significantly misaligned with the outer disk. The near-infrared excess is uncorrelated with the mm dust mass of the inner disk. The size-luminosity correlation known for protoplanetary disks is recovered for the inner disks as well, consistent with radial drift. The inner disks are depleted in dust relative to the outer disk and their dust mass is uncorrelated with the accretion rates. This is interpreted as the result of radial drift and trapping by planets in a low $α$ ($\sim 10^{-3}$) disk, or a failure of the $α$-disk model to describe angular momentum transport and accretion. The only disk in our sample with confirmed planets in the gap, PDS 70, has an inner disk with a significantly larger radius and lower inferred gas-to-dust ratio than other disks in the sample. We hypothesize that these inner disk properties and the detection of planets are due to the gap having only been opened recently by young, actively accreting planets.

preprint2020arXiv

GW Ori: Interactions Between a Triple-star System and its Circumtriple Disk in Action

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

preprint2020arXiv

Protoplanetary disk masses in NGC 2024: Evidence for two populations

Protoplanetary disks in dense, massive star-forming regions are strongly affected by their environment. How this environmental impact changes over time is an important constraint on disk evolution and external photoevaporation models. We characterize the dust emission from 179 disks in the core of the young (0.5 Myr) NGC 2024 cluster. By studying how the disk mass varies within the cluster, and comparing these disks to those in other regions, we aim to determine how external photoevaporation influences disk properties over time. Using the Atacama Large Millimeter/submillimeter Array, a 2.9' x 2.9' mosaic centered on NGC 2024 FIR 3 was observed at 225 GHz with a resolution of 0.25'', or ~100 AU. The imaged region contains 179 disks identified at IR wavelengths, seven new disk candidates, and several protostars. The overall detection rate of disks is $32 \pm 4\%$. Few of the disks are resolved, with the exception of a giant (R = 300 AU) transition disk. Serendipitously, we observe a millimeter flare from an X-ray bright young stellar object (YSO), and resolve continuum emission from a Class 0 YSO in the FIR 3 core. Two distinct disk populations are present: a more massive one in the east, along the dense molecular ridge hosting the FIR 1-5 YSOs, with a detection rate of $45 \pm 7\%$. In the western population, towards IRS 1, only $15 \pm 4\%$ of disks are detected. NGC 2024 hosts two distinct disk populations. Disks along the dense molecular ridge are young (0.2 - 0.5 Myr) and partly shielded from the far ultraviolet radiation of IRS 2b; their masses are similar to isolated 1 - 3 Myr old SFRs. The western population is older and at lower extinctions, and may be affected by external photoevaporation from both IRS 1 and IRS 2b. However, it is possible these disks had lower masses to begin with.

preprint2020arXiv

Solving grain size inconsistency between ALMA polarization and VLA continuum in the Ophiuchus IRS 48 protoplanetary disk

The protoplanetary disk around Ophiuchus IRS 48 shows an azimuthally asymmetric dust distribution in (sub-)millimeter observations, which is interpreted as a vortex, where millimeter/centimeter-sized particles are trapped at the location of the continuum peak. In this paper, we present 860 $μ$m ALMA observations of polarized dust emission of this disk. The polarized emission was detected toward a part of the disk. The polarization vectors are parallel to the disk minor axis, and the polarization fraction was derived to be $1-2$\%. These characteristics are consistent with models of self-scattering of submillimeter-wave emission, which indicate a maximum grain size of $\sim100$ $μ$m. However, this is inconsistent with the previous interpretation of millimeter/centimeter dust particles being trapped by a vortex. To explain both, ALMA polarization and previous ALMA and VLA observations, we suggest that the thermal emission at 860 $μ$m wavelength is optically thick ($τ_{\rm abs}\sim7.3$) at the dust trap with the maximum observable grain size of $\sim100$ $μ$m rather than an optically thin case with $\sim$ cm dust grains. We note that we cannot rule out that larger dust grains are accumulated near the midplane if the 860 $μ$m thermal emission is optically thick.

preprint2019arXiv

Canada and the SKA from 2020-2030

This white paper submitted for the 2020 Canadian Long-Range Planning process (LRP2020) presents the prospects for Canada and the Square Kilometre Array (SKA) from 2020-2030, focussing on the first phase of the project (SKA1) scheduled to begin construction early in the next decade. SKA1 will make transformational advances in our understanding of the Universe across a wide range of fields, and Canadians are poised to play leadership roles in several. Canadian key SKA technologies will ensure a good return on capital investment in addition to strong scientific returns, positioning Canadian astronomy for future opportunities well beyond 2030. We therefore advocate for Canada's continued scientific and technological engagement in the SKA from 2020-2030 through participation in the construction and operations phases of SKA1.

preprint2019arXiv

Development Plans for the Atacama Large Millimeter/submillimeter Array (ALMA)

(abridged) The Atacama Large Millimeter/submillimeter Array (ALMA) was the top-ranked priority for a new ground-based facility in the 2000 Canadian Long Range Plan. Ten years later, at the time of LRP2010, ALMA construction was well underway, with first science observations anticipated for 2011. In the past 8 years, ALMA has proved itself to be a high-impact, high-demand observatory, with record numbers of proposals submitted to the annual calls and large numbers of highly cited scientific papers across fields from protoplanetary disks to high-redshift galaxies and quasars. The LRP2010 ALMA white paper laid out 8 specific metrics that could be used to judge the success of Canada's participation in ALMA. Among these metrics were publications (number; impact), collaborations (international; multi-wavelength), and student training. To call out one particular metric, Canadians are making excellent use of ALMA in training graduate students and postdocs: as of June 2018, 12 of 23 Canadian first-author papers were led by a graduate student, and a further 4 papers were led by postdocs. All 8 metrics argue for Canada's involvement in ALMA over the past decade to be judged a success. The successful achievement of these wide-ranging goals argues strongly for Canada's continuing participation in ALMA over the next decade and beyond. Looking forward, our community needs to: (1) maintain Canadian access to ALMA and our competitiveness in using ALMA; (2) preserve full Canadian funding for our share of ALMA operations; (3) identify components of ALMA development in which Canada can play a significant role, including stimulating expertise in submillimetre instrumentation to capitalize on future opportunities; and (4) keep Canadians fully trained and engaged in ALMA, as new capabilities become available, reaching the widest possible community of potential users.

preprint2019arXiv

LRP2020: Signposts of planet formation in protoplanetary disks

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

preprint2019arXiv

The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

preprint2019arXiv

The Next Generation Very Large Array

The next generation Very Large Array (ngVLA) is a transformational radio observatory being designed by the U.S. National Radio Astronomy Observatory (NRAO). It will provide order of magnitude improvements in sensitivity, resolution, and uv coverage over the current Jansky Very Large Array (VLA) at ~1.2-50 GHz and extend the frequency range up to 70-115 GHz. This document is a white paper written by members of the Canadian community for the 2020 Long Range Plan panel, which will be making recommendations on Canada's future directions in astronomy. Since Canadians have been historically major users of the VLA and have been valued partners with NRAO for ALMA, Canada's participation in ngVLA is welcome. Canadians have been actually involved in ngVLA discussions for the past five years, and have played leadership roles in the ngVLA Science and Technical Advisory Councils. Canadian technologies are also very attractive for the ngVLA, in particular our designs for radio antennas, receivers, correlates, and data archives, and our industrial capacities to realize them. Indeed, the Canadian designs for the ngVLA antennas and correlator/beamformer are presently the baseline models for the project. Given the size of Canada's radio community and earlier use of the VLA (and ALMA), we recommend Canadian participation in the ngVLA at the 7% level. Such participation would be significant enough to allow Canadian leadership in gVLA's construction and usage. Canada's participation in ngVLA should not preclude its participation in SKA; access to both facilities is necessary to meet Canada's radio astronomy needs. Indeed, ngVLA will fill the gap between those radio frequencies observable with the SKA and ALMA at high sensitivities and resolutions. Canada's partnership in ngVLA will give it access to cutting-edge facilities together covering approximately three orders of magnitude in frequency.

preprint2016arXiv

The (w)hole survey: an unbiased sample study of transition disk candidates based on Spitzer catalogs

Understanding disk evolution and dissipation is essential for studies of planet formation. Transition disks, i.e., disks with large dust cavities and gaps, are promising candidates of active evolution. About two dozen SED-selected candidates have been confirmed to have dust cavities through millimeter interferometric imaging, but this sample is biased towards the brightest disks. The Spitzer surveys of nearby low-mass star forming regions have resulted in more than 4000 Young Stellar Objects (YSOs). Using color criteria we have selected a sample of ~150 candidates, and an additional 40 candidates and known transition disks from the literature. The Spitzer data were complemented by new observations at longer wavelengths, including new JCMT and APEX submillimeter photometry, and WISE and Herschel-PACS mid and far-infrared photometry. Furthermore, optical spectroscopy was obtained and stellar types were derived for 85% of the sample, including information from the literature. The SEDs were fit to a grid of RADMC-3D disk models with a limited number of parameters: disk mass, inner disk mass, scale height and flaring, and disk cavity radius, where the latter is the main parameter of interest. A large fraction of the targets possibly have dust cavities based on the SED. The derived cavity sizes are consistent with imaging/modeling results in the literature, where available. Trends are found with Ldisk/Lstar and stellar mass and a possible connection with exoplanet orbital radii. A comparison with a previous study where color observables are used (Cieza et al. 2010) reveals large overlap between their category of planet-forming disks and our transition disks with cavities. A large number of the new transition disk candidates are suitable for follow-up observations with ALMA.

preprint2016arXiv

The Matryoshka Disk: Keck/NIRC2 Discovery of a Solar System-Scale, Radially Segregated Residual Protoplanetary Disk Around HD 141569A

Using Keck/NIRC2 $L^\prime$ (3.78 $μm$) data, we report the direct imaging discovery of a scattered light-resolved, solar system-scale residual protoplanetary disk around the young A-type star HD 141569A, interior to and concentric with the two ring-like structures at wider separations. The disk is resolved down to $\sim$ 0\farcs{}25 and appears as an arc-like rim with attached hook-like features. It is located at an angular separation intermediate between that of warm CO gas identified from spatially-resolved mid-infrared spectroscopy and diffuse dust emission recently discovered with the \textit{Hubble Space Telescope}. The inner disk has a radius of $\sim$ 39 AU, a position angle consistent with north-up, an inclination of $i$ $\sim$ 56$^{o}$, and has a center offset from the star. Forward-modeling of the disk favors a thick torus-like emission sharply truncated at separations beyond the torus' photocenter and heavily depleted at smaller separations. In particular, the best-fit density power law for the dust suggests that the inner disk dust and gas (as probed by CO) are radially segregated, a feature consistent with the dust trapping mechanism inferred from observations of "canonical" transitional disks. However, the inner disk component may instead be explained by radiation pressure-induced migration in optically-thin conditions, in contrast to the two stellar companion/planet-influenced ring-like structures at wider separations. HD 141569A's circumstellar environment --- with three nested, gapped, concentric dust populations --- is an excellent laboratory for understanding the relationship between planet formation and the evolution of both dust grains and disk architecture.

preprint2016arXiv

Vortices and spirals in the HD135344B transition disk

In recent years spiral structures have been seen in scattered light observations and signs of vortices in millimeter images of protoplanetary disks, both probably linked with the presence of planets. We present ALMA Band 7 (335 GHz or 0.89 mm) continuum observations of the transition disk HD135344B at unprecedented spatial resolution of 0.16", using superuniform weighting. The data show that the asymmetric millimeter dust ring seen in previous work actually consists of an inner ring and an outer asymmetric structure. The outer feature is cospatial with the end of one of the spiral arms seen in scattered light, but the feature itself is not consistent with a spiral arm due to its coradiance. We propose a new possible scenario to explain the observed structures at both wavelengths. Hydrodynamical simulations show that a massive planet can generate a primary vortex (which dissipates at longer timescales, becoming an axisymmetric ring) and trigger the formation of a second generation vortex further out. Within this scenario the two spiral arms observed at scattered light originate from a planet at ~30 AU and from the secondary vortex at ~75 AU rather than a planet further out as previously reported.

preprint2015arXiv

A Combined Spitzer and Herschel Infrared Study of Gas and Dust in the Circumbinary Disk Orbiting V4046 Sgr

We present results from a spectroscopic Spitzer and Herschel mid-to-far-infrared study of the circumbinary disk orbiting the evolved (age ~12-23 Myr) close binary T Tauri system V4046 Sgr. Spitzer IRS spectra show emission lines of [Ne II], H_2 S(1), CO_2 and HCN, while Herschel PACS and SPIRE spectra reveal emission from [O I], OH, and tentative detections of H_2O and high-J transitions of CO. We measure [Ne III]/[Ne II] < 0.13, which is comparable to other X-ray/EUV luminous T Tauri stars that lack jets. We use the H_2 S(1) line luminosity to estimate the gas mass in the relatively warm surface layers of the inner disk. The presence of [O I] emission suggests that CO, H_2O, and/or OH is being photodissociated, and the lack of [C I] emission suggests any excess C may be locked up in HCN, CN and other organic molecules. Modeling of silicate dust grain emission features in the mid-infrared indicates that the inner disk is composed mainly of large (r~5 um) amorphous pyroxene and olivine grains (~86% by mass) with a relatively large proportion of crystalline silicates. These results are consistent with other lines of evidence indicating that planet building is ongoing in regions of the disk within ~30 AU of the central, close binary.

preprint2015arXiv

A concentration of centimeter-sized grains in the Oph IRS 48 dust trap

Azimuthally asymmetric dust distributions observed with ALMA in transition disks have been interpreted as dust traps. We present VLA Ka band (34 GHz or 0.9 cm) and ALMA Cycle 2 Band 9 (680 GHz or 0.45 mm) observations at 0.2" resolution of the Oph IRS 48 disk, which suggest that larger particles could be more azimuthally concentrated than smaller dust grains, assuming an axisymmetric temperature field or optically thin 680 GHz emission. Fitting an intensity model to both data demonstrates that the azimuthal extent of the millimeter emission is 2.3 $\pm0.9$ times as wide as the centimeter emission, marginally consistent with the particle trapping mechanism under the above assumptions. The 34 GHz continuum image also reveals evidence for ionized gas emission from the star. Both the morphology and the spectral index variations are consistent with an increase of large particles in the center of the trap, but uncertainties remain due to the continuum optical depth at 680 GHz. Particle trapping has been proposed in planet formation models to allow dust particles to grow beyond millimeter sizes in the outer regions of protoplanetary disks. The new observations in the Oph IRS 48 disk provide support for the dust trapping mechanism for centimeter-sized grains, although additional data is required for definitive confirmation.

preprint2015arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: IV. Mechanical and radiative feedback

During the embedded stage of star formation, bipolar molecular outflows and UV radiation from the protostar are important feedback processes. Our aim is to quantify the feedback, mechanical and radiative, for a large sample of low-mass sources. The outflow activity is compared to radiative feedback in the form of UV heating by the accreting protostar to search for correlations and evolutionary trends. Large-scale maps of 26 young stellar objects, which are part of the Herschel WISH key program are obtained using the CHAMP+ instrument on the APEX (12CO and 13CO 6-5), and the HARP-B instrument on the JCMT (12CO and 13CO 3-2). Maps are used to determine outflow parameters and envelope models are used to quantify the amount of UV-heated gas and its temperature from 13CO 6-5 observations. All sources in our sample show outflow activity and the outflow force, F_CO, is larger for Class 0 sources than for Class I sources, even if their luminosities are comparable. The outflowing gas typically extends to much greater distances than the power-law envelope and therefore influences the surrounding cloud material directly. Comparison of the CO 6-5 results with Herschel-HIFI H2O and PACS high-J CO lines, both tracing currently shocked gas, shows that the two components are linked, even though the transitions do not probe the same gas. The link does not extend down to CO 3-2. The conclusion is that CO 6-5 depends on the shock characteristics (density and velocity), whereas CO 3-2 is more sensitive to conditions in the surrounding environment (density). The radiative feedback is responsible for increasing the gas temperature by a factor of two, up to 30-50 K, on scales of a few thousand AU, particularly along the direction of the outflow. The mass of the UV heated gas exceeds the mass contained in the entrained outflow in the inner ~3000 AU and is therefore at least as important on small scales.

preprint2015arXiv

Gas density drops inside dust cavities of transitional disks around young stars observed with ALMA

Transitional disks with large dust cavities are important laboratories to study planet formation and disk evolution. Cold gas may still be present inside these cavities, but the quantification of this gas is challenging. The gas content is important to constrain the origin of the dust cavity. We use Atacama Large Millimeter/submillimeter Array (ALMA) observations of 12CO 6--5 and 690 GHz (Band 9) continuum of five well-studied transitional disks. In addition, we analyze previously published Band 7 observations of a disk in 12CO 3--2 line and 345 GHz continuum. The observations are used to set constraints on the gas and dust surface density profiles, in particular the drop delta-gas of the gas density inside the dust cavity. The physical-chemical modeling code DALI is used to analyze the gas and dust images simultaneously. We model SR21, HD135344B, LkCa15, SR24S and RXJ1615-3255 (Band 9) and J1604-2130 (Band 7). The SED and continuum visibility curve constrain the dust surface density. Subsequently, the same model is used to calculate the 12CO emission, which is compared with the observations through spectra and intensity cuts. The amount of gas inside the cavity is quantified by varying the delta-gas parameter. Model fits to the dust and gas indicate that gas is still present inside the dust cavity for all disks but at a reduced level. The gas surface density drops inside the cavity by at least a factor 10, whereas the dust density drops by at least a factor 1000. Disk masses are comparable with previous estimates from the literature, cavity radii are found to be smaller than in the 345 GHz SubMillimeter Array (SMA) data. The derived gas surface density profiles suggest clearing of the cavity by one or more companions in all cases, trapping the millimeter-sized dust at the edge of the cavity.

preprint2015arXiv

Resolved gas cavities in transitional disks inferred from CO isotopologues with ALMA

Transitional disks around young stars are promising candidates to look for recently formed, embedded planets. Planet-disk interaction models predict that planets clear a gap in the gas while trapping dust at larger radii. Other physical mechanisms could be responsible for cavities as well. Previous observations have revealed that gas is still present inside these cavities, but the spatial distribution of this gas remains uncertain. We present high spatial resolution observations with the Atacama Large Millimeter/submillimeter Array (ALMA) of 13CO and C18O lines of four well-studied transitional disks. The observations are used to set constraints on the gas surface density, specifically cavity size and density drop inside the cavity. The physical-chemical model DALI is used to analyze the gas images of SR21, HD135344B, DoAr44 and IRS48. The main parameters of interest are the size, depth and shape of the gas cavity. CO isotope-selective photodissociation is included to properly constrain the surface density in the outer disk from C18O emission. The gas cavities are up to 3 times smaller than those of the dust in all four disks. Model fits indicate that the surface density inside the gas cavities decreases by a factor of 100-10000 compared with the surface density profile derived from the outer disk. A comparison with an analytical model of gap depths by planet-disk interaction shows that the disk viscosities are likely low, with a<1E-3 for planet masses <10 MJup. The resolved measurements of the gas and dust in transition disk cavities support the predictions of models that describe how planet-disk interactions sculpt gas disk structures and influence the evolution of dust grains. These observed structures strongly suggest the presence of giant planetary companions in transition disk cavities, although at smaller orbital radii than is typically indicated from the dust cavity radii alone.

preprint2014arXiv

SEEDS Adaptive Optics Imaging of the Asymmetric Transition Disk Oph IRS 48 in Scattered Light

We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $\sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$μ$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$\pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{\odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.

preprint2014arXiv

Warm formaldehyde in the Oph IRS 48 transitional disk

Simple molecules like H2CO and CH3OH in protoplanetary disks are the starting point for the production of more complex organic molecules. So far, the observed chemical complexity in disks has been limited due to freeze out of molecules onto grains in the bulk of the cold outer disk. Complex molecules can be studied more directly in transitional disks with large inner holes, as these have a higher potential of detection, through UV heating of the outer disk and the directly exposed midplane at the wall. We use Atacama Large Millimeter/submillimeter Array (ALMA) Band 9 (~680 GHz) line data of the transitional disk Oph IRS 48, previously shown to have a large dust trap, to search for complex molecules in regions where planetesimals are forming. We report the detection of the H2CO 9(1,8)-8(1,7) line at 674 GHz, which is spatially resolved as a semi-ring at ~60 AU radius centered south from the star. The inferred H2CO abundance is ~10^{-8} derived by combining a physical disk model of the source with a non-LTE excitation calculation. Upper limits for CH3OH lines in the same disk give an abundance ratio H2CO/CH3OH>0.3, which points to both ice formation and gas-phase routes playing a role in the H2CO production. Upper limits on the abundances of H13CO+, CN and several other molecules in the disk are also derived and found to be consistent with full chemical models. The detection of the H2CO line demonstrates the start of complex organic molecules in a planet-forming disk. Future ALMA observations should be able to push down the abundance detection limits of other molecules by 1-2 orders of magnitude and test chemical models of organic molecules in (transitional) disks.

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.

preprint2013arXiv

Gas structure inside dust cavities of transition disks: Oph IRS 48 observed by ALMA

(Abridged) Transition disks are recognized by the absence of emission of small dust grains inside a radius of up to several 10s of AUs. Due to the lack of angular resolution and sensitivity, the gas content of such dust holes has not yet been determined, but is of importance to constrain the mechanism leading to the dust holes. Transition disks are thought to currently undergo the process of dispersal, setting an end to the giant planet formation process. We present new high-resolution observations with the Atacama Large Millimeter/ submillimeter Array (ALMA) of gas lines towards the transition disk Oph IRS 48 previously shown to host a large dust trap. ALMA has detected the $J=6-5$ line of $^{12}$CO and C$^{17}$O around 690 GHz (434 $μ$m) at a resolution of $\sim$0.25$''$ corresponding to $\sim$30 AU (FWHM). The observed gas lines are used to set constraints on the gas surface density profile. New models of the physical-chemical structure of gas and dust in Oph IRS 48 are developed to reproduce the CO line emission together with the spectral energy distribution (SED) and the VLT-VISIR 18.7 $μ$m dust continuum images. Integrated intensity cuts and the total spectrum from models having different trial gas surface density profiles are compared to observations. Using the derived surface density profiles, predictions for other CO isotopologues are made, which can be tested by future ALMA observations of the object. The derived gas surface density profile points to the clearing of the cavity by one or more massive planet/companion rather than just photoevaporation or grain-growth.

preprint2013arXiv

Outflow forces of low mass embedded objects in Ophiuchus: a quantitative comparison of analysis methods

The outflow force of molecular bipolar outflows is a key parameter in theories of young stellar feedback on their surroundings. The focus of many outflow studies is the correlation between the outflow force, bolometric luminosity and envelope mass. However, it is difficult to combine the results of different studies in large evolutionary plots over many orders of magnitude due to the range of data quality, analysis methods and corrections for observational effects such as opacity and inclination. We aim to determine the outflow force for a sample of low luminosity embedded sources. We will quantify the influence of the analysis method and the assumptions entering the calculation of the outflow force. We use the James Clerk Maxwell Telescope to map 12CO J=3-2 over 2'x2' regions around 16 Class I sources of a well-defined sample in Ophiuchus at 15" resolution. The outflow force is then calculated using seven different methods differing e.g. in the use of intensity-weighted emission and correction factors for inclination. The results from the analysis methods differ from each other by up to a factor of 6, whereas observational properties and choices in the analysis procedure affect the outflow force by up to a factor of 4. For the sample of Class I objects, bipolar outflows are detected around 13 sources including 5 new detections, where the three non-detections are confused by nearby outflows from other sources. When combining outflow forces from different studies, a scatter by up to a factor of 5 can be expected. Although the true outflow force remains unknown, the separation method (separate calculation of dynamical time and momentum) is least affected by the uncertain observational parameters. The correlations between outflow force, bolometric luminosity and envelope mass are further confirmed down to low luminosity sources.

preprint2012arXiv

APEX-CHAMP+ high-J CO observations of low-mass young stellar objects: III. NGC 1333 IRAS 4A/4B envelope, outflow and UV heating

NGC 1333 IRAS 4A and IRAS 4B sources are among the best studied Stage 0 low-mass protostars which are driving prominent bipolar outflows. Most studies have so far concentrated on the colder parts (T<30K) of these regions. The aim is to characterize the warmer parts of the protostellar envelope in order to quantify the feedback of the protostars on their surroundings in terms of shocks, UV heating, photodissociation and outflow dispersal. Fully sampled large scale maps of the region were obtained; APEX-CHAMP+ was used for 12CO 6-5, 13CO 6-5 and [CI] 2-1, and JCMT-HARP-B for 12CO 3-2 emissions. Complementary Herschel-HIFI and ground-based lines of CO and its isotopologs, from 1-0 upto 10-9 (Eu/k 300K), are collected at the source positions. Radiative-transfer models of the dust and lines are used to determine temperatures and masses of the outflowing and UV-heated gas and infer the CO abundance structure. Broad CO emission line profiles trace entrained shocked gas along the outflow walls, with typical temperatures of ~100K. At other positions surrounding the outflow and the protostar, the 6-5 line profiles are narrow indicating UV excitation. The narrow 13CO 6-5 data directly reveal the UV heated gas distribution for the first time. The amount of UV-heated and outflowing gas are found to be comparable from the 12CO and 13CO 6-5 maps, implying that UV photons can affect the gas as much as the outflows. Weak [CI] emission throughout the region indicates a lack of CO dissociating photons. Modeling of the C18O lines indicates the necessity of a "drop" abundance profile throughout the envelope where the CO freezes out and is reloaded back into the gas phase, thus providing quantitative evidence for the CO ice evaporation zone around the protostars. The inner abundances are less than the canonical value of CO/H_2=2.7x10^-4, indicating some processing of CO into other species on the grains.

preprint2011arXiv

Complex molecules toward low-mass protostars: the Serpens core

Gas-phase complex organic molecules are commonly detected toward high-mass protostellar hot cores. Detections toward low-mass protostars and outflows are comparatively rare, and a larger sample is key to investigate how the chemistry responds to its environment. Guided by the prediction that complex organic molecules form in CH3OH-rich ices and thermally or non-thermally evaporate with CH3OH, we have identified three sight-lines in the Serpens core - SMM1, SMM4 and SMM4-W - which are likely to be rich in complex organics. Using the IRAM 30m telescope, narrow lines (FWHM of 1-2 km s-1) of CH3CHO and CH3OCH3 are detected toward all sources, HCOOCH3 toward SMM1 and SMM4-W, and C2H5OH not at all. Beam-averaged abundances of individual complex organics range between 0.6 and 10% with respect to CH3OH when the CH3OH rotational temperature is applied. The summed complex organic abundances also vary by an order of magnitude, with the richest chemistry toward the most luminous protostar SMM1. The range of abundances compare well with other beam-averaged observations of low-mass sources. Complex organic abundances are of the same order of magnitude toward low-mass protostars and high-mass hot cores, but HCOOCH3 is relatively more important toward low-mass protostars. This is consistent with a sequential ice photochemistry, dominated by CHO-containing products at low temperatures and early times.