No Arabic abstract
The validity of the widely used linear mixing approximation for the equations of state (EOS) of planetary ices is investigated at pressure-temperature conditions typical for the interior of Uranus and Neptune. The basis of this study are ab initio data ranging up to 1000 GPa and 20 000 K calculated via density functional theory molecular dynamics simulations. In particular, we calculate a new EOS for methane and EOS data for the 1:1 binary mixtures of methane, ammonia, and water, as well as their 2:1:4 ternary mixture. Additionally, the self-diffusion coefficients in the ternary mixture are calculated along three different Uranus interior profiles and compared to the values of the pure compounds. We find that deviations of the linear mixing approximation from the results of the real mixture are generally small; for the thermal EOS they amount to 4% or less. The diffusion coefficients in the mixture agree with those of the pure compounds within 20% or better. Finally, a new adiabatic model of Uranus with an inner layer of almost pure ices is developed. The model is consistent with the gravity field data and results in a rather cold interior ($mathrm{T_{core}} mathtt{sim}$ 4000 K).
We discuss in a thermodynamic, geologically empirical way the long-term nature of the stable majority ices that could be present in Kuiper Belt Object 2014 MU69 after its 4.6 Gyr residence in the EKB as a cold classical object. Considering the stability versus sublimation into vacuum for the suite of ices commonly found on comets, Centaurs, and KBOs at the average ~40K sunlit surface temperature of MU69 over Myr to Gyr, we find only 3 common ices that are truly refractory: HCN, CH3OH, and H2O (in order of increasing stability). NH3 and H2CO ices are marginally stable and may be removed by any positive temperature excursions in the EKB, as produced every 1e8 - 1e9 yrs by nearby supernovae and passing O/B stars. To date the NH team has reported the presence of abundant CH3OH and evidence for H2O on MU69s surface (Lisse et al. 2017, Grundy et al. 2020). NH3 has been searched for, but not found. We predict that future absorption feature detections will be due to an HCN or poly-H2CO based species. Consideration of the conditions present in the EKB region during the formation era of MU69 lead us to infer that it formed in the dark, in an optically thick mid-plane, unable to see the nascent, variable, highly luminous Young Stellar Object-TTauri Sun, and that KBOs contain HCN and CH3OH ice phases in addition to the H2O ice phases found in their Short Period comet descendants. Finally, when we apply our ice thermal stability analysis to bodies/populations related to MU69, we find that methanol ice may be ubiquitous in the outer solar system; that if Pluto is not a fully differentiated body, then it must have gained its hypervolatile ices from proto-planetary disk sources in the first few Myr of the solar systems existence; and that hypervolatile rich, highly primordial comet C/2016 R2 was placed onto an Oort Cloud orbit on a similar timescale.
Context. The study of linear waves and instabilities is necessary to understand the physical evolution of an atmosphere, and can provide physical interpretation of the complex flows found in simulations performed using Global Circulation Models (GCM). In particular, the acceleration of superrotating flow at the equator of hot Jupiters has mostly been studied under several simplifying assumptions, the relaxing of which may impact final results. Aims. We develop and benchmark a publicly available algorithm to identify the eigenmodes of an atmosphere around any initial steady state. We also solve for linear steady states. Methods. We linearise the hydrodynamical equations of a planetary atmosphere in a steady state with arbitrary velocities and thermal profile. We then discretise the linearised equations on an appropriate staggered grid, and solve for eigenvectors and linear steady solutions with the use of a parallel library for linear algebra: ScaLAPACK. We also implement a posteriori calculation of an energy equation in order to obtain more information on the underlying physics of the mode. Results. Our code is benchmarked against classical wave and instability test cases in multiple geometries. The steady linear circulation calculations also reproduce expected results for the atmosphere of hot Jupiters. We finally show the robustness of our energy equation, and its power to obtain physical insight into the modes. Conclusions. We have developed and benchmarked a code for the study of linear processes in planetary atmospheres, with an arbitrary steady state. The calculation of an a posteriori energy equation provides both increased robustness and physical meaning to the obtained eigenmodes. This code can be applied to various problems, and notably to further study the initial spin up of superrotation of GCM simulations of hot Jupiter.
Molecular oxygen, nitrogen, and ozone have been detected in the Solar System. They are also expected to be present in ice-grain mantles within star-forming regions. Laboratory experiments that simulate energetic processing (ions, photons, and electrons) of ices are essential for interpreting and directing future astronomical observations. We provide VUV photoabsorption spectroscopic data of energetically processed nitrogen- and oxygen-rich ices that will help to identify absorption bands and/or spectral slopes observed on icy objects in the Solar System and on ice-grain mantles of the interstellar medium. We present VUV photoabsorption spectra of frozen O2 and N2, a 1:1 mixture of both, and a new systematic set of pure and mixed nitrogen oxide ices. Spectra were obtained at 22 K before and after 1 keV electron bombardment of the ice sample. Ices were then annealed to higher temperatures to study their thermal evolution. In addition, Fourier-transform infrared spectroscopy was used as a secondary probe of molecular synthesis to better identify the physical and chemical processes at play. Our VUV data show that ozone and the azide radical (N3) are observed in our experiments after electron irradiation of pure O2 and N2 ices, respectively. Energetic processing of an O2:N2 = 1:1 ice mixture leads to the formation of ozone along with a series of nitrogen oxides. The electron irradiation of solid nitrogen oxides, pure and in mixtures, induces the formation of new species such as O2, N2 , and other nitrogen oxides not present in the initial ice. Results are discussed here in light of their relevance to various astrophysical environments. Finally, we show that VUV spectra of solid NO2 and water can reproduce the observational VUV profile of the cold surface of Enceladus, Dione, and Rhea, strongly suggesting the presence of nitrogen oxides on the surface of the icy Saturn moons.
Oxygen atom addition and insertion reactions may provide a pathway to chemical complexity in ices that are too cold for radicals to diffuse and react. We have studied the ice-phase reactions of photo-produced oxygen atoms with C2 hydrocarbons under ISM-like conditions. The main products of oxygen atom reactions with ethane are ethanol and acetaldehyde; with ethylene are ethylene oxide and acetaldehyde; and with acetylene is ketene. The derived branching ratio from ethane to ethanol is ~0.74 and from ethylene to ethylene oxide is ~0.47. For all three hydrocarbons there is evidence of an effectively barrierless reaction with O(^1D) to form oxygen-bearing organic products; in the case of ethylene, there may be an additional barriered contribution of the ground-state O(^3P) atom. Thus, oxygen atom reactions with saturated and unsaturated hydrocarbons are a promising pathway to chemical complexity even at very low temperatures where the diffusion of radical species is thermally inaccessible.
Ices are an important constituent of protoplanetary disks. New observational facilities, notably JWST, will greatly enhance our view of disk ices by measuring their infrared spectral features. We present a suite of models to complement these upcoming observations. Our models use a kinetics-based gas-grain chemical evolution code to simulate the distribution of ices in a disk, followed by a radiative transfer code using a subset of key ice species to simulate the observations. We present models reflecting both molecular inheritance and chemical reset initial conditions. We find that H$_2$O, CO$_2$, and CH$_3$OH near-to-mid-IR absorption features are readily observable in disk-integrated spectra of highly-inclined disks while CO, NH$_3$, and CH$_4$ ice do not show prominent features. CH$_3$OH ice has low abundance and is not observable in the reset model, making this species an excellent diagnostic of initial chemical conditions. CO$_2$ ice features exhibit the greatest change over disk lifetime: decreasing and increasing for the inheritance and reset models, respectively. Spatially-resolved spectra of edge-on disks, possible with JWSTs integral field unit observing modes, are ideal for constraining the vertical distribution of ices and may be able to isolate features from ices closer to the midplane (e.g., CO) given sufficient sensitivity. Spatially-resolved spectra of face-on disks can trace scattered-light features from H$_2$O, CO$_2$, and CH$_3$OH, plus CO and CH$_4$ from the outermost regions. We additionally simulate far-IR H$_2$O ice emission features and find they are strongest for disks viewed face-on.