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

Farzana Meru

Farzana Meru 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-ph.IMphysics.comp-phphysics.data-anphysics.space-ph

Catalog footprint

What is connected

29works
7topics
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

29 published item(s)

preprint2026arXiv

Disc fragmentation. II. Ejection of low mass Free Floating Planets from growing binary systems

Observations indicate that disc fragmentation due to Gravitational Instability (GI) is the likely origin of massive companions to stars, such as giant planets orbiting M-dwarf stars, Brown Dwarf (BD) companions to FGK stars, and binary stars with separations smaller than 100 au. Additionally, we have recently showed that disc fragmentation in young rapidly evolving binary systems ejects an abundant population of massive Jupiter-mass Free-Floating Planets (FFPs). In this model, a massive disc around an initially single protostar undergoes GI and hatches a number of fragments; the most massive oligarch grows by runaway accretion into the secondary star. As the system rearranges itself from a single to a binary star configuration, a dramatic "pincer movement" by the binary ejects planets through dynamical interactions with the stars. Here we propose that the same scenario applies to an even more abundant population of smaller FFPs discovered by the microlensing surveys. Although disc fragmentation is usually believed to form only massive objects, several pathways for forming small core-dominated planets at separations of tens of au exist. We present results from three complementary simulation approaches, all of which confirm planet ejection efficiency as high as 0.5 for secondaries more massive than $\sim 10$\% of the primary star mass. On the other hand, Jovian mass planets migrate through the region of tens of au too rapidly to eject planets from that region. We discuss implications of this scenario, concluding that microlensing FFPs may be the most convincing evidence yet that disc fragmentation forms planets much less massive than Jupiter.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories

ALMA observations of dust ring/gap structures in a minority but growing sample of protoplanetary disks can be explained by the presence of planets at large disk radii - yet the origins of these planets remains debated. We perform planet formation simulations using a semi-analytic model of the HL Tau disk to follow the growth and migration of hundreds of planetary embryos initially distributed throughout the disk, assuming either a high or low turbulent $α$ viscosity. We have discovered that there is a bifurcation in the migration history of forming planets as a consequence of varying the disk viscosity. In our high viscosity disks, inward migration prevails and yields compact planetary systems, tempered only by planet trapping at the water iceline around 5 au. In our lower viscosity models however, low mass planets can migrate outward to twice their initial orbital radii, driven by a radially extended region of strong outward-directed corotation torques located near the heat transition (where radiative heating of the disk by the star is comparable to viscous heating) - before eventually migrating inwards. We derive analytic expressions for the planet mass at which the corotation torque dominates, and find that this "corotation mass" scales as $M_{\rm p, corot} \sim α^{2/3}$. If disk winds dominate the corotation torque, the corotation mass scales linearly with wind strength. We propose that the observed bifurcation in disk demographics into a majority of compact dust disks and a minority of extended ring/gap systems is a consequence of a distribution of viscosity across the disk population.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Warping Away Gravitational Instabilities in Protoplanetary Discs

We perform 3D SPH simulations of warped, non-coplanar gravitationally unstable discs to show that as the warp propagates through the self-gravitating disc, it heats up the disc rendering it gravitationally stable. Thus losing their spiral structure and appearing completely axisymmetric. In their youth, protoplanetary discs are expected to be massive and self-gravitating, which results in non-axisymmetric spiral structures. However recent observations of young protoplanetary discs with ALMA have revealed that discs with large-scale spiral structure are rarely observed in the midplane. Instead, axisymmetric discs with some also having ring & gap structures are more commonly observed. Our work invloving warps, non-coplanar disc structures that are expected to commonly occur in young discs, potentially resolves this discrepancy between observations and theoretical predictions. We demonstrate that they are able to suppress the large-scale spiral structure of self-gravitating protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

A remnant planetary core in the hot-Neptune desert

