Showing 1–50 of 53 results for author: Meru, F

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

The orbit of HD 142527 B is too compact to explain many of the disc features

Authors: M. Nowak, S. Rowther, S. Lacour, F. Meru, R. Nealon, D. J. Price

Abstract: HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13… ▽ More HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13 $\pm$ 0.03 $M_\odot{}$), which is thought to be responsible for most of the structures observed in the surrounding disc. We gathered VLTI/GRAVITY observations of HD 142527, either from our own programmes or from the ESO archive. We used this inhomogeneous set of data to extract a total of seven high-precision measurements of the relative astrometry between HD 142527 A and B, spread from mid-2017 to early 2021. Combined with what is available in the literature, we now have 9 yr of astrometric monitoring on HD 142527. We used orbit fitting tools to determine the orbital parameters of HD 142527 B, and used them as inputs for a 3D hydrodynamical model of the disc to determine whether or not the binary is able to create the structures observed in the disc. Our VLTI/GRAVITY astrometry gives excellent constraints on the orbit of HD 142527 B. We show that the secondary is following an orbit of semi-major axis a = 10.80 $\pm$ 0.22 au, with moderate eccentricity (e = 0.47 $\pm$ 0.01). With such a compact orbit, we show that HD 142527 B can only generate a gap and spiral arms of $\sim$30 au in the disc, which is much smaller than what is revealed by observations. Even from a theoretical standpoint, the observed cavity size of $\sim$100 au far exceeds even the most generous predictions for a companion like HD 142527 B on such a compact orbit. Thus, we conclude that the low-mass companion cannot be solely responsible for the observed morphology of the disc surrounding the system. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 11 pages, 6 figures. Accepted for publication in A&A

The Role of Drag and Gravity on Dust Concentration in a Gravitationally Unstable Disc

Authors: Sahl Rowther, Rebecca Nealon, Farzana Meru, James Wurster, Hossam Aly, Richard Alexander, Ken Rice, Richard A. Booth

Abstract: We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If w… ▽ More We carry out three dimensional smoothed particle hydrodynamics simulations to study the role of gravitational and drag forces on the concentration of large dust grains (St > 1) in the spiral arms of gravitationally unstable protoplanetary discs, and the resulting implications for planet formation. We find that both drag and gravity play an important role in the evolution of large dust grains. If we include both, grains that would otherwise be partially decoupled will become well coupled and trace the spirals. For the dust grains most influenced by drag (with Stokes numbers near unity), the dust disc quickly becomes gravitationally unstable and rapidly forms clumps with masses between 0.15 - 6 Earth masses. A large fraction of clumps are below the threshold where runaway gas accretion can occur. However, if dust self-gravity is neglected, the dust is unable to form clumps, despite still becoming trapped in the gas spirals. When large dust grains are unable to feel either gas gravity or drag, the dust is unable to trace the gas spirals. Hence, full physics is needed to properly simulate dust in gravitationally unstable discs. Dust trapping of large grains in spiral arms of discs stable to gas fragmentation could explain planet formation in very young discs by a population of planetesimals formed due to the combined roles of drag and gravity in the earliest stages of a disc's evolution. Furthermore, it highlights that gravitationally unstable discs are not just important for forming gas giants quickly, it can also rapidly form Earth mass bodies. △ Less

Submitted 24 January, 2024; v1 submitted 17 January, 2024; originally announced January 2024.

Comments: Accepted for publication in MNRAS. Updated references

Sandwiched planet formation: restricting the mass of a middle planet

Authors: Matthew Pritchard, Farzana Meru, Sahl Rowther, David Armstrong, Kaleb Randall

Abstract: We conduct gas and dust hydrodynamical simulations of protoplanetary discs with one and two embedded planets to determine the impact that a second planet located further out in the disc has on the potential for subsequent planet formation in the region locally exterior to the inner planet. We show how the presence of a second planet has a strong influence on the collection of solid material near t… ▽ More We conduct gas and dust hydrodynamical simulations of protoplanetary discs with one and two embedded planets to determine the impact that a second planet located further out in the disc has on the potential for subsequent planet formation in the region locally exterior to the inner planet. We show how the presence of a second planet has a strong influence on the collection of solid material near the inner planet, particularly when the outer planet is massive enough to generate a maximum in the disc's pressure profile. This effect in general acts to reduce the amount of material that can collect in a pressure bump generated by the inner planet. When viewing the inner pressure bump as a location for potential subsequent planet formation of a third planet, we therefore expect that the mass of such a planet will be smaller than it would be in the case without the outer planet, resulting in a small planet being sandwiched between its neighbours - this is in contrast to the expected trend of increasing planet mass with radial distance from the host star. We show that several planetary systems have been observed that do not show this trend but instead have a smaller planet sandwiched in between two more massive planets. We present the idea that such an architecture could be the result of the subsequent formation of a middle planet after its two neighbours formed at some earlier stage. △ Less

Submitted 9 October, 2023; v1 submitted 2 October, 2023; originally announced October 2023.

Comments: 13 pages, 10 figures, accepted by MNRAS (in press)

Continuing to Hide Signatures of Gravitational Instability in Protoplanetary Discs with Planets

Authors: Sahl Rowther, Rebecca Nealon, Farzana Meru

Abstract: We carry out three dimensional smoothed particle hydrodynamics simulations to study the impact of planet-disc interactions on a gravitationally unstable protoplanetary disc. We find that the impact of a planet on the disc's evolution can be described by three scenarios. If the planet is sufficiently massive, the spiral wakes generated by the planet dominate the evolution of the disc and gravitatio… ▽ More We carry out three dimensional smoothed particle hydrodynamics simulations to study the impact of planet-disc interactions on a gravitationally unstable protoplanetary disc. We find that the impact of a planet on the disc's evolution can be described by three scenarios. If the planet is sufficiently massive, the spiral wakes generated by the planet dominate the evolution of the disc and gravitational instabilities are completely suppressed. If the planet's mass is too small, then gravitational instabilities are unaffected. If the planet's mass lies between these extremes, gravitational instabilities are weakened. We present mock Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations showing that the observability of large-scale spiral structures is diminished or completely suppressed when the planet is massive enough to influence the disc's evolution. Our results show that massive discs that would be expected to be gravitationally unstable can appear axisymmetric in the presence of a planet. Thus, the absence of observed large-scale spiral structures alone is not enough to place upper limits on the disc's mass, which could have implications on observations of young Class I discs with rings & gaps. △ Less

