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

Gijs D. Mulders

Gijs D. Mulders 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.IMastro-ph.GA

Catalog footprint

What is connected

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

16 published item(s)

preprint2022arXiv

Pterodactyls: A Tool to Uniformly Search and Vet for Young Transiting Planets In TESS Primary Mission Photometry

Kepler's short-period exoplanet population has revealed evolutionary features such as the Radius Valley and the Hot Neptune desert that are likely sculpted by atmospheric loss over time. These findings suggest that the primordial planet population is different from the Gyr-old Kepler population, and motivates exoplanet searches around young stars. Here, we present pterodactyls , a data reduction pipeline specifically built to address the challenges in discovering exoplanets around young stars and to work with TESS Primary Mission 30-min cadence photometry, since most young stars were not pre-selected TESS 2-min cadence targets. pterodactyls builds on publicly available and tested tools in order to extract, detrend, search, and vet transiting young planet candidates. We search five clusters with known transiting planets: Tucana-Horologium Association, IC 2602, Upper Centaurus Lupus, Ursa Major and Pisces Eridani. We show that pterodactyls recovers seven out of the eight confirmed planets and one out of the two planet candidates, most of which were initially detected in 2-min cadence data. For these clusters, we conduct injection-recovery tests to characterize our detection efficiency, and compute an intrinsic planet occurrence rate of 49+-20% for sub-Neptunes and Neptunes (1.8-6 Re) within 12.5 days, which is higher than Kepler's Gyr-old occurrence rates of 6.8+-0.3%. This potentially implies that these planets have shrunk with time due to atmospheric mass loss. However, a proper assessment of the occurrence of transiting young planets will require a larger sample unbiased to planets already detected. As such, pterodactyls will be used in future work to search and vet for planet candidates in nearby clusters and moving groups.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Earths in Other Solar Systems N-body simulations: the Role of Orbital Damping in Reproducing the Kepler Planetary Systems

The population of exoplanetary systems detected by Kepler provides opportunities to refine our understanding of planet formation. Unraveling the conditions needed to produce the observed exoplanets will sallow us to make informed predictions as to where habitable worlds exist within the galaxy. In this paper, we examine using N-body simulations how the properties of planetary systems are determined during the final stages of assembly. While accretion is a chaotic process, trends in the ensemble properties of planetary systems provide a memory of the initial distribution of solid mass around a star prior to accretion. We also use EPOS, the Exoplanet Population Observation Simulator, to account for detection biases and show that different accretion scenarios can be distinguished from observations of the Kepler systems. We show that the period of the innermost planet, the ratio of orbital periods of adjacent planets, and masses of the planets are determined by the total mass and radial distribution of embryos and planetesimals at the beginning of accretion. In general, some amount of orbital damping, either via planetesimals or gas, during accretion is needed to match the whole population of exoplanets. Surprisingly, all simulated planetary systems have planets that are similar in size, showing that the "peas in a pod" pattern can be consistent with both a giant impact scenario and a planet migration scenario. The inclusion of material at distances larger than what Kepler observes has a profound impact on the observed planetary architectures, and thus on the formation and delivery of volatiles to possible habitable worlds.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Joint Radial Velocity and Direct Imaging Planet Yield Calculations: I. Self-consistent Planet Populations

Planet yield calculations may be used to inform the target selection strategy and science operations of space observatories. Forthcoming and proposed NASA missions, such as the Wide-Field Infrared Survey Telescope (WFIRST), the Habitable Exoplanet Imaging Mission (HabEx), and the Large UV/Optical/IR Surveyor (LUVOIR), are expected to be equipped with sensitive coronagraphs and/or starshades. We are developing a suite of numerical simulations to quantify the extent to which ground-based radial velocity (RV) surveys could boost the detection efficiency of direct imaging missions. In this paper, we discuss the first step in the process of estimating planet yields: generating synthetic planetary systems consistent with observed occurrence rates from multiple detection methods. In an attempt to self-consistently populate stars with orbiting planets, it is found that naive extrapolation of occurrence rates (mass, semi-major axis) results in an unrealistically large number-density of Neptune-mass planets beyond the ice-line ($a \gtrsim 5$au), causing dynamic interactions that would destabilize orbits. We impose a stability criterion for multi-planet systems based on mutual Hill radii separation. Considering the influence of compact configurations containing Jovian-mass and Neptune-mass planets results in a marked suppression in the number of terrestrial planets that can exist at large radii. This result has a pronounced impact on planet yield calculations particularly in regions accessible to high-contrast imaging and microlensing. The dynamically compact configurations and occurrence rates that we develop may be incorporated as input into joint RV and direct imaging yield calculations to place meaningful limits on the number of detectable planets with future missions.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

