Showing 1–50 of 193 results for author: Adams, F C

Exploring the Complex Ionization Environment of the Turbulent DM Tau Disk

Authors: Deryl E. Long, L. Ilsedore Cleeves, Fred C. Adams, Sean Andrews, Edwin A. Bergin, Viviana V. Guzmán, Jane Huang, A. Meredith Hughes, Chunhua Qi, Kamber Schwarz, Jacob B. Simon, David Wilner

Abstract: Ionization drives important chemical and dynamical processes within protoplanetary disks, including the formation of organics and water in the cold midplane and the transportation of material via accretion and magneto-hydrodynamic (MHD) flows. Understanding these ionization-driven processes is crucial for understanding disk evolution and planet formation. We use new and archival ALMA observations… ▽ More Ionization drives important chemical and dynamical processes within protoplanetary disks, including the formation of organics and water in the cold midplane and the transportation of material via accretion and magneto-hydrodynamic (MHD) flows. Understanding these ionization-driven processes is crucial for understanding disk evolution and planet formation. We use new and archival ALMA observations of HCO+, H13CO+, and N2H+ to produce the first forward-modeled 2D ionization constraints for the DM Tau protoplanetary disk. We include ionization from multiple sources and explore the disk chemistry under a range of ionizing conditions. Abundances from our 2D chemical models are post-processed using non-LTE radiative transfer, visibility sampling, and imaging, and are compared directly to the observed radial emission profiles. The observations are best fit by a modestly reduced CR ionization rate ($ζ_{CR}$ ~ 10$^{-18}$ s$^{-1}$) and a hard X-ray spectrum (hardness ratio [HR] = 0.3), which we associate with stellar flaring conditions. Our best-fit model under-produces emission in the inner disk, suggesting that there may be an additional mechanism enhancing ionization in DM Tau's inner disk. Overall, our findings highlight the complexity of ionization in protoplanetary disks and the need for high resolution multi-line studies. △ Less

Submitted 26 June, 2024; originally announced June 2024.

Comments: 18 pages, 12 figures, accepted to be published in The Astrophysical Journal (June 25, 2024)

Formation and Structure of Circumplanetary Disks and Envelopes during the Late Stages of Giant Planet Formation

Authors: Aster G. Taylor, Fred C. Adams

Abstract: Giant planets are expected to form within circumstellar 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 ab… ▽ More Giant planets are expected to form within circumstellar 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. △ Less

Submitted 12 March, 2024; v1 submitted 8 March, 2024; originally announced March 2024.

Comments: 19 pages, 9 figures. Published at Icarus

Journal ref: Icarus (2024) 116044

Can Neptune's Distant Mean-Motion Resonances Constrain Undiscovered Planets in the Solar System? Lessons from a Case Study of the 9:1

Authors: Matthew W. Porter, David W. Gerdes, Kevin J. Napier, Hsing Wen Lin, Fred C. Adams

Abstract: Recent observational surveys of the outer Solar System provide evidence that Neptune's distant $n$:1 mean-motion resonances may harbor relatively large reservoirs of trans-Neptunian objects (TNOs). In particular, the discovery of two securely classified 9:1 resonators, 2015 KE$_{172}$ and 2007 TC$_{434}$, by the Outer Solar System Origins Survey is consistent with a population of order $10^4$ such… ▽ More Recent observational surveys of the outer Solar System provide evidence that Neptune's distant $n$:1 mean-motion resonances may harbor relatively large reservoirs of trans-Neptunian objects (TNOs). In particular, the discovery of two securely classified 9:1 resonators, 2015 KE$_{172}$ and 2007 TC$_{434}$, by the Outer Solar System Origins Survey is consistent with a population of order $10^4$ such objects in the 9:1 resonance with absolute magnitude $H_r < 8.66$. This work investigates whether the long-term stability of such populations in Neptune's $n$:1 resonances can be used to constrain the existence of distant $5-10M_{\oplus}$ planets orbiting at hundreds of AU. The existence of such a planet has been proposed to explain a reported clustering in the orbits of highly eccentric "extreme" trans-Neptunian objects (eTNOs), although this hypothesis remains controversial. We engage in a focused computational case-study of the 9:1 resonance, generating synthetic populations and integrating them for 1 Gyr in the presence of 81 different test planets with various masses, perihelion distances, eccentricities, and inclinations. While none of the tested planets are incompatible with the existence of 9:1 resonators, our integrations shed light on the character of the interaction between such planets and nearby $n$:1 resonances, and we use this knowledge to construct a simple, heuristic method for determining whether or not a given planet could destabilize a given resonant population. We apply this method to the currently estimated properties of Planet 9, and find that a large primordial population in the 15:1 resonance (or beyond), if discovered in the future, could potentially constrain the existence of this planet. △ Less

Submitted 31 January, 2024; originally announced February 2024.

Comments: 20 pages, 16 figures

Generalized Models for Inflationary Preheating: Oscillations and Symmetries

Authors: Leia Barrowes, Fred C. Adams, Anthony M. Bloch, Scott Watson

Abstract: The paradigm of the inflationary universe provides a possible explanation for several observed cosmological properties. In order for such solutions to be successful, the universe must convert the energy stored in the inflaton potential into standard model particles through a process known as reheating. In this paper, we reconsider the reheating process for the case where the inflaton potential res… ▽ More The paradigm of the inflationary universe provides a possible explanation for several observed cosmological properties. In order for such solutions to be successful, the universe must convert the energy stored in the inflaton potential into standard model particles through a process known as reheating. In this paper, we reconsider the reheating process for the case where the inflaton potential respects an approximate (but spontaneously broken) conformal symmetry during the reheating epoch. After reviewing the Effective Field Theory of Reheating, we present solutions for the nonlinear oscillations of the inflaton field, derive the corresponding Hill's equation for the coupled reheating field, and determine the stability diagram for parametric resonance. For this class of models -- the simplest realization being a scalar field with a quartic term -- the expansion of the universe drives the coupled field toward a more unstable part of parameter space, in contrast to the standard case. We also generalize this class of models to include quadratic breaking terms in the potential during the reheating epoch and address the process of stability in that universality class of models. △ Less

Submitted 23 January, 2024; originally announced January 2024.

Comments: 35 pages, 15 figures

Potential Melting of Extrasolar Planets by Tidal Dissipation

Authors: Darryl Z. Seligman, Adina D. Feinstein, Dong Lai, Luis Welbanks, Aster G. Taylor, Juliette Becker, Fred C. Adams, Marvin Morgan, Jennifer B. Bergner

