Earth and Planetary Astrophysics (astro-ph.EP)

Abstract: Context: Polycyclic Aromatic Hydrocarbons, largely known as PAHs, are widespread in the universe and have been identified in a vast array of astronomical observations from the interstellar medium to protoplanetary discs. They are likely to be associated with the chemical history of the universe and the emergence of life on Earth. However, their abundance on exoplanets remains unknown. Aims: We aim to investigate the feasibility of PAH formation in the thermalized atmospheres of irradiated and non-irradiated hot Jupiters around Sun-like stars. Methods: To this aim, we introduced PAHs in the 1-D self-consistent forward modeling code petitCODE. We simulated a large number of planet atmospheres with different parameters (e.g. carbon to oxygen ratio, metallicity, and effective planetary temperature) to study PAH formation. By coupling the thermochemical equilibrium solution from petitCODE with the 1-D radiative transfer code, petitRADTRANS, we calculated the synthetic transmission and emission spectra for irradiated and non-irradiated planets, respectively, and explored the role of PAHs on planet spectra. Results: Our models show strong correlations between PAH abundance and the aforementioned parameters. In thermochemical equilibrium scenarios, an optimal temperature, elevated carbon to oxygen ratio, and increased metallicity values are conducive to the formation of PAHs, with the carbon to oxygen ratio having the largest effect.

Abstract: This paper is written as a follow-up observations to reinterpret the radial velocity (RV) of HD 36384, where the existence of planetary systems is known to be ambiguous. In giants, it is, in general, difficult to distinguish the signals of planetary companions from those of stellar activities. Thus, known exoplanetary giant hosts are relatively rare. We, for many years, have obtained RV data in evolved stars using the high-resolution, fiber-fed Bohyunsan Observatory Echelle Spectrograph (BOES) at the Bohyunsan Optical Astronomy Observatory (BOAO). Here, we report the results of RV variations in the M giant HD 36384. We have found two significant periods of 586d and 490d. Considering the orbital stability, it is impossible to have two planets at so close orbits. To determine the nature of the RV variability variations, we analyze the HIPPARCOS photometric data, some indicators of stellar activities, and line profiles. A significant period of 580d was revealed in the HIPPARCOS photometry. H{\alpha} EW variations also show a meaningful period of 582d. Thus, the period of 586d may be closely related to the rotational modulations and/or stellar pulsations. On the other hand, the other significant period of 490d is interpreted as the result of the orbiting companion. Our orbital fit suggests that the companion was a planetary mass of 6.6 MJ and is located at 1.3 AU from the host.

Abstract: The interior of Enceladus, a medium sized icy moon of Saturn hosts hydrothermal activity and exhibits tidal heating and related geyser-like activity. There are major disagreements in the existing literature on the porosity of the interior, due to the different theoretical assumptions on which porosity related calculations were based. We present an application of experimental equations - derived for Earth - for icy planetary objects and Enceladus in particular. We chose a set of boundary values for our initial parameters from measured porosity values of chondrite samples as references, and calculated the porosity related values of Enceladus using various approaches. We present a comprehensive investigation of the effects of using these different porosity calculation methods on icy moons. With our most realistic approach we also calculated the same values for Earth and Mars for comparison. Our result for Enceladus is a minimum porosity of about 5\% at the centre of the body. For the total pore volume we estimated $1.51*10^7 km^3$ for Enceladus, $2.11*10^8 km^3$ for Earth and $1.62*10^8 km^3$ for Mars. Using the same method, we estimated the total pore surface area. From this we derived that the pore surface under a given $1 km^2$ area of the surface on Enceladus is about $1.37*10^9 km^2$, while for Earth this value is only $5.07*10^7 km^2$.

Abstract: Context. The long-term evolution of an atmosphere and the remote detectability of its chemical constituents are susceptible to how the atmospheric gas responds to stellar irradiation. The response remains poorly characterized for water and its dissociation products, however, this knowledge is relevant to our understanding of hypothetical water-rich exoplanets. Aims: Our work investigates the effect of photoelectrons, namely, the non-thermal electrons produced by photoionizing stellar radiation on the heating and ionization of extended atmospheres dominated by the dissociation products of water. Methods: We used a Monte Carlo model and up-to-date collision cross sections to simulate the slowing down of photoelectrons in O-H mixtures for a range of fractional ionizations and photoelectron energies. Results: We find that that the fraction of energy of a photoelectron that goes into heating is similar in a pure H gas and in O-H mixtures, except for very low fractional ionizations, whereby the O atom remains an efficient sink of energy. The O-H mixtures will go on to produce more electrons because the O atom is particularly susceptible to ionization. We quantified all that information and present it in a way that can be easily incorporated into photochemical-hydrodynamical models. Conclusions: Neglecting the role of photoelectrons in models of water-rich atmospheres will result in overestimations of the atmospheric heating and, foreseeably, the mass-loss rates as well. It will also underestimate the rate at which the atmospheric gas becomes ionized, which may have implications for the detection of extended atmospheres with Lyman-{\alpha} transmission spectroscopy. Our simulations for the small exoplanets {\pi} Men c and TRAPPIST-1 b reveal that they respond very differently to irradiation from their host stars, with water remaining in molecular form at lower pressures in the latter case.