Submitted 31 October, 2022; originally announced October 2022.

Comments: Accepted for publication in MNRAS. 11 pages, 9 figures

Turbulent Disk Viscosity and the Bifurcation of Planet Formation Histories

Authors: Jessica Speedie, Ralph E. Pudritz, Alex J. Cridland, Farzana Meru, Richard A. Booth

Abstract: 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 thro… ▽ More 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. △ Less

Submitted 4 February, 2022; v1 submitted 2 November, 2021; originally announced November 2021.

Comments: 25 pages, 15 figures, 10 movies (links embedded). Published in MNRAS

Journal ref: MNRAS, Volume 510, Issue 4, March 2022, Pages 6059-6084

Warping Away Gravitational Instabilities in Protoplanetary Discs

Authors: Sahl Rowther, Rebecca Nealon, Farzana Meru

Abstract: 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-axisy… ▽ More 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. △ Less

Submitted 5 December, 2021; v1 submitted 12 October, 2021; originally announced October 2021.

Comments: Accepted for publication in ApJ

Hiding Signatures of Gravitational Instability in Protoplanetary Discs with Planets

Authors: Sahl Rowther, Farzana Meru, Grant M. Kennedy, Rebecca Nealon, Christophe Pinte

Abstract: We carry out three dimensional SPH simulations to show that a migrating giant planet strongly suppresses the spiral structure in self-gravitating discs. We present mock ALMA continuum observations which show that in the absence of a planet, spiral arms due to gravitational instability are easily observed. Whereas in the presence of a giant planet, the spiral structures are suppressed by the migrat… ▽ More We carry out three dimensional SPH simulations to show that a migrating giant planet strongly suppresses the spiral structure in self-gravitating discs. We present mock ALMA continuum observations which show that in the absence of a planet, spiral arms due to gravitational instability are easily observed. Whereas in the presence of a giant planet, the spiral structures are suppressed by the migrating planet resulting in a largely axisymmetric disc with a ring and gap structure. Our modelling of the gas kinematics shows that the planet's presence could be inferred, for example, using optically thin 13C16O. Our results show that it is not necessary to limit the gas mass of discs by assuming high dust-to-gas mass ratios in order to explain a lack of spiral features that would otherwise be expected in high mass discs. △ Less

Submitted 9 November, 2020; originally announced November 2020.

Comments: Accepted to ApJL. 9 pages, 5 figures

Planet Migration in Self-Gravitating Discs: Survival of Planets

Authors: Sahl Rowther, Farzana Meru

Abstract: 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… ▽ More 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. △ Less

Submitted 4 June, 2020; originally announced June 2020.

Comments: 13 pages, 14 figures, accepted for publication in MNRAS

TW Hya: an old protoplanetary disc revived by its planet

Authors: Sergei Nayakshin, Takashi Tsukagoshi, Cassandra Hall, Allona Vazan, Ravit Helled, Jack Humphries, Farzana Meru, Patrick Neunteufel, Olja Panic

Abstract: 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 ca… ▽ More 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. △ Less

Submitted 14 May, 2020; v1 submitted 21 April, 2020; originally announced April 2020.

Comments: Typos fixed, references and author list updated

A remnant planetary core in the hot-Neptune desert

Authors: David J. Armstrong, Théo A. Lopez, Vardan Adibekyan, Richard A. Booth, Edward M. Bryant, Karen A. Collins, Alexandre Emsenhuber, Chelsea X. Huang, George W. King, Jorge Lillo-box, Jack J. Lissauer, Elisabeth C. Matthews, Olivier Mousis, Louise D. Nielsen, Hugh Osborn, Jon Otegi, Nuno C. Santos, Sérgio G. Sousa, Keivan G. Stassun, Dimitri Veras, Carl Ziegler, Jack S. Acton, Jose M. Almenara, David R. Anderson, David Barrado , et al. (69 additional authors not shown)

Abstract: 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 s… ▽ More 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. △ Less

Submitted 16 July, 2020; v1 submitted 23 March, 2020; originally announced March 2020.

Comments: Published in Nature. This is a preprint of the article, before minor changes made during the refereeing and editing process. The published PDF is at https://www.nature.com/articles/s41586-020-2421-7 and can be accessed for free by following this link: https://rdcu.be/b5miB . Abstract updated to match published version

Mass determinations of the three mini-Neptunes transiting TOI-125

Authors: L. D. Nielsen, D. Gandolfi, D. J. Armstrong, J. S. Jenkins, M. Fridlund, N. C. Santos, F. Dai, V. Adibekyan, R. Luque, J. H. Steffen, M. Esposito, F. Meru, S. Sabotta, E. Bolmont, D. Kossakowski, J. F. Otegi, F. Murgas, M. Stalport, F. ~Rodler, M. R. Díaz, N. T. ~Kurtovic, G. Ricker, R. Vanderspek, D. W. Latham, S. Seager , et al. (55 additional authors not shown)

