The Early Instability Scenario: Terrestrial Planet Formation During the Giant Planet Instability, and the Effect of Collisional Fragmentation (Icarus, 2019)

The solar system’s dynamical state can be explained by an orbital instability among the giant planets. A recent model has proposed that the giant planet instability happened during terrestrial planet formation. This scenario has been shown to match the inner solar system by stunting Mars’ growth and preventing planet formation in the asteroid belt. Here we present a large sample of new simulations of the “Early Instability” scenario. We use an N-body integration scheme that accounts for collisional fragmentation, and also perform a large set of control simulations that do not include an early giant planet instability. Since the total particle number decreases slower when collisional fragmentation is accounted for, the growing planets’ orbits are damped more strongly via dynamical friction and encounters with small bodies that dissipate angular momentum (eg: hit-and-run impacts). Compared with simulations without collisional fragmentation, our fully evolved systems provide better matches to the solar system’s terrestrial planets in terms of their compact mass distribution and dynamically cold orbits. Collisional processes also tend to lengthen the dynamical accretion timescales of Earth analogs, and shorten those of Mars analogs. This yields systems with relative growth timescales more consistent with those inferred from isotopic dating. Accounting for fragmentation is thus supremely important for any successful evolutionary model of the inner solar system.

Continue Reading →

How to Characterize the Atmosphere of a Transiting Exoplanet (PASP, 2018)

This tutorial is an introduction to techniques used to characterize the atmospheres of transiting exoplanets. We intend it to be a useful guide for the undergraduate, graduate student, or postdoctoral scholar who wants to begin research in this field, but who has no prior experience with transiting exoplanets. We begin with a discussion of the properties of exoplanetary systems that allow us to measure exoplanetary spectra, and the principles that underlie transit techniques. Subsequently, we discuss the most favorable wavelengths for observing, and explain the specific techniques of secondary eclipses and eclipse mapping, phase curves, transit spectroscopy, and convolution with spectral templates. Our discussion includes factors that affect the data acquisition, and also a separate discussion of how the results are interpreted. Other important topics that we cover include statistical methods to characterize atmospheres such as stacking, and the effects of stellar activity. We conclude by projecting the future utility of large-aperture observatories such as the James Webb Space Telescope and the forthcoming generation of extremely large ground-based telescopes.

Continue Reading →

The Rotationally Modulated Polarization of ? Boo a (MNRAS, 2018)

We have observed the active star ? Boo A (HD 131156A) with high precision broadband linear polarimetry contemporaneously with circular spectropolarimetry. We find both signals are modulated by the 6.43 d rotation period of ? Boo A. The signals from the two techniques are 0.25 out of phase, consistent with the broadband linear polarization resulting from differential saturation of spectral lines in the global transverse magnetic field. The mean magnitude of the linear polarization signal is ?4 ppm?G–1 but its structure is complex and the amplitude of the variations suppressed relative to the longitudinal magnetic field. The result has important implications for current attempts to detect polarized light from hot Jupiters orbiting active stars in the combined light of the star and planet. In such work stellar activity will manifest as noise, both on the time-scale of stellar rotation, and on longer time-scales – where changes in activity level will manifest as a baseline shift between observing runs.

Continue Reading →

Can Moons Have Moons? (MNRAS Letters, 2019)

Each of the giant planets within the Solar system has large moons but none of these moons have their own moons (which we call submoons). By analogy with studies of moons around short-period exoplanets, we investigate the tidal-dynamical stability of submoons. We find that 10?km-scale submoons can only survive around large (1000?km-scale) moons on wide-separation orbits. Tidal dissipation destabilizes the orbits of submoons around moons that are small or too close to their host planet; this is the case for most of the Solar system’s moons. A handful of known moons are, however, capable of hosting long-lived submoons: Saturn’s moons Titan and Iapetus, Jupiter’s moon Callisto, and Earth’s Moon. Based on its inferred mass and orbital separation, the newly discovered exomoon candidate Kepler-1625b-I can in principle host a large submoon, although its stability depends on a number of unknown parameters. We discuss the possible habitability of submoons and the potential for subsubmoons. The existence, or lack thereof, of submoons may yield important constraints on satellite formation and evolution in planetary systems.

Continue Reading →

Nitrogen Isotope Evidence for Anoxic Deep Marine Environments From the Mesoarchean Mosquito Creek Formation, Australia (Precambrian Research, 2019)