The interiors of giant planets remain poorly understood. Even for the planets in the Solar System, difficulties in observation lead to large uncertainties in the properties of planetary cores. Exoplanets that have undergone rare evolutionary processes provide a route to understanding planetary interiors. Planets found in and near the typically barren hot-Neptune 'desert' (a region in mass-radius space that contains few planets) have proved to be particularly valuable in this regard. These planets include HD149026b, which is thought to have an unusually massive core, and recent discoveries such as LTT9779b and NGTS-4b, on which photoevaporation has removed a substantial part of their outer atmospheres. Here we report observations of the planet TOI-849b, which has a radius smaller than Neptune's but an anomalously large mass of $39.1^{+2.7}_{-2.6}$ Earth masses and a density of $5.2^{+0.7}_{-0.8}$ grams per cubic centimetre, similar to Earth's. Interior structure models suggest that any gaseous envelope of pure hydrogen and helium consists of no more than $3.9^{+0.8}_{-0.9}$ per cent of the total planetary mass. The planet could have been a gas giant before undergoing extreme mass loss via thermal self-disruption or giant planet collisions, or it could have avoided substantial gas accretion, perhaps through gap opening or late formation. Although photoevaporation rates cannot account for the mass loss required to reduce a Jupiter-like gas giant, they can remove a small (a few Earth masses) hydrogen and helium envelope on timescales of several billion years, implying that any remaining atmosphere on TOI-849b is likely to be enriched by water or other volatiles from the planetary interior. We conclude that TOI-849b is the remnant core of a giant planet.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Constraining planet formation around 6$M_{\odot}$-8$M_{\odot}$ stars

Identifying planets around O-type and B-type stars is inherently difficult; the most massive known planet host has a mass of only about $3M_{\odot}$. However, planetary systems which survive the transformation of their host stars into white dwarfs can be detected via photospheric trace metals, circumstellar dusty and gaseous discs, and transits of planetary debris crossing our line-of-sight. These signatures offer the potential to explore the efficiency of planet formation for host stars with masses up to the core-collapse boundary at $\approx 8M_{\odot}$, a mass regime rarely investigated in planet formation theory. Here, we establish limits on where both major and minor planets must reside around $\approx 6M_{\odot}-8M_{\odot}$ stars in order to survive into the white dwarf phase. For this mass range, we find that intact terrestrial or giant planets need to leave the main sequence beyond approximate minimum star-planet separations of respectively about 3 and 6 au. In these systems, rubble pile minor planets of radii 10, 1.0, and 0.1 km would have been shorn apart by giant branch radiative YORP spin-up if they formed and remained within, respectively, tens, hundreds and thousands of au. These boundary values would help distinguish the nature of the progenitor of metal-pollution in white dwarf atmospheres. We find that planet formation around the highest mass white dwarf progenitors may be feasible, and hence encourage both dedicated planet formation investigations for these systems and spectroscopic analyses of the highest mass white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Massive discs around low-mass stars

We use a suite of SPH simulations to investigate the susceptibility of protoplanetary discs to the effects of self-gravity as a function of star-disc properties. We also include passive irradiation from the host star using different models for the stellar luminosities. The critical disc-to-star mass ratio for axisymmetry (for which we produce criteria) increases significantly for low-mass stars. This could have important consequences for increasing the potential mass reservoir in a proto Trappist-1 system, since even the efficient Ormel et al. (2017) formation model will be influenced by processes like external photoevaporation, which can rapidly and dramatically deplete the dust reservoir. The aforementioned scaling of the critical $M_d/M_*$ for axisymmetry occurs in part because the Toomre $Q$ parameter has a linear dependence on surface density (which promotes instability) and only an $M_*^{1/2}$ dependence on shear (which reduces instability), but also occurs because, for a given $M_d/M_*$, the thermal evolution depends on the host star mass. The early phase stellar irradiation of the disc (for which the luminosity is much higher than at the zero age main sequence, particularly at low stellar masses) can also play a key role in significantly reducing the role of self-gravity, meaning that even Solar mass stars could support axisymmetric discs a factor two higher in mass than usually considered possible. We apply our criteria to the DSHARP discs with spirals, finding that self-gravity can explain the observed spirals so long as the discs are optically thick to the host star irradiation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planet Migration in Self-Gravitating Discs: Survival of Planets

