Earth and Planetary Astrophysics (astro-ph.EP)

Abstract: In the conventional theory of planet formation, it is assumed that protoplanetary disks are axisymmetric and have a smooth radial profile. However, recent radio observations of protoplanetary disks have revealed that many of them have complex radial structures. In this study, we perform a series of $N$-body simulations to investigate how planets are formed in protoplanetary disks with radial structures. For this purpose, we consider the effect of continuous pebble accretion onto the discontinuity boundary within the terrestrial planet-forming region ($\sim0.6$ AU). We found that protoplanets grow efficiently at the discontinuity boundary, reaching the Earth mass within $\sim10^4$ years. We confirmed that giant collisions of protoplanets occur universally in our model. Moreover, we found that multiple planet-sized bodies form at regular intervals in the vicinity of the discontinuity boundary. These results indicate the possibility of the formation of solar system-like planetary systems in radially structured protoplanetary disks.

Abstract: We study the response of hot Jupiters to a static tidal perturbation using the Concentric MacLaurin Spheroid (CMS) method. For strongly irradiated planets, we first performed radiative transfer calculations to relate the planet's equilibrium temperature, T_eq, to its interior entropy. We then determined the gravity harmonics, shape, moment of inertia, and the static Love numbers for a range of two-layer interior models that assume a rocky core plus a homogeneous and isentropic envelope composed of hydrogen, helium, and heavier elements. We identify general trends and then study HAT-P-13b, the WASP planets 4b, 12b, 18b, 103b, and 121b, as well as Kepler-75b and CoRot-3b. We compute the Love numbers, k_nm, and transit radius correction, Delta R, which we compare with predictions in the literature. We find that the Love number, k_22, of tidally locked giant planets cannot exceed the value 0.6, and that the high T_eq consistent with strongly irradiated hot Jupiters tend %lead to further lower k_22. While most tidally locked planets are well described by a linear-regime response of k_22 = 3 J_2/q_0 (where q_0 is the rotation parameter of the gravitational potential), for extreme cases such as WASP-12b, WASP-103b and WASP-121b, nonlinear effects can account for over 10% of the predicted k_22. k_22 values larger than 0.6, as they have been reported for planets WASP-4b and HAT-P13B, cannot result from a static tidal response without extremely rapid rotation, and thus are inconsistent with their expected tidally-locked state.

Abstract: Far ultraviolet (FUV) emission lines from dwarf stars are important driving sources of photochemistry in planetary atmospheres. Properly interpreting spectral features of planetary atmospheres critically depends on the emission of its host star. While the spectral energy distributions (SEDs) of K- and M-type stars have been extensively characterized by previous observational programs, the full X-ray to infrared SED of F-type stars has not been assembled to support atmospheric modeling. On the second flight of the Suborbital Imaging Spectrograph for Transition-region Irradiance from Nearby Exoplanet host stars (SISTINE-2) rocket-borne spectrograph, we successfully captured the FUV spectrum of Procyon A (F5 IV-V) and made the first simultaneous observation of several emission features across the FUV bandpass (1010 - 1270 and 1300 - 1565 \r{A}) of any cool star. We combine flight data with stellar models and archival observations to develop the first SED of a mid-F star. We model the response of a modern Earth-like exoplanet's upper atmosphere to the heightened X-ray and extreme ultraviolet radiation within the habitable zone of Procyon A. These models indicate that this planet would not experience significant atmospheric escape. We simulate observations of the Ly$\alpha$ transit signal of this exoplanet with the Hubble Space Telescope (HST) and the Habitable Worlds Observatory (HWO). While marginally detectable with HST, we find that H I Ly$\alpha$ transits of potentially habitable exoplanets orbiting high radial velocity F-type stars could be observed with HWO for targets up to 150 pc away.