Abstract: Tidal heating on Io due to its finite eccentricity was predicted to drive surface volcanic activity, which was subsequently confirmed by the $\textit{Voyager}$ spacecrafts. 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 g… ▽ More Tidal heating on Io due to its finite eccentricity was predicted to drive surface volcanic activity, which was subsequently confirmed by the $\textit{Voyager}$ spacecrafts. 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. $μ/ (ρg R_{\rm P})\gtrsim0.1$ for a body with density $ρ$, surface gravitational acceleration $g$, and radius $R_{\rm 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 $μ/ (ρg R_{\rm P})$, and find that more than $90\%$ of the core will melt due to this runaway for $μ/ (ρg R_{\rm P})\gtrsim1$. From all currently confirmed exoplanets, we find that the terrestrial planets in the L98-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\%$ SO$_2$ component, which is a tracer of volcanic activity. We find a $\gtrsim 3 σ$ preference for a model with SO$_2$ with 5-10 transits with $\textit{JWST}$ for L98-59bcd. △ Less

Submitted 2 November, 2023; originally announced November 2023.

Comments: 16 pages, 8 Figures, accepted for publication in ApJ

The DECam Ecliptic Exploration Project (DEEP) II. Observational Strategy and Design

Authors: Chadwick A. Trujillo, Cesar Fuentes, David W. Gerdes, Larissa Markwardt, Scott S. Sheppard, Ryder Strauss, Colin Orion Chandler, William J. Oldroyd, David E. Trilling, Hsing Wen Lin, Fred C. Adams, Pedro H. Bernardinelli, Matthew J. Holman, Mario Juric, Andrew McNeill, Michael Mommert, Kevin J. Napier, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Hayden Smotherman

Abstract: We present the DECam Ecliptic Exploration Project (DEEP) survey strategy including observing cadence for orbit determination, exposure times, field pointings and filter choices. The overall goal of the survey is to discover and characterize the orbits of a few thousand Trans-Neptunian Objects (TNOs) using the Dark Energy Camera (DECam) on the Cerro Tololo Inter-American Observatory (CTIO) Blanco 4… ▽ More We present the DECam Ecliptic Exploration Project (DEEP) survey strategy including observing cadence for orbit determination, exposure times, field pointings and filter choices. The overall goal of the survey is to discover and characterize the orbits of a few thousand Trans-Neptunian Objects (TNOs) using the Dark Energy Camera (DECam) on the Cerro Tololo Inter-American Observatory (CTIO) Blanco 4 meter telescope. The experiment is designed to collect a very deep series of exposures totaling a few hours on sky for each of several 2.7 square degree DECam fields-of-view to achieve a magnitude of about 26.2 using a wide VR filter which encompasses both the V and R bandpasses. In the first year, several nights were combined to achieve a sky area of about 34 square degrees. In subsequent years, the fields have been re-visited to allow TNOs to be tracked for orbit determination. When complete, DEEP will be the largest survey of the outer solar system ever undertaken in terms of newly discovered object numbers, and the most prolific at producing multi-year orbital information for the population of minor planets beyond Neptune at 30 au. △ Less

Submitted 30 October, 2023; originally announced October 2023.

Comments: 29 pages, 4 figures and 4 tables

First Near-IR Spectroscopic Survey of Neptune Trojans with JWST: Distinct Surface Compositions of Red vs Ultra-Red Neptune Trojans

Authors: Larissa Markwardt, Bryan J. Holler, Hsing Wen Lin, David W. Gerdes, Fred C. Adams, Renu Malhotra, Kevin J. Napier

Abstract: Neptune's Trojan asteroids have been observed to have a variety of optical colors, most notably red (g $-$ r < 0.75) vs. ultra-red (g $-$ r > 0.75), but the underlying cause of these different color classifications is unknown. Near-IR spectroscopy can be used as a probe of the surface composition of these objects, as broad ice bands for a variety of materials are present in the near-IR. Here, we p… ▽ More Neptune's Trojan asteroids have been observed to have a variety of optical colors, most notably red (g $-$ r < 0.75) vs. ultra-red (g $-$ r > 0.75), but the underlying cause of these different color classifications is unknown. Near-IR spectroscopy can be used as a probe of the surface composition of these objects, as broad ice bands for a variety of materials are present in the near-IR. Here, we present the first results of a spectroscopic survey of Neptune's Trojan asteroids using the NIRSpec instrument on JWST. We compare the near-IR spectra of eight Neptune Trojans (NTs) based on different optical color classifications and with model spectra of different ices. We find that most of our targets are consistent with a surface covered in a thin layer of H$_2$O and CO$_2$ ices, while the only NT to reliably be classified as ultra-red is covered in ice tholins in addition to CO$_2$. Ice tholins are a known reddening agent when subjected to irradiation, so these results support the hypothesis that differences in optical color are due to differences in irradiation of the surfaces of these bodies. Since NTs have very similar orbits and therefore generally similar levels of irradiation at the current time, our results suggest that these objects have unique origins or there is ongoing processing of the surfaces of these objects through stochastic disturbances such as impacts. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 10 pages, 2 figures, Submitted to ApJL

The DECam Ecliptic Exploration Project (DEEP) VI: first multi-year observations of trans-Neptunian objects

Authors: Hayden Smotherman, Pedro H. Bernardinelli, Stephen K. N. Portillo, Andrew J. Connolly, J. Bryce Kalmbach, Steven Stetzler, Mario Juric, Dino Bektesvic, Zachary Langford, Fred C. Adams, William J. Oldroyd, Matthew J. Holman, Colin Orion Chandler, Cesar Fuentes, David W. Gerdes, Hsing Wen Lin, Larissa Markwardt, Andrew McNeill, Michael Mommert, Kevin J. Napier, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Scott S. Sheppard , et al. (3 additional authors not shown)

Abstract: We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and… ▽ More We present the first set of trans-Neptunian objects (TNOs) observed on multiple nights in data taken from the DECam Ecliptic Exploration Project (DEEP). Of these 110 TNOs, 105 do not coincide with previously known TNOs and appear to be new discoveries. Each individual detection for our objects resulted from a digital tracking search at TNO rates of motion, using two to four hour exposure sets, and the detections were subsequently linked across multiple observing seasons. This procedure allows us to find objects with magnitudes $m_{VR} \approx 26$. The object discovery processing also included a comprehensive population of objects injected into the images, with a recovery and linking rate of at least $94\%$. The final orbits were obtained using a specialized orbit fitting procedure that accounts for the positional errors derived from the digital tracking procedure. Our results include robust orbits and magnitudes for classical TNOs with absolute magnitudes $H \sim 10$, as well as a dynamically detached object found at 76 au (semi-major axis $a\approx 77 \, \mathrm{au}$). We find a disagreement between our population of classical TNOs and the CFEPS-L7 three component model for the Kuiper belt. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted to AJ, companion paper do DEEP III. Objects will be released in the journal version (or contacting the authors)

The DECam Ecliptic Exploration Project (DEEP) III: Survey characterization and simulation methods

Authors: Pedro H. Bernardinelli, Hayden Smotherman, Zachary Langford, Stephen K. N. Portillo, Andrew J. Connolly, J. Bryce Kalmbach, Steven Stetzler, Mario Juric, William J. Oldroyd, Hsing Wen Lin, Fred C. Adams, Colin Orion Chandler, Cesar Fuentes, David W. Gerdes, Matthew J. Holman, Larissa Markwardt, Andrew McNeill, Michael Mommert, Kevin J. Napier, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Scott S. Sheppard, Ryder Strauss , et al. (2 additional authors not shown)

Abstract: We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characteriz… ▽ More We present a detailed study of the observational biases of the DECam Ecliptic Exploration Project's (DEEP) B1 data release and survey simulation software that enables direct statistical comparisons between models and our data. We inject a synthetic population of objects into the images, and then subsequently recover them in the same processing as our real detections. This enables us to characterize the survey's completeness as a function of apparent magnitudes and on-sky rates of motion. We study the statistically optimal functional form for the magnitude, and develop a methodology that can estimate the magnitude and rate efficiencies for all survey's pointing groups simultaneously. We have determined that our peak completeness is on average 80\% in each pointing group, and our magnitude drops to $25\%$ of this value at $m_{25} = 26.22$. We describe the freely available survey simulation software and its methodology. We conclude by using it to infer that our effective search area for objects at 40 au is $14.8°^2$, and that our lack of dynamically cold distant objects means that there at most $8\times 10^3$ objects with $60 < a < 80$ au and absolute magnitudes $H \leq 8$. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: Accepted to AJ, companion paper to DEEP VI

The DECam Ecliptic Exploration Project (DEEP): V. The Absolute Magnitude Distribution of the Cold Classical Kuiper Belt

Authors: Kevin J. Napier, Hsing-Wen Lin, David W. Gerdes, Fred C. Adams, Anna M. Simpson, Matthew W. Porter, Katherine G. Weber, Larissa Markwardt, Gabriel Gowman, Hayden Smotherman, Pedro H. Bernardinelli, Mario Jurić, Andrew J. Connolly, J. Bryce Kalmbach, Stephen K. N. Portillo, David E. Trilling, Ryder Strauss, William J. Oldroyd, Chadwick A. Trujillo, Colin Orion Chandler, Matthew J. Holman, Hilke E. Schlichting, Andrew McNeill, the DEEP Collaboration

Abstract: The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth,… ▽ More The DECam Ecliptic Exploration Project (DEEP) is a deep survey of the trans-Neptunian solar system being carried out on the 4-meter Blanco telescope at Cerro Tololo Inter-American Observatory in Chile using the Dark Energy Camera (DECam). By using a shift-and-stack technique to achieve a mean limiting magnitude of $r \sim 26.2$, DEEP achieves an unprecedented combination of survey area and depth, enabling quantitative leaps forward in our understanding of the Kuiper Belt populations. This work reports results from an analysis of twenty 3 sq.\ deg.\ DECam fields along the invariable plane. We characterize the efficiency and false-positive rates for our moving-object detection pipeline, and use this information to construct a Bayesian signal probability for each detected source. This procedure allows us to treat all of our Kuiper Belt Object (KBO) detections statistically, simultaneously accounting for efficiency and false positives. We detect approximately 2300 candidate sources with KBO-like motion at S/N $>6.5$. We use a subset of these objects to compute the luminosity function of the Kuiper Belt as a whole, as well as the Cold Classical (CC) population. We also investigate the absolute magnitude ($H$) distribution of the CCs, and find consistency with both an exponentially tapered power-law, which is predicted by streaming instability models of planetesimal formation, and a rolling power law. Finally, we provide an updated mass estimate for the Cold Classical Kuiper Belt of $M_{CC}(H_r < 12) = 0.0017^{+0.0010}_{-0.0004} M_{\oplus}$, assuming albedo $p = 0.15$ and density $ρ= 1$ g cm$^{-3}$. △ Less

Submitted 18 September, 2023; originally announced September 2023.

Comments: Accepted by PSJ

The DECam Ecliptic Exploration Project (DEEP) IV: Constraints on the shape distribution of bright TNOs

Authors: R. Strauss, D. E. Trilling, P. H. Bernardinelli, C. Beach, W. J. Oldroyd, S. S. Sheppard, H. E. Schlichting, D. W. Gerdes, F. C. Adams, C. O. Chandler, C. Fuentes, M. J. Holman, M. Jurić, H. W. Lin, L. Markwardt, A. McNeill, M. Mommert, K. J. Napier, M. J. Payne, D. Ragozzine, A. S. Rivkin, H. Smotherman, C. A. Trujillo

Abstract: We present the methods and results from the discovery and photometric measurement of 26 bright (VR $>$ 24 trans-Neptunian objects (TNOs) during the first year (2019-20) of the DECam Ecliptic Exploration Project (DEEP). The DEEP survey is an observational TNO survey with wide sky coverage, high sensitivity, and a fast photometric cadence. We apply a computer vision technique known as a progressive… ▽ More We present the methods and results from the discovery and photometric measurement of 26 bright (VR $>$ 24 trans-Neptunian objects (TNOs) during the first year (2019-20) of the DECam Ecliptic Exploration Project (DEEP). The DEEP survey is an observational TNO survey with wide sky coverage, high sensitivity, and a fast photometric cadence. We apply a computer vision technique known as a progressive probabilistic Hough transform to identify linearly-moving transient sources within DEEP photometric catalogs. After subsequent visual vetting, we provide a photometric and astrometric catalog of our TNOs. By modeling the partial lightcurve amplitude distribution of the DEEP TNOs using Monte Carlo techniques, we find our data to be most consistent with an average TNO axis ratio b/a $<$ 0.5, implying a population dominated by non-spherical objects. Based on ellipsoidal gravitational stability arguments, we find our data to be consistent with a TNO population containing a high fraction of contact binaries or other extremely non-spherical objects. We also discuss our data as evidence that the expected binarity fraction of TNOs may be size-dependent. △ Less

Submitted 7 September, 2023; originally announced September 2023.

The DECam Ecliptic Exploration Project (DEEP): I. Survey description, science questions, and technical demonstration

Authors: David E. Trilling, David W. Gerdes, Mario Juric, Chadwick A. Trujillo, Pedro H. Bernardinelli, Kevin J. Napier, Hayden Smotherman, Ryder Strauss, Cesar Fuentes, Matthew J. Holman, Hsing Wen Lin, Larissa Markwardt, Andrew McNeill, Michael Mommert, William J. Oldroyd, Matthew J. Payne, Darin Ragozzine, Andrew S. Rivkin, Hilke Schlichting, Scott S. Sheppard, Fred C. Adams, Colin Orion Chandler

Abstract: We present here the DECam Ecliptic Exploration Project (DEEP), a three year NOAO/NOIRLab Survey that was allocated 46.5 nights to discover and measure the properties of thousands of trans-Neptunian objects (TNOs) to magnitudes as faint as VR~27, corresponding to sizes as small as 20 km diameter. In this paper we present the science goals of this project, the experimental design of our survey, and… ▽ More We present here the DECam Ecliptic Exploration Project (DEEP), a three year NOAO/NOIRLab Survey that was allocated 46.5 nights to discover and measure the properties of thousands of trans-Neptunian objects (TNOs) to magnitudes as faint as VR~27, corresponding to sizes as small as 20 km diameter. In this paper we present the science goals of this project, the experimental design of our survey, and a technical demonstration of our approach. The core of our project is "digital tracking," in which all collected images are combined at a range of motion vectors to detect unknown TNOs that are fainter than the single exposure depth of VR~23 mag. Through this approach we reach a depth that is approximately 2.5 magnitudes fainter than the standard LSST "wide fast deep" nominal survey depth of 24.5 mag. DEEP will more than double the number of known TNOs with observational arcs of 24 hours or more, and increase by a factor of 10 or more the number of known small (<50 km) TNOs. We also describe our ancillary science goals, including measuring the mean shape distribution of very small main belt asteroids, and briefly outline a set of forthcoming papers that present further aspects of and preliminary results from the DEEP program. △ Less

Submitted 6 September, 2023; originally announced September 2023.

Comments: AJ, in press. First in a series of papers

Cluster Cosmology Redux: A Compact Model of the Halo Mass Function

Authors: Cameron E. Norton, Fred C. Adams, August E. Evrard

Abstract: Massive halos hosting groups and clusters of galaxies imprint coherent, arcminute-scale features across the spectrophotometric sky, especially optical-IR clusters of galaxies, distortions in the sub-mm CMB, and extended sources of X-ray emission. Statistical modeling of such features often rely upon the evolving space-time density of dark matter halos -- the halo mass function (HMF) -- as a common… ▽ More Massive halos hosting groups and clusters of galaxies imprint coherent, arcminute-scale features across the spectrophotometric sky, especially optical-IR clusters of galaxies, distortions in the sub-mm CMB, and extended sources of X-ray emission. Statistical modeling of such features often rely upon the evolving space-time density of dark matter halos -- the halo mass function (HMF) -- as a common theoretical ground for cosmological, astrophysical and fundamental physics studies. We propose a compact (eight parameter) representation of the HMF with readily interpretable parameters that stem from polynomial expansions, first in terms of log-mass, then expanding those coefficients similarly in redshift. We demonstrate good ($\sim \! 5\%$) agreement of this form, referred to as the dual-quadratic (DQ-HMF), with Mira-Titan N-body emulator estimates for halo masses above $10^{13.7} h^{-1} {\rm M}_\odot$ over the redshift range $0.1 < z < 1.5$, present best-fit parameters for a Planck 2018 cosmology, and present parameter variation in the $σ_8 - Ω_{\rm m}$ plane. Convolving with a minimal mass-observable relation (MOR) yields closed-form expressions for counts, mean mass, and mass variance of cluster samples characterized by some observable property. Performing information-matrix forecasts of potential parameter constraints from existing and future surveys under different levels of systematic uncertainties, we demonstrate the potential for percent-level constraints on model parameters by an LSST-like optical cluster survey of 300,000 clusters and a richness-mass variance of $0.3^2$. Even better constraints could potentially be achieved by a survey with one-tenth the sample size but with a reduced selection property variance of $0.1^2$. Potential benefits and extensions to the basic MOR parameterization are discussed. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: 18 pages, 9 figures

Avalanches and the Distribution of Reconnection Events in Magnetized Circumstellar Disks

Authors: Marco Fatuzzo, Fred C. Adams, Adina D. Feinstein, Darryl Z. Seligman

Abstract: Cosmic rays produced by young stellar objects can potentially alter the ionization structure, heating budget, chemical composition, and accretion activity in circumstellar disks. The inner edges of these disks are truncated by strong magnetic fields, which can reconnect and produce flaring activity that accelerates cosmic radiation. The resulting cosmic rays can provide a source of ionization and… ▽ More Cosmic rays produced by young stellar objects can potentially alter the ionization structure, heating budget, chemical composition, and accretion activity in circumstellar disks. The inner edges of these disks are truncated by strong magnetic fields, which can reconnect and produce flaring activity that accelerates cosmic radiation. The resulting cosmic rays can provide a source of ionization and produce spallation reactions that alter the composition of planetesimals. This reconnection and particle acceleration are analogous to the physical processes that produce flaring in and heating of stellar coronae. Flaring events on the surface of the Sun exhibit a power-law distribution of energy, reminiscent of those measured for Earthquakes and avalanches. Numerical lattice-reconnection models are capable of reproducing the observed power-law behavior of solar flares under the paradigm of self-organized criticality. One interpretation of these experiments is that the solar corona maintains a nonlinear attractor -- or ``critical'' -- state by balancing energy input via braided magnetic fields and output via reconnection events. Motivated by these results, we generalize the lattice-reconnection formalism for applications in the truncation region of magnetized disks. Our numerical experiments demonstrate that these nonlinear dynamical systems are capable of both attaining and maintaining criticality in the presence of Keplerian shear and other complications. The resulting power-law spectrum of flare energies in the equilibrium attractor state is found to be nearly universal in magnetized disks. This finding indicates that magnetic reconnection and flaring in the inner regions of circumstellar disks occur in a manner similar to activity on stellar surfaces. △ Less

Submitted 12 August, 2023; originally announced August 2023.

Photometric Survey of Neptune's Trojan Asteroids I: The Color Distribution

Authors: Larissa Markwardt, Hsing Wen Lin, David Gerdes, Fred C. Adams

Abstract: In 2018, Jewitt identified the "The Trojan Color Conundrum", namely that Neptune's Trojan asteroids (NTs) had no ultra-red members, unlike the the nearby Kuiper Belt. Since then, numerous ultra-red NTs have been discovered, seemingly resolving this conundrum (Lin et al. 2019; Bolin et al.12 2023). However, it is still unclear whether or not the Kuiper Belt has a color distribution consistent with… ▽ More In 2018, Jewitt identified the "The Trojan Color Conundrum", namely that Neptune's Trojan asteroids (NTs) had no ultra-red members, unlike the the nearby Kuiper Belt. Since then, numerous ultra-red NTs have been discovered, seemingly resolving this conundrum (Lin et al. 2019; Bolin et al.12 2023). However, it is still unclear whether or not the Kuiper Belt has a color distribution consistent with the NT population, as would be expected if it were the source population. In this work, we present a new photometric survey of 15 out of 31 NTs. We utilized the Sloan g'r'i'z' filters on the IMACS f/4 instrument which is mounted on the 6.5m Baade telescope. In this survey, we identify four NTs as being ultra-red using a Principal Component Analysis (PCA). This result brings the ratio of red to ultra-red NTs to 7.75:1, more consistent with the corresponding Trans-Neptunian Object (TNO) ratio of 4-11:1. We also identify three targets as being blue (nearly Solar) in color. Such objects may be C-type surfaces, but we see more of these blue NTs than has been observed in the Kuiper Belt (Seccull et al. 2018). Finally, we show that there are hints of a color-absolute magnitude (H) correlation, with larger H (smaller sized, lower albedo) tending to be more red, but more data is needed to confirm this result. The origin of such a correlation remains an open question which will be addressed by future observations of the surface composition of these targets and their rotational properties. △ Less

Submitted 19 July, 2023; originally announced July 2023.

Comments: 15 pages, 8 figures, accepted to PSJ

The Origin of Universality in the Inner Edges of Planetary Systems

Authors: Konstantin Batygin, Fred C. Adams, Juliette Becker

Abstract: The characteristic orbital period of the inner-most objects within the galactic census of planetary and satellite systems appears to be nearly universal, with $P$ on the order of a few days. This paper presents a theoretical framework that provides a simple explanation for this phenomenon. By considering the interplay between disk accretion, magnetic field generation by convective dynamos, and Kel… ▽ More The characteristic orbital period of the inner-most objects within the galactic census of planetary and satellite systems appears to be nearly universal, with $P$ on the order of a few days. This paper presents a theoretical framework that provides a simple explanation for this phenomenon. By considering the interplay between disk accretion, magnetic field generation by convective dynamos, and Kelvin-Helmholtz contraction, we derive an expression for the magnetospheric truncation radius in astrophysical disks, and find that the corresponding orbital frequency is independent of the mass of the host body. Our analysis demonstrates that this characteristic frequency corresponds to a period of $P\sim3$ days, although intrinsic variations in system parameters are expected to introduce a factor of $\sim2-3$ spread in this result. Standard theory of orbital migration further suggests that planets should stabilize at an orbital period that exceeds disk truncation by a small margin. Cumulatively, our findings predict that the periods of close-in bodies should span $P\sim2-12$ days - a range that is consistent with observations. △ Less

Submitted 14 June, 2023; originally announced June 2023.

Comments: 9 pages, 4 figures, accepted to ApJL

Two Warm Super-Earths Transiting the Nearby M Dwarf TOI-2095

Authors: Elisa V. Quintana, Emily A. Gilbert, Thomas Barclay, Michele L. Silverstein, Joshua E. Schlieder, Ryan Cloutier, Samuel N. Quinn, Joseph E. Rodriguez, Andrew Vanderburg, Benjamin J. Hord, Dana R. Louie, Colby Ostberg, Stephen R. Kane, Kelsey Hoffman, Jason F. Rowe, Giada N. Arney, Prabal Saxena, Taran Richardson, Matthew S. Clement, Nicholas M. Kartvedt, Fred C. Adams, Marcus Alfred, Travis Berger, Allyson Bieryla, Paul Bonney , et al. (33 additional authors not shown)

Abstract: We report the detection and validation of two planets orbiting TOI-2095 (TIC 235678745). The host star is a 3700K M1V dwarf with a high proper motion. The star lies at a distance of 42 pc in a sparsely populated portion of the sky and is bright in the infrared (K=9). With data from 24 Sectors of observation during TESS's Cycles 2 and 4, TOI-2095 exhibits two sets of transits associated with super-… ▽ More We report the detection and validation of two planets orbiting TOI-2095 (TIC 235678745). The host star is a 3700K M1V dwarf with a high proper motion. The star lies at a distance of 42 pc in a sparsely populated portion of the sky and is bright in the infrared (K=9). With data from 24 Sectors of observation during TESS's Cycles 2 and 4, TOI-2095 exhibits two sets of transits associated with super-Earth-sized planets. The planets have orbital periods of 17.7 days and 28.2 days and radii of 1.30 and 1.39 Earth radii, respectively. Archival data, preliminary follow-up observations, and vetting analyses support the planetary interpretation of the detected transit signals. The pair of planets have estimated equilibrium temperatures of approximately 400 K, with stellar insolations of 3.23 and 1.73 times that of Earth, placing them in the Venus zone. The planets also lie in a radius regime signaling the transition between rock-dominated and volatile-rich compositions. They are thus prime targets for follow-up mass measurements to better understand the properties of warm, transition radius planets. The relatively long orbital periods of these two planets provide crucial data that can help shed light on the processes that shape the composition of small planets orbiting M dwarfs. △ Less

Submitted 18 April, 2023; originally announced April 2023.

Comments: Submitted to AAS Journals

The Influence of Tidal Heating on the Habitability of Planets Orbiting White Dwarfs

Authors: Juliette Becker, Darryl Z. Seligman, Fred C. Adams, Marshall J. Styczinski

Abstract: In recent years, there have been a growing number of observations indicating the presence of rocky material in short-period orbits around white dwarfs. In this Letter, we revisit the prospects for habitability around these post-main-sequence star systems. In addition to the typically considered radiative input luminosity, potentially habitable planets around white dwarfs are also subjected to sign… ▽ More In recent years, there have been a growing number of observations indicating the presence of rocky material in short-period orbits around white dwarfs. In this Letter, we revisit the prospects for habitability around these post-main-sequence star systems. In addition to the typically considered radiative input luminosity, potentially habitable planets around white dwarfs are also subjected to significant tidal heating. The combination of these two heating sources can, for a narrow range of planetary properties and orbital parameters, continuously maintain surface temperatures amenable for habitability for planets around white dwarfs over time scales up to 10 Gyr. We show that for a specific locus of orbital parameter space, tidal heating can substantially extend the timescale of continuous habitability for a planet around a white dwarf. △ Less

Submitted 3 March, 2023; originally announced March 2023.

Comments: Accepted to ApJL

Formation History of HD106906 and the Vertical Warping of Debris Disks by an External Inclined Companion

Authors: Nathaniel W. H. Moore, Gongjie Li, Lee Hassenzahl, Erika R. Nesvold, Smadar Naoz, Fred C. Adams

Abstract: HD106906 is a planetary system that hosts a wide-orbit companion, as well as an eccentric and flat debris disk, which hold important constraints on its formation and subsequent evolution. The recent observations of the companion constrain its orbit to be eccentric and inclined relative to the plane of the debris disk. Here, we show that, in the presence of the inclined companion, the debris disk q… ▽ More HD106906 is a planetary system that hosts a wide-orbit companion, as well as an eccentric and flat debris disk, which hold important constraints on its formation and subsequent evolution. The recent observations of the companion constrain its orbit to be eccentric and inclined relative to the plane of the debris disk. Here, we show that, in the presence of the inclined companion, the debris disk quickly ($\lesssim5$ Myr) becomes warped and puffy. This suggests that the current configuration of the system is relatively recent. We explore the possibility that a recent close encounter with a free floating planet could produce a companion with orbital parameters that agree with observations of HD106906b. We find that this scenario is able to recreate the structure of the debris disk while producing a companion in agreement with observation. △ Less

Submitted 28 November, 2022; v1 submitted 23 September, 2022; originally announced September 2022.

Comments: 13 pages, 7 figures. Accepted for publication in ApJ

ARMADA II: Further Detections of Inner Companions to Intermediate Mass Binaries with Micro-Arcsecond Astrometry at CHARA and VLTI

Authors: Tyler Gardner, John D. Monnier, Francis C. Fekel, Jean-Baptiste Le Bouquin, Adam Scovera, Gail Schaefer, Stefan Kraus, Fred C. Adams, Narsireddy Anugu, Jean-Philippe Berger, Theo Ten Brummelaar, Claire L. Davies, Jacob Ennis, Douglas R. Gies, Keith J. C. Johnson, Pierre Kervella, Kaitlin M. Kratter, Aaron Labdon, Cyprien Lanthermann, Johannes Sahlmann, Benjamin R. Setterholm

Abstract: We started a survey with CHARA/MIRC-X and VLTI/GRAVITY to search for low mass companions orbiting individual components of intermediate mass binary systems. With the incredible precision of these instruments, we can detect astrometric "wobbles" from companions down to a few tens of micro-arcseconds. This allows us to detect any previously unseen triple systems in our list of binaries. We present t… ▽ More We started a survey with CHARA/MIRC-X and VLTI/GRAVITY to search for low mass companions orbiting individual components of intermediate mass binary systems. With the incredible precision of these instruments, we can detect astrometric "wobbles" from companions down to a few tens of micro-arcseconds. This allows us to detect any previously unseen triple systems in our list of binaries. We present the orbits of 12 companions around early F to B-type binaries, 9 of which are new detections and 3 of which are first astrometric detections of known RV companions. The masses of these newly detected components range from 0.45-1.3 solar masses. Our orbits constrain these systems to a high astrometric precision, with median residuals to the orbital fit of 20-50 micro-arcseconds in most cases. For 7 of these systems we include newly obtained radial velocity data, which help us to identify the system configuration and to solve for masses of individual components in some cases. Although additional RV measurements are needed to break degeneracy in the mutual inclination, we find that the majority of these inner triples are not well-aligned with the wide binary orbit. This hints that higher mass triples are more misaligned compared to solar and lower mass triples, though a thorough study of survey biases is needed. We show that the ARMADA survey is extremely successful at uncovering previously unseen companions in binaries. This method will be used in upcoming papers to constrain companion demographics in intermediate mass binary systems down to the planetary mass regime. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 25 pages, 14 figures, accepted for publication in AJ

Analytic Approach to the Late Stages of Giant Planet Formation

Authors: Fred C Adams, Konstantin Batygin

Abstract: This paper constructs an analytic description for the late stages of giant planet formation. During this phase of evolution, the planet gains the majority of its final mass through gas accretion at a rapid rate. This work determines the density and velocity fields for material falling onto the central planet and its circumplanetary disk, and finds the corresponding column density of this infalling… ▽ More This paper constructs an analytic description for the late stages of giant planet formation. During this phase of evolution, the planet gains the majority of its final mass through gas accretion at a rapid rate. This work determines the density and velocity fields for material falling onto the central planet and its circumplanetary disk, and finds the corresponding column density of this infalling envelope. We derive a steady-state solution for the surface density of the disk as a function of its viscosity (including the limiting case where no disk accretion occurs). Planetary magnetic fields truncate the inner edge of the disk and determine the boundary conditions for mass accretion onto the planet from both direct infall and from the disk. The properties of the forming planet and its circumplanetary disk are determined, including the luminosity contributions from infall onto the planet and disk surfaces, and from disk viscosity. The radiative signature of the planet formation process is explored using a quasi-spherical treatment of the emergent spectral energy distributions. The analytic solutions developed herein show how the protoplanet properties (envelope density distribution, velocity field, column density, disk surface density, luminosity, and radiative signatures) vary with input parameters (instantaneous mass, orbital location, accretion rate, and planetary magnetic field strength). △ Less

Submitted 20 June, 2022; originally announced June 2022.

Comments: 52 pages, 9 figures, accepted to The Astrophysical Journal

Gemini-LIGHTS: Herbig Ae/Be and massive T-Tauri protoplanetary disks imaged with Gemini Planet Imager

Authors: Evan A. Rich, John D. Monnier, Alicia Aarnio, Anna S. E. Laws, Benjamin R. Setterholm, David J. Wilner, Nuria Calvet, Tim Harries, Chris Miller, Claire L. Davies, Fred C. Adams, Sean M. Andrews, Jaehan Bae, Catherine Espaillat, Alexandra Z. Greenbaum, Sasha Hinkley, Stefan Kraus, Lee Hartmann, Andrea Isella, Melissa McClure, Rebecca Oppenheimer, Laura M. Pérez, Zhaohuan Zhu

Abstract: We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of… ▽ More We present the complete sample of protoplanetary disks from the Gemini- Large Imaging with GPI Herbig/T-tauri Survey (Gemini-LIGHTS) which observed bright Herbig Ae/Be stars and T-Tauri stars in near-infrared polarized light to search for signatures of disk evolution and ongoing planet formation. The 44 targets were chosen based on their near- and mid-infrared colors, with roughly equal numbers of transitional, pre-transitional, and full disks. Our approach explicitly did not favor well-known, "famous" disks or those observed by ALMA, resulting in a less-biased sample suitable to probe the major stages of disk evolution during planet formation. Our optimized data reduction allowed polarized flux as low as 0.002% of the stellar light to be detected, and we report polarized scattered light around 80% of our targets. We detected point-like companions for 47% of the targets, including 3 brown dwarfs (2 confirmed, 1 new), and a new super-Jupiter mass candidate around V1295 Aql. We searched for correlations between the polarized flux and system parameters, finding a few clear trends: presence of a companion drastically reduces the polarized flux levels, far-IR excess correlates with polarized flux for non-binary systems, and systems hosting disks with ring structures have stellar masses $<$ 3 Msun. Our sample also included four hot, dusty "FS CMa" systems and we detected large-scale ($>100$ au) scattered light around each, signs of extreme youth for these enigmatic systems. Science-ready images are publicly available through multiple distribution channels using a new FITS file standard jointly developed with members of the VLT/SPHERE team. △ Less

Submitted 12 June, 2022; originally announced June 2022.

Comments: 51 pages, 31 figures, 7 tables, accepted to AJ

Transfer of Rocks between Planetary Systems: Panspermia Revisited

Authors: Fred C Adams, Kevin J Napier

Abstract: Motivated by the recent discovery of interstellar objects passing through the solar system, and by recent developments in dynamical simulations, this paper reconsiders the likelihood for life bearing rocks to be transferred from one planetary system to another. The astronomical aspects of this lithopanspermia process can now be estimated, including the cross sections for rock capture, the velocity… ▽ More Motivated by the recent discovery of interstellar objects passing through the solar system, and by recent developments in dynamical simulations, this paper reconsiders the likelihood for life bearing rocks to be transferred from one planetary system to another. The astronomical aspects of this lithopanspermia process can now be estimated, including the cross sections for rock capture, the velocity distributions of rocky ejecta, the survival times for captured objects, and the dynamics of the solar system in both its birth cluster and in the field. The remaining uncertainties are primarily biological, i.e., the probability of life developing on a planet, the time required for such an event, and the efficiency with which life becomes seeded in a new environment. Using current estimates for the input quantities, we find that the transfer rates are enhanced in the birth cluster, but the resulting odds for success are too low for panspermia to be a likely occurrence. In contrast, the expected inventory of alien rocks in the solar system is predicted to be substantial (where the vast majority of such bodies are not biologically active and do not interact with Earth). △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 18 pages, 3 figures, to appear in Astrobiology

Possible Solution to the Triple Alpha Fine-Tuning Problem: Spallation Reactions during Planet Formation

Authors: Fred C Adams

Abstract: Carbon is produced during the helium burning phase of sufficiently massive stars through the triple alpha process. The $0^+$ energy level of the carbon nucleus allows for resonant nuclear reactions, which act to greatly increase the carbon yields compared to the non-resonant case. Many authors have argued that small changes to the energy level of this resonance would lead to a significantly lower… ▽ More Carbon is produced during the helium burning phase of sufficiently massive stars through the triple alpha process. The $0^+$ energy level of the carbon nucleus allows for resonant nuclear reactions, which act to greatly increase the carbon yields compared to the non-resonant case. Many authors have argued that small changes to the energy level of this resonance would lead to a significantly lower carbon abundance in the universe, and this sensitivity is often considered an example of fine-tuning. By considering spallation reactions occuring during the process of planet formation, this paper presents a partial solution to this triple alpha fine-tuning problem. Young stellar objects generate substantial luminosities of particle radiation (cosmic rays) that can drive nuclear reactions through spallation. If the standard triple alpha process is inoperative, stars tend to synthesize oxygen (and other alpha elements) rather than carbon. Cosmic rays can interact with oxygen nuclei to produce carbon while planets are forming. The resulting carbon abundances are significant, but much smaller than those observed in our universe. However, for a range of conditions -- as delineated herein -- spallation reactions can result in carbon-to-oxygen ratios roughly comparable to those found on Earth and thereby obviate the triple alpha fine-tuning problem. △ Less

Submitted 7 May, 2022; originally announced May 2022.

Comments: 25 pages, 2 figures, accepted for publication in Astroparticle Physics

Inferring Late Stage Enrichment of Exoplanet Atmospheres from Observed Interstellar Comets

Authors: Darryl Z. Seligman, Juliette Becker, Fred C. Adams, Adina D. Feinstein, Leslie A. Rogers

Abstract: The discovery of the first two interstellar objects implies that, on average, every star contributes a substantial amount of material to the galactic population by ejecting such bodies from the host system. Since scattering is a chaotic process, a comparable amount of material should be injected into the inner regions of each system that ejects comets. For comets that are transported inwards and i… ▽ More The discovery of the first two interstellar objects implies that, on average, every star contributes a substantial amount of material to the galactic population by ejecting such bodies from the host system. Since scattering is a chaotic process, a comparable amount of material should be injected into the inner regions of each system that ejects comets. For comets that are transported inwards and interact with planets, this Letter estimates the fraction of material that is accreted or outward-scattered as a function of planetary masses and orbital parameters. These calculations indicate that planets with escape velocities smaller than their current day orbital velocities will efficiently accrete comets. We estimate the accretion efficiency for members of the current census of extrasolar planets, and find that planetary populations including but not limited to hot and warm Jupiters, sub-Neptunes and super-Earths can efficiently capture incoming comets. This cometary enrichment may have important ramifications for post-formation atmospheric composition and chemistry. As a result, future detections and compositional measurements of interstellar comets will provide direct measurements of material that potentially enriched a sub-population of the extrasolar planets. Finally, we estimate the efficiency of this enrichment mechanism for extrasolar planets that will be observed with the $\textit{James Webb Space Telescope}$ (JWST). With JWST currently operational and these observations imminently forthcoming, it is of critical importance to investigate how enrichment from interstellar comet analogues may affect the interpretations of exoplanet atmospheric compositions. △ Less

Submitted 15 June, 2022; v1 submitted 26 April, 2022; originally announced April 2022.

Comments: Accepted for publication at ApJL, 9 pages, 3 figures, preprint for reference at Exoplanets IV Program Number 405.03

A Collision Mechanism for the Removal of Earth's Trojan Asteroids

Authors: Kevin J. Napier, Larissa Markwardt, Fred C. Adams, David W. Gerdes, Hsing Wen Lin

Abstract: Due to their strong resonances with their host planet, Trojan asteroids can remain in stable orbits for billions of years. As a result, they are powerful probes for constraining the dynamical and chemical history of the solar system. Although we have detected thousands of Jupiter Trojans and dozens of Neptune Trojans, there are currently no known long-term stable Earth Trojans. Dynamical simulatio… ▽ More Due to their strong resonances with their host planet, Trojan asteroids can remain in stable orbits for billions of years. As a result, they are powerful probes for constraining the dynamical and chemical history of the solar system. Although we have detected thousands of Jupiter Trojans and dozens of Neptune Trojans, there are currently no known long-term stable Earth Trojans. Dynamical simulations show that the parameter space for stable Earth Trojans in substantial, so their apparent absence poses a mystery. This work uses a large ensemble of N-body simulations to explore how the Trojan population dynamically responds if Earth suffers large collisions, such as those thought to have occurred to form the Moon and/or to have given Earth its Late Veneer. We show that such collisions can be highly disruptive to the primordial Trojan population, and could have eliminated it altogether. More specifically, if Earth acquired the final 1\% of its mass through ${\cal O}(10)$ collisions, then only $\sim1\%$ of the previously bound Trojan population would remain. △ Less

Submitted 21 April, 2022; originally announced April 2022.

Comments: Accepted for publication in the Planetary Science Journal

Architectures of Compact Multi-planet Systems: Diversity and Uniformity

Authors: Lauren M. Weiss, Sarah C. Millholland, Erik A. Petigura, Fred C. Adams, Konstantin Batygin, Anthony M. Bloch, Christoph Mordasini

Abstract: One of the most important developments in exoplanet science in the past decade is the discovery of multi-planet systems with sub-Neptune-sized planets interior to 1~AU. This chapter explores the architectures of these planetary systems, which often display a remarkable degree of uniformity: the planets have nearly equal sizes, regular orbital spacing, low eccentricities, and small mutual inclinati… ▽ More One of the most important developments in exoplanet science in the past decade is the discovery of multi-planet systems with sub-Neptune-sized planets interior to 1~AU. This chapter explores the architectures of these planetary systems, which often display a remarkable degree of uniformity: the planets have nearly equal sizes, regular orbital spacing, low eccentricities, and small mutual inclinations. This uniformity stands in sharp contrast to the diverse nature of the exoplanet sample considered as a whole (as well as our inner solar system). We begin with a critical review of the observations -- including possible biases -- and find that these peas-in-a-pod planetary systems are apparently a common outcome of the planet formation process. Modest departures from exact uniformity suggest additional patterns, such as the planet mass slowly increasing with semi-major axis. The star formation process naturally produces circumstellar disks with the properties required to produce these planetary systems, although the solid material must move inward from its initial location. We discuss primary modes of planetary assembly, the role of orbital migration, and post-nebular atmospheric loss. Mature planetary systems are found to be near their minimum energy (tidal equilibrium) configurations; this finding provides a partial explanation for their observed properties and indicates that efficient energy dissipation must occur. Finally, we consider population synthesis models and show that peas-in-a-pod patterns emerge with reasonable choices for the input parameters. Nonetheless, interesting observational and theoretical challenges remain in order to understand how these surprisingly organized planetary systems arise from the disorder of their formation processes. △ Less

Submitted 18 March, 2022; originally announced March 2022.

Comments: Chapter for publication in Protostars and Planets VII, 42 pages, 17 Figures

Theoretical and Observational Evidence for Coriolis Effects in Coronal Magnetic Fields Via Direct Current Driven Flaring Events

Authors: Darryl Z. Seligman, Leslie A. Rogers, Adina D. Feinstein, Mark R. Krumholz, James R. Beattie, Christoph Federrath, Fred C. Adams, Marco Fatuzzo, Maximilian N. Günther

Abstract: All stars produce explosive surface events such as flares and coronal mass ejections. These events are driven by the release of energy stored in coronal magnetic fields, generated by the stellar dynamo. However, it remains unclear if the energy deposition in the magnetic fields is driven by direct or alternating currents. Recently, we presented observational measurements of the flare intensity dis… ▽ More All stars produce explosive surface events such as flares and coronal mass ejections. These events are driven by the release of energy stored in coronal magnetic fields, generated by the stellar dynamo. However, it remains unclear if the energy deposition in the magnetic fields is driven by direct or alternating currents. Recently, we presented observational measurements of the flare intensity distributions for a sample of $\sim10^5$ stars across the main sequence observed by $\textit{TESS}$, all of which exhibited power-law distributions similar to those observed in the Sun, albeit with varying slopes. Here we investigate the mechanisms required to produce such a distribution of flaring events via direct current energy deposition, in which coronal magnetic fields braid, reconnect, and produce flares. We adopt a topological model for this process which produces a power-law distribution of energetic flaring events. We expand this model to include the Coriolis effect, which we demonstrate produces a shallower distribution of flare energies in stars that rotate more rapidly (corresponding to a weaker decline in occurrence rates toward increasing flare energies). We present tentative evidence for the predicted rotation-power-law index correlation in the observations. We advocate for future observations of stellar flares that would improve our measurements of the power-law exponents, and yield key insights into the underlying dynamo mechanisms that underpin the self-similar flare intensity distributions. △ Less

Submitted 5 March, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: 19 pages, 11 figures, accepted for publication in ApJ

Theoretical Distributions of Short-Lived Radionuclides for Star Formation in Molecular Clouds

Authors: M. Fatuzzo, F. C. Adams

Abstract: Short-lived radioactive nulcei (half-life $τ_{1/2}\sim1$ Myr) influence the formation of stars and planetary systems by providing sources of heating and ionization. Whereas many previous studies have focused on the possible nuclear enrichment of our own Solar System, the goal of this paper is to estimate the distributions of short-lived radionuclides (SLRs) for the entire population of stars formi… ▽ More Short-lived radioactive nulcei (half-life $τ_{1/2}\sim1$ Myr) influence the formation of stars and planetary systems by providing sources of heating and ionization. Whereas many previous studies have focused on the possible nuclear enrichment of our own Solar System, the goal of this paper is to estimate the distributions of short-lived radionuclides (SLRs) for the entire population of stars forming within a molecular cloud. Here we focus on the nuclear species $^{60}$Fe and $^{26}$Al, which have the largest impact due to their relatively high abundances. We construct molecular cloud models and include nuclear contributions from both supernovae and stellar winds. The resulting distributions of SLRs are time dependent with widths of $\sim3$ orders of magnitude and mass fractions $ρ_{\scriptstyle SLR}/ρ_\ast\sim10^{-11}-10^{-8}$. Over the range of scenarios explored herein, the SLR distributions show only modest variations with the choice of cloud structure (fractal dimension), star formation history, and cluster distribution. The most important variation arises from the diffusion length scale for the transport of SLRs within the cloud. The expected SLR distributions are wide enough to include values inferred for the abundances in our Solar System, although most of the stars are predicted to have smaller enrichment levels. In addition, the ratio of $^{60}$Fe/$^{26}$Al is predicted to be greater than unity, on average, in contrast to Solar System results. One explanation for this finding is the presence of an additional source for the $^{26}$Al isotope. △ Less

Submitted 18 November, 2021; originally announced November 2021.

Comments: accepted for publication in ApJ

On the Fate of Interstellar Objects Captured by our Solar System

Authors: Kevin J Napier, Fred C Adams, Konstantin Batygin

Abstract: With the recent discoveries of interstellar objects `Oumuamua and Borisov traversing the solar system, understanding the dynamics of interstellar objects is more pressing than ever. These detections have highlighted the possibility that captured interstellar material could be trapped in our solar system. The first step in rigorously investigating this question is to calculate a capture cross secti… ▽ More With the recent discoveries of interstellar objects `Oumuamua and Borisov traversing the solar system, understanding the dynamics of interstellar objects is more pressing than ever. These detections have highlighted the possibility that captured interstellar material could be trapped in our solar system. The first step in rigorously investigating this question is to calculate a capture cross section for interstellar objects as a function of hyperbolic excess velocity, which can be convolved with any velocity dispersion to compute a capture rate (Napier et. al. 2021). Although the cross section provides the first step toward calculating the mass of alien rocks residing in our solar system, we also need to know the lifetime of captured objects. We use an ensemble of N-body simulations to characterize a dynamical lifetime for captured interstellar objects and determines the fraction of surviving objects as a function of time (since capture). We also illuminate the primary effects driving their secular evolution. Finally, we use the resulting dynamical lifetime function to estimate the current inventory of captured interstellar material in the solar system. We find that capture from the field yields a steady state mass of only $\sim 10^{-13} M_{\oplus}$, whereas the mass remaining from capture events in the birth cluster is roughly $10^{-9} M_{\oplus}$. △ Less

Submitted 22 September, 2021; originally announced September 2021.

Comments: Accepted for publication in the Planetary Science Journal. 12 pages, 5 figures

Testing Self-Organized Criticality Across the Main Sequence using Stellar Flares from TESS

Authors: Adina D. Feinstein, Darryl Z. Seligman, Maximilian N. Günther, Fred C. Adams

Abstract: Self-organized criticality describes a class of dynamical systems that maintain themselves in an attractor state with no intrinsic length or time scale. Fundamentally, this theoretical construct requires a mechanism for instability that may trigger additional instabilities locally via dissipative processes. This concept has been invoked to explain nonlinear dynamical phenomena such as featureless… ▽ More Self-organized criticality describes a class of dynamical systems that maintain themselves in an attractor state with no intrinsic length or time scale. Fundamentally, this theoretical construct requires a mechanism for instability that may trigger additional instabilities locally via dissipative processes. This concept has been invoked to explain nonlinear dynamical phenomena such as featureless energy spectra that have been observed empirically for earthquakes, avalanches, and solar flares. If this interpretation proves correct, it implies that the solar coronal magnetic field maintains itself in a critical state via a delicate balance between the dynamo-driven injection of magnetic energy and the release of that energy via flaring events. All-sky high-cadence surveys like the Transiting Exoplanet Survey Satellite (TESS) provide the necessary data to compare the energy distribution of flaring events in stars of different spectral types to that observed in the Sun. We identified $\sim 10^6$ flaring events on $\sim 10^5$ stars observed by TESS at 2-minute cadence. By fitting the flare frequency distribution for different mass bins, we find that all main sequence stars exhibit distributions of flaring events similar to that observed in the Sun, independent of their mass or age. This may suggest that stars universally maintain a critical state in their coronal topologies via magnetic reconnection events. If this interpretation proves correct, we may be able to infer properties of magnetic fields, interior structure, and dynamo mechanisms for stars that are otherwise unresolved point sources. △ Less

Submitted 12 January, 2022; v1 submitted 14 September, 2021; originally announced September 2021.

Comments: 6 pages, 3 figures, Accepted to ApJL

Radioactive Planet Formation

Authors: Fred C Adams

Abstract: Young stellar objects are observed to have large X-ray fluxes and are thought to produce commensurate luminosities in energetic particles (cosmic rays). This particle radiation, in turn, can synthesize short-lived radioactive nuclei through spallation. With a focus on $^{26}$Al, this paper estimates the expected abundances of radioactive nulcei produced by spallation during the epoch of planet for… ▽ More Young stellar objects are observed to have large X-ray fluxes and are thought to produce commensurate luminosities in energetic particles (cosmic rays). This particle radiation, in turn, can synthesize short-lived radioactive nuclei through spallation. With a focus on $^{26}$Al, this paper estimates the expected abundances of radioactive nulcei produced by spallation during the epoch of planet formation. In this model, cosmic rays are accelerated near the inner truncation radii of circumstellar disks, $r_{\scriptstyle X}\approx0.1$ AU, where intense magnetic activity takes place. For planets forming in this region, radioactive abundances can be enhanced over the values inferred for the early solar system (from meteoritic measurements) by factors of $\sim10-20$. These short-lived radioactive nuclei influence the process of planet formation and the properties of planets in several ways. The minimum size required for planetesimals to become fully molten decreases with increasing levels of radioactive enrichment, and such melting leads to loss of volatile components including water. Planets produced with an enhanced radioactive inventory have significant internal luminosity which can be comparable to that provided by the host star; this additional heating affects both atmospheric mass loss and chemical composition. Finally, the habitable zone of red dwarf stars is coincident with the magnetic reconnection region, so that planets forming at those locations will experience maximum exposure to particle radiation, and subsequent depletion of volatiles. △ Less

Submitted 7 July, 2021; originally announced July 2021.

Comments: 29 pages, 4 figures, accepted to ApJ

Investigating the Relative Gas and Small Dust Grain Surface Heights in Protoplanetary Disks

Authors: Evan A. Rich, Richard Teague, John D. Monnier, Claire L. Davies, Arthur Bosman, Tim J. Harries, Nuria Calvet, Fred C. Adams, Zhaohuan Zhu

Abstract: Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($\sim$1 $μm$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $β$ of the small dust grain sc… ▽ More Dust evolution in protoplanetary disks from small dust grains to pebbles is key to the planet formation process. The gas in protoplanetary disks should influence the vertical distribution of small dust grains ($\sim$1 $μm$) in the disk.Utilizing archival near-infrared polarized light and millimeter observations, we can measure the scale height and the flare parameter $β$ of the small dust grain scattering surface and $^{12}$CO gas emission surface for three protoplanetary disks IM Lup, HD 163296, and HD 97048 (CU Cha). For two systems, IM Lup and HD 163296, the $^{12}$CO gas and small dust grains at small radii from the star have similar heights but at larger radii ($>$100 au) the dust grain scattering surface height is lower than the $^{12}$CO gas emission surface height. In the case of HD 97048, the small dust grain scattering surface has similar heights to the $^{12}$CO gas emission surface at all radii. We ran a protoplanetary disk radiative transfer model of a generic protoplanetary disk with TORUS and showed that there is no difference between the observed scattering surface and $^{12}$CO emission surface. We also performed analytical modeling of the system and found that gas-to-dust ratios larger than 100 could explain the observed difference in IM Lup and HD 163296. This is the first direct comparison of observations of gas and small dust grain heights distribution in protoplanetary disks. Future observations of gas emission and near-infrared scattered light instruments are needed to look for similar trends in other protoplanetary disks. △ Less

Submitted 15 April, 2021; originally announced April 2021.

Comments: Accepted ApJ, 13 pages, 3 figures, 4 tables

Effects of Bound Diprotons and Enhanced Nuclear Reaction Rates on Stellar Evolution

Authors: Fred C. Adams, Alex R. Howe, Evan Grohs, George M. Fuller

Abstract: Deuterium represents the only bound isotope in the universe with atomic mass number $A=2$. Motivated by the possibility of other universes, where the strong force could be stronger, this paper considers the effects of bound diprotons and dineutrons on stars. We find that the existence of additional stable nuclei with $A=2$ has relatively modest effects on the universe. Previous work indicates that… ▽ More Deuterium represents the only bound isotope in the universe with atomic mass number $A=2$. Motivated by the possibility of other universes, where the strong force could be stronger, this paper considers the effects of bound diprotons and dineutrons on stars. We find that the existence of additional stable nuclei with $A=2$ has relatively modest effects on the universe. Previous work indicates that Big Bang Nucleosynthesis (BBN) produces more deuterium, but does not lead to catastrophic heavy element production. This paper revisits BBN considerations and confirms that the universe is left with an ample supply of hydrogen and other light nuclei for typical cosmological parameters. Using the $MESA$ numerical package, we carry out stellar evolution calculations for universes with stable diprotons, with nuclear cross sections enhanced by large factors $X$. This work focuses on $X=10^{15}-10^{18}$, but explores the wider range $X$ = $10^{-3}-10^{18}$. For a given stellar mass, the presence of stable diprotons leads to somewhat brighter stars, with the radii and photospheric temperatures roughly comparable to thoese of red giants. The central temperature decreases from the characteristic value of $T_c\approx1.5\times10^7$ K for hydrogen burning down to the value of $T_c\approx10^6$ K characteristic of deuterium burning. The stellar lifetimes are smaller for a given mass, but with the extended possible mass range, the smallest stars live for trillions of years, far longer than the current cosmic age. Finally, the enhanced cross sections allow for small, partially degenerate objects with mass $M_\ast=1-10M_J$ to produce significant steady-state luminosity and thereby function as stars. △ Less

Submitted 29 March, 2021; originally announced March 2021.

Comments: 54 pages, 14 figures, accepted to Astroparticle Physics

Evidence for a Cosmic Ray Gradient in the IM Lup Protoplanetary Disk

Authors: Richard A. Seifert, L. Ilsedore. Cleeves, Fred C. Adams, Zhi-Yun Li

Abstract: Protoplanetary disk evolution is strongly impacted by ionization from the central star and local environment, which collectively have been shown to drive chemical complexity and are expected to impact the transport of disk material. Nonetheless, ionization remains a poorly constrained input to many detailed modeling efforts. We use new and archival ALMA observations of N$_2$H$^+$ 3--2 and H… ▽ More Protoplanetary disk evolution is strongly impacted by ionization from the central star and local environment, which collectively have been shown to drive chemical complexity and are expected to impact the transport of disk material. Nonetheless, ionization remains a poorly constrained input to many detailed modeling efforts. We use new and archival ALMA observations of N$_2$H$^+$ 3--2 and H$^{13}$CO$^+$ 3--2 to derive the first observationally-motivated ionization model for the IM Lup protoplanetary disk. Incorporating ionization from multiple internal and external sources, we model N$_2$H$^+$ and H$^{13}$CO$^+$ abundances under varying ionization environments, and compare these directly to the imaged ALMA observations by performing non-LTE radiative transfer, visibility sampling, and imaging. We find that the observations are best reproduced using a radially increasing cosmic ray (CR) gradient, with low CR ionization in the inner disk, high CR ionization in the outer disk, and a transition at $\sim 80 - 100$ au. This location is approximately coincident with the edge of spiral structure identified in millimeter emission. We also find that IM Lup shows evidence for enhanced UV-driven formation of HCO$^+$, which we attribute to the disk's high flaring angle. In summary, IM Lup represents the first protoplanetary disk with observational evidence for a CR gradient, which may have important implications for IM Lup's on-going evolution, especially given the disk's young age and large size. △ Less

Submitted 19 March, 2021; originally announced March 2021.

Comments: 19 pages, 13 figures, accepted to ApJ

On the Capture of Interstellar Objects by our Solar System

Authors: K. J. Napier, F. C. Adams, K. Batygin

Abstract: Motivated by recent visits from interstellar comets, along with continuing discoveries of minor bodies in orbit of the Sun, this paper studies the capture of objects on initially hyperbolic orbits by our solar system. Using an ensemble of $\sim500$ million numerical experiments, this work generalizes previous treatments by calculating the capture cross section as a function of asymptotic speed. Th… ▽ More Motivated by recent visits from interstellar comets, along with continuing discoveries of minor bodies in orbit of the Sun, this paper studies the capture of objects on initially hyperbolic orbits by our solar system. Using an ensemble of $\sim500$ million numerical experiments, this work generalizes previous treatments by calculating the capture cross section as a function of asymptotic speed. The resulting velocity-dependent cross section can then be convolved with any distribution of relative speeds to determine the capture rate for incoming bodies. This convolution is carried out for the usual Maxwellian distribution, as well as the velocity distribution expected for rocky debris ejected from planetary systems. We also construct an analytic description of the capture process that provides an explanation for the functional form of the capture cross section in both the high velocity and low velocity limits. △ Less

Submitted 16 February, 2021; originally announced February 2021.

Comments: 18 pages, 5 figures. Accepted to the Planetary Science Journal

A Theoretical Framework for the Mass Distribution of Gas Giant Planets forming through the Core Accretion Paradigm

Authors: Fred C Adams, Michael R Meyer, Arthur D Adams

Abstract: This paper constructs a theoretical framework for calculating the distribution of masses for gas giant planets forming via the core accretion paradigm. Starting with known properties of circumstellar disks, we present models for the planetary mass distribution over the range $0.1M_J < M_{\rm p} < 10M_J$. If the circumstellar disk lifetime is solely responsible for the end of planetary mass accreti… ▽ More This paper constructs a theoretical framework for calculating the distribution of masses for gas giant planets forming via the core accretion paradigm. Starting with known properties of circumstellar disks, we present models for the planetary mass distribution over the range $0.1M_J < M_{\rm p} < 10M_J$. If the circumstellar disk lifetime is solely responsible for the end of planetary mass accretion, the observed (nearly) exponential distribution of disk lifetime would imprint an exponential fall-off in the planetary mass function. This result is in apparent conflict with observations, which suggest that the mass distribution has a (nearly) power-law form $dF/dM_{\rm p}\sim M_{\rm p}^{-p}$, with index $p\approx1.3$, over the relevant planetary mass range (and for stellar masses $\sim0.5-2M_\odot$). The mass accretion rate onto the planet depends on the fraction of the (circumstellar) disk accretion flow that enters the Hill sphere, and on the efficiency with which the planet captures the incoming material. Models for the planetary mass function that include distributions for these efficiencies, with uninformed priors, can produce nearly power-law behavior, consistent with current observations. The disk lifetimes, accretion rates, and other input parameters depend on the mass of the host star. We show how these variations lead to different forms for the planetary mass function for different stellar masses. Compared to stars with masses $M_\ast$ = $0.5-2M_\odot$, stars with smaller masses are predicted to have a steeper planetary mass function (fewer large planets). △ Less

Submitted 17 January, 2021; originally announced January 2021.

Comments: 29 pages, 4 figures, accepted to ApJ

ARMADA I: Triple Companions Detected in B-Type Binaries alpha Del and nu Gem

Authors: Tyler Gardner, John D. Monnier, Francis C. Fekel, Gail Schaefer, Keith J. C. Johnson, Jean-Baptiste Le Bouquin, Stefan Kraus, Narsireddy Anugu, Benjamin R. Setterholm, Aaron Labdon, Claire L. Davies, Cyprien Lanthermann, Jacob Ennis, Michael Ireland, Kaitlin M. Kratter, Theo Ten Brummelaar, Judit Sturmann, Laszlo Sturmann, Chris Farrington, Douglas R. Gies, Robert Klement, Fred C. Adams

Abstract: Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ~10 micro-arcsecond precision. This precision makes it possible to detect "wobbles" in the binary motion due to the gravitational pull from additional short period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the MIRC-X instrume… ▽ More Ground-based optical long-baseline interferometry has the power to measure the orbits of close binary systems at ~10 micro-arcsecond precision. This precision makes it possible to detect "wobbles" in the binary motion due to the gravitational pull from additional short period companions. We started the ARrangement for Micro-Arcsecond Differential Astrometry (ARMADA) survey with the MIRC-X instrument at the CHARA array for the purpose of detecting giant planets and stellar companions orbiting individual stars in binary systems. We describe our observations for the survey, and introduce the wavelength calibration scheme that delivers precision at the tens of micro-arcseconds level for <0.2 arcsecond binaries. We test our instrument performance on a known triple system kappa Peg, and show that our survey is delivering a factor of 10 better precision than previous similar surveys. We present astrometric detections of tertiary components to two B-type binaries: a 30-day companion to alpha Del, and a 50-day companion to nu Gem. We also collected radial velocity data for alpha Del with the Tennessee State University Automated Spectroscopic Telescope at Fairborn Observatory. We are able to measure the orbits and masses of all three components in these systems. We find that the previously published RV orbit for the inner pair of nu Gem is not consistent with our visual orbit. The precision achieved for these orbits suggests that our ARMADA survey will be successful at discovering new compact triple systems to A/B-type binary systems, leading to better statistics of hierarchical system architectures and formation history. △ Less

Submitted 1 December, 2020; originally announced December 2020.

Comments: 25 pages, 15 figures, accepted for publication in AJ

The Resonance Hopping Effect in the Neptune-Planet Nine System

Authors: Tali Khain, Juliette Becker, Fred C. Adams

Abstract: The observed physical clustering of the orbits of small bodies in the distant Kuiper Belt (TNOs) has recently prompted the prediction of an additional planet in the outer solar system. Since the initial posing of the hypothesis, the effects of Planet Nine on the dynamics of the main cluster of TNOs - the objects anti-aligned with its orbit - have been well-studied. In particular, numerical simulat… ▽ More The observed physical clustering of the orbits of small bodies in the distant Kuiper Belt (TNOs) has recently prompted the prediction of an additional planet in the outer solar system. Since the initial posing of the hypothesis, the effects of Planet Nine on the dynamics of the main cluster of TNOs - the objects anti-aligned with its orbit - have been well-studied. In particular, numerical simulations have revealed a fascinating phenomenon, referred to as "resonance hopping", in which these objects abruptly transition between different mean-motion commensurabilities with Planet Nine. In this work, we explore this effect in greater detail, with the goal of understanding what mechanism prompts the hopping events to occur. In the process, we elucidate the often underestimated role of Neptune scattering interactions, which leads to diffusion in the semi-major axes of these distant TNOs. In addition, we demonstrate that although some resonant interactions with Planet Nine do occur, the anti-aligned objects are able to survive without the resonances, confirming that the dynamics of the TNOs are predominantly driven by secular, rather than resonant, interactions with Planet Nine. △ Less

Submitted 5 October, 2020; originally announced October 2020.

Comments: accepted to PASP

The Origin of Systems of Tightly Packed Inner Planets with Misaligned, Ultra-Short-Period Companions

Authors: Juliette Becker, Konstantin Batygin, Daniel C. Fabrycky, Fred C. Adams, Gongjie Li, Andrew Vanderburg, Joseph E. Rodriguez

Abstract: Ultra-short period planets provide a window into the inner edge of the parameter space occupied by planetary orbits. In one particularly intriguing class of multi-planet systems, the ultra-short period planet is flanked by short-period companions, and the outer planets occupy a discernibly distinct dynamical state. In the observational database, this phenomenon is represented by a small number of… ▽ More Ultra-short period planets provide a window into the inner edge of the parameter space occupied by planetary orbits. In one particularly intriguing class of multi-planet systems, the ultra-short period planet is flanked by short-period companions, and the outer planets occupy a discernibly distinct dynamical state. In the observational database, this phenomenon is represented by a small number of stars hosting systems of tightly packed co-planar planets as well as an ultra-short period planet, whose orbit of is misaligned relative to the mutual plane of the former. In this work, we explore two different mechanisms that can produce an ultra-short period planet that is misaligned with the rest of its compact planetary system: natural decoupling between the inner and outer system via the stellar quadrupole moment, and decoupling forced by an external companion with finely-tuned orbital parameters. These two processes operate with different timescales, and can thus occur simultaneously. In this work, we use the K2-266 system as an illustrative example to elucidate the dynamics of these two processes, and highlight the types of constraints that may arise regarding the dynamical histories of systems hosting ultra-short period planets. △ Less

Submitted 22 September, 2020; originally announced September 2020.

Comments: Accepted to AJ

The Great Inequality and the Dynamical Disintegration of the Outer Solar System

Authors: Jon K. Zink, Konstantin Batygin, Fred C. Adams

Abstract: Using an ensemble of N-body simulations, this paper considers the fate of the outer gas giants (Jupiter, Saturn, Uranus, and Neptune) after the Sun leaves the main sequence and completes its stellar evolution. Due to solar mass-loss -- which is expected to remove roughly half of the star's mass -- the orbits of the giant planets expand. This adiabatic process maintains the orbital period ratios, b… ▽ More Using an ensemble of N-body simulations, this paper considers the fate of the outer gas giants (Jupiter, Saturn, Uranus, and Neptune) after the Sun leaves the main sequence and completes its stellar evolution. Due to solar mass-loss -- which is expected to remove roughly half of the star's mass -- the orbits of the giant planets expand. This adiabatic process maintains the orbital period ratios, but the mutual interactions between planets and the width of mean-motion resonances (MMR) increase, leading to the capture of Jupiter and Saturn into a stable 5:2 resonant configuration. The expanded orbits, coupled with the large-amplitude librations of the critical MMR angle, make the system more susceptible to perturbations from stellar flyby interactions. Accordingly, within about 30 Gyr, stellar encounters perturb the planets onto the chaotic sub-domain of the 5:2 resonance, triggering a large-scale instability, which culminates in the ejections of all but one planet over the subsequent $\sim10$ Gyr. After an additional $\sim50$ Gyr, a close stellar encounter (with a perihelion distance less than $\sim200$ AU) liberates the final planet. Through this sequence of events, the characteristic timescale over which the solar system will be completely dissolved is roughly 100 Gyr. Our analysis thus indicates that the expected dynamical lifetime of the solar system is much longer than the current age of the universe, but is significantly shorter than previous estimates. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 12 pages, 5 figures; accepted for publication in AJ

A Giant Planet Candidate Transiting a White Dwarf

Authors: Andrew Vanderburg, Saul A. Rappaport, Siyi Xu, Ian Crossfield, Juliette C. Becker, Bruce Gary, Felipe Murgas, Simon Blouin, Thomas G. Kaye, Enric Palle, Carl Melis, Brett Morris, Laura Kreidberg, Varoujan Gorjian, Caroline V. Morley, Andrew W. Mann, Hannu Parviainen, Logan A. Pearce, Elisabeth R. Newton, Andreia Carrillo, Ben Zuckerman, Lorne Nelson, Greg Zeimann, Warren R. Brown, René Tronsgaard , et al. (39 additional authors not shown)

Abstract: Astronomers have discovered thousands of planets outside the solar system, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star, but more distant planets can survive this phase and remain in orbit around the white dwarf. Some white dwarfs show evidence for rocky material… ▽ More Astronomers have discovered thousands of planets outside the solar system, most of which orbit stars that will eventually evolve into red giants and then into white dwarfs. During the red giant phase, any close-orbiting planets will be engulfed by the star, but more distant planets can survive this phase and remain in orbit around the white dwarf. Some white dwarfs show evidence for rocky material floating in their atmospheres, in warm debris disks, or orbiting very closely, which has been interpreted as the debris of rocky planets that were scattered inward and tidally disrupted. Recently, the discovery of a gaseous debris disk with a composition similar to ice giant planets demonstrated that massive planets might also find their way into tight orbits around white dwarfs, but it is unclear whether the planets can survive the journey. So far, the detection of intact planets in close orbits around white dwarfs has remained elusive. Here, we report the discovery of a giant planet candidate transiting the white dwarf WD 1856+534 (TIC 267574918) every 1.4 days. The planet candidate is roughly the same size as Jupiter and is no more than 14 times as massive (with 95% confidence). Other cases of white dwarfs with close brown dwarf or stellar companions are explained as the consequence of common-envelope evolution, wherein the original orbit is enveloped during the red-giant phase and shrinks due to friction. In this case, though, the low mass and relatively long orbital period of the planet candidate make common-envelope evolution less likely. Instead, the WD 1856+534 system seems to demonstrate that giant planets can be scattered into tight orbits without being tidally disrupted, and motivates searches for smaller transiting planets around white dwarfs. △ Less

Submitted 15 September, 2020; originally announced September 2020.

Comments: 50 pages, 12 figures, 2 tables. Published in Nature on Sept. 17, 2020. The final authenticated version is available online at: https://www.nature.com/articles/s41586-020-2713-y

Inclination Excitation of Solar System Debris Disk due to Stellar Flybys

Authors: Nathaniel W. H. Moore, Gongjie Li, Fred C. Adams

Abstract: Most stars form in clusters where relatively close encounters with other stars are common and can leave imprints on the orbital architecture of planetary systems. In this paper, we investigate the inclination excitation of debris disk particles due to such stellar encounters. We derive an analytical expression that describes inclination excitation in the hierarchical limit where the stellar flyby… ▽ More Most stars form in clusters where relatively close encounters with other stars are common and can leave imprints on the orbital architecture of planetary systems. In this paper, we investigate the inclination excitation of debris disk particles due to such stellar encounters. We derive an analytical expression that describes inclination excitation in the hierarchical limit where the stellar flyby is distant. We then obtain numerical results for the corresponding particle inclination distribution in the non-hierarchical regime using a large ensemble of N-body simulations. For encounters with expected parameters, we find that the bulk inclination of the disk particles remains low. However, a distinct high inclination population is produced by prograde stellar encounters for particles with final pericenter distances above $50$AU. The maximum extent $i_t$ of the inclination distribution scales with the inclination of the encounter $\sin(i_s)$ for massive star flybys with low incoming velocity. The inclination distribution of observed trans-Neptunian objects places constraints on the dynamical history of our Solar System. For example, these results imply an upper limit on product of the number density $n$ of the solar birth cluster and the Sun's residence time $τ$ of the form $nτ\lesssim8\times10^4$ Myr pc$^{-3}$. Stronger constraints can be derived with future observational surveys of the outer Solar System. △ Less

Submitted 18 August, 2020; v1 submitted 30 July, 2020; originally announced July 2020.

Comments: 15 pages, 12 figures. This paper has been accepted for publication in The Astrophysical Journal, and comments are welcome

A Coupled Analysis of Atmospheric Mass Loss and Tidal Evolution in XUV Irradiated Exoplanets: the TRAPPIST-1 Case Study

Authors: Juliette Becker, Elena Gallo, Edmund Hodges-Kluck, Fred C. Adams, Rory Barnes

Abstract: Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in XUV irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system, and place const… ▽ More Exoplanets residing close to their stars can experience evolution of both their physical structures and their orbits due to the influence of their host stars. In this work, we present a coupled analysis of dynamical tidal dissipation and atmospheric mass loss for exoplanets in XUV irradiated environments. As our primary application, we use this model to study the TRAPPIST-1 system, and place constraints on the interior structure and orbital evolution of the planets. We start by reporting on a UV continuum flux measurement (centered around $\sim1900$ Angstroms) for the star TRAPPIST-1, based on 300 ks of Neil Gehrels Swift Observatory data, and which enables an estimate of the XUV-driven thermal escape arising from XUV photo-dissociation for each planet. We find that the X-ray flaring luminosity, measured from our X-ray detections, of TRAPPIST-1 is 5.6 $\times$10$^{-4} L_{*}$, while the full flux including non-flaring periods is 6.1 $\times$10$^{-5} L_{*}$, when $L_{*}$ is TRAPPIST-1's bolometric luminosity. We then construct a model that includes both atmospheric mass-loss and tidal evolution, and requires the planets to attain their present-day orbital elements during this coupled evolution. We use this model to constrain the ratio $Q'=3Q/2k_{2}$ for each planet. Finally, we use additional numerical models implemented with the Virtual Planet Simulator \texttt{VPLanet} to study ocean retention for these planets using our derived system parameters. △ Less

Submitted 4 May, 2020; originally announced May 2020.

Comments: accepted to AJ. See extra figures, data, and scripts at https://github.com/jxcbecker/trappist1

A triple star system with a misaligned and warped circumstellar disk shaped by disk tearing

Authors: Stefan Kraus, Alexander Kreplin, Alison K. Young, Matthew R. Bate, John D. Monnier, Tim J. Harries, Henning Avenhaus, Jacques Kluska, Anna S. E. Laws, Evan A. Rich, Matthew Willson, Alicia N. Aarnio, Fred C. Adams, Sean M. Andrews, Narsireddy Anugu, Jaehan Bae, Theo ten Brummelaar, Nuria Calvet, Michel Curé, Claire L. Davies, Jacob Ennis, Catherine Espaillat, Tyler Gardner, Lee Hartmann, Sasha Hinkley , et al. (7 additional authors not shown)

Abstract: Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple… ▽ More Young stars are surrounded by a circumstellar disk of gas and dust, within which planet formation can occur. Gravitational forces in multiple star systems can disrupt the disk. Theoretical models predict that if the disk is misaligned with the orbital plane of the stars, the disk should warp and break into precessing rings, a phenomenon known as disk tearing. We present observations of the triple star system GWOrionis, finding evidence for disk tearing. Our images show an eccentric ring that is misaligned with the orbital planes and the outer disk. The ring casts shadows on a strongly warped intermediate region of the disk. If planets can form within the warped disk, disk tearing could provide a mechanism for forming wide-separation planets on oblique orbits. △ Less

Submitted 4 September, 2020; v1 submitted 2 April, 2020; originally announced April 2020.

Comments: 63 pages, 4+13 Figures, 6 Tables, published at Science

Energy Optimization in Binary Star Systems: Explanation for Equal Mass Members in Close Orbits

Authors: Fred C. Adams, Konstantin Batygin, Anthony M. Bloch

Abstract: Observations indicate that members of close stellar binaries often have mass ratios close to unity, while longer-period systems exhibit a more uniform mass-ratio distribution. This paper provides a theoretical explanation for this finding by determining the tidal equilibrium states for binary star systems --- subject to the constraints of conservation of angular momentum and constant total mass. T… ▽ More Observations indicate that members of close stellar binaries often have mass ratios close to unity, while longer-period systems exhibit a more uniform mass-ratio distribution. This paper provides a theoretical explanation for this finding by determining the tidal equilibrium states for binary star systems --- subject to the constraints of conservation of angular momentum and constant total mass. This work generalizes previous treatments by including the mass fraction as a variable in the optimization problem. The results show that the lowest energy state accessible to the system corresponds to equal mass stars on a circular orbit, where the stellar spin angular velocities are both synchronized and aligned with the orbit. These features are roughly consistent with observed properties of close binary systems. We also find the conditions required for this minimum energy state to exist: [1] The total angular momentum must exceed a critical value, [2] the orbital angular momentum must be three times greater than the total spin angular momentum, and [3] the semimajor axis is bounded from above. The last condition implies that sufficiently wide binaries are not optimized with equal mass stars, where the limiting binary separation occurs near $a_0\approx16R_\ast$. △ Less

Submitted 23 March, 2020; originally announced March 2020.

Comments: 13 pages, 2 figures, accepted for publication in MNRAS

Testing the isotropy of the Dark Energy Survey's extreme trans-Neptunian objects

Authors: Pedro H. Bernardinelli, Gary M. Bernstein, Masao Sako, Stephanie Hamilton, David W. Gerdes, Fred C. Adams, William R. Saunders, M. Aguena, S. Allam, S. Avila, D. Brooks, H. T. Diehl, P. Doel, S. Everett, J. García-Bellido, E. Gaztanaga, R. A. Gruendl, K. Honscheid, R. L. C. Ogando, A. Palmese, D. L. Tucker, A. R. Walker

Abstract: We test whether the population of "extreme" trans-Neptunian objects (eTNOs) detected in the Y4 Dark Energy Survey (DES) data exhibit azimuthal asymmetries which might be evidence of gravitational perturbations from an unseen super-Earth in a distant orbit. By rotating the orbits of the detected eTNOs, we construct a synthetic population which, when subject to the DES selection function, reproduces… ▽ More We test whether the population of "extreme" trans-Neptunian objects (eTNOs) detected in the Y4 Dark Energy Survey (DES) data exhibit azimuthal asymmetries which might be evidence of gravitational perturbations from an unseen super-Earth in a distant orbit. By rotating the orbits of the detected eTNOs, we construct a synthetic population which, when subject to the DES selection function, reproduces the detected distribution of eTNOs in the orbital elements $a,e,$ and $i$ as well as absolute magnitude $H$, but has uniform distributions in mean anomaly $M$, longitude of ascending node $Ω,$ and argument of perihelion $ω.$ We then compare the detected distributions in each of $Ω, ω,$ and $\varpi\equivΩ+ω$ to those expected from the isotropic population, using Kuiper's variant of the Kolmogorov-Smirnov test. The three angles are tested for each of 4 definitions of the eTNO population, choosing among $a>(150,250)$ AU and perihelion $q>(30,37)$ AU. These choices yield 3--7 eTNOs in the DES Y4 sample. Among the twelve total tests, two have the likelihood of drawing the observed angles from the isotropic population at $p<0.05.$ The 3 detections at $a>250, q>37$ AU, and the 4 detections at $a>250, q>30$ AU, have $Ω$ distribution with $p=0.03$ of coming from the isotropic construction, but this is not strong evidence of anisotropy given the 12 different tests. The DES data taken on their own are thus consistent with azimuthal isotropy and do not require a "Planet 9" hypothesis. The limited sky coverage and object count mean, however, that the DES data by no means falsify this hypothesis. △ Less

Submitted 17 July, 2020; v1 submitted 19 March, 2020; originally announced March 2020.

Comments: Accepted on PSJ

Report number: FERMILAB-PUB-20-113-AE

Energy Optimization in Extrasolar Planetary Systems: The Transition from Peas-in-a-Pod to Runaway Growth

Authors: Fred C. Adams, Konstantin Batygin, Anthony M. Bloch, Gregory Laughlin

Abstract: Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimizatio… ▽ More Motivated by the trends found in the observed sample of extrasolar planets, this paper determines tidal equilibrium states for forming planetary systems --- subject to conservation of angular momentum, constant total mass, and fixed orbital spacing. In the low-mass limit, valid for superearth-class planets with masses of order $m_{\rm p}\sim10M_\oplus$, previous work showed that energy optimization leads to nearly equal mass planets, with circular orbits confined to a plane. The present treatment generalizes previous results by including the self-gravity of the planetary bodies. For systems with sufficiently large total mass $m_{\scriptstyle T}$ in planets, the optimized energy state switches over from the case of nearly equal mass planets to a configuration where one planet contains most of the material. This transition occurs for a critical mass threshold of approximately $m_{\scriptstyle T} \ge m_{\scriptstyle C}\sim40M_\oplus$ (where the value depends on the semimajor axes of the planetary orbits, the stellar mass, and other system properties). These considerations of energy optimization apply over a wide range of mass scales, from binary stars to planetary systems to the collection of moons orbiting the giant planets in our solar system. △ Less

Submitted 24 February, 2020; originally announced February 2020.

Comments: 16 pages, 5 figures, accepted for publication in MNRAS

The Solar wind prevents re-accretion of debris after Mercury's giant impact

Authors: Christopher Spalding, Fred C. Adams

Abstract: The planet Mercury possesses an anomalously large iron core, and a correspondingly high bulk density. Numerous hypotheses have been proposed in order to explain such a large iron content. A long-standing idea holds that Mercury once possessed a larger silicate mantle which was removed by a giant impact early in the the Solar system's history. A central problem with this idea has been that material… ▽ More The planet Mercury possesses an anomalously large iron core, and a correspondingly high bulk density. Numerous hypotheses have been proposed in order to explain such a large iron content. A long-standing idea holds that Mercury once possessed a larger silicate mantle which was removed by a giant impact early in the the Solar system's history. A central problem with this idea has been that material ejected from Mercury is typically re-accreted onto the planet after a short (~Myr) timescale. Here, we show that the primordial Solar wind would have provided sufficient drag upon ejected debris to remove them from Mercury-crossing trajectories before re-impacting the planet's surface. Specifically, the young Sun likely possessed a stronger wind, fast rotation and strong magnetic field. Depending upon the time of the giant impact, the ram pressure associated with this wind would push particles outward into the Solar system, or inward toward the Sun, on sub-Myr timescales, depending upon the size of ejected debris. Accordingly, the giant impact hypothesis remains a viable pathway toward the removal of planetary mantles, both on Mercury and extrasolar planets, particularly those close to young stars with strong winds. △ Less

Submitted 20 February, 2020; v1 submitted 18 February, 2020; originally announced February 2020.

Comments: 14 pages, 6 Figures, Accepted for publication in The Planetary Science Journal (AAS journals)

Dynamics of Planetary Systems Within Star Clusters: Aspects of the Solar System's Early Evolution

Authors: Konstantin Batygin, Fred C. Adams, Yuri K. Batygin, Erik A. Petigura

Abstract: Most planetary systems -- including our own -- are born within stellar clusters, where interactions with neighboring stars can help shape the system architecture. This paper develops an orbit-averaged formalism to characterize the cluster's mean-field effects as well as the physics of long-period stellar encounters. Our secular approach allows for an analytic description of the dynamical consequen… ▽ More Most planetary systems -- including our own -- are born within stellar clusters, where interactions with neighboring stars can help shape the system architecture. This paper develops an orbit-averaged formalism to characterize the cluster's mean-field effects as well as the physics of long-period stellar encounters. Our secular approach allows for an analytic description of the dynamical consequences of the cluster environment on its constituent planetary systems. We analyze special cases of the resulting Hamiltonian, corresponding to eccentricity evolution driven by planar encounters, as well as hyperbolic perturbations upon dissipative disks. We subsequently apply our results to the early evolution of our solar system, where the cluster's collective potential perturbs the solar system's plane, and stellar encounters act to increase the velocity dispersion of the Kuiper belt. Our results are two-fold: first, we find that cluster effects can alter the mean plane of the solar system by $\lesssim1°$, and are thus insufficient to explain the $ψ\approx6°$ obliquity of the sun. Second, we delineate the extent to which stellar flybys excite the orbital dispersion of the cold classical Kuiper belt, and show that while stellar flybys may grow the cold belt's inclination by the observed amount, the resulting distribution is incompatible with the data. Correspondingly, our calculations place an upper limit on the product of the stellar number density and residence time of the sun in its birth cluster, $η\,τ\lesssim2\times10^4\,$Myr/pc$^3$. △ Less

Submitted 13 February, 2020; originally announced February 2020.

Comments: 30 pages, 16 figures, published in The Astronomical Journal

Journal ref: The Astronomical Journal (2020) volume 159, number 3