Planet-disk Interactions Research Articles

Mind the Kinematics Simulation of Planet–Disk Interactions: Time Evolution and Numerical Resolution

Planet–disk interactions can produce kinematic signatures in protoplanetary disks. While recent observations have detected non-Keplerian gas motions in disks, their origins are still being debated. To explore this, we conduct 3D hydrodynamic simulations using the code FARGO3D to study nonaxisymmetric kinematic perturbations at two scale heights induced by Jovian planets in protoplanetary disks, followed by examinations of detectable signals in synthetic CO emission line observations at millimeter wavelengths. We advocate for using residual velocity or channel maps, generated by subtracting an azimuthally averaged background of the disk, to identify planet-induced kinematic perturbations. We investigate the effects of two basic simulation parameters, simulation duration and numerical resolution, on the simulation results. Our findings suggest that a short simulation (e.g., 100 orbits) is insufficient to establish a steady velocity pattern given our chosen viscosity (α = 10−3) and displays plenty of fluctuations on an orbital timescale. Such transient features could be detected in observations. By contrast, a long simulation (e.g., 1000 orbits) is required to reach steady state in kinematic structures. At 1000 orbits, the strongest detectable velocity structures are found in the spiral wakes close to the planet. Through numerical convergence tests, we find hydrodynamics results converge in spiral regions at a resolution of 14 cells per disk scale height or higher. Meanwhile, synthetic observations produced from hydrodynamic simulations at different resolutions are indistinguishable with 0.″1 angular resolution and 10 hr of integration time on Atacama Large Millimeter/submillimeter Array.

The Search for Disk Perturbing Planets Around the Asymmetrical Debris Disk System HD 111520 Using REBOUND

Debris disks, which are optically thin, dusty disks around main-sequence stars, are often found to have structures and/or asymmetries associated with planet–disk interactions. Debris disk morphologies can hence be used as probes for planets in these systems, which are unlikely to be detected with other current exoplanet detection methods. In this study we take a look at the very asymmetrical debris disk around HD 111520, which harbours several signs of perturbation such as a “fork”-like structure in the NW, as well as a 4° warp from the midplane on either side of the disk. We simulate the complicated disk morphology using the code REBOUND, with the goal of constraining the possible mass and orbit of the planet responsible for the observed structures. We find that an ∼1 M jup, eccentric planet that is inclined relative to the disk and is orbiting outside the warp location is able to reproduce the majority of disk features including the warp, fork, and radial extent asymmetry. To create the surface brightness asymmetry, a second eccentric planet is required inside the disk inner edge (50 au), although we are unable to produce the 2:1 brightness asymmetry observed, suggesting that a second mechanism may be required. Our work demonstrates how debris disk morphologies alone can be used to learn more about the architecture and evolution of a system as a whole, and can provide planet constraints to determine potential targets for current/future instruments such as JWST/NIRCam and the Gemini Planet Imager 2.0.

A Search for Collisions and Planet–Disk Interactions in the Beta Pictoris Disk with 26 Years of High-precision HST/STIS Imaging

β Pictoris's well-studied debris disk and two known giant planets, in combination with the stability of the Hubble Space Telescope’s Space Telescope Imaging Spectrograph (HST/STIS) (and now also JWST), offers a unique opportunity to test planet–disk interaction models and to observe recent planetesimal collisions. We present HST/STIS coronagraphic imaging from two new epochs of data taken between 2021 and 2023, complementing earlier data taken in 1997 and 2012. This data set enables a temporal comparison with the longest baseline and highest precision of any debris disk to date, with sensitivity to variations in temporal surface brightness of sub-percent levels in the midplane of the disk. While no localized changes in surface brightness are detected, which would be indicative of a recent planetesimal collision, there is a tentative brightening of the southeast side of the disk over the past decade. We link the constraints on surface brightness variations to dynamical models of the planetary system’s evolution and to the collisional history of planetesimals. Using a coupled collisional model and injection/recovery framework, we estimate sensitivity to expanding collisional debris down to a Ceres mass per progenitor in the most sensitive regions of the disk midplane. These results demonstrate the capabilities of long-baseline, temporal studies with HST (and also soon with JWST) for constraining the physical processes occurring within debris disks.

Spiral excitation in protoplanetary disks through gap-edge illumination