We carry out three-dimensional SPH simulations to study whether planets can survive in self-gravitating protoplanetary discs. The discs modelled here use a cooling prescription that mimics a real disc which is only gravitationally unstable in the outer regions. We do this by modelling the cooling using a simplified method such that the cooling time in the outer parts of the disc is shorter than in the inner regions, as expected in real discs. We find that both giant (> M_Sat) and low mass (< M_Nep) planets initially migrate inwards very rapidly, but are able to slow down in the inner gravitationally stable regions of the disc without needing to open up a gap. This is in contrast to previous studies where the cooling was modelled in a more simplified manner where regardless of mass, the planets were unable to slow down their inward migration. This shows the important effect the thermodynamics has on planet migration. In a broader context, these results show that planets that form in the early stages of the discs' evolution, when they are still quite massive and self-gravitating, can survive.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

TW Hya: an old protoplanetary disc revived by its planet

Dark rings with bright rims are the indirect signposts of planets embedded in protoplanetary discs. In a recent first, an azimuthally elongated AU-scale blob, possibly a planet, was resolved with ALMA in TW Hya. The blob is at the edge of a cliff-like rollover in the dust disc rather than inside a dark ring. Here we build time-dependent models of TW Hya disc. We find that the classical paradigm cannot account for the morphology of the disc and the blob. We propose that ALMA-discovered blob hides a Neptune mass planet losing gas and dust. We show that radial drift of mm-sized dust particles naturally explains why the blob is located on the edge of the dust disc. Dust particles leaving the planet perform a characteristic U-turn relative to it, producing an azimuthally elongated blob-like emission feature. This scenario also explains why a 10 Myr old disc is so bright in dust continuum. Two scenarios for the dust-losing planet are presented. In the first, a dusty pre-runaway gas envelope of about 40 Earth mass Core Accretion planet is disrupted, e.g., as a result of a catastrophic encounter. In the second, a massive dusty pre-collapse gas giant planet formed by Gravitational Instability is disrupted by the energy released in its massive core. Future modelling may discriminate between these scenarios and allow us to study planet formation in an entirely new way -- by analysing the flows of dust and gas recently belonging to planets, informing us about the structure of pre-disruption planetary envelopes.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Grand challenges in protoplanetary disc modelling

The Protoplanetary Discussions conference --- held in Edinburgh, UK, from 7th --11th March 2016 --- included several open sessions led by participants. This paper reports on the discussions collectively concerned with the multiphysics modelling of protoplanetary discs, including the self-consistent calculation of gas and dust dynamics, radiative transfer and chemistry. After a short introduction to each of these disciplines in isolation, we identify a series of burning questions and grand challenges associated with their continuing development and integration. We then discuss potential pathways towards solving these challenges, grouped by strategical, technical and collaborative developments. This paper is not intended to be a review, but rather to motivate and direct future research and collaboration across typically distinct fields based on \textit{community driven input}, to encourage further progress in our understanding of circumstellar and protoplanetary discs.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Planet Formation Imager (PFI): science vision and key requirements

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On the gap-opening criterion of migrating planets in protoplanetary disks

We perform two-dimensional hydrodynamical simulations to quantitatively explore the torque balance criterion for gap-opening (as formulated by Crida et al. 2006) in a variety of disks when considering a migrating planet. We find that even when the criterion is satisfied, there are instances when planets still do not open gaps. We stress that gap-opening is not only dependent on whether a planet has the ability to open a gap, but whether it can do so quickly enough. This can be expressed as an additional condition on the gap-opening timescale versus the crossing time, i.e. the time it takes the planet to cross the region which it is carving out. While this point has been briefly made in the previous literature, our results quantify it for a range of protoplanetary disk properties and planetary masses, demonstrating how crucial it is for gap-opening. This additional condition has important implications for the survival of planets formed by core accretion in low mass disks as well as giant planets or brown dwarfs formed by gravitational instability in massive disks. It is particularly important for planets with intermediate masses susceptible to Type III-like migration. For some observed transition disks or disks with gaps, we expect that estimates on the potential planet masses based on the torque balance gap-opening criterion alone may not be sufficient. With consideration of this additional timescale criterion theoretical studies may find a reduced planet survivability or that planets may migrate further inwards before opening a gap.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Triggered fragmentation in self-gravitating discs: forming fragments at small radii