The Exoplanet Population Observation Simulator. II -- Population Synthesis in the Era of Kepler

The collection of planetary system properties derived from large surveys such as Kepler provides critical constraints on planet formation and evolution. These constraints can only be applied to planet formation models, however, if the observational biases and selection effects are properly accounted for. Here we show how epos, the Exoplanet Population Observation Simulator, can be used to constrain planet formation models by comparing the Bern planet population synthesis models to the Kepler exoplanetary systems. We compile a series of diagnostics, based on occurrence rates of different classes of planets and the architectures of multi-planet systems, that can be used as benchmarks for future and current modeling efforts. Overall, we find that a model with 100 seed planetary cores per protoplanetary disk provides a reasonable match to most diagnostics. Based on these diagnostics we identify physical properties and processes that would result in the Bern model more closely matching the known planetary systems. These are: moving the planet trap at the inner disk edge outward; increasing the formation efficiency of mini-Neptunes; and reducing the fraction of stars that form observable planets. We conclude with an outlook on the composition of planets in the habitable zone, and highlight that the majority of simulated planets smaller than 1.7 Earth radii have substantial hydrogen atmospheres. The software used in this paper is available online for public scrutiny at https://github.com/GijsMulders/epos

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Eruption of the Candidate Young Star ASASSN-15qi

Outbursts on young stars are usually interpreted as accretion bursts caused by instabilities in the disk or the star-disk connection. However, some protostellar outbursts may not fit into this framework. In this paper, we analyze optical and near-infrared spectra and photometry to characterize the 2015 outburst of the probable young star ASASSN-15qi. The $\sim 3.5$ mag brightening in the $V$ band was sudden, with an unresolved rise time of less than one day. The outburst decayed exponentially by 1 mag for 6 days and then gradually back to the pre-outburst level after 200 days. The outburst is dominated by emission from $\sim10,000$ K gas. An explosive release of energy accelerated matter from the star in all directions, seen in a spectacular cool, spherical wind with a maximum velocity of 1000 km/s. The wind and hot gas both disappeared as the outburst faded and the source returned to its quiescent F-star spectrum. Nebulosity near the star brightened with a delay of 10-20 days. Fluorescent excitation of H$_2$ is detected in emission from vibrational levels as high as $v=11$, also with a possible time delay in flux increase. The mid-infrared spectral energy distribution does not indicate the presence of warm dust emission, although the optical photospheric absorption and CO overtone emission could be related to a gaseous disk. Archival photometry reveals a prior outburst in 1976. Although we speculate about possible causes for this outburst, none of the explanations are compelling.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The SPHERE view of the planet-forming disk around HD100546

We image with unprecedented spatial resolution and sensitivity disk features that could be potential signs of planet-disk interaction. Two companion candidates have been claimed in the disk around the young Herbig Ae/Be star HD100546. Thus, this object serves as an excellent target for our investigation of the natal environment of giant planets. We exploit the power of extreme adaptive optics operating in conjunction with the new high-contrast imager SPHERE to image HD100546 in scattered light. We obtain the first polarized light observations of this source in the visible (with resolution as fine as 2 AU) and new H and K band total intensity images that we analyze with the Pynpoint package. The disk shows a complex azimuthal morphology, where multiple scattering of photons most likely plays an important role. High brightness contrasts and arm-like structures are ubiquitous in the disk. A double-wing structure (partly due to ADI processing) resembles a morphology newly observed in inclined disks. Given the cavity size in the visible (11 AU), the CO emission associated to the planet candidate 'c' might arise from within the circumstellar disk. We find an extended emission in the K band at the expected location of 'b'. The surrounding large-scale region is the brightest in scattered light. There is no sign of any disk gap associated to 'b'.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

A stellar-mass-dependent drop in planet occurrence rates