The advent of high-resolution, near-infrared (NIR) instruments such as VLT/SPHERE and Gemini/GPI has helped uncover a wealth of substructure in planet-forming disks, including large, prominent spiral arms in MWC 75 8, SAO 206462, and V1247 Ori. In the classical theory of disk-planet interaction, these arms are consistent with Lindblad-resonance driving by companions of multiple Jupiter masses. Despite improved detection limits, evidence for massive bodies like this in connection with spiral substructure has been inconclusive. In search of an alternative explanation, we used the PLUTO code to run 3D hydrodynamical simulations with two comparatively low planet masses (Saturn mass and Jupiter mass) and two thermodynamic prescriptions (three-temperature radiation hydrodynamics, and the more traditional β-cooling) in a low-mass disk. In the radiative cases, an m = 2 mode, potentially attributable to the interaction of stellar radiation with gap-edge asymmetries, creates an azimuthal pressure gradient, which in turn gives rise to prominent spiral arms in the upper layers of the disk. Monte Carlo radiative transfer post-processing with RADMC3D revealed that in NIR scattered light, these gap-edge spirals are significantly more prominent than the traditional Lindblad spirals for planets in the mass range we tested. Our results demonstrate that even intermediate-mass protoplanets, which are less detectable, but more ubiquitous than super-Jupiters, are capable of indirectly inducing large-scale spiral disk features, and underscore the importance of including radiation physics in any efforts to reproduce observations.

Horseshoes and spiral waves: capturing the 3D flow induced by a low-mass planet analytically

ABSTRACT The key difficulty faced by 2D models for planet–disc interaction is in appropriately accounting for the impact of the disc’s vertical structure on the dynamics. 3D effects are often mimicked via softening of the planet’s potential; however, the planet-induced flow and torques often depend strongly on the choice of softening length. We show that for a linear adiabatic flow perturbing a vertically isothermal disc, there is a particular vertical average of the 3D equations of motion that exactly reproduces 2D fluid equations for arbitrary adiabatic index. There is a strong connection here with the Lubow–Pringle 2D mode of the disc. Correspondingly, we find a simple, general prescription for the consistent treatment of planetary potentials embedded within ‘2D’ discs. The flow induced by a low-mass planet involves large-scale excited spiral density waves that transport angular momentum radially away from the planet and ‘horseshoe streamlines’ within the coorbital region. We derive simple linear equations governing the flow that locally capture both effects faithfully simultaneously. We present an accurate coorbital flow solution allowing for inexpensive future study of corotation torques, and predict the vertical structure of the coorbital flow and horseshoe region width for different values of adiabatic index, as well as the vertical dependence of the initial shock location. We find strong agreement with the flow computed in 3D numerical simulations, and with 3D one-sided Lindblad torque estimates, which are a factor of 2–3 lower than values from previous 2D simulations.

Rocky planet formation in compact disks around M dwarfs