We carry out three dimensional radiation hydrodynamical simulations of gravitationally unstable discs to explore the movement of mass in a disc following its initial fragmentation. We find that the radial velocity of the gas in some parts of the disc increases by up to a factor of approximately 10 after the disc fragments, compared to before. While the movement of mass occurs in both the inward and outward directions, the inwards movement can cause the inner spirals of a self-gravitating disc to become sufficiently dense such that they can potentially fragment. This suggests that the dynamical behaviour of fragmented discs may cause subsequent fragmentation to occur at smaller radii than initially expected, but only after an initial fragment has formed in the outer disc.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Discovery of a Companion Candidate in the HD169142 Transition Disk and the Possibility of Multiple Planet Formation

We present L' and J-band high-contrast observations of HD169142, obtained with the VLT/NACO AGPM vector vortex coronagraph and the Gemini Planet Imager, respectively. A source located at 0".156+/-0".032 north of the host star (PA=7.4+/-11.3 degrees) appears in the final reduced L' image. At the distance of the star (~145 pc), this angular separation corresponds to a physical separation of 22.7+/-4.7 AU, locating the source within the recently resolved inner cavity of the transition disk. The source has a brightness of L'=12.2+/-0.5 mag, whereas it is not detected in the J band (J>13.8 mag). If its L' brightness arose solely from the photosphere of a companion and given the J-L' color constraints, it would correspond to a 28-32 MJupiter object at the age of the star, according to the COND models. Ongoing accretion activity of the star suggests, however, that gas is left in the inner disk cavity from which the companion could also be accreting. In this case the object could be lower in mass and its luminosity enhanced by the accretion process and by a circumplanetary disk. A lower mass object is more consistent with the observed cavity width. Finally, the observations enable us to place an upper limit on the L'-band flux of a second companion candidate orbiting in the disk annular gap at ~50 AU, as suggested by millimeter observations. If the second companion is also confirmed, HD169142 might be forming a planetary system, with at least two companions opening gaps and possibly interacting with each other.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Long-lasting dust rings in gas-rich disks: sculpting by single and multiple planets

We propose a mechanism by which dust rings in protoplanetary disks can form and be long-lasting compared to gas rings. This involves the existence of a pressure maximum which traps dust either in between two gap-opening planets or at the outermost gap edge of a single or multiple planet system, combined with the decoupling of large dust particles from the gas. We perform 2D gas hydrodynamical simulations of disks with one and two giant planets which may open deep or partial gaps. A gas ring forms in between two planets such that the surface mass density is higher than on either side of it. This ring is a region of pressure maximum where we expect large grains, which are marginally coupled to the gas and would otherwise be subject to radial drift, to collect. Such a pressure maximum also occurs at the outermost gap edge in a disk with one or more planets. We infer the dust evolution in these regions as the gas disk evolves, to understand the longer term behavior of the resulting dust rings. Over time the gas surface density in the ring(s) decreases, which may cause the larger trapped particles to decouple. Consequently, these particles are expected to stay in ring structure(s) longer than the gas. For a Minimum Mass Solar Nebula model, we expect that millimeter and centimeter-sized grains in the outer O(10) au would most likely undergo this trapping and decoupling process.

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

Growth and fragmentation of centimetre-sized dust aggregates: the dependence on aggregate size and porosity

We carry out three-dimensional Smoothed Particle Hydrodynamics simulations of spherical homogeneous SiO2 dust aggregates to investigate how the mass and the porosity of the aggregates affects their ability to survive an impact at various different collision velocities (between 1 - 27.5m/s). We explore how the threshold velocities for fragmentation vary with these parameters. Crucially, we find that the porosity plays a part of utmost importance in determining the outcome of collisions. In particular, we find that aggregates with filling factors >37% are significantly weakened and that the velocity regime in which the aggregates grow is reduced or even non-existent (instead, the aggregates either rebound off each other or break apart). At filling factors less than ~37% we find that more porous objects are weaker but not as weak as highly compact objects with filling factors >37%. In addition we find that (for a given aggregate density) collisions between very different mass objects have higher threshold velocities than those between very similar mass objects. We find that fragmentation velocities are higher than the typical values of 1m/s and that growth can even occur for velocities as high as 27.5m/s. Therefore, while the growth of aggregates is more likely if collisions between different sized objects occurs or if the aggregates are porous with filling factor <37%, it may also be hindered if the aggregates become too compact.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Growth of grains in Brown Dwarf disks

