Search

WASP-49Ab, a low-density, Saturn-like planet in a tight orbit around a Sun-like star within a wide binary system, is a compelling candidate for hosting a volcanic moon, as suggested by the detection of Doppler-shifted sodium.This study evaluates the stability of potential satellites around WASP-49Ab under the influence of planetary oblateness, relativistic effects, and perturbations from a close companion star, focusing on their impact on light curve parameters such as transit duration and impact parameter variations, driven by the evolution of the planet's orbit in this extreme environment. Using N-body simulations and semi-analytical methods, we analysed moon's dynamics across varied initial conditions and gravitational frameworks including the potential of an oblate planet and the effects of the general relativity. We find that `selenity', a moon survival indicator, is high in close orbits with low eccentricity, near the Roche limit, especially for masses greater than Io's. Stability decreases as eccentricity or distance from the planet increases. Additionally, We observe a strong destabilising resonance near 1.4 planetary radii when planetary eccentricities are slightly greater than zero. This study confirms the potential for stable exomoons around WASP-49Ab despite its hostile environment, emphasizing the importance of incorporating diverse physical effects in stability analyses, aiding future detection efforts., Comment: Accepted for publication in Astronomy & Astrophysics Letters

Stars and planets form in regions of enhanced stellar density, subjecting protoplanetary discs to gravitational perturbations from neighbouring stars. Observations in the Taurus star-forming have uncovered evidence of at least three recent, star-disc encounters that have truncated discs (HV/DO Tau, RW Aurigae, UX Tau), raising questions about the frequency of such events. We aim to assess the probability of observing truncating star-disc encounters in Taurus. We generate a physically motivated dynamical model including binaries and spatial-kinematic substructure to follow the historical dynamical evolution and stellar encounters in the Taurus star forming region. We track the star-disc encounters and outer disc radius evolution over the lifetime of Taurus. A quarter of discs are truncated below 30 au by dynamical encounters, but this truncation mostly occurs in binaries over the course of a few orbital periods, on a time-scale $\lesssim 0.1$ Myr. Nonetheless, some truncating encounters still occur up to the present age of Taurus. Strongly truncating encounters (ejecting $\gtrsim 10$ percent of the disc mass) occur at a rate $\sim 10$ Myr$^{-1}$, sufficient to explain the encounter between HV and DO Tau $\sim 0.1$ Myr ago. If encounters that eject only $\sim 1$ percent of the disc mass are responsible for RW Aurigae and UX Tau, then they are also expected with encounter rate $\Gamma_\mathrm{enc} \sim 100{-}200$ Myr$^{-1}$. However, the observed sample of recent encounters is probably incomplete, since these examples occurred in systems that are not consistent with random drawing from the mass function. One more observed example would statistically imply additional physics, such as replenishment of the outer disc material., Comment: 20 pages, 20 figures, accepted for publication in A&A

Molecular deuteration is a powerful diagnostic tool for probing the physical conditions and chemical processes in astrophysical environments. In this work, we focus on formaldehyde deuteration in the protobinary system NGC\,1333 IRAS\,4A, located in the Perseus molecular cloud. Using high-resolution ($\sim$\,100\,au) ALMA observations, we investigate the [D$_2$CO]/[HDCO] ratio along the cavity walls of the outflows emanating from IRAS\,4A1. Our analysis reveals a consistent decrease in the deuteration ratio (from $\sim$\,60-20\% to $\sim$\,10\%) with increasing distance from the protostar (from $\sim$\,2000\,au to $\sim$\,4000\,au). Given the large measured [D$_2$CO]/[HDCO], both HDCO and D$_2$CO are likely injected by the shocks along the cavity walls into the gas-phase from the dust mantles, formed in the previous prestellar phase. We propose that the observed [D$_2$CO]/[HDCO] decrease is due to the density profile of the prestellar core from which NGC\,1333 IRAS\,4A was born. When considering the chemical processes at the base of formaldehyde deuteration, the IRAS\,4A's prestellar precursor had a predominantly flat density profile within 3000\,au and a decrease of density beyond this radius.

The accretion/ejection processes in T-Tauri stars are fundamental to their physical evolution, while also impacting the properties and evolution of the circumstellar material at a time when planet formation takes place. To this date, characterization of ongoing accretion processes in stellar pairs at 5-50\,au scales has been challenging, high angular resolution spectrographs are required to extract the spectral features of each component. We present the analysis of spectroscopic observations of the tight (160mas, 25au) T-Tauri system HT Lup A/B, obtained with MUSE at VLT in March and July of 2021. We focus on constraining the accretion/ejection processes and variability of the secondary component HT Lup B, by searching for accretion tracers applying High-Resolution Spectral Differential Imaging techniques. We retrieve strong (SNR $>$ 5) $H\alpha, H\beta$ and [OI]$\lambda6300$ emission in both epochs. The $H\alpha$ and $H\beta$ line fluxes showcase high variability, with variations up to 400-500\% between epochs. The fluxes are consistent with accretion rates of $8\times10^{-9} M_\odot \, yr^{-1}$ and $2\times10^{-9} M_\odot \, yr^{-1}$ for the first and second epoch, respectively. We attribute the increased accretion activity during the first night to a "burst" like event, followed by a relaxation period more representative of the common accretion activity of the system. The [OI]$\lambda6300$ line profiles remain relatively similar between epochs and suggest ejection rates on the order of $10^{-9}-10^{-10} M_\odot \, yr^{-1}$, compatible with moderate disk winds emission. Our results also indicate that the accretion processes of HT Lup B are compatible with Classical T Tauri Stars, unlike previous classifications, Comment: 28 pages, 13 fgures, Accepted by ApJ

Rings are complex structures surrounding giant planets and some minor bodies in the Solar System. While some formation mechanisms could also potentially foster their existence around (regular or irregular) satellites, none of these bodies currently bear these structures. We aim to understand the underlying mechanisms that govern the potential formation, stability, and/or decay of hypothetical circumsatellital rings (CSRs), orbiting the largest moons in the Solar System. This extends to the exploration of short-term morphological features within these rings, providing insights into the ring survival time-scales and the interactions that drive their evolution. To conduct this study, we use numerical N-body simulations under the perturbing influence of the host planet and other moon companions. We found that moons with a lower Roche-to-Hill radius can preserve their rings over extended periods. Moreover, the gravitational environment in which these rings are immersed influences the system's morphological evolution, inducing gaps through the excitation of eccentricity and inclination of constituent particles. Specifically, our results show that Iapetus' and Rhea's rings experience minimal variations in their orbital parameters, enhancing their long-term stability. This agrees with the hypothesis that some of the features of Iapetus and Rhea were produced by ancient ring systems, for example, the huge ridge in Iapetus equator as a result of a decaying ring. From a dynamical perspective, we found that there are no mechanisms that preclude the existence of CSRs and we attribute their current absence to non-gravitational phenomena. Effects such as stellar radiation, magnetic fields, and the influence of magnetospheric plasma can significantly impact the dynamics of constituent particles and trigger their decay, highlighing the importance of future studies on these effects., Comment: Accepted for publication in Astronomy & Astrophysics

The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.

V892 Tau is a young binary star surrounded by a circumbinary disc which show hints of interaction with the low-mass nearby star V892 Tau NE. The goal of this paper is to constrain the orbit of V892 Tau NE and to determine the resulting circumbinary disc dynamics. We present new ALMA observations of the V892 Tau circumbinary disc at a twice higher angular and spectral resolution. We model the data with V892 Tau as a triple system and perform a grid of hydrodynamical simulations testing several orbits of the companion. The simulation outputs are then post-processed to build synthetic maps that we compare to the observations. The 12CO emission of the disc shows clear non-Keplerian features such as spiral arms. When comparing the data with our synthetic observations, we interpret these features as ongoing interactions with the companion. Our simulations indicate that an eccentricity of 0.5 of the companion is needed to reproduce the observed disc extent and that a mutual inclination of approximately 60{\deg} with the inner binary reproduces the measured disc tilt. In order to explain most of the features of the circumbinary disc, we propose that V892 Tau NE follows a misaligned eccentric orbit, with an eccentricity between 0.2 and 0.5 and a mutual inclination between 30{\deg} and 60{\deg}. Such a misaligned companion suggests the disc is oscillating and precessing with time, stabilising in an intermediate plane with a non-zero mutual inclination with the inner binary. Given that orbital configuration, we show that the stability of future planets is compromised in the second half of the disc once the gas has dissipated., Comment: Published in A&A, 29 pages, 15 figures

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 $\mu$m ALMA dust continuum emission and 1.67 $\mu$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, $\tau_{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 ($\Delta\theta= 30^\circ$) on timescales shorter than $\tau_{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 $\tau_{ecc}$ and the models do not reproduce the observed morphology anymore. This implies that either the age of the system is shorter than $\tau_{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., Comment: 14 pages, 12 figures

High-order multiple (triple and beyond) systems are relatively common. Their interaction with circumstellar and circumbinary material can have a large impact on the formation and evolution of planetary systems and depends on their orbital properties. GG\,Tau and UX\,Tau are two pre-main sequence high-order multiple systems in which the tightest pair has a projected separation of $\approx5$--20\,au. Characterizing precisely their orbits is crucial to establish their long-term stability, to predict the dynamics and evolution of circumstellar matter, and to evaluate the potential for planet formation in such systems. We combine existing astrometric measurements with previously unpublished high-resolution observations of the GG\,Tau\,Ab and UX\,Tau\,B pairs and perform Keplerian orbital fits. For GG\,Tau\,Ab the data presented here represent the first detection of orbital motion. For both systems they yield dramatic increases in orbital coverage ($\gtrsim60\%$ and $\approx100\%$ for UX\,Tau\,B and GG\,Tau\,Ab, for orbital periods of $\approx32$ and $\approx8$\,yr, respectively) and allow us to obtain well-constrained orbital fits, including dynamical masses with $\lesssim10\%$ and $\lesssim7\%$ random and systematic uncertainties. We find that both GG\,Tau\,A and UX\,Tau\,A--B likely form stable hierarchical systems, although one possible deprojection solution for GG\,Tau is strongly misaligned and could experience von Zeipel-Lidov-Kozai oscillations. We further find that the UX\,Tau\,B orbit is much more eccentric than the GG\,Tau\,Ab one, possibly explaining the lack of circumstellar material in the former. The newly-determined orbits revive the question of the dynamical fate of gas and dust in these two hierarchical systems and should spur new dedicated simulations to assess the long-term evolution of the systems and the dynamical perturbations imposed by the close binaries they host., Comment: Accepted for publication in Astronomy & Astrophysics

Alpha Centauri is a triple stellar system, and it contains the closest star to Earth (Proxima Centauri). Over the last decades, the stars in Alpha Cen and their orbits have been investigated in great detail. However, the possible scenarios for planet formation and evolution in this triple stellar system remain to be explored further. First, we present a 3D hydrodynamical simulation of the circumstellar discs in the binary Alpha Cen AB. Then, we compute stability maps for the planets within Alpha Cen obtained through N-body integrations. Last, we estimate the radial velocity (RV) signals of such planets. We find that the circumstellar discs within the binary cannot exceed 3 au in radius and that the available dust mass to form planets is about 30 $M_\oplus$. Planets around A and B are stable if their semimajor axes are below 3 au, while those around C are stable and remain unperturbed by the binary AB. For rocky planets, the planetary mass has only a mild effect on the stability. Therefore, Alpha Cen could have formed and hosted rocky planets around each star, which may be detected with RV methods in the future. The exoplanetary hunt in this triple stellar system must continue., Comment: 17 pages, 12 figures, accepted for publication in Universe

Some post-common-envelope binaries are binary stars with short periods that exhibit significant period variations over long observational time spans. These eclipse timing variations (ETVs) are most likely to be accounted for by the presence of an unseen massive companion, potentially of planetary or substellar nature, and the light-travel time (LTT) effect. In this study, our main objective is to describe the diversity of compatible nontransit companions around PCEBs and explore the robustness of the solutions by employing tools for uncertainty estimation. We select the controversial data of the QS Vir binary star, which previous studies have suggested hosts a planet. We employ a minimizing strategy, using genetic algorithms to explore the global parameter space followed by refinement of the solution using the simplex method. We evaluate errors through the classical MCMC approach and discuss the error range for parameters. Our results highlight the strong dependence of ETV models for close binaries on the dataset used, which leads to relatively loose constraints on the parameters of the unseen companion. We find that the shape of the $O-C$ curve is influenced by the dataset employed. We propose an alternative method to evaluate errors on the orbital fits based on a grid search surrounding the best-fit values, obtaining a wider range of plausible solutions that are compatible with goodness-of-fit statistics. We also analyze how the parameter solutions are affected by the choice of the dataset, and find that this system continuously changes the compatible solutions as new data are obtained from eclipses. The best-fit parameters for QS Vir correspond to a low-mass stellar companion (57.71 $M_{jup}$ ranging from 40 to 64 $M_{jup}$) on an eccentric orbit ($e=0.91^{+0.07}_{-0.17}$) with a variety of potential periods ($P = 16.69 ^{+0.47}_{-0.42}$ yr.), Comment: 13 Pages. 14 Figures. A&A A&A approved. article aa47030-23

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., Comment: accepted for publication in A&A

Observations reveal protoplanetary discs being perturbed by flyby candidates. We simulate a scenario where an unbound perturber, i.e., a flyby, undergoes an inclined grazing encounter, capturing material and forming a second-generation protoplanetary disc. We run $N$--body and three-dimensional hydrodynamical simulations of a parabolic flyby grazing a particle disc and a gas-rich protoplanetary disc, respectively. In both our $N$--body and hydrodynamic simulations, we find that the captured, second-generation disc forms at a tilt twice the initial flyby tilt. This relationship is robust to variations in the flyby's tilt, position angle, periastron, and mass. We extend this concept by also simulating the case where the flyby has a disc of material prior to the encounter but we do not find the same trend. An inclined disc with respect to the primary disc around a misaligned flyby is tilted by a few degrees, remaining close to its initial disc tilt. Therefore, if a disc is present around the flyby before the encounter, the disc may not tilt up to twice the perturber tilt depending on the balance between the angular momentum of the circumsecondary disc and captured particles. In the case where the perturber has no initial disc, analyzing the orientation of these second-generation discs can give information about the orbital properties of the flyby encounter., Comment: 17 pages, 15 figures, accepted to MNRAS

Episodic accretion is a low-mass pre-main sequence phenomenon characterized by sudden outbursts of enhanced accretion. These objects are classified into two: protostars with elevated levels of accretion that lasts for decades or more, called FUors, and protostars with shorter and repetitive bursts, called EXors. HBC 494 is a FUor object embedded in the Orion Molecular Cloud. Earlier Atacama Large (sub-)Millimeter Array (ALMA) continuum observations showed an asymmetry in the disk at 0.''2 resolution. Here, we present follow-up observations at ~0.''03, resolving the system into two components: HBC 494 N (primary) and HBC 494 S (secondary). No circumbinary disk was detected. Both disks are resolved with a projected separation of ~0.''18 (75 au). Their projected dimensions are 84+/-1.8 x 66.9+/-1.5 mas for HBC 494 N and 64.6+/-2.5 x 46.0+/-1.9 mas for HBC 494 S. The disks are almost aligned and with similar inclinations. The observations show that the primary is ~5 times brighter/more massive and ~2 times bigger than the secondary. We notice that the northern component has a similar mass to the FUors, while the southern has to EXors. The HBC 494 disks show individual sizes that are smaller than single eruptive YSOs. In this work, we also report 12CO, 13CO, and C18O molecular line observations. At large scale, the 12CO emission shows bipolar outflows, while the 13CO and C18O maps show a rotating and infalling envelope around the system. At a smaller scale, the 12CO and 13CO moment zero maps show cavities within the continuum disks' area, which may indicate continuum over-subtraction or slow-moving jets and chemical destruction along the line-of-sight., Comment: Published in MNRAS. Main text (11 pages, 8 figures). Appendix (11 pages, 9 figures). Total (22 pages, 17 figures)

Spiral arms are observed in numerous protoplanetary discs. These spiral arms can be excited by companions, either on bound or unbound orbits. We simulate a scenario where an unbound perturber, i.e. a flyby, excites spiral arms during a periastron passage. We run three-dimensional hydrodynamical simulations of a parabolic flyby encountering a gaseous protoplanetary disc. The perturber mass ranges from $10\, \rm M_J$ to $1\, \rm M_{\odot}$. The perturber excites a two-armed spiral structure, with a more prominent spiral feature for higher mass perturbers. The two arms evolve over time, eventually winding up, consistent with previous works. We focus on analysing the pattern speed and pitch angle of these spirals during the whole process. The initial pattern speed of the two arms are close to the angular velocity of the perturber at periastron, and then it decreases over time. The pitch angle also decreases over time as the spiral winds up. The spirals disappear after several local orbital times. An inclined prograde orbit flyby induces similar disc substructures as a coplanar flyby. A solar-mass flyby event causes increased eccentricity growth in the protoplanetary disc, leading to an eccentric disc structure which dampens over time. The spirals' morphology and the disc eccentricity can be used to search for potential unbound stars or planets around discs where a flyby is suspected. Future disc observations at high resolution and dedicated surveys will help to constrain the frequency of such stellar encounters in nearby star-forming regions., Comment: 17 pages, 18 figures, accepted to MNRAS

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., Comment: 11 pages, 9 figures, accepted for publication on MNRAS

Dense stellar environments as hosts of ongoing star formation increase the probability of gravitational encounters among stellar systems during the early stages of evolution. Stellar interaction may occur through non-recurring, hyperbolic or parabolic passages (a so-called 'fly-by'), through secular binary evolution, or through binary capture. In all three scenarios, the strong gravitational perturbation is expected to manifest itself in the disc structures around the individual stars. Here, we present near-infrared polarised light observations that were taken with the SPHERE/IRDIS instrument of three known interacting twin-disc systems: AS 205, EM* SR 24, and FU Orionis. The scattered light exposes spirals likely caused by the gravitational interaction. On a larger scale, we observe connecting filaments between the stars. We analyse their very complex polarised intensity and put particular attention to the presence of multiple light sources in these systems. The local angle of linear polarisation indicates the source whose light dominates the scattering process from the bridging region between the two stars. Further, we show that the polarised intensity from scattering with multiple relevant light sources results from an incoherent summation of the individuals' contribution. This can produce nulls of polarised intensity in an image, as potentially observed in AS 205. We discuss the geometry and content of the systems by comparing the polarised light observations with other data at similar resolution, namely with ALMA continuum and gas emission. Collective observational data can constrain the systems' geometry and stellar trajectories, with the important potential to differentiate between dynamical scenarios of stellar interaction., Comment: 23 pages, 17 figures, accepted for publication in MNRAS

Most of the planetary systems discovered around binary stars are located at approximately three semi-major axes from the barycentre of their system, curiously close to low-order mean-motion resonances (MMRs). The formation mechanism of these circumbinary planets is not yet fully understood. In situ formation is extremely challenging because of the strong interaction with the binary. One possible explanation is that, after their formation, the interactions between these planets and the surrounding protoplanetary disc cause them to migrate at velocities dependent on the nature of the disc and the mass of the exoplanet. Although extensive data can be obtained with direct hydrodynamical simulations, their computational cost remains too high. On the other hand, the direct n-body simulations approach allows us to model a large variety of parameters at much lower cost. We analyse the planetary migration around a wide variety of binary stars using Stokes-like forces that mimic planetary migration at a constant rate. Our goal is to identify the main parameters responsible for the ejection of planets at different resonances with the inner binary. We performed 4200 n-body simulations with Stokes-like forces and analysed their evolution and outcome as a function of the properties of each system. For each simulated exoplanet, we applied an ensemble learning method for classification in order to clarify the relationship between the inspected parameters and the process of MMR capture. We identify the capture probability for different N/1 MMRs, 4/1 being the most prone to capture exoplanets, with 37% probability, followed by MMR 5/1 with $\sim$ 23% of probability. The eccentricity of the binary is found to be the most important parameter in determining the MMR capture of each circumbinary exoplanet, followed by the mass ratio of the binary and the initial eccentricity of the planet., Comment: 13 pages, 10 figures, A\&A accepted

We perform 3D hydrodynamics simulations of disc-disc stellar flybys with on-the-fly Monte Carlo radiative transfer. We show that pre-existing circumstellar discs around both stars result in fast rising ($\sim$yrs) outbursts lasting 2-5 times longer than for a star-disc flyby. The perturber always goes into outburst ($\dot{M}>10^{-5}~{\rm M_{\odot}~ yr^{-1}}$). Whereas we find that the primary goes into a decades long outburst only when the flyby is retrograde to the circumprimary disc rotation. High accretion rates during the outburst are triggered by angular momentum cancellation in misaligned material generated by the encounter. A large fraction of accreted material is alien., Comment: 11 pages, 9 figures, accepted for publication in MNRAS

We review the role of stellar flybys and encounters in shaping planet-forming discs around young stars, based on the published literature on this topic in the last 30 years. Since most stars $\leq~2$ Myr old harbour protoplanetary discs, tidal perturbations affect planet formation. First, we examine the probability of experiencing flybys or encounters: More than 50\% of stars with planet-forming discs in a typical star forming environment should experience a close stellar encounter or flyby within 1000 au. Second, we detail the dynamical effects of flybys on planet-forming discs. Prograde, parabolic, disc-penetrating flybys are the most destructive. Grazing and penetrating flybys in particular lead to the capture of disc material by the secondary to form a highly misaligned circumsecondary disc with respect to the disc around the primary. One or both discs may undergo extreme accretion and outburst events, similar to the ones observed in FU Orionis-type stars. Warps and broken discs are distinct signatures of retrograde flybys. Third, we review some recently observed stellar systems with discs where a stellar flyby or an encounter is suspected -- including UX Tau, RW Aur, AS 205, Z CMa, and FU Ori. Finally, we discuss the implications of stellar flybys for planet formation and exoplanet demographics, including possible imprints of a flyby in the Solar System in the orbits of trans-Neptunian objects and the Sun's obliquity., Comment: 16 pages, 7 figures, 1 table, 255 references. Invited review article accepted in EPJ+

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., Comment: Accepted for publication in MNRAS, 15 pages, 13 figures, 1 tables, 2 appendices

We report a study of the low-mass Class-0 multiple system VLA 1623AB in the Ophiuchus star-forming region, using H$^{13}$CO$^+$ ($J=3-2$), CS ($J=5-4$), and CCH ($N=3-2$) lines as part of the ALMA Large Program FAUST. The analysis of the velocity fields revealed the rotation motion in the envelope and the velocity gradients in the outflows (about 2000 au down to 50 au). We further investigated the rotation of the circum-binary VLA 1623A disk as well as the VLA 1623B disk. We found that the minor axis of the circum-binary disk of VLA 1623A is misaligned by about 12 degrees with respect to the large-scale outflow and the rotation axis of the envelope. In contrast, the minor axis of the circum-binary disk is parallel to the large-scale magnetic field according to previous dust polarization observations, suggesting that the misalignment may be caused by the different directions of the envelope rotation and the magnetic field. If the velocity gradient of the outflow is caused by rotation, the outflow has a constant angular momentum and the launching radius is estimated to be $5-16$ au, although it cannot be ruled out that the velocity gradient is driven by entrainments of the two high-velocity outflows. Furthermore, we detected for the first time a velocity gradient associated with rotation toward the VLA 16293B disk. The velocity gradient is opposite to the one from the large-scale envelope, outflow, and circum-binary disk. The origin of its opposite gradient is also discussed., Comment: 27 pages, 21 figures, 2 Tables, Accepted for publication in ApJ

Close encounters between young stellar objects in star forming clusters are expected to dramatically perturb circumstellar disks. Such events are witnessed in numerical simulations of star formation, but few direct observations of ongoing encounters have been made. Here we report sub-0".1 resolution Atacama Large Millimeter Array (ALMA) and Jansky Very Large Array (JVLA) observations towards the million year old binary protostar Z CMa in dust continuum and molecular line emission. A point source ~4700 au from the binary has been discovered at both millimeter and centimeter wavelengths. It is located along the extension of a ~2000 au streamer structure previously found in scattered light imaging, whose counterpart in dust and gas emission is also newly identified. Comparison with simulations shows signposts of a rare flyby event in action. Z CMa is a "double burster", as both binary components undergo accretion outbursts, which may be facilitated by perturbations to the host disk by flybys., Comment: Published in Nature Astronomy. Here is the authors' version with the Supplementary Information integrated into the Methods section

In recent years, technical and theoretical work to detect moons and rings around exoplanets has been attempted. The small mass/size ratios between moons and planets means this is very challenging, having only one exoplanetary system where spotting an exomoon might be feasible (i.e. Kepler-1625b i). In this work, we study the dynamical evolution of ringed exomoons, dubbed "cronomoons" after their similarity with Cronus (Greek for Saturn), and after Chronos (the epitome of time), following the Transit Timing Variations (TTV) and Transit Duration Variation (TDV) that they produce on their host planet. Cronomoons have extended systems of rings that make them appear bigger than they actually are when transiting in front of their host star. We explore different possible scenarios that could lead to the formation of such circumsatellital rings, and through the study of the dynamical/thermodynamic stability and lifespan of their dust and ice ring particles, we found that an isolated cronomoon can survive for time-scales long enough to be detected and followed up. If these objects exist, cronomoons' rings will exhibit gaps similar to Saturn's Cassini Division and analogous to the asteroid belt's Kirkwood gaps, but instead raised due to resonances induced by the host planet. Finally, we analyse the case of Kepler-1625b i under the scope of this work, finding that the controversial giant moon could instead be an Earth-mass cronomoon. From a theoretical perspective, this scenario can contribute to a better interpretation of the underlying phenomenology in current and future observations., Comment: Accepted to MNRAS. See a Twitter thread at @MarioSucerquia. Press release: http://www.npf.cl/ 14 pages, 11 figures

We examine whether stellar flyby simulations can initiate FU Orionis outbursts using 3D hydrodynamics simulations coupled to live Monte Carlo radiative transfer. We find that a flyby where the secondary penetrates the circumprimary disc triggers a 1-2 year rise in the mass accretion rate to $10^{-4}~{\rm M_{\odot}~ yr^{-1}}$ that remains high ($\gtrsim 10^{-5}~{\rm M_{\odot}~yr^{-1}}$) for more than a hundred years, similar to the outburst observed in FU Ori. Importantly, we find that the less massive star becomes the dominant accretor, as observed., Comment: 5 pages, 4 figures, accepted for publication in MNRAS, movies available at https://emborchert.github.io/videos, press release available at https://www.monash.edu/science/news/current/astronomers-find-clue-to-solar-system-formation-through-well-known-star

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$\alpha$ 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 $\theta\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., Comment: 19 pages, 14 figures, 1 table, accepted for publication in MNRAS

The gas dissipation from a protoplanetary disk is one of the key processes affecting planet formation, and it is widely accepted that it happens on timescales of a few million years for disks around single stars. Over the last years, several protoplanetary disks have been discovered in multiple star systems, and despite the complex environment in which they find themselves, some of them seem to be quite old, a situation that may favor planet formation. A clear example of this is the disk around HD 98800 B, a binary in a hierarchical quadruple stellar system, which at a $\sim$10 Myr age seems to still be holding significant amounts of gas. Here we present a 1D+1D model to compute the vertical structure and gas evolution of circumbinary disks in hierarchical triple star systems considering different stellar and disk parameters. We show that tidal torques due to the inner binary together with the truncation of the disk due to the external companion strongly reduce the viscous accretion and expansion of the disk. Even allowing viscous accretion by tidal streams, disks in these kind of environments can survive for more than 10 Myr, depending on their properties, with photoevaporation being the main gas dissipation mechanism. We particularly apply our model to the circumbinary disk around HD 98800 B and confirm that its longevity, along with the current non-existence of a disk around the companion binary HD 98800 A, can be explained with our model and by this mechanism., Comment: 20 pages, 9 figures. Accepted for publication in ApJ

We present sensitive ALMA observations of TWA 3, a nearby, young ($\sim$10 Myr) hierarchical system composed of three pre-main sequence M3--M4.5 stars. For the first time, we detected ${}^{12}$CO and ${}^{13}$CO $J$=2-1 emission from the circumbinary protoplanetary disk around TWA 3A. We jointly fit the protoplanetary disk velocity field, stellar astrometric positions, and stellar radial velocities to infer the architecture of the system. The Aa and Ab stars ($0.29\pm0.01\,M_\odot$ and $0.24\pm0.01\,M_\odot$, respectively) comprising the tight ($P=35$ days) eccentric ($e=0.63\pm0.01$) spectroscopic binary are coplanar with their circumbinary disk (misalignment $< 6^{\circ}$ with 68% confidence), similar to other short-period binary systems. From models of the spectral energy distribution, we found the inner radius of the circumbinary disk ($r_\mathrm{inner} = 0.50 - 0.75$ au) to be consistent with theoretical predictions of dynamical truncation $r_\mathrm{cav}/a_\mathrm{inner} \approx 3$. The outer orbit of the tertiary star B ($0.40\pm0.28\,M_\odot$, $a\sim65 \pm 18$ au, $e=0.3\pm0.2$) is not as well constrained as the inner orbit, however, orbits coplanar with the A system are still preferred (misalignment $ < 20^{\circ}$). To better understand the influence of the B orbit on the TWA 3A circumbinary disk, we performed SPH simulations of the system and found that the outer edge of the gas disk ($r_\mathrm{outer}=8.5\pm0.2$ au) is most consistent with truncation from a coplanar, circular or moderately eccentric orbit, supporting the preference from the joint orbital fit., Comment: 19 pages, 9 figures, accepted to ApJ

Planets form in young circumstellar disks called protoplanetary disks. However, it is still difficult to catch planet formation in-situ. Nevertheless, from recent ALMA/SPHERE data, encouraging evidence of the direct and indirect presence of embedded planets has been identified in disks around young stars: co-moving point sources, gravitational perturbations, rings, cavities, and emission dips or shadows cast on disks. The interpretation of these observations needs a robust physical framework to deduce the complex disk geometry. In particular, protoplanetary disk models usually assume the gas pressure scale-height given by the ratio of the sound speed over the azimuthal velocity $H/r = c_{s\rm }/v_{\rm k}$. By doing so, \textit{radiative} pressure fields are often ignored, which could lead to a misinterpretation of the real vertical structure of such disks. We follow the evolution of a gaseous disk with an embedded Jupiter mass planet through hydrodynamical simulations, computing the disk scale-height including radiative pressure, which was derived from a generalization of the stellar atmosphere theory. We focus on the vertical impact of the radiative pressure in the vicinity of circumplanetary disks, where temperatures can reach $\gtrsim 1000$ K for an accreting planet, and radiative forces can overcome gravitational forces from the planet. The radiation-pressure effects create a vertical optically thick column of gas and dust at the proto-planet location, casting a shadow in scattered light. This mechanism could explain the peculiar illumination patterns observed in some disks around young stars such as HD 169142 where a moving shadow has been detected, or the extremely high aspect-ratio $H/r \sim 0.2$ observed in systems like AB Aur and CT Cha., Comment: 12 pages, 7 figures. Accepted for publication by ApJ

The protoplanetary disc HD 100453 exhibits a curious combination of spirals, shadows and a relative misalignment between the observed outer disc and inferred inner disc. This disc is accompanied by a secondary star on a bound orbit exterior to the disc. Recent observations have suggested there may be an additional low-mass companion residing within the disc inner cavity. In our companion paper the orbit of the secondary was shown to be misaligned by 61 degrees to the plane of the outer disc. Here we investigate the properties of the inner companion and the origin of the misalignment between the inner and outer disc. Using numerical simulations and synthetic observations, we show that the disc structure and kinematics are consistent with a $\lesssim$ 5 Jupiter mass planet located at 15-20au. We find that the disc evolution over around 50 binary orbits (about 10$^5$ yrs) is governed by differential precession and to a lesser extent, the Kozai-Lidov effect. In our proposed model the misalignment observed between the outer and inner disc arises naturally as a result of the misaligned outer companion driving the outer disc to precess more rapidly than the inner disc., Comment: Accepted for publication in MNRAS, 12 pages, 6 figures. Movie of simulations available at https://youtu.be/MvH8TO21W_4

In recent years, several protoplanetary discs have been observed to exhibit spirals, both in scattered light and (sub)millimetre continuum data. The HD 100453 binary star system hosts such a disc around its primary. Previous work has argued that the spirals were caused by the gravitational interaction of the secondary, which was assumed to be on a circular orbit, coplanar with the disc (meaning here the large outer disc, as opposed to the very small inner disc). However, recent observations of the CO gas emission were found incompatible with this assumption. In this paper, we run SPH simulations of the gas and dust disc for seven orbital configurations taken from astrometric fits and compute synthetic observations from their results. Comparing to high-resolution ALMA $^{12}$CO data, we find that the best agreement is obtained for an orbit with eccentricity $e=0.32$ and semi-major axis $a=207$ au, inclined by $61^\circ$ relative to the disc plane. The large misalignment between the disc and orbit planes is compatible with the tidal evolution of a circumprimary disc in an eccentric, unequal-mass binary star., Comment: Accepted for publication in MNRAS. 21 pages, 28 figures

Trojans are defined as objects that share the orbit of a planet at the stable Lagrangian points $L_4$ and $L_5$. In the Solar System, these bodies show a broad size distribution ranging from micrometer($\mu$m) to centimeter(cm) particles (Trojan dust) and up to kilometer (km) rocks (Trojan asteroids). It has also been theorized that earth-like Trojans may be formed in extra-solar systems. The Trojan formation mechanism is still under debate, especially theories involving the effects of dissipative forces from a viscous gaseous environment. We perform hydro-simulations to follow the evolution of a protoplanetary disk with an embedded 1--10 Jupiter-mass planet. On top of the gaseous disk, we set a distribution of $\mu$m--cm dust particles interacting with the gas. This allows us to follow dust dynamics as solids get trapped around the Lagrangian points of the planet. We show that large vortices generated at the Lagrangian points are responsible for dust accumulation, where the leading Lagrangian point $L_4$ traps a larger amount of submillimeter (submm) particles than the trailing $L_5$, which traps mostly mm--cm particles. However, the total bulk mass, with typical values of $\sim M_{\rm moon}$, is more significant in $L_5$ than in $L_4$, in contrast to what is observed in the current Solar System a few gigayears later. Furthermore, the migration of the planet does not seem to affect the reported asymmetry between $L_4$ and $L_5$. The main initial mass reservoir for Trojan dust lies in the same co-orbital path of the planet, while dust migrating from the outer region (due to drag) contributes very little to its final mass, imposing strong mass constraints for the in situ formation scenario of Trojan planets., Comment: 18 pages, 18 figures. Accepted for publication in Astronomy & Astrophysics

Spiral arms in protoplanetary discs are thought to be linked to the presence of companions. We test the hypothesis that the double spiral arm morphology observed in the transition disc MWC 758 can be generated by an $\approx 10$ M$_{\rm Jup}$ companion on an eccentric orbit internal to the spiral arms. Previous studies on MWC 758 have assumed an external companion. We compare simulated observations from three dimensional hydrodynamics simulations of disc-companion interaction to scattered light, infrared and CO molecular line observations, taking into account observational biases. The inner companion hypothesis is found to explain the double spiral arms, as well as several additional features seen in MWC 758 -- the arc in the northwest, substructures inside the spiral arms, the cavity in CO isotopologues, and the twist in the kinematics. Testable predictions include detection of fainter spiral structure, detection of a point source south-southeast of the primary, and proper motion of the spiral arms., Comment: 13 pages, 6 figures, resubmitted to MNRAS

In this work we demonstrate that the inner spiral structure observed in AB Aurigae can be created by a binary star orbiting inside the dust cavity. We find that a companion with a mass-ratio of 0.25, semi-major axis of 40 au, eccentricity of 0.5, and inclination of 90{\deg} produces gaseous spirals closely matching the ones observed in $^{12}$CO (2-1) line emission. Based on dust dynamics in circumbinary discs (Poblete, Cuello, and Cuadra 2019), we constrain the inclination of the binary with respect to the circumbinary disc to range between 60{\deg} and 90{\deg}. We predict that the stellar companion is located roughly 0.18 arcsec from the central star towards the east-southeast, above the plane of the disc. Should this companion be detected in the near future, our model indicates that it should be moving away from the primary star at a rate of 6 mas/yr on the plane of the sky. Since our companion is inclined, we also predict that the spiral structure will appear to change with time, and not simply co-rotate with the companion., Comment: Accepted for publication in MNRAS. 11 Pages, 7 Figures

FU Orionis objects are low-mass pre-main sequence stars characterized by dramatic outbursts of several magnitudes in brightness. These outbursts are linked to episodic accretion events in which stars gain a significant portion of their mass. The physical processes behind these accretion events are not yet well understood. The archetypical FU Ori system, FU Orionis, is composed of two young stars with detected gas and dust emission. The continuum emitting regions have not been resolved until now. Here, we present 1.3 mm observations of the FU Ori binary system with ALMA. The disks are resolved at 40 mas resolution. Radiative transfer modeling shows that the emission from FU Ori north (primary) is consistent with a dust disk with a characteristic radius of $\sim$11 au. The ratio between major and minor axes shows that the inclination of the disk is $\sim$37 deg. FU Ori south is consistent with a dust disk of similar inclination and size. Assuming the binary orbit shares the same inclination angle as the disks, the deprojected distance between north and south components is 0.6'', i.e. $\sim$250 au. Maps of $^{12}$CO emission show a complex kinematic environment with signatures disk rotation at the location of the northern component, and also (to a lesser extent) for FU Ori south. The revised disk geometry allows us to update FU Ori accretion models (Zhu et al.), yielding a stellar mass and mass accretion rate of FU Ori north of 0.6 M$_{\odot}$ and 3.8$\times10^{-5}$ M$_{\odot}$ yr$^{-1}$, respectively., Comment: 9 pages, 4 figures, accepted for publication in ApJ

We now have several observational examples of misaligned broken protoplanetary discs, where the disc inner regions are strongly misaligned with respect to the outer disc. Current models suggest that this disc structure can be generated with an internal misaligned companion (stellar or planetary), but the occurrence rate of these currently unobserved companions remains unknown. Here we explore whether a strong misalignment between the inner and outer disc can be formed without such a companion. We consider a disc that has an existing gap --- essentially separating the disc into two regions --- and use a flyby to disturb the discs, leading to a misalignment. Despite considering the most optimistic parameters for this scenario, we find maximum misalignments between the inner and outer disc of $\sim$45$^{\circ}$ and that these misalignments are short-lived. We thus conclude that the currently observed misaligned discs must harbour internal, misaligned companions., Comment: Accepted for publication in MNRAS. 9 pages with 6 figures

Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa)., Comment: 11 pages, 5 figures, accepted for publication in MNRAS

We test the hypothesis that the disc cavity in the `transition disc' Oph IRS 48 is carved by an unseen binary companion. We use 3D dust-gas smoothed-particle hydrodynamics simulations to demonstrate that marginally coupled dust grains concentrate in the gas over-density that forms in in the cavity around a low binary mass ratio binary. This produces high contrast ratio dust asymmetries at the cavity edge similar to those observed in the disc around IRS 48 and other transition discs. This structure was previously assumed to be a vortex. However, we show that the observed velocity map of IRS 48 displays a peculiar asymmetry that is not predicted by the vortex hypothesis. We show the unusual kinematics are naturally explained by the non-Keplerian flow of gas in an eccentric circumbinary cavity. We further show that perturbations observed in the isovelocity curves of IRS 48 may be explained as the product of the dynamical interaction between the companion and the disc. The presence of a $\sim$0.4 M$_{\odot}$ companion at a $\sim$10 au separation can qualitatively explain these observations. High spatial resolution line and continuum imaging should be able to confirm this hypothesis., Comment: 9 pages, 7 figures, accepted for publication in MNRAS