Abstract: The Transiting Exoplanet Survey Satellite, TESS, is currently carrying out an all-sky search for small planets transiting bright stars. In the first year of the TESS survey, steady progress was made in achieving the mission's primary science goal of establishing bulk densities for 50 planets smaller than Neptune. During that year, TESS's observations were focused on the southern ecliptic hemispher… ▽ More The Transiting Exoplanet Survey Satellite, TESS, is currently carrying out an all-sky search for small planets transiting bright stars. In the first year of the TESS survey, steady progress was made in achieving the mission's primary science goal of establishing bulk densities for 50 planets smaller than Neptune. During that year, TESS's observations were focused on the southern ecliptic hemisphere, resulting in the discovery of three mini-Neptunes orbiting the star TOI-125, a V=11.0 K0 dwarf. We present intensive HARPS radial velocity observations, yielding precise mass measurements for TOI-125b, TOI-125c and TOI-125d. TOI-125b has an orbital period of 4.65 days, a radius of $2.726 \pm 0.075 ~\mathrm{R_{\rm E}}$, a mass of $ 9.50 \pm 0.88 ~\mathrm{M_{\rm E}}$ and is near the 2:1 mean motion resonance with TOI-125c at 9.15 days. TOI-125c has a similar radius of $2.759 \pm 0.10 ~\mathrm{R_{\rm E}}$ and a mass of $ 6.63 \pm 0.99 ~\mathrm{M_{\rm E}}$, being the puffiest of the three planets. TOI-125d, has an orbital period of 19.98 days and a radius of $2.93 \pm 0.17~\mathrm{R_{\rm E}}$ and mass $13.6 \pm 1.2 ~\mathrm{M_{\rm E}}$. For TOI-125b and TOI-125d we find unusual high eccentricities of $0.19\pm 0.04$ and $0.17^{+0.08}_{-0.06}$, respectively. Our analysis also provides upper mass limits for the two low-SNR planet candidates in the system; for TOI-125.04 ($R_P=1.36 ~\mathrm{R_{\rm E}}$, $P=$0.53 days) we find a $2σ$ upper mass limit of $1.6~\mathrm{M_{\rm E}}$, whereas TOI-125.05 ( $R_P=4.2^{+2.4}_{-1.4} ~\mathrm{R_{\rm E}}$, $P=$ 13.28 days) is unlikely a viable planet candidate with upper mass limit $2.7~\mathrm{M_{\rm E}}$. We discuss the internal structure of the three confirmed planets, as well as dynamical stability and system architecture for this intriguing exoplanet system. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS

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

Authors: Dimitri Veras, Pier-Emmanuel Tremblay, J. J. Hermes, Catriona H. McDonald, Grant M. Kennedy, Farzana Meru, Boris T. Gänsicke

Abstract: 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… ▽ More 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. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: Accepted for publication in MNRAS

Massive discs around low-mass stars

Authors: Thomas J. Haworth, James Cadman, Farzana Meru, Cassandra Hall, Emma Albertini, Duncan Forgan, Ken Rice, James E. Owen

Abstract: 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. T… ▽ More 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. △ Less

Submitted 27 March, 2020; v1 submitted 17 January, 2020; originally announced January 2020.

Comments: 21 pages, Accepted for publication MNRAS

A Gap in the Mass Distribution for Warm Neptune and Terrestrial Planets

Authors: David J. Armstrong, Farzana Meru, Daniel Bayliss, Grant M. Kennedy, Dimitri Veras

Abstract: Structure in the planet distribution provides an insight into the processes that shape the formation and evolution of planets. The Kepler mission has led to an abundance of statistical discoveries in regards to planetary radius, but the number of observed planets with measured masses is much smaller. By incorporating results from recent mass determination programs, we have discovered a new gap eme… ▽ More Structure in the planet distribution provides an insight into the processes that shape the formation and evolution of planets. The Kepler mission has led to an abundance of statistical discoveries in regards to planetary radius, but the number of observed planets with measured masses is much smaller. By incorporating results from recent mass determination programs, we have discovered a new gap emerging in the planet population for sub-Neptune mass planets with orbital periods less than 20 days. The gap follows a slope of decreasing mass with increasing orbital period, has a width of a few $M_\oplus$, and is potentially completely devoid of planets. Fitting gaussian mixture models to the planet population in this region favours a bimodel distribution over a unimodel one with a reduction in Bayesian Information Criterion (BIC) of 19.9, highlighting the gap significance. We discuss several processes which could generate such a feature in the planet distribution, including a pileup of planets above the gap region, tidal interactions with the host star, dynamical interactions with the disk, with other planets, or with accreting material during the formation process. △ Less

Submitted 27 June, 2019; originally announced June 2019.

Comments: 8 pages, 4 figures, accepted for publication in ApJ Letters

A planetesimal orbiting within the debris disc around a white dwarf star

Authors: Christopher J. Manser, Boris T. Gänsicke, Siegfried Eggl, Mark Hollands, Paula Izquierdo, Detlev Koester, John D. Landstreet, Wladimir Lyra, Thomas R. Marsh, Farzana Meru, Alexander J. Mustill, Pablo Rodríguez-Gil, Odette Toloza, Dimitri Veras, David J. Wilson, Matthew R. Burleigh, Melvyn B. Davies, Jay Farihi, Nicola Gentile Fusillo, Domitilla de Martino, Steven G. Parsons, Andreas Quirrenbach, Roberto Raddi, Sabine Reffert, Melania Del Santo , et al. (7 additional authors not shown)

Abstract: Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. W… ▽ More Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength. △ Less

Submitted 3 April, 2019; originally announced April 2019.

Comments: 11 Figures, 3 Tables

Revealing signatures of planets migrating in protoplanetary discs with ALMA multi-wavelength observations

Authors: Pooneh Nazari, Richard A. Booth, Cathie J. Clarke, Giovanni P. Rosotti, Marco Tazzari, Attila Juhasz, Farzana Meru