The Kepler Spacecraft has discovered a large number of planets up to one-year periods and down to terrestrial sizes. While the majority of the target stars are main-sequence dwarfs of spectral type F, G, and K, Kepler covers stars with effective temperature as low as 2500 K, which corresponds to M stars. These cooler stars allow characterization of small planets near the habitable zone, yet it is not clear if this population is representative of that around FGK stars. In this paper, we calculate the occurrence of planets around stars of different spectral types as a function of planet radius and distance from the star, and show that they are significantly different from each other. We further identify two trends: First, the occurrence of Earth to Neptune-sized planets is successively higher toward later spectral types at all orbital periods probed by Kepler; Planets around M stars occur twice as frequently as around G stars, and thrice as frequently as around F stars. Second, a drop in planet occurrence is evident at all spectral types inward of a 10 day orbital period, with a plateau further out. By assigning to each spectral type a median stellar mass, we show that the distance from the star where this drop occurs is stellar mass dependent, and scales with semi-major axis as the cube root of stellar mass. By comparing different mechanisms of planet formation, trapping and destruction, we find that this scaling best matches the location of the pre-main-sequence co-rotation radius, indicating efficient trapping of migrating planets or planetary building blocks. These results demonstrate the stellar-mass dependence of the planet population, both in terms of occurrence rate and of orbital distribution. The prominent stellar-mass dependence of the inner boundary of the planet population shows that the formation or migration of planets is sensitive to the stellar parameters.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

An Increase in the Mass of Planetary Systems around Lower-Mass Stars

Trends in the planet population with host star mass provide an avenue to constrain planet formation theories. We derive the planet radius distribution function for Kepler stars of different spectral types, sampling a range in host star masses. We find that M dwarf stars have 3.5 times more small planets (1.0-2.8 R_Earth) than main-sequence FGK stars, but two times fewer Neptune-sized and larger planets (>2.8 R_Earth). We find no systematic trend in the planet size distribution between spectral types F, G, and K to explain the increasing occurrence rates. Taking into account the mass-radius relationship and heavy-element mass of observed exoplanets, and assuming those are independent of spectral type, we derive the inventory of the heavy-element mass locked up in exoplanets at short orbits. The overall higher planet occurrence rates around M stars are not consistent with the redistribution of the same mass into more, smaller planets. At the orbital periods and planet radii where Kepler observations are complete for all spectral types, the average heavy-element mass locked up in exoplanets increases roughly inversely with stellar mass from 4 M_Earth in F stars to 5 M_Earth in G and K stars to 7 M_Earth in M stars. This trend stands in stark contrast with observed protoplanetary disk masses that decrease towards lower mass stars, and provides a challenge for current planet formation models. Neither models of in situ formation nor migration of fully-formed planets are consistent with these results. Instead, these results are indicative of large-scale inward migration of planetary building blocks --- either through type-I migration or radial drift of dust grains --- that is more efficient for lower mass stars, but does not result in significantly larger or smaller planets.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Optical imaging polarimetry of the LkCa 15 protoplanetary disk with SPHERE ZIMPOL

We present the first optical (590--890 nm) imaging polarimetry observations of the pre-transitional protoplanetary disk around the young solar analog LkCa 15, addressing a number of open questions raised by previous studies. We detect the previously unseen far side of the disk gap, confirm the highly eccentric scattered-light gap shape that was postulated from near-infrared imaging, at odds with the symmetric gap inferred from millimeter interferometry. Furthermore, we resolve the inner disk for the first time and trace it out to 30 AU. This new source of scattered light may contribute to the near-infrared interferometric signal attributed to the protoplanet candidate LkCa 15 b, which lies embedded in the outer regions of the inner disk. Finally, we present a new model for the system architecture of LkCa 15 that ties these new findings together. These observations were taken during science verification of SPHERE ZIMPOL and demonstrate this facility's performance for faint guide stars under adverse observing conditions.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The snow line in viscous disks around low-mass stars: implications for water delivery to terrestrial planets in the habitable zone