Recent observations have revealed that protoplanetary discs often exhibit cavities and azimuthal asymmetries such as dust traps and clumps. The presence of a stellar binary system in the inner disc regions has been proposed to explain the formation of these structures. Here, we study the dust and gas dynamics in circumbinary discs around eccentric and inclined binaries. This is done through two-fluid simulations of circumbinary discs, considering different values of the binary eccentricity and inclination. We find that two kinds of dust structures can form in the disc: a single horseshoe-shaped clump, on top of a similar gaseous over-density; or numerous clumps, distributed along the inner disc rim. The latter features form through the complex interplay between the dust particles and the gaseous spirals caused by the binary. All these clumps survive between one and several tens of orbital periods at the feature location. We show that their evolution strongly depends on the gas-dust coupling and the binary parameters. Interestingly, these asymmetric features could in principle be used to infer or constrain the orbital parameters of a stellar companion - potentially unseen - inside the inner disc cavity. Finally, we apply our findings to the disc around AB Aurigae., Comment: 12 pages, 13 figures, accepted for publication in MNRAS

Dynamical studies suggest that most of the circumbinary discs (CBDs) should be coplanar. However, under certain initial conditions, the CBD can evolve toward polar orientation. Here we extend the parametric study of polar configurations around detached close-in binaries through $N$-body simulations. For polar configurations around binaries with mass ratios $q$ below $0.7$, the nominal location of the mean motion resonance (MMR) $1~:~4$ predicts the limit of stability for $e_{B} > 0.1$. Alternatively, for $e_{B} < 0.1$ or $q \sim 1$, the nominal location of the MMR $1~:~3$ is the closest stable region. The presence of a} giant planet increases the region of forbidden polar configurations around low mass ratio binaries with eccentricities $e_B\sim0.4$ with respect to rocky earth-like planets. For equal mass stars, the eccentricity excitation $\Delta e$ of polar orbits smoothly increases with decreasing distance to the binary. For $q<1$, $\Delta e$ can reach values as high as $0.4$. Finally, we studied polar configurations around $HD~98800BaBb$ and show that the region of stability is strongly affected by the relative positions of the nodes. The most stable configurations in the system correspond to polar particles, which are not expected to survive on longer time-scales due to the presence of the external perturber HD~$98800AaAb$., Comment: 9 pages. 7 figures. Accepted for publication in the special issue of the Journal of Physics: Conference Series aimed to the works presented at the Brazilian Colloquium on Orbital Dynamics - CBDO 2018