Abstract: Recent observations show that rings and gaps are ubiquitous in protoplanetary discs. These features are often interpreted as being due to the presence of planets; however, the effect of planetary migration on the observed morphology has not been investigated hitherto. In this work we investigate whether multiwavelength mm/submm observations can detect signatures of planet migration, using 2D dusty… ▽ More Recent observations show that rings and gaps are ubiquitous in protoplanetary discs. These features are often interpreted as being due to the presence of planets; however, the effect of planetary migration on the observed morphology has not been investigated hitherto. In this work we investigate whether multiwavelength mm/submm observations can detect signatures of planet migration, using 2D dusty hydrodynamic simulations to model the structures generated by migrating planets and synthesising ALMA continuum observations at 0.85 and 3 mm. We identify three possible morphologies for a migrating planet: a slowly migrating planet is associated with a single ring outside the planet's orbit, a rapidly migrating planet is associated with a single ring inside the planet's orbit while a planet migrating at intermediate speed generates one ring on each side of the planet's orbit. We argue that multiwavelength data can distinguish multiple rings produced by a migrating planet from other scenarios for creating multiple rings, such as multiple planets or discs with low viscosity. The signature of migration is that the outer ring has a lower spectral index, due to larger dust grains being trapped there. Of the recent ALMA observations revealing protoplanetary discs with multiple rings and gaps, we suggest that Elias 24 is the best candidate for a planet migrating in the intermediate speed regime. △ Less

Submitted 21 March, 2019; v1 submitted 7 March, 2019; originally announced March 2019.

Comments: Accepted for publication in MNRAS

Giant planets and brown dwarfs on wide orbits: a code comparison project

Authors: Mark Fletcher, Sergei Nayakshin, Dimitris Stamatellos, Walter Dehnen, Farzana Meru, Lucio Mayer, Hongping Deng, Ken Rice

Abstract: Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion wit… ▽ More Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them. △ Less

Submitted 23 January, 2019; originally announced January 2019.

Comments: submitted to MNRAS. This is version 2 of the paper after considering referees comments and changes

HD219666b: A hot-Neptune from TESS Sector 1

Authors: M. Esposito, D. J. Armstrong, D. Gandolfi, V. Adibekyan, M. Fridlund, N. C. Santos, J. H. Livingston, E. Delgado Mena, L. Fossati, J. Lillo-Box, O. Barragán, D. Barrado, P. E. Cubillos, B. Cooke, A. B. Justesen, F. Meru, R. F. Díaz, F. Dai, L. D. Nielsen, C. M. Persson, P. J. Wheatley, A. P. Hatzes, V. Van Eylen, M. M. Musso, R. Alonso , et al. (51 additional authors not shown)

