[en] We analyze the emission spectrum of the hot Jupiter WASP-12b using our HELIOS-R retrieval code and HELIOS-K opacity calculator. When interpreting Hubble and Spitzer data, the retrieval outcomes are found to be prior-dominated. When the prior distributions of the molecular abundances are assumed to be log-uniform, the volume mixing ratio of HCN is found to be implausibly high. A VULCAN chemical kinetics model of WASP-12b suggests that chemical equilibrium is a reasonable assumption even when atmospheric mixing is implausibly rigorous. Guided by (exo)planet formation theory, we set Gaussian priors on the elemental abundances of carbon, oxygen, and nitrogen with the Gaussian peaks being centered on the measured C/H, O/H, and N/H values of the star. By enforcing chemical equilibrium, we find substellar O/H and stellar to slightly superstellar C/H for the dayside atmosphere of WASP-12b. The superstellar carbon-to-oxygen ratio is just above unity, regardless of whether clouds are included in the retrieval analysis, consistent with Madhusudhan et al. Furthermore, whether a temperature inversion exists in the atmosphere depends on one’s assumption for the Gaussian width of the priors. Our retrieved posterior distributions are consistent with the formation of WASP-12b in a solar-composition protoplanetary disk, beyond the water iceline, via gravitational instability or pebble accretion (without core erosion) and migration inward to its present orbital location via a disk-free mechanism, and are inconsistent with both in situ formation and core accretion with disk migration, as predicted by Madhusudhan et al. We predict that the interpretation of James Webb Space Telescope WASP-12b data will not be prior-dominated.
Journal
Astrophysical Journal Letters; ISSN 2041-8205; 
; v. 847(1); [7 p.]
; v. 847(1); [7 p.]
[en] We develop empirical relationships for the accretion and erosion of colliding gravity-dominated bodies of various compositions under conditions expected in late-stage solar system formation. These are fast, easily coded relationships based on a large database of smoothed particle hydrodynamics (SPH) simulations of collisions between bodies of different compositions, including those that are water rich. The accuracy of these relations is also comparable to the deviations of results between different SPH codes and initial thermal/rotational conditions. We illustrate the paucity of disruptive collisions between major bodies, as compared to collisions between less massive planetesimals in late-stage planet formation, and thus focus on more probable, low-velocity collisions, though our relations remain relevant to disruptive collisions as well. We also pay particular attention to the transition zone between merging collisions and those where the impactor does not merge with the target, but continues downrange, a “hit-and-run” collision. We find that hit-and-run collisions likely occur more often in density-stratified bodies and across a wider range of impact angles than suggested by the most commonly used analytic approximation. We also identify a possible transitional zone in gravity-dominated collisions where larger bodies may undergo more disruptive collisions when the impact velocity exceeds the sound speed, though understanding this transition warrants further study. Our results are contrary to the commonly assumed invariance of total mass (scale), density structure, and material composition on the largest remnants of giant impacts. We provide an algorithm for adopting our model into N-body planet formation simulations, so that the mass of growing planets and debris can be tracked.
Journal