Publications

2024

1. Radiative Signatures of Circumplanetary Disks and Envelopes During the Late Stages of Giant Planet Formation

   Aster G. Taylor, and Fred C. Adams

   Icarus, 2024

   Abs

   During the late stages of giant planet formation, protoplanets are surrounded by a circumplanetary disk and an infalling envelope of gas and dust. For systems with sufficient cooling, material entering the sphere of influence of the planet falls inward and approaches ballistic conditions. Due to conservation of angular momentum, most of the incoming material falls onto the disk rather than directly onto the planet. This paper determines the spectral energy distributions of forming planets in this stage of evolution. Generalizing previous work, we consider a range of possible geometries for the boundary conditions of the infall and determine the two-dimensional structure of the envelope, as well as the surface density of the disk. After specifying the luminosity sources for the planet and disk, we calculate the corresponding radiative signatures for the system, including the emergent spectral energy distributions and emission maps. These results show how the observational appearance of forming planets depend on the input parameters, including the instantaneous mass, mass accretion rate, semimajor axis of the orbit, and the planetary magnetic field strength (which sets the inner boundary condition for the disk). We also consider different choices for the form of the opacity law and attenuation due to the background circumstellar disk. Although observing forming planets will be challenging, these results show how the observational signatures depend on the underlying properties of the planet/disk/envelope system.

2. The Dynamical Origins of the Dark Comets and a Proposed Evolutionary Track

   Aster G. Taylor, Jordan K. Steckloff, Darryl Z. Seligman, and 5 more authors

   Icarus, Jul 2024

   Abs

   So-called ‘dark comets’ are small, morphologically inactive near-Earth objects (NEOs) that exhibit nongravitational accelerations inconsistent with radiative effects. These objects exhibit short rotational periods (minutes to hours), where measured. We find that the strengths required to prevent catastrophic disintegration are consistent with those measured in cometary nuclei and expected in rubble pile objects. We hypothesize that these dark comets are the end result of a rotational fragmentation cascade, which is consistent with their measured physical properties. We calculate the predicted size-frequency distribution for objects evolving under this model. Using dynamical simulations, we further demonstrate that the majority of these bodies originated from the ν6 resonance, implying the existence of volatiles in the current inner main belt. Moreover, one of the dark comets, (523599) 2003 RM, likely originated from the outer main belt, although a JFC origin is also plausible. These results provide strong evidence that volatiles from a reservoir in the inner main belt are present in the near-Earth environment.

3. Formation and Structure of Circumplanetary Disks and Envelopes During the Late Stages of Giant Planet Formation

   Aster G. Taylor, and Fred C. Adams

   Icarus, Jun 2024

   Abs

   Giant planets are expected to form within circumplanetary disks, which shape their formation history and the local environment. Here, we consider the formation and structure of circumplanetary disks that arise during the late stages of giant planet formation. During this phase, when most of the final mass is accumulated, incoming material enters the Hill sphere and falls toward the planet. In the absence of torques, the falling parcels of gas conserve their specific angular momentum and collect into a circumplanetary disk. Generalizing previous work, we consider a range of possible geometries for the flow entering the sphere of influence of the planet. Specifically, we consider five geometric patterns for the inward flow, ranging from concentration toward the rotational poles of the system to isotropic flow to concentration along the equatorial plane. For each case, we derive analytic descriptions for the density field of the infall region, the disk surface density in the absence of viscosity, and steady-state solutions for viscous disks. These results, in turn, specify the luminosity contributions of the planet, the circumplanetary disk, and the envelope. These power sources, in conjunction with the surrounding material, collectively determine the observational appearance of the forming planet. We conclude with an approximate determination of these radiative signatures.