Current evidence for oxygenated environments in the Mesoarchean is limited to the shallowest marine and fluvio-lacustrine settings. It is not until the Neoarchean that signs of oxygenated surface waters above outer shelf and basinal depositional environments become evident. In order to further explore the Mesoarchean redox landscape for signs of basinward surface water oxygenation, we present nitrogen and carbon isotope ratios from the turbiditic Mosquito Creek Formation of the Nullagine Group (?2.9?Ga). The ?15N and ?13Corg values are invariant around ?1.8‰ and ?32‰ respectively throughout a 70?m section of drill-core, suggesting an ecosystem dominated by nitrogen fixers (anaerobic nitrogen cycling) and CO2 fixation by the Calvin Cycle. When compared with other Archean isotopic data, these results (i) provide further evidence that the Mosquito Creek Formation was deposited in a marine basin, and (ii) contain ?15N values that highlight the prevalence of nitrogen fixation by Mo-nitrogenase and the dearth of aerobic nitrogen metabolisms in the Mesoarchean.

Continue Reading →

Evolved Climates and Observational Discriminants for the TRAPPIST-1 Planetary System (The Astrophysical Journal, 2018)

The TRAPPIST-1 planetary system provides an unprecedented opportunity to study terrestrial exoplanet evolution with the James Webb Space Telescope (JWST) and ground-based observatories. Since M dwarf planets likely experience extreme volatile loss, the TRAPPIST-1 planets may have highly evolved, possibly uninhabitable atmospheres. We used a versatile, 1D terrestrial planet climate model with line-by-line radiative transfer and mixing length convection (VPL Climate) coupled to a terrestrial photochemistry model to simulate environmental states for the TRAPPIST-1 planets. We present equilibrium climates with self-consistent atmospheric compositions and observational discriminants of postrunaway, desiccated, 10–100 bar O2- and CO2-dominated atmospheres, including interior outgassing, as well as for water-rich compositions

Continue Reading →

Evaluating Climate Sensitivity to CO2 Across Earth’s History (JGR: Atmospheres, 2018)

CO2?driven changes to climate have occurred during many epochs of Earth’s history when the solar insolation, atmospheric CO2 concentration, and surface temperature of the planet were all significantly different than today. Each of these aspects affects the implied radiative forcings, climate feedbacks, and resultant changes in global mean surface temperature. Here we use a three?dimensional climate system model to study the effects of increasing CO2 on Earth’s climate, across many orders of magnitude of variation, and under solar inputs relevant for paleo, present, and future Earth scenarios. We find that the change in global mean surface temperature from doubling CO2 (i.e., the equilibrium climate sensitivity) may vary between 2.6 and 21.6 K over the course of Earth’s history. In agreement with previous studies, we find that the adjusted radiative forcing from doubling CO2 increases at high concentrations up to about 1.5 bars partial pressure, generally resulting in larger changes in the surface temperature. We also find that the cloud albedo feedback causes an abrupt transition in climate for warming atmospheres that depends both on the mean surface temperature and the total solar insolation. Climate sensitivity to atmospheric CO2 has probably varied considerably across Earth’s history.

Continue Reading →

Identifying Inflated Super-Earths and Photo-Evaporated Cores (The Astrophysical Journal, 2018)

We present empirical evidence, supported by a planet formation model, to show that the curve $R/{R}_{oplus }=1.05{(F/{F}_{oplus })}^{0.11}$ approximates the location of the so-called photo-evaporation valley. Planets below that curve are likely to have experienced complete photo-evaporation, and planets just above it appear to have inflated radii; thus we identify a new population of inflated super-Earths and mini-Neptunes. Our N-body simulations are set within an evolving protoplanetary disk and include prescriptions for orbital migration, gas accretion, and atmospheric loss due to giant impacts. Our simulated systems broadly match the sizes and periods of super-Earths in the Kepler catalog. They also reproduce the relative sizes of adjacent planets in the same system, with the exception of planet pairs that straddle the photo-evaporation valley. This latter group is populated by planet pairs with either very large or very small size ratios (R out/R in Gt 1 or R out/R in Lt 1) and a dearth of size ratios near unity. It appears that this feature could be reproduced if the planet outside the photo-evaporation valley (typically the outer planet, but sometimes not) has its atmosphere significantly expanded by stellar irradiation. This new population of planets may be ideal targets for future transit spectroscopy observations with the upcoming James Webb Space Telescope.

Continue Reading →

Migration-driven Diversity of Super-Earth Compositions (MNRAS Letters, 2018)

