The effect of radiation pressure on the dispersal of photoevaporating
  discs

By: Alfie Robinson, James E. Owen, Richard A. Booth

Observed IR excesses indicate that protoplanetary discs evolve slowly for the majority of their lifetime before losing their near- and mid-IR excesses on short timescales. Photoevaporation models can explain this "two-timescale" nature of disc evolution through the removal of inner regions of discs after a few million years. However, they also predict the existence of a population of non-accreting discs with large cavities. Such discs are s... more
Observed IR excesses indicate that protoplanetary discs evolve slowly for the majority of their lifetime before losing their near- and mid-IR excesses on short timescales. Photoevaporation models can explain this "two-timescale" nature of disc evolution through the removal of inner regions of discs after a few million years. However, they also predict the existence of a population of non-accreting discs with large cavities. Such discs are scarce within the observed population, suggesting the models are incomplete. We explore whether radiation-pressure-driven outflows are able to remove enough dust to fit observations. We simulate these outflows using cuDisc, including dust dynamics, growth/fragmentation, radiative transfer and a parameterisation of internal photoevaporation. We find that, in most cases, dust mass-loss rates are around 5-10 times too small to meet observational constraints. Particles are launched from the disc inner rim, however grains larger than around a micron do not escape in the outflow, meaning mass-loss rates are too low for the initial dust masses at gap-opening. Only systems that have smooth photoevaporation profiles with gas mass-loss rates $>\sim 5 \times 10^{-9}$ $M_\odot$ yr$^{-1}$ and disc dust masses $<\sim$1 $M_\oplus$ at the time of gap opening can meet observational constraints; in the current models these manifest as EUV winds driven by atypically large high-energy photon fluxes. We also find that the height of the disc's photosphere is controlled by small grains in the outflow as opposed to shadowing from a hot inner rim; the effect of this can be seen in synthetic scattered light observations. less
C/O ratios in self-gravitating protoplanetary discs with dust evolution

By: Tamara Molyarova, Eduard Vorobyov, Vitaly Akimkin

Elemental abundances, particularly the C/O ratio, are seen as a way to connect the composition of planetary atmospheres with planet formation scenario and the disc chemical environment. We model the chemical composition of gas and ices in a self-gravitating disc on timescales of 0.5\,Myr since its formation to study the evolution of C/O ratio due to dust dynamics and growth, and phase transitions of the volatile species. We use the thin-dis... more
Elemental abundances, particularly the C/O ratio, are seen as a way to connect the composition of planetary atmospheres with planet formation scenario and the disc chemical environment. We model the chemical composition of gas and ices in a self-gravitating disc on timescales of 0.5\,Myr since its formation to study the evolution of C/O ratio due to dust dynamics and growth, and phase transitions of the volatile species. We use the thin-disc hydrodynamic code FEOSAD, which includes disc self-gravity, thermal balance, dust evolution and turbulent diffusion, and treats dust as a dynamically different and evolving component interacting with the gas. It also describes freeze-out, sublimation and advection of four volatile species: H$_2$O, CO$_2$, CH$_4$ and CO. We demonstrate the effect of gas and dust substructures on the distribution of volatiles and C/O ratios, including the formation of multiple snowlines of one species, and point out the anticorrelation between dust-to-gas ratio and total C/O ratio emerging due to the contribution of oxygen-rich ice mantles. We identify time and spatial locations where two distinct trigger mechanisms for planet formation are operating and differentiate them by C/O ratio range: wide range of the C/O ratios of $0-1.4$ for streaming instability, and a much narrower range $0.3-0.6$ for gravitational instability (with the initial value of 0.34). This conclusion is corroborated by observations, showing that transiting exoplanets, which possibly experienced migration through a variety of disc conditions, have significantly larger spread of C/O in comparison with directly imaged exoplanets likely formed in gravitationally unstable outer disk regions. We show that the ice-phase C/O$\approx0.2-0.3$ between the CO, CO$_2$ and CH$_4$ snowlines corresponds to the composition of the Solar system comets, that represent primordial planetesimals. less
Novel Physics of Escaping Secondary Atmospheres May Shape the Cosmic
  Shoreline

By: Richard D. Chatterjee, Raymond T. Pierrehumbert