Abstract: The enigmatic nature of 55 Cancri e has defied theoretical explanation. Any explanation needs to account for the observed variability of its secondary eclipse depth, which is at times consistent with zero in the visible/optical range of wavelengths -- a phenomenon that does not occur with its also variable infrared eclipses. Yet, despite this variability its transit depth remains somewhat constant in time and is inconsistent with opaque material filling its Hill sphere. The current study explores the possibility of a thin, transient, secondary atmosphere on 55 Cancri e that is sourced by geochemical outgassing. Its transient nature derives from the inability of outgassing to be balanced by atmospheric escape. As the outgassed atmosphere escapes and is replenished, it rapidly adjusts to radiative equilibrium and the temperature fluctuations cause the infrared eclipse depths to vary. Atmospheres of pure carbon dioxide or carbon monoxide produce sufficient Rayleigh scattering to explain the observed optical/visible eclipse depths, which vanish in the absence of an atmosphere and the presence of a dark rocky surface. Atmospheres of pure methane are ruled out because they produce insufficient Rayleigh scattering. Upcoming observations by the James Webb Space Telescope will potentially allow the atmospheric temperature and surface pressure, as well as the surface temperature, to be measured.

Abstract: How multiple close-in super-Earths form around stars with masses lower than that of the Sun is still an open issue. Several recent modeling studies have focused on planet formation around M-dwarf stars, but so far no studies have focused specifically on K dwarfs, which are of particular interest in the search for extraterrestrial life. We aim to reproduce the currently known population of close-in super-Earths observed around K-dwarf stars and their system characteristics. We performed 48 high-resolution N-body simulations of planet formation via planetesimal accretion using the existing GENGA software running on GPUs. In the simulations we varied the initial disk mass and the solid and gas surface density profiles. Each simulation began with 12000 bodies with radii of between 200 and 2000 km around two different stars, with masses of 0.6 and 0.8 $M_{\odot}$. Most simulations ran for 20 Myr, with several simulations extended to 40 or 100 Myr. The mass distributions for the planets with masses between 2 and 12 $M_\oplus$ show a strong preference for planets with masses $M_p<6$ $M_\oplus$ and a lesser preference for planets with larger masses, whereas the mass distribution for the observed sample increases almost linearly. However, we managed to reproduce the main characteristics and architectures of the known planetary systems and produce mostly long-term angular-momentum-deficit-stable, nonresonant systems, but we require an initial disk mass of 15 $M_\oplus$ or higher and a gas surface density value at 1 AU of 1500 g cm$^{-2}$ or higher. Our simulations also produce many low-mass planets with $M<2$ $M_\oplus$, which are not yet found in the observed population, probably due to the observational biases. The final systems contain only a small number of planets, which could possibly accrete substantial amounts of gas, and these formed after the gas had mostly dissipated.

Abstract: We present new optical transmission spectra for two hot Jupiters: WASP-25b (M = 0.56~M$_J$; R = 1.23 R$_J$; P =~3.76 days) and WASP-124b (M = 0.58~M$_J$; R = 1.34 R$_J$; P = 3.37 days), with wavelength coverages of 4200 - 9100\AA\ and 4570 - 9940\AA, respectively. These spectra are from the ESO Faint Object Spectrograph and Camera (v.2) mounted on the New Technology Telescope (NTT) and Inamori-Magellan Areal Camera & Spectrograph on Magellan Baade. No strong spectral features were found in either spectra, with the data probing 4 and 6 scale heights, respectively. \texttt{Exoretrievals} and \texttt{PLATON} retrievals favor stellar activity for WASP-25b, while the data for WASP-124b did not favor one model over another. For both planets the retrievals found a wide range in the depths where the atmosphere could be optically thick ($\sim0.4\mu$ - 0.2 bars for WASP-25b and 1.6 $\mu$ -- 32 bars for WASP-124b) and recovered a temperature that is consistent with the planets' equilibrium temperatures, but with wide uncertainties (up to $\pm$430$^\circ$K). For WASP-25b, the models also favor stellar spots that are $\sim$500-3000$^\circ$K cooler than the surrounding photosphere. The fairly weak constraints on parameters are owing to the relatively low precision of the data, with an average precision of 840 and 1240 ppm per bin for WASP-25b and WASP-124b, respectively. However, some contribution might still be due to an inherent absence of absorption or scattering in the planets' upper atmospheres, possibly because of aerosols. We attempt to fit the strength of the sodium signals to the aerosol-metallicity trend proposed by McGruder et al. 2023, and find WASP-25b and WASP-124b are consistent with the prediction, though their uncertainties are too large to confidently confirm the trend.