Showing 1–50 of 188 results for author: Lodato, G

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.

Comments: Accepted for publication in A&A letters, 6 pages, 3 images

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.

Comments: 13 pages, 13 figures, accepted for publication on A&A

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.

Comments: 10 pages, 7 figures. Accepted for publication in A&A May 27, 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.

Comments: 12 pages + appendices, 15 figures, accepted to ApJL. Data+movies on zenodo at https://doi.org/10.5281/zenodo.11438154 . Movies on YouTube... watch this one first: https://www.youtube.com/watch?v=KlffYCtYxFU then: https://www.youtube.com/watch?v=bPurbeQNgvI&t=19s and finally watch the expanding "bubble" with rho and T at photosphere: https://www.youtube.com/watch?v=D4dyZeUbfvM

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.

Comments: 14 pages, 12 figures

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.

Comments: 16 pages, 12 figures, accepted for publication on A&A

Journal ref: A&A 686, A9 (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.

Comments: accepted for publication in A&A

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.

Comments: Accepted for publication in ApJ, 12 pages, 7 figures with appendix

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.

Comments: Accepted for publication on Astronomy & Astrophysics. 13 pages, 7 figures + appendix. Abstract abridged to meet arXiv requirements

Journal ref: A&A 679, A82 (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.

Comments: 11 pages, 6 figures, 1 table. Published in ApJL

Journal ref: ApJL 954 (2023) L13

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.

Comments: Accepted for publication in A&A; 13 pages, 13 figures, 1 table

Journal ref: A&A 677, A118 (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.

Comments: Accepted for publication in MNRAS, 20 pages

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.

Comments: 15 pages, 5 figures. Accepted for publication in ApJ

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.

Comments: Accepted by A&A

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.

Comments: 11 pages, 9 figures, accepted for publication on MNRAS

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.

Comments: 17 pages, 14 figures

Journal ref: A&A 671, A82 (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.

Comments: 14 pages, 8 figures, 1 appendix, accepted for publication on MNRAS

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.

Comments: 14 pages, 10 figures, MNRAS in press

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.

Comments: 5 pages, 3 figures

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.

Comments: Accepted for publication on A&A

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.

Comments: 10 figures, 14 pages. Accepted for publication in MNRAS on 25 July 2022

Report number: MN-22-0627-MJ.R1

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.

Comments: Accepted for publication in MNRAS. 14 pages, 8 figures

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.

Comments: Accepted for publication in A&A, 22 pages, 21 figures, 3 tables

Journal ref: A&A 665, A25 (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.

Comments: Accepted for publication in MNRAS, 15 pages, 13 figures, 1 tables, 2 appendices

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.

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

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.

Comments: Review Chapter for Protostars and Planets VII, Editors: Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Final accepted version after interaction with the referees and community feedback. 21 pages (24 with references), 8 figures. Data table available at http://ppvii.org/chapter/15/

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.

Journal ref: Living Reviews in Relativity, Volume 26, Article number: 2 (2023)

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.

Comments: 12 pages, 10 figures, accepted for publication in A&A

Journal ref: A&A 661, A68 (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.

Comments: A&A in press

Journal ref: A&A 663, A98 (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.

Comments: 25 pages, 29 figures, 9 tables. Accepted for publication in A&A

Journal ref: A&A 662, A121 (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.

Comments: Accepted for publication in MNRAS Letters

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.

Comments: Accepted for publication in MNRAS

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.

Comments: 10 Pages, 6 figures. Accepted in ApJL

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.

Comments: Accepted for publication on MNRAS

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.

Comments: Accepted for publication in MNRAS, 15 pages, 13 figures

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.

Comments: Accepted for publication in APJL, 13 pages, 4 figures

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.

Comments: 9 pages, 8 figures, 4 tables

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.

Comments: 19 pages, 14 figures, 1 table, accepted for publication in MNRAS

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.

Comments: 17 pages, 14 figures

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.

Comments: 20 pages, 10 figures; accepted for publication in MNRAS

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.

Comments: Accepted for publication

Journal ref: A&A 651, A90 (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.

Comments: 14 pages, 9 figures, 2 tables, accepted for publication on ApJL

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.

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

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.

Comments: Accepted for publication in ApJ, 38 pages, 23 Figures

PENELLOPE: the ESO data legacy program to complement the Hubble UV Legacy Library of Young Stars (ULLYSES) I. Survey presentation and accretion properties of Orion OB1 and $σ$-Orionis

Authors: C. F. Manara, A. Frasca, L. Venuti, M. Siwak, G. J. Herczeg, N. Calvet, J. Hernandez, Ł. Tychoniec, M. Gangi, J. M. Alcalá, H. M. J. Boffin, B. Nisini, M. Robberto, C. Briceno, J. Campbell-White, A. Sicilia-Aguilar, P. McGinnis, D. Fedele, Á. Kóspál, P. Ábrahám, J. Alonso-Santiago, S. Antoniucci, N. Arulanantham, F. Bacciotti, A. Banzatti , et al. (47 additional authors not shown)

Abstract: The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msu… ▽ More The evolution of young stars and disks is driven by the interplay of several processes, notably accretion and ejection of material. Critical to correctly describe the conditions of planet formation, these processes are best probed spectroscopically. About five-hundred orbits of the Hubble Space Telescope (HST) are being devoted in 2020-2022 to the ULLYSES public survey of about 70 low-mass (M<2Msun) young (age<10 Myr) stars at UV wavelengths. Here we present the PENELLOPE Large Program that is being carried out at the ESO Very Large Telescope (VLT) to acquire, contemporaneous to HST, optical ESPRESSO/UVES high-resolution spectra to investigate the kinematics of the emitting gas, and UV-to-NIR X-Shooter medium-resolution flux-calibrated spectra to provide the fundamental parameters that HST data alone cannot provide, such as extinction and stellar properties. The data obtained by PENELLOPE have no proprietary time, and the fully reduced spectra are made available to the whole community. Here, we describe the data and the first scientific analysis of the accretion properties for the sample of thirteen targets located in the Orion OB1 association and in the sigma-Orionis cluster, observed in Nov-Dec 2020. We find that the accretion rates are in line with those observed previously in similarly young star-forming regions, with a variability on a timescale of days of <3. The comparison of the fits to the continuum excess emission obtained with a slab model on the X-Shooter spectra and the HST/STIS spectra shows a shortcoming in the X-Shooter estimates of <10%, well within the assumed uncertainty. Its origin can be either a wrong UV extinction curve or due to the simplicity of this modelling, and will be investigated in the course of the PENELLOPE program. The combined ULLYSES and PENELLOPE data will be key for a better understanding of the accretion/ejection mechanisms in young stars. △ Less

Submitted 6 April, 2021; v1 submitted 23 March, 2021; originally announced March 2021.

Comments: Accepted for publication on Astronomy & Astrophysics. 15 pages + appendix, language edited version

Journal ref: A&A 650, A196 (2021)

Distinguishing Tidal Disruption Events from Impostors

Authors: Ann Zabludoff, Iair Arcavi, Stephanie La Massa, Hagai B. Perets, Benny Trakhtenbrot, B. Ashley Zauderer, Katie Auchettl, Jane L. Dai, K. Decker French, Tiara Hung, Erin Kara, Giuseppe Lodato, W. Peter Maksym, Yujing Qin, Enrico Ramirez-Ruiz, Nathaniel Roth, Jessie C. Runnoe, Thomas Wevers

Abstract: Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transi… ▽ More Recent claimed detections of tidal disruption events (TDEs) in multi-wavelength data have opened potential new windows into the evolution and properties of otherwise dormant supermassive black holes (SMBHs) in the centres of galaxies. At present, there are several dozen TDE candidates, which share some properties and differ in others. The range in properties is broad enough to overlap other transient types, such as active galactic nuclei (AGN) and supernovae (SNe), which can make TDE classification ambiguous. A further complication is that "TDE signatures" have not been uniformly observed to similar sensitivities or even targeted across all candidates. This chapter reviews those events that are unusual relative to other TDEs, including the possibility of TDEs in pre-existing AGN, and summarises those characteristics thought to best distinguish TDEs from continuously accreting AGN, strongly flaring AGN, SNe, and Gamma-Ray Bursts (GRBs), as well as other potential impostors like stellar collisions, "micro-TDEs," and circumbinary accretion flows. We conclude that multiple observables should be used to classify any one event as a TDE. We also consider the TDE candidate population as a whole, which, for certain host galaxy or SMBH characteristics, is distinguishable statistically from non-TDEs, suggesting that at least some TDE candidates do in fact arise from SMBH-disrupted stars. △ Less

Submitted 8 April, 2021; v1 submitted 22 March, 2021; originally announced March 2021.

Comments: 57 pages, 17 figures. Accepted for publication in Springer Space Science Reviews. Chapter in ISSI review "The Tidal Disruption of Stars by Massive Black Holes" vol. 79. Corrected several typos from previous submission

Dynamical dust traps in misaligned circumbinary discs: analytical theory and numerical simulations

Authors: Cristiano Longarini, Giuseppe Lodato, Claudia Toci, Hossam Aly

Abstract: Recent observations have shown that circumbinary discs can be misaligned with respect to the binary orbital plane.The lack of spherical symmetry, together with the non-planar geometry of these systems, causes differential precession which might induce the propagation of warps. While gas dynamics in such environments is well understood, little is known about dusty discs. In this work, we analytical… ▽ More Recent observations have shown that circumbinary discs can be misaligned with respect to the binary orbital plane.The lack of spherical symmetry, together with the non-planar geometry of these systems, causes differential precession which might induce the propagation of warps. While gas dynamics in such environments is well understood, little is known about dusty discs. In this work, we analytically study the problem of dust traps formation in misaligned circumbinary discs. We find that pile-ups may be induced not by pressure maxima, as the usual dust traps, but by a difference in precession rates between the gas and dust. Indeed, this difference makes the radial drift inefficient in two locations, leading to the formation of two dust rings whose position depends on the system parameters. This phenomenon is likely to occur to marginally coupled dust particles $(\text{St}\gtrsim1)$ as both the effect of gravitational and drag force are considerable. We then perform a suite of three-dimensional SPH numerical simulations to compare the results with our theoretical predictions. We explore the parameter space, varying stellar mass ratio, disc thickness, radial extension, and we find a general agreement with the analytical expectations. Such dust pile-up prevents radial drift, fosters dust growth and may thus promote the planet formation in circumbinary discs. △ Less

Submitted 22 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS, 12 pages, 12 figures, 2 tables, 2 appendices

Observable gravitational waves from tidal disruption events and their electromagnetic counterpart

Authors: Hugo Pfister, Martina Toscani, Thomas Hong Tsun Wong, Jane Lixin Dai, Giuseppe Lodato, Elena M. Rossi

Abstract: We estimate the rate of tidal disruption events (TDEs) that will be detectable with future gravitational wave detectors as well as the most probable properties of these events and their possible electromagnetic counterpart. To this purpose we combine standard gravitational waves and electromagnetic results with detailed rates estimates. We find that the \emph{Laser Interferometer Space Antenna} (L… ▽ More We estimate the rate of tidal disruption events (TDEs) that will be detectable with future gravitational wave detectors as well as the most probable properties of these events and their possible electromagnetic counterpart. To this purpose we combine standard gravitational waves and electromagnetic results with detailed rates estimates. We find that the \emph{Laser Interferometer Space Antenna} (LISA) should not detect any TDEs, unless black holes (BHs) are typically embedded by a young stellar population which, in this situation, could lead up to few 10 events during the duration of the mission. If there are gravitational wave observations, these events should also be observable in the X-ray or the optical/UV part of the electromagnetic spectrum, which may open up the multi-messenger era for TDEs. The generation of detectors following LISA will at least yearly observe $10^4$ TDEs at cosmological distances, allowing to do population studies and constrain the black hole mass function. In all cases, most probable events should be around black holes with a mass such that the Keplerian frequency at the Schwarzschild radius is similar to the optimal frequency of the detector and with a large penetration factor. △ Less

Submitted 17 November, 2021; v1 submitted 10 March, 2021; originally announced March 2021.

Comments: accepted to MNRAS

The theory of kinks -- I. A semi-analytic model of velocity perturbations due to planet-disc interaction

Authors: Francesco Bollati, Giuseppe Lodato, Daniel J. Price, Christophe Pinte

Abstract: A new technique to detect protoplanets is by observing the kinematics of the surrounding gas. Gravitational perturbations from a planet produce peculiar `kinks' in channel maps of different gas species. In this paper, we show that such kinks can be reproduced using semi-analytic models for the velocity perturbation induced by a planet. In doing so we i) confirm that the observed kinks are consiste… ▽ More A new technique to detect protoplanets is by observing the kinematics of the surrounding gas. Gravitational perturbations from a planet produce peculiar `kinks' in channel maps of different gas species. In this paper, we show that such kinks can be reproduced using semi-analytic models for the velocity perturbation induced by a planet. In doing so we i) confirm that the observed kinks are consistent with the planet-induced wake; ii) show how to quantify the planet mass from the kink amplitude; in particular, we show that the kink amplitude scales with the square root of the planet mass for channels far from the planet velocity, steepening to linear as the channels approach the planet; iii) show how to extend the theory to include the effect of damping, which may be needed in order to have localized kinks. △ Less

Submitted 19 April, 2021; v1 submitted 26 February, 2021; originally announced March 2021.

Comments: Submission to MNRAS Accepted