Close-in giant planets represent the most significant evidence of planetary migration. If large exomoons form around migrating giant planets which are more stable (e.g. those in the Solar System), what happens to these moons after migration is still under intense research. This paper explores the scenario where large regular exomoons escape after tidal-interchange of angular momentum with its parent planet, becoming small planets by themselves. We name this hypothetical type of object a \textit{ploonet}. By performing semi-analytical simulations of tidal interactions between a large moon with a close-in giant, and integrating numerically their orbits for several Myr, we found that in $\sim$50 per cent of the cases a young ploonet may survive ejection from the planetary system, or collision with its parent planet and host star, being in principle detectable. Volatile-rich ploonets are dramatically affected by stellar radiation during both planetocentric and siderocentric orbital evolution, and their radius and mass change significantly due to the sublimation of most of their material during time-scales of hundred of Myr. We estimate the photometric signatures that ploonets may produce if they transit the star during the phase of evaporation, and compare them with noisy lightcurves of known objects ("Kronian" stars and non-periodical dips in dusty lightcurves). Additionally, the typical transit timing variations (TTV) induced by the interaction of a ploonet with its planet are computed. We find that present and future photometric surveys' capabilities can detect these effects and distinguish them from those produced by other nearby planetary encounters., Comment: Published in Monthly Notices of the Royal Astronomical Society. 11 pages, 7 figures