The water ice or snow line is one of the key properties of protoplanetary disks that determines the water content of terrestrial planets in the habitable zone. Its location is determined by the properties of the star, the mass accretion rate through the disk, and the size distribution of dust suspended in the disk. We calculate the snow line location from recent observations of mass accretion rates and as a function of stellar mass. By taking the observed dispersion in mass accretion rates as a measure of the dispersion in initial disk mass, we find that stars of a given mass will exhibit a range of snow line locations. At a given age and stellar mass, the observed dispersion in mass accretion rates of 0.4 dex naturally leads to a dispersion in snow line locations of 0.2 dex. For ISM-like dust sizes, the one-sigma snow line location among solar mass stars of the same age ranges from 2 to 5 au. For more realistic dust opacities that include larger grains, the snow line is located up to two times closer to the star. We use these locations and the outcome of N-body simulations to predict the amount of water delivered to terrestrial planets that formed in situ in the habitable zone. We find that the dispersion in snow line locations leads to a large range in water content. For ISM-like dust sizes, a significant fraction of habitable-zone terrestrial planets around sun-like stars remain dry, and no water is delivered to the habitable zones of low-mass M stars (less than half a solar mass) as in previous works. The closer-in snow line in disks with larger grains enables water delivery to the habitable zone for a significant fraction of M stars and all FGK stars. Considering their larger numbers and higher planet occurrence, M stars may host most of the water-rich terrestrial planets in the galaxy if these planets are able to hold on to their water in their subsequent evolution.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Volatile Delivery to Planets from Water-rich Planetesimals around Low Mass Stars

Most models of volatile delivery to accreting terrestrial planets assume that the carriers for water are similar in water content to the carbonaceous chondrites in our Solar System. Here we suggest that the water content of primitive bodies in many planetary systems may actually be much higher, as carbonaceous chondrites have lost some of their original water due to heating from short-lived radioisotopes that drove parent body alteration. Using N-body simulations, we explore how planetary accretion would be different if bodies beyond the water line contained a water mass fraction consistent with chemical equilibrium calculations, and more similar to comets, as opposed to the more traditional water-depleted values. We apply this model to consider planet formation around stars of different masses and identify trends in the properties of Habitable Zone planets and planetary system architecture which could be tested by ongoing exoplanet census data collection. Comparison of such data with the model predicted trends will serve to evaluate how well the N-body simulations and the initial conditions used in studies of planetary accretion can be used to understand this stage of planet formation.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Resolved images of the protoplanetary disk around HD 100546 with ALMA

The disk around the Herbig Ae/Be star HD 100546 has been extensively studied and it is one of the systems for which there are observational indications of ongoing and/or recent planet formation. However, up until now no resolved image of the millimeter dust emission or the gas has been published. We present the first resolved images of the disk around HD 100546 obtained in Band 7 with the ALMA observatory. The CO (3-2) image reveals a gas disk that extends out to 350 au radius at the 3-sigma level. Surprisingly, the 870um dust continuum emission is compact (radius <60 au) and asymmetric. The dust emission is well matched by a truncated disk with outer radius of $\approx$50 au. The lack of millimeter-sized particles outside the 60 au is consistent with radial drift of particles of this size. The protoplanet candidate, identified in previous high-contrast NACO/VLT L' observations, could be related to the sharp outer edge of the millimeter-sized particles. Future higher angular resolution ALMA observations are needed to determine the detailed properties of the millimeter emission and the gas kinematics in the inner region (<2arcsec). Such observations could also reveal the presence of a planet through the detection of circumplanetary disk material.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Planet or Brown Dwarf? Inferring the Companion Mass in HD 100546 from the Wall Shape using Mid-Infrared Interferometry

Giant planets form in protoplanetary disks while these disks are still gas-rich, and can reveal their presence through the annular gaps they carve out. HD 100546 is a gas-rich disk with a wide gap between between a radius of ~1 and 13 AU, possibly cleared out by a planetary companion or planetary system. We want to identify the nature of the unseen companion near the far end of the disk gap. We use mid-infrared interferometry at multiple baselines to constrain the curvature of the disk wall at the far end of the gap. We use 2D hydrodynamical simulations of embedded planets and brown dwarfs to estimate viscosity of the disk and the mass of a companion close to the disk wall. We find that the disk wall at the far end of the gap is not vertical, but rounded-off by a gradient in the surface density. Such a gradient can be reproduced in hydrodynamical simulations with a single, heavy companion (>=30...80 M_Jup) while the disk has viscosity of at least α >~ 5 10^-3. Taking into account the changes in the temperature structure after gap opening reduces the lower limit on the planet mass and disk viscosity to 20 M_Jup and α = 2 10^-3. The object in the disk gap of HD 100546 that shapes the disk wall is most likely a 60(+20)(-40) M_Jup brown dwarf, while the disk viscosity is estimated to be at least α = 2 10^-3. The disk viscosity is an important factor in estimating planetary masses from disk morphologies: more viscous disks need heavier planets to open an equally deep gap.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Probing the turbulent mixing strength in protoplanetary disks across the stellar mass range: no significant variations

