-
Gravitational instability in a planet-forming disk
Authors:
Jessica Speedie,
Ruobing Dong,
Cassandra Hall,
Cristiano Longarini,
Benedetta Veronesi,
Teresa Paneque-Carreño,
Giuseppe Lodato,
Ya-Wen Tang,
Richard Teague,
Jun Hashimoto
Abstract:
The canonical theory for planet formation in circumstellar disks proposes that planets are grown from initially much smaller seeds. The long-considered alternative theory proposes that giant protoplanets can be formed directly from collapsing fragments of vast spiral arms induced by gravitational instability -- if the disk is gravitationally unstable. For this to be possible, the disk must be mass…
▽ More
The canonical theory for planet formation in circumstellar disks proposes that planets are grown from initially much smaller seeds. The long-considered alternative theory proposes that giant protoplanets can be formed directly from collapsing fragments of vast spiral arms induced by gravitational instability -- if the disk is gravitationally unstable. For this to be possible, the disk must be massive compared to the central star: a disk-to-star mass ratio of 1/10 is widely held as the rough threshold for triggering gravitational instability, inciting significant non-Keplerian dynamics and generating prominent spiral arms. While estimating disk masses has historically been challenging, the motion of the gas can reveal the presence of gravitational instability through its effect on the disk velocity structure. Here we present kinematic evidence of gravitational instability in the disk around AB Aurigae, using deep observations of 13CO and C18O line emission with the Atacama Large Millimeter/submillimeter Array (ALMA). The observed kinematic signals strongly resemble predictions from simulations and analytic modelling. From quantitative comparisons, we infer a disk mass of up to 1/3 the stellar mass enclosed within 1" to 5" on the sky.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
X-ray and optical polarization aligned with the radio jet ejecta in GX 339-4
Authors:
G. Mastroserio,
B. De Marco,
M. C. Baglio,
F. Carotenuto,
S. Fabiani,
T. D. Russell,
F. Capitanio,
Y. Cavecchi,
S. Motta,
D. M. Russell,
M. Dovciak,
M. Del Santo,
K. Alabarta,
A. Ambrifi,
S. Campana,
P. Casella,
S. Covino,
G. Illiano,
E. Kara,
E. V. Lai,
G. Lodato,
A. Manca,
I. Mariani,
A. Marino,
C. Miceli
, et al. (5 additional authors not shown)
Abstract:
We present the first X-ray polarization measurements of GX 339-4. IXPE observed this source twice during its 2023-2024 outburst, once in the soft-intermediate state and again during a soft state. The observation taken during the intermediate state shows significant ($4σ$) polarization degree P = $1.3\% \pm 0.3\%$ and polarization angle $θ$ = -74\degree $\pm$ 7\degree only in the 3 - 8 keV band. FO…
▽ More
We present the first X-ray polarization measurements of GX 339-4. IXPE observed this source twice during its 2023-2024 outburst, once in the soft-intermediate state and again during a soft state. The observation taken during the intermediate state shows significant ($4σ$) polarization degree P = $1.3\% \pm 0.3\%$ and polarization angle $θ$ = -74\degree $\pm$ 7\degree only in the 3 - 8 keV band. FORS2 at VLT observed the source simultaneously detecting optical polarization in the B, V, R, I bands (between $0.1%$ and $0.7\%$), all roughly aligned with the X-ray polarization. We also detect a discrete jet knot from radio observations taken later in time; this knot would have been ejected from the system around the same time as the hard-to-soft X-ray state transition and a bright radio flare occurred $\sim$3 months earlier. The proper motion of the jet knot provides a direct measurement of the jet orientation angle on the plane of the sky at the time of the ejection. We find that both the X-ray and optical polarization angles are aligned with the direction of the ballistic jet.
△ Less
Submitted 9 August, 2024;
originally announced August 2024.
-
The evolution of the $M_{\mathrm{d}}-M_{\star}$ and $\dot M-M_{\star}$ correlations traces protoplanetary disc dispersal
Authors:
Alice Somigliana,
Leonardo Testi,
Giovanni Rosotti,
Claudia Toci,
Giuseppe Lodato,
Rossella Anania,
Benoît Tabone,
Marco Tazzari,
Ralf Klessen,
Ugo Lebreuilly,
Patrick Hennebelle,
Sergio Molinari
Abstract:
(Abridged) Observational surveys of entire star-forming regions have provided evidence of power-law correlations between the disc properties and the stellar mass, especially the disc mass (${M_d \propto M_*}^{λ_m}$) and the accretion rate ($\dot M \propto {M_*}^{λ_{acc}}$). Whether the secular disc evolution affects said correlations is still debated: while the purely viscous scenario has been pro…
▽ More
(Abridged) Observational surveys of entire star-forming regions have provided evidence of power-law correlations between the disc properties and the stellar mass, especially the disc mass (${M_d \propto M_*}^{λ_m}$) and the accretion rate ($\dot M \propto {M_*}^{λ_{acc}}$). Whether the secular disc evolution affects said correlations is still debated: while the purely viscous scenario has been probed, other mechanisms could impact differently. We study the evolution of the slopes $λ_m$ and $λ_{acc}$ in the wind-driven and hybrid case and compare it to the viscous prediction, using a combination of analytical calculations and numerical simulations (performed with the 1D population synthesis code Diskpop, that we also present and release). Assuming $M_d(0) \propto {M_*}^{λ_{m, 0}}$ and $\dot M(0) \propto {M_*}^{λ_{acc, 0}}$ as initial conditions, we find that viscous and hybrid accretion preserve the shape of the correlations and evolve their slope; on the other hand, MHD winds change the shape of the correlations, bending them according to the scaling of the accretion timescale with the stellar mass. We also show how a spread in the initial conditions conceals this behaviour. We then analyse the impact of disc dispersal, and find that the currently available sample sizes ($\sim 30$ discs at 5 Myr) introduce stochastic oscillations in the slopes evolution, which dominate over the physical signatures. Increasing the sample size could mitigate this issue: $\sim 140$ discs at 5 Myr, corresponding to the complete Upper Sco sample, would give small enough error bars to use the evolution of the slopes as a proxy for the driving mechanism of disc evolution. Finally, we discuss how the observational claim of steepening slopes necessarily leads to an initially steeper $M_d - M_*$ correlation with respect to $\dot M - M_*$.
△ Less
Submitted 30 July, 2024;
originally announced July 2024.
-
Angular momentum transport via gravitational instability in the Elias 2-27 disc
Authors:
Cristiano Longarini,
Giuseppe Lodato,
Cathie J. Clarke,
Jessica Speedie,
Teresa Paneque-Carreno,
Edoardo Arrigoni,
Pietro Curone,
Claudia Toci,
Cassandra Hall
Abstract:
Gravitational instability is thought to be one of the main drivers of angular momentum transport in young protoplanetary discs. The disc around Elias 2-27 offers a unique example of gravitational instability at work. It is young and massive, displaying two prominent spiral arms in dust continuum emission and global non-axisymmetric kinematic signatures in molecular line data. In this work, we used…
▽ More
Gravitational instability is thought to be one of the main drivers of angular momentum transport in young protoplanetary discs. The disc around Elias 2-27 offers a unique example of gravitational instability at work. It is young and massive, displaying two prominent spiral arms in dust continuum emission and global non-axisymmetric kinematic signatures in molecular line data. In this work, we used archival ALMA observations of $^{13}$CO line emission to measure the efficiency of angular momentum transport in the Elias 2-27 system through the kinematic signatures generated by gravitational instability, known as 'GI wiggles'. Assuming the angular momentum is transported by the observed spiral structure and leveraging previously-derived dynamical disc mass measurements, the amount of angular momentum transport we found corresponds to an $α-$viscosity of $α=0.038\pm0.018$. This value implies an accretion rate onto the central star of $\log_{10}\dot{M}_\star=-6.99\pm0.17\text{M}_\odot/\text{yr, which}$ reproduces the one observed value of $\log_{10}\dot{M}_{\star,\text{obs}}=-7.2\pm0.5\text{M}_\odot/\text{yr }$ very well. The excellent agreement we have found serves as further proof that gravitational instability is the main driver of angular momentum transport acting in this system.
△ Less
Submitted 11 June, 2024; v1 submitted 9 June, 2024;
originally announced June 2024.
-
Weighing protoplanetary discs with kinematics: physical model, method and benchmark
Authors:
Benedetta Veronesi,
Cristiano Longarini,
Giuseppe Lodato,
Guillaume Laibe,
Cassandra Hall,
Stefano Facchini,
Leonardo Testi
Abstract:
The mass of protoplanetary discs sets the amount of material available for planet formation, determines the level of coupling between gas and dust, and possibly sets gravitational instabilities. Measuring mass of discs is challenging, since it is not possible to directly detect H$_2$, and CO-based estimates remain poorly constrained. An alternative method that does not rely on tracers-to-H$_2$ rat…
▽ More
The mass of protoplanetary discs sets the amount of material available for planet formation, determines the level of coupling between gas and dust, and possibly sets gravitational instabilities. Measuring mass of discs is challenging, since it is not possible to directly detect H$_2$, and CO-based estimates remain poorly constrained. An alternative method that does not rely on tracers-to-H$_2$ ratios has recently been proposed to dynamically measure the disc mass altogether with the star mass and the disc critical radius by looking at deviations from Keplerian rotation induced by the self-gravity of the disc. So far, this method has been applied to weigh three protoplanetary discs: Elias 2-27, IM Lup and GM Aurigae. We provide here a numerical benchmark of the method by simulating isothermal self-gravitating discs with a range of masses from 0.01 to $0.2 \,M_{\odot}$ with the phantom code and post-process them with radiative transfer (mcfost) to obtain synthetic observations. We find that dynamical weighing allows to retrieve the expected value of disc masses as long as the disc-to-star mass ratio is larger than $M_d/M_\star=0.05$. The estimated uncertainty for the disc mass measurement is $\sim 25\%$.
△ Less
Submitted 24 May, 2024;
originally announced May 2024.
-
Planet-driven spirals in protoplanetary discs: limitations of the semi-analytical theory for observations
Authors:
D. Fasano,
A. J. Winter,
M. Benisty,
G. Rosotti,
A. Ruzza,
G. Lodato,
C. Toci,
T. Hilder,
A. Izquierdo,
D. Price
Abstract:
Detecting protoplanets during their formation stage is an important but elusive goal of modern astronomy. Kinematic detections via the spiral wakes in the gaseous disc are a promising avenue to achieve this goal. We aim to test the applicability to observations in the low and intermediate planet mass regimes of a commonly used semi-analytical model for planet induced spiral waves. In contrast with…
▽ More
Detecting protoplanets during their formation stage is an important but elusive goal of modern astronomy. Kinematic detections via the spiral wakes in the gaseous disc are a promising avenue to achieve this goal. We aim to test the applicability to observations in the low and intermediate planet mass regimes of a commonly used semi-analytical model for planet induced spiral waves. In contrast with previous works which proposed to use the semi-analytical model to interpret observations, in this study we analyse for the first time both the structure of the velocity and density perturbations. We run a set of FARGO3D hydrodynamic simulations and compare them with the output of the semi-analytic model in the code wakeflow, which is obtained by solving Burgers' equation using the simulations as an initial condition. We find that the velocity field derived from the analytic theory is discontinuous at the interface between the linear and nonlinear regions. After 0.2 r$_p$ from the planet, the behaviour of the velocity field closely follows that of the density perturbations. In the low mass limit, the analytical model is in qualitative agreement with the simulations, although it underestimates the azimuthal width and the amplitude of the perturbations, predicting a stronger decay but a slower azimuthal advance of the shock fronts. In the intermediate regime, the discrepancy increases, resulting in a different pitch angle between the spirals of the simulations and the analytic model. The implementation of a fitting procedure based on the minimisation of intensity residuals is bound to fail due to the deviation in pitch angle between the analytic model and the simulations. In order to apply this model to observations, it needs to be revisited accounting also for higher planet masses.
△ Less
Submitted 19 June, 2024; v1 submitted 24 May, 2024;
originally announced May 2024.
-
Eddington envelopes: The fate of stars on parabolic orbits tidally disrupted by supermassive black holes
Authors:
Daniel J. Price,
David Liptai,
Ilya Mandel,
Joanna Shepherd,
Giuseppe Lodato,
Yuri Levin
Abstract:
Stars falling too close to massive black holes in the centres of galaxies can be torn apart by the strong tidal forces. Simulating the subsequent feeding of the black hole with disrupted material has proved challenging because of the range of timescales involved. Here we report a set of simulations that capture the relativistic disruption of the star, followed by one year of evolution of the retur…
▽ More
Stars falling too close to massive black holes in the centres of galaxies can be torn apart by the strong tidal forces. Simulating the subsequent feeding of the black hole with disrupted material has proved challenging because of the range of timescales involved. Here we report a set of simulations that capture the relativistic disruption of the star, followed by one year of evolution of the returning debris stream. These reveal the formation of an expanding asymmetric bubble of material extending to hundreds of astronomical units -- an outflowing Eddington envelope with an optically thick inner region. Such envelopes have been hypothesised as the reprocessing layer needed to explain optical/UV emission in tidal disruption events, but never produced self-consistently in a simulation. Our model broadly matches the observed light curves with low temperatures, faint luminosities, and line widths of 10,000--20,000 km/s.
△ Less
Submitted 9 July, 2024; v1 submitted 14 April, 2024;
originally announced April 2024.
-
Orbital dynamics in the GG Tau A system: investigating its enigmatic disc
Authors:
Claudia Toci,
Simone Ceppi,
Nicolás Cuello,
Gaspard Duchêne,
Enrico Ragusa,
Giuseppe Lodato,
Francesca Farina,
François Ménard,
Hossam Aly
Abstract:
GG Tau is one of the most studied multiple young stellar systems: GG Tau A is a hierarchical triple surrounded by a massive disc and its companion, GG Tau B, is also a binary. Despite numerous observational attempts, an understanding of the geometry of the GG Tau A system is still elusive. We provide new astrometric measures of the system and we run a set of hydrodynamical simulations with two rep…
▽ More
GG Tau is one of the most studied multiple young stellar systems: GG Tau A is a hierarchical triple surrounded by a massive disc and its companion, GG Tau B, is also a binary. Despite numerous observational attempts, an understanding of the geometry of the GG Tau A system is still elusive. We provide new astrometric measures of the system and we run a set of hydrodynamical simulations with two representative orbits to test how they impact a disc composed of dust and gas. We test the dynamical evolution of the two scenarios on short and long timescales. We obtain synthetic flux emission from our simulations and we compare them with 1300 $μ$m ALMA dust continuum emission and 1.67 $μ$m SPHERE dust scattering images to infer the most likely orbital arrangement. We extend the analysis of the binary orbital parameters using six new epochs from archival data, showing that the current measurements alone are not capable of breaking the degeneracy between families of coplanar and misaligned orbits. We found that the time-scale for the onset of the disc eccentricity growth, $τ_{ecc}$, is a fundamental time-scale for the morphology of the system. Results from numerical simulations show that the best match between is obtained with the misaligned configuration ($Δθ= 30^\circ$) on timescales shorter than $τ_{ecc}$. The results exhibit an almost circular cavity and dust ring. However, for both scenarios, the cavity size and its eccentricity quickly grow for timescales longer than $τ_{ecc}$ and the models do not reproduce the observed morphology anymore. This implies that either the age of the system is shorter than $τ_{ecc}$ or that the disc eccentricity growth is not triggered or dissipated. This finding raises questions on the future evolution of the GG Tau A system and, more in general, on the time evolution of eccentric binaries and their circumbinary discs.
△ Less
Submitted 11 April, 2024;
originally announced April 2024.
-
Generalisation of the Spectral Difference scheme for the diffused-interface five equation model
Authors:
Niccolò Tonicello,
Guido Lodato,
Matthias Ihme
Abstract:
The present work focuses on the generalisation of the Spectral Difference (SD) scheme to the reduced Baer-Nunziato system known as five-equation model for the simulation of two immiscible compressible fluids. This five equation model is considered with the additional Allen-Cahn regularisation to avoid both over-diffusion and over-thinning of the phase field representing the interface. Finally, in…
▽ More
The present work focuses on the generalisation of the Spectral Difference (SD) scheme to the reduced Baer-Nunziato system known as five-equation model for the simulation of two immiscible compressible fluids. This five equation model is considered with the additional Allen-Cahn regularisation to avoid both over-diffusion and over-thinning of the phase field representing the interface. Finally, in order to preserve contact discontinuities, in the reconstruction step of the spectral difference scheme, a change of variables from conservative to primitive is used. This approach is shown to be beneficial in avoiding pressure oscillations at material interfaces. An extensive series of numerical tests are proposed to assess accuracy and robustness of the present method. Both kinematic (Rider-Kothe vortex) and two-phase flow problems (Rayleigh-Taylor instability, shock-droplet interaction, Taylor-Green vortex) are considered.
△ Less
Submitted 29 March, 2024; v1 submitted 28 March, 2024;
originally announced March 2024.
-
DBNets: A publicly available deep learning tool to measure the masses of young planets in dusty protoplanetary discs
Authors:
Alessandro Ruzza,
Giuseppe Lodato,
Giovanni Pietro Rosotti
Abstract:
Current methods to characterize embedded planets in protoplanetary disc observations are severely limited either in their ability to fully account for the observed complex physics or in their computational and time costs. To address this shortcoming, we developed DBNets: a deep learning tool, based on convolutional neural networks, that analyses substructures observed in the dust continuum emissio…
▽ More
Current methods to characterize embedded planets in protoplanetary disc observations are severely limited either in their ability to fully account for the observed complex physics or in their computational and time costs. To address this shortcoming, we developed DBNets: a deep learning tool, based on convolutional neural networks, that analyses substructures observed in the dust continuum emission of protoplanetary discs to quickly infer the mass of allegedly embedded planets. We focussed on developing a method to reliably quantify not only the planet mass, but also the associated uncertainty introduced by our modelling and adopted techniques. Our tests gave promising results achieving an 87% reduction of the log Mp mean squared error with respect to an analytical formula fitted on the same data (DBNets metrics: lmse 0.016, r2-score 97%). With the goal of providing the final user of DBNets with all the tools needed to interpret their measurements and decide on their significance, we extensively tested our tool on out-of-distribution data. We found that DBNets can identify inputs strongly outside its training scope returning an uncertainty above a specific threshold and we thus provided a rejection criterion that helps determine the significance of the results obtained. Additionally, we outlined some limitations of our tool: it can be reliably applied only on discs observed with inclinations below approximately 60°, in the optically thin regime, with a resolution 8 times better than the gap radial location and with a signal-to-noise ratio higher than approximately ten. Finally, we applied DBNets to 33 actual observations of protoplanetary discs measuring the mass of 48 proposed planets and comparing our results with the available literature. We confirmed that most of the observed gaps imply planets in the sub-Jupiter regime. DBNets is publicly available at dbnets.fisica.unimi.it.
△ Less
Submitted 19 February, 2024;
originally announced February 2024.
-
Rotation curves in protoplanetary disks with thermal stratification
Authors:
Paola Martire,
Cristiano Longarini,
Giuseppe Lodato,
Giovanni P. Rosotti,
Andrew Winter,
Stefano Facchini,
Caitlyn Hardiman,
Myriam Benisty,
Jochen Stadler,
Andrés F. Izquierdo,
Leonardo Testi
Abstract:
In recent years the gas kinematics probed by molecular lines detected with ALMA has opened a new window to study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure and correctly interpreting these data allow us to gain a better comprehension of the planet formation process. We investigate the impact of thermal strat…
▽ More
In recent years the gas kinematics probed by molecular lines detected with ALMA has opened a new window to study protoplanetary disks. High spatial and spectral resolution observations have revealed the complexity of protoplanetary disk structure and correctly interpreting these data allow us to gain a better comprehension of the planet formation process. We investigate the impact of thermal stratification on the azimuthal velocity of protoplanetary disks. High resolution gas observations are showing velocity differences between CO isotopologues, which cannot be adequately explained with vertically isothermal models. The aim of this work is to determine whether a stratified model can explain this discrepancy. We analytically solve the hydrostatic equilibrium for a stratified disk and we derive the azimuthal velocity. We test the model with SPH numerical simulations and then we use it to fit for star mass, disk mass and scale radius of the sources in the MAPS sample. In particular, we use 12CO and 13CO datacubes.
△ Less
Submitted 11 March, 2024; v1 submitted 19 February, 2024;
originally announced February 2024.
-
Probing initial distributions of orbital eccentricity and disc misalignment via polar discs
Authors:
Simone Ceppi,
Nicolás Cuello,
Giuseppe Lodato,
Cristiano Longarini,
Daniel J. Price,
Daniel Elsender,
Matthew R. Bate
Abstract:
In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the di…
▽ More
In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the distributions of initial eccentricity and mutual inclination in multiple stellar systems at birth.
△ Less
Submitted 8 January, 2024;
originally announced January 2024.
-
Dual-Band Observations of the Asymmetric Ring around CIDA 9A: Dead or Alive?
Authors:
Daniel Harsono,
Feng Long,
Paola Pinilla,
Alessia A. Rota,
Carlo F. Manara,
Gregory J. Herczeg,
Doug Johnstone,
Giovanni Rosotti,
Giuseppe Lodato,
Francois Menard,
Marco Tazzari,
Yangfan Shi
Abstract:
While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data.…
▽ More
While the most exciting explanation of the observed dust asymmetries in protoplanetary disks is the presence of protoplanets, other mechanisms can also form the dust features. This paper presents dual-wavelength Atacama Large Millimeter/submillimeter Array (ALMA) observations of a large asymmetric dusty ring around the M-type star CIDA 9A. We detect a dust asymmetry in both 1.3 mm and 3.1 mm data. To characterize the asymmetric structure, a parametric model is used to fit the observed visibilities. We report a tentative azimuthal shift of the dust emission peaks between the observations at the two wavelengths. This shift is consistent with a dust trap caused by a vortex, which may be formed by an embedded protoplanet or other hydrodynamical instabilities, such as a dead zone. Deep high-spatial observations of dust and molecular gas are needed to constrain the mechanisms that formed the observed millimeter cavity and dust asymmetry in the protoplanetary disk around CIDA 9A.
△ Less
Submitted 17 October, 2023;
originally announced October 2023.
-
Testing external photoevaporation in the $σ$-Orionis cluster with spectroscopy and disk mass measurements
Authors:
K. Maucó,
C. F. Manara,
M. Ansdell,
G. Bettoni,
R. Claes,
J. Alcala,
A. Miotello,
S. Facchini,
T. J. Haworth,
G. Lodato,
J. P. Williams
Abstract:
The evolution of protoplanetary disks is regulated by an interplay of several processes, either internal to the system or related to the environment. As most of the stars and planets have formed in massive stellar clusters, studying the effects of UV radiation on disk evolution is of paramount importance. Here we test the impact of external photoevaporation on the evolution of disks in the $σ$ Ori…
▽ More
The evolution of protoplanetary disks is regulated by an interplay of several processes, either internal to the system or related to the environment. As most of the stars and planets have formed in massive stellar clusters, studying the effects of UV radiation on disk evolution is of paramount importance. Here we test the impact of external photoevaporation on the evolution of disks in the $σ$ Orionis cluster by conducting the first combined large-scale UV to IR spectroscopic and mm-continuum survey of this region. We study a sample of 50 targets located at increasing distances from the central, OB system $σ$ Ori. We combine new VLT/X-Shooter spectra with new and previously published ALMA measurements of disk dust and gas fluxes and masses. We confirm the previously found decrease of $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. This is particularly evident when considering the disks around the more massive stars ($\ge$ 0.4 $M_{\odot}$), where those located in the inner part ($<$ 0.5 pc) have $M_{\rm dust}$ about an order of magnitude lower than the more distant ones. About half of the sample is located in the region of the $\dot{M}_{\rm acc}$ vs $M_{\rm disk}$ expected by models of external photoevaporation, namely showing shorter disk lifetimes. These are observed for all targets with projected separation from $σ$ Ori $<$ 0.5 pc, proving that the presence of a massive stellar system affects disk evolution. External photoevaporation is a viable mechanism to explain the observed shorter disk lifetimes and lower $M_{\rm dust}$ in the inner $\sim$0.5 pc of the cluster. Follow-up observations of the low stellar mass targets are crucial to confirm the dependence of the external photoevaporation process with stellar host mass. This work confirms that the effects of external photoevaporation are significant down to impinging radiation as low as $\sim 10^{4}$ G$_0$.
△ Less
Submitted 11 September, 2023;
originally announced September 2023.
-
The time evolution of $M_{\mathrm{d}}/\dot M$ in protoplanetary discs as a way to disentangle between viscosity and MHD winds
Authors:
Alice Somigliana,
Leonardo Testi,
Giovanni Rosotti,
Claudia Toci,
Giuseppe Lodato,
Benoît Tabone,
Carlo Manara,
Marco Tazzari
Abstract:
As the classic viscous paradigm for protoplanetary disk accretion is challenged by the observational evidence of low turbulence, the alternative scenario of MHD disk winds is being explored as potentially able to reproduce the same observed features traditionally explained with viscosity. Although the two models lead to different disk properties, none of them has been ruled out by observations - m…
▽ More
As the classic viscous paradigm for protoplanetary disk accretion is challenged by the observational evidence of low turbulence, the alternative scenario of MHD disk winds is being explored as potentially able to reproduce the same observed features traditionally explained with viscosity. Although the two models lead to different disk properties, none of them has been ruled out by observations - mainly due to instrumental limitations. In this work, we present a viable method to distinguish between the viscous and MHD framework based on the different evolution of the distribution in the disk mass ($M_{\mathrm{d}}$) - accretion rate ($\dot M$) plane of a disk population. With a synergy of analytical calculations and 1D numerical simulations, performed with the population synthesis code \texttt{Diskpop}, we find that both mechanisms predict the spread of the observed ratio $M_{\mathrm{d}}/\dot M$ in a disk population to decrease over time; however, this effect is much less pronounced in MHD-dominated populations as compared to purely viscous populations. Furthermore, we demonstrate that this difference is detectable with the current observational facilities: we show that convolving the intrinsic spread with the observational uncertainties does not affect our result, as the observed spread in the MHD case remains significantly larger than in the viscous scenario. While the most recent data available show a better agreement with the wind model, ongoing and future efforts to obtain direct gas mass measurements with ALMA and ngVLA will cause a reassessment of this comparison in the near future.
△ Less
Submitted 7 September, 2023;
originally announced September 2023.
-
Radio multiwavelength analysis of the compact disk CX Tau: Presence of strong free-free variability or anomalous microwave emission
Authors:
Pietro Curone,
Leonardo Testi,
Enrique Macias,
Marco Tazzari,
Stefano Facchini,
Jonathan P. Williams,
Cathie J. Clarke,
Antonella Natta,
Giovanni Rosotti,
Claudia Toci,
Giuseppe Lodato
Abstract:
Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosp…
▽ More
Protoplanetary disks emit radiation across a broad range of wavelengths, requiring a multiwavelength approach to fully understand their physical mechanisms and how they form planets. Observations at sub-millimeter to centimeter wavelengths can provide insights into the thermal emission from dust, free-free emission from ionized gas, and possible gyro-synchrotron emission from the stellar magnetosphere. This work is focused on CX Tau, a ${\sim}0.4\,M_\odot$ star with an extended gas emission and a compact and apparently structureless dust disk, with an average millimeter flux when compared to Class II sources in Taurus. We present observations from the Karl G. Jansky Very Large Array (VLA) across four bands (between 9.0 mm and 6.0 cm) and combine them with archival data from the Atacama Large Millimeter/submillimeter Array (ALMA), the Submillimeter Array (SMA) and the Plateau de Bure Interferometer (PdBI). This multiwavelength approach allows us to separate the dust continuum from other emissions. After isolating the dust thermal emission, we derived an upper limit of the dust disk extent at 1.3 cm which is consistent with theoretical predictions of a radial drift-dominated disk. Centimeter data show a peculiar behavior: deep observations at 6.0 cm did not detect the source, while at 1.3 cm the flux density is anomalously higher than adjacent bands. Intraband spectral indices suggest a dominant contribution from free-free emission, whereas gyro-synchrotron emission is excluded. To explain these observations, we propose a strong variability among the free-free emission with timescales shorter than a month. Another possible interpretation is the presence of anomalous microwave emission from spinning dust grains.
△ Less
Submitted 22 August, 2023; v1 submitted 20 July, 2023;
originally announced July 2023.
-
The role of the drag force in the gravitational stability of dusty planet-forming disc -- II. Numerical simulations
Authors:
Cristiano Longarini,
Philip J. Armitage,
Giuseppe Lodato,
Daniel J. Price,
Simone Ceppi
Abstract:
Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed p…
▽ More
Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St>1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time, and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio, and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1-10 Earth masses. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity.
△ Less
Submitted 5 May, 2023;
originally announced May 2023.
-
A Large Double-ring Disk around the Taurus M Dwarf J04124068+2438157
Authors:
Feng Long,
Bin B. Ren,
Nicole L. Wallack,
Daniel Harsono,
Gregory J. Herczeg,
Paola Pinilla,
Dimitri Mawet,
Michael C. Liu,
Sean M. Andrews,
Xue-Ning Bai,
Sylvie Cabrit,
Lucas A. Cieza,
Doug Johnstone,
Jarron M. Leisenring,
Giuseppe Lodato,
Yao Liu,
Carlo F. Manara,
Gijs D. Mulders,
Enrico Ragusa,
Steph Sallum,
Yangfan Shi,
Marco Tazzari,
Taichi Uyama,
Kevin Wagner,
David J. Wilner
, et al. (1 additional authors not shown)
Abstract:
Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcse…
▽ More
Planet formation imprints signatures on the physical structures of disks. In this paper, we present high-resolution ($\sim$50 mas, 8 au) Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and CO line emission toward the disk around the M3.5 star 2MASS J04124068+2438157. The dust disk consists only of two narrow rings at radial distances of 0.47 and 0.78 arcsec ($\sim$70 and 116 au), with Gaussian $σ$ widths of 5.6 and 8.5 au, respectively. The width of the outer ring is smaller than the estimated pressure scale height by $\sim25\%$, suggesting dust trapping in a radial pressure bump. The dust disk size, set by the location of the outermost ring, is significantly larger (by $3σ$) than other disks with similar millimeter luminosity, which can be explained by an early formation of local pressure bump to stop radial drift of millimeter dust grains. After considering the disk's physical structure and accretion properties, we prefer planet--disk interaction over dead zone or photoevaporation models to explain the observed dust disk morphology. We carry out high-contrast imaging at $L'$ band using Keck/NIRC2 to search for potential young planets, but do not identify any source above $5σ$. Within the dust gap between the two rings, we reach a contrast level of $\sim$7 mag, constraining the possible planet below $\sim$2--4 $M_{\rm Jup}$. Analyses of the gap/ring properties suggest a $\sim$Saturn mass planet at $\sim$90 au is likely responsible for the formation of the outer ring, which can be potentially revealed with JWST.
△ Less
Submitted 25 March, 2023;
originally announced March 2023.
-
Kinematics signature of a giant planet in the disk of AS 209
Authors:
D. Fedele,
F. Bollati,
G. Lodato
Abstract:
[abridged] ALMA observations of dust in protoplanetary disks are revealing the existence of sub-structures such as rings, gaps and cavities. Such morphology are expected to be the outcome of dynamical interaction between the disk and planets. However, other mechanisms are able to produce similar dust sub-structures. A solution is to look at the perturbation induced by the planet to the gas surface…
▽ More
[abridged] ALMA observations of dust in protoplanetary disks are revealing the existence of sub-structures such as rings, gaps and cavities. Such morphology are expected to be the outcome of dynamical interaction between the disk and planets. However, other mechanisms are able to produce similar dust sub-structures. A solution is to look at the perturbation induced by the planet to the gas surface density and/or to the kinematics. In the case of the disk around AS 209, a prominent gap has been reported in the surface density of CO at $r \sim 100\,$au. Recently, Bae et al. (2022) detected a localized velocity perturbation in the $^{12}$CO $J=2-1$ emission along with a clump in $^{13}$CO $J=2-1$ at nearly 200 au, interpreted as a gaseous circumplanetary disk. We report a new analysis of ALMA archival observations of $^{12}$CO and $^{13}$CO J=2-1. A clear kinematics perturbation (kink) is detected in multiple channels and over a wide azimuth range in both dataset. We compared the observed perturbation with a semi-analytic model of velocity perturbations due to planet-disk interaction. The observed kink is not consistent with a planet at 200\,au as this would require a low gas disk scale height ($< 0.05$) in contradiction with previous estimate ($h/r \sim 0.118$ at $r = 100$ au). When we fix the disk scale height to 0.118 (at $r = 100$ au) we find instead that a planet of 3-5 M$_{\rm Jup}$ at 100 au induces a kinematics perturbation similar to the observed one. Thus, we conclude that a giant protoplanet orbiting at $r \sim 100\,$au is responsible of the large scale kink as well as of the perturbed dust and gas surface density previously detected. The position angle of the planet is constrained to be between 60$^{\circ}$-100$^{\circ}$. Future observations with high contrast imaging technique in the near- and mid- infrared are needed to confirm the presence and position of such a planet.
△ Less
Submitted 21 February, 2023;
originally announced February 2023.
-
Precession and polar alignment of accretion discs in triple (or multiple) stellar systems
Authors:
Simone Ceppi,
Cristiano Longarini,
Giuseppe Lodato,
Nicolás Cuello,
Stephen H. Lubow
Abstract:
We investigate the mechanism of polar alignment for accretion discs in hierarchical systems (HSs) with more than two stars. In eccentric binary systems, low mass discs that are sufficiently tilted to the binary orbit align in a polar configuration with respect to the binary plane by aligning their angular momentum to the binary eccentricity vector. In HSs, secular evolution of the orbital paramete…
▽ More
We investigate the mechanism of polar alignment for accretion discs in hierarchical systems (HSs) with more than two stars. In eccentric binary systems, low mass discs that are sufficiently tilted to the binary orbit align in a polar configuration with respect to the binary plane by aligning their angular momentum to the binary eccentricity vector. In HSs, secular evolution of the orbital parameters makes the eccentricity vector of the system precess with time. This precession undermines the stability of the polar orbit for accretion discs hosted in HSs. We analytically show that the binary criteria for polar alignment derived in the literature are necessary but not sufficient conditions for polar alignment in HSs. Then, we derive an analytical criterion for polar alignment in HSs. In general, we find that discs orbiting the innermost level of a HS can go polar. Conversely, radially extended discs orbiting the outer levels of a HS cannot polarly align and evolve as orbiting around a circular binary. We confirm our findings through detailed numerical simulations. Also, our results are compatible with the observed distribution of disc-orbit mutual inclination. Finally, we compare the observed distribution of disc inclinations in the binary and in the HS populations. Binaries host mainly coplanar discs, while HSs show a wide range of disc inclinations. We suggest that the wider range of inclinations in HSs results from the secular oscillation of their orbital parameters (such as Kozai-Lidov oscillations), rather than from a different initial condition or evolution between HSs and binaries.
△ Less
Submitted 7 February, 2023;
originally announced February 2023.
-
Disk Evolution Study Through Imaging of Nearby Young Stars (DESTINYS): Characterization of the young star T CrA and its circumstellar environment
Authors:
E. Rigliaco,
R. Gratton,
S. Ceppi,
C. Ginski,
M. Hogerheijde,
M. Benisty,
T. Birnstiel,
M. Dima,
S. Facchini,
A. Garufi,
J. Bae,
M. Langlois,
G. Lodato,
E. Mamajek,
C. F. Manara,
F. Ménard,
Á. Ribas,
A. Zurlo
Abstract:
Birth environments of young stars have strong imprints on the star itself and their surroundings. We present a detailed analysis of the wealthy circumstellar environment around the young Herbig Ae/Be star TCrA. Our aim is to understand the nature of the stellar system and the extended circumstellar structures as seen in scattered light images. We conduct our analysis combining archival data, and n…
▽ More
Birth environments of young stars have strong imprints on the star itself and their surroundings. We present a detailed analysis of the wealthy circumstellar environment around the young Herbig Ae/Be star TCrA. Our aim is to understand the nature of the stellar system and the extended circumstellar structures as seen in scattered light images. We conduct our analysis combining archival data, and new adaptive optics high-contrast and high-resolution images. The scattered light images reveal the presence of a complex environment composed of a bright forward scattering rim of the disk's surface that is seen at very high inclination, a dark lane of the disk midplane, bipolar outflows, and streamer features likely tracing infalling material from the surrounding birth cloud onto the disk. The analysis of the light curve suggests the star is a binary with a period of 29.6yrs. The comparison of the scattered light images with ALMA continuum and 12CO line emission shows the disk is in keplerian rotation, with the northern side of the outflowing material receding, while the southern side approaching the observer. The disk is itself seen edge-on. The direction of the outflows seen in scattered light is in agreement with the direction of the more distant molecular hydrogen emission-line objects (MHOs) associated to the star. Modeling of the SED using a radiative transfer scheme well agrees with the proposed configuration, as well as the hydrodynamical simulation performed using a Smoothed Particle Hydrodynamics code. We find evidence of streamers of accreting material around TCrA. These streamers connect the filament along which TCrA is forming with the outer parts of the disk, suggesting that the strong misalignment between the inner and outer disk is due to a change in the direction of the angular momentum of the material accreting on the disk during the late phase of star formation.
△ Less
Submitted 4 January, 2023;
originally announced January 2023.
-
The role of the drag force in the gravitational stability of dusty planet forming disc -- I. Analytical theory
Authors:
Cristiano Longarini,
Giuseppe Lodato,
Giuseppe Bertin,
Philip J. Armitage
Abstract:
Recent observations show that planet formation is already underway in young systems, when the protostar is still embedded into the molecular cloud and the accretion disc is massive. In such environments, the role of self gravity (SG) and gravitational instability (GI) is crucial in determining the dynamical evolution of the disc. In this work, we study the dynamical role of drag force in self-grav…
▽ More
Recent observations show that planet formation is already underway in young systems, when the protostar is still embedded into the molecular cloud and the accretion disc is massive. In such environments, the role of self gravity (SG) and gravitational instability (GI) is crucial in determining the dynamical evolution of the disc. In this work, we study the dynamical role of drag force in self-gravitating discs as a way to form planetesimals in early protoplanetary stages. We obtain the dispersion relation for density-wave perturbations on a fluid composed of two phases (gas and dust) interacting through the common gravitation field and the mutual drag force, and we find that the stability threshold is determined by three parameters: the local dust-to-gas density ratio, the dust relative temperature and the relevant Stokes number. In a region of parameters space, where young protoplanetary discs are likely to be found, the instability can be dust driven, occurring at small wavelengths. In this regime, the Jeans mass is much smaller than the one predicted by the standard gravitational instability model. This mechanism can be a viable way to form planetary cores in protostellar discs, since their predicted mass is about 10 Earth masses.
△ Less
Submitted 9 December, 2022;
originally announced December 2022.
-
Dynamical mass measurements of two protoplanetary discs
Authors:
G. Lodato,
L. Rampinelli,
E. Viscardi,
C. Longarini,
A. Inzquierdo,
T. Paneque-Carreno,
L. Testi,
S. Facchini,
A. Miotello,
B. Veronesi,
C. Hall
Abstract:
ALMA observations of line emission from planet forming discs have demonstrated to be an excellent tool to probe the internal disc kinematics, often revealing subtle effects related to important dynamical processes occurring in them, such as turbulence, or the presence of planets, that can be inferred from pressure bumps perturbing the gas motion, or from detection of the planetary wake. In particu…
▽ More
ALMA observations of line emission from planet forming discs have demonstrated to be an excellent tool to probe the internal disc kinematics, often revealing subtle effects related to important dynamical processes occurring in them, such as turbulence, or the presence of planets, that can be inferred from pressure bumps perturbing the gas motion, or from detection of the planetary wake. In particular, we have recently shown for the case of the massive disc in Elias 2-27 how one can use such kind of observations to measure deviations from Keplerianity induced by the disc self-gravity, thus constraining the total disc mass with good accuracy and independently on mass conversion factors between the tracer used and the total mass. Here, we refine our methodology and extend it to two additional sources, GM Aur and IM Lup, for which archival line observations are available for both the 12CO and the 13CO line. For IM Lup, we are able to obtain a consistent disc mass of Mdisc=0.1 Msun, implying a disc-star mass ratio of 0.1 (consistent with the observed spiral structure in the continuum emission) and a gas/dust ratio of ~ 65 (consistent with standard assumptions), with a systematic uncertainty by a factor ~2 due to the different methods to extract the rotation curve. For GM Aur, the two lines we use provide slightly inconsistent rotation curves, that cannot be attributed only to a difference in the height of the emitting layer, nor to a vertical temperature stratification. Our best fit disc mass measurement is Mdisc=0.26Msun, implying a disc-star mass ratio of ~0.35 and a gas/dust ratio of ~130... ABRIDGED
△ Less
Submitted 7 November, 2022;
originally announced November 2022.
-
An analytic solution to measure the gas size in protoplanetary discs in the viscous self-similar scenario
Authors:
Claudia Toci,
Giuseppe Lodato,
Francesco Gerardo Livio,
Giovanni Rosotti,
Leon Trapman
Abstract:
In order to understand which mechanism is responsible for accretion in protoplanetary discs, a robust knowledge of the observed disc radius using gas tracers such as $^{12}$CO and other CO isotopologues is pivotal. Indeed, the two main theories proposed, viscous accretion and wind-driven accretion, predict different time evolution for the disc radii. In this Letter, we present an analytical soluti…
▽ More
In order to understand which mechanism is responsible for accretion in protoplanetary discs, a robust knowledge of the observed disc radius using gas tracers such as $^{12}$CO and other CO isotopologues is pivotal. Indeed, the two main theories proposed, viscous accretion and wind-driven accretion, predict different time evolution for the disc radii. In this Letter, we present an analytical solution for the evolution of the disc radii in viscously evolving protoplanetary discs using $^{12}$CO as a tracer, under the assumption that the $^{12}$CO radius is the radius where the surface density of the disc is equal to the threshold for CO photo-dissociation. We discuss the properties of the solution and the limits of its applicability as a simple numerical prescription to evaluate the observed disc radii of populations of discs. Our results suggest that, in addition to photo-dissociation, also freeze out plays an important role in setting the disc size. We find an effective reduction of the CO abundance by about two orders of magnitude at the location of CO photo-dissociation, which however should not be interpreted as the bulk abundance of CO in the disc. The use of our analytical solution will allow to compute disc sizes for large quantities of models without using expensive computational resources such as radiative transfer calculations.
△ Less
Submitted 10 November, 2022; v1 submitted 26 October, 2022;
originally announced October 2022.
-
H$_2$O distribution in the disc of HD 100546 and HD 163296: the role of dust dynamics and planet--disc interaction
Authors:
L. M. Pirovano,
D. Fedele,
E. F. van Dishoeck,
M. R. Hogerheijde,
G. Lodato,
S. Bruderer
Abstract:
[Abridged] Far-infrared observations with Herschel revealed a surprisingly low abundance of cold-water reservoirs in protoplanetary discs. On the other hand, a handful of discs show emission of hot water transitions excited at temperatures above a few hundred Kelvin. In particular, the protoplanetary discs around the Herbig Ae stars HD 100546 and HD 163296 show opposite trends in terms of cold ver…
▽ More
[Abridged] Far-infrared observations with Herschel revealed a surprisingly low abundance of cold-water reservoirs in protoplanetary discs. On the other hand, a handful of discs show emission of hot water transitions excited at temperatures above a few hundred Kelvin. In particular, the protoplanetary discs around the Herbig Ae stars HD 100546 and HD 163296 show opposite trends in terms of cold versus hot water emission: in the first case, the ground-state transitions are detected and the high-J lines are undetected, while the trend is opposite in HD 163296. We performed a spectral analysis using the thermo-chemical model DALI. We find that HD 163296 is characterised by a water-rich (abundance $\gtrsim 10^{-5}$) hot inner disc (within the snowline) and a water-poor ($< 10^{-10}$) outer disc: the relative abundance may be due to the thermal desorption of icy grains that have migrated inward. Remarkably, the size of the H$_2$O emitting region corresponds to a narrow dust gap visible in the millimeter continuum at $r=10\,$au with ALMA. The low-J lines detected in HD 100546 instead imply an abundance of a few $10^{-9}$ in the cold outer disc ($> 40$ au). The emitting region of the cold H$_2$O transitions is spatially coincident with that of the H$_2$O ice previously seen in the near-infrared. Notably, millimetre observations with ALMA reveal the presence of a large dust gap between nearly 40 and 150 au, likely opened by a massive embedded protoplanet. In both discs, we find that the warm molecular layer in the outer region (beyond the snow line) is highly depleted of water molecules, implying an oxygen-poor chemical composition of the gas. We speculate that gas-phase oxygen in the outer disc is readily depleted and its distribution in the disc is tightly coupled to the dynamics of the dust grains.
△ Less
Submitted 21 July, 2022;
originally announced July 2022.
-
External or internal companion exciting the spiral arms in CQ Tau?
Authors:
Iain Hammond,
Valentin Christiaens,
Daniel J. Price,
Maria Giulia Ubeira-Gabellini,
Jennifer Baird,
Josh Calcino,
Myriam Benisty,
Giuseppe Lodato,
Leonardo Testi,
Christophe Pinte,
Claudia Toci,
Davide Fedele
Abstract:
We present new high-contrast images in near-infrared wavelengths ($λ$ = 1.04, 1.24, 1.62, 2.18 and 3.78$μ$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M…
▽ More
We present new high-contrast images in near-infrared wavelengths ($λ$ = 1.04, 1.24, 1.62, 2.18 and 3.78$μ$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M$_{\rm Jup}$ immediately outside the spiral arms at $\sim$60 au and above 2-3 M$_{\rm Jup}$ beyond 100 au to 5$σ$ confidence assuming hot-start models. We do, however, detect four spiral arms, a double-arc and evidence for shadows in scattered light cast by a misaligned inner disc. Our observations may be explained by an unseen close-in companion on an inclined and eccentric orbit. Such a hypothesis would also account for the disc CO cavity and disturbed kinematics.
△ Less
Submitted 31 July, 2022; v1 submitted 18 July, 2022;
originally announced July 2022.
-
On the time evolution of the $M_{\rm d} - M_\star$ and $\dot M - M_\star$ correlations for protoplanetary discs: the viscous timescale increases with stellar mass
Authors:
Alice Somigliana,
Claudia Toci,
Giovanni Rosotti,
Giuseppe Lodato,
Marco Tazzari,
Carlo Manara,
Leonardo Testi,
Federico Lepri
Abstract:
Large surveys of star-forming regions have unveiled power-law correlations between the stellar mass and the disc parameters, such as the disc mass $M_{\mathrm{d}} \propto {M_{\star}}^{λ_{\mathrm{m}}}$ and the accretion rate $\dot M \propto {M_{\star}}^{λ_{\mathrm{acc}}}$. The observed slopes appear to be increasing with time, but the reason behind the establishment of these correlations and their…
▽ More
Large surveys of star-forming regions have unveiled power-law correlations between the stellar mass and the disc parameters, such as the disc mass $M_{\mathrm{d}} \propto {M_{\star}}^{λ_{\mathrm{m}}}$ and the accretion rate $\dot M \propto {M_{\star}}^{λ_{\mathrm{acc}}}$. The observed slopes appear to be increasing with time, but the reason behind the establishment of these correlations and their subsequent evolution is still uncertain. We conduct a theoretical analysis of the impact of viscous evolution on power-law initial conditions for a population of protoplanetary discs. We find that, for evolved populations, viscous evolution enforces the two correlations to have the same slope, $λ_{\mathrm{m}}$ = $λ_{\mathrm{acc}}$, and that this limit is uniquely determined by the initial slopes $λ_{\mathrm{m}, 0}$ and $λ_{\mathrm{acc}, 0}$. We recover the increasing trend claimed from the observations when the difference in the initial values, $δ_0 = λ_{\mathrm{m}, 0} - λ_{\mathrm{acc}, 0}$, is larger than $1/2$; moreover, we find that this increasing trend is a consequence of a positive correlation between the viscous timescale and the stellar mass. We also present the results of disc population synthesis numerical simulations, that allow us to introduce a spread and analyse the effect of sampling, which show a good agreement with our analytical predictions. Finally, we perform a preliminary comparison of our numerical results with observational data, which allows us to constrain the parameter space of the initial conditions to $λ_{\mathrm{m}, 0} \in [1.2, 2.1]$, $λ_{\mathrm{acc}, 0} \in [0.7, 1.5]$.
△ Less
Submitted 8 June, 2022;
originally announced June 2022.
-
A giant planet shaping the disk around the very low-mass star CIDA 1
Authors:
P. Curone,
A. F. Izquierdo,
L. Testi,
G. Lodato,
S. Facchini,
A. Natta,
P. Pinilla,
N. T. Kurtovic,
C. Toci,
M. Benisty,
M. Tazzari,
F. Borsa,
M. Lombardi,
C. F. Manara,
E. Sanchis,
L. Ricci
Abstract:
(Abridged) Exoplanetary research has provided us with exciting discoveries of planets around very low-mass (VLM) stars (e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from ALMA of the disk around CIDA 1, a VLM star in Taurus, s…
▽ More
(Abridged) Exoplanetary research has provided us with exciting discoveries of planets around very low-mass (VLM) stars (e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from ALMA of the disk around CIDA 1, a VLM star in Taurus, show substructures hinting at the presence of a massive planet. We aim to reproduce the dust ring of CIDA 1, observed in the dust continuum emission in ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm), along with its $^{12}$CO (J=3-2) and $^{13}$CO (J=3-2) channel maps, assuming the structures are shaped by the interaction of the disk with a massive planet. We seek to retrieve the mass and position of the putative planet. We model the protoplanetary disk with a set of hydrodynamical simulations, varying the mass and locations of the embedded planet. We compute the dust and gas emission using radiative transfer simulations, and, finally, we obtain the synthetic observations treating the images as the actual ALMA observations. Our models indicate that a planet with a minimum mass of $\sim1.4\,\text{M}_\text{Jup}$ orbiting at a distance of $\sim 9-10$ au can explain the morphology and location of the observed dust ring at Band 7 and Band 4. We can reproduce the low spectral index ($\sim 2$) observed where the dust ring is detected. Our synthetic images reproduce the morphology of the $^{12}$CO and $^{13}$CO observed channel maps where the cloud absorption allowed a detection. Applying an empirical relation between planet mass and gap width in the dust, we predict a maximum planet mass of $\sim4 - 8\,\text{M}_\text{Jup}$. Our results suggest the presence of a massive planet orbiting CIDA 1, thus challenging our understanding of planet formation around VLM stars.
△ Less
Submitted 17 August, 2022; v1 submitted 20 May, 2022;
originally announced May 2022.
-
Accretion rates in hierarchical triple systems with discs
Authors:
Simone Ceppi,
Nicolás Cuello,
Giuseppe Lodato,
Cathie Clarke,
Claudia Toci,
Daniel J. Price
Abstract:
Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical…
▽ More
Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical configurations. Despite its astrophysical relevance, differential accretion in hierarchical systems remains to be understood. In this work, we investigate whether the accretion trends expected in binaries are valid for higher order multiples. We performed a set of 3D Smoothed Particle Hydrodynamics simulations of binaries and of hierarchical triples (HTs) embedded in an accretion disc, with the code Phantom. We identify for the first time accretion trends in HTs and their deviations compared to binaries. These deviations, due to the interaction of the small binary with the infalling material from the circum-triple disc, can be described with a semi-analytical prescription. Generally, the smaller binary of a HT accretes more mass than a single star of the same mass as the smaller binary. We found that in a HT, if the small binary is heavier than the third body, the standard differential accretion scenario (whereby the secondary accretes more of the mass) is hampered. Reciprocally, if the small binary is lighter than the third body, the standard differential accretion scenario is enhanced. The peculiar differential accretion mechanism we find in HTs is expected to affect their mass ratio distribution.
△ Less
Submitted 18 May, 2022;
originally announced May 2022.
-
Extreme Mass Ratio Inspirals triggered by Massive Black Hole Binaries: from Relativistic Dynamics to Cosmological Rates
Authors:
Giovanni Mazzolari,
Matteo Bonetti,
Alberto Sesana,
Riccardo M. Colombo,
Massimo Dotti,
Giuseppe Lodato,
David Izquierdo-Villalba
Abstract:
Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass-ratio $q=m/M$ in the range $~10^{-9}-10^{-4}$ and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the centra…
▽ More
Extreme mass ratio inspirals (EMRIs) are compact binary systems characterized by a mass-ratio $q=m/M$ in the range $~10^{-9}-10^{-4}$ and represent primary gravitational wave (GW) sources for the forthcoming Laser Interferometer Space Antenna (LISA). While their standard formation channel involves relaxation processes deflecting compact objects on very low angular momentum orbits around the central massive black hole, a number of alternative formation channels has been proposed, including binary tidal break-up, migration in accretion disks and secular and chaotic dynamics around a massive black hole binary (MBHB). In this work, we take an extensive closer look at this latter scenario, investigating how EMRIs can be triggered by a MBHBs, formed in the aftermath of galaxy mergers. By employing a suite of relativistic three-body simulations, we evaluate the efficiency of EMRI formation for different parameters of the MBHB, assessing the importance of both secular and chaotic dynamics. By modelling the distribution of compact objects in galaxy nuclei, we estimate the resulting EMRI formation rate, finding that EMRI are produced in a sharp burst, with peak rates that are 10-100 times higher than the standard two-body relaxation channel, lasting for 10$^6$--10$^8$ years. By coupling our results with an estimate of the cosmic MBHB merger rate, we finally forecast that LISA could observe ${\cal O}(10)$ EMRIs per year formed by this channel.
△ Less
Submitted 8 August, 2022; v1 submitted 11 April, 2022;
originally announced April 2022.
-
Demographics of young stars and their protoplanetary disks: lessons learned on disk evolution and its connection to planet formation
Authors:
C. F. Manara,
M. Ansdell,
G. P. Rosotti,
A. M. Hughes,
P. J. Armitage,
G. Lodato,
J. P. Williams
Abstract:
Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys…
▽ More
Since Protostars and Planets VI (PPVI), our knowledge of the global properties of protoplanetary and debris disks, as well as of young stars, has dramatically improved. At the time of PPVI, mm-observations and optical to near-infrared spectroscopic surveys were largely limited to the Taurus star-forming region, especially of its most massive disk and stellar population. Now, near-complete surveys of multiple star-forming regions cover both spectroscopy of young stars and mm interferometry of their protoplanetary disks. This provides an unprecedented statistical sample of stellar masses and mass accretion rates, as well as disk masses and radii, for almost 1000 young stellar objects within 300 pc from us, while also sampling different evolutionary stages, ages, and environments. At the same time, surveys of debris disks are revealing the bulk properties of this class of more evolved objects. This chapter reviews the statistics of these measured global star and disk properties and discusses their constraints on theoretical models describing global disk evolution. Our comparisons of observations to theoretical model predictions extends beyond the traditional viscous evolution framework to include analytical descriptions of magnetic wind effects. Finally, we discuss how recent observational results can provide a framework for models of planet population synthesis and planet formation.
△ Less
Submitted 1 March, 2023; v1 submitted 18 March, 2022;
originally announced March 2022.
-
Astrophysics with the Laser Interferometer Space Antenna
Authors:
Pau Amaro Seoane,
Jeff Andrews,
Manuel Arca Sedda,
Abbas Askar,
Quentin Baghi,
Razvan Balasov,
Imre Bartos,
Simone S. Bavera,
Jillian Bellovary,
Christopher P. L. Berry,
Emanuele Berti,
Stefano Bianchi,
Laura Blecha,
Stephane Blondin,
Tamara Bogdanović,
Samuel Boissier,
Matteo Bonetti,
Silvia Bonoli,
Elisa Bortolas,
Katelyn Breivik,
Pedro R. Capelo,
Laurentiu Caramete,
Federico Cattorini,
Maria Charisi,
Sylvain Chaty
, et al. (134 additional authors not shown)
Abstract:
The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery…
▽ More
The Laser Interferometer Space Antenna (LISA) will be a transformative experiment for gravitational wave astronomy, and, as such, it will offer unique opportunities to address many key astrophysical questions in a completely novel way. The synergy with ground-based and space-born instruments in the electromagnetic domain, by enabling multi-messenger observations, will add further to the discovery potential of LISA. The next decade is crucial to prepare the astrophysical community for LISA's first observations. This review outlines the extensive landscape of astrophysical theory, numerical simulations, and astronomical observations that are instrumental for modeling and interpreting the upcoming LISA datastream. To this aim, the current knowledge in three main source classes for LISA is reviewed; ultracompact stellar-mass binaries, massive black hole binaries, and extreme or intermediate mass ratio inspirals. The relevant astrophysical processes and the established modeling techniques are summarized. Likewise, open issues and gaps in our understanding of these sources are highlighted, along with an indication of how LISA could help making progress in the different areas. New research avenues that LISA itself, or its joint exploitation with upcoming studies in the electromagnetic domain, will enable, are also illustrated. Improvements in modeling and analysis approaches, such as the combination of numerical simulations and modern data science techniques, are discussed. This review is intended to be a starting point for using LISA as a new discovery tool for understanding our Universe.
△ Less
Submitted 25 May, 2023; v1 submitted 11 March, 2022;
originally announced March 2022.
-
Recurrent X-ray flares of the black hole candidate in the globular cluster RZ 2109 in NGC 4472
Authors:
A. Tiengo,
P. Esposito,
M. Toscani,
G. Lodato,
M. Arca Sedda,
S. E. Motta,
F. Contato,
M. Marelli,
R. Salvaterra,
A. De Luca
Abstract:
We report on the systematic analysis of the X-ray observations of the ultra-luminous X-ray source XMMU J122939.7+075333 located in the globular cluster RZ 2109 in the Virgo galaxy NGC 4472. The inclusion of observations and time intervals ignored in previous works and the careful selection of extraction regions and energy bands have allowed us to identify new flaring episodes, in addition to the o…
▽ More
We report on the systematic analysis of the X-ray observations of the ultra-luminous X-ray source XMMU J122939.7+075333 located in the globular cluster RZ 2109 in the Virgo galaxy NGC 4472. The inclusion of observations and time intervals ignored in previous works and the careful selection of extraction regions and energy bands have allowed us to identify new flaring episodes, in addition to the ones that made it one of the best black hole candidates in globular clusters. Although most observations are too short and sparse to recognize a regular pattern, the spacing of the three most recent X-ray flares is compatible with a ~34 hours recurrence time. If confirmed by future observations, such behavior, together with the soft spectrum of the X-ray flares, would be strikingly similar to the quasi-periodic eruptions recently discovered in galactic nuclei. Following one of the possible interpretations of these systems and of a peculiar class of extra-galactic X-ray transients, we explore the possibility that XMMU J122939.7+075333 might be powered by the partial disruption of a white dwarf by an intermediate mass (M~700 Msun) black hole.
△ Less
Submitted 17 February, 2022;
originally announced February 2022.
-
The protoplanetary disk population in the rho-Ophiuchi region L1688 and the time evolution of Class II YSOs
Authors:
L. Testi,
A. Natta,
C. F. Manara,
I. de Gregorio Monsalvo,
G. Lodato,
C. Lopez,
K. Muzic,
I. Pascucci,
E. Sanchis,
A. Santamaria Miranda,
A. Scholz,
M. De Simone,
J. P. Williams
Abstract:
(Abridged) We present a study of the disk population in L1688, the densest and youngest region in Ophiuchus, and we compare it with other nearby regions of different age, namely Lupus, Chamaeleon I, Corona Australis, Taurus and Upper Scorpius. We select our L1688 sample using a combination of criteria (ALMA data, Gaia, optical/near-IR spectroscopy) and determine stellar and disk properties, specif…
▽ More
(Abridged) We present a study of the disk population in L1688, the densest and youngest region in Ophiuchus, and we compare it with other nearby regions of different age, namely Lupus, Chamaeleon I, Corona Australis, Taurus and Upper Scorpius. We select our L1688 sample using a combination of criteria (ALMA data, Gaia, optical/near-IR spectroscopy) and determine stellar and disk properties, specifically stellar mass (Mstar), average population age, mass accretion rate (Macc) and disk dust mass (Mdust). a) In L1688 the relations between Macc and Mstar, Mdust and Mstar, and Macc and Mdust have a roughly linear trend with slopes 1.8-1.9 for the first two relations and ~1 for the third, similarly to what found in the other regions. b) When ordered according to the characteristic age of each region, Macc decreases as 1/t, when corrected for the different stellar mass content; Mdust follows roughly the same trend between 0.5 and 5 Myr, but has an increase of a factor ~3 at ages of 2-3 Myr. We suggest that this could result from an earlier planet formation, followed by collisional fragmentation that temporarily replenishes the millimeter-size grain population. c) The dispersion of Macc and Mdust around the best-fitting relation with Mstar, as well as that of Macc versus Mdust are large: we find that the dispersions have continuous distributions with a log-normal shape and similar width (~0.8 dex). The amount of dust observed at ~1 Myr does not appear to be sufficient to assemble the majority of planetary systems, which suggests an earlier planetary cores formation. The dust mass traces to a large extent the disk gas mass evolution. Two properties remain puzzling: the steep dependence of Macc and Mdust on Mstar and the cause of the large dispersion in the three relations analyzed in this paper, in particular the one of the Macc versus Mdust relation.
△ Less
Submitted 17 May, 2022; v1 submitted 11 January, 2022;
originally announced January 2022.
-
Observational constraints on disc sizes in protoplanetary discs in multiple systems in the Taurus region. II. Gas disc sizes
Authors:
A. A. Rota,
C. F. Manara,
A. Miotello,
G. Lodato,
S. Facchini,
M. Koutoulaki,
G. Herczeg,
F. Long,
M. Tazzari,
S. Cabrit,
D. Harsono,
F. Menard,
P. Pinilla,
G. van der Plas,
E. Ragusa,
H. -W. Yen
Abstract:
The formation of multiple stellar systems is a natural by-product of the star-formation process, and its impact on the properties of protoplanetary discs and on the formation of planets is still to be fully understood. To date, no detailed uniform study of the gas emission from a sample of protoplanetary discs around multiple stellar systems has been performed. Here we analyse new ALMA observation…
▽ More
The formation of multiple stellar systems is a natural by-product of the star-formation process, and its impact on the properties of protoplanetary discs and on the formation of planets is still to be fully understood. To date, no detailed uniform study of the gas emission from a sample of protoplanetary discs around multiple stellar systems has been performed. Here we analyse new ALMA observations at a $\sim$21 au resolution of the molecular CO gas emission targeting discs in eight multiple stellar systems in the Taurus star-forming regions. $^{12}$CO gas emission is detected around all primaries and in seven companions. With these data, we estimate the inclination and the position angle for all primary discs and for five secondary or tertiary discs, and measure the gas disc radii of these objects with a cumulative flux technique on the spatially resolved zeroth moment images. When considering the radius including 95\% of the flux as a metric, the estimated gas disc size in multiple stellar systems is found to be on average $\sim 4.2$ times larger than the dust disc size. This ratio is higher than what was recently found in a population of more isolated and single systems. On the contrary, when considering the radius including 68\% of the flux, no difference between multiple and single discs is found in the distribution of ratios. This discrepancy is due to the sharp truncation of the outer dusty disc observed in multiple stellar systems. The measured gas disc sizes are consistent with tidal truncation models in multiple stellar systems assuming eccentricities of $\sim0.15$-$0.5$, as expected in typical binary systems.
△ Less
Submitted 10 January, 2022;
originally announced January 2022.
-
MHD disc winds can reproduce fast disc dispersal and the correlation between accretion rate and disc mass in Lupus
Authors:
Benoît Tabone,
Giovanni P. Rosotti,
Giuseppe Lodato,
Philip J. Armitage,
Alexander J. Cridland,
Ewine F. van Dishoeck
Abstract:
The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solut…
▽ More
The final architecture of planetary systems depends on the extraction of angular momentum and mass-loss processes of the discs in which they form. Theoretical studies proposed that magnetized winds launched from the discs (MHD disc winds) could govern accretion and disc dispersal. In this work, we revisit the observed disc demographics in the framework of MHD disc winds, combining analytical solutions of disc evolution and a disc population synthesis approach. We show that MHD disc winds alone can account for both disc dispersal and accretion properties. The decline of disc fraction over time is reproduced by assuming that the initial accretion timescale (a generalization of the viscous timescale) varies from disc to disc and that the decline of the magnetic field strength is slower than that of the gas. The correlation between accretion rate and disc mass, and the dispersion of the data around the mean trend as observed in Lupus is then naturally reproduced. The model also accounts for the rapidity of the disc dispersal. This paves the way for planet formation models in the paradigm of wind-driven accretion.
△ Less
Submitted 29 November, 2021;
originally announced November 2021.
-
Fully-discrete spatial eigenanalysis of discontinuous spectral element methods: insights into well-resolved and under-resolved vortical flows
Authors:
Niccolò Tonicello,
Rodrigo C Moura,
Guido Lodato,
Gianmarco Mengaldo
Abstract:
This study presents a comprehensive spatial eigenanalysis of fully-discrete discontinuous spectral element methods, now generalizing previous spatial eigenanalysis that did not include time integration errors. The influence of discrete time integration is discussed in detail for different explicit Runge-Kutta (1st to 4th order accurate) schemes combined with either Discontinuous Galerkin (DG) or S…
▽ More
This study presents a comprehensive spatial eigenanalysis of fully-discrete discontinuous spectral element methods, now generalizing previous spatial eigenanalysis that did not include time integration errors. The influence of discrete time integration is discussed in detail for different explicit Runge-Kutta (1st to 4th order accurate) schemes combined with either Discontinuous Galerkin (DG) or Spectral Difference (SD) methods, both here recovered from the Flux Reconstruction (FR) scheme. Selected numerical experiments using the improved SD method by Liang and Jameson [1] are performed to quantify the influence of time integration errors on actual simulations. These involve test cases of varied complexity, from one-dimensional linear advection equation studies to well-resolved and under-resolved inviscid vortical flows. It is shown that, while both well-resolved and under-resolved simulations of linear problems correlate well with the eigenanalysis prediction of time integration errors, the correlation can be much worse for under-resolved nonlinear problems. The effect of mesh regularity is also considered, where time integration errors are found to be, in the case of irregular grids, less pronounced than those of the spatial discretisation. In fact, for the under-resolved vortical flows considered, the predominance of spatial errors made it practically impossible for time integration errors to be distinctly identified. Nevertheless, for well-resolved nonlinear simulations, the effect of time integration errors could still be recognized. This highlights that the interaction between space and time discretisation errors is more complex than otherwise anticipated, contributing to the current understanding about when eigenanalysis can effectively predict the behaviour of numerical errors in practical under-resolved nonlinear problems, including under-resolved turbulence computations.
△ Less
Submitted 27 November, 2021;
originally announced November 2021.
-
Secular evolution of MHD wind-driven discs: analytical solutions in the expanded α-framework
Authors:
Benoît Tabone,
Giovanni P. Rosotti,
Alexander J. Cridland,
Philip J. Armitage,
Giuseppe Lodato
Abstract:
The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscou…
▽ More
The evolution of protoplanetary discs and the related process of planet formation is regulated by angular momentum transport and mass-loss processes. Over the past decade, the paradigm of viscosity has been challenged and MHD disc winds appear as a compelling scenario to account for disc accretion. In this work, we aim to construct the equivalent of the widely used analytical description of viscous evolution for the MHD wind case. The transport of angular momentum and mass induced by the wind is parameterized by an $α$-like parameter and by the magnetic lever arm parameter $λ$. Extensions of the paradigmatic Lynden-Bell and Pringle similarity solutions to the wind case are presented. We show that wind-driven accretion leads to a steeper decrease in the disc mass and accretion rate than in viscous models due to the absence of disc spreading. If the decline of the magnetic field strength is slower than that of the gas surface density, the disc is dispersed after a finite time. The evolution of the disc in the $\dot{M}_*-M_D$ plane is sensitive to the wind and turbulence parameters. A disc population evolving under the action of winds can exhibit a correlation between $\dot{M}_*$ and $M_D$ depending on the initial conditions. The simplified framework proposed in this work opens to a new avenue to test the effectiveness of wind-driven accretion from the observed disc demographics and constitutes an important step to include wind-driven accretion in planet population synthesis models.
△ Less
Submitted 19 November, 2021;
originally announced November 2021.
-
Mapping the Planetary Wake in HD 163296 with Kinematics
Authors:
Josh Calcino,
Thomas Hilder,
Daniel J. Price,
Christophe Pinte,
Francesco Bollati,
Giuseppe Lodato,
Brodie J. Norfolk
Abstract:
We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD~163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk major axis confirm that the wake induces predominantly radial motion, as predicted by models of planet-disk interaction. These resu…
▽ More
We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD~163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk major axis confirm that the wake induces predominantly radial motion, as predicted by models of planet-disk interaction. These results provide the first direct confirmation of planet wakes generated by Lindblad resonances. Mapping the wake provides a constraint on the disk aspect ratio, which is required to measure the mass of the planet.
△ Less
Submitted 6 April, 2022; v1 submitted 14 November, 2021;
originally announced November 2021.
-
Gravitational waves from tidal disruption events: an open and comprehensive catalogue
Authors:
Martina Toscani,
Giuseppe Lodato,
Daniel J. Price,
David Liptai
Abstract:
We present an online, open and comprehensive template library of gravitational waveforms produced during the tidal disruptions of stars by massive black holes, spanning a broad space of parameters. We build this library thanks to a new feature that we implement in the general relativistic version of PHANTOM, a smoothed particle hydrodynamics code for three dimensional simulations in general relati…
▽ More
We present an online, open and comprehensive template library of gravitational waveforms produced during the tidal disruptions of stars by massive black holes, spanning a broad space of parameters. We build this library thanks to a new feature that we implement in the general relativistic version of PHANTOM, a smoothed particle hydrodynamics code for three dimensional simulations in general relativity. We first perform a series of numerical tests to show that the gravitational wave (GW) signal obtained is in excellent agreement with the one expected from theory. This benchmark is done for well studied scenarios (such as binary stellar systems). We then apply our code to calculate the GW signals from tidal disruption events (TDEs), finding that our results are consistent with the theoretical estimates obtained in previous studies for selected parameters. We illustrate interesting results from the catalogue, where we stress how the gravitational signal is affected by variations of some parameters (like black hole spin, stellar orbital eccentricity and inclination). The full catalogue is available online. It is intended to be a living catalogue.
△ Less
Submitted 9 November, 2021;
originally announced November 2021.
-
Dust Traffic Jams in Inclined Circumbinary Protoplanetary Discs I. Morphology and Formation Theory
Authors:
Hossam Aly,
Jean-François Gonzalez,
Rebecca Nealon,
Cristiano Longarini,
Giuseppe Lodato,
Daniel J. Price
Abstract:
Gas and dust in inclined orbits around binaries experience precession induced by the binary gravitational torque. The difference in precession between gas and dust alters the radial drift of weakly coupled dust and leads to density enhancements where the radial drift is minimised. We explore this phenomenon using 3D hydrodynamical simulations to investigate the prominence of these `dust traffic ja…
▽ More
Gas and dust in inclined orbits around binaries experience precession induced by the binary gravitational torque. The difference in precession between gas and dust alters the radial drift of weakly coupled dust and leads to density enhancements where the radial drift is minimised. We explore this phenomenon using 3D hydrodynamical simulations to investigate the prominence of these `dust traffic jams' and the evolution of the resulting dust sub-structures at different disc inclinations and binary eccentricities. We then derive evolution equations for the angular momentum of warped dust discs and implement them in a 1D code and present calculations to further explain these traffic jams. We find that dust traffic jams in inclined circumbinary discs provide significant dust density enhancements that are long lived and can have important consequences for planetesimal formation.
△ Less
Submitted 27 September, 2021;
originally announced September 2021.
-
Investigating Protoplanetary Disk Cooling through Kinematics: Analytical GI Wiggle
Authors:
Cristiano Longarini,
Giuseppe Lodato,
Claudia Toci,
Benedetta Veronesi,
Cassandra Hall,
Ruobing Dong,
Jason Patrick Terry
Abstract:
It is likely that young protostellar disks undergo a self-gravitating phase. Such systems are characterized by the presence of a spiral pattern that can be either in a quasi-steady state or in a nonlinear unstable condition. This spiral wave affects both the gas dynamics and kinematics, resulting in deviations from the Keplerian rotation. Recently, a lot of attention has been devoted to kinematic…
▽ More
It is likely that young protostellar disks undergo a self-gravitating phase. Such systems are characterized by the presence of a spiral pattern that can be either in a quasi-steady state or in a nonlinear unstable condition. This spiral wave affects both the gas dynamics and kinematics, resulting in deviations from the Keplerian rotation. Recently, a lot of attention has been devoted to kinematic studies of planet-forming environments, and we are now able to measure even small perturbations of velocity field ($\lesssim 1 \%$ of the Keplerian speed) thanks to high spatial and spectral resolution observations of protostellar disks. In this work, we investigate the kinematic signatures of gravitational instability: we perform an analytical study of the linear response of a self-gravitating disk to a spiral-like perturbation, focusing our attention on the velocity field perturbations. We show that unstable disks have clear kinematic imprints into the gas component across the entire disk extent, due to the GI spiral wave perturbation, resulting in deviations from Keplerian rotation. The shape of these signatures depends on several parameters, but they are significantly affected by the cooling factor: by detecting these features, we can put constraints on protoplanetary disk cooling.
△ Less
Submitted 11 January, 2022; v1 submitted 25 August, 2021;
originally announced August 2021.
-
Constraining protoplanetary disc mass using the GI wiggle
Authors:
J. P. Terry,
C. Hall,
Cristiano Longarini,
Giuseppe Lodato,
Claudia Toci,
B. Veronesi,
T. Paneque-Carreño,
C. Pinte
Abstract:
Exoplanets form in protoplanetary accretion discs. The total protoplanetary disc mass is the most fundamental parameter, since it sets the mass budget for planet formation. Although observations with the Atacama Large Millimeter/Submillimeter array (ALMA) have dramatically increased our understanding of these discs, total protoplanetary disc mass remains difficult to measure. If a disc is sufficie…
▽ More
Exoplanets form in protoplanetary accretion discs. The total protoplanetary disc mass is the most fundamental parameter, since it sets the mass budget for planet formation. Although observations with the Atacama Large Millimeter/Submillimeter array (ALMA) have dramatically increased our understanding of these discs, total protoplanetary disc mass remains difficult to measure. If a disc is sufficiency massive ($\gtrsim$ 10\% of the host star mass), it can excite gravitational instability (GI). Recently, it has been revealed that GI leaves kinematic imprints of its presence known as the "GI Wiggle." In this work, we use numerical simulations to empirically determine an approximately linear relationship between the amplitude of the wiggle and the host disc-to-star mass ratio, and show that measurements of the amplitude are possible with the spatial and spectral capabilities of ALMA. These measurements can therefore be used to constrain disc-to-star mass ratio.
△ Less
Submitted 30 November, 2021; v1 submitted 25 August, 2021;
originally announced August 2021.
-
Circumbinary and circumstellar discs around the eccentric binary IRAS 04158+2805 -- a testbed for binary-disc interaction
Authors:
Enrico Ragusa,
Daniele Fasano,
Claudia Toci,
Gaspard Duchêne,
Nicolás Cuello,
Marion Villenave,
Gerrit van der Plas,
Giuseppe Lodato,
François Ménard,
Daniel J. Price,
Christophe Pinte,
Karl Stapelfeldt,
Schuyler Wolff
Abstract:
IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$α$ imaging and Keck ast…
▽ More
IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$α$ imaging and Keck astrometric measures of the binary to develop a coherent dynamical model of this system. The system features an azimuthal asymmetry detected at the western edge of the cavity in Band~7 observations and a wiggling outflow. Dust emission in ALMA Band 4 from the proximity of the individual stars suggests the presence of marginally resolved circumstellar discs. We estimate the binary orbital parameters from the measured arc of the orbit from Keck and ALMA astrometry. We further constrain these estimates using considerations from binary-disc interaction theory. We finally perform three SPH gas + dust simulations based on the theoretical constraints; we post-process the hydrodynamic output using radiative transfer Monte Carlo methods and directly compare the models with observations. Our results suggest that a highly eccentric $e\sim 0.5\textrm{--}0.7$ equal mass binary, with a semi-major axis of $\sim 55$ au, and small/moderate orbital plane vs. circumbinary disc inclination $θ\lesssim 30^\circ$ provides a good match with observations. A dust mass of $\sim 1.5\times 10^{-4} {\rm M_\odot}$ best reproduces the flux in Band 7 continuum observations. Synthetic CO line emission maps qualitatively capture both the emission from the central region and the non-Keplerian nature of the gas motion in the binary proximity.
△ Less
Submitted 28 July, 2021;
originally announced July 2021.
-
On the secular evolution of the ratio between gas and dust radii in protoplanetary discs
Authors:
Claudia Toci,
Giovanni Rosotti,
Giuseppe Lodato,
Leonardo Testi,
Leon Trapman
Abstract:
A key problem in protoplanetary disc evolution is understanding the efficiency of dust radial drift. This process makes the observed dust disc sizes shrink on relatively short timescales, implying that discs started much larger than what we see now. In this paper we use an independent constraint, the gas radius (as probed by CO rotational emission), to test disc evolution models. In particular, we…
▽ More
A key problem in protoplanetary disc evolution is understanding the efficiency of dust radial drift. This process makes the observed dust disc sizes shrink on relatively short timescales, implying that discs started much larger than what we see now. In this paper we use an independent constraint, the gas radius (as probed by CO rotational emission), to test disc evolution models. In particular, we consider the ratio between the dust and gas radius, $R_{\rm CO}/R_{\rm dust}$. We model the time evolution of protoplanetary discs under the influence of viscous evolution, grain growth, and radial drift. Then, using the radiative transfer code RADMC with approximate chemistry, we compute the dust and gas radii of the models and investigate how $R_{\rm CO}/R_{\rm dust}$ evolves. Our main finding is that, for a broad range of values of disc mass, initial radius, and viscosity, $R_{\rm CO}/R_{\rm dust}$ becomes large (>5) after only a short time (<1 Myr) due to radial drift. This is at odds with measurements in young star forming regions such as Lupus, which find much smaller values, implying that dust radial drift is too efficient in these models. Substructures, commonly invoked to stop radial drift in large, bright discs, must then be present, although currently unresolved, in most discs.
△ Less
Submitted 21 July, 2021;
originally announced July 2021.
-
On dust evolution in planet-forming discs in binary systems. II -- Comparison with Taurus and $ρ$ Ophiuchus (sub-)millimetre observations: discs in binaries have small dust sizes
Authors:
Francesco Zagaria,
Giovanni P. Rosotti,
Giuseppe Lodato
Abstract:
The recently discovered exoplanets in binary or higher-order multiple stellar systems sparked a new interest in the study of proto-planetary discs in stellar aggregations. Here we focus on disc solids, as they make up the reservoir out of which exoplanets are assembled and dominate (sub-)millimetre disc observations. These observations suggest that discs in binary systems are fainter and smaller t…
▽ More
The recently discovered exoplanets in binary or higher-order multiple stellar systems sparked a new interest in the study of proto-planetary discs in stellar aggregations. Here we focus on disc solids, as they make up the reservoir out of which exoplanets are assembled and dominate (sub-)millimetre disc observations. These observations suggest that discs in binary systems are fainter and smaller than in isolated systems. In addition, disc dust sizes are consistent with tidal truncation only if they orbit very eccentric binaries. In a previous study we showed that the presence of a stellar companion hastens the radial migration of solids, shortening disc lifetime and challenging planet formation. In this paper we confront our theoretical and numerical results with observations: disc dust fluxes and sizes from our models are computed at ALMA wavelengths and compared with Taurus and $ρ$ Ophiuchus data. A general agreement between theory and observations is found. In particular, we show that the dust disc sizes are generally smaller than the binary truncation radius due to the combined effect of grain growth and radial drift: therefore, small disc sizes do not require implausibly high eccentricities to be explained. Furthermore, the observed binary discs are compatible within $1σ$ with a quadratic flux-radius correlation similar to that found for single-star discs and show a close match with the models. However, the observational sample of resolved binary discs is still small and additional data are required to draw more robust conclusions on the flux-radius correlation and how it depends on the binary properties.
△ Less
Submitted 12 July, 2021;
originally announced July 2021.
-
ALMA 870 $μ$m continuum observations of HD 100546. Evidence of a giant planet on a wide orbit
Authors:
D. Fedele,
C. Toci,
L. Maud,
G. Lodato
Abstract:
This paper reports on a new analysis of archival ALMA $870\,μ$m dust continuum observations. Along with the previously observed bright inner ring ($r \sim 20-40\,$au), two addition substructures are evident in the new continuum image: a wide dust gap, $r \sim 40-150\,$au, and a faint outer ring ranging from $r \sim 150\,$au to $r \sim 250\,$au and whose presence was formerly postulated in low-angu…
▽ More
This paper reports on a new analysis of archival ALMA $870\,μ$m dust continuum observations. Along with the previously observed bright inner ring ($r \sim 20-40\,$au), two addition substructures are evident in the new continuum image: a wide dust gap, $r \sim 40-150\,$au, and a faint outer ring ranging from $r \sim 150\,$au to $r \sim 250\,$au and whose presence was formerly postulated in low-angular-resolution ALMA cycle 0 observations but never before observed. Notably, the dust emission of the outer ring is not homogeneous, and it shows two prominent azimuthal asymmetries that resemble an eccentric ring with eccentricity $e = 0.07 $. The characteristic double-ring dust structure of HD 100546 is likely produced by the interaction of the disk with multiple giant protoplanets. This paper includes new smoothed-particle-hydrodynamic simulations with two giant protoplanets, one inside of the inner dust cavity and one in the dust gap. The simulations qualitatively reproduce the observations, and the final masses and orbital distances of the two planets in the simulations are 3.1 $M_{J}$ at 15 au and 8.5 $M_{J}$ at 110 au, respectively. The massive outer protoplanet substantially perturbs the disk surface density distribution and gas dynamics, producing multiple spiral arms both inward and outward of its orbit. This can explain the observed perturbed gas dynamics inward of 100 au as revealed by ALMA observations of CO.
Finally, the reduced dust surface density in the $\sim 40-150\,$au dust gap can nicely clarify the origin of the previously detected H$_2$O gas and ice emission.
△ Less
Submitted 14 June, 2021;
originally announced June 2021.
-
A dynamical measurement of the disk mass in Elias 2-27
Authors:
Benedetta Veronesi,
Teresa Paneque-Carreno,
Giuseppe Lodato,
Leonardo Testi,
Laura M. Pérez,
Giuseppe Bertin,
Cassandra Hall
Abstract:
Recent multi-wavelength ALMA observations of the protoplanetary disk orbiting around Elias 2-27 revealed a two armed spiral structure. The observed morphology together with the young age of the star and the disk-to-star mass ratio estimated from dust continuum emission make this system a perfect laboratory to investigate the role of self-gravity in the early phases of star formation. This is parti…
▽ More
Recent multi-wavelength ALMA observations of the protoplanetary disk orbiting around Elias 2-27 revealed a two armed spiral structure. The observed morphology together with the young age of the star and the disk-to-star mass ratio estimated from dust continuum emission make this system a perfect laboratory to investigate the role of self-gravity in the early phases of star formation. This is particularly interesting if we consider that gravitational instabilities could be a fundamental first step for the formation of planetesimals and planets. In this Letter, we model the rotation curve obtained by CO data of Elias 2-27 with a theoretical rotation curve including both the disk self-gravity and the star contribution to the gravitational potential. We compare this model with a purely Keplerian one and with a simple power-law function. We find that (especially for the $^{13}$CO isotopologue) the rotation curve is better described by considering not only the star, but also the disk self-gravity. We are thus able to obtain for the first time a dynamical estimate of the disk mass of $0.08\pm0.04\,M_{\odot}$ and the star mass of $0.46\pm0.03\,M_{\odot}$ (in the more general case), the latter being comparable with previous estimates. From these values, we derive that the disk is 17$\%$ of the star mass, meaning that it could be prone to gravitational instabilities. This result would strongly support the hypothesis that the two spiral arms are generated by gravitational instabilities.
△ Less
Submitted 7 May, 2021; v1 submitted 19 April, 2021;
originally announced April 2021.
-
On dust evolution in planet-forming discs in binary systems. I -- Theoretical and numerical modelling: radial drift is faster in binary discs
Authors:
Francesco Zagaria,
Giovanni P. Rosotti,
Giuseppe Lodato
Abstract:
Many stars are in binaries or higher-order multiple stellar systems. Although in recent years a large number of binaries have been proven to host exoplanets, how planet formation proceeds in multiple stellar systems has not been studied much yet from the theoretical standpoint. In this paper we focus on the evolution of the dust grains in planet-forming discs in binaries. We take into account the…
▽ More
Many stars are in binaries or higher-order multiple stellar systems. Although in recent years a large number of binaries have been proven to host exoplanets, how planet formation proceeds in multiple stellar systems has not been studied much yet from the theoretical standpoint. In this paper we focus on the evolution of the dust grains in planet-forming discs in binaries. We take into account the dynamics of gas and dust in discs around each component of a binary system under the hypothesis that the evolution of the circumprimary and the circumsecondary discs is independent. It is known from previous studies that the secular evolution of the gas in binary discs is hastened due to the tidal interactions with their hosting stars. Here we prove that binarity affects dust dynamics too, possibly in a more dramatic way than the gas. In particular, the presence of a stellar companion significantly reduces the amount of solids retained in binary discs because of a faster, more efficient radial drift, ultimately shortening their lifetime. We prove that how rapidly discs disperse depends both on the binary separation, with discs in wider binaries living longer, and on the disc viscosity. Although the less-viscous discs lose high amounts of solids in the earliest stages of their evolution, they are dissipated slowly, while those with higher viscosities show an opposite behaviour. The faster radial migration of dust in binary discs has a striking impact on planet formation, which seems to be inhibited in this hostile environment, unless other disc substructures halt radial drift further in. We conclude that if planetesimal formation were viable in binary discs, this process would take place on very short time scales.
△ Less
Submitted 7 April, 2021;
originally announced April 2021.
-
Spiral Arms and a Massive Dust Disk with non-Keplerian Kinematics: Possible Evidence for Gravitational Instability in the Disk of Elias 2-27
Authors:
T. Paneque-Carreño,
L. M. Perez,
M. Benisty,
C. Hall,
B. Veronesi,
G. Lodato,
A. Sierra,
J. M. Carpenter,
S. M. Andrews,
Jaehan Bae,
Th. Henning,
W. Kwon,
H. Linz,
L. Loinard,
C. Pinte,
L. Ricci,
M. Tazzari,
L. Testi,
D. Wilner
Abstract:
To determine the origin of the spiral structure observed in the dust continuum emission of Elias 2-27 we analyze multi-wavelength continuum ALMA data with a resolution of $\sim$0.2 arcsec ($\sim$23au) at 0.89, 1.3 and 3.3mm. We also study the kinematics of the disk with $^{13}$CO and C$^{18}$O ALMA observations in the $J=$3-2 transition. The spiral arm morphology is recovered at all wavelengths in…
▽ More
To determine the origin of the spiral structure observed in the dust continuum emission of Elias 2-27 we analyze multi-wavelength continuum ALMA data with a resolution of $\sim$0.2 arcsec ($\sim$23au) at 0.89, 1.3 and 3.3mm. We also study the kinematics of the disk with $^{13}$CO and C$^{18}$O ALMA observations in the $J=$3-2 transition. The spiral arm morphology is recovered at all wavelengths in the dust continuum observations, where we measure contrast and spectral index variations along the spiral arms and detect subtle dust-trapping signatures. We determine that the emission from the midplane is cold and interpret the optical depth results as signatures of a higher disk mass than previous constraints. From the gas data, we search for deviations from Keplerian motion and trace the morphology of the emitting surfaces and the velocity profiles. We find an azimuthally varying emission layer height in the system, large-scale emission surrounding the disk, and strong perturbations in the channel maps, co-located with the spirals. Additionally, we develop multigrain dust and gas SPH simulations of a gravitationally unstable disk and compare them to the observations. Given the large scale emission and highly perturbed gas structure, together with the comparison of continuum observations to theoretical predictions, we propose infall-triggered gravitational instabilities as origin for the observed spiral structure.
△ Less
Submitted 24 May, 2021; v1 submitted 25 March, 2021;
originally announced March 2021.