A leading model for the origin of super-Earths proposes that planetary embryos migrate inward and pile up on close-in orbits. As large embryos are thought to preferentially form beyond the snowline, this naively predicts that most super-Earths should be very water-rich. Here we show that the shortest period planets formed in the migration model are often purely rocky. The inward migration of icy embryos through the terrestrial zone accelerates the growth of rocky planets via resonant shepherding. We illustrate this process with a simulation that provided a match to the Kepler-36 system of two planets on close orbits with very different densities. In the simulation, two super-Earths formed in a Kepler-36-like configuration; the inner planet was pure rock while the outer one was ice-rich. We conclude from a suite of simulations that the feeding zones of close-in super-Earths are likely to be broad and disconnected from their final orbital radii.

Continue Reading →

Excitation of a Primordial Cold Asteroid Belt as an Outcome of Planetary Instability (The Astrophysical Journal, 2018)

The main asteroid belt (MB) is low in mass but dynamically excited. Here we propose a new mechanism to excite the MB during the giant planet (the “Nice model”) instability, which is expected to feature repeated close encounters between Jupiter and one or more ice giants (“jumping Jupiter” or JJ). We show that, when Jupiter temporarily reaches a high-enough level of excitation, both in eccentricity and inclination, it induces strong forced vectors of eccentricity and inclination across the MB region. Because during the JJ instability Jupiter’s orbit “jumps” around, the forced vectors keep changing both in magnitude and phase throughout the whole MB region. The entire cold primordial MB is thus excited as a natural outcome of the JJ instability. The level of such an excitation, however, is typically larger than the current orbital excitation observed in the MB. We show that the subsequent evolution of the solar system is capable of reshaping the resultant overexcited MB to its present-day orbital state, and that a strong mass depletion (~90%) is associated with the JJ instability phase and its subsequent evolution throughout the age of the solar system.

Continue Reading →

A Search for Refraction in Kepler Photometry of Gas Giants (Research Notes of the AAS, 2018)

Refraction can lead to a brightening just before ingress and just after egress of a transit, as light passes through the exoplanet’s atmosphere and is refracted into our line of sight (Sidis & Sari 2010; Misra & Meadows 2014; Misra et al. 2014; Dalba 2017; Alp & Demory 2018). Refraction just outside of transit has been seen and modeled in our own solar system during transits of Venus (Pasachoff et al. 2011; García Muñoz & Mills 2012; Tanga et al. 2012). For short-period planets, the model of (Sidis & Sari 2010, hereafter S&S) implies refraction peaks typically under 100 parts per million (ppm) and comparable in duration to ingress and egress. Kepler photometry (Borucki et al. 2010) currently provides the best opportunity for detecting refraction. We search for the signature of refraction just outside of transit in Kepler photometry of 45 gas giants and firmly rule out the S&S model for four candidates.

We select Kepler Objects of Interest (KOIs) with radii at least twice that of Earth for which the S&S Equation (30) implies a peak effect greater than 10 parts per million (ppm), adjusted for Rayleigh scattering using their Equations (40)–(45). We eliminate KOIs with grazing transits as well as those identified in Ford et al. (2012), Mazeh et al. (2013), and Holczer et al. (2016) as having significant transit timing variations. We also eliminate a few KOIs identified by Holczer et al. (2016) as likely planetary false positives based on the behavior of the light curves, leaving 45 planet candidates. To calculate the expected effect, we adopt the masses predicted in Chen & Kipping (2018).

Continue Reading →

Formation of Terrestrial Planets in Eccentric and Inclined Giant Planet Systems (Astronomy & Astrophysics, 2018)

Aims. Evidence of mutually inclined planetary orbits has been reported for giant planets in recent years. Here we aim to study the impact of eccentric and inclined massive giant planets on the terrestrial planet formation process, and investigate whether it can possibly lead to the formation of inclined terrestrial planets.

Methods. We performed 126 simulations of the late-stage planetary accretion in eccentric and inclined giant planet systems. The physical and orbital parameters of the giant planet systems result from n-body simulations of three giant planets in the late stage of the gas disc, under the combined action of Type II migration and planet-planet scattering. Fourteen two- and three-planet configurations were selected, with diversified masses, semi-major axes (resonant configurations or not), eccentricities, and inclinations (including coplanar systems) at the dispersal of the gas disc. We then followed the gravitational interactions of these systems with an inner disc of planetesimals and embryos (nine runs per system), studying in detail the final configurations of the formed terrestrial planets.

Results. In addition to the well-known secular and resonant interactions between the giant planets and the outer part of the disc, giant planets on inclined orbits also strongly excite the planetesimals and embryos in the inner part of the disc through the combined action of nodal resonance and the Lidov–Kozai mechanism. This has deep consequences on the formation of terrestrial planets. While coplanar giant systems harbour several terrestrial planets, generally as massive as the Earth and mainly on low-eccentric and low-inclined orbits, terrestrial planets formed in systems with mutually inclined giant planets are usually fewer, less massive (<0.5 M?), and with higher eccentricities and inclinations. This work shows that terrestrial planets can form on stable inclined orbits through the classical accretion theory, even in coplanar giant planet systems emerging from the disc phase.