Dust settling and grain growth are the first steps in the planet-formation process in protoplanetary disks. These disks are observed around stars with different spectral types, and there are indications that the disks around lower mass stars are significantly flatter, which could indicate that they settle and evolve faster, or in a different way. We aim to test this assumption by modeling the median spectral energy distributions (SEDs) of three samples of protoplanetary disks: around Herbig stars, T Tauri stars and brown dwarfs. We focus on the turbulent mixing strength to avoid a strong observational bias from disk and stellar properties that depend on stellar mass. We generated SEDs with the radiative transfer code MCMax, using a hydrostatic disk structure and settling the dust in a self-consistent way with the alpha-prescription to probe the turbulent mixing strength. We are able to fit all three samples with a disk with the same input parameters, scaling the inner edge to the dust evaporation radius and disk mass to millimeter photometry. The Herbig stars require a special treatment for the inner rim regions, while the T-Tauri stars require viscous heating, and the brown dwarfs lack a good estimate of the disk mass because only few millimeter detections exist. We find that the turbulent mixing strength does not vary across the stellar mass range for a fixed grain size distribution and gas-to-dust ratio. Regions with the same temperature have a self-similar vertical structure independent of stellar mass, but regions at the same distance from the central star appear more settled in disks around lower mass stars. We find a relatively low turbulent mixing strength of alpha = 10^(-4) for a standard grain size distribution, but our results are also consistent with alpha = 0.01 for a grain size distribution with fewer small grains or a lower gas-to-dust ratio.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Why circumstellar disks are so faint in scattered light: The case of HD 100546

Scattered light images of circumstellar disks play an important role in characterizing the planet forming environments around young stars. The characteristic size of the scattering dust grains can be estimated from the observed brightness asymmetry between the front and back side of the disk, for example using standard Mie theory. However such models often overpredict their brightness by one or two orders of magnitude, and have difficulty explaining very red disk colors. We aim to develop a dust model that explains simultaneously the observed disk surface brightness, colors and asymmetry in scattered light, focusing on constraining grain sizes. We use the 2D radiative transfer code MCMax with anisotropic scattering to explore the effects of grain size on synthetic scattered light images of circumstellar disks. We compare the results with surface brightness profiles of the protoplanetary disk HD 100546 in scattered light at wavelengths from 0.4 to 2.2 micron. We find that extreme forward scattering by micron sized particles lowers the /effective/ dust albedo and creates a faint and red disk that /appears/ only slightly forward scattering. For the outer (>~100 AU) disk of HD 100546 we derive a minimum grain size of 2.5 micron, likely present in the form of aggregates. Intermediate sized grains are too bright, whereas smaller grains are faint and scatter more isotropically, but also produce disk colors that are too blue. Observed surface brightness asymmetries alone are not sufficient to constrain the grain size in circumstellar disks. Additional information, such as the brightness and colors of the disk are needed to provide additional constraints.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The effects of disk and dust structure on observed polarimetric images of protoplanetary disks

Imaging polarimetry is a powerful tool for imaging faint circumstellar material. For a correct analysis of observations we need to fully understand the effects of dust particle parameters, as well as the effects of the telescope, atmospheric seeing, and assumptions about the data reduction and processing of the observed signal. Here we study the major effects of dust particle structure, size-dependent grain settling, and instrumental properties. We performed radiative transfer modeling using different dust particle models and disk structures. To study the influence of seeing and telescope diffraction we ran the models through an instrument simulator for the ExPo dual-beam imaging polarimeter mounted at the 4.2m William Herschel Telescope (WHT). Particle shape and size have a strong influence on the brightness and detectability of the disks. In the simulated observations, the central resolution element also contains contributions from the inner regions of the protoplanetary disk besides the unpolarized central star. This causes the central resolution element to be polarized, making simple corrections for instrumental polarization difficult. This effect strongly depends on the spatial resolution, so adaptive optics systems are needed for proper polarization calibration. We find that the commonly employed homogeneous sphere model gives results that differ significantly from more realistic models. For a proper analysis of the wealth of data available now or in the near future, one must properly take the effects of particle types and disk structure into account. The observed signal depends strongly on the properties of these more realistic models, thus providing a potentially powerful diagnostic. We conclude that it is important to correctly understand telescope depolarization and calibration effects for a correct interpretation of the degree of polarization.

0 citations0 reviewsNo code linkedOpen access