Abstract: We study the relationship of zonal gravity coefficients, J_2n, zonal winds, and axial moment of inertia (MoI) by constructing models for the interiors of giant planets. We employ the nonperturbative concentric Maclaurin spheroid (CMS) method to construct both physical (realistic equation of state and barotropes) and abstract (small number of constant-density spheroids) interior models. We find that accurate gravity measurements of Jupiter's and Saturn's J_2, J_4, and J_6 by Juno and Cassini spacecrafts do not uniquely determine the MoI of either planet but do constrain it to better than 1%. Zonal winds (or differential rotation, DR) then emerge as the leading source of uncertainty. For Saturn, they are predicted to decrease the MoI by 0.4% because they reach a depth of ~9000 km while on Jupiter, they appear to reach only ~3000 km. We thus predict DR to affect Jupiter's MoI by only 0.01%, too small by one order of magnitude to be detectable by the Juno spacecraft. We find winds primarily affect the MoI indirectly via the gravity harmonic J_6 while direct contributions are much smaller because the effects of pro- and retrograde winds cancel. DR contributes +6% and -0.8% to Saturn's and Jupiter's J_6 value, respectively. This changes the J_6 contribution that comes from the uniformly rotating bulk of the planet that correlates most strongly with the predicted MoI. With our physical models, we predict Jupiter's MoI to be 0.26393+-0.00001. For Saturn, we predict 0.2181+-0.0002, assuming a rotation period of 10:33:34 h that matches the observed polar radius.

Abstract: We construct models for Jupiter's interior that match the gravity data obtained by the Juno and Galileo spacecrafts. To generate ensembles of models, we introduce a novel quadratic Monte Carlo technique that is more efficient in confining fitness landscapes than affine invariant method that relies on linear stretch moves. We compare how long it takes the ensembles of walkers in both methods to travel to the most relevant parameter region. Once there, we compare the autocorrelation time and error bars of the two methods. For a ring potential and the 2d Rosenbrock function, we find that our quadratic Monte Carlo technique is significantly more efficient. Furthermore we modified the walk moves by adding a scaling factor. We provide the source code and examples so that this method can be applied elsewhere. Here we employ our method to generate five-layer models for Jupiter's interior that include winds and a prominent dilute core, which allows us to match the planet's even and odd gravity harmonics. We compare predictions from the different model ensembles and analyze how much an increase of the temperature at 1 bar and ad hoc change to the equation of state affects the inferred amount of heavy elements in atmosphere and in the planet overall.

Abstract: We investigate the evolution of deuterium to hydrogen mass ratio (D/H) driven by EUV photoevaporation of hydrogen-rich atmospheres of close-in sub-Neptunes around solar-type stars. For the first time, the diffusion-limited approach in conjunction with energy-limited photoevaporation is considered in evaluating deuterium escape from evolving exoplanet H/He envelopes. We find that the planets with smaller initial gas envelopes and thus smaller sizes can lead to weaker atmospheric escape, which facilitates hydrogen-deuterium fractionation. Specifically, in our grid of simulations with a low envelope mass fraction less than 0.005, a low-mass sub-Neptune (4-$5M_\oplus$) at about 0.25-0.4 au or a high-mass sub-Neptune (10-$15M_\oplus$) at about 0.1-0.25 au can increase the D/H values by greater than 20% over 7.5 Gyrs. Akin to the helium-enhanced envelopes of sub-Neptunes due to photoevaporating escape, the planets along the upper boundary of the radius valley are the best targets to detect high D/H ratios. The ratio can rise by a factor of $\lesssim$ 1.65 within 7.5 Gyrs in our grid of evolutionary calculations. D/H is expected to be higher in thinner envelopes so long as the planets do not become bare rocky cores.

Abstract: When observing transmission spectra produced by atmospheres of ultra-hot Jupiters, large telescopes are typically the instrument of choice due to the very weak signal of the planet's atmosphere. This study aims to alleviate the desire for large telescopes by illustrating that the same science is possible with smaller telescope classes. We use the cross-correlation technique to showcase the potential of the high-resolution spectrograph FOCES at Wendelstein Observatory and demonstrate its potential to resolve the atmosphere of the ultra-hot Jupiter, KELT-9 b. A performance comparison is conducted between FOCES and HARPS-N spectrographs, considering both single transit and combined observations over three nights. With FOCES, we have detected seven species in KELT-9 b's atmosphere: Ti II, Fe I, Fe II, Na I, Mg I, Na II, Cr II, Sc II. Although HARPS-N surpasses FOCES in performance, our results reveal that smaller telescope classes are capable of resolving ultra-hot Jupiter atmospheres. This broadens the scope of potential studies, allowing for investigations into phenomena like temporal variations in atmospheric signals and the atmospheric loss characteristics of these close-in planets.