We perform coagulation and fragmentation simulations using the new physically-motivated model by Garaud et al. (2013) to determine growth locally in brown dwarf disks. We show that large grains can grow and that if brown dwarf disks are scaled down versions of T Tauri disks (in terms of stellar mass, disk mass and disk radius) growth at an equivalent location with respect to the disk truncation radius can occur to the same size in both disks. We show that similar growth occurs because the collisional timescales in the two disks are comparable. Our model may therefore potentially explain the recent observations of grain growth to millimetre sizes in brown dwarf disks, as seen in T Tauri disks.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Large grains can grow in circumstellar discs

We perform coagulation & fragmentation simulations to understand grain growth in T Tauri & brown dwarf discs. We present a physically-motivated approach using a probability distribution function for the collision velocities and separating the deterministic & stochastic velocities. We find growth to larger sizes compared to other models. Furthermore, if brown dwarf discs are scaled-down versions of T Tauri discs (in terms of stellar & disc mass, and disc radius), growth at the same location with respect to the outer edge occurs to similar sizes in both discs.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Small vs large dust grains in transitional disks: do different cavity sizes indicate a planet?

Transitional disks represent a short stage of the evolution of circumstellar material. Studies of dust grains in these objects can provide pivotal information on the mechanisms of planet formation. Dissimilarities in the spatial distribution of small (micron-size) and large (millimeter-size) dust grains have recently been pointed out. Constraints on the small dust grains can be obtained by imaging the distribution of scattered light at near-infrared wavelengths. We aim at resolving structures in the surface layer of transitional disks (with particular emphasis on the inner 10 - 50 AU), thus increasing the scarce sample of high resolution images of these objects. We obtained VLT/NACO near-IR high-resolution polarimetric differential imaging observations of SAO 206462 (HD135344B). This technique allows one to image the polarized scattered light from the disk without any occulting mask and to reach an inner working angle of 0.1''. A face-on disk is detected in H and Ks bands between 0.1'' and 0.9''. No significant differences are seen between the H and Ks images. In addition to the spiral arms, these new data allow us to resolve for the first time an inner cavity for small dust grains. The cavity size (about 28 AU) is much smaller than what is inferred for large dust grains from (sub)mm observations (39 to 50 AU). The interaction between the disk and potential orbiting companion(s) can explain both the spiral arm structure and the discrepant cavity sizes for small and large dust grains. One planet may be carving out the gas (and, thus, the small grains) at 28 AU, and generating a pressure bump at larger radii (39 AU), which holds back the large grains. We analytically estimate that, in this scenario, a single giant planet (with a mass between 5 and 15 Jupiter masses) at 17 to 20 AU from the star is consistent with the observed cavity sizes.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Triggered fragmentation in gravitationally unstable discs: forming fragments at small radii

We carry out three dimensional radiation hydrodynamical simulations of gravitationally unstable discs using to explore the movement of mass in a disc following its fragmentation. Compared to a more quiescent state before it fragments, the radial velocity of the gas increases by up to a factor of ~2-3 after fragmentation. While the mass movement occurs both inwards and outwards, the inwards motion can cause the inner spirals to be sufficiently dense that they may become unstable and potentially fragment. Consequently, the dynamical behaviour of fragmented discs may cause subsequent fragmentation at smaller radii after an initial fragment has formed in the outer disc.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

From dust to planetesimals: an improved model for collisional growth in protoplanetary disks

Planet formation occurs within the gas and dust rich environments of protoplanetary disks. Observations of these objects show that the growth of primordial sub micron sized particles into larger aggregates occurs at the earliest stages of the disks. However, theoretical models of particle growth that use the Smoluchowski equation to describe collisional coagulation and fragmentation have so far failed to produce large particles while maintaining a significant populations of small grains. This has been generally attributed to the existence of two barriers impeding growth due to bouncing and fragmentation of colliding particles. In this paper, we demonstrate that the importance of these barriers has been artificially inflated through the use of simplified models that do not take into account the stochastic nature of the particle motions within the gas disk. We present a new approach in which the relative velocities between two particles is described by a probability distribution function that models both deterministic motion and stochastic motion. Taking both into account can give quite different results to what has been considered recently in other studies. We demonstrate the vital effect of two "ingredients" for particle growth: the proper implementation of a velocity distribution function that overcomes the bouncing barrier and, in combination with mass transfer in high-mass-ratio collisions, boosts the growth of larger particles beyond the fragmentation barrier. A robust result of our simulations is the emergence of two particle populations (small and large), potentially explaining simultaneously a number of long-standing problems in protoplanetary disks, including planetesimal formation close to the central star, the presence of mm to cm size particles far out in the disk, and the persistence of micron-size grains for millions of years.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the convergence of the critical cooling timescale for the fragmentation of self-gravitating discs