Abstract: We report on the confirmation and mass determination of a transiting planet orbiting the old and inactive G7 dwarf star HD219666 (Mstar = 0.92 +/- 0.03 MSun, Rstar = 1.03 +/- 0.03 RSun, tau_star = 10 +/- 2 Gyr). With a mass of Mb = 16.6 +/- 1.3 MEarth, a radius of Rb = 4.71 +/- 0.17 REarth, and an orbital period of P ~ 6 days, HD219666b is a new member of a rare class of exoplanets: the hot-Neptun… ▽ More We report on the confirmation and mass determination of a transiting planet orbiting the old and inactive G7 dwarf star HD219666 (Mstar = 0.92 +/- 0.03 MSun, Rstar = 1.03 +/- 0.03 RSun, tau_star = 10 +/- 2 Gyr). With a mass of Mb = 16.6 +/- 1.3 MEarth, a radius of Rb = 4.71 +/- 0.17 REarth, and an orbital period of P ~ 6 days, HD219666b is a new member of a rare class of exoplanets: the hot-Neptunes. The Transiting Exoplanet Survey Satellite (TESS) observed HD219666 (also known as TOI-118) in its Sector 1 and the light curve shows four transit-like events, equally spaced in time. We confirmed the planetary nature of the candidate by gathering precise radial-velocity measurements with HARPS@ESO3.6m. We used the co-added HARPS spectrum to derive the host star fundamental parameters (Teff = 5527 +/- 65 K, log g = 4.40 +/- 0.11 (cgs), [Fe/H]= 0.04 +/- 0.04 dex, log R'HK = -5.07 +/- 0.03), as well as the abundances of many volatile and refractory elements. The host star brightness (V = 9.9) makes it suitable for further characterisation by means of in-transit spectroscopy. The determination of the planet orbital obliquity, along with the atmospheric metal-to-hydrogen content and thermal structure could provide us with important clues on the formation mechanisms of this class of objects. △ Less

Submitted 18 February, 2019; v1 submitted 14 December, 2018; originally announced December 2018.

Comments: Accepted for publication in A&A

Journal ref: A&A 623, A165 (2019)

High-Resolution ALMA Observations of HD100546: Asymmetric Circumstellar Ring, and Circumplanetary Disk Upper Limits

Authors: J. E. Pineda, J. Szulágyi, S. P. Quanz, E. F. van Dishoeck, A. Garufi, F. Meru, G. D. Mulders, L. Testi, M. R. Meyer, M. Reggiani

Abstract: We present long baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870$\,μ$m dust continuum emission and CO (3-2) from the protoplanetary disk around the Herbig Ae/Be star HD 100546, which is one of the few systems claimed to have two young embedded planets. These observations achieve a resolution of 4 au (3.8 mas), an rms noise of 66$μ$Jy/beam, and reveal an asymmetr… ▽ More We present long baseline Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 870$\,μ$m dust continuum emission and CO (3-2) from the protoplanetary disk around the Herbig Ae/Be star HD 100546, which is one of the few systems claimed to have two young embedded planets. These observations achieve a resolution of 4 au (3.8 mas), an rms noise of 66$μ$Jy/beam, and reveal an asymmetric ring between $\sim$20-40 au with largely optically thin dust continuum emission. This ring is well fit by two concentric and overlapping Gaussian rings of different widths and a Vortex. In addition, an unresolved component is detected at a position consistent with the central star, which may trace the central inner disk ($<$2au in radius). We report a lack of compact continuum emission at the positions of both claimed protoplanets. We use this result to constrain the circumplanetary disk (CPD) mass and size of 1.44M$_{\rm Earth}$ and 0.44au in the optically thin and thick regime, respectively, for the case of the previously directly imaged protoplanet candidate at $\sim$55 au (HD100546 b). We compare these empirical CPD constraints to previous numerical simulations. This suggests that HD100546 b is inconsistent with several planet accretion models, while gas-starved models are also still compatible. We estimate the planetary mass as 1.65 M$_J$ by using the relation between planet, circumstellar, and circumplanetary masses derived from numerical simulations. Finally, the CO integrated intensity map shows a possible spiral arm feature that would match the spiral features identified in Near-Infrared scattered light polarized emission, which suggests a real spiral feature in the disk surface that needs to be confirmed with further observations. △ Less

Submitted 26 November, 2018; originally announced November 2018.

Comments: 14 pages, accepted by ApJ

Is the ring inside or outside the planet?: The effect of planet migration on dust rings

Authors: Farzana Meru, Giovanni P. Rosotti, Richard A. Booth, Pooneh Nazari, Cathie J. Clarke

Abstract: Planet migration in protoplanetary discs plays an important role in the longer term evolution of planetary systems, yet we currently have no direct observational test to determine if a planet is migrating in its gaseous disc. We explore the formation and evolution of dust rings - now commonly observed in protoplanetary discs by ALMA - in the presence of relatively low mass (12-60 Earth masses) mig… ▽ More Planet migration in protoplanetary discs plays an important role in the longer term evolution of planetary systems, yet we currently have no direct observational test to determine if a planet is migrating in its gaseous disc. We explore the formation and evolution of dust rings - now commonly observed in protoplanetary discs by ALMA - in the presence of relatively low mass (12-60 Earth masses) migrating planets. Through two dimensional hydrodynamical simulations using gas and dust we find that the importance of perturbations in the pressure profile interior and exterior to the planet varies for different particle sizes. For small sizes a dust enhancement occurs interior to the planet, whereas it is exterior to it for large particles. The transition between these two behaviours happens when the dust drift velocity is comparable to the planet migration velocity. We predict that an observational signature of a migrating planet consists of a significant outwards shift of an observed midplane dust ring as the wavelength is increased. △ Less

Submitted 15 October, 2018; originally announced October 2018.

Comments: 18 pages, 18 figures, accepted for publication in MNRAS

High resolution millimetre imaging of the CI Tau protoplanetary disc - a massive ensemble of protoplanets from 0.1 - 100 au

Authors: Cathie J. Clarke, Marco Tazzari, Attila Juhász, Giovanni Rosotti, Richard Booth, Stefano Facchini, John D. Ilee, Christopher M. Johns-Krull, Mihkel Kama, Farzana Meru, Lisa Prato

Abstract: We present high resolution millimeter continuum imaging of the disc surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disc system. The system has extended mm emission on which are superposed three prominent annular gaps at radii ~ 13, 39 and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including t… ▽ More We present high resolution millimeter continuum imaging of the disc surrounding the young star CI Tau, a system hosting the first hot Jupiter candidate in a protoplanetary disc system. The system has extended mm emission on which are superposed three prominent annular gaps at radii ~ 13, 39 and 100 au. We argue that these gaps are most likely to be generated by massive planets so that, including the hot Jupiter, the system contains four gas giant planets at an age of only 2 Myr. Two of the new planets are similarly located to those inferred in the famous HL Tau protoplanetary disc; in CI Tau, additional observational data enables a more complete analysis of the system properties than was possible for HL Tau. Our dust and gas dynamical modeling satisfies every available observational constraint and points to the most massive ensemble of exo-planets ever detected at this age, with its four planets spanning a factor 1000 in orbital radius.Our results show that the association between hot Jupiters and gas giants on wider orbits, observed in older stars, is apparently in place at an early evolutionary stage. △ Less

Submitted 21 September, 2018; originally announced September 2018.

Comments: Accepted on Astrophysical Journal Letters; 7 pages, 2 figures

Are Elias 2-27's spiral arms driven by self-gravity, or by a companion? A comparative spiral morphology study

Authors: Duncan H. Forgan, John D. Ilee, Farzana Meru

Abstract: The spiral waves detected in the protostellar disc surrounding Elias 2-27 have been suggested as evidence of the disc being gravitationally unstable. However, previous work has shown that a massive, stable disc undergoing an encounter with a massive companion are also consistent with the observations. We compare the spiral morphology of smoothed particle hydrodynamic simulations modelling both cas… ▽ More The spiral waves detected in the protostellar disc surrounding Elias 2-27 have been suggested as evidence of the disc being gravitationally unstable. However, previous work has shown that a massive, stable disc undergoing an encounter with a massive companion are also consistent with the observations. We compare the spiral morphology of smoothed particle hydrodynamic simulations modelling both cases. The gravitationally unstable disc produces symmetric, tightly wound spiral arms with constant pitch angle, as predicted by the literature. The companion disc's arms are asymmetric, with pitch angles that increase with radius. However, these arms are not well-fitted by standard analytic expressions, due to the high disc mass and relatively low companion mass. We note that differences (or indeed similarities) in morphology between pairs of spirals is a crucial discriminant between scenarios for Elias 2-27, and hence future studies must fit spiral arms individually. If Elias 2-27 continues to show symmetric tightly wound spiral arms in future observations, then we posit that it is the first observed example of a gravitationally unstable protostellar disc. △ Less

Submitted 21 May, 2018; originally announced May 2018.

Comments: 9 pages, 3 figures, accepted for publication in ApJL

Breakthrough revisited: investigating the requirements for growth of dust beyond the bouncing barrier

Authors: Richard A. Booth, Farzana Meru, Man Hoi Lee, Cathie J. Clarke

Abstract: For grain growth to proceed effectively and lead to planet formation a number of barriers to growth must be overcome. One such barrier, relevant for compact grains in the inner regions of the disc, is the `bouncing barrier' in which large grains ($\sim$ mm size) tend to bounce off each other rather than sticking. However, by maintaining a population of small grains it has been suggested that cm-si… ▽ More For grain growth to proceed effectively and lead to planet formation a number of barriers to growth must be overcome. One such barrier, relevant for compact grains in the inner regions of the disc, is the `bouncing barrier' in which large grains ($\sim$ mm size) tend to bounce off each other rather than sticking. However, by maintaining a population of small grains it has been suggested that cm-size particles may grow rapidly by sweeping up these small grains. We present the first numerically resolved investigation into the conditions under which grains may be lucky enough to grow beyond the bouncing barrier by a series of rare collisions leading to growth (so-called `breakthrough'). Our models support previous results, and show that in simple models breakthrough requires the mass ratio at which high velocity collisions transition to growth instead of causing fragmentation to be low, $φ\lesssim 50$. However, in models that take into account the dependence of the fragmentation threshold on mass-ratio, we find breakthrough occurs more readily, even if mass transfer is relatively inefficient. This suggests that bouncing may only slow down growth, rather than preventing growth beyond a threshold barrier. However, even when growth beyond the bouncing barrier is possible, radial drift will usually prevent growth to arbitrarily large sizes. △ Less