Context. Dynamical studies suggest that most circumbinary discs (CBDs) should be coplanar. However, some theoretical works show that under certain conditions a CBD can become polar, which means that its rotation vector is orthogonal with respect to the binary orbital plane. Interestingly, very recent observations show that polar CBDs exist in nature (e.g. HD 98800). Aims. We test the predictions of CBD alignment around eccentric binaries based on linear theory. In particular, we compare prograde and retrograde CBD configurations. Then, assuming planets form in these systems, we thoroughly characterise the orbital behaviour and stability of misaligned (P-type) particles. This is done for massless and massive particles. Methods. The evolution of the CBD alignment for various configurations was modelled through three-dimensional hydrodynamical simulations. For the orbital characterisation and the analysis stability, we relied on long-term N-body integrations and structure and chaos indicators, such as $\Delta e$ and MEGNO. Results. We confirm previous analytical predictions on CBD alignment, but find an unexpected symmetry breaking between prograde and retrograde configurations. Therefore, the likelihood of becoming polar for a highly misaligned CBD is higher than previously thought. Regarding the stability of circumbinary P-type planets (also know as Tatooines), polar orbits are stable over a wide range of binary parameters. In particular, for binary eccentricities below 0.4 the orbits are stable for any value of the binary mass ratio. In the absence of gas, planets with masses below $10^{-5}\,M_{\odot}$ have negligible effects on the binary orbit. Finally, we suggest that mildly eccentric equal-mass binaries should be searched for polar Tatooines., Comment: 9 pages, 9 figures; Astronomy & Astrophysics (in press)

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade., Comment: White paper submitted to the Astro2020 decadal survey