4. Potential Melting of Extrasolar Planets by Tidal Dissipation

   Darryl Z. Seligman, Adina D. Feinstein, Dong Lai, and 6 more authors

   Astrophysical Journal, Jan 2024

   Abs

   Tidal heating on Io due to its finite eccentricity was predicted to drive surface volcanic activity, which was subsequently confirmed by the Voyager spacecraft. Although the volcanic activity in Io is more complex, in theory volcanism can be driven by runaway melting in which the tidal heating increases as the mantle thickness decreases. We show that this runaway melting mechanism is generic for a composite planetary body with liquid core and solid mantle, provided that (i) the mantle rigidity, μ, is comparable to the central pressure, i.e., μ/(ρ gR P) ≳ 0.1 for a body with density ρ, surface gravitational acceleration g, and radius R P; (ii) the surface is not molten; (iii) tides deposit sufficient energy; and (iv) the planet has nonzero eccentricity. We calculate the approximate liquid core radius as a function of μ/(ρ gR P), and find that more than 90% of the core will melt due to this runaway for μ/(ρ gR P) ≳ 1. From all currently confirmed exoplanets, we find that the terrestrial planets in the L 98-59 system are the most promising candidates for sustaining active volcanism. However, uncertainties regarding the quality factors and the details of tidal heating and cooling mechanisms prohibit definitive claims of volcanism on any of these planets. We generate synthetic transmission spectra of these planets assuming Venus-like atmospheric compositions with an additional 5%, 50%, and 98% SO2 component, which is a tracer of volcanic activity. We find a ≳3σ preference for a model with SO2 with 5–10 transits with JWST for L 98-59bcd.

5. Seasonally Varying Outgassing as an Explanation for Dark Comet Accelerations

   Aster G. Taylor, Davide Farnocchia, David Vokrouhlický, and 3 more authors

   Icarus, Jan 2024

   Abs

   Significant nonradial, nongravitational accelerations with magnitudes incompatible with radiation-driven effects have been reported in seven photometrically inactive near-Earth objects. Two of these objects exhibit large transverse accelerations (i.e., within the orbital plane but orthogonal to the radial direction), and six exhibit significant out-of-plane accelerations. Here, we find that anisotropic outgassing resulting from differential heating on a nucleus with nonzero spin-pole obliquity, averaged over an eccentric orbit, can explain these accelerations for most of the objects. This balanced outgassing model depends on three parameters — the spin pole orientation (R.A. and Dec.) and an acceleration magnitude. For these “dark comets” (excepting 2003 RM), we obtain parameter values that reproduce the observed nongravitational accelerations. We derive formulae for the component accelerations under certain assumptions for the acceleration scaling over heliocentric distance. Although we lack estimates of these objects’ spin axes to confirm our values, this mechanism is nevertheless a plausible explanation for the observed accelerations, and produces accurate perturbations to the heliocentric motions of most of these objects. This model may also be applied to active objects outside of the dark comets group.

2023

1. Fitting the Light Curve of 1I/’Oumuamua with a Nonprincipal Axis Rotational Model and Outgassing Torques

   Aster G. Taylor, Darryl Z. Seligman, Olivier R. Hainaut, and 1 more author

   Planetary Science Journal, Oct 2023

   Abs

   In this paper, we investigate the nonprincipal axis (NPA) rotational state of 1I/‘Oumuamua—the first interstellar object discovered traversing the inner solar system—from its photometric light curve. Building upon Mashchenko, we develop a model which incorporates NPA rotation and Sun-induced, time-varying outgassing torques to generate synthetic light curves of the object. The model neglects tidal forces, which are negligible compared to outgassing torques over the distances at which ‘Oumuamua was observed. We implement an optimization scheme that incorporates the NPA rotation model to calculate the initial rotation state of the object. We find that an NPA rotation state with an average period of 〈P〉 ≃ 7.34 hr best reproduces the photometric data. The discrepancy between this period and previous estimates is due to continuous period modulation induced by outgassing torques in the rotational model, as well as different periods being used. The best fit to the 2017 October data does not reproduce the 2017 November data (although the later measurements are too sparse to fit). The light curve is consistent with there being no secular evolution of the angular momentum, which is somewhat in tension with the empirical correlations between nuclear spin-up and cometary outgassing. The complex rotation of ‘Oumuamua may be the result of primordial rotation about the smallest principal axis if (i) the object experienced hypervolatile outgassing and (ii) our idealized outgassing model is accurate.