Recent James Webb Space Telescope observations of cool, rocky exoplanets reveal a probable lack of thick atmospheres, suggesting prevalent escape of the secondary atmospheres formed after losing primordial hydrogen. Yet, simulations indicate that hydrodynamic escape of secondary atmospheres, composed of nitrogen and carbon dioxide, requires intense fluxes of ionizing radiation (XUV) to overcome the effects of high molecular weight and effic... more
Recent James Webb Space Telescope observations of cool, rocky exoplanets reveal a probable lack of thick atmospheres, suggesting prevalent escape of the secondary atmospheres formed after losing primordial hydrogen. Yet, simulations indicate that hydrodynamic escape of secondary atmospheres, composed of nitrogen and carbon dioxide, requires intense fluxes of ionizing radiation (XUV) to overcome the effects of high molecular weight and efficient line cooling. This transonic outflow of hot, ionized metals (not hydrogen) presents a novel astrophysical regime ripe for exploration. We introduce an analytic framework to determine which planets retain or lose their atmospheres, positioning them on either side of the cosmic shoreline. We model the radial structure of escaping atmospheres as polytropic expansions - power-law relationships between density and temperature driven by local XUV heating. Our approach diagnoses line cooling with a three-level atom model and incorporates how ion-electron interactions reduce mean molecular weight. Crucially, hydrodynamic escape onsets for a threshold XUV flux dependent upon the atmosphere's gravitational binding. Ensuing escape rates either scale linearly with XUV flux when weakly ionized (energy-limited) or are controlled by a collisional-radiative thermostat when strongly ionized. Thus, airlessness is determined by whether the XUV flux surpasses the critical threshold during the star's active periods, accounting for expendable primordial hydrogen and revival by volcanism. We explore atmospheric escape from Young-Sun Mars and Earth, LHS-1140 b and c, and TRAPPIST-1 b. Our modeling characterizes the bottleneck of atmospheric loss on the occurrence of observable Earth-like habitats and offers analytic tools for future studies. less
Using the helium triplet as a tracer of the physics of giant planet
  outflows

By: Matthäus Schulik, James Owen

Hydrodynamic outflows, such as those observed escaping close-in gas giant planets, are not isothermal in structure. Their highly ionized nature allows them to cool adiabatically at distances beyond several planetary radii. The contrast between the hottest gas temperatures at around 10,000K and the coldest at around 1,000K triggers an excess population of the observable helium triplet. This excess is caused by the suppression of collisional ... more
Hydrodynamic outflows, such as those observed escaping close-in gas giant planets, are not isothermal in structure. Their highly ionized nature allows them to cool adiabatically at distances beyond several planetary radii. The contrast between the hottest gas temperatures at around 10,000K and the coldest at around 1,000K triggers an excess population of the observable helium triplet. This excess is caused by the suppression of collisional de-excitation from the triplet state at cool temperatures. Using radiation-hydrodynamic simulations, we show that this helium triplet excess may explain the excess broadening seen in HD 189733b's observed transmission spectrum, demonstrating adiabatic cooling of its outflow, confirming its hydrodynamic nature on scales of several planetary radii. However, further observations are required to confirm this conclusion. Furthermore, we explore a range of electron transitions for neutral helium which were not considered in the previous literature. We find that the He$2^1$S state is unavailable as a potential reservoir for He$2^3$S electrons. Additionally, the de-excitation to the ground state must be considered for stellar spectra later than K2 in predicting the correct helium triplet population. Importantly, since triplet helium inherits momentum from ionized helium as it is generated by recombination, it is significantly less prone to fractionation than ground-state neutral helium. However at separations of $\gtrsim 0.05$~au, ionization at the flow base and drag on helium weaken, leading to significant fractionation of the then mostly neutral helium. This in turn, can cause a suppression of the Helium transit depth, even though the helium line width remains large. less
The First JWST View of a 30-Myr-old Protoplanetary Disk Reveals a
  Late-stage Carbon-rich Phase

By: Feng Long, Ilaria Pascucci, Adrien Houge, Andrea Banzatti, Klaus M. Pontoppidan, Joan Najita, Sebastiaan Krijt, Chengyan Xie, Joe Williams, Gregory J. Herczeg, Sean M. Andrews, Edwin Bergin, Geoffrey A. Blake, María José Colmenares, Daniel Harsono, Carlos E. Romero-Mirza, Rixin Li, Cicero X. Lu, Paola Pinilla, David J. Wilner, Miguel Vioque, Ke Zhang, the JDISCS collaboration