Due to the improvements in radial velocity and transit techniques, we know that rocky or rocky-icy planets, in particular close-in super-Earths in compact configurations, are the most common ones around M dwarfs. On the other hand, thanks to the high angular resolution of ALMA we know that many disks around very low-mass stars (between 0.1 and 0.5 M$_ are rather compact and small (without observable substructures and radius less than 20 au), which favors the idea of an efficient radial drift that could enhance planet formation in the terrestrial zone. Our aim was to investigate the potential formation paths of the observed close-in rocky exoplanet population around M dwarfs, especially close-in super-Earths, assuming that planet formation could take place in compact disks with an efficient dust radial drift. We developed N-body simulations that include a sample of embryos growing by pebble accretion exposed to planet-disk interactions, star-planet tidal interactions, and general relativistic corrections that include the evolution of the luminosity, radius, and rotational period of the star. For a star of 0.1 M$_ we considered different gas disk viscosities and initial embryo distributions. We also explored planet formation by pebble accretion around stars of 0.3 M$_ and 0.5 M$_ Lastly, for each stellar mass, we ran simulations that include a sample of embryos growing by planetesimal accretion instead of pebble accretion. Our main result is that the sample of simulated planets that grow by pebble accretion in a gas disk with low viscosity ($ $) can reproduce the close-in low-mass exoplanet population around M dwarfs in terms of multiplicity, mass, and semi-major axis. Furthermore, we found that a gas disk with high viscosity ($ $), and thus lower pebble accretion rates, cannot reproduce the observed planet masses as no planet more massive than 0.5 M$_ could be formed in our simulations. In addition, we show that planetesimal accretion favors the formation of smaller planets than pebble accretion does. Whether this planet population truly exists remains unknown with the current instrumental sensitivity. Rocky planet formation around M dwarfs can take place in compact and small dust disks driven by an efficient radial drift in a gas disk with low viscosity ($ $). This result points toward a new approach in the direction of the disk conditions needed for rocky planet formation around very low-mass stars.

Shoulder of Dust Rings Formed by Planet–Disk Interactions

Recent analyses of millimeter-wavelength protoplanetary disk observations have revealed several emission excesses on the previously identified dust rings, referred to as dust shoulders. The prevalence of dust shoulders suggests that they trace a common but unclear mechanism. In this work, we combine 3D, multifluid hydrodynamic simulations with radiative transfer calculations to explain the formation of dust shoulders. We find that the ring-shoulder pairs can result from the 3D planet–disk interactions with massive, gap-opening planets. The key driver is the dust filtration effect at the local pressure maximum due to planet-driven outward gas flows. Our work provides a possible explanation for the outer dust shoulders in recent super-resolution analyses of Atacama Large Millimeter/submillimeter Array observations. It also provides insights into the formation of the inner dust shoulder in the PDS 70 disk and highlights the role of 3D effects in planet–disk interaction studies.

Three-temperature radiation hydrodynamics with PLUTO: Thermal and kinematic signatures of accreting protoplanets

In circumstellar disks around young stars, the gravitational influence of nascent planets produces telltale patterns in density, temperature, and kinematics. To better understand these signatures, we first performed 3D hydrodynamical simulations of a 0.012 M⊙ disk with a Saturn-mass planet orbiting circularly in-plane at 40 au. We tested four different disk thermodynamic prescriptions (in increasing order of complexity: local isothermality, β cooling, two-temperature radiation hydrodynamics, and three-temperature radiation hydrodynamics), finding that β cooling offers a reasonable approximation for the three-temperature approach when the planet is not massive or luminous enough to substantially alter the background temperature and density structure. Thereafter, using the three-temperature scheme, we relaxed this assumption, simulating a range of different planet masses (Neptune-mass, Saturn-mass, and Jupiter-mass) and accretion luminosities (0 and 10−3 L⊙) in the same disk. Our investigation revealed that signatures of disk–planet interaction strengthen with increasing planet mass, with circumplanetary flows becoming prominent in the high-planet-mass regime. Accretion luminosity, which adds pressure support around the planet, was found to weaken the midplane Doppler flip, which is potentially visible in optically thin tracers such as C18O, while strengthening the spiral signature, particularly in upper disk layers sensitive to thicker lines, such as those of 12CO.

Buoyancy torques prevent low-mass planets from stalling in low-turbulence radiative discs

ABSTRACT Low-mass planets migrating inwards in laminar protoplanetary discs (PPDs) experience a dynamical corotation torque (DCT), which is expected to slow down migration to a stall. However, baroclinic effects can reduce or even reverse this effect, leading to rapid inward migration. In the radiatively inefficient inner disc, one such mechanism is the buoyancy response of the disc to an embedded planet. Recent work has suggested that radiative cooling can quench this response, but for parameters that are not necessarily representative of the inner regions of PPDs. We perform global 3D inviscid radiation hydrodynamics simulations of planet–disc interaction to investigate the effect of radiative cooling on the buoyancy-driven torque in a more realistic disc model. We find that the buoyancy response exerts a negative DCT – albeit partially damped due to radiative cooling – resulting in sustained, rapid inward migration. Models that adopt a local cooling prescription significantly overestimate the impact of the buoyancy response, highlighting the importance of a realistic treatment of radiation transport that includes radiative diffusion. Our results suggest that low-mass planets should migrate inwards faster than has been previously expected in radiative discs, with implications for the formation and orbital distribution of super-Earths and sub-Neptunes at intermediate distances from their host stars, unless additional physical processes that can slow down migration are considered.

The interplay between forming planets and photoevaporating discs

Context. Disc winds and planet–disc interactions are two crucial mechanisms that define the structure, evolution, and dispersal of protoplanetary discs. While winds are capable of removing material from discs, eventually leading to their dispersal, massive planets can shape their disc by creating sub-structures such as gaps and spiral arms. Aims. We studied the interplay between an X-ray photoevaporative disc wind and the sub-structures generated due to planet–disc interactions to determine how their mutual interactions affect the disc’s and the planet’s evolution. Methods. We performed 3D hydrodynamic simulations of viscous discs (α = 6.9 × 10−4) that host a Jupiter-like planet and undergo X-ray photoevaporation. We traced the gas flows within the disc and wind and measured the rate of accretion onto the planet, as well as the gravitational torque that is acting on it. Results. Our results show that the planetary gap removes the wind’s pressure support, allowing wind material to fall back into the gap. This opens new pathways for material from the inner disc (and part of the outer disc) to be redistributed through the wind towards the gap. Consequently, the gap becomes shallower and the flow of mass across the gap in both directions is significantly increased, as is the planet’s mass-accretion rate (by factors of ≈5 and ≈2, respectively). Moreover, the wind-driven redistribution results in a denser inner disc and a less dense outer disc, which, combined with the recycling of a significant portion of the inner wind, leads to longer lifetimes for the inner disc, contrary to the expectation in a planet-induced photoevaporation scenario that has been proposed in the past.

Increasing the raw contrast of VLT/SPHERE with the dark hole technique

Context. Imaging exoplanetary systems is essential to characterizing exoplanet atmospheres and orbits as well as circumstellar disks and to studying planet-disk interactions to understand the planet formation processes. Imaging exoplanets or circumstellar disks in the visible and near-infrared is challenging, however, because these objects are very faint relative to their star, even though only fractions of an arcsecond away. Coronagraphic instruments have already allowed the imaging of a few exoplanets, but their performance is limited by wavefront aberrations. Adaptive optics systems partly compensate for the Earth’s atmosphere turbulence, but they cannot fully control the wavefront. Some of the starlight leaks through the coronagraph and forms speckles in the astrophysical image. Focal plane wavefront control, used as a second stage after the adaptive optics system, has been proposed to minimize the speckle intensity within an area called the dark hole. Aims. We previously demonstrated the on-sky performance of dark hole techniques, pairwise probing coupled with electric field conjugation, using the apodized pupil Lyot coronagraph of the VLT/SPHERE instrument. In this paper, we probe their performance using the SPHERE four-quadrant phase mask coronagraph, and we demonstrate the interest of combining dark hole techniques and reference differential imaging. Methods. We used these dark hole techniques on-sky to create a dark hole in the narrow band around 1.7 |j.m observing HR 4796. We then recorded broadband images of HR 4796 and a reference star at the H band. Results. The dark hole techniques improved the H-band detection limit by a factor of three. The dark hole was stable from one star to a nearby star enabling reference differential imaging. Conclusions. This stability offers two new strategies of observation. First, one can quickly create a dark hole observing a bright star before pointing to a faint target star. Furthermore, one can couple dark hole techniques and reference differential imaging. A very interesting point is that the performance of these methods does not depend on the astrophysical signal.

A Recipe for Eccentricity and Inclination Damping for Partial-gap Opening Planets in 3D Disks

In a previous paper, we showed that, like the migration speed, the eccentricity damping efficiency is modulated linearly by the depth of the partial gap a planet carves in the disk surface density profile, resulting in less efficient e-damping compared to the prescription commonly used in population synthesis works. Here, we extend our analysis to 3D, refining our e-damping formula and studying how the inclination damping efficiency is also affected. We perform high-resolution 3D locally isothermal hydrodynamical simulations of planets with varying masses embedded in disks with varying aspect ratios and viscosities. We extract the gap profile and orbital damping timescales for fixed eccentricities and inclinations up to the disk scale height. The limit in gap depths below which vortices appear, in the low-viscosity case, happens roughly at the transition between classical type-I and type-II migration regimes. The orbital damping timescales can be described by two linear trends with a break around gap depths ∼80% and with slopes and intercepts depending on the eccentricity and inclination. These trends are understood on physical grounds and are reproduced by simple fitting formulas whose error is within the typical uncertainty of type-I torque formulas. Thus, our recipes for the gap depth and orbital damping efficiencies yield a simple description for planet–disk interactions to use in N-body codes in the case of partial-gap opening planets that is consistent with high-resolution 3D hydrosimulations. Finally, we show examples of how our novel orbital damping prescription can affect the outcome of population synthesis experiments.

ALMA 2D super-resolution imaging of Taurus–Auriga protoplanetary disks: Probing statistical properties of disk substructures

Abstract In the past decade, ALMA observations of protoplanetary disks revealed various substructures including gaps and rings. Their origin of substructures may be probed through statistical studies of their physical properties. We present the analyses of archival ALMA Band 6 continuum data of 43 disks (39 Class II and four Herbig Ae) in the Taurus–Auriga region. We employ a novel 2D super-resolution imaging technique based on sparse modeling to obtain images with high fidelity and spatial resolution. As a result, we have obtained images with spatial resolutions comparable to a few au (${0_{.}^{\prime\prime}02}$–${0_{.}^{\prime\prime}1}$), which is two to three times better than conventional CLEAN methods. All dust disks are spatially resolved, with the radii ranging from 8 to 238 au with a median radius of 45 au. Half of the disks harbor clear gap structures, the radial locations of which show a bimodal distribution with peaks at ≲20 au and ≳30 au. We also see structures indicating weak gaps at all the radii in the disk. We find that the widths of these gaps increase with their depths, which is consistent with the model of planet–disk interactions. The inferred planet mass–orbital radius distribution indicates that the planet distribution is analogous to our solar system. However, planets with Neptune mass or lower may exist in all the radii.

MINDS: JWST/NIRCam imaging of the protoplanetary disk PDS 70

Context. Two protoplanets have recently been discovered within the PDS 70 protoplanetary disk. JWST/NIRCam offers a unique opportunity to characterize them and their birth environment at wavelengths that are difficult to access from the ground. Aims. We image the circumstellar environment of PDS 70 at 1.87 μm and 4.83 μm, assess the presence of Pa-α emission due to accretion onto the protoplanets, and probe any IR excess indicative of heated circumplanetary material. Methods. We obtained noncoronagraphic JWST/NIRCam images of PDS 70 within the MIRI mid-INfrared Disk Survey (MINDS) program. We leveraged the Vortex Image Processing (VIP) package for data reduction, and we developed dedicated routines for optimal stellar point spread function subtraction, unbiased imaging of the disk, and protoplanet flux measurement in this type of dataset. A radiative transfer model of the disk was used to separate the contributions from the disk and the protoplanets. Results. We redetect both protoplanets and identify extended emission after subtracting a disk model, including a large-scale spiral-like feature. We interpret its signal in the direct vicinity of planet c as tracing the accretion stream that feeds its circumplanetary disk, while the outer part of the feature may rather reflect asymmetric illumination of the outer disk. We also report a bright signal that is consistent with a previously proposed protoplanet candidate enshrouded in dust near the 1:2:4 mean-motion resonance with planets b and c. The 1.87 μm flux of planet b is consistent with atmospheric model predictions, but the flux of planet c is not. We discuss potential origins for this discrepancy, including significant Pa-α line emission. The 4.83 μm fluxes of planets b and c suggest enshrouding dust or heated CO emission from their circumplanetary environment. Conclusions. The use of image-processing methods that are optimized for extended disk signals on high-sensitivity and high-stability from JWST can uniquely identify signatures of planet–disk interactions and enable accurate photometry of protoplanets at wavelengths that are difficult to probe from the ground. Our results indicate that more protoplanets can be identified and characterized in other JWST datasets.

Forming localized dust concentrations in a dust ring: DM Tau case study

Context. Previous high-angular-resolution 225 GHz (~1.3 mm) continuum observations of the transitional disk DM Tau have resolved an outer ring at 20–120 au radii that is weakly azimuthally asymmetric. Aims. We aim to examine dust growth and filtration in the outer ring of DM Tau. Methods. We performed ~0″.06 (~8.7 au) resolution Karl G. Jansky Very Large Array (JVLA) 40–48 GHz (~7 mm; Q band) continuum observations, along with complementary observations at lower frequencies. In addition, we analyzed the archival JVLA observations undertaken since 2010. Results. Intriguingly, the Q band image resolved the azimuthally highly asymmetric, knotty dust emission sources close to the inner edge of the outer ring. Fitting the 8–700 GHz spectral energy distribution (SED) with two dust components indicates that the maximum grain size (amax) in these knotty dust emission sources is likely ≳300 µm, whereas it is ≳50 µm in the rest of the ring. These results may be explained by a trapping of inwardly migrating “grown” dust close to the ring inner edge. The exact mechanism for developing the azimuthal asymmetry has not yet been identified, which may be due to planet-disk interaction that might also be responsible for the creation of the dust cavity and pressure bump. Otherwise, it may be due to the fluid instabilities and vortex formation as a result of shear motions. Finally, we remark that the asymmetries in DM Tau are difficult to diagnose from the ≳225 GHz observations, owing to a high optical depth at the ring. In other words, the apparent symmetric or asymmetric morphology of the transitional disks may be related to the optical depths of those disks at the observing frequency.

Small and Large Dust Cavities in Disks around Mid-M Stars in Taurus

High angular resolution imaging by Atacama Large Millimeter/submillimeter Array (ALMA) has revealed the near universality and diversity of substructures in protoplanetary disks. However, disks around M-type pre-main-sequence stars are still poorly sampled, despite the prevalence of M dwarfs in the Galaxy. Here we present high-resolution (∼50 mas, 8 au) ALMA Band 6 observations of six disks around mid-M stars in Taurus. We detect dust continuum emission in all six disks, 12CO in five disks, and 13CO line in two disks. The size ratios between gas and dust disks range from 1.6 to 5.1. The ratio of about 5 for 2M0436 and 2M0450 indicates efficient dust radial drift. Four disks show rings and cavities, and two disks are smooth. The cavity sizes occupy a wide range: 60 au for 2M0412, and ∼10 au for 2M0434, 2M0436, and 2M0508. Detailed visibility modeling indicates that small cavities of 1.7 and 5.7 au may hide in the two smooth disks 2M0450 and CIDA 12. We perform radiative transfer fitting of the infrared spectral energy distributions to constrain the cavity sizes, finding that micron-sized dust grains may have smaller cavities than millimeter grains. Planet–disk interactions are the preferred explanation to produce the large 60 au cavity, while other physics could be responsible for the three ∼10 au cavities under current observations and theories. Currently, disks around mid- to late M stars in Taurus show a higher detection frequency of cavities than earlier-type stars, although a more complete sample is needed to evaluate any dependence of substructure on stellar mass.

PDS 70 unveiled by star-hopping: Total intensity, polarimetry, and millimeter imaging modeled in concert

Most ground-based direct-imaging planet search campaigns use angular differential imaging that distorts the signal from extended sources, such as protoplanetary disks. In the case of the young system PDS 70, for which two planets were detected within the cavity of a protoplanetary disk, obtaining a reliable image of both planets and the disk is essential to understanding planet-disk interactions. Our goals are to reveal the true intensity of the planets and disk without self-subtraction effects for the first time, search for new giant planets beyond separations of 0.1$''$, and to study the morphology of the disk shaped by two massive planets. We present $YJHK$-band imaging, polarimetry, and spatially resolved spectroscopy of PDS 70 using near-simultaneous reference star differential imaging, also known as star-hopping. We created a radiative transfer model of the system to try to match the near-infrared imaging and polarimetric data within measurement errors. Sub-millimeter imaging data from ALMA were also modeled. Furthermore, we extracted the spectra of the planets and the disk and compared them. With strong constraints, we find that the disk is quite flared, with a scale height of sim 15<!PCT!> at the outer edge of the disk at sim 90 au, similar to some disks in the literature. The gap inside sim 50 au is estimated to have sim 1<!PCT!> of the dust density of the outer disk. The northeast outer disk arc seen in previous observations is likely the outer lip of the flared disk. Abundance ratios of sub-micron, micron, and grains estimated by the modeling indicate a shallow grain-size index greater than $-2.7$, instead of the canonical -3.5. There is both vertical and radial segregation of grains. Planet $c$ is well separated from the disk and has a spectrum similar to planet $b$, and it is clearly redder than the disk spectra. Planet $c$ is possibly associated with the sudden flaring of the disk starting at sim 50 au. We found no new planets in the system. If we assume DUSTY models and an age of 5 Myrs, this indicates no new planets more massive than 5 outside a 12 au separation.

FARGOCPT: 2D Multiphysics code for simulating disk interactions with stars, planets, and particles

Context. Planet-disk interactions play a crucial role in the understanding of planet formation and disk evolution. There are multiple numerical tools available to simulate these interactions, including the commonly used FARGO code and its variants. Many of the codes have been extended over time to include additional physical processes, with a focus on their accurate modeling. Aims. We introduce FARGOCPT, an updated version of FARGO that incorporates other previous enhancements to the code, to provide a simulation environment tailored to studies of the interactions between stars, planets, and disks. It is meant to ensure an accurate representation of planet systems, hydrodynamics, and dust dynamics, with a focus on usability. Methods. The radiation-hydrodynamics part of FARGOCPT uses a second-order upwind scheme in 2D polar coordinates, supporting multiple equations of state, radiation transport, heating and cooling, and self-gravity. Shocks are considered using artificial viscosity. The integration of the N-body system is achieved by leveraging the REBOUND code. The dust module utilizes massless tracer particles, adapted to drag laws for the Stokes and Epstein regimes. Moreover, FARGOCPT provides mechanisms to simulate accretion onto stars and planets. Results. The code has been tested in practice in the context of multiple studies. Additionally, it comes with an automated test suite for checking the physics modules. It is available online. Conclusions. FARGOCPT offers a unique set of simulation capabilities within the current landscape of publicly available planet-disk interaction simulation tools. Its structured interface and underlying technical updates are intended to assist researchers in ongoing explorations of planet formation.

Polarimetric differential imaging with VLT/NACO

Context. The observed diversity of exoplanets can possibly be traced back to the planet formation processes. Planet-disk interactions induce sub-structures in the circumstellar disk that can be revealed via scattered light observations. However, a high-contrast imaging technique such as polarimetric differential imaging (PDI) must first be applied to suppress the stellar diffraction halo. Aims. In this work we present the PDI PiPelIne for NACO data (PIPPIN), which reduces the archival polarimetric observations made with the NACO instrument at the Very Large Telescope. Prior to this work, such a comprehensive pipeline to reduce polarimetric NACO data did not exist. We identify a total of 243 datasets of 57 potentially young stellar objects observed before NACO’s decommissioning. Methods. The PIPPIN pipeline applies various levels of instrumental polarisation correction and is capable of reducing multiple observing setups, including half-wave plate or de-rotator usage and wire-grid observations. A novel template-matching method is applied to assess the detection significance of polarised signals in the reduced data. Results. In 22 of the 57 observed targets, we detect polarised light resulting from a scattering of circumstellar dust. The detections exhibit a collection of known sub-structures, including rings, gaps, spirals, shadows, and in- or outflows of material. Since NACO was equipped with a near-infrared wavefront sensor, it made unique polarimetric observations of a number of embedded protostars. This is the first time detections of the Class I objects Elia 2-21 and YLW 16A have been published. Alongside the outlined PIPPIN pipeline, we publish an archive of the reduced data products, thereby improving the accessibility of these data for future studies.

The PFS View of TOI-677 b: A Spin–Orbit Aligned Warm Jupiter in a Dynamically Hot System* *This paper includes data gathered with the 6.5 meter Magellan Telescopes located at Las Campanas Observatory, Chile.

TOI-677 b is part of an emerging class of “tidally detached” gas giants (a/R ⋆ ≳ 11) that exhibit large orbital eccentricities and yet low stellar obliquities. Such sources pose a challenge for models of giant planet formation, which must account for the excitation of high eccentricities without large changes in the orbital inclination. In this work, we present a new Rossiter–McLaughlin measurement of the tidally detached warm Jupiter TOI-677 b, obtained using high-precision radial velocity observations with Magellan’s Planet Finder Spectrograph (PFS). Combined with previously published observations from the Very Large Telescope’s ESPRESSO spectrograph, we derive one of the most precisely constrained sky-projected spin–orbit angle measurements to date for an exoplanet. The combined fit offers a refined set of self-consistent parameters, including a low sky-projected stellar obliquity of λ=3.°2−1.∘5+1.∘6 and a moderately high eccentricity of e=0.460−0.018+0.019 , which further constrain the puzzling architecture of this system. We examine several potential scenarios that may have produced the current TOI-677 orbital configuration, ultimately concluding that TOI-677 b most likely had its eccentricity excited through disk–planet interactions. This system adds to a growing population of aligned warm Jupiters on eccentric orbits around hot (T eff > 6100 K) stars.