2. Interstellar Comets from Post-Main-Sequence Systems as Tracers of Extrasolar Oort Clouds

   W. Garrett Levine, Aster G. Taylor, Darryl Z. Seligman, and 4 more authors

   Planetary Science Journal, Jul 2023

   Abs

   Interstellar small bodies are unique probes into the histories of exoplanetary systems. One hypothesized class of interlopers are "Jurads," exocomets released into the Milky Way during the post-main-sequence as the thermally pulsing asymptotic giant branch (AGB) host stars lose mass. In this study, we assess the prospects for the Legacy Survey of Space and Time (LSST) to detect a Jurad and examine whether such an interloper would be observationally distinguishable from exocomets ejected during the (pre-)main-sequence. Using analytic and numerical methods, we estimate the fraction of exo-Oort Cloud objects that are released from 1-8 M ⊙ stars during post-main-sequence evolution. We quantify the extent to which small bodies are altered by the increased luminosity and stellar outflows during the AGB, finding that some Jurads may lack hypervolatiles and that stellar winds could deposit dust that covers the entire exocomet surface. Next, we construct models of the interstellar small body reservoir for various size-frequency distributions and examine the LSST’s ability to detect members of those hypothesized populations. Combining these analyses, we highlight the joint constraints that the LSST will place on power-law size-frequency distribution slopes, characteristic sizes, and the total mass sequestered in the minor planets of exo-Oort Clouds. Even with the LSST’s increased search volume compared to contemporary surveys, we find that detecting a Jurad is unlikely but not infeasible given the current understanding of (exo)planet formation.

3. Numerical Simulations of Tidal Deformation and Resulting Light Curves of Small Bodies: Material Constraints of 99942 Apophis and 1I/’Oumuamua

   Aster G. Taylor, Darryl Z. Seligman, Douglas R. MacAyeal, and 2 more authors

   Planetary Science Journal, May 2023

   Abs

   In this paper, we present an open-source software (Simulator of Asteroid Malformation Under Stress, SAMUS) that simulates constant-density, constant-viscosity liquid bodies subject to tidal forces for a range of assumed viscosities and sizes. This software solves the Navier-Stokes equations on a finite-element mesh, incorporating the centrifugal, Coriolis, self-gravitational, and tidal forces. The primary functionality is to simulate the deformation of minor bodies under the influence of tidal forces. It may therefore be used to constrain the composition and physical structure of bodies experiencing significant tidal forces, such as 99942 Apophis and 1I/’Oumuamua. We demonstrate that SAMUS will be useful to constrain the material properties of Apophis during its near-Earth flyby in 2029. Depending on the material properties, Apophis may experience an area change of up to 0.5%, with similar effects on the photometric brightness. We also apply SAMUS to constrain the material dynamic viscosity of 1I/’Oumuamua, the first interstellar object discovered traversing the inner solar system. ’Oumuamua experienced a close approach to the Sun at perihelion (q ≃ 0.25 au) during which there were significant tidal forces that may have caused deformation of the body. This deformation could have lead to observable changes in the photometric light curve based on the material properties. The application of SAMUS to produce synthetic observations which incorporate tidal deformation effects demonstrates that no deformation-an infinite dynamic viscosity-best reproduces the photometric data. While these results indicate that ’Oumuamua did not experience significant tidal deformation, a sophisticated model incorporating nonprincipal axis rotation is necessary to conclusively analyze both ’Oumuamua and Apophis.

4. Dark Comets? Unexpectedly Large Nongravitational Accelerations on a Sample of Small Asteroids

   Darryl Z. Seligman, Davide Farnocchia, Marco Micheli, and 15 more authors

   Planetary Science Journal, Feb 2023

   Abs

   We report statistically significant detections of nonradial, nongravitational accelerations based on astrometric data in the photometrically inactive objects 1998 KY26, 2005 VL1, 2016 NJ33, 2010 VL65, 2016 RH120, and 2010 RF12. The magnitudes of the nongravitational accelerations are greater than those typically induced by the Yarkovsky effect, and there is no radiation-based, nonradial effect that can be so large. Therefore, we hypothesize that the accelerations are driven by outgassing and calculate implied H2O production rates for each object. We attempt to reconcile outgassing-induced acceleration with the lack of visible comae or photometric activity via the absence of surface dust and low levels of gas production. Although these objects are small, and some are rapidly rotating, the surface cohesive forces are stronger than the rotational forces, and rapid rotation alone cannot explain the lack of surface debris. It is possible that surface dust was removed previously, perhaps via outgassing activity that increased the rotation rates to their present-day value. We calculate dust production rates of order  10-4 g s-1 in each object, assuming that the nuclei are bare, within the upper limits of dust production from a sample stacked image of 1998 KY26 of M˙Dust<0.2 g s-1. This production corresponds to brightness variations of order  0.0025%, which are undetectable in extant photometric data. We assess the future observability of each of these targets and find that the orbit of 1998 KY26-which is also the target of the extended Hayabusa2 mission-exhibits favorable viewing geometry before 2025.