We carry out simulations of gravitationally unstable discs using a Smoothed Particle Hydrodynamics (SPH) code and a grid-based hydrodynamics code, FARGO, to understand the previous non-convergent results reported by Meru & Bate (2011a). We obtain evidence that convergence with increasing resolution occurs with both SPH and FARGO and in both cases we find that the critical cooling timescale is larger than previously thought. We show that SPH has a first-order convergence rate while FARGO converges with a second-order rate. We show that the convergence of the critical cooling timescale for fragmentation depends largely on the numerical viscosity employed in both SPH and FARGO. With SPH, particle velocity dispersion may also play a role. We show that reducing the dissipation from the numerical viscosity leads to larger values of the critical cooling time at a given resolution. For SPH, we find that the effect of the dissipation due to the numerical viscosity is somewhat larger than had previously been appreciated. In particular, we show that using a quadratic term in the SPH artificial viscosity (beta_{SPH}) that is too low appears to lead to excess dissipation in gravitationally unstable discs, which may affect any results that sensitively depend on the thermodynamics, such as disc fragmentation. We show that the two codes converge to values of the critical cooling timescale, beta_{crit} > 20 (for a ratio of specific heats of gamma=5/3), and perhaps even as large as beta_{crit} \approx 30. These are approximately 3-5 times larger than has been found by most previous studies. This is equivalent to a maximum gravitational stress that a disc can withstand without fragmenting of alpha_{GI,crit} \approx 0.013-0.02, which is much smaller than the values typically used in the literature. It is therefore easier for self-gravitating discs to fragment than has been concluded from most past studies.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Numerical convergence in self-gravitating disc simulations: initial conditions and edge effects

We study the numerical convergence of hydrodynamical simulations of self-gravitating accretion discs, in which a simple cooling law is balanced by shock heating. It is well-known that there exists a critical cooling time scale for which shock heating can no longer compensate for the energy losses, at which point the disc fragments. The numerical convergence of previous results of this critical cooling time scale was questioned recently using Smoothed Particle Hydrodynamics (SPH). We employ a two-dimensional grid-based code to study this problem, and find that for smooth initial conditions, fragmentation is possible for slower cooling as the resolution is increased, in agreement with recent SPH results. We show that this non-convergence is at least partly due to the creation of a special location in the disc, the boundary between the turbulent and the laminar region, when cooling towards a gravito-turbulent state. Converged results appear to be obtained in setups where no such sharp edges appear, and we then find a critical cooling time scale of ~ 4 $Ω^{-1}$, where $Ω$ is the local angular velocity.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Rapid inward migration of planets formed by gravitational instability

The observation of massive exoplanets at large separation from their host star, like in the HR 8799 system, challenges theories of planet formation. A possible formation mechanism involves the fragmentation of massive self-gravitating discs into clumps. While the conditions for fragmentation have been extensively studied, little is known of the subsequent evolution of these giant planet embryos, in particular their expected orbital migration. Assuming a single planet has formed by fragmentation, we investigate its interaction with the gravitoturbulent disc it is embedded in. Two-dimensional hydrodynamical simulations are used with a simple prescription for the disc cooling. A steady gravitoturbulent disc is first set up, after which simulations are restarted including a planet with a range of masses approximately equal to the clump's initial mass expected in fragmenting discs. Planets rapidly migrate inwards, despite the stochastic kicks due to the turbulent density fluctuations. We show that the migration timescale is essentially that of type I migration, with the planets having no time to open a gap. In discs with aspect ratio ~ 0.1 at their forming location, planets with a mass comparable to, or larger than Jupiter's can migrate in as short as 10000 years, that is, about 10 orbits at 100 AU. Massive planets formed at large separation from their star by gravitational instability are thus unlikely to stay in place, and should rapidly migrate towards the inner parts of protoplanetary discs, regardless of the planet mass.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The four-populations model: a new classification scheme for pre-planetesimal collisions

