2025-11-17 | | Total: 16
TOI-1422 is a G2 V star ($V = 10.6$ mag) known to host a warm Neptune-sized planet, TOI-1422 b, with a mass and radius of about $9M_{\oplus}$ and $4R_{\oplus}$, on a circular orbit with a period of $12.997$ days. An outer planetary candidate in this system had previously been suggested on the basis of a residual signal in the radial velocity (RV) data with a tentative period of $\sim$29 days, along with a possible single transit-like event, although it was not clear at the time whether the two signals belonged to the same companion. In this work, we confirm the presence of a second transiting planet, TOI-1422 c, a sub-Neptune ($R=2.61\pm0.14 R_{\oplus}$) that orbits with a longer period of 34.563 days. This confirmation is based on the detection of three TESS transits, two from newly available sectors, combined with new and archival RV measurements. The sub-Neptune ($ρ_{\rm c}=4.3^{+1.3}_{-1.0}$ g cm$^{-3}$) is more massive than the inner Neptune ($ρ_{\rm b}=0.93^{+0.21}_{-0.20}$ g cm$^{-3}$), having a mass of $M_{\rm c}=14\pm3 M_{\oplus}$, making TOI-1422 a rare anti-ordered system. Furthermore, we detect transit timing variations (TTVs) on the inner planet, with amplitudes of up to 5 hours, suggesting ongoing dynamical interactions. A dynamical analysis that combined TTVs and RVs indicates that planet c alone is unlikely to account for the full TTV amplitude observed on TOI-1422 b. We investigated whether an additional, as yet undetected companion could account for the observed signal, exploring a range of plausible orbital configurations and finding that a low-mass planet located between the two known orbits may be responsible.
We present pkdgrav3, a high-performance, fully parallel tree-SPH code designed for large-scale hydrodynamic simulations including self-gravity. Building upon the long development history of pkdgrav, the code combines an efficient hierarchical tree algorithm for gravity and neighbor finding with a modern implementation of Smoothed Particle Hydrodynamics (SPH) optimized for massively parallel hybrid CPU/GPU architectures. Its hybrid shared/distributed memory model, combined with an asynchronous communication scheme, allows pkdgrav3 to scale efficiently to thousands of CPU cores and GPUs. We validate the numerical accuracy of pkdgrav3 using a suite of standard tests, demonstrating excellent agreement with analytic or reference solutions. The code was already used in several peer-reviewed publications to model planetary-scale impacts, where SPH's Lagrangian nature allows accurate tracking of material origin and thermodynamic evolution. These examples highlight pkdgrav3's robustness and efficiency in simulating highly dynamical, self-gravitating systems. pkdgrav3 thus provides a powerful, flexible, and scalable platform for astrophysical and planetary applications, capable of exploiting the full potential of modern heterogeneous high-performance computing systems.
Shadows are commonly observed in protoplanetary disks in near-infrared and (sub)millimeter images, often cast by misaligned inner disks or other obscuring material. While recent studies show that shadows can alter disk dynamics, only the case symmetric across the midplane (e.g., from a polar-aligned inner disk) has been studied. Here we study shadows cast by an inner disk with a $30^\circ$ mutual inclination using 3D radiation-hydrodynamical simulations. Given the same shadow shape and amplitude, the $30^\circ$ inclined shadow leads to a much stronger accretion compared with the polar case, reaching $α\sim$ 1, because the disk is squeezed twice in one azimuth, leading to shocks and strong radial flows near the midplane. The outer disk develops a warp: the inner disk region tilts toward alignment with the shadow, while the outer, exponentially tapered disk tilts and twists in a different direction, inclined $\sim$ 32$^\circ$ relative to the inner region. Locally isothermal simulations with a prescribed temperature structure reproduce the effect, confirming that it is thermally driven. Fourier-Hermite analysis shows that it is the m=1, n=1 temperature perturbation that drives the warp by launching bending waves, with the tilting response of the disk approximately proportional to the modal amplitude. This mode always exists unless the shadow is coplanar or polar. Given a fixed temperature contrast, the m=1,n=1 mode peaks at $\sim$15$^\circ$ mutual inclination, but still contributes substantially across 3$^\circ$ to 30$^\circ$. Shadows cause disk warps--they are not only a consequence of them. We discuss testable predictions for current and future ALMA and NIR observations.
Planet formation simulations consistently predict compact systems of numerous small planets in chains of mean motion resonances formed by planet-disk interaction, but transiting planet surveys have found most systems to be non-resonant and somewhat dynamically excited. A scenario in which nearly all of the primordial resonant chains undergo dynamical instabilities and collisions has previously been found to closely match many features of the observed planet sample. However, existing models have not been tested against new observations that show a steep decline in the resonant fraction as a function of stellar age on a timescale of ~100 Myr. We construct a simplified model incorporating Type I migration, growth from embryos, and N-body integrations continued to 500 Myr and use it to generate a synthetic planet population. Nearly all systems exit the disk phase in a resonant configuration but begin slowly diffusing away from the center of the resonance. Dynamical instabilities can arise on timescales of tens or hundreds of Myr, especially when systems formed in disks with a convergent migration trap. In this case, a secondary chain of smaller planets that remained at their birth location eventually breaks, destabilizing the inner resonant chain. We also show that the instability statistics are well modeled by a Weibull distribution, and use this to extrapolate our population to Gyr ages. The close match of our modeled systems to the observed population implies that the high resonance fraction predicted by this class of models is in fact consistent with the data, and the previously-reported overabundance of resonant systems was a consequence of comparing simulations of early evolution to mature Gyr-old systems. This result also suggests that instabilities triggered by disk dissipation or other very early mechanisms are unlikely to be consistent with observed young systems.
Recent observations have revealed the distribution of orbital period ratios of adjacent planets in multiple super-Earth systems and how these distributions change with time. The aim of this study is to clarify under what conditions the observed features of orbital period ratios of super-Earths can be explained, and to identify what causes the dynamical instability of super-Earths captured into resonant chains. We perform N-body simulations for 100 Myr that follow the formation and orbital evolution of super-Earths originating from a ring of planetary embryos at 1 au from the star. The simulations show that super-Earths undergo inward migration in the disk and are captured into mean-motion resonances with their neighbors. As a result, several resonant pairs form a resonant chain. After disk dispersal, some of these chains become dynamically unstable. In such cases, the final distribution of orbital period ratios and their time evolution can be consistent with recent observations. The instabilities of resonant chains are likely triggered by secular perturbations from embryos that remain on outer orbits beyond 1 au, indicating that not only giant planets but also small embryos can disrupt the resonances among inner super-Earths. We therefore further investigate the secular perturbations from outer embryos using analytic formulas and additional orbital calculations. We discuss the conditions required to excite the eccentricities of inner super-Earths on a timescale of about 100 Myr. These conditions include the need for large eccentricities of the outer embryos, as well as constraints on their masses and semimajor axes.
China is planing to launch the Tianwen-4 mission around the year 2030, with its aim being the exploration of Jupiter and its moon, Callisto. Within the realm of deep space exploration, the accuracy of ephemerides is of great importance. Current ephemerides employ a simplified rotation model for Callisto, which this study addresses by proposing a novel dynamical model. This model enhancesthe existing orbital dynamics by integrating Callisto's rotational motions influenced by gravitational torques from the Sun, Jupiter, and other Galilean moons within an inertial frame, capturing the intricate coupling between Callisto's orbital and rotational dynamics. The study establishes a full dynamical model by deriving analytical expressions for this coupling and developing an adjustment model for data fitting using precise orbit determination methods. Furthermore, the influence of tidal effects on Callisto's motion is investigated, considering its multi-layered internal structure. Results demonstrate that the difference between the newly established full model and the model in current ephemerides is on the order of tens of meters. When calculating the impact of different internal structures of Callisto on its orbit, the influence of three-layered and two-layered structures is on the order of meters, suggesting that the development of a high-precision dynamical model requires additional constraints on the internal structure of Callisto. This research provides a novel alternative for a new generation of precise numerical ephemerides for Callisto. Additionally, these findings provide a testing platform for the data from the Tianwen-4 mission.
The physical and chemical conditions within a protoplanetary disk play a crucial role in determining its chemical composition, which is subsequently inherited by any forming planets. To probe these conditions, high-resolution molecular line observations, coupled with modelling, are essential. In this study, we investigate the chemistry of the nearby, massive, and relatively line-rich protoplanetary disk around HD 163296 using high-resolution observations from ALMA across Bands 3, 4, 6, and 7. We constrain the disk-averaged and radial distributions of column density and excitation temperature for the detected molecules using the new retrieval code DRive. The disk chemistry is modelled using the astrochemical code PEGASIS, with variations in the initial elemental C/O ratio. Our modelling, informed by molecular observations of HCO+, DCO+, HCN, DCN, CS, HC3N, H2CO, CH3OH, HNCO, and NH2CHO, allows us to place strong constraints on the C/O ratio, with a best-fit value of 1.1 that is broadly consistent with previous estimates. We present the highest-resolution DCO+ emission map of this disk to date, revealing triple-ringed chemical substructures that closely align with the dust continuum rings. Additionally, our results provide the first and most stringent upper limits on the column densities of NH2CHO and HNCO in this protoplanetary disk, measured at < 7e11 cm-2 and < 1e11 cm-2, respectively. Our chemical models suggest that NH2CHO and HNCO predominantly form on grain surfaces within the disk. However, physico-chemical desorption mechanisms are inefficient at releasing these species into detectable gas-phase abundances, yet they remain promising targets for future ALMA observations.
The TRAPPIST-1 planets have become prime targets for studying the atmospheric and geophysical properties of planets around M-dwarf stars. To effectively identify their atmospheric composition, we first must understand their geological evolution. For this study, we focus on enhancing an existing atmosphere-interior exchange model by incorporating additional geological processes relevant to rocky planets. We have extended the model to include the carbon cycle, which enables the model to track four key gas species - CO$_2$, CO, H$_2$O, and H$_2$ - across four planetary reservoirs: the mantle, plate, ocean, and atmosphere. Major features added include surface temperature calculations which are crucial for the carbon cycle, oxygen fugacity as a planetary interior parameter in the model, and oxidation reactions and diffusion-limited escape calculations to the atmosphere portion of the model. We successfully validated the model for Earth and applied this model to study the effect of oxygen fugacity and initial water abundance on TRAPPIST-1 d, e and f. Our results for present-day abundances show that oxygen fugacity significantly affects the partial pressures of H$_2$ and CO$_2$ for all three planets with minor effects for CO on two of the planets. We also found that H$_2$ is strongly dependent on water mass fraction (WMF). The addition of atmospheric processes produced a significant difference in the H$_2$ and CO abundances at present-day. These results highlight the importance of considering interior parameters to be able to further constrain the geological evolution of these planets and effectively put atmosphere observations into context.
We present a comprehensive photometric study of transNeptunian objects (TNOs) by combining data from SDSS, Col-OSSOS, DES, and the recent Rubin First Look (RFL) data. Our database comprises 43 878 measurements in the u, g, r, i, z, and J filters, from which we derived 8 738 phase curves for 1 921 unique objects. From these data, we computed 12 852 absolute color measurements and spectral slope differences for 1 761 objects, allowing a statistical characterization of phase coloring effects. The colors show no strong bimodality or correlation with orbital parameters, emphasizing the importance of phase correction even for small phase angles. The increase in sample size and application of phase corrections fill previously empty regions in color magnitude space likely affected by observational biases, as redder (and thus darker) objects are preferentially lost near detection limits. Notably, our dataset includes the first photometric measurements from Rubin Observatory during RFL, covering eight objects (five newly discovered TNOs and three previously known). These early LSST observations occupy sparsely sampled regions of parameter space, particularly at faint magnitudes, highlighting the discovery and characterization potential of the full survey. We confirm previous results showing that TNO colors vary with phase angle, exhibiting both reddening and bluening trends. Correlations between (dS'/dalpha) and (alpha) strengthen with increasing (Delta lambda), except for Hi - Hz, which tends to neutralize, consistent with the spectral flattening previously reported in visible wavelengths.
Few gas giants have been observed in their accretion phase. ELT instruments will have a higher sensitivity and a smaller inner working angle than instruments up to now, allowing detailed characterisation. We study the observability of accreting gas giants with METIS, the first-generation ELT spectrograph with a resolution R = 1e5. We focus on the accretion-tracing hydrogen recombination lines accessible to METIS, mainly Brackett alpha and Pfund-series lines. Our approach is general but we take PDS 70 b as a fiducial case. To calculate high-resolution line profiles, we combine a semianalytical multidimensional description of the flow onto an accreting planet and its circumplanetary disc (CPD) with local NLTE shock-emission calculations. We assume the limiting scenario of no extinction, appropriate for gas giants in gaps, and negligible contribution from magnetospheric accretion. We use resolved photospheric models, and include detector sensitivities to compute required observing times. Both the planet surface and the CPD surface shocks contribute to the total line profile, which is non-Gaussian and much narrower than the free-fall velocity. When the accretion onto PDS 70 b is at its maximum observed strength, the Br a line peak is equal to the photospheric continuum, which is modulated mostly by water features. However, the rotation of the planet broadens the features, making the shock excess stand out. At Br a, already only the continuum of PDS 70 b should yield SNR = 12 in 4 h, and for the fiducial accretion rate, the line peak will have a per-bin SNR = 28 in 4 h. The peak excess should require only about 10 min to reach SNR = 3. Br a is a potent planet formation tracer accessible to ELT/METIS in little integration time. Resolved line profiles will place independent constraints especially on the mass and radius of accreting planets, and help identify the accretion mechanism(s) at work.
Gravitational coupling between a protoplanetary disc and an embedded planet is often studied in a frame attached to a central star. This frame is non-inertial because of the stellar reflex motion, leading to indirect forces arising in the star-planet-disc system. Here we examine the impact produced by these forces on several aspects of disc-planet coupling using analytical and numerical means. We explore how neglecting indirect forces changes (1) the spatial pattern of the surface density perturbation in the disc, (2) the calculation of the torque exerted on the disc by the planet, and (3) the torque on the planet exerted by the disc. For low-mass planets, in the linear regime, the differences in the perturbation pattern are only in its $m=1$ azimuthal harmonic, with an amplitude increasing with the distance from the star. In this regime both the torque on the planet and the deposition torque density in the disc are only weakly affected by non-inclusion of indirect forces, corroborating some results of studies neglecting indirect forces altogether. For higher mass planets, a broader range of azimuthal harmonics of the perturbation are affected. Also, indirect forces have a stronger effect on the planetary torque and on planet migration in the Type II regime. We highlight the importance of including the planetary indirect force in the calculation of the torque on the disc (if disc evolution accounts for indirect force) to ensure conservation of angular momentum carried by the planet-driven density waves. The corresponding indirect torque has an oscillatory, radially-diverging character.
The discovery of infant (< 50 Myr), close-in (<30-day period) planets is vital in understanding the formation mechanisms that lead to the distribution of mature transiting planets as discovered by Kepler. Despite several discoveries in this age bin, the sample is still too small for a robust statistical comparison to older planets. Here we report the validation of TOI-6448 b, an 8.8 +/- 0.8 Re planet on a 14.8 day orbit. TOI-6448 was previously identified to be a likely member of Vela Population IV. We confirm the star's membership and re-derive the age of the cluster using isochrones, variability, and gyrochronology. We find the star, and thus planet, to be 34 +/- 3 Myr. Like other young planets, TOI-6448 b lands in a region of parameter space with few older planets. While just one data point, this fits with prior findings of an excess of 5-11Re planets around young stars far beyond what can be explained by reduced sensitivity at young ages. Our ongoing search of Vela, Taurus-Auriga, Sco-Cen, and Orion are expected to reveal dozens more < 50 Myr transiting planets.
A widely explored pathway for hot Jupiter (HJ) formation is high-eccentricity migration driven by von Zeipel-Lidov-Kozai cycles induced by an exterior companion. However, for a distant or low-mass companion, this mechanism typically demands that the planet's initial orbit be very nearly perpendicular to that of the companion. In previous work (Yang et al. 2025), we demonstrated that such fine-tuning can be circumvented in the HAT-P-7 system due to the presence of an intermediate body that efficiently couples the orbits of the planet and the distant companion -- a mechanism we termed the eccentricity cascade (EC). In this work, we analytically characterize the dynamics governing the EC and delineate the parameter space within which it effectively operates. Our qualitative results are as follows: (i) The proto-HJ's eccentricity is most efficiently excited when the inner triple is on the verge of dynamical instability, (ii) the addition of a distant fourth body allows this instability to be approached gradually, and (iii) the instability mechanism is closely related to the stability of circumstellar (S-type) planets in binaries. By deriving an analytic criterion for S-type instability, we obtain closed-form expressions describing the onset of the EC. Our results show that efficient HJ formation via the EC occurs across a broad range of intermediate perturbers, highlighting its potential as a robust migration channel.
Stably stratified layers are thought to develop at the top of the liquid metallic cores of many terrestrial planets. We consider the case where the thermal gradient is stable but the compositional gradient is unstable, a situation particularly relevant to Mercury. The strong contrast between molecular diffusivities of temperature and composition leads to fingering convection. We investigate this process using hydrodynamical simulations in a rotating spherical shell, systematically varying the stratification strength N relative to the rotation rate $Ω$. In all regimes, the primary fingering mode forms narrow, elongated structures that shift orientation from the rotation axis to the direction of gravity as $N^2/Ω^2$ exceeds 10. The fingers remain laminar, with transverse scales proportional to thermal stratification but independent of rotation. Fingering convection also drives secondary large-scale flows across most of the explored parameter space, producing diverse dynamics including zonal flows, hemispherical convection, axisymmetric poloidal bands, finger clusters, and toroidal gyres. In the rapidly-rotating regime, laterally inhomogeneous mixing generates zonal flows in thermo-compositional wind balance; zonal flow direction and amplitude depend on $N^2/Ω^2$, with amplitude weakening for strong stratification $N^2/Ω^2>10$. In the intermediate regime ($N^2/Ω^2\sim 1$), axisymmetric or spiraling poloidal bands emerge within the tangent cylinder, gradually overtaking the primary fingers. For stronger stratification, finger clusters and weak, large-scale density anomalies surrounded by toroidal gyres form in the upper domain. These diverse large-scale flows may interact with the dynamo-generated magnetic field in the deeper core, potentially influencing surface magnetic fields.
This PRobe far-Infrared Mission for Astrophysics (PRIMA) mission concept is a proposed mission to NASA's Astrophysics Probe Explorer (APEX) call. The concept features a cryogenically cooled 1.8 m diameter telescope, and is designed to carry two science instruments covering the 24 to 264 $μ$m wavelength range: an imaging polarimeter (PRIMAger) and a spectrometer (FIRESS). The majority of PRIMA's time (75%) will be open to observations proposed by the community (General Observer science / GO), and all of data will be publicly available for archival research (Guest Investigator science / GI). Following up on the successful community engagement created by the first volume of the GO PRIMA Science Book (arXiv:2310.20572), Volume 2 gathers 120 new and updated contributed science cases which could be performed within the context of the PRIMA GO/GI program. This volume reflects the strong development of the community interest, awareness and involvement in PRIMA, and further develops how PRIMA's unprecedented capabilities can be leveraged for an impactful and innovative GO/GI program covering most areas of astrophysics and over 90% of the scientific questions and discovery areas in the Astro2020 decadal survey.
While the influence of supermassive black hole (SMBH) activity on habitability has garnered attention, the specific effects of active galactic nuclei (AGN) winds, particularly ultrafast outflows (UFOs), on planetary atmospheres remain largely unexplored. This study aims to fill this gap by investigating the relationship between SMBH mass at the galactic center and exoplanetary habitability, given that SMBH masses are empirically confirmed to span approximately 5 orders of magnitude in galaxies. Through simplified models, we account for various results involving the relationships between the distance from the planet to the central SMBH and the mass of the SMBH. Specifically, we show that increased SMBH mass leads to higher atmospheric heating and elevated temperatures, greater molecular thermal velocities, and enhanced mass loss, all of which diminish with distance from the galactic center. Energy-driven winds consistently have a stronger impact than momentum-driven ones. Crucially, ozone depletion is shown to rise with SMBH mass and decrease with distance from the galactic center, with nearly complete ozone loss ($\sim100\%$) occurring across galactic scales for SMBHs $\geq 10^8 M_\odot$ in the energy-driven case. This study emphasizes that SMBH growth over cosmic time may have produced markedly different impacts on galactic habitability, depending on both the mass of the central black hole (BH) and the location of planetary systems within their host galaxies.