Submitted 8 December, 2017; originally announced December 2017.

Comments: 15 pages, 12 Figures. Accepted for publication in MNRAS

Towards a population synthesis model of self-gravitating disc fragmentation and tidal downsizing II: The effect of fragment-fragment interactions

Authors: D. H. Forgan, C. Hall, F. Meru, W. K. M. Rice

Abstract: It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the f… ▽ More It is likely that most protostellar systems undergo a brief phase where the protostellar disc is self-gravitating. If these discs are prone to fragmentation, then they are able to rapidly form objects that are initially of several Jupiter masses and larger. The fate of these disc fragments (and the fate of planetary bodies formed afterwards via core accretion) depends sensitively not only on the fragment's interaction with the disc, but with its neighbouring fragments. We return to and revise our population synthesis model of self-gravitating disc fragmentation and tidal downsizing. Amongst other improvements, the model now directly incorporates fragment-fragment interactions while the disc is still present. We find that fragment-fragment scattering dominates the orbital evolution, even when we enforce rapid migration and inefficient gap formation. Compared to our previous model, we see a small increase in the number of terrestrial-type objects being formed, although their survival under tidal evolution is at best unclear. We also see evidence for disrupted fragments with evolved grain populations - this is circumstantial evidence for the formation of planetesimal belts, a phenomenon not seen in runs where fragment-fragment interactions are ignored. In spite of intense dynamical evolution, our population is dominated by massive giant planets and brown dwarfs at large semimajor axis, which direct imaging surveys should, but only rarely, detect. Finally, disc fragmentation is shown to be an efficient manufacturer of free floating planetary mass objects, and the typical multiplicity of systems formed via gravitational instability will be low. △ Less

Submitted 7 November, 2017; v1 submitted 3 November, 2017; originally announced November 2017.

Comments: 14 pages, 17 figures, accepted for publication in MNRAS

Convergence of the critical cooling rate for protoplanetary disk fragmentation achieved; the key role of numerical dissipation of angular momentum

Authors: Hongping Deng, Lucio Mayer, Farzana Meru