Within the collision growth scenario for planetesimal formation, the growth step from centimetre sized pre-planetesimals to kilometre sized planetesimals is still unclear. The formation of larger objects from the highly porous pre-planetesimals may be halted by a combination of fragmentation in disruptive collisions and mutual rebound with compaction. However, the right amount of fragmentation is necessary to explain the observed dust features in late T Tauri discs. Therefore, detailed data on the outcome of pre-planetesimal collisions is required and has to be presented in a suitable and precise format. We propose and apply a new classification scheme for pre-planetesimal collisions based on the quantitative aspects of four fragment populations: the largest and second largest fragment, a power-law population, and a sub-resolution population. For the simulations of pre-planetesimal collisions, we adopt the SPH numerical scheme with extensions for the simulation of porous solid bodies. By means of laboratory benchmark experiments, this model was previously calibrated and tested for the correct simulation of the compaction, bouncing, and fragmentation behaviour of macroscopic highly porous silica dust aggregates. It is shown that previous attempts to map collision data were much too oriented on qualitatively categorising into sticking, bouncing, and fragmentation events. We show that the four-populations model encompasses all previous categorisations and in addition allows for transitions. This is because it is based on quantitative characteristic attributes of each population such as the mass, kinetic energy, and filling factor. As a demonstration of the applicability and the power of the four-populations model, we utilise it to present the results of a study on the influence of collision velocity in head-on collisions of intermediate porosity aggregates.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Exploring the conditions required to form giant planets via gravitational instability in massive protoplanetary discs

We carry out global three-dimensional radiation hydrodynamical simulations of self-gravitating accretion discs to determine if, and under what conditions, a disc may fragment to form giant planets. We explore the parameter space (in terms of the disc opacity, temperature and size) and include the effect of stellar irradiation. We find that the disc opacity plays a vital role in determining whether a disc fragments. Specifically, opacities that are smaller than interstellar Rosseland mean values promote fragmentation (even at small radii, R < 25AU) since low opacities allow a disc to cool quickly. This may occur if a disc has a low metallicity or if grain growth has occurred. With specific reference to the HR 8799 planetary system, given its star is metal-poor, our results suggest that the formation of its imaged planetary system could potentially have occurred by gravitational instability. We also find that the presence of stellar irradiation generally acts to inhibit fragmentation (since the discs can only cool to the temperature defined by stellar irradiation). However, fragmentation may occur if the irradiation is sufficiently weak that it allows the disc to attain a low Toomre stability parameter.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-convergence of the critical cooling timescale for fragmentation of self-gravitating discs

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results suggest that at the very least, the critical value of the cooling timescale is longer than previously thought. However, the absence of convergence also raises the question of whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Non-convergence of the critical cooling timescale for fragmentation of self-gravitating discs

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results certainly suggest that beta_{crit} is larger than previously thought. However, the absence of convergence also questions whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

On the fragmentation criteria of self-gravitating protoplanetary discs

We investigate the fragmentation criterion in massive self-gravitating discs. We present new analysis of the fragmentation conditions which we test by carrying out global three-dimensional numerical simulations. Whilst previous work has placed emphasis on the cooling timescale in units of the orbital timescale, β, we find that at a given radius the surface mass density (i.e. disc mass and profile) and star mass also play a crucial role in determining whether a disc fragments or not as well as where in the disc fragments form. We find that for shallow surface mass density profiles (p<2, where Σ\propto R^{-p}), fragments form in the outer regions of the disc. However for steep surface mass density profiles (p \gtrsim 2), fragments form in the inner regions of a disc. In addition, we also find that the critical value of the cooling timescale in units of the orbital timescale found in previous simulations is only applicable to certain disc surface mass density profiles and for particular disc radii and is not a general rule for all discs. We find an empirical fragmentation criteria between the cooling timescale in units of the orbital timescale, β, the surface mass density, the star mass and the radius.

0 citations0 reviewsNo code linkedOpen access