We present a JWST MIRI/MRS spectrum of the inner disk of WISE J044634.16$-$262756.1B (hereafter J0446B), an old ($\sim$34 Myr) M4.5 star but with hints of ongoing accretion. The spectrum is molecule-rich and dominated by hydrocarbons. We detect 14 molecular species (H$_2$, CH$_3$, CH$_4$, C$_2$H$_2$, $^{13}$CCH$_2$, C$_2$H$_4$, C$_2$H$_6$, C$_3$H$_4$, C$_4$H$_2$, C$_6$H$_6$, HCN, HC$_3$N, CO$_2$ and $^{13}$CO$_2$) and 2 atomic lines ([Ne II... more
We present a JWST MIRI/MRS spectrum of the inner disk of WISE J044634.16$-$262756.1B (hereafter J0446B), an old ($\sim$34 Myr) M4.5 star but with hints of ongoing accretion. The spectrum is molecule-rich and dominated by hydrocarbons. We detect 14 molecular species (H$_2$, CH$_3$, CH$_4$, C$_2$H$_2$, $^{13}$CCH$_2$, C$_2$H$_4$, C$_2$H$_6$, C$_3$H$_4$, C$_4$H$_2$, C$_6$H$_6$, HCN, HC$_3$N, CO$_2$ and $^{13}$CO$_2$) and 2 atomic lines ([Ne II] and [Ar II]), all observed for the first time in a disk at this age. The detection of spatially unresolved H$_2$ and Ne gas strongly supports that J0446B hosts a long-lived primordial disk, rather than a debris disk. The marginal H$_2$O detection and the high C$_2$H$_2$/CO$_2$ column density ratio indicate that the inner disk of J0446B has a very carbon-rich chemistry, with a gas-phase C/O ratio $\gtrsim$2, consistent with what have been found in most primordial disks around similarly low-mass stars. In the absence of significant outer disk dust substructures, inner disks are expected to first become water-rich due to the rapid inward drift of icy pebbles, and evolve into carbon-rich as outer disk gas flows inward on longer timescales. The faint millimeter emission in such low-mass star disks implies that they may have depleted their outer icy pebble reservoir early and already passed the water-rich phase. Models with pebble drift and volatile transport suggest that maintaining a carbon-rich chemistry for tens of Myr likely requires a slowly evolving disk with $\alpha-$viscosity $\lesssim10^{-4}$. This study represents the first detailed characterization of disk gas at $\sim$30 Myr, strongly motivating further studies into the final stages of disk evolution. less
Predictions of Dust Continuum Emission from a Potential Circumplanetary
  Disk: A Case Study of the Planet Candidate AB Aurigae b

By: Yuhito Shibaike, Jun Hashimoto, Ruobing Dong, Christoph Mordasini, Misato Fukagawa, Takayuki Muto

Gas accreting planets embedded in protoplanetary disks are expected to show dust thermal emission from their circumplanetary disks (CPDs). However, a recently reported gas accreting planet candidate, AB Aurigae b, has not been detected in (sub)millimeter continuum observations. We calculate the evolution of dust in the potential CPD of AB Aurigae b and predict its thermal emission at 1.3 mm wavelength as a case study, where the obtained fea... more
Gas accreting planets embedded in protoplanetary disks are expected to show dust thermal emission from their circumplanetary disks (CPDs). However, a recently reported gas accreting planet candidate, AB Aurigae b, has not been detected in (sub)millimeter continuum observations. We calculate the evolution of dust in the potential CPD of AB Aurigae b and predict its thermal emission at 1.3 mm wavelength as a case study, where the obtained features may also be applied to other gas accreting planets. We find that the expected flux density from the CPD is lower than the 3-sigma level of the previous continuum observation by ALMA with broad ranges of parameters, consistent with the non-detection. However, the expected planet mass and gas accretion rate are higher if the reduction of the observed near-infrared continuum and H-alpha line emission due to the extinction by small grains is considered, resulting in higher flux density of the dust emission from the CPD at (sub)millimeter wavelength. We find that the corrected predictions of the dust emission are stronger than the 3-sigma level of the previous observation with the typical dust-to-gas mass ratio of the inflow to the CPD. This result suggests that the dust supply to the vicinity of AB Aurigae b is small if the planet candidate is not the scattered light of the star but is a planet and has a CPD. Future continuum observations at shorter wavelength are preferable to obtain more robust clues to the question whether the candidate is a planet or not. less
Electron-Induced Radiolysis of Water Ice and the Buildup of Oxygen

By: Chantal Tinner, André Galli, Fiona Bär, Antoine Pommerol, Martin Rubin, Audrey Vorburger, Peter Wurz

Irradiation by energetic ions, electrons, and UV photons induces sputtering and chemical processes (radiolysis) in the surfaces of icy moons, comets, and icy grains. Laboratory experiments, both of ideal surfaces and of more complex and realistic analog samples, are crucial to understand the interaction of surfaces of icy moons and comets with their space environment. This study shows the first results of mass spectrometry measurements from... more
Irradiation by energetic ions, electrons, and UV photons induces sputtering and chemical processes (radiolysis) in the surfaces of icy moons, comets, and icy grains. Laboratory experiments, both of ideal surfaces and of more complex and realistic analog samples, are crucial to understand the interaction of surfaces of icy moons and comets with their space environment. This study shows the first results of mass spectrometry measurements from porous water ice regolith samples irradiated with electrons as a representative analogy to water-ice rich surfaces in the solar system. Previous studies have shown that most electron-induced H2O radiolysis products leave the ice as H2 and O2 and that O2 can be trapped under certain conditions in the irradiated ice. Our new laboratory experiments confirm these findings. Moreover, they quantify residence times and saturation levels of O2 in originally pure water ice. H2O may also be released from the water ice by irradiation, but the quantification of the released H2O is more difficult and the total amount is sensitive to the electron flux and energy. less
Insights on the Rotational State and Shape of Asteroid (203) Pompeja
  from TESS Photometry

By: Oriel A. Humes, Josef Hanuš

The Main Belt asteroid (203) Pompeja shows evidence of extreme variability in visible and near-infrared spectral slope with time. The observed spectral variability has been hypothesized to be attributed to spatial variations across Pompeja's surface. In this scenario, the observed spectrum of Pompeja is dependent on the geometry of the Sun and the observer relative to the asteroid's spin pole and surface features. Knowledge of the rotationa... more
The Main Belt asteroid (203) Pompeja shows evidence of extreme variability in visible and near-infrared spectral slope with time. The observed spectral variability has been hypothesized to be attributed to spatial variations across Pompeja's surface. In this scenario, the observed spectrum of Pompeja is dependent on the geometry of the Sun and the observer relative to the asteroid's spin pole and surface features. Knowledge of the rotational spin pole and shape can be gleaned from light curves and photometric measurements. However, dense light curves of Pompeja are only available from two apparitions. Further, previous estimates of Pompeja's sidereal period are close to being Earth-commensurate, making ground-based light curves difficult to obtain. To overcome these difficulties, we implement a pipeline to extract a dense light curve of Pompeja from cutouts of TESS Full Frame Images. We succeeded in obtaining a dense light curve of Pompeja covering $\sim$22 complete rotations. We measure a synodic period of $P_{syn} =24.092 \pm 0.005$ hours and amplitude of 0.073 $\pm$ 0.002 magnitudes during Pompeja's 2021 apparition in the TESS field of view. We use this light curve to refine models of Pompeja's shape and spin pole orientation, yielding two spin pole solutions with sidereal periods and spin pole ecliptic coordinates of $P_{\mathrm{sid}, 1} = 24.0485 \pm 0.0001$ hours, $\lambda_1 = 132^{\circ}$, $\beta_1 = +41^{\circ}$ and $P_{\mathrm{sid}, 2} = 24.0484 \pm 0.0001$ hours, $\lambda_2 =307^{\circ}$, $\beta_2 =+34^{\circ}$. Finally, we discuss the implications of the derived shape and spin models for spectral variability on Pompeja. less
Deep imaging of three accelerating stars using SHARK-NIR and LMIRCam at
  LBT

By: D. Mesa, R. Gratton, V. D'Orazi, E. Carolo, D. Vassallo, J. Farinato, L. Marafatto, K. Wagner, J. Hom, S. Ertel, Th. Henning, C. Desgrange, D. Barbato, M. Bergomi, P. Cerpelloni, S. Desidera, S. Di Filippo, D. Doelman, T. S. Gomes Machado, D. Greggio, P. Grenz, M. Kenworthy, F. Laudisio, C. Lazzoni, J. Leisenring, A. Lorenzetto, K. K. R. Santhakumari, D. Ricci, F. Snik, G. Umbriaco, M. C. Vega Pallauta, V. Viotto

The combination of detection techniques enhances our ability to identify companions orbiting nearby stars. We employed high-contrast imaging to constrain mass and separation of possible companions responsible for the significant proper motion anomalies of the nearby stars HIP 11696, HIP 47110 and HIP 36277. These targets were observed using the LBT's high-contrast camera, SHARK-NIR, in H-band using a Gaussian coronagraph, and with the LMIRC... more
The combination of detection techniques enhances our ability to identify companions orbiting nearby stars. We employed high-contrast imaging to constrain mass and separation of possible companions responsible for the significant proper motion anomalies of the nearby stars HIP 11696, HIP 47110 and HIP 36277. These targets were observed using the LBT's high-contrast camera, SHARK-NIR, in H-band using a Gaussian coronagraph, and with the LMIRCam instrument in the L'-band and using a vAPP coronagraph. Both observations were conducted simultaneously. Additionally, constraints at short separations from the host star are derived analyzing the renormalized unit weight error (RUWE) values from the Gaia catalogue. We find that the companion responsible for the anomaly signal of HIP 11696 is likely positioned at a distance from 2.5 to 28 astronomical units from its host. Its mass is estimated to be between 4 and 16 Jupiter masses, with the greater mass possible only at the upper end of the separation range. Similar limits were obtained for HIP 47110 where the companion should reside between 3 and 30 au with a mass between 3 and 10 MJup. For HIP 36277, we identified a faint stellar companion at large separation, though it might be substellar depending on the assumed age for the star. Considering the older age, this object accounts for the absolute value of the PMa vector but not for its direction. Additionally, we found a substellar candidate companion at a closer separation that could explain the PMa signal, considering a younger age for the system. less
Ionospheric conductances at the giant planets of the Solar System:a
  comparative study of ionization sources and the impact of meteoric ions

By: Noé Clément, Yuki Nakamura, Michel Blanc, Yuxian Wang, Sariah Al Saati

The dynamics of giant planet magnetospheres is controlled by a complex interplay between their fast rotation, their interaction with the solar wind, and their diverse internal plasma and momentum sources. In the ionosphere, the Hall and Pedersen conductances are two key parameters that regulate the intensity of currents coupling the magnetosphere and the ionosphere, and the rate of angular momentum transfer and power carried by these curren... more
The dynamics of giant planet magnetospheres is controlled by a complex interplay between their fast rotation, their interaction with the solar wind, and their diverse internal plasma and momentum sources. In the ionosphere, the Hall and Pedersen conductances are two key parameters that regulate the intensity of currents coupling the magnetosphere and the ionosphere, and the rate of angular momentum transfer and power carried by these currents. We perform a comparative study of Hall and Pedersen conductivities and conductances in the four giant planets of our Solar System - Jupiter, Saturn, Uranus and Neptune. We use a generic ionospheric model (restraining the studied ions to H3+, CH5+, and meteoric ions) to study the dependence of conductances on the structure and composition of these planets' upper atmospheres and on the main ionization sources (photoionization, ionization by precipitating electrons, and meteoroid ablation). After checking that our model reproduces the conclusions of Nakamura et al. (2022, https://doi.org/10.1029/2022JA030312) at Jupiter, i.e. the contribution of meteoric ions to the height-integrated conductances is non-negligible, we show that this contribution could also be non-negligible at Saturn, Uranus and Neptune, compared with ionization processes caused by precipitating electrons of energies lower than a few keV (typical energies on these planets). However, because of their weaker magnetic field, the conductive layer of these planets is higher than the layer where meteoric ions are mainly produced, limiting their role in magnetosphere-ionosphere coupling. less