Abstract: We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics(SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code (Hopkins 2015). Our aim is to understand the cause of the non-convergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementation… ▽ More We carry out simulations of gravitationally unstable disks using smoothed particle hydrodynamics(SPH) and the novel Lagrangian meshless finite mass (MFM) scheme in the GIZMO code (Hopkins 2015). Our aim is to understand the cause of the non-convergence of the cooling boundary for fragmentation reported in the literature. We run SPH simulations with two different artificial viscosity implementations, and compare them with MFM, which does not employ any artificial viscosity. With MFM we demonstrate convergence of the critical cooling time scale for fragmentation at β_{crit} =3.. Non-convergence persists in SPH codes, although it is significantly mitigated with schemes having reduced artificial viscosity such as inviscid SPH (ISPH) (Cullen & Dehnen 2010). We show how the non-convergence problem is caused by artificial fragmentation triggered by excessive dissipation of angular momentum in domains with large velocity derivatives. With increased resolution such domains become more prominent. Vorticity lags behind density due to numerical viscous dissipation in these regions, promoting collapse with longer cooling times. Such effect is shown to be dominant over the competing tendency of artificial viscosity to diminish with increasing resolution. When the initial conditions are first relaxed for several orbits, the flow is more regular, with lower shear and vorticity in non-axisymmetric regions, aiding convergence. Yet MFM is the only method that converges exactly. Our findings are of general interest as numerical dissipation via artificial viscosity or advection errors can also occur in grid-based codes. Indeed for the FARGO code values of β_{crit} significantly higher than our converged estimate have been reported in the literature. Finally, we discuss implications for giant planet formation via disk instability. △ Less

Submitted 9 November, 2017; v1 submitted 1 June, 2017; originally announced June 2017.

Comments: please check the ApJ version

Journal ref: Hongping Deng et al 2017 ApJ 847 43

On the origin of the spiral morphology in the Elias 2-27 circumstellar disc

Authors: Farzana Meru, Attila Juhasz, John D. Ilee, Cathie J. Clarke, Giovanni P. Rosotti, Richard A. Booth

Abstract: The young star Elias 2-27 has recently been observed to posses a massive circumstellar disc with two prominent large-scale spiral arms. In this Letter we perform three-dimensional Smoothed Particle Hydrodynamics simulations, radiative transfer modelling, synthetic ALMA imaging and an unsharped masking technique to explore three possibilities for the origin of the observed structures -- an undetect… ▽ More The young star Elias 2-27 has recently been observed to posses a massive circumstellar disc with two prominent large-scale spiral arms. In this Letter we perform three-dimensional Smoothed Particle Hydrodynamics simulations, radiative transfer modelling, synthetic ALMA imaging and an unsharped masking technique to explore three possibilities for the origin of the observed structures -- an undetected companion either internal or external to the spirals, and a self-gravitating disc. We find that a gravitationally unstable disc and a disc with an external companion can produce morphology that is consistent with the observations. In addition, for the latter, we find that the companion could be a relatively massive planetary mass companion (less than approximately 10 - 13 MJup) and located at large radial distances (between approximately 300 - 700 au). We therefore suggest that Elias 2-27 may be one of the first detections of a disc undergoing gravitational instabilities, or a disc that has recently undergone fragmentation to produce a massive companion. △ Less

Submitted 15 March, 2017; originally announced March 2017.

Comments: 6 pages, 5 figures, accepted for publication in ApJ Letters

Grand challenges in protoplanetary disc modelling

Authors: Thomas J. Haworth, John D. Ilee, Duncan H. Forgan, Stefano Facchini, Daniel J. Price, Community authors, Dominika M. Boneberg, Richard A. Booth, Cathie J. Clarke, Jean-François Gonzalez, Mark A. Hutchison, Inga Kamp, Guillaume Laibe, Wladimir Lyra, Farzana Meru, Subhanjoy Mohanty, Olja Panić, Ken Rice, Takeru Suzuki, Catherine Walsh, Peter Woitke

Abstract: 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 t… ▽ More 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. △ Less

Submitted 20 September, 2016; v1 submitted 3 August, 2016; originally announced August 2016.

Comments: Accepted for publication in PASA. 25 pages

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

Authors: Stefan Kraus, John D. Monnier, Michael J. Ireland, Gaspard Duchene, Catherine Espaillat, Sebastian Hoenig, Attila Juhasz, Chris Mordasini, Johan Olofsson, Claudia Paladini, Keivan Stassun, Neal Turner, Gautam Vasisht, Tim J. Harries, Matthew R. Bate, Jean-Francois Gonzalez, Alexis Matter, Zhaohuan Zhu, Olja Panic, Zsolt Regaly, Alessandro Morbidelli, Farzana Meru, Sebastian Wolf, John Ilee, Jean-Philippe Berger , et al. (53 additional authors not shown)

Abstract: 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 prot… ▽ More 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. △ Less

Submitted 16 August, 2016; v1 submitted 1 August, 2016; originally announced August 2016.

Comments: 13 pages, 5 figures, Proceedings of SPIE 2016

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

Authors: Farzana Meru

Abstract: 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 an… ▽ More 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. △ Less

Submitted 16 September, 2015; v1 submitted 11 September, 2015; originally announced September 2015.

Comments: 11 pages, 11 figures, accepted for publication by MNRAS. Movies of Simulation 1 available at www.ast.cam.ac.uk/~fmeru/Movies/mass_movement_sigma.mov and www.ast.cam.ac.uk/~fmeru/Movies/mass_movement_vR.mov

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

Authors: Matej Malik, Farzana Meru, Lucio Mayer, Michael R. Meyer

Abstract: 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 ha… ▽ More 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. △ Less

Submitted 23 February, 2015; originally announced February 2015.

Comments: Accepted by ApJ, 22 pages, 13 figures, 6 tables

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

Authors: Farzana Meru, Sascha P. Quanz, Maddalena Reggiani, Clement Baruteau, Jaime E. Pineda

Abstract: 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 simu… ▽ More 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. △ Less

Submitted 19 November, 2014; originally announced November 2014.

Comments: Submitted to ApJ. Comments very welcome

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

Authors: Maddalena Reggiani, Sascha P. Quanz, Michael R. Meyer, Laurent Pueyo, Olivier Absil, Adam Amara, Guillem Anglada, Henning Avenhaus, Julien H. Girard, Carlos Carrasco Gonzalez, Graham James, Dimitri Mawet, Farzana Meru, Julien Milli, Mayra Osorio, Schuyler Wolff, Jose-Maria Torrelles

Abstract: 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… ▽ More 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. △ Less

Submitted 4 August, 2014; originally announced August 2014.

Comments: Accepted to ApJL, see also Biller et al. 2014

Resolved images of the protoplanetary disk around HD 100546 with ALMA

Authors: Jaime E. Pineda, Sascha P. Quanz, Farzana Meru, Gijs D. Mulders, Michael R. Meyer, Olja Panić, Henning Avenhaus

Abstract: 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 obs… ▽ More 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. △ Less

Submitted 22 May, 2014; originally announced May 2014.

Comments: 6 pages, 4 figures. Accepted in ApJL

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

Authors: Antonio Garufi, Sascha P Quanz, Henning Avenhaus, Esther Buenzli, Carsten Dominik, Farzana Meru, Michael R Meyer, Paola Pinilla, Hans Martin Schmid, Sebastian Wolf

Abstract: 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… ▽ More 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. △ Less

Submitted 17 November, 2013; originally announced November 2013.

Comments: 11 pages, 6 figures; accepted for publication in A&A

Giant Planet Formation, Evolution, and Internal Structure

Authors: Ravit Helled, Peter Bodenheimer, Morris Podolak, Aaron Boley, Farzana Meru, Sergei Nayakshin, Jonathan J. Fortney, Lucio Mayer, Yann Alibert, Alan P. Boss

Abstract: The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ide… ▽ More The large number of detected giant exoplanets offers the opportunity to improve our understanding of the formation mechanism, evolution, and interior structure of gas giant planets. The two main models for giant planet formation are core accretion and disk instability. There are substantial differences between these formation models, including formation timescale, favorable formation location, ideal disk properties for planetary formation, early evolution, planetary composition, etc. First, we summarize the two models including their substantial differences, advantages, and disadvantages, and suggest how theoretical models should be connected to available (and future) data. We next summarize current knowledge of the internal structures of solar- and extrasolar- giant planets. Finally, we suggest the next steps to be taken in giant planet exploration. △ Less

Submitted 5 November, 2013; originally announced November 2013.

Comments: Accepted for publication as a chapter in Protostars and Planets VI, to be published in 2014 by University of Arizona Press

Large grains can grow in circumstellar discs

Authors: Farzana Meru, Marina Galvagni, Christoph Olczak, Pascale Garaud

Abstract: 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… ▽ More 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. △ Less

Submitted 23 August, 2013; originally announced August 2013.

Comments: Submitted to the conference proceedings of the IAU Symposium 299 - Exploring the formation and evolution of planetary systems. 2 pages; 2 figures

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

Authors: Farzana Meru, Ralf J. Geretshauser, Christoph Schaefer, Roland Speith, Wilhelm Kley

Abstract: 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 tha… ▽ More 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. △ Less

Submitted 4 August, 2013; originally announced August 2013.

Comments: 21 pages, 18 figures, accepted by MNRAS

Growth of grains in Brown Dwarf disks

Authors: Farzana Meru, Marina Galvagni, Christoph Olczak

Abstract: 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 ra… ▽ More 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. △ Less

Submitted 14 July, 2013; originally announced July 2013.

Comments: 6 pages, 4 figures, accepted by ApJL

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

Authors: Farzana Meru

Abstract: 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… ▽ More 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. △ Less

Submitted 30 May, 2013; originally announced May 2013.

Comments: Submitted to the conference proceedings of: Instabilities and Structures in Proto-Planetary Disks. 5 pages; 4 figures

Journal ref: EPJ Web of Conferences; Volume 46, 2013; Instabilities and Structures in Proto-Planetary Disks

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

Authors: Farzana Meru, Matthew R. Bate

Abstract: 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 larg… ▽ More 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. △ Less

Submitted 5 September, 2012; originally announced September 2012.

Comments: Accepted for publication by MNRAS. 26 pages, 17 figures

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

Authors: Pascale Garaud, Farzana Meru, Marina Galvagni, Christoph Olczak

Abstract: 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 fa… ▽ More 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. △ Less

Submitted 15 December, 2012; v1 submitted 31 August, 2012; originally announced September 2012.

Comments: Accepted for publication in ApJ. Additional appendix included. Minor changes from previous versions. 46 pages, 10 figures

Treating gravity in thin disk simulations

Authors: T. W. A. Müller, W. Kley, F. Meru

Abstract: In 2D simulations of thin gaseous disks with embedded planets or self-gravity the gravitational potential needs to be smoothed to avoid singularities in the numerical evaluation of the gravitational potential or force. In order to correctly resemble the realistic case of vertically extended 3D disks the softening prescription used in 2D needs to be adjusted properly. In this paper we analyze the e… ▽ More In 2D simulations of thin gaseous disks with embedded planets or self-gravity the gravitational potential needs to be smoothed to avoid singularities in the numerical evaluation of the gravitational potential or force. In order to correctly resemble the realistic case of vertically extended 3D disks the softening prescription used in 2D needs to be adjusted properly. In this paper we analyze the embedded planet and the self-gravity case and provide a method to evaluate the required smoothing in 2D simulations of thin disks. Starting from the averaged hydrodynamic equations and using a vertically isothermal disk model, we calculate the force. We compare our results to the often used Plummer form of the potential which runs as \propto 1/(r^2+ε^2)^{1/2}. For that purpose we compute the required smoothing length εas a function of distance r to the planet or to a disk element within a self-gravitating disk. We find that for larger distances εis determined solely by the vertical disk thickness H. For the planet case we find that outside r\simH a value of ε= 0.7H describes the averaged force very well, while in the self-gravitating disk the value needs to be larger, ε= 1.2H. For smaller distances the smoothing needs to be reduced significantly. Comparing torque densities of 3D and 2D simulations we show that the modification to the vertical density stratification as induced by an embedded planet needs to be taken into account to obtain agreeing results. In disk fragmentation simulations the choice of εcan determine whether a disk will fragment or not. Because a wrong smoothing length can change even the direction of migration, it is very important to include the effect of the planet on the local scale height in 2D planet-disk simulations. We provide an approximate and fast method for this purpose which gives very good agreement to full 3D simulations. △ Less

Submitted 4 April, 2012; v1 submitted 7 March, 2012; originally announced March 2012.

Comments: 13 pages, 19 figures, accepted for publication in A&A

Journal ref: A&A 541, A123 (2012)

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

Authors: Sijme-Jan Paardekooper, Clement Baruteau, Farzana Meru

Abstract: 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 coo… ▽ More 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. △ Less

Submitted 22 June, 2011; originally announced June 2011.

Comments: 5 pages, 5 figures, accepted for publication in MNRAS

Rapid inward migration of planets formed by gravitational instability

Authors: Clément Baruteau, Farzana Meru, Sijme-Jan Paardekooper

Abstract: 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, i… ▽ More 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. △ Less

Submitted 2 June, 2011; originally announced June 2011.

Comments: 13 pages, 10 figures, accepted for publication in MNRAS

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

Authors: Ralf J. Geretshauser, Farzana Meru, Roland Speith, Wilhelm Kley

Abstract: 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… ▽ More 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. △ Less

Submitted 19 May, 2011; v1 submitted 4 May, 2011; originally announced May 2011.

Comments: 14 pages, 11 figures, 5 tables, to be published in Astronomy and Astrophysics

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

Authors: Farzana Meru, Matthew R. Bate

Abstract: 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 re… ▽ More 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. △ Less

Submitted 25 November, 2010; originally announced November 2010.

Comments: Submitted to the conference proceedings of the IAU Symposium 276 - The Astrophysics of Planetary Systems: Formation, Structure, and Dynamical Evolution. 3 pages; 2 figures

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

Authors: Farzana Meru, Matthew R. Bate

Abstract: 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 resoluti… ▽ More 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. △ Less

Submitted 3 November, 2010; originally announced November 2010.

Comments: Accepted for publication by MNRAS Letters. 6 pages, 3 figures

On the fragmentation criteria of self-gravitating protoplanetary discs

Authors: Farzana Meru, Matthew R. Bate

Abstract: 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)… ▽ More 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. △ Less

Submitted 3 August, 2010; originally announced August 2010.

Comments: Accepted for publication by MNRAS. 15 pages, 18 figures

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

Authors: Farzana Meru, Matthew R. Bate

Abstract: 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 wheth… ▽ More 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. △ Less

Submitted 21 April, 2010; originally announced April 2010.

Comments: Accepted for publication by MNRAS. 11 pages, 12 figures.