We present 3D smoothed particle hydrodynamics simulations of protoplanetary discs undergoing a flyby by a stellar perturber on a parabolic orbit lying in a plane inclined relative to the disc mid-plane. We model the disc as a mixture of gas and dust, with grains ranging from 1 {\mu}m to 10 cm in size. Exploring different orbital inclinations, periastron distances and mass ratios, we investigate the disc dynamical response during and after the flyby. We find that flybys induce evolving spiral structure in both gas and dust which can persist for thousands of years after periastron. Gas and dust structures induced by the flyby differ because of drag-induced effects on the dust grains. Variations in the accretion rate by up to an order of magnitude occur over a time-scale of order 10 years or less, inducing FU Orionis-like outbursts. The remnant discs are truncated and warped. The dust disc is left more compact than the gas disc, both because of disc truncation and accelerated radial drift of grains induced by the flyby., Comment: 27 pages, 24 figures, accepted for publication in MNRAS

Context. Despite the recent discovery of spiral-shaped features in protoplanetary discs in the near-infrared and millimetric wavelengths, there is still an active discussion to understand how they formed. In fact, the spiral waves observed in discs around young stars can be due to different physical mechanisms: planet/companion torques, gravitational perturbations or illumination effects. Aims. We study the spirals formed in the gaseous phase due to two diametrically opposed shadows cast at fixed disc locations. The shadows are created by an inclined non-precessing disc inside the cavity, which is assumed to be optically thick. In particular, we analyse the effect of these spirals on the dynamics of the dust particles and discuss their detectability in transition discs. Methods. We perform gaseous hydrodynamical simulations with shadows, then we compute the dust evolution on top of the gaseous distribution, and finally we produce synthetic ALMA observations of the dust emission based on radiative transfer calculations. Results. Our main finding is that mm- to cm-sized dust particles are efficiently trapped inside the shadow-triggered spirals. We also observe that particles of various sizes starting at different stellocentric distances are well mixed inside these pressure maxima. This dynamical effect would favour grain growth and affect the resulting composition of planetesimals in the disc. In addition, our radiative transfer calculations show spiral patterns in the disc at 1.6 {\mu}m and 1.3 mm. Due to their faint thermal emission (compared to the bright inner regions of the disc) the spirals cannot be detected with ALMA. Our synthetic observations prove however that shadows are observable as dips in the thermal emission., Comment: 15 pages, 11 figures, accepted for publication in A&A

We present a methodology to interpret observations of protoplanetary discs where a flyby, also called a tidal encounter, is suspected. In case of a flyby, protoplanetary discs can be significantly disturbed. The resulting dynamical and kinematical signatures can last for several thousands of years after the flyby and hence deeply affect the evolution of the disc. These effects are stronger for closer encounters and more massive perturbers. For the very same flyby parameters, varying the inclination of the perturber's orbit produces a broad range of disc structures: spirals, bridges, warps and cavities. We study this kind of features both in the gas and in the dust for grains ranging from 1 {\mu}m to 10 cm in size. Interestingly, the dust exhibits a different dynamical behaviour compared to the gas because of gas-drag effects. Finally, flybys can also trigger high accretion events in the disc-hosting star, readily similar to FU Orionis-type outbursts. All this information can be used to infer the flyby parameters from an incomplete set of observations at different wavelengths. Therefore, the main scope of our flyby scene investigation (FSI) methodology is to help to interpret recent "puzzling" disc observations., Comment: Proceedings of the SF2A 2018 Conference; 7 pages, 5 figures

We present 3D hydrodynamical models of the HD142527 protoplanetary disc, a bright and well studied disc that shows spirals and shadows in scattered light around a 100 au gas cavity, a large horseshoe dust structure in mm continuum emission, together with mysterious fast radial flows and streamers seen in gas kinematics. By considering several possible orbits consistent with the observed arc, we show that all of the main observational features can be explained by one mechanism - the interaction between the disc and the observed binary companion. We find that the spirals, shadows and horseshoe are only produced in the correct position angles by a companion on an inclined and eccentric orbit approaching periastron - the 'red' family from Lacour et al. (2016). Dust-gas simulations show radial and azimuthal concentration of dust around the cavity, consistent with the observed horseshoe. The success of this model in the HD142527 disc suggests other mm-bright transition discs showing cavities, spirals and dust asymmetries may also be explained by the interaction with central companions., Comment: 16 pages, 12 figures, accepted to MNRAS. Movies at http://users.monash.edu.au/~dprice/pubs/HD142527/

Shadows and spirals seem to be common features of transition discs. Among the spiral-triggering mechanisms proposed, only one establishes a causal link between shadows and spirals so far. In fact, provided the presence of shadows in the disc, the combined effect of temperature gradient and differential disc rotation, creates strong azimuthal pressure gradients. After several thousand years, grand-design spirals develop in the gas phase. Previous works have only considered static shadows caused by an inclined inner disc. However, in some cases, the inner regions of circumbinary discs can break and precess. Thus, it is more realistic to consider moving shadow patterns in the disc. In this configuration, the intersection between the inner and the outer discs defines the line of nodes at which the shadows are cast. Here, we consider moving shadows and study the resulting circumbinary disc structure. We find that only static and prograde shadows trigger spirals, in contrast to retrograde ones. Interestingly, if a region of the disc corotates with the shadow, a planet-like signature develops at the co-rotation position. The resulting spirals resemble those caused by a planet embedded in the disc, with similar pitch angles., Comment: 6 pages, 4 figures. Accepted for publication in MNRAS letters