Continue Reading →

Linearized Flux Evolution (LiFE): A technique for rapidly adapting fluxes from full-physics radiative transfer models (Journal of Quantitative Spectroscopy and Radiative Transfer, 2018)

Solar and thermal radiation are critical aspects of planetary climate, with gradients in radiative energy fluxes driving heating and cooling. Climate models require that radiative transfer tools be versatile, computationally efficient, and accurate. Here, we describe a technique that uses an accurate full-physics radiative transfer model to generate a set of atmospheric radiative quantities which can be used to linearly adapt radiative flux profiles to changes in the atmospheric and surface state—the Linearized Flux Evolution (LiFE) approach. These radiative quantities describe how each model layer in a plane-parallel atmosphere reflects and transmits light, as well as how the layer generates diffuse radiation by thermal emission and by scattering light from the direct solar beam. By computing derivatives of these layer radiative properties with respect to dynamic elements of the atmospheric state, we can then efficiently adapt the flux profiles computed by the full-physics model to new atmospheric states. We validate the LiFE approach, and then apply this approach to Mars, Earth, and Venus, demonstrating the information contained in the layer radiative properties and their derivatives, as well as how the LiFE approach can be used to determine the thermal structure of radiative and radiative-convective equilibrium states in one-dimensional atmospheric models.

Continue Reading →

Mars’ Growth Stunted by an Early Giant Planet Instability (Icarus, 2018)

Many dynamical aspects of the solar system can be explained by the outer planets experiencing a period of orbital instability sometimes called the Nice Model. Though often correlated with a perceived delayed spike in the lunar cratering record known as the Late Heavy Bombardment (LHB), recent work suggests that this event may have occurred much earlier; perhaps during the epoch of terrestrial planet formation. While current simulations of terrestrial accretion can reproduce many observed qualities of the solar system, replicating the small mass of Mars requires modification to standard planet formation models. Here we use 800 dynamical simulations to show that an early instability in the outer solar system strongly influences terrestrial planet formation and regularly yields properly sized Mars analogs. Our most successful outcomes occur when the terrestrial planets evolve an additional 1–10 million years (Myr) following the dispersal of the gas disk, before the onset of the giant planet instability. In these simulations, accretion has begun in the Mars region before the instability, but the dynamical perturbation induced by the giant planets’ scattering removes large embryos from Mars’ vicinity. Large embryos are either ejected or scattered inward toward Earth and Venus (in some cases to deliver water), and Mars is left behind as a stranded embryo. An early giant planet instability can thus replicate both the inner and outer solar system in a single model.

Continue Reading →

The Habitability of Proxima Centauri b: Environmental States and Observational Discriminants (Astrobiology, 2018)

Proxima Centauri b provides an unprecedented opportunity to understand the evolution and nature of terrestrial planets orbiting M dwarfs. Although Proxima Cen b orbits within its star’s habitable zone, multiple plausible evolutionary paths could have generated different environments that may or may not be habitable. Here, we use 1-D coupled climate-photochemical models to generate self-consistent atmospheres for several evolutionary scenarios, including high-O2, high-CO2, and more Earth-like atmospheres, with both oxic and anoxic compositions. We show that these modeled environments can be habitable or uninhabitable at Proxima Cen b’s position in the habitable zone. We use radiative transfer models to generate synthetic spectra and thermal phase curves for these simulated environments, and use instrument models to explore our ability to discriminate between possible planetary states.

Continue Reading →

Exo-Milankovitch Cycles. I. Orbits and Rotation States (The Astronomical Journal, 2018)

The obliquity of the Earth, which controls our seasons, varies by only ~2fdg5 over ~40,000 years, and its eccentricity varies by only ~0.05 over 100,000 years. Nonetheless, these small variations influence Earth’s ice ages. For exoplanets, however, variations can be significantly larger. Previous studies of the habitability of moonless Earth-like exoplanets have found that high obliquities, high eccentricities, and dynamical variations can extend the outer edge of the habitable zone by preventing runaway glaciation (snowball states). We expand upon these studies by exploring the orbital dynamics with a semianalytic model that allows us to map broad regions of parameter space. We find that, in general, the largest drivers of obliquity variations are secular spin–orbit resonances. We show how the obliquity varies in several test cases, including Kepler-62 f, across a wide range of orbital and spin parameters. These obliquity variations, alongside orbital variations, will have a dramatic impact on the climates of such